Difference between revisions of "MD WildeType"

From Bioinformatikpedia
(energy calculations for pressure, temperature, potential and total energy)
(Blanked the page)
 
Line 1: Line 1:
=== check the trajectory ===
 
 
We checked the trajectory with following command:
 
 
gmxcheck -f 2GJX_A_md.xtc
 
 
With the command we got following results:
 
 
Reading frame 0 time 0.000
 
# Atoms 96543
 
Precision 0.001 (nm)
 
Last frame 2000 time 10000.000
 
 
Furthermore, we got some detailed results about the different items during the simulation.
 
 
{| border="1" style="text-align:center; border-spacing:0;"
 
|Item
 
|#frames
 
|Timestep (ps)
 
|-
 
|Step
 
|2001
 
|5
 
|-
 
|Time
 
|2001
 
|5
 
|-
 
|Lambda
 
|0
 
| -
 
|-
 
|Coords
 
|2001
 
|5
 
|-
 
|Velocities
 
|0
 
| -
 
|-
 
|Forces
 
|0
 
| -
 
|-
 
|Box
 
|2001
 
|5
 
|-
 
|}
 
 
The simulation finished on node 0 Thu Sep 15 23:45:08 2011
 
 
{| border="1" style="text-align:center; border-spacing:0;"
 
|colspan="3" | Time
 
|-
 
|Node (s)
 
|Real (s)
 
|%
 
|-
 
|22438.875
 
|22438.875
 
|1oo%
 
|-
 
|colspan="3" | 6h13:58
 
|-
 
|}
 
 
The complete simulation needs 6 hours and 13 minutes to finishing.
 
{| border="1" style="text-align:center; border-spacing:0;"
 
|colspan="4" | Performance
 
|-
 
|Mnbf/s
 
|GFlops
 
|ns/day
 
|hour/ns
 
|-
 
|1271.745
 
|93.383
 
|38.505
 
|0.623
 
|-
 
|}
 
 
As you can see in the table above, it takes about half an hour to simulate 1ns of the system. So therefore, it would be possible to simulate about 40ns in one complete day calculation time.
 
 
=== Visualize in pymol ===
 
 
First of all, we visualized the simulation with with ngmx, because it draws bonds based on the topology file. ngmx gave the user the possibility to choose different parameters. Therefore, we decided to visualize the system with following parameters:
 
 
{| border="1" style="text-align:center; border-spacing:0;"
 
|Group 1
 
|Group 2
 
|-
 
|System
 
|Water
 
|-
 
|Protein
 
|Ion
 
|-
 
|Backbone
 
|NA
 
|-
 
|MainChain+H
 
|CL
 
|-
 
|SideChain
 
|
 
|-
 
|}
 
 
Figure 1 shows the visualization with ngmx:
 
 
[[Image:ngmx_2GJX.png|thumb|center|Figure 1: Visualisation of the MD simulation for the wildtype with ngmx]]
 
 
=== Create a movie ===
 
 
Next, we want to visualize the protein with pymol. Therefore, we extracted 1000 frames from the trajectory, leaving out the water and jump over the boundaries to make continouse trajectories. Therefore, we used following command:
 
 
trjconv -s fole.tpr -f file.xtc -o output_file.pdb -pbc nojump -dt 10
 
 
The program asks for the a group as output. We only want to see the protein, therefore we decided to choose group 1.
 
 
Todo: Filem erstellen und die filtered machen
 
 
 
=== energy calculations for pressure, temperature, potential and total energy ===
 
 
==== Temperature ====
 
 
{| border="1" style="text-align:center; border-spacing:0;"
 
|Average (in K)
 
|297.94
 
|-
 
|Error Estimation
 
|0.0029
 
|-
 
|RMSD
 
|0.942857
 
|-
 
|Tot-Drift
 
|0.0066475
 
|-
 
|}
 
 
The plot with the temperature distribution of the system can be seen here:
 
 
[[Image:wt_md_temperatur.png|thumb|center|Figure 2: Plot of the temperature distribution of the MD system.]]
 
 
As you can see on Figure 2, most of the time the system has a temperature about 298K. The maximal difference between this average temperature and the minimum/maxmimum temperature is only about 4 K, which is not that high. But we have to keep in mind, that only some degree difference can destroy the function of a protein. 298 K is about 25°C, which is relatively cold for a protein to work, because the temperature in our bodies is about 36°C.
 
 
==== Potential ====
 
 
{| border="1" style="text-align:center; border-spacing:0;"
 
|Average (in kJ/mol)
 
| -1.2815e+06
 
|-
 
|Error Estimation
 
|100
 
|-
 
|RMSD
 
|1078.55
 
|-
 
|Tot-Drift
 
| -661.902
 
|-
 
|}
 
 
The plot with the potential energy distribution of the system can be seen here:
 
 
[[Image:wt_potential.png|thumb|center|Figure 3: Plot of the potential energy distribution of the MD system.]]
 
 
As can be seen on Figure 3, the potential energy of the system is between -1.285e+06 and -1.278e+06, which is a relatively low energy. Therefore this means that the protein is stable. So we can suggest, that the protein with such a low energy is able to function and is stable and therefore, our simulation could be true. Otherwise, if the energy of the simulated system is too high, we can not trust the results, because the protein is too instable to work.
 
 
==== Total energy ====
 
 
{| border="1" style="text-align:center; border-spacing:0;"
 
|Average (in kJ/mol)
 
| -1.05177e+06
 
|-
 
|Error Estimation
 
|100
 
|-
 
|RMSD
 
|1321.31
 
|-
 
|Tot-Drift
 
| -656.777
 
|-
 
|}
 
 
The plot with the total energy distribution of the system can be seen here:
 
 
[[Image:wt_total.png|thumb|center|Figure 4: Plot of the total energy distribution of the MD system.]]
 
 
As we can see on Figure 4 above, the total energy of the protein is a little bit higher than the potential energy of the protein. In this case, the energy is between -1.055e+06 and -1.048e+06. But these values are already in a range, where we can suggest that the energy of the protein is low enough so that this one can work.
 
 
==== Pressure ====
 
 
{| border="1" style="text-align:center; border-spacing:0;"
 
|Average (in bar)
 
|1.00711
 
|-
 
|Error Estimation
 
|0.0087
 
|-
 
|RMSD
 
|71.2473
 
|-
 
|Tot-Drift
 
| -0.0454746
 
|-
 
|}
 
 
The plot with the pressure distribution of the system can be seen here:
 
 
[[Image:wt_pressure.png|thumb|center|Figure 5: Plot of the pressure distribution of the MD system.]]
 
 
As you can see on Figure 5, the pressure in the system is most of the time about 0, but there a big outlier with 250 and about -200 bar. So therefore we are not sure, if a protein can work with such a pressure.
 

Latest revision as of 09:37, 19 September 2011