We attempt to control the Lipschitz constant of the solution vector-field to certain drift-diffusion equations. Utilizing techniques introduced by Kiselev, Nazarov, Volberg and Shterenberg to analyze certain active scalars, we track the evolution of moduli of continuity by drift-diffusion systems. As a consequence, we obtain a criterion for the preservation of such moduli of continuity by regular solutions to the incompressible Navier-Stokes system in any dimension d ≥ 3 and in the absence of physical boundaries. We also consider two related systems: a nonlinear drift-diffusion equation, forced nonlocally by a singular integral operator of order zero, and a linear drift-diffusion equation with pressure and Hölder drift velocity b. For the former, we show that the Lipschitz constant grows at most double exponentially in time, leading to global regularity. For the latter, under a supercritical assumption on b ∈ L 1 t C 0,β x , where β ∈ (0, 1) is arbitrary, we prove that the Lipschitz constant of the solution is logarithmically integrable in time. An immediate corollary is a partial regularity result for L 1 t C 0,βx solutions to the incompressible NSE (any β ∈ (0, 1)). On the other hand, critical or subcritical assumptions on b lead to integrability of the Lipschitz constant. All estimates are homogenous involving only the Hölder semi-norm of b and W 1,∞ norm of initial data.