Difference between revisions of "Molecular Dynamics Simulations Analysis Hemochromatosis"

From Bioinformatikpedia
(internal RMSD)
(Pymol analysis of average and bfactor)
Line 254: Line 254:
   
 
=== Pymol analysis of average and bfactor ===
 
=== Pymol analysis of average and bfactor ===
  +
  +
  +
<figtable id="comparison">
  +
{| class="wikitable" style="float: right; margin: 0 0 0 1em; border: 2px solid darkgray;" cellpadding="0"
  +
! scope="row" align="left" |
  +
| align="right" | [[File:Hemo_MD_1a6zCProtAvg.png|thumb|300px]]
  +
| align="right" | [[File:Hemo_MD_R224WProtAvg.png|thumb|300px]]
  +
| align="right" | [[File:Hemo_MD_C282SProtAvg.png|thumb|300px]]
  +
|-
  +
|+ style="caption-side: bottom; text-align: left" |<font size=1>'''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>
  +
 
<br style="clear:both;">
 
<br style="clear:both;">
   

Revision as of 11:44, 29 July 2012

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 1: 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 1: 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>


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


RMSF for protein and C-alpha

<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>


<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>


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/>