Difference between revisions of "Molecular Dynamics Simulations Analysis (PKU)"

From Bioinformatikpedia
(Short Introduction)
Line 3: Line 3:
   
 
==Initial Checks==
 
==Initial Checks==
  +
  +
All three simulations run for the desired 10 ns, the trajectories contain 2000 frames in 5 ps steps each. The wildtype simulation took significantly longer, since we used only 16 cores for the widtype, 32 for the mutants). Almost half of the calculation time, 44.2% in each run, is spent on calculating Coulomb interactions and the Lennard-Jones potential of the solvent molecules.
  +
  +
<figtable id="tab:modelling_scores">
  +
<caption>Key Values of the different SNP databases</caption>
  +
{| style="border-collapse: separate; border-spacing: 0; border-width: 1px; border-style: solid; padding-left:5px; padding-right:5px; border-color: #000; padding: 0"
  +
! style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| Mutation
  +
! style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| Sim. time
  +
! style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| Sim. speed
  +
! style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"| time to reach 1 s
  +
! style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"|
  +
! style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"|
  +
! style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0;"|
  +
! style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 0px 1px 0;"|
  +
|-
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0; font-style:italic;"| Wildtype
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"| 11:32 h
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"| 20.8 ns/day
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"| 131,621 years
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"|
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"|
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"|
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 0px 1px 0;"|
  +
|-
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0; font-style:italic;"| ALA322GLY
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"| 4:20 h
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"| 55.3 ns/day
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"| 49,543 years
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"|
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"|
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"|
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 0px 1px 0;"|
  +
|-
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; border-width: 0 1px 1px 0; font-style:italic;"| ARG408TRP
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"| 4:26 h
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"| 54.1 ns/day
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"| 50,685 years
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"|
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"|
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 1px 1px 0;"|
  +
| style="border-style: solid; padding-left:5px; padding-right:5px; text-align: right; border-width: 0 0px 1px 0;"|
  +
|-
  +
|}
  +
</figtable>

Revision as of 11:46, 11 July 2012

Short Introduction

We will analyze our completed molecular dynamics simulations, following the task description and the tutorial of the Utrecht University Molecular Modeling Practical. We have completed one run for the wildtype protein and for the mutations ALA322GLY and ARG408TRP, a second run of the wildtype is pending. The second run for the wildtype might be necessary as the trajectory of the wildtype differs significantly from both the mutants. The commands used to generate plots, images etc. can be found in our journal.

Initial Checks

All three simulations run for the desired 10 ns, the trajectories contain 2000 frames in 5 ps steps each. The wildtype simulation took significantly longer, since we used only 16 cores for the widtype, 32 for the mutants). Almost half of the calculation time, 44.2% in each run, is spent on calculating Coulomb interactions and the Lennard-Jones potential of the solvent molecules.

<figtable id="tab:modelling_scores"> Key Values of the different SNP databases

Mutation Sim. time Sim. speed time to reach 1 s
Wildtype 11:32 h 20.8 ns/day 131,621 years
ALA322GLY 4:20 h 55.3 ns/day 49,543 years
ARG408TRP 4:26 h 54.1 ns/day 50,685 years

</figtable>