Task 10 - Molecular Dynamics Simulations 2011

From Bioinformatikpedia

By now you should all have finished your MD simulations. Two files that were produced are the most important ones, the XTC file (trajectory) and the EDR file (energy).

As there is a very good tutorial already available you will focus on following the instructions found here: http://md.chem.rug.nl/~mdcourse/analysis1.html[1] Mind, there is a second page of the analysis tutorial.

In this tutorial you will first find steps to assure quality of the simulations itself and secondly analysis of the molecules. In the second part the tutorial addresses a specific protein - skip the specific analysis and questions.

The steps should be applied for WT and the mutants simulated. The results should be compared. Most XVG files used for the graphs contain the raw data and could be further analysed (stdev, average...). Whenever an ECHO is preceding the command this means that these are the selections chosen in the interactive part of the command.

The following steps are done:

  • check the trajectory
  • Visualize in pymol only: the intrafit command is correct this way: intra_fit [filename without filetype ending]
  • create a movie and skip the g_filter step
  • energy calculations for pressure, temperature, potential and total energy
  • minimum distance between periodic boundary cells
  • RMSF for protein and C-alpha
  • Pymol analysis of average and bfactor
  • Radius of gyration
  • solvent accesible surface area
  • hydrogen-bonds between protein and protein / protein and water
  • salt bridges: is a mutant in the list of files or a close residue that could be affected? any big differences between WT and mutants
  • Ramachandran plot: compare it to the standard ramachandran plots. also compare WT and mutations
  • RMSD matrix
  • cluster analysis
  • internal RMSD