mol_model_course_2023

Workshop: Analyzing MD simulations in Gromacs

Description

• We will outline basic approaches to analyze MD simulations trajectories.
• Duration: 90 minutes
• Objectives: be able to
  • 6 seminar
  • Open and analyze trajectories using MDAnalysis
  • Properly analyze equillibration, autocorellation and uncertainty.
  • Analyze protein RMSD
  • Analyze distances, angles, contacts, hydrogen bonds, etc.
  • Analyze radial distribution functions.
  • Analyze kintetic parameters (diffusion, etc.)
  • Analyze thermodynamic fluctuations
  • RMSD matrix calculation and cluster analysis using Gromacs
  • Principal component analysis using gromacs.

Jupyter notebook

• MD_analysis.ipynb

Required software and resources

• Access to a Jupyter notebook evironment with Python 3, MDanalysis, nglview libraries
• Access to Newton cluster with Gromacs installed or install Gromacs at your local workstation. See http://www.gromacs.org

Learning resources

• MD Analysis tutorials

Assignments

From previous assignments you should have a simulation for a protein of your choice set up in Gromacs. You may need to adjust the simulation parameters if needed (trajectory write parameters, duration, temperature etc). Following analyses of this system will be needed as an assignment:

Seminar 6

1. Plot RMSD matrix for MD frames; perform cluster analysis for mainchain+H atoms; visualize clusters and dynamics of systems transitions between them
2. Perform PCA: plot atom covariance matrix, eigenvalues plot, calculate fraction of explained variability by first 6 vectors, make vector implementations and make representing first eigenvectors snapshots.

Troubleshooting

• Consult with the seminar protocol/recording
• Ask questions in Slack