The Yukawa Lagrangian
$$\mathcal{L}=-\dfrac{1}{4}F_{\mu\nu}F^{\mu\nu}+\dfrac{1}{2}m^2A_\mu A^\mu-A_\mu J^\mu$$
can be used to derive the equations of motion
$$\Box A_\mu=J_\mu-m^2A_\mu$$
for the four-potential $A_\mu$. It turns out that if we have $J_\mu$ the current of one charge $e$ at rest at the origin, that is, $J_0(\mathbf{x})=e\delta^{(3)}(\mathbf{x})$ and $J_i=0$, then we can solve the equation by Fourier transform and get
$$A_0(r)=\dfrac{e}{4\pi r} e^{-mr}.$$
It is that then that this potential could explain how protons are bound in the nucleus.
How can we see that? Because honestly I don't see it. I mean, is it by computing the associated force?
In that case I computed it and we would have (assuming the force works as in electrostatics $F_Y = qA_0(r)$.
$$F_Y=\dfrac{e^2}{4\pi r^2}e^{-mr}+\dfrac{me^2}{4\pi r}e^{-mr}.$$
Now, why would this force explain how the protons are bound? I believe this has to do with some specific value for $m$ that allows this force to overcome the repulsion due to the coulomb force
$$F_C=\dfrac{e^2}{4\pi r^2},$$
but I don't know how to make this precise.