Molecular dynamics (MD) simulation is a computational technique that models the physical movements of atoms and molecules over time. By numerically solving Newton’s equations of motion for a system of interacting particles, MD provides atomic-level insight into protein structure, dynamics, and function.
The fundamental input to an MD simulation is a force field, which describes the potential energy of the system as a function of atomic coordinates. Force fields combine bonded terms (bond stretching, angle bending, dihedral torsion) with non-bonded terms (electrostatic and van der Waals interactions). Common force fields include AMBER, CHARMM, OPLS, and GROMOS, each parametrized for specific classes of biomolecules.
The simulation proceeds by integrating the equations of motion using algorithms such as Verlet or leapfrog integration. The time step is typically 1–2 fs, limited by the fastest vibrational motions. Periodic boundary conditions simulate bulk behavior by replicating the system in all directions. Long-range electrostatics are calculated using particle mesh Ewald summation.
Solvent is modeled explicitly with thousands of water molecules or implicitly using continuum dielectric models. Ions are added to neutralize the system and achieve physiological ionic strength. Temperature and pressure are controlled using thermostats and barostats.
Analysis of MD trajectories reveals structural fluctuations, root-mean-square deviation, and root-mean-square fluctuation of atomic positions. Principal component analysis identifies dominant collective motions. Free energy calculations using umbrella sampling, thermodynamic integration, or metadynamics determine binding affinities and conformational free energy landscapes.
MD simulations are used to study protein folding, molecular docking, ligand binding, enzyme catalysis, membrane dynamics, and protein-protein interactions. The simulation timescale has increased from nanoseconds to microseconds and beyond with specialized hardware such as Anton and GPU-accelerated computing. Coarse-grained MD reduces computational cost by grouping atoms into beads, enabling millisecond-scale simulations of large systems.
