Protocol BCKDHA MD
Contents
This Protocol show how to analyse the data from molecular dynamics simulation
A brief check of results
To verified that the simulations finished properly we first use the command
gmxcheck -f wt.xtc
Quality assurance
Energy calculations
To calculate the different energies we used the command:
g_energy -f wtMD.edr -o energy.xvg
After submitting this command we had to choose the energy which should calculated.
- Pressure: 13
- Temperature: 12
- Potential: 9
- Total Energy: 11
Minimum distance between periodic boundary cells
It is important to calculate the minimal distances to find out if there are direct interactions. Such interactions could appear if the distances are shorter than the cut off value of electrostatic interactions.
To calculate the minimum distance we used the command
g_mindist -f wtMD.xtc -s wtMD.tpr -od minimal-periodic-distance.xvg -pi
After submitting this command we chose group 1 to calculate the minimum distance for the whole protein.
Root mean square fluctuations
For calculating the RMSF of a protein each atom of this protein is compared with the calculated average stucture of the protein. By comparing them it is possible to find out how much it varies from its average position and so the flexibility of this region can be calculated. Regions with a high fluctuation are more flexible than regions with a low one.
To calculate the minimum distance we used the command
g_rmsf -f wtMD.xtc -s wtMD.tpr -o rmsf-per-residue.xvg -ox average.pdb -oq bfactors.pdb -res
After submitting this command we had to choose the group we want the RMSF to be calculated for:
- Protein: 1
- C-alpha: 3
Radius of gyration
The radius of gyration reflects how the structure changes during the simulation and how the shape changes during the time.
To calculate the radius of gyration we used the command
g_gyrate -f wtMD.xtc -s wtMD.tpr -o radius-of-gyration.xvg
After submitting this command we chose group 1 to calculate the radius of gyration for the whole protein
