{"id":4853,"date":"2026-08-11T12:27:43","date_gmt":"2026-08-11T03:27:43","guid":{"rendered":"https:\/\/bdsl.jbnu.ac.kr\/blog\/?p=4853"},"modified":"2026-08-11T14:36:15","modified_gmt":"2026-08-11T05:36:15","slug":"introduction-to-molecular-dynamics-simulation","status":"publish","type":"post","link":"https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/","title":{"rendered":"Introduction to Molecular Dynamics Simulation"},"content":{"rendered":"\n<h2 class=\"wp-block-heading\">Introduction<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Molecular dynamics (MD) simulation is a computational method for investigating the dynamic behavior of molecules, particularly biomolecules such as proteins, nucleic acids, lipids, and their complexes.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">In an MD simulation, a molecular system is represented as a collection of atoms whose interactions are described using a predefined mathematical model called a <strong>force field<\/strong>. The biomolecule is usually placed in a simulation box containing solvent molecules, such as water, together with appropriate ions.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The forces acting on each atom are calculated from the force field, and the motion of the atoms is simulated by numerically solving Newton&#8217;s equations of motion. By repeatedly updating the positions and velocities of atoms over very short time intervals, the simulation generates a <strong>trajectory<\/strong> describing how the molecular system changes over time.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Experimental methods such as X-ray crystallography and cryo-electron microscopy often provide structural snapshots of biomolecules. MD simulation complements these experiments by providing information about molecular motion, flexibility, conformational transitions, interactions, and possible mechanisms of biological function.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">However, an MD trajectory should not simply be regarded as a molecular movie. Under appropriate conditions and with sufficient sampling, it can also be used to characterize the <strong>distribution of conformational states<\/strong> accessible to the molecular system.<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>Discussion: What does an experimentally determined structure represent?<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Is the structure deposited in the PDB necessarily the equilibrium structure of the molecule in solution? Experimental structures are obtained under particular experimental conditions and may represent one or several conformational states.<\/p>\n<\/blockquote>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Classical Molecular Dynamics<\/h1>\n\n\n\n<h2 class=\"wp-block-heading\">Classical Mechanics<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Conventional biomolecular MD simulations are based primarily on <strong>classical mechanics<\/strong>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Atoms are treated as classical particles with properties such as mass, position, velocity, and partial electric charge. Their interactions are described using mathematical functions rather than by explicitly calculating electrons and molecular orbitals from quantum mechanics.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The force field therefore provides an <strong>effective approximation<\/strong> of the underlying quantum mechanical interactions.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">This approximation makes it possible to simulate systems containing hundreds of thousands or even millions of atoms over relatively long timescales.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">However, conventional classical MD generally assumes a fixed chemical topology. Consequently, processes involving substantial changes in electronic structure are usually not represented directly.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Examples include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>breaking chemical bonds<\/li>\n\n\n\n<li>forming new chemical bonds<\/li>\n\n\n\n<li>proton transfer<\/li>\n\n\n\n<li>electron transfer<\/li>\n\n\n\n<li>chemical reactions<\/li>\n\n\n\n<li>substantial changes in electronic polarization<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">Such problems may require quantum mechanical calculations, reactive force fields, or hybrid quantum mechanics\/molecular mechanics approaches such as QM\/MM.<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>Discussion: What cannot conventional classical MD simulate?<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Consider enzymatic reactions, proton transfer, metal coordination changes, electron transfer, and changes in protonation state.<\/p>\n<\/blockquote>\n\n\n\n<h3 class=\"wp-block-heading\">Material to study<\/h3>\n\n\n\n<ul class=\"wp-block-list\">\n<li><a href=\"Classical mechanics\">Classical mechanics<\/a><\/li>\n\n\n\n<li><a href=\"Newtonian mechanics\">Newtonian mechanics<\/a><\/li>\n\n\n\n<li><a href=\"Force field\">Force field<\/a><\/li>\n\n\n\n<li><a href=\"Quantum mechanics\">Quantum mechanics<\/a><\/li>\n\n\n\n<li><a href=\"QM-MM\">QM-MM<\/a><br><\/li>\n<\/ul>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Molecular Representation<\/h1>\n\n\n\n<h2 class=\"wp-block-heading\">Atoms<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">In classical MD, each atom is represented using a set of parameters.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">These commonly include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>atomic mass<\/li>\n\n\n\n<li>Cartesian coordinates<\/li>\n\n\n\n<li>velocity<\/li>\n\n\n\n<li>partial electric charge<\/li>\n\n\n\n<li>Lennard-Jones parameters<\/li>\n\n\n\n<li>atom type<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">An atom in a force field should therefore not be interpreted as a literal hard sphere containing an explicit nucleus and electrons.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Instead, the parameters assigned to each atom are designed so that the collection of atoms reproduces important molecular properties observed experimentally or obtained from quantum mechanical calculations.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The apparent size of an atom in an MD simulation arises mainly from its non-bonded interaction parameters rather than from a rigid physical boundary.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Molecular Topology<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">A biomolecule is represented as a set of atoms together with information describing how those atoms are chemically connected.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">This information is called the <strong>molecular topology<\/strong>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The topology defines properties such as:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>which atoms are bonded<\/li>\n\n\n\n<li>bond types<\/li>\n\n\n\n<li>bond angles<\/li>\n\n\n\n<li>dihedral relationships<\/li>\n\n\n\n<li>atom types<\/li>\n\n\n\n<li>partial charges<\/li>\n\n\n\n<li>parameters describing the interactions<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">The topology usually remains fixed during a conventional MD simulation.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">The Force Field<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">A <strong>force field<\/strong> is a collection of mathematical functions and parameters used to calculate the potential energy of a molecular system.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Examples of commonly used biomolecular force fields include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>CHARMM<\/li>\n\n\n\n<li>AMBER<\/li>\n\n\n\n<li>OPLS<\/li>\n\n\n\n<li>GROMOS<br><\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">Different force fields use somewhat different functional forms and parameterization strategies, but they generally describe similar classes of molecular interactions.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The potential energy of the system can broadly be divided into:<\/p>\n\n\n\n<ol class=\"wp-block-list\">\n<li>bonded interactions<\/li>\n\n\n\n<li>non-bonded interactions<\/li>\n<\/ol>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Bonded Interactions<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Atoms that are chemically connected require potential-energy terms that maintain realistic molecular geometry.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Bond stretching<\/h3>\n\n\n\n<p class=\"wp-block-paragraph\">A chemical bond is commonly approximated as a harmonic oscillator.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A simplified bond potential can be written as:<\/p>\n\n\n\n<div class=\"wp-block-math\"><math display=\"block\"><semantics><mrow><msub><mi>V<\/mi><mrow><mtext><\/mtext><mi>bond<\/mi><\/mrow><\/msub><mo>=<\/mo><mfrac><mn>1<\/mn><mn>2<\/mn><\/mfrac><msub><mi>k<\/mi><mi>b<\/mi><\/msub><mo form=\"prefix\" stretchy=\"false\">(<\/mo><mi>r<\/mi><mo>\u2212<\/mo><msub><mi>r<\/mi><mn>0<\/mn><\/msub><msup><mo form=\"postfix\" stretchy=\"false\">)<\/mo><mn>2<\/mn><\/msup><\/mrow><annotation encoding=\"application\/x-tex\">V_{\\mathrm{bond}} = \\frac{1}{2}k_b(r-r_0)^2<\/annotation><\/semantics><\/math><\/div>\n\n\n\n<p class=\"wp-block-paragraph\">where:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>(r) is the current bond length,<\/li>\n\n\n\n<li>(r_0) is the equilibrium bond length,<\/li>\n\n\n\n<li>(k_b) is the bond force constant.<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">The energy therefore increases when a bond is stretched or compressed away from its preferred length.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Angle Bending<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Three connected atoms define a bond angle.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Deviation from the preferred angle is also commonly represented using a harmonic potential:<\/p>\n\n\n\n<div class=\"wp-block-math\"><math display=\"block\"><semantics><mrow><msub><mi>V<\/mi><mrow><mtext><\/mtext><mi>angle<\/mi><\/mrow><\/msub><mo>=<\/mo><mfrac><mn>1<\/mn><mn>2<\/mn><\/mfrac><msub><mi>k<\/mi><mi>\u03b8<\/mi><\/msub><mo form=\"prefix\" stretchy=\"false\">(<\/mo><mi>\u03b8<\/mi><mo>\u2212<\/mo><msub><mi>\u03b8<\/mi><mn>0<\/mn><\/msub><msup><mo form=\"postfix\" stretchy=\"false\">)<\/mo><mn>2<\/mn><\/msup><\/mrow><annotation encoding=\"application\/x-tex\">V_{\\mathrm{angle}} = \\frac{1}{2}k_\\theta(\\theta-\\theta_0)^2<\/annotation><\/semantics><\/math><\/div>\n\n\n\n<p class=\"wp-block-paragraph\">This term helps maintain realistic local molecular geometry.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Dihedral Angles<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Four sequentially connected atoms define a <strong>dihedral angle<\/strong>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Rotation around chemical bonds is extremely important because it allows biomolecules to adopt different conformations.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Dihedral potentials describe the energetic preference for particular rotational states.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">For example, rotation around many single bonds is possible but not equally favorable at all angles.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Double bonds, aromatic systems, and other chemically restricted structures are represented using parameters that strongly favor particular geometries.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Classical force fields do not normally describe explicit sigma and pi electrons. Instead, the underlying electronic properties of these chemical bonds are incorporated indirectly into the force-field parameters.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Improper Dihedrals<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Improper dihedral terms are commonly used to maintain properties such as:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>planarity<\/li>\n\n\n\n<li>chirality<\/li>\n\n\n\n<li>appropriate geometry around particular functional groups<\/li>\n<\/ul>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">CHARMM36 Bonded Terms<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">For example, the CHARMM36 force field includes several bonded interaction terms such as:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>bond stretching<\/li>\n\n\n\n<li>angle bending<\/li>\n\n\n\n<li>Urey-Bradley interactions<\/li>\n\n\n\n<li>proper dihedrals<\/li>\n\n\n\n<li>improper dihedrals<\/li>\n<\/ul>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>Discussion: Why are bonds, angles, and dihedrals represented separately?<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">What structural property is controlled by each term, and why would bond lengths alone be insufficient to maintain a realistic molecular structure?<\/p>\n<\/blockquote>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Non-Bonded Interactions<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">Atoms also interact with atoms to which they are not directly bonded.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The two major non-bonded interactions in conventional biomolecular MD are:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>van der Waals interactions<\/li>\n\n\n\n<li>electrostatic interactions<\/li>\n<\/ul>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Lennard-Jones Potential<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Van der Waals interactions are commonly modeled using the <strong>Lennard-Jones potential<\/strong>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A typical form is:<\/p>\n\n\n\n<div class=\"wp-block-math\"><math display=\"block\"><semantics><mrow><msub><mi>V<\/mi><mrow><mtext><\/mtext><mi>LJ<\/mi><\/mrow><\/msub><mo form=\"prefix\" stretchy=\"false\">(<\/mo><mi>r<\/mi><mo form=\"postfix\" stretchy=\"false\">)<\/mo><mo>=<\/mo><mn>4<\/mn><mi>\u03f5<\/mi><mrow><mo fence=\"true\" form=\"prefix\">[<\/mo><msup><mrow><mo fence=\"true\" form=\"prefix\">(<\/mo><mfrac><mi>\u03c3<\/mi><mi>r<\/mi><\/mfrac><mo fence=\"true\" form=\"postfix\">)<\/mo><\/mrow><mn>12<\/mn><\/msup><mo>+<\/mo><msup><mrow><mo fence=\"true\" form=\"prefix\">(<\/mo><mfrac><mi>\u03c3<\/mi><mi>r<\/mi><\/mfrac><mo fence=\"true\" form=\"postfix\">)<\/mo><\/mrow><mn>6<\/mn><\/msup><mo fence=\"true\" form=\"postfix\">]<\/mo><\/mrow><\/mrow><annotation encoding=\"application\/x-tex\">V_{\\mathrm{LJ}}(r) = 4\\epsilon \\left[ \\left(\\frac{\\sigma}{r}\\right)^{12} +\\left(\\frac{\\sigma}{r}\\right)^6 \\right] <\/annotation><\/semantics><\/math><\/div>\n\n\n\n<p class=\"wp-block-paragraph\">At intermediate distances, the attractive term approximates dispersion interactions.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">At very short distances, the potential becomes strongly repulsive.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">This repulsive term prevents atoms from unrealistically occupying the same region of space and approximates short-range quantum mechanical effects associated with overlapping electron clouds.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">It should not be interpreted simply as an elastic collision between rigid spheres.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Electrostatic Interactions<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Atoms in classical force fields are assigned partial electric charges.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The electrostatic interaction between two charged atoms is approximately described by Coulomb&#8217;s law:<\/p>\n\n\n\n<div class=\"wp-block-math\"><math display=\"block\"><semantics><mrow><msub><mi>V<\/mi><mrow><mtext><\/mtext><mi>Coulomb<\/mi><\/mrow><\/msub><mo form=\"prefix\" stretchy=\"false\">(<\/mo><mi>r<\/mi><mo form=\"postfix\" stretchy=\"false\">)<\/mo><mo>=<\/mo><mfrac><mn>1<\/mn><mrow><mn>4<\/mn><mi>\u03c0<\/mi><msub><mi>\u03f5<\/mi><mn>0<\/mn><\/msub><\/mrow><\/mfrac><mfrac><mrow><msub><mi>q<\/mi><mi>i<\/mi><\/msub><msub><mi>q<\/mi><mi>j<\/mi><\/msub><\/mrow><mi>r<\/mi><\/mfrac><\/mrow><annotation encoding=\"application\/x-tex\">V_{\\mathrm{Coulomb}}(r) = \\frac{1}{4\\pi\\epsilon_0} \\frac{q_iq_j}{r} <\/annotation><\/semantics><\/math><\/div>\n\n\n\n<p class=\"wp-block-paragraph\">Electrostatic interactions are particularly important for:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>salt bridges<\/li>\n\n\n\n<li>hydrogen-bonding environments<\/li>\n\n\n\n<li>protein-ligand recognition<\/li>\n\n\n\n<li>protein-DNA interactions<\/li>\n\n\n\n<li>ion interactions<\/li>\n\n\n\n<li>solvent organization<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">Because electrostatic interactions are long-range, their calculation requires special treatment in periodic MD simulations.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Particle Mesh Ewald (PME) is commonly used for this purpose.<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>Discussion: How are long-range interactions calculated efficiently?<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Topics to explore:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>interaction cutoffs<\/li>\n\n\n\n<li>neighbor lists<\/li>\n\n\n\n<li>periodic boundaries<\/li>\n\n\n\n<li>Ewald summation<\/li>\n\n\n\n<li>Particle Mesh Ewald<\/li>\n<\/ul>\n<\/blockquote>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Potential Energy, Force, and Motion<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">The force field defines the <strong>potential energy<\/strong> (V) of the molecular system.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The force acting on atom (i) is calculated from the gradient of this potential energy:<\/p>\n\n\n\n<div class=\"wp-block-math\"><math display=\"block\"><semantics><mrow><msub><mi>\ud835\udc05<\/mi><mi>i<\/mi><\/msub><mo>=<\/mo><mo form=\"prefix\" stretchy=\"false\">\u2212<\/mo><msub><mo form=\"prefix\" stretchy=\"false\">\u2207<\/mo><mi>i<\/mi><\/msub><mi>V<\/mi><\/mrow><annotation encoding=\"application\/x-tex\">\\mathbf{F}_i = -\\nabla_i V<\/annotation><\/semantics><\/math><\/div>\n\n\n\n<p class=\"wp-block-paragraph\">The force is then related to acceleration through Newton&#8217;s second law:<\/p>\n\n\n\n<div class=\"wp-block-math\"><math display=\"block\"><semantics><mrow><msub><mi>\ud835\udc05<\/mi><mi>i<\/mi><\/msub><mo>=<\/mo><msub><mi>m<\/mi><mi>i<\/mi><\/msub><msub><mi>\ud835\udc1a<\/mi><mi>i<\/mi><\/msub><\/mrow><annotation encoding=\"application\/x-tex\">\\mathbf{F}_i = m_i\\mathbf{a}_i<\/annotation><\/semantics><\/math><\/div>\n\n\n\n<p class=\"wp-block-paragraph\">Thus, the conceptual sequence of an MD simulation is:<\/p>\n\n\n\n<div class=\"wp-block-math\"><math display=\"block\"><semantics><mrow><mtext>atomic&nbsp;coordinates<\/mtext><mo stretchy=\"false\">\u2192<\/mo><mtext>potential&nbsp;energy<\/mtext><mo stretchy=\"false\">\u2192<\/mo><mtext>forces<\/mtext><mo stretchy=\"false\">\u2192<\/mo><mtext>acceleration<\/mtext><mo stretchy=\"false\">\u2192<\/mo><mtext>new&nbsp;positions&nbsp;and&nbsp;velocities<\/mtext><\/mrow><annotation encoding=\"application\/x-tex\">\\text{atomic coordinates} \\rightarrow \\text{potential energy} \\rightarrow \\text{forces} \\rightarrow \\text{acceleration} \\rightarrow \\text{new positions and velocities}<\/annotation><\/semantics><\/math><\/div>\n\n\n\n<p class=\"wp-block-paragraph\">This procedure is repeated many times throughout the simulation.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Kinetic Energy and Temperature<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">Atoms possess velocity and therefore have <strong>kinetic energy<\/strong>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">For a single atom:<\/p>\n\n\n\n<div class=\"wp-block-math\"><math display=\"block\"><semantics><mrow><msub><mi>E<\/mi><mi>k<\/mi><\/msub><mo>=<\/mo><mfrac><mn>1<\/mn><mn>2<\/mn><\/mfrac><mi>m<\/mi><msup><mi>v<\/mi><mn>2<\/mn><\/msup><\/mrow><annotation encoding=\"application\/x-tex\">E_k = \\frac{1}{2}mv^2<\/annotation><\/semantics><\/math><\/div>\n\n\n\n<p class=\"wp-block-paragraph\">The total kinetic energy of the molecular system is related statistically to its temperature.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">At finite temperature, atoms continuously move and exchange energy through their interactions with surrounding atoms.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The total energy of a system can broadly be expressed as:<\/p>\n\n\n\n<div class=\"wp-block-math\"><math display=\"block\"><semantics><mrow><msub><mi>E<\/mi><mrow><mtext><\/mtext><mi>total<\/mi><\/mrow><\/msub><mo>=<\/mo><msub><mi>E<\/mi><mrow><mtext><\/mtext><mi>kinetic<\/mi><\/mrow><\/msub><mo>+<\/mo><msub><mi>E<\/mi><mrow><mtext><\/mtext><mi>potential<\/mi><\/mrow><\/msub><\/mrow><annotation encoding=\"application\/x-tex\">E_{\\mathrm{total}} =E_{\\mathrm{kinetic}} + E_{\\mathrm{potential}}<\/annotation><\/semantics><\/math><\/div>\n\n\n\n<p class=\"wp-block-paragraph\">In an isolated idealized system, total energy should be conserved.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">However, many biomolecular simulations intentionally exchange energy or volume with external algorithms in order to maintain a desired temperature or pressure.<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>Discussion: Potential energy, kinetic energy, and total energy<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">How does temperature relate to kinetic energy? What should happen to total energy in an NVE simulation compared with an NVT simulation?<\/p>\n<\/blockquote>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Numerical Integration and Time Steps<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">The equations of motion cannot generally be solved analytically for a large molecular system.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Instead, MD programs calculate the trajectory using <strong>numerical integration<\/strong>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Common integration algorithms include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>Verlet<\/li>\n\n\n\n<li>velocity Verlet<\/li>\n\n\n\n<li>leap-frog<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">The system is advanced in very small discrete time intervals called <strong>time steps<\/strong>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">For atomistic biomolecular simulations, a typical time step is approximately:<\/p>\n\n\n\n<div class=\"wp-block-math\"><math display=\"block\"><semantics><mrow><mn>1<\/mn><mo>\u2212<\/mo><mn>2<\/mn><mtext>&nbsp;<\/mtext><mrow><mtext><\/mtext><mi>fs<\/mi><\/mrow><\/mrow><annotation encoding=\"application\/x-tex\">1-2\\ \\mathrm{fs} <\/annotation><\/semantics><\/math><\/div>\n\n\n\n<p class=\"wp-block-paragraph\">where:<\/p>\n\n\n\n<div class=\"wp-block-math\"><math display=\"block\"><semantics><mrow><mn>1<\/mn><mtext>&nbsp;<\/mtext><mrow><mtext><\/mtext><mi>fs<\/mi><\/mrow><mo>=<\/mo><msup><mn>10<\/mn><mrow><mo lspace=\"0em\" rspace=\"0em\">\u2212<\/mo><mn>15<\/mn><\/mrow><\/msup><mtext>&nbsp;<\/mtext><mrow><mi mathvariant=\"normal\">s<\/mi><\/mrow><\/mrow><annotation encoding=\"application\/x-tex\">1\\ \\mathrm{fs}=10^{-15}\\ \\mathrm{s}<\/annotation><\/semantics><\/math><\/div>\n\n\n\n<p class=\"wp-block-paragraph\">The time step must be short enough to accurately describe the fastest important molecular motions.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Bond vibrations involving hydrogen atoms are particularly fast. Constraints such as <strong>LINCS<\/strong> can be used to restrict selected bond vibrations, allowing somewhat longer time steps.<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>Discussion: Why are MD time steps so short?<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Explore the relationship between molecular vibration frequencies, hydrogen-containing bonds, numerical stability, and constraint algorithms such as LINCS.<\/p>\n<\/blockquote>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Simulation Environment<\/h1>\n\n\n\n<h2 class=\"wp-block-heading\">Solvent<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Biomolecules normally function in an aqueous environment.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A protein is therefore typically placed inside a simulation box containing explicit water molecules.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Common water models include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>TIP3P<\/li>\n\n\n\n<li>TIP4P<\/li>\n\n\n\n<li>SPC<\/li>\n\n\n\n<li>SPC\/E<br><\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">The choice of water model is usually made together with the selected force field.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Ions may also be added to:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>neutralize the total charge of the system<\/li>\n\n\n\n<li>approximate physiological ionic strength<br><\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">For example, a simulation may contain approximately 0.15 M NaCl.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Periodic Boundary Conditions<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">A simulation can contain only a finite number of water molecules.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">If the edge of the simulation box represented a physical wall, molecules near that boundary would experience an artificial environment.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">This problem is commonly avoided using <strong>periodic boundary conditions<\/strong>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The simulation box is conceptually repeated in all directions.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">When a molecule leaves one side of the box, an equivalent image appears from the opposite side.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">This creates an approximately bulk-like environment without requiring an infinitely large number of solvent molecules.<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>Discussion: Periodic boundary conditions<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">What does it mean for a simulation box to be periodically repeated? What artifacts can occur if the box is too small?<\/p>\n<\/blockquote>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Temperature Control<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">Most biomolecular simulations are performed at a specified temperature, such as:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>approximately 298 K for room temperature<\/li>\n\n\n\n<li>approximately 310 K for physiological temperature<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">A <strong>thermostat<\/strong> is used to generate or maintain the desired temperature distribution.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A thermostat should not simply be regarded as an external environment that manually adds energy whenever the system becomes cold.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Rather, it is a mathematical algorithm designed to control the statistical distribution of kinetic energy.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Common thermostats include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>velocity-rescaling thermostat<\/li>\n\n\n\n<li>Nos\u00e9-Hoover thermostat<\/li>\n<\/ul>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>Discussion: What does a thermostat really do?<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A thermostat does not necessarily force the instantaneous temperature to equal the reference temperature. Temperature naturally fluctuates. The important question is whether the thermostat generates the appropriate statistical ensemble.<\/p>\n<\/blockquote>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Pressure Control<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">Biomolecules are commonly studied at approximately atmospheric pressure.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A <strong>barostat<\/strong> is used to control pressure.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A barostat typically changes the volume and dimensions of the simulation box so that the system samples states consistent with the specified pressure.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Examples include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>Berendsen barostat<\/li>\n\n\n\n<li>C-rescale barostat<\/li>\n\n\n\n<li>Parrinello-Rahman barostat<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">Different barostats have different properties and are not necessarily appropriate for every phase of a simulation.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Statistical Ensembles<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">MD simulations can be performed under different thermodynamic conditions.<\/p>\n\n\n\n<h2 class=\"wp-block-heading\">NVE Ensemble<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">The NVE ensemble maintains:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>N: number of particles<\/li>\n\n\n\n<li>V: volume<\/li>\n\n\n\n<li>E: total energy<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">No thermostat or barostat is used in an ideal NVE simulation.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">NVT Ensemble<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">The NVT ensemble maintains:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>N: number of particles<\/li>\n\n\n\n<li>V: volume<\/li>\n\n\n\n<li>T: temperature<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">A thermostat is used to control temperature while the box volume remains fixed.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">NPT Ensemble<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">The NPT ensemble maintains:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>N: number of particles<\/li>\n\n\n\n<li>P: pressure<\/li>\n\n\n\n<li>T: temperature<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">A thermostat and barostat are used.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">For biomolecules in solution, NPT conditions are commonly used because they allow the solvent density and simulation-box volume to adjust at the desired temperature and pressure.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Typical MD Workflow<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">A conventional protein MD simulation commonly includes the following stages.<\/p>\n\n\n\n<h2 class=\"wp-block-heading\">1. Structure Preparation<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">The starting molecular structure is prepared.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">This may include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>removing or retaining crystallographic molecules<\/li>\n\n\n\n<li>adding missing atoms<\/li>\n\n\n\n<li>resolving alternate conformations<\/li>\n\n\n\n<li>assigning protonation states<\/li>\n\n\n\n<li>adding hydrogen atoms<\/li>\n\n\n\n<li>parameterizing ligands<\/li>\n\n\n\n<li>checking unusual residues<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">The quality of the starting structure strongly affects the simulation.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">2. Solvation<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">The biomolecule is placed in a simulation box and surrounded by water molecules.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">3. Ion Addition<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Counterions are added to neutralize the system.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Additional ions may be introduced to reproduce an approximate experimental ionic strength.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">4. Energy Minimization<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">The initial molecular system may contain unfavorable contacts or steric clashes.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Energy minimization adjusts atomic coordinates to reduce these severe interactions before dynamic simulation begins.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Energy minimization is not itself an MD simulation because it does not represent physical motion over time.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">5. NVT Equilibration<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">The system is gradually brought to the desired temperature while maintaining a fixed simulation-box volume.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Position restraints may initially be applied to the biomolecule while the solvent reorganizes around it.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">6. NPT Equilibration<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Pressure control is introduced.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The simulation box and solvent density can then adjust toward conditions appropriate for the target temperature and pressure.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">7. Production MD<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">After the system has been sufficiently prepared and equilibrated, a production simulation is performed.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The resulting trajectory is used for subsequent scientific analysis.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A common simplified workflow is therefore:<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">\\text{Preparation} \\rightarrow \\text{Energy minimization} \\rightarrow \\text{NVT} \\rightarrow \\text{NPT} \\rightarrow \\text{Production MD}<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Trajectories and Molecular Ensembles<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">During a simulation, atomic coordinates are periodically saved.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The resulting collection of structures is called a <strong>trajectory<\/strong>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A trajectory can be visualized as a molecular movie, but this interpretation is incomplete.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">At finite temperature, a biomolecule does not possess only one structure. Instead, it occupies an <strong>ensemble of conformational states<\/strong>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Some states may be highly populated, whereas others may occur only rarely.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">MD simulation attempts to sample these states over time.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Consequently, many biologically meaningful quantities are better described as distributions or ensemble averages rather than as properties of one particular structure.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Examples include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>average distances<\/li>\n\n\n\n<li>probability of a hydrogen bond<\/li>\n\n\n\n<li>distribution of dihedral angles<\/li>\n\n\n\n<li>frequency of loop opening<\/li>\n\n\n\n<li>conformational populations<\/li>\n\n\n\n<li>residue flexibility<\/li>\n<\/ul>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>Discussion: Trajectory versus ensemble<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A trajectory is one realization of molecular motion over time. Under what assumptions can observations from a finite trajectory be interpreted as properties of the molecular ensemble?<\/p>\n<\/blockquote>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Sampling and Timescale<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">One of the major limitations of MD simulation is <strong>sampling<\/strong>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Molecular processes occur over a wide range of timescales.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Examples may include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>bond vibrations: femtoseconds<\/li>\n\n\n\n<li>side-chain rotations: picoseconds to nanoseconds<\/li>\n\n\n\n<li>loop rearrangements: nanoseconds to microseconds or longer<\/li>\n\n\n\n<li>domain motion: microseconds to milliseconds<\/li>\n\n\n\n<li>protein folding: microseconds to seconds or longer<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">A simulation can only observe processes that occur within the accessible simulation time and are sampled sufficiently frequently.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Therefore, the absence of a conformational transition in a simulation does not demonstrate that the transition is biologically impossible.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Similarly, observing a transition only once may not be sufficient to estimate its probability.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Multiple independent simulations, longer trajectories, or enhanced-sampling methods may therefore be required.<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>Discussion: Simulation time versus biological time<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Does a 100-ns simulation adequately sample a protein&#8217;s conformational behavior? The answer depends on the molecular process being investigated.<\/p>\n<\/blockquote>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Assessing a Simulation<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">Before interpreting biological mechanisms from an MD trajectory, the quality and behavior of the simulation should be examined.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Several different concepts should be distinguished:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>numerical stability<\/li>\n\n\n\n<li>thermodynamic stability<\/li>\n\n\n\n<li>structural equilibration<\/li>\n\n\n\n<li>sampling<\/li>\n\n\n\n<li>convergence<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">These are related but are not equivalent.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Temperature, Pressure, and Energy<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Basic simulation properties should first be inspected.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">These may include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>temperature<\/li>\n\n\n\n<li>pressure<\/li>\n\n\n\n<li>density<\/li>\n\n\n\n<li>potential energy<\/li>\n\n\n\n<li>kinetic energy<\/li>\n\n\n\n<li>total energy<\/li>\n\n\n\n<li>simulation-box volume<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">These measurements help identify major problems such as numerical instability or inappropriate equilibration.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">However, stable temperature and pressure do not necessarily mean that the biomolecule itself has fully equilibrated.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Root Mean Square Deviation<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">One of the most widely used structural measurements is the <strong>root mean square deviation (RMSD)<\/strong>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">RMSD measures the structural difference between two sets of atomic coordinates, generally after appropriate structural alignment.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">For example, the protein structure at each simulation time point can be compared with the initial structure.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">RMSD can indicate whether the protein:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>remains close to its initial conformation<\/li>\n\n\n\n<li>relaxes into a different structural state<\/li>\n\n\n\n<li>undergoes large conformational transitions<\/li>\n\n\n\n<li>shows possible structural instability<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">However, RMSD must be interpreted carefully.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A plateau in RMSD does <strong>not<\/strong> prove that the simulation has reached thermodynamic equilibrium or fully converged.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Likewise, a large RMSD does not necessarily indicate a failed simulation.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A large change may represent:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>normal relaxation from the experimental starting structure<\/li>\n\n\n\n<li>loop rearrangement<\/li>\n\n\n\n<li>domain motion<\/li>\n\n\n\n<li>ligand-induced conformational change<\/li>\n\n\n\n<li>transition between functional states<\/li>\n\n\n\n<li>unfolding<\/li>\n\n\n\n<li>simulation artifact<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">The biological context and structural details should therefore be examined.<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>Discussion: Does a stable RMSD mean convergence?<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Usually not. A stable RMSD shows only that the overall deviation from the selected reference remains within a certain range.<\/p>\n<\/blockquote>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Stability, Equilibration, and Convergence<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">These concepts should be distinguished.<\/p>\n\n\n\n<h2 class=\"wp-block-heading\">Stability<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">A simulation is stable when it does not show obvious numerical or structural failure.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Examples of instability include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>rapidly increasing energy<\/li>\n\n\n\n<li>unrealistic temperature<\/li>\n\n\n\n<li>severe structural distortion<\/li>\n\n\n\n<li>numerical integration failure<\/li>\n<\/ul>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Equilibration<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">A system can be considered equilibrated with respect to a particular property when that property no longer shows a strong systematic dependence on the initial preparation.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Different components may equilibrate on different timescales.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">For example:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>temperature may equilibrate rapidly<\/li>\n\n\n\n<li>solvent density may equilibrate relatively rapidly<\/li>\n\n\n\n<li>side-chain conformations may take longer<\/li>\n\n\n\n<li>large protein domains may require much longer sampling<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">Thus, there is no single universal equilibration time.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Convergence<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Convergence is a stronger concept.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A measured property can be considered reasonably converged when additional simulation does not substantially change its estimated value or distribution.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Examples include convergence of:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>conformational populations<\/li>\n\n\n\n<li>average distances<\/li>\n\n\n\n<li>free-energy profiles<\/li>\n\n\n\n<li>principal components<\/li>\n\n\n\n<li>contact probabilities<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">Convergence is often difficult to establish from a single trajectory.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Independent replicate simulations can provide important evidence about reproducibility and sampling.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Structural and Dynamic Analyses<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">Once the simulation has been evaluated sufficiently, the trajectory can be analyzed in greater detail.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The appropriate analysis depends on the biological question.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">RMSD<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>Root mean square deviation<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Measures overall structural deviation relative to a reference structure.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Useful for:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>monitoring global conformational changes<\/li>\n\n\n\n<li>comparing different trajectories<\/li>\n\n\n\n<li>identifying large structural transitions<\/li>\n<\/ul>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">RMSF<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>Root mean square fluctuation<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Measures how strongly individual atoms or residues fluctuate around their average positions.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Useful for identifying:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>flexible loops<\/li>\n\n\n\n<li>rigid structural regions<\/li>\n\n\n\n<li>ligand-induced changes in flexibility<\/li>\n<\/ul>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Radius of Gyration<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Measures the overall spatial compactness of a molecule.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">It can help identify:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>expansion<\/li>\n\n\n\n<li>compaction<\/li>\n\n\n\n<li>unfolding-related changes<\/li>\n<\/ul>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Distance and Angle Analysis<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Distances or angles between selected structural features can be monitored directly.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Examples include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>distance between two domains<\/li>\n\n\n\n<li>catalytic residue distance<\/li>\n\n\n\n<li>ligand-protein distance<\/li>\n\n\n\n<li>opening angle of a protein domain<\/li>\n\n\n\n<li>loop opening and closing<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">These measurements are often more biologically interpretable than global RMSD.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Hydrogen-Bond Analysis<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Hydrogen bonds can be monitored over time.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Questions may include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>How frequently is a hydrogen bond formed?<\/li>\n\n\n\n<li>Is an interaction stable throughout the trajectory?<\/li>\n\n\n\n<li>Does ligand binding introduce new hydrogen bonds?<\/li>\n<\/ul>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Contact Analysis<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Residue-residue or protein-ligand contacts can be monitored.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A contact map can reveal how interactions change across different conformational states.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Secondary Structure<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">The formation and loss of:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>alpha helices<\/li>\n\n\n\n<li>beta sheets<\/li>\n\n\n\n<li>turns<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">can be monitored throughout the trajectory.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Solvent Accessible Surface Area<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Solvent-accessible surface area (SASA) can be used to assess:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>protein exposure to solvent<\/li>\n\n\n\n<li>burial of hydrophobic regions<\/li>\n\n\n\n<li>accessibility of ligand-binding regions<\/li>\n<\/ul>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Clustering<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Trajectory structures can be grouped according to structural similarity.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Clustering can help identify representative conformational states rather than examining thousands of individual frames.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Essential Dynamics<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">Proteins often contain many atomic fluctuations, but only a smaller number of collective motions may dominate their behavior.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>Principal component analysis (PCA)<\/strong>, also called <strong>essential dynamics<\/strong> in the context of molecular simulation, can be used to identify these dominant motions.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">PCA can reveal patterns such as:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>domain opening and closing<\/li>\n\n\n\n<li>collective loop movements<\/li>\n\n\n\n<li>large-scale twisting<\/li>\n\n\n\n<li>hinge motion<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">These motions may be related to biological function.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Dynamic Correlation<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">Different regions of a protein may move together or in opposite directions.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Dynamic cross-correlation analysis can be used to examine relationships between these motions.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">It may help identify:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>mechanically coupled residues<\/li>\n\n\n\n<li>communication between distant domains<\/li>\n\n\n\n<li>residues associated with allosteric regulation<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">However, correlated motion alone does not prove direct causality or molecular communication.<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>Discussion: Correlation versus causation<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">If two residues move together, does one residue control the other? Additional structural and experimental evidence is generally required.<\/p>\n<\/blockquote>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">From Motion to Biological Mechanism<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">The objective of MD analysis should not simply be to calculate as many structural measurements as possible.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Instead, analyses should be selected according to a biological hypothesis.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">For example, suppose a protein is believed to become active through opening of an active-site loop.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A possible analysis strategy could be:<\/p>\n\n\n\n<ol class=\"wp-block-list\">\n<li>Measure the distance or angle describing loop opening.<\/li>\n\n\n\n<li>Determine the distribution of open and closed conformations.<\/li>\n\n\n\n<li>Identify residues whose interactions change during opening.<\/li>\n\n\n\n<li>Use PCA to determine whether loop opening represents a dominant collective motion.<\/li>\n\n\n\n<li>Use correlation analysis to identify residues whose motions are associated with the loop.<\/li>\n\n\n\n<li>Compare wild-type and mutant simulations.<\/li>\n\n\n\n<li>Test whether the computational prediction is consistent with experimental observations.<\/li>\n<\/ol>\n\n\n\n<p class=\"wp-block-paragraph\">In this way, MD becomes a tool for testing mechanistic hypotheses rather than simply generating molecular animations.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Limitations of Classical MD<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">Although molecular dynamics is powerful, several important limitations should always be considered.<\/p>\n\n\n\n<h2 class=\"wp-block-heading\">Force-Field Approximation<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">The accuracy of the simulation depends on the quality of the force field.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Force fields are approximations of molecular interactions and may not reproduce every physical property equally well.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Fixed Chemical Topology<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Conventional MD generally does not allow bonds to form or break.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Chemical reactions therefore require other approaches.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Electronic Polarization<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Many commonly used biomolecular force fields use fixed atomic partial charges.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Real electron distributions respond dynamically to their environment.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">This polarization is only approximated indirectly in standard fixed-charge force fields.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Sampling Limitations<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">The simulation may be too short to sample important conformational transitions.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A trajectory that appears stable may therefore represent only one region of a much larger conformational landscape.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Initial Structure Dependence<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Simulation results may depend on the starting conformation.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Independent replicas can help determine whether observations are reproducible.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h2 class=\"wp-block-heading\">Model Dependence<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">Results can be affected by choices such as:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>force field<\/li>\n\n\n\n<li>water model<\/li>\n\n\n\n<li>ion parameters<\/li>\n\n\n\n<li>protonation states<\/li>\n\n\n\n<li>simulation-box dimensions<\/li>\n\n\n\n<li>thermostat<\/li>\n\n\n\n<li>barostat<\/li>\n\n\n\n<li>simulation length<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\">These choices should therefore be documented and justified.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Key Conceptual Summary<\/h1>\n\n\n\n<p class=\"wp-block-paragraph\">A molecular dynamics simulation can be summarized conceptually as:<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">\\boxed{ \\text{Molecular structure} \\rightarrow \\text{Force field} \\rightarrow \\text{Potential energy} \\rightarrow \\text{Forces} \\rightarrow \\text{Atomic motion} \\rightarrow \\text{Trajectory} \\rightarrow \\text{Conformational ensemble} \\rightarrow \\text{Biological interpretation} }<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The force field describes how the molecular system stores potential energy.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The gradient of this energy determines the forces acting on atoms.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Newtonian mechanics determines how those forces change atomic positions and velocities.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Repeated calculations generate a molecular trajectory.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The trajectory can then be analyzed to understand the conformational states and motions of the molecule.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The major scientific challenge is not simply generating a trajectory, but determining whether the simulation has sampled the molecular behavior sufficiently to support a reliable biological interpretation.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Topics for Further Study<\/h1>\n\n\n\n<ul class=\"wp-block-list\">\n<li><a href=\"Classical mechanics\">Classical mechanics<\/a><\/li>\n\n\n\n<li><a href=\"Newtonian mechanics\">Newtonian mechanics<\/a><\/li>\n\n\n\n<li><a href=\"Force field\">Force field<\/a><\/li>\n\n\n\n<li><a href=\"CHARMM36\">CHARMM36<\/a><\/li>\n\n\n\n<li><a href=\"Potential energy\">Potential energy<\/a><\/li>\n\n\n\n<li><a href=\"Lennard-Jones potential\">Lennard-Jones potential<\/a><\/li>\n\n\n\n<li><a href=\"Coulomb potential\">Coulomb potential<\/a><\/li>\n\n\n\n<li><a href=\"Molecular topology\">Molecular topology<\/a><\/li>\n\n\n\n<li><a href=\"Periodic boundary conditions\">Periodic boundary conditions<\/a><\/li>\n\n\n\n<li><a href=\"Particle Mesh Ewald\">Particle Mesh Ewald<\/a><\/li>\n\n\n\n<li><a href=\"Numerical integration\">Numerical integration<\/a><\/li>\n\n\n\n<li><a href=\"LINCS\">LINCS<\/a><\/li>\n\n\n\n<li><a href=\"Thermostat\">Thermostat<\/a><\/li>\n\n\n\n<li><a href=\"Barostat\">Barostat<\/a><\/li>\n\n\n\n<li><a href=\"NVE ensemble\">NVE ensemble<\/a><\/li>\n\n\n\n<li><a href=\"NVT ensemble\">NVT ensemble<\/a><\/li>\n\n\n\n<li><a href=\"NPT ensemble\">NPT ensemble<\/a><\/li>\n\n\n\n<li><a href=\"Energy minimization\">Energy minimization<\/a><\/li>\n\n\n\n<li><a href=\"Equilibration\">Equilibration<\/a><\/li>\n\n\n\n<li><a href=\"Convergence\">Convergence<\/a><\/li>\n\n\n\n<li><a href=\"RMSD\">RMSD<\/a><\/li>\n\n\n\n<li><a href=\"RMSF\">RMSF<\/a><\/li>\n\n\n\n<li><a href=\"Principal component analysis\">Principal component analysis<\/a><\/li>\n\n\n\n<li><a href=\"Essential dynamics\">Essential dynamics<\/a><\/li>\n\n\n\n<li><a href=\"Dynamic cross-correlation\">Dynamic cross-correlation<\/a><\/li>\n\n\n\n<li><a href=\"Molecular ensemble\">Molecular ensemble<\/a><\/li>\n\n\n\n<li><a href=\"Enhanced sampling\">Enhanced sampling<\/a><\/li>\n<\/ul>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h1 class=\"wp-block-heading\">Resources<\/h1>\n\n\n\n<h2 class=\"wp-block-heading\">MD Tutorials<\/h2>\n\n\n\n<ul class=\"wp-block-list\">\n<li><a href=\"http:\/\/www.mdtutorials.com\/\" target=\"_blank\" rel=\"noreferrer noopener\">http:\/\/www.mdtutorials.com\/<\/a><\/li>\n\n\n\n<li>GROMACS tutorials developed by Justin A. Lemkul<\/li>\n\n\n\n<li>A useful practical introduction to preparing, running, and analyzing molecular dynamics simulations with GROMACS.<\/li>\n<\/ul>\n\n\n\n<h2 class=\"wp-block-heading\">GROMACS Documentation<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">The official GROMACS documentation is particularly useful when connecting these theoretical concepts to practical simulation parameters, including:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>integrators<\/li>\n\n\n\n<li>constraints<\/li>\n\n\n\n<li>thermostats<\/li>\n\n\n\n<li>barostats<\/li>\n\n\n\n<li>electrostatics<\/li>\n\n\n\n<li>cutoffs<\/li>\n\n\n\n<li>trajectory output<\/li>\n\n\n\n<li>analysis methods<\/li>\n<\/ul>\n\n\n\n<p class=\"wp-block-paragraph\"><\/p>\n\n\n\n<h2 class=\"wp-block-heading\">ChatGPT<\/h2>\n\n\n\n<p class=\"wp-block-paragraph\">This document was originally written by Seungpyo Hong and subsequently substantially edited with the assistance of ChatGPT.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Introduction Molecular dynamics (MD) simulation is a computational method for investigating the dynamic behavior of molecules, particularly biomolecules such as proteins, nucleic acids, lipids, and their complexes. In an MD simulation, a molecular system is represented as a collection of atoms whose interactions are described using a predefined mathematical model called a force field. The [&hellip;]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"_uag_custom_page_level_css":"","site-sidebar-layout":"default","site-content-layout":"","ast-site-content-layout":"default","site-content-style":"default","site-sidebar-style":"default","ast-global-header-display":"","ast-banner-title-visibility":"","ast-main-header-display":"","ast-hfb-above-header-display":"","ast-hfb-below-header-display":"","ast-hfb-mobile-header-display":"","site-post-title":"","ast-breadcrumbs-content":"","ast-featured-img":"","footer-sml-layout":"","ast-disable-related-posts":"","theme-transparent-header-meta":"","adv-header-id-meta":"","stick-header-meta":"","header-above-stick-meta":"","header-main-stick-meta":"","header-below-stick-meta":"","astra-migrate-meta-layouts":"set","ast-page-background-enabled":"default","ast-page-background-meta":{"desktop":{"background-color":"var(--ast-global-color-4)","background-image":"","background-repeat":"repeat","background-position":"center center","background-size":"auto","background-attachment":"scroll","background-type":"","background-media":"","overlay-type":"","overlay-color":"","overlay-opacity":"","overlay-gradient":""},"tablet":{"background-color":"","background-image":"","background-repeat":"repeat","background-position":"center center","background-size":"auto","background-attachment":"scroll","background-type":"","background-media":"","overlay-type":"","overlay-color":"","overlay-opacity":"","overlay-gradient":""},"mobile":{"background-color":"","background-image":"","background-repeat":"repeat","background-position":"center center","background-size":"auto","background-attachment":"scroll","background-type":"","background-media":"","overlay-type":"","overlay-color":"","overlay-opacity":"","overlay-gradient":""}},"ast-content-background-meta":{"desktop":{"background-color":"var(--ast-global-color-5)","background-image":"","background-repeat":"repeat","background-position":"center center","background-size":"auto","background-attachment":"scroll","background-type":"","background-media":"","overlay-type":"","overlay-color":"","overlay-opacity":"","overlay-gradient":""},"tablet":{"background-color":"var(--ast-global-color-5)","background-image":"","background-repeat":"repeat","background-position":"center center","background-size":"auto","background-attachment":"scroll","background-type":"","background-media":"","overlay-type":"","overlay-color":"","overlay-opacity":"","overlay-gradient":""},"mobile":{"background-color":"var(--ast-global-color-5)","background-image":"","background-repeat":"repeat","background-position":"center center","background-size":"auto","background-attachment":"scroll","background-type":"","background-media":"","overlay-type":"","overlay-color":"","overlay-opacity":"","overlay-gradient":""}},"footnotes":"","_members_access_role":[],"_members_access_error":""},"categories":[48,43],"tags":[],"class_list":["post-4853","post","type-post","status-publish","format-standard","hentry","category-computation-biochemistry","category-structural-bioinformatics-bdsl"],"yoast_head":"<!-- This site is optimized with the Yoast SEO plugin v28.4 - https:\/\/yoast.com\/product\/yoast-seo-wordpress\/ -->\n<title>Introduction to Molecular Dynamics Simulation - Biomedical Data Science Laboratory<\/title>\n<meta name=\"robots\" content=\"index, follow, max-snippet:-1, max-image-preview:large, max-video-preview:-1\" \/>\n<link rel=\"canonical\" href=\"https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/\" \/>\n<meta property=\"og:locale\" content=\"en_US\" \/>\n<meta property=\"og:type\" content=\"article\" \/>\n<meta property=\"og:title\" content=\"Introduction to Molecular Dynamics Simulation - Biomedical Data Science Laboratory\" \/>\n<meta property=\"og:description\" content=\"Introduction Molecular dynamics (MD) simulation is a computational method for investigating the dynamic behavior of molecules, particularly biomolecules such as proteins, nucleic acids, lipids, and their complexes. In an MD simulation, a molecular system is represented as a collection of atoms whose interactions are described using a predefined mathematical model called a force field. The [&hellip;]\" \/>\n<meta property=\"og:url\" content=\"https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/\" \/>\n<meta property=\"og:site_name\" content=\"Biomedical Data Science Laboratory\" \/>\n<meta property=\"article:published_time\" content=\"2026-08-11T03:27:43+00:00\" \/>\n<meta property=\"article:modified_time\" content=\"2026-08-11T05:36:15+00:00\" \/>\n<meta name=\"author\" content=\"sphong\" \/>\n<meta name=\"twitter:card\" content=\"summary_large_image\" \/>\n<meta name=\"twitter:label1\" content=\"Written by\" \/>\n\t<meta name=\"twitter:data1\" content=\"sphong\" \/>\n\t<meta name=\"twitter:label2\" content=\"Est. reading time\" \/>\n\t<meta name=\"twitter:data2\" content=\"17 minutes\" \/>\n<script type=\"application\/ld+json\" class=\"yoast-schema-graph\">{\"@context\":\"https:\\\/\\\/schema.org\",\"@graph\":[{\"@type\":\"Article\",\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/introduction-to-molecular-dynamics-simulation\\\/#article\",\"isPartOf\":{\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/introduction-to-molecular-dynamics-simulation\\\/\"},\"author\":{\"name\":\"sphong\",\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/#\\\/schema\\\/person\\\/8d22f775fcc218b1184ec0d890034e9b\"},\"headline\":\"Introduction to Molecular Dynamics Simulation\",\"datePublished\":\"2026-08-11T03:27:43+00:00\",\"dateModified\":\"2026-08-11T05:36:15+00:00\",\"mainEntityOfPage\":{\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/introduction-to-molecular-dynamics-simulation\\\/\"},\"wordCount\":4014,\"commentCount\":0,\"publisher\":{\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/#organization\"},\"articleSection\":[\"Computation Biochemistry\",\"Structural Bioinformatics (BDSL)\"],\"inLanguage\":\"en-US\",\"potentialAction\":[{\"@type\":\"CommentAction\",\"name\":\"Comment\",\"target\":[\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/introduction-to-molecular-dynamics-simulation\\\/#respond\"]}]},{\"@type\":\"WebPage\",\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/introduction-to-molecular-dynamics-simulation\\\/\",\"url\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/introduction-to-molecular-dynamics-simulation\\\/\",\"name\":\"Introduction to Molecular Dynamics Simulation - Biomedical Data Science Laboratory\",\"isPartOf\":{\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/#website\"},\"datePublished\":\"2026-08-11T03:27:43+00:00\",\"dateModified\":\"2026-08-11T05:36:15+00:00\",\"breadcrumb\":{\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/introduction-to-molecular-dynamics-simulation\\\/#breadcrumb\"},\"inLanguage\":\"en-US\",\"potentialAction\":[{\"@type\":\"ReadAction\",\"target\":[\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/introduction-to-molecular-dynamics-simulation\\\/\"]}]},{\"@type\":\"BreadcrumbList\",\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/introduction-to-molecular-dynamics-simulation\\\/#breadcrumb\",\"itemListElement\":[{\"@type\":\"ListItem\",\"position\":1,\"name\":\"Home\",\"item\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/\"},{\"@type\":\"ListItem\",\"position\":2,\"name\":\"Introduction to Molecular Dynamics Simulation\"}]},{\"@type\":\"WebSite\",\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/#website\",\"url\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/\",\"name\":\"Biomedical Data Science Laboratory\",\"description\":\"Home page for Biomedical Data Science Laboratory.\",\"publisher\":{\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/#organization\"},\"potentialAction\":[{\"@type\":\"SearchAction\",\"target\":{\"@type\":\"EntryPoint\",\"urlTemplate\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/?s={search_term_string}\"},\"query-input\":{\"@type\":\"PropertyValueSpecification\",\"valueRequired\":true,\"valueName\":\"search_term_string\"}}],\"inLanguage\":\"en-US\"},{\"@type\":\"Organization\",\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/#organization\",\"name\":\"Biomedical Data Science Laboratory\",\"url\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/\",\"logo\":{\"@type\":\"ImageObject\",\"inLanguage\":\"en-US\",\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/#\\\/schema\\\/logo\\\/image\\\/\",\"url\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/wp-content\\\/uploads\\\/2022\\\/11\\\/logo.png\",\"contentUrl\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/wp-content\\\/uploads\\\/2022\\\/11\\\/logo.png\",\"width\":1792,\"height\":638,\"caption\":\"Biomedical Data Science Laboratory\"},\"image\":{\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/#\\\/schema\\\/logo\\\/image\\\/\"}},{\"@type\":\"Person\",\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/#\\\/schema\\\/person\\\/8d22f775fcc218b1184ec0d890034e9b\",\"name\":\"sphong\",\"image\":{\"@type\":\"ImageObject\",\"inLanguage\":\"en-US\",\"@id\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/wp-content\\\/uploads\\\/2022\\\/04\\\/cropped-\uc99d\uba85\uc0ac\uc9c4-96x96.jpg\",\"url\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/wp-content\\\/uploads\\\/2022\\\/04\\\/cropped-\uc99d\uba85\uc0ac\uc9c4-96x96.jpg\",\"contentUrl\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/wp-content\\\/uploads\\\/2022\\\/04\\\/cropped-\uc99d\uba85\uc0ac\uc9c4-96x96.jpg\",\"caption\":\"sphong\"},\"description\":\"Assistant professor Department of Molecular Biology Jeonbuk National University\",\"sameAs\":[\"https:\\\/\\\/bdsl.jbnu.ac.kr\"],\"url\":\"https:\\\/\\\/bdsl.jbnu.ac.kr\\\/blog\\\/author\\\/sphong\\\/\"}]}<\/script>\n<!-- \/ Yoast SEO plugin. -->","yoast_head_json":{"title":"Introduction to Molecular Dynamics Simulation - Biomedical Data Science Laboratory","robots":{"index":"index","follow":"follow","max-snippet":"max-snippet:-1","max-image-preview":"max-image-preview:large","max-video-preview":"max-video-preview:-1"},"canonical":"https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/","og_locale":"en_US","og_type":"article","og_title":"Introduction to Molecular Dynamics Simulation - Biomedical Data Science Laboratory","og_description":"Introduction Molecular dynamics (MD) simulation is a computational method for investigating the dynamic behavior of molecules, particularly biomolecules such as proteins, nucleic acids, lipids, and their complexes. In an MD simulation, a molecular system is represented as a collection of atoms whose interactions are described using a predefined mathematical model called a force field. The [&hellip;]","og_url":"https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/","og_site_name":"Biomedical Data Science Laboratory","article_published_time":"2026-08-11T03:27:43+00:00","article_modified_time":"2026-08-11T05:36:15+00:00","author":"sphong","twitter_card":"summary_large_image","twitter_misc":{"Written by":"sphong","Est. reading time":"17 minutes"},"schema":{"@context":"https:\/\/schema.org","@graph":[{"@type":"Article","@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/#article","isPartOf":{"@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/"},"author":{"name":"sphong","@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/#\/schema\/person\/8d22f775fcc218b1184ec0d890034e9b"},"headline":"Introduction to Molecular Dynamics Simulation","datePublished":"2026-08-11T03:27:43+00:00","dateModified":"2026-08-11T05:36:15+00:00","mainEntityOfPage":{"@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/"},"wordCount":4014,"commentCount":0,"publisher":{"@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/#organization"},"articleSection":["Computation Biochemistry","Structural Bioinformatics (BDSL)"],"inLanguage":"en-US","potentialAction":[{"@type":"CommentAction","name":"Comment","target":["https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/#respond"]}]},{"@type":"WebPage","@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/","url":"https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/","name":"Introduction to Molecular Dynamics Simulation - Biomedical Data Science Laboratory","isPartOf":{"@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/#website"},"datePublished":"2026-08-11T03:27:43+00:00","dateModified":"2026-08-11T05:36:15+00:00","breadcrumb":{"@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/#breadcrumb"},"inLanguage":"en-US","potentialAction":[{"@type":"ReadAction","target":["https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/"]}]},{"@type":"BreadcrumbList","@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/introduction-to-molecular-dynamics-simulation\/#breadcrumb","itemListElement":[{"@type":"ListItem","position":1,"name":"Home","item":"https:\/\/bdsl.jbnu.ac.kr\/blog\/"},{"@type":"ListItem","position":2,"name":"Introduction to Molecular Dynamics Simulation"}]},{"@type":"WebSite","@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/#website","url":"https:\/\/bdsl.jbnu.ac.kr\/blog\/","name":"Biomedical Data Science Laboratory","description":"Home page for Biomedical Data Science Laboratory.","publisher":{"@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/#organization"},"potentialAction":[{"@type":"SearchAction","target":{"@type":"EntryPoint","urlTemplate":"https:\/\/bdsl.jbnu.ac.kr\/blog\/?s={search_term_string}"},"query-input":{"@type":"PropertyValueSpecification","valueRequired":true,"valueName":"search_term_string"}}],"inLanguage":"en-US"},{"@type":"Organization","@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/#organization","name":"Biomedical Data Science Laboratory","url":"https:\/\/bdsl.jbnu.ac.kr\/blog\/","logo":{"@type":"ImageObject","inLanguage":"en-US","@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/#\/schema\/logo\/image\/","url":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-content\/uploads\/2022\/11\/logo.png","contentUrl":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-content\/uploads\/2022\/11\/logo.png","width":1792,"height":638,"caption":"Biomedical Data Science Laboratory"},"image":{"@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/#\/schema\/logo\/image\/"}},{"@type":"Person","@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/#\/schema\/person\/8d22f775fcc218b1184ec0d890034e9b","name":"sphong","image":{"@type":"ImageObject","inLanguage":"en-US","@id":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-content\/uploads\/2022\/04\/cropped-\uc99d\uba85\uc0ac\uc9c4-96x96.jpg","url":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-content\/uploads\/2022\/04\/cropped-\uc99d\uba85\uc0ac\uc9c4-96x96.jpg","contentUrl":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-content\/uploads\/2022\/04\/cropped-\uc99d\uba85\uc0ac\uc9c4-96x96.jpg","caption":"sphong"},"description":"Assistant professor Department of Molecular Biology Jeonbuk National University","sameAs":["https:\/\/bdsl.jbnu.ac.kr"],"url":"https:\/\/bdsl.jbnu.ac.kr\/blog\/author\/sphong\/"}]}},"uagb_featured_image_src":{"full":false,"thumbnail":false,"medium":false,"medium_large":false,"large":false,"1536x1536":false,"2048x2048":false},"uagb_author_info":{"display_name":"sphong","author_link":"https:\/\/bdsl.jbnu.ac.kr\/blog\/author\/sphong\/"},"uagb_comment_info":1,"uagb_excerpt":"Introduction Molecular dynamics (MD) simulation is a computational method for investigating the dynamic behavior of molecules, particularly biomolecules such as proteins, nucleic acids, lipids, and their complexes. In an MD simulation, a molecular system is represented as a collection of atoms whose interactions are described using a predefined mathematical model called a force field. The&hellip;","_links":{"self":[{"href":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-json\/wp\/v2\/posts\/4853","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-json\/wp\/v2\/comments?post=4853"}],"version-history":[{"count":5,"href":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-json\/wp\/v2\/posts\/4853\/revisions"}],"predecessor-version":[{"id":4878,"href":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-json\/wp\/v2\/posts\/4853\/revisions\/4878"}],"wp:attachment":[{"href":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-json\/wp\/v2\/media?parent=4853"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-json\/wp\/v2\/categories?post=4853"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/bdsl.jbnu.ac.kr\/blog\/wp-json\/wp\/v2\/tags?post=4853"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}