Molecular Dynamics Simulations BCKDHA
For the molecular dynamics simulation two mutations should be chosen. We decided to take one mutation, which leads to a drastic change in the amino acid properties and may therefore be crucial. The other mutation changes methionine to lysine, which have similar physiochemical properties and structures. The mutations we chose are the following:
- M82L: (in PDB position 37)
- C264W: (in PDB position 219)
- Crystal water Extraction
repairPDB bckdha_mod.pdb -ssw 15 > crystalWater.pdb
- bckdha_mod.pdb is the manipulated pdb from last tast, where the missing residues in the 1U5B pdb file were handled.
The new version of repairPDB (found in Dropbox -> scrips) was used to execute all repairPDB commands of this task.
- When trying to extract the crystal water we received an empty file, indicating that there is no crystal water below 15 Å in our pdb file.
- Extract Protein
repairPDB bckdha_mod -nosol >wt.pdb
- first, extract sequence from pdb file:
repairPDB bckdha_mod.pdb -seq > wt.seq
- run SCWRL
scwrl -i wt.pdb -s wt.seq -o wt.pdb
- wt.seq is the extracted protein sequence in all lower case letters. In the case of mutants the mutated residue is given in an upper case letter.
- remove hydrogens
repairPDB wt.pdb -noh > wt.pdb
- Concatenate protein and crystal water
- As we didn't have crystal water, the concatenation step could be skipped.
Create important GROMACS files
We executed all commands as given in the task description here.
- Genion neutralized the charge by adding 56 Cl- and 57 Na+ ions in all three cases (wt, mutant 1 and mutant 2).
- The corresponding .top files were changed, s.t.
... SOL x-(56-57) NA 57 CL 56
- Crystal water overlap
- As we didn't have crystal water in our inital pdb file, there was no overlap. We still executed the following command to make sure there will be no problem.
repairPDB wt_solv.pdb -cleansol > wt_solv2.pdb
- The REMARK tag in the last line of the corresponding output files always showed that zero water molecules were removed.
- Therefore we didn't have our .top files this time.
Equilibration of system
To do the final molecular dynamics situation, we submitted the following script to the LRZ, which calls mdrun for the wildtype and the two mutated structures:
- $-o $HOME/LRZ/mdrun.WT_all.out -j y
- $-N MapleSyrup_ALL
- $-S /bin/bash
- $-M firstname.lastname@example.org
- $-l h_rt=32:00:00
- $-l march=x86_64
- $-pe mpi_32 32
mpirun -np 32 mdrun_mpi -v -deffnm $HOME/LRZ/wtMD
mpirun -np 32 mdrun_mpi -v -deffnm $HOME/LRZ/mut1MD
mpirun -np 32 mdrun_mpi -v -deffnm $HOME/LRZ/mut5MD