🧬 Molecular Dynamics Simulations — Theoretical Summary

Source:

🎯 Why Molecular Dynamics (MD) simulations exist

In experiments, proteins exist as ensembles of many conformations in solution. But when we look at a structure on a screen (e.g., PDB), we see only one static conformation.

MD simulations aim to:

✅ Mimic how atoms move in real life ✅ Explore all possible conformations over time ✅ Understand thermodynamic averages ✅ Reveal motions that are experimentally hard to observe

You can think of MD as a “computational microscope” — it lets us observe atomic motion with time resolution that experiments often cannot achieve.

🎲 Monte Carlo vs Molecular Dynamics

Two major simulation approaches exist:

🎲 Monte Carlo (MC)

• Samples conformations randomly
• Gives statistical distribution of states
• No time information

👉 Example idea: Like roulette statistics — you only know how often something happens.

⏱ Molecular Dynamics (MD)

• Simulates real motion over time
• Conformations are connected through physical trajectories

👉 Therefore MD can reveal mechanisms, not just probabilities.

⛰ Potential Energy Surface (Energy Landscape)

A key theoretical concept.

Proteins move on an energy landscape.

🔽 Energy minima

• Low-energy conformations
• Most stable → most likely in nature

✔ Yes — in thermodynamics, systems tend toward minimum free energy states (but not always global minimum due to kinetic barriers).

🌍 Global minimum

• Absolute lowest energy state
• Often considered the native fold

🪜 Local minima

• Metastable states
• Protein can get “trapped” there

🪑 Saddle Point (Important!)

A saddle point is:

• Maximum in one direction
• Minimum in another direction

Meaning:

➡ It represents a transition pathway between two conformations. ➡ It is related to activation barriers in conformational change.

So MD helps identify:

• Which conformational change route is easiest
• Which barrier must be crossed

This is fundamental for:

• Folding
• Channel gating
• Ligand binding pathways

⚙️ How MD simulation actually works

Basic algorithm:

1. Start with structure
2. Let atoms move for a tiny time step (~femtoseconds)
3. Stop
4. Calculate all forces
5. Update velocities and positions
6. Repeat

Thus:

👉 Motion becomes a time-resolved trajectory.

This is based on:

• Newton’s laws
• Potential energy function

🌊 Do atoms need to be bonded to interact?

❗ No.

Atoms only need to be close enough.

This means MD includes:

• Protein-protein interactions
• Protein-solvent interactions
• Lipid interactions
• Ion interactions

So yes — protein solvent effects are explicitly modeled.

This is extremely important because:

• Water stabilizes structure
• Solvation influences folding
• Channels allow transport

🧪 Why can ions not pass membranes freely?

Because lipid bilayers are:

• Hydrophobic
• Energetically unfavorable for charged species

Thus:

✔ Ions must pass via channel proteins

Experimentally:

• You can detect current
• But cannot easily track atomic path

To locate ions structurally:

• You would need to freeze motion
• Then crystallize

Very difficult → MD provides dynamic insight.

💊 Drug binding theory in MD

Docking vs MD difference:

Docking

• Often treats receptor as rigid
• Ligand flexible
• Cannot capture induced fit well

MD

• Protein and ligand both flexible
• Can observe binding pathways

🧭 Two strategies for ligand binding simulations

Strategy 1 — Start ligand inside

• Run MD → observe how it leaves
• Reverse trajectory → infer entry path

Strategy 2 — Start ligand outside

• Apply pulling force
• Guide ligand into active site

These approaches help understand:

• Binding kinetics
• Entry channels
• Energy barriers

⏳ Binding strength interpretation

If ligand:

• Leaves quickly → weak binding
• Stays long → strong binding

Residence time correlates with:

👉 Binding affinity

This is crucial in drug design.

🔁 Allosteric regulation

MD can reveal:

• Which ligand binds first
• How binding changes shape elsewhere
• Cooperative mechanisms

This is difficult with static docking.

🔓 Protein functional mechanisms

Example:

Active → inactive transition.

If ligand removed:

• Protein may relax to inactive form
• MD shows structural pathway

Thus MD explains:

✔ Signal transduction ✔ Channel gating ✔ Receptor activation

🧬 Protein folding theory

Huge challenge.

Reasons:

• Unknown unfolded starting structure
• Folding occurs over microseconds–seconds
• Requires enormous computing power

Historically:

• Distributed computing (Folding@home)
• Now GPU computing accelerates simulations

💡 Why static structures are misleading

A crystal/NMR structure is:

❗ A snapshot of a dynamic system

Like photographing a running horse mid-air and concluding it can fly.

Also:

• Structure has no temperature
• No vibrational motion
• Effectively resembles 0 K condition

Thus simulations must:

🔥 Heat system → equilibration phase

• Adjust velocities
• Stabilize pressure and temperature

📊 Production phase

• Collect meaningful data

🧮 Force fields — core theoretical idea

Force field =

👉 Large lookup table of parameters:

• Bond lengths
• Angles
• Charges
• van der Waals terms

Speeds up simulation because:

• Bonded interactions pre-parameterized
• Only non-bonded interactions calculated dynamically

⚠️ Limitation

Force fields are not universal.

Each has:

• Specific strengths
• Specific failures

Therefore:

✔ Always check what force field cannot do before using it.

🐢 Time-scale problem in MD

Different motions:

• Bond vibrations → very fast
• Side chain rotations → medium
• Loop movements → very slow

Thus:

To observe large conformational change:

👉 Need long simulation time.

💊 Drug “on-off” behavior — why dynamic binding can be better

Important theoretical pharmacology insight:

A drug that:

• Binds → unbinds → rebinds

may be better than one that binds permanently.

Why?

Because:

• Real biological channels open and close
• Dynamic drugs mimic physiology
• Permanent blockers may cause toxicity

MD helps evaluate such behavior.

⭐ Key Takeaway

Molecular dynamics is fundamentally about:

Understanding how biomolecules move on an energy landscape over time.

It connects:

• Thermodynamics
• Kinetics
• Structure
• Function
• Drug action

into one theoretical framework.