Workshop: Analyzing MD simulations in Gromacs
Description
- We will outline basic approaches to analyze MD simulations trajectories.
- Duration: 90 minutes
- Objectives: be able to
- 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.)
- RMSD matrix calculation and cluster analysis using Gromacs
- Principal component analysis using gromacs.
Jupyter notebook
Required software and resources
- Access to a Jupyter notebook evironment with Python 3, MDanalysis, nglview and Gromacs
Learning resources
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:
- Calculate system density, perform proper etimate of statistical uncertainty.
- Calcuate RMSD of the protein with time. Make conclusions if the system has reached local equilibrium state.
- Calcuate distance between N- and C-ends of protein. Make conclusions if the system has reached local equilibrium state.
- Calculate radial distribution function between water molecules.
- Plot RMSD matrix for MD frames; perform cluster analysis for mainchain+H atoms; visualize clusters and dynamics of systems transitions between them.
- Perform PCA: plot atom covariance matrix, eigenvalues plot, make vector implementations and make representing first eigenvectors snapshots.
(*) Calcualte the number of hydrogen bonds within your protein and with water. Estimate average number (with uncertainty). Plot variation with time.
(*) Estimate water diffusion coefficients from mean square displacement.
Troubleshooting
- Consult with the seminar protocol/recording
- Ask questions in Telegram