Molecular dynamics

Molecular dynamics (MD) is a computer simulation method for studying the physical movements of atoms and molecules over time. By numerically solving Newton's equations of motion for a system of interacting particles, MD generates trajectories that reveal how a molecular system evolves from a given initial configuration under the influence of interatomic forces.

In structural biology, MD complements experimental methods like X-ray crystallography and cryo-EM by adding the dimension of time that these static methods omit. A crystal structure captures a single conformational snapshot; MD reveals the conformational ensemble, the breathing motions of binding pockets, and the allosteric pathways that transmit signals across a protein. The method has become essential for understanding how proteins fold, how drugs bind, and how mutations alter molecular function.

The Achilles' heel of MD is the timescale problem. The characteristic time step in an all-atom simulation is a femtosecond (10⁻¹⁵ seconds), while many biologically relevant processes — protein folding, drug unbinding, allosteric transitions — occur on millisecond to second timescales. Bridging this gap of twelve orders of magnitude requires specialized hardware (Anton), coarse-grained models, or enhanced sampling algorithms that sacrifice atomic detail for temporal reach.

Molecular dynamics is the closest thing structural biology has to a time machine, but it is a time machine built on approximations. The force fields that drive these simulations are parameterized against experimental data, meaning MD does not discover new physics — it extrapolates from known physics. The question is whether extrapolation across twelve orders of magnitude is insight or sophisticated interpolation.