Lets say we are working in a classical scalar field theory and we have two functional $ F[\phi, \pi](x)$ and $G[\phi, \pi](x)$. In most of the references, starting with two functional the Poisson bracket is defined as $$\{F(x),G(y)\} = \int d^3z \left( \frac{\delta F(x)}{\delta \phi(z)}\frac{\delta G(y)}{\delta \pi(z)} - \frac{\delta F(x)}{\delta \pi(z)}\frac{\delta G(y)}{\delta \phi(z)}\right) . $$
But as explained here the functional derivative $\frac{\delta F}{\delta \phi} $ is a distribution rather than a function, so the previous definition does not make much sense. I was wondering then, if the Poisson bracket can be interpreted as the convolution calculated in $(x-y)$ (in the sense of distributions) between the functional derivatives. This works in case of interest such as $\{\phi(x), \pi(y) \}$ but I'm not sure it can be applied for two generic functional (the dependence $(x-y)$ is not explicit). Is there a proof that the Poisson bracket is a convolution? More in general, can field theories be formulated in a formal way in the sense of distributions?