Glucocerebrosidase Molecular Dynamics Simulation

Introduction

In this task we used Gromacs and performed a Molecular Dynamics Simulation with the wildtype structure 2NT0 of Glucocerebrosidase as well as with two mutated structures. The task description is available here.

Selection of the mutations

In addition to the wildtype structure we used two mutated structures for the Molecular Dynamics Simulation. Therefore we decided to take mutation number 7 and 10 as numbered in the previous tasks (6 and 7).

Mutation 7: Asn - Ser (Pos. 409/370)

We chose mutation number seven because it is the most common mutation in Gaucher Disease Type 1. In the sequence based analysis we had problems to clearly classify the mutation as damaging. We hope to get more insight into the damaging effect by molecular dynamics.

Mutation 10: Leu - Pro (Pos. 509/470)

We chose mutation number ten because in the sequence based analysis we predicted it as damaging. Although all our observations led us to this prediction the mutation is not listed in HGMD. So we want to analyse this mutation closer with molecular dynamics.

Workflow

We followed all the steps described here. Therefore we used one chain. Our molecule had no crystal water, so we had not to create a new file with the crystal water. After adding the number of Natrium and Chlorid ions we had to check, if there are clashes with the crystal water. As we had none in our structure there were no clashes and we got REMARK 0 as a result.

For the production run we copied the .mdp, .gro, .top and .tpr files to our LRZ directory of di69dub. There we created a "molecular dynamics" directory with one directory for each run (wt, mut7 and mut10).