La simulation de dynamique moléculaire (MD) est une technique computationnelle qui modélise les mouvements physiques des atomes et des molécules dans le temps. En résolvant numériquement les équations du mouvement de Newton pour un système de particules en interaction, la MD fournit un aperçu au niveau atomique de la structure, de la dynamique et de la fonction des biomolécules.
L’entrée fondamentale d’une simulation MD est un champ de force, qui décrit l’énergie potentielle du système en fonction des coordonnées atomiques. Les champs de force combinent des termes liés (élongation de liaison, flexion d’angle, torsion dièdre) avec des termes non liés (interactions électrostatiques et de van der Waals). Les champs de force courants incluent AMBER, CHARMM, OPLS et GROMOS, chacun paramétré pour des classes spécifiques de biomolécules.
La simulation procède par intégration des équations du mouvement en utilisant des algorithmes tels que l’intégration de Verlet ou leapfrog. Le pas de temps est typiquement de 1 à 2 fs, limité par les mouvements vibrationnels les plus rapides. Les conditions aux limites périodiques simulent le comportement volumique en répliquant le système dans toutes les directions. Les électrostatiques à longue portée sont calculées en utilisant la sommation d’Ewald par maillage de particules.
Le solvant est modélisé explicitement avec des milliers de molécules d’eau ou implicitement en utilisant des modèles diélectriques continus. Des ions sont ajoutés pour neutraliser le système et atteindre la force ionique physiologique. La température et la pression sont contrôlées à l’aide de thermostats et de barostats.
L’analyse des trajectoires MD révèle les fluctuations structurales, l’écart quadratique moyen et la fluctuation quadratique moyenne des positions atomiques. L’analyse en composantes principales identifie les mouvements collectifs dominants. Les calculs d’énergie libre utilisant l’échantillonnage par ombrelle, l’intégration thermodynamique ou la métadynamique déterminent les affinités de liaison et les paysages d’énergie libre conformationnels.
Les simulations MD sont utilisées pour étudier le repliement des protéines, la liaison de ligands, la catalyse enzymatique, la dynamique membranaire et les interactions protéine-protéine. L’échelle de temps des simulations est passée de la nanoseconde à la microseconde et au-delà avec du matériel spécialisé comme Anton et le calcul accéléré par GPU. La MD à grain grossier réduit le coût computationnel en regroupant les atomes en billes, permettant des simulations à l’échelle de la milliseconde de grands systèmes.
