We introduce dynamical versions of loop (or Dyson-Schwinger) equations for large families of two--dimensional interacting particle systems, including Dyson Brownian motion, Nonintersecting Bernoulli/Poisson random walks, $β$--corners processes, uniform and Jack-deformed measures on Gelfand-Tsetlin patterns, Macdonald processes, and $(q,κ)$-distributions on lozenge tilings. Under technical assumptions, we show that the dynamical loop equations lead to Gaussian field type fluctuations. As an application, we compute the limit shape for $(q,κ)$--distributions on lozenge tilings and prove that their height fluctuations converge to the Gaussian Free Field in an appropriate complex structure.