Molecular Dynamics Simulations GLA
by Benjamin Drexler and Fabian Grandke
The aim of the task was to simulate the molecular dynamics of our protein. For that we used the tool GROMACS that has been used in the previous task, as well. Molecular dynamics are very useful to visualize the atomic movements within a protein. In our case we want to investigate, if there are differences between the molecular dynamics of the wildtype of the protein and the two mutations I117S and Q283P.
According to the task description we had to prepare some files and did a lot preprocessing steps, before we submit our final calculation to the LRZ.
At first we had to prepare our files. Thereto we completed missing sidechains and removed hydrogen atoms.
In this step we created a box around the protein and filled it with water. Then we added the particular Na and Cl ions, made sure that the charge is neutral ,and that there are no clashes between the initial crystal water and the box-fill-water.
In this step first the solvent was minimized and then the solvent and the side chains of the protein were minimized in combination.
Equilibration of the system
The equilibration process consists of two steps of heating up the system. First the particle number and volume have been kept constant, while the pressure changed. Secondly the particle number and pressure were kept constant, while the volume changed. In both steps the temperature was increased during many single steps.
The heated system then got finished for the final molecular dynamics simulation.
The previous prepared files were submitted to the LRZ Cluster and are calculated there, due to the fact that huge calculation powers are required.
The results of the molecular dynamics simulation can be found here.