Molecular Dynamics Simulations Analysis Hemochromatosis

From Bioinformatikpedia
Revision as of 09:17, 30 July 2012 by Joerdensv (talk | contribs) (Potential)

Hemochromatosis>>Task 10: Molecular dynamics simulations analysis


Short task description

Detailed description: Molecular dynamics simulations analysis


Protocol

A protocol with a description of the data acquisition and other scripts used for this task is available here.


Dummy

Calculation statistics

<figtable id="tab:simulation_stats"> Statistics of the MD simulations

Input Calc. time Calc. speed time to reach 1 s
Wildtype 13h31:15 17.750 ns/day 154350,8 years
C282S 13h35:05 17.667 ns/day 155075,9 years
R224W 13h35:02 17.668 ns/day 155067,1 years

</figtable>

GMXcheck revealed for all calculations that all 2001 frames were calculated, resulting in a 10ns model.

Energies

Pressure

<figtable id="comparison">

Hemo MD 1a6zC pressure.png
Hemo MD R224W pressure.png
Hemo MD C282S pressure.png
Table 1: different pressures of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutatuin at position 282 (C282S)

</figtable>


Temperature

<figtable id="comparison">

Hemo MD 1a6zC temperature.png
Hemo MD R224W temperature.png
Hemo MD C282S temperature.png
Table 2: different temperature energies of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>


Potential

<figtable id="comparison">

Hemo MD 1a6zC potential.png
Hemo MD R224W potential.png
Hemo MD C282S potential.png
Table 3: different potential energies of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>


Total energy

<figtable id="comparison">

Hemo MD 1a6zC totalEnergy.png
Hemo MD R224W totalEnergy.png
Hemo MD C282S totalEnergy.png
Table 1: different total energies of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>


Minimum distance between periodic boundary cells

<figtable id="comparison">

Hemo MD 1a6zC minPeriodicDist.png
Hemo MD R224W minPeriodicDist.png
Hemo MD C282S minPeriodicDist.png
Table 1: different total energies of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>


RMSF for protein and C-alpha

Protein based

<figtable id="comparison">

Hemo MD 1a6zC prot rmsf.png
Hemo MD R224W prot rmsf.png
Hemo MD C282S prot rmsf.png
Table 1: different RMS fluctuations (based on the whole protein) of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>

C-Alpha based

<figtable id="comparison">

Hemo MD 1a6zC ca rmsf.png
Hemo MD R224W ca rmsf.png
Hemo MD C282S ca rmsf.png
Table 1: different tRMS fluctuations (based on the the C-alpha of the backbone of the protein) of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>


Statistical values

Average StdDeviation Ratio significance
Wildtype 0.183621323529412 0.0605331780782099 50.0281092824371
R224W 0.232876838235294 0.0895836329605749 42.87282263833
C282S 0.268558823529412 0.0982244780281883 45.0924827835932
Table 1:

Pymol analysis of average and bfactor

<figtable id="comparison">

Hemo MD 1a6zCProtAvg.png
Hemo MD R224WProtAvg.png
Hemo MD C282SProtAvg.png
Table 1: different t gyrations (based on the the C-alphas of the backbone of the protein) of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>


Radius of gyration

<figtable id="comparison">

Hemo MD 1a6zC ca gyration.png
Hemo MD R224W ca gyration.png
Hemo MD C282S ca gyration.png
Table 1: different t gyrations (based on the the C-alphas of the backbone of the protein) of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>


<figtable id="comparison">

Hemo MD 1a6zC prot gyration.png
Hemo MD R224W prot gyration.png
Hemo MD C282S prot gyration.png
Table 1: different t gyrations (based on the the whole protein) of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>

solvent accessible surface area

<figtable id="comparison">

Hemo MD 1a6zC SAS.png
Hemo MD R224W SAS.png
Hemo MD C282S SAS.png
Table 1: display of the different solvent accessible surface values of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>

<figtable id="comparison">

Hemo MD 1a6zC resSAS.png
Hemo MD R224W resSAS.png
Hemo MD C282S resSAS.png
Table 1: display of the different solvent accessible surface values (normalized to per residue values) of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>


hydrogen-bonds between protein and protein / protein and water

Protein-Protein

<figtable id="comparison">

Hemo MD 1a6zC hBondsP2P.png
Hemo MD R224W hBondsP2P.png
Hemo MD C282S hBondsP2P xvg.png
Table 1: the number of hydrogen bonds inside the protein of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>

Protein-Water

<figtable id="comparison">

Hemo MD 1a6zC hBondsP2W.png
Hemo MD R224W hBondsP2W.png
Hemo MD C282S hBondsP2W.png
Table 1: the number of hydrogen bonds of the protein with water of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>



Ramachandran plots

<figtable id="comparison">

Hemo MD 1a6zC ramaPlot.png
Hemo MD R224W ramaPlot.png
Hemo MD C282S ramaPlot.png
Table 1: tramachandran Plots of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>



RMSD matrix

<figtable id="comparison">

Hemo MD 1a6zC PP rmsd-matrix.png
Hemo MD R224W PP rmsd-matrix.png
Hemo MD C282S PP rmsd-matrix.png
Table 1: rmsd matrices of the three calculated models over time (based on the whole protein). From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>


<figtable id="comparison">

Hemo MD 1a6zC MCB rmsd-matrix.png
Hemo MD R224W MCB rmsd-matrix.png
Hemo MD C282S MCB rmsd-matrix.png
Table 1: rmsd matrices of the three calculated models over time (based on the mainchain and C-betas). From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>



cluster analysis

whole protein based

<figtable id="comparison">

Hemo MD 1a6zC PP cluster-id-over-time.png
Hemo MD R224W PP cluster-id-over-time.png
Hemo MD C282S PP cluster-id-over-time.png
Hemo MD 1a6zC PPcluster-transitions.png
Hemo MD R224W PP cluster-transitions.png
Hemo MD C282S PP cluster-transitions.png
Hemo MD 1a6zC PP cluster-sizes.png
Hemo MD R224W PP cluster-sizes.png
Hemo MD C282S PP cluster-sizes.png
Table 1: rmsd matrices of the three calculated models over time (based on the whole protein). From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>



C-alpha based

<figtable id="comparison">

Hemo MD 1a6zC MCB cluster-id-over-time.png
Hemo MD R224W MCB cluster-id-over-time.png
Hemo MD C282S MCB cluster-id-over-time.png
Hemo MD 1a6zC MCB cluster-transitions.png
Hemo MD R224W MCB cluster-transitions.png
Hemo MD C282S MCB cluster-transitions.png
Hemo MD 1a6zC MCB cluster-sizes.png
Hemo MD R224W MCB cluster-sizes.png
Hemo MD C282S MCB cluster-sizes.png
Table 1: rmsd matrices of the three calculated models over time (based on the whole protein). From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>



internal RMSD

against starting structure

<figtable id="comparison">

Hemo MD 1a6zC rmsd-ca-vs-start.png
Hemo MD R224W rmsd-ca-vs-start.png
Hemo MD C282S rmsd-ca-vs-start.png
Hemo MD 1a6zC rmsd-allProt-atom-vs-start.png
Hemo MD R224W rmsd-allProt-atom-vs-start.png
Hemo MD C282S rmsd-allProt-atom-vs-start.png
Table 1: tramachandran Plots of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>




against average structure

<figtable id="comparison">

Hemo MD 1a6zC rmsd-ca-vs-average.png
Hemo MD R224W rmsd-ca-vs-average.png
Hemo MD C282S rmsd-ca-vs-average.png
Hemo MD 1a6zC rmsd-all-atom-vs-average.png
Hemo MD R224W rmsd-all-atom-vs-average.png
Hemo MD C282S rmsd-all-atom-vs-average.png
Table 1: tramachandran Plots of the three calculated models over time. From left to right: 1a6zC (wildtype), mutation at position 224 (R224W) and mutation at position 282 (C282S)

</figtable>




References

<references/>