Team:EPF-Lausanne/Modeling/RunSimulation

From 2009.igem.org

Revision as of 09:37, 23 July 2009 by Guittet (Talk | contribs)

How to run a complete simulation

Theory

The .pdb is generated from X-ray diffraction crystallography. This process requires a crystal of the protein, which occurs only at low temperature (~5°K). Taking this in consideration, we have to go through different steps to bring our protein to lab conditions.

Initial minimization

We start with a few minimization rounds, to reach a minmimum in the sense of potential energy.

Heating

When the protein is stable, we have to add heat (kinetic energy) to reach a higher temperature (~300°K). We have to take a special care not to add heat too fast, which would result in protein burst and explosion of our system.

First NPT

This is a relaxation step, with the number of atoms N, the pressure P and the temperature T all kept constant (NPT step). This is a kind of homogenization of the distribution of atoms inside our box.

NVT

This is also a relaxation step, with atom number, volume and temperature constant.

Second NPT

Now we perform another relaxation NPT to reach lab conditions.

Final NPT

This is the final NPT, which last much longer and gives us the ouput of the simulation.


.conf file

Initial minimization

This is quite easy to achieve, as NAMD has a specific keyword to run a simulation:

minimize #OF_MINIMIZ_STEPS

Heating

We tried different protocol, but this one is the best:

reassignTemp $temp_step ;# starting temp
reassignIncr $temp_step ;# temp increment
reassignHold 300 ;# final temp
reassignFreq 100 ;# num timesteps to temp incr.