The reason a Langevin diffusion leaves $\nu(x)=e^{-U(x)}$ invariant is because it is *[symmetric or reversible with respect to $\nu$][1]*.  In comparison to general diffusion processes, the ergodic properties of symmetric diffusions are not only easier to analyze, but it is often possible to obtain an explicit formula for the stationary density of a symmetric diffusion.  This symmetry property is a bit more transparent when the Langevin equation is written as $$
dX_t = \nabla \log \nu(X_t) dt + \sqrt{2} dW_t \;. \tag{$\star$}
$$ The symmetry property can then be easily verified by checking that the infinitesimal generator of the Langevin diffusion is symmetric in an $L^2$-inner product weighted by $\nu$.  This property of the Langevin diffusion is analogous to the [detailed balance condition][2] of a Markov chain.  It is also the basis of [standard MCMC methods based on Langevin diffusions][3]. In fact, many random-walk based MCMC methods (like random walk Metropolis) with target probability density proportional to $\nu(x)$, can weakly approximate the corresponding Langevin diffusion ($\star$).  This is remarkable because, e.g., the random walk Metropolis algorithm does not involve the gradient of $U$.  

To read more about this connection, see [the introduction and Theorem 5.2 of this paper][4].


  [1]: https://www.cambridge.org/core/journals/advances-in-applied-probability/article/abs/time-reversible-diffusions/6CEF3A3069A15D58C55FBB06D6F467D8
  [2]: https://cims.nyu.edu/~holmes/teaching/asa19/handout_Lecture3_2019.pdf
  [3]: https://projecteuclid.org/journals/bernoulli/volume-2/issue-4/Exponential-convergence-of-Langevin-distributions-and-their-discrete-approximations/bj/1178291835.full
  [4]: https://epubs.siam.org/doi/pdf/10.1137/130937470