Molecular Dynamics Simulation

‎Molecular dynamics (MD) simulation is a computer simulation technique used in drug design to model the time-dependent behaviour of biomolecular systems. In this method, the atoms (all atoms simulation) are allowed to interact for a fixed period of time, giving a view of the dynamic evolution of protein–ligand complexes under near-physiological conditions.

‎MD generates trajectories of atoms and molecules, which are determined by numerically solving Newton's equations of motion for a system of interacting particles. The forces governing these movements are calculated using molecular force fields, which account for bonded interactions (bond stretching, angle bending, torsions) and non-bonded interactions (electrostatic and van der Waals forces).

‎

‎In computational drug design and discovery, MD simulations help assess binding stability, conformational flexibility, induced fit effects, allosteric modulation, and binding free energy estimations.

‎

‎In MD simulations, energies can be calculated either with first-principle approaches (quantum mechanics) or by treating the molecules as classical objects resembling the 'ball and stick' model (molecular mechanics). 

‎

‎The quantum simulations incorporates the quantum nature of the chemical bond. The electron density function for the valence electrons that determine bonding interactions is computed using quantum mechanical equations, while the nuclear motion is typically propagated using classical dynamics. ‎Quantum molecular dynamics (MD) simulations are computationally expensive and demanding because they require the calculation of electron density at every simulation step. 

‎‎

‎Classical mechanics calculations are much less computationally expensive than quantum mechanical methods. ‎Classical MD simulates the dynamics of the nuclei, but electrons are not present explicitly in the treatment. The nuclei move on a potential energy surface, and electrons are indirectly introduced through this surface that is solely a function of the atomic positions. In this technique, atoms are treated as classical particles, often conceptualized as soft spheres connected by elastic bonds, with parameters carefully optimized to reproduce experimentally observed or quantum-derived properties. 

‎

‎In computational drug design, we usually deal with large biomolecular systems such as protein–ligand complexes, sometimes containing thousands of atoms. Running fully quantum mechanical simulations on systems of that size is simply not practical. That is why classical molecular dynamics remains the preferred method. Its efficiency and scalability make it possible to monitor binding stability, conformational flexibility, and dynamic molecular interactions over biologically meaningful timescales.

‎

‎Although quantum MD gives deeper insight into electronic structure and chemical reactivity, its computational demands restrict its routine use in large biomolecular simulations. Classical MD, through carefully parameterised force fields, provides a realistic and computationally feasible way to explore molecular behaviour in drug discovery. 

‎

References.

Yan, D., Ma, Y., Chen, X., Deng, S., & Wang, Q. (2025). Molecular dynamics-driven drug discovery: Applications and challenges in early compound discovery. Physical Chemistry Chemical Physics, 27, 12633–12651. https://doi.org/10.1039/D5CP00380F

Liu, H. (2026). Role of molecular dynamics simulations and related methods in drug discovery [Special issue]. International Journal of Molecular Sciences, 27(1). 

AlRawashdeh, S., & Barakat, K. H. (2024). Applications of molecular dynamics simulations in drug discovery. In Methods in Molecular Biology (Vol. 2714) (pp. 127-141). 

Koirala, K., Joshi, K., Adediwura, V., Wang, J., & Miao, Y. (2024). Accelerating molecular dynamics simulations for drug discovery. In Methods in Molecular Biology (Vol. 2714) (pp. 187-202). https://pubmed.ncbi.nlm.nih.gov/37676600/

Chen, M., Jiang, X., Zhang, L., Chen, X., Wen, Y., Gu, Z., … Li, X. (2024). The emergence of machine learning force fields in drug design. Medicinal Research Reviews, 44(3), 1147-1182. https://pubmed.ncbi.nlm.nih.gov/38173298/

Farhadi, B., Beygisangchin, M., Ghamari, N., Jakmunee, J., & Tang, T. (2025). Molecular dynamics simulations of proteins: Computational strategies and roles in medicinal chemistry and drug development. Biological Cybernetics, 119. https://pubmed.ncbi.nlm.nih.gov/41003729/

Guo, J., & Liu, H. (2024). The applications of molecular dynamics simulation in studying protein structure and dynamics. Current Medicinal Chemistry, 31(20), 2839–2840. https://rjpbr.com/0929-8673/article/view/645205