GROMACS (GROningen MAchine for Chemical Simulations) is a widely-used open-source molecular dynamics simulation engine designed for biochemical molecules like proteins, lipids, and nucleic acids[1]. It is one of the fastest and most versatile MD engines available, making it essential for computational studies of neurodegenerative disease mechanisms.
GROMACS supports multiple force fields relevant to neurodegeneration research:
| Force Field | Applications | Notes |
|---|---|---|
| AMBER | Proteins, nucleic acids | Widely used |
| CHARMM | Lipids, proteins | Good for membranes |
| OPLS-AA | Small molecules | Solvation |
| GROMOS | General purpose | Less accurate |
| Method | Use Case |
|---|---|
| Metadynamics | Free energy landscapes |
| Umbrella Sampling | Binding affinities |
| Accelerated MD | Conformational changes |
| Steered MD | Pulling experiments |
GROMACS enables detailed studies of amyloid-beta and alpha-synuclein aggregation:
1. Protein Structure (PDB)
→ 2. pdb2gmx (force field)
→ 3. GROMACS (minimization/equilibration)
→ 4. Production MD
→ 5. Analysis (RMSD, contacts, etc.)
# From source
cmake ..
make -j$(nproc)
make install
# Conda
conda install -c conda-forge gromacs
# Energy minimization
gmx mdrun -v -deffnm em
# Production run
gmx mdrun -deffnm md_0_1
# Analysis
gmx rms -s md_0_1.tpr -f md_0_1.xtc -o rmsd.xvg
GROMACS excels at simulating membrane proteins relevant to neurodegeneration[2]:
Alchemical free energy methods in GROMACS enable precise binding affinity calculations[3]:
The Martini force field allows simulation of larger systems[2:1]:
Abraham MJ, et al. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX. 2015. ↩︎
Wang Y, et al. Molecular dynamics simulations of amyloid-beta and alpha-synuclein membrane interactions. Journal of Molecular Biology. 2021. ↩︎ ↩︎
Jorgensen WL, et al. Free energy calculations for drug discovery: Applications to neurodegeneration targets. Drug Discovery Today: Technologies. 2022. ↩︎