Let $$\begin{align*}dX_t&=-\nabla V_1(X_t)dt+\sqrt2 dB_t \\ dY_t&=-\nabla V_2(Y_t)dt+\sqrt2 dB_t.\end{align*}$$ It is well-known that under suitable conditions of $V_1,V_2$ (e.g. $\lambda$-convexity), the processes $X,Y$ converge to equilibrium $\mu_1,\mu_2$ resepectively, where $$\begin{align*}\mu_i(dx)&=\frac{\exp(-V_i(x))}{\int\exp(-V_i(y))dy}dx\end{align*}$$
Hence, the densities of the stationary distributions are known. However, in general it is quite difficult to bound the Wasserstein distance of measures in terms of their densities, see e.g. here.
In this particular case, i doubt that a meaningful bound can be achieved as you describe it. Consider the case where $$V_i(x)=\theta_i(x-m_i)^2$$
leading to $\mu_i$ being Gaussian with mean $m_i$ and variance $\theta_i^{-1}$.
Now note that, unless $\theta_1=\theta_2$ as well as $m_1=m_2$ (in which case the Wasserstein distance would be $0$ anyway), we have both $$\Vert V_1-V_2\Vert_2^2=\infty, \hspace{1cm} \Vert \nabla V_1-\nabla V_2\Vert_2^2=\infty.$$
Hence i doubt that one would be able to describe the "similarity" between $\mu_1$ and $\mu_2$ using such quantities.
Perhaps another approach is the following: One can view the evolution described by the SDEs as a Gradient Flow of a certain energy functional, namely the relative entropy functional $\text{Ent}_{\mu_i}$, on the Wasserstein space. Perhaps a more promising approach would be to compare $\mu_1$ and $\mu_2$ by some comparison between $\text{Ent}_{\mu_1}$ and $\text{Ent}_{\mu_2}$. For more information on viewing this problem as a Gradient Flow, see the book Lectures on Optimal Transport, specifically Chapter 18.