Difference between revisions of "MD WildeType"

From Bioinformatikpedia
(Create a movie)
(Create a movie)
Line 122: Line 122:
   
 
Todo: Filem erstellen und die filtered machen
 
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.944618
  +
|-
  +
|Tot-Drift
  +
|0.00834573
  +
|-
  +
|}
  +
  +
The plot with the temperature distribution of the system can be seen here:
  +
  +
[[Image:mut436_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 relativly 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.28165e+06
  +
|-
  +
|Error Estimation
  +
|100
  +
|-
  +
|RMSD
  +
|1080.9
  +
|-
  +
|Tot-Drift
  +
| -714.814
  +
|-
  +
|}
  +
  +
The plot with the potential energy distribution of the system can be seen here:
  +
  +
[[Image:mut436_md_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 engery of the system is between -1.282e+06 and -1.281e+06, which is a relativly 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.0519e+06
  +
|-
  +
|Error Estimation
  +
|100
  +
|-
  +
|RMSD
  +
|1322.68
  +
|-
  +
|Tot-Drift
  +
| -708.38
  +
|-
  +
|}
  +
  +
The plot with the total energy distribution of the system can be seen here:
  +
  +
[[Image:mut436_md_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.05e+06 and -1.051e+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.0066
  +
|-
  +
|Error Estimation
  +
|0.014
  +
|-
  +
|RMSD
  +
|71.218
  +
|-
  +
|Tot-Drift
  +
| -0.083422
  +
|-
  +
|}
  +
  +
The plot with the pressure distribution of the system can be seen here:
  +
  +
[[Image:mut436_md_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 1, but there a big outlier with 250 and -250 bar. So therefore we are not sure, if a protein can work with such a pressure.

Revision as of 09:25, 19 September 2011

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.

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

Time
Node (s) Real (s) %
22438.875 22438.875 1oo%
6h13:58

The complete simulation needs 6 hours and 13 minutes to finishing.

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:

Group 1 Group 2
System Water
Protein Ion
Backbone NA
MainChain+H CL
SideChain

Figure 1 shows the visualization with ngmx:

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

Average (in K) 297.94
Error Estimation 0.0029
RMSD 0.944618
Tot-Drift 0.00834573

The plot with the temperature distribution of the system can be seen here:

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 relativly cold for a protein to work, because the temperature in our bodies is about 36°C.

Potential

Average (in kJ/mol) -1.28165e+06
Error Estimation 100
RMSD 1080.9
Tot-Drift -714.814

The plot with the potential energy distribution of the system can be seen here:

Figure 3: Plot of the potential energy distribution of the MD system.

As can be seen on Figure 3, the potential engery of the system is between -1.282e+06 and -1.281e+06, which is a relativly 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

Average (in kJ/mol) -1.0519e+06
Error Estimation 100
RMSD 1322.68
Tot-Drift -708.38

The plot with the total energy distribution of the system can be seen here:

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.05e+06 and -1.051e+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

Average (in bar) 1.0066
Error Estimation 0.014
RMSD 71.218
Tot-Drift -0.083422

The plot with the pressure distribution of the system can be seen here:

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 1, but there a big outlier with 250 and -250 bar. So therefore we are not sure, if a protein can work with such a pressure.