https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Task_8_-_Molecular_Dynamics_Simulations&feed=atom&action=history
Task 8 - Molecular Dynamics Simulations - Revision history
2024-03-28T13:10:11Z
Revision history for this page on the wiki
MediaWiki 1.31.16
https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Task_8_-_Molecular_Dynamics_Simulations&diff=25527&oldid=prev
Boidolj: /* Intro */
2012-06-27T15:21:48Z
<p><span dir="auto"><span class="autocomment">Intro</span></span></p>
<table class="diff diff-contentalign-left" data-mw="interface">
<col class="diff-marker" />
<col class="diff-content" />
<col class="diff-marker" />
<col class="diff-content" />
<tr class="diff-title" lang="en">
<td colspan="2" style="background-color: #fff; color: #222; text-align: center;">← Older revision</td>
<td colspan="2" style="background-color: #fff; color: #222; text-align: center;">Revision as of 15:21, 27 June 2012</td>
</tr><tr>
<td colspan="2" class="diff-lineno">Line 3:</td>
<td colspan="2" class="diff-lineno">Line 3:</td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>In this section we will simulate the wildtype protein and two interesting mutants with MD, e.g. the gromacs package. For this we will use an automatic pipeline. As the final simulations will take a while, we will post the analysis part at a later point.</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>In this section we will simulate the wildtype protein and two interesting mutants with MD, e.g. the gromacs package. For this we will use an automatic pipeline. As the final simulations will take a while, we will post the analysis part at a later point.</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>The pipeline is available as a git repository. All the work needs to be done on the LRZ now.</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>The pipeline is available as a git repository. All the work needs to be done on the LRZ now.</div></td>
</tr>
<tr>
<td colspan="2" class="diff-empty"> </td>
<td class="diff-marker">+</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"></td>
</tr>
<tr>
<td colspan="2" class="diff-empty"> </td>
<td class="diff-marker">+</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>The slides of the task: [[File:MD talk.pdf]]</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>== LRZ ==</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>== LRZ ==</div></td>
</tr>
</table>
Boidolj
https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Task_8_-_Molecular_Dynamics_Simulations&diff=25443&oldid=prev
Offman: /* Submit Job */
2012-06-26T11:19:11Z
<p><span dir="auto"><span class="autocomment">Submit Job</span></span></p>
<table class="diff diff-contentalign-left" data-mw="interface">
<col class="diff-marker" />
<col class="diff-content" />
<col class="diff-marker" />
<col class="diff-content" />
<tr class="diff-title" lang="en">
<td colspan="2" style="background-color: #fff; color: #222; text-align: center;">← Older revision</td>
<td colspan="2" style="background-color: #fff; color: #222; text-align: center;">Revision as of 11:19, 26 June 2012</td>
</tr><tr>
<td colspan="2" class="diff-lineno">Line 49:</td>
<td colspan="2" class="diff-lineno">Line 49:</td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>If the test simulation fails due to a gromacs problem try to use only 16 cores and change that also for the commandline call of AGroS.</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>If the test simulation fails due to a gromacs problem try to use only 16 cores and change that also for the commandline call of AGroS.</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
</tr>
<tr>
<td class="diff-marker">−</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>In the real script you choose the standard cluster and instead of only <del class="diffchange diffchange-inline">4</del> hours (limit) you set something like 16-32 hours depending on the size of your protein.</div></td>
<td class="diff-marker">+</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>In the real script you choose the standard cluster and instead of only <ins class="diffchange diffchange-inline">2</ins> hours (limit) you set something like 16-32 hours depending on the size of your protein.</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>=Waiting=</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>=Waiting=</div></td>
</tr>
</table>
Offman
https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Task_8_-_Molecular_Dynamics_Simulations&diff=25442&oldid=prev
Offman: /* Prepare Job Scripts */
2012-06-26T11:18:54Z
<p><span dir="auto"><span class="autocomment">Prepare Job Scripts</span></span></p>
<table class="diff diff-contentalign-left" data-mw="interface">
<col class="diff-marker" />
<col class="diff-content" />
<col class="diff-marker" />
<col class="diff-content" />
<tr class="diff-title" lang="en">
<td colspan="2" style="background-color: #fff; color: #222; text-align: center;">← Older revision</td>
<td colspan="2" style="background-color: #fff; color: #222; text-align: center;">Revision as of 11:18, 26 June 2012</td>
</tr><tr>
<td colspan="2" class="diff-lineno">Line 31:</td>
<td colspan="2" class="diff-lineno">Line 31:</td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div> #SBATCH --ntasks=32</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div> #SBATCH --ntasks=32</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div> #SBATCH --mail-type=end</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div> #SBATCH --mail-type=end</div></td>
</tr>
<tr>
<td class="diff-marker">−</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div> #SBATCH --mail-user=<del class="diffchange diffchange-inline"><</del>offman@lrz.de<del class="diffchange diffchange-inline">></del></div></td>
<td class="diff-marker">+</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div> #SBATCH --mail-user=offman@lrz.de</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div> #SBATCH --export=NONE</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div> #SBATCH --export=NONE</div></td>
</tr>
<tr>
<td colspan="2" class="diff-empty"> </td>
<td class="diff-marker">+</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div> #SBATCH --time=02:00:00</div></td>
</tr>
<tr>
<td class="diff-marker">−</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div> #SBATCH --export=PATH="$HOME/test/AGroS:$PATH",PATH="$HOME/apps/bin/:$PATH"</div></td>
<td colspan="2" class="diff-empty"> </td>
</tr>
<tr>
<td class="diff-marker">−</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div> #SBATCH --time=04:00:00</div></td>
<td colspan="2" class="diff-empty"> </td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div> source /etc/profile.d/modules.sh</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div> source /etc/profile.d/modules.sh</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div> module load gromacs</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div> module load gromacs</div></td>
</tr>
<tr>
<td colspan="2" class="diff-empty"> </td>
<td class="diff-marker">+</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div> export PATH="$HOME/test/AGroS:$PATH"</div></td>
</tr>
<tr>
<td colspan="2" class="diff-empty"> </td>
<td class="diff-marker">+</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div> export PATH="$HOME/apps/bin/:$PATH"</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div> AGroS 1whz_new.pdb -dir /home/hpc/pr32fi/lu32xul/test -threads 32</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div> AGroS 1whz_new.pdb -dir /home/hpc/pr32fi/lu32xul/test -threads 32</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div></code></div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div></code></div></td>
</tr>
</table>
Offman
https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Task_8_-_Molecular_Dynamics_Simulations&diff=25440&oldid=prev
Offman at 08:18, 26 June 2012
2012-06-26T08:18:09Z
<p></p>
<table class="diff diff-contentalign-left" data-mw="interface">
<col class="diff-marker" />
<col class="diff-content" />
<col class="diff-marker" />
<col class="diff-content" />
<tr class="diff-title" lang="en">
<td colspan="2" style="background-color: #fff; color: #222; text-align: center;">← Older revision</td>
<td colspan="2" style="background-color: #fff; color: #222; text-align: center;">Revision as of 08:18, 26 June 2012</td>
</tr><tr>
<td colspan="2" class="diff-lineno">Line 41:</td>
<td colspan="2" class="diff-lineno">Line 41:</td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>In this script we do not use the standard cluster <code>--clusters=mpp1</code> but a test queue to get a quicker answer whether the simulation works at all.</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>In this script we do not use the standard cluster <code>--clusters=mpp1</code> but a test queue to get a quicker answer whether the simulation works at all.</div></td>
</tr>
<tr>
<td colspan="2" class="diff-empty"> </td>
<td class="diff-marker">+</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"></td>
</tr>
<tr>
<td colspan="2" class="diff-empty"> </td>
<td class="diff-marker">+</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>=Submit Job=</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>Submission is done using the following command <code>sbatch job.script</code></div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>Submission is done using the following command <code>sbatch job.script</code></div></td>
</tr>
<tr>
<td colspan="2" class="diff-lineno">Line 48:</td>
<td colspan="2" class="diff-lineno">Line 50:</td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>In the real script you choose the standard cluster and instead of only 4 hours (limit) you set something like 16-32 hours depending on the size of your protein.</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>In the real script you choose the standard cluster and instead of only 4 hours (limit) you set something like 16-32 hours depending on the size of your protein.</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
</tr>
<tr>
<td colspan="2" class="diff-empty"> </td>
<td class="diff-marker">+</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>=Waiting=</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>The state of the job and whether it really sits in the queue can be checked with the command <code>squeue -u <username> <queue></code> where the queue can either be <code>--clusters=mpp1</code> or <code>--partition=mpp1_inter</code>.</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>The state of the job and whether it really sits in the queue can be checked with the command <code>squeue -u <username> <queue></code> where the queue can either be <code>--clusters=mpp1</code> or <code>--partition=mpp1_inter</code>.</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
</tr>
</table>
Offman
https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Task_8_-_Molecular_Dynamics_Simulations&diff=25439&oldid=prev
Offman at 08:17, 26 June 2012
2012-06-26T08:17:09Z
<p></p>
<a href="https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Task_8_-_Molecular_Dynamics_Simulations&diff=25439&oldid=25432">Show changes</a>
Offman
https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Task_8_-_Molecular_Dynamics_Simulations&diff=25432&oldid=prev
Offman: /* Intro */
2012-06-26T07:30:58Z
<p><span dir="auto"><span class="autocomment">Intro</span></span></p>
<table class="diff diff-contentalign-left" data-mw="interface">
<col class="diff-marker" />
<col class="diff-content" />
<col class="diff-marker" />
<col class="diff-content" />
<tr class="diff-title" lang="en">
<td colspan="2" style="background-color: #fff; color: #222; text-align: center;">← Older revision</td>
<td colspan="2" style="background-color: #fff; color: #222; text-align: center;">Revision as of 07:30, 26 June 2012</td>
</tr><tr>
<td colspan="2" class="diff-lineno">Line 1:</td>
<td colspan="2" class="diff-lineno">Line 1:</td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>== Intro ==</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>== Intro ==</div></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
</tr>
<tr>
<td class="diff-marker">−</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>In this section we will simulate the wildtype protein and two interesting mutants with MD, e.g. the gromacs package. <del class="diffchange diffchange-inline">In</del> <del class="diffchange diffchange-inline">the following</del> we <del class="diffchange diffchange-inline">supply</del> <del class="diffchange diffchange-inline">you</del> <del class="diffchange diffchange-inline">with</del> <del class="diffchange diffchange-inline">the</del> pipeline<del class="diffchange diffchange-inline"> to simulate the proteins</del>. As the final simulations will take a while, we will post the analysis part at a later point.</div></td>
<td class="diff-marker">+</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #a3d3ff; vertical-align: top; white-space: pre-wrap;"><div>In this section we will simulate the wildtype protein and two interesting mutants with MD, e.g. the gromacs package. <ins class="diffchange diffchange-inline">For</ins> <ins class="diffchange diffchange-inline">this</ins> we <ins class="diffchange diffchange-inline">will</ins> <ins class="diffchange diffchange-inline">use</ins> <ins class="diffchange diffchange-inline">an</ins> <ins class="diffchange diffchange-inline">automatic</ins> pipeline. As the final simulations will take a while, we will post the analysis part at a later point.</div></td>
</tr>
<tr>
<td class="diff-marker">−</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"></td>
<td colspan="2" class="diff-empty"> </td>
</tr>
<tr>
<td class="diff-marker">−</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>A good tutorial can be found here:</div></td>
<td colspan="2" class="diff-empty"> </td>
</tr>
<tr>
<td class="diff-marker">−</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>http://md.chem.rug.nl/~mdcourse/index.html#introduction[http://md.chem.rug.nl/~mdcourse/index.html#introduction]</div></td>
<td colspan="2" class="diff-empty"> </td>
</tr>
<tr>
<td class="diff-marker">−</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"></td>
<td colspan="2" class="diff-empty"> </td>
</tr>
<tr>
<td class="diff-marker">−</td>
<td style="color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #ffe49c; vertical-align: top; white-space: pre-wrap;"><div>The related talk can be found here: [[File:Molecular_dynamics_talk.pdf]]</div></td>
<td colspan="2" class="diff-empty"> </td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"></td>
</tr>
<tr>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>== Preparation ==</div></td>
<td class="diff-marker"> </td>
<td style="background-color: #f8f9fa; color: #222; font-size: 88%; border-style: solid; border-width: 1px 1px 1px 4px; border-radius: 0.33em; border-color: #eaecf0; vertical-align: top; white-space: pre-wrap;"><div>== Preparation ==</div></td>
</tr>
</table>
Offman
https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Task_8_-_Molecular_Dynamics_Simulations&diff=25431&oldid=prev
Offman at 07:30, 26 June 2012
2012-06-26T07:30:09Z
<p></p>
<a href="https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Task_8_-_Molecular_Dynamics_Simulations&diff=25431&oldid=25430">Show changes</a>
Offman
https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Task_8_-_Molecular_Dynamics_Simulations&diff=25430&oldid=prev
Offman: Blanked the page
2012-06-26T07:27:16Z
<p>Blanked the page</p>
<a href="https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Task_8_-_Molecular_Dynamics_Simulations&diff=25430&oldid=25429">Show changes</a>
Offman
https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php?title=Task_8_-_Molecular_Dynamics_Simulations&diff=25429&oldid=prev
Offman: Created page with "== Task description == A detailed task description can be found here. == Intro and selection of two mutants == The description …"
2012-06-26T07:26:54Z
<p>Created page with "== Task description == A detailed task description can be found <a href="/wiki/index.php/Task_7_-_Structure-based_mutation_analysis" title="Task 7 - Structure-based mutation analysis"> here</a>. == Intro and selection of two mutants == The description …"</p>
<p><b>New page</b></p><div>== Task description ==<br />
<br />
A detailed task description can be found [[Task_7_-_Structure-based_mutation_analysis| here]].<br />
<br />
== Intro and selection of two mutants ==<br />
<br />
The description given here was applied for our wild type (1J8U). However, the steps done for our two mutants R408W and P281L are the same. <br />
<br />
We selected R408 because it is the most abundant mutation of PAH which is associated with phenylketonuria. It has a total frequency of 6.67% (see [http://www.pahdb.mcgill.ca/cgi-bin/pahdb/mutation_statistics-1.cgi]).<br />
<br />
The second mutation, P281L, was selected because it is the closed mutation to the binding site (HIS 285) and we hope that we are able to see some interesting things here which give us an explanation why this mutation is disease causing.<br />
<br />
== Preperation ==<br />
<br />
Before we can start to generate all necessary files for our MD simulation we have to prepare our structure first. We employed the following steps:<br />
<br />
'''1. Extracting the crystal water below 15 Å''' <br />
<br />
<code><br />
repairPDB 1J8U.pdb -ssw 15 > water_below_15.out<br />
</code><br />
<br />
We received the following output:<br />
<br />
<code><br />
HETATM 2559 O HOH A1008 2.996 9.738 16.094 1.00 14.91 O <br />
HETATM 2561 O HOH A1010 -3.667 12.788 8.538 1.00 14.80 O <br />
HETATM 2572 O HOH A1021 3.557 27.911 10.913 1.00 14.69 O <br />
HETATM 2879 O HOH A1328 -5.335 24.271 14.490 1.00 14.88 O <br />
TER<br />
</code><br />
<br />
<br />
'''2. Extract only the protein'''<br />
<br />
<code><br />
repairPDB 1J8U.pdb -nosol > 1J8U_nosol.pdb<br />
</code><br />
<br />
<br />
'''3. Extract the amino acid sequence and turn it into lower case letter to give it as a input for SCWRL'''<br />
<br />
<code><br />
repairPDB 1J8U.pdb -seq > 1J8U_seq.txt<br />
<br />
vim 1J8U_seq.txt<br />
:%s/.*/\L&/g<br />
</code><br />
<br />
<br />
'''4. Use SCWRL to complete the side chains of all amino acids'''<br />
<br />
<code><br />
/apps/scwrl4/Scwrl4 -i 1J8U_nosol.pdb -s 1J8U_seq.txt -o 1J8U_nosol_after_SCWRL.pdb | tee scwrl_1J8U_nosol_after_scwrl.out<br />
</code><br />
<br />
<br />
'''5. Remove the hydrogen atoms from the SCWRL output'''<br />
<br />
<code><br />
repairPDB 1J8U_nosol_after_SCWRL.pdb -noh > 1J8U_nosol_after_SCWRL_no_h.pdb<br />
</code><br />
<br />
<br />
'''6. concatenate the protein and the crystal water into one file'''<br />
<br />
We just added the output of step 1 to the end of the PDB file 1J8U_nosol_after_SCWRL_no_h.pdb and named this file ''1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.pdb''<br />
<br />
<br />
== Create important GROMACS files ==<br />
<br />
'''1. We have to create a TOP (topology) and a PAR (parameter) file with the following command:'''<br />
<br />
<code><br />
pdb2gmx -f 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.pdb -o 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.gro -p 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.top -water tip3p -ff amber03 -vsite hydrogens | tee pdb2gmx.out<br />
</code><br />
<br />
<br />
'''2. Next we create a box around the protein and fill it with water'''<br />
<br />
Create the box:<br />
<br />
<code><br />
editconf -f 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.gro -o 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_box.gro -bt dodecahedron -d 1.0 | tee editconf.out<br />
</code><br />
<br />
Fill it with water:<br />
<br />
<code><br />
genbox -cp 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_box.gro -cs spc216.gro -p 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.top -o 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_water.gro<br />
</code><br />
<br />
<br />
'''3. In the next step we have to add the ions. To do so we need to use the grompp command, which is always used to prepare a system for more complicated steps. This command needs an mdp file where we can define several additional parameters.'''<br />
<br />
Create parameter file addions.mdp:<br />
<br />
<code><br />
vim addions.mdp <br />
</code><br />
<br />
And fill it with the following content:<br />
<br />
<code><br />
integrator = steep<br />
emtol = 1.0<br />
nsteps = 500<br />
nstenergy = 1<br />
energygrps = System<br />
coulombtype = PME<br />
rcoulomb = 0.9<br />
rvdw = 0.9<br />
rlist = 0.9<br />
fourierspacing = 0.12<br />
pme_order = 4<br />
ewald_rtol = 1e-5<br />
</code><br />
<br />
Now, we execute it with this command:<br />
<br />
<code><br />
grompp -v -f addions.mdp -c 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_water.gro -p 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.top -o 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_water.tpr | tee grommp.out<br />
</code><br />
<br />
<br />
'''4. Now we add the solvent to the system and neutralize the charge that is part of the system due to charged amino acids. To do so we use genion:'''<br />
<br />
<code><br />
genion -s 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_water.tpr -o 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv.pdb -conc 0.1 -neutral -pname NA+ -nname CL- | tee genion.out<br />
</code><br />
<br />
During the call we were asked for a number we selected 13.<br />
<br />
<br />
'''5. Next we have to adjust the SOL value and add the number of NA+ and CL- in our 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.top'''<br />
<br />
We had 25 NA+ and 21 CL- ions generated from our last step. So we had to reduce the number in the SOL field to 10170 - 46 = '''10124''' the .top file had in the end the following values:<br />
<br />
<code><br />
[ molecules ]<br />
; Compound #mols<br />
Protein_chain_A 1<br />
SOL 4<br />
SOL 10124<br />
NA 25<br />
CL 21<br />
</code><br />
<br />
<br />
'''6. To make sure our initial crystal water does not overlap with the added water we use repairPDB again.'''<br />
<br />
<code><br />
repairPDB 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv.pdb -cleansol > 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2.pdb<br />
</code><br />
<br />
<br />
'''7. In the last line there is a REMARK tag with the number of removed water molecules. This has to be changed in the TOP file if larger than 0.'''<br />
<br />
<code><br />
REMARK RM 0<br />
</code><br />
<br />
Value is 0 so we did not have to change anything.<br />
<br />
<br />
'''8. We need to extract restraints from the structures. These restraints are useful to reduce the simulation time by disallowing very fast vibrations such as seen for hydrogen atoms. The command we use is called genrestr.'''<br />
<br />
<code><br />
genrestr -f 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2.pdb -o 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2.itp | tee genrestr.out<br />
</code><br />
<br />
In the command line menu we choose number 1 for protein.<br />
<br />
<br />
== Minimization solvent ==<br />
<br />
'''1. Creating a .mdp file for grompp to prepare the system for the solvent minimization'''<br />
<br />
We created an file with the name minimize_solvent.mdp by using vim:<br />
<br />
<code><br />
vim minimize_solvent.mdp<br />
</code><br />
<br />
The content of the file is as follows:<br />
<br />
<code><br />
define = -DPOSRES<br />
integrator = steep<br />
emtol = 1.0<br />
nsteps = 500<br />
nstenergy = 1<br />
energygrps = System<br />
coulombtype = PME<br />
rcoulomb = 0.9<br />
rvdw = 0.9<br />
rlist = 0.9<br />
fourierspacing = 0.12<br />
pme_order = 4<br />
ewald_rtol = 1e-5<br />
pbc = xyz<br />
</code><br />
<br />
<br />
'''2. execution of grompp to prepare the system for the MD run'''<br />
<br />
<code><br />
grompp -v -f minimize_solvent.mdp -c 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2.pdb -p 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.top -o 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2_min.tpr | tee grompp_minimize_solvent.out<br />
</code><br />
<br />
<br />
'''3. Executing mdrun to minimize the solvent'''<br />
<br />
<code><br />
mdrun -v -deffnm 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2_min -c 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2_min.pdb | tee mdrun_minimize_solvent.out<br />
</code><br />
<br />
<br />
== Minimization system ==<br />
<br />
In this part we minimize the solvent AND the side chains of the protein.<br />
<br />
'''1. Creating the position restraint file.'''<br />
<br />
<code><br />
genrestr -f 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2_min.pdb -o 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2_min.itp | tee genrestr_minimize_system.out<br />
</code><br />
<br />
We selected option 4 which says that we only put restraints on the backbone of the protein.<br />
<br />
<br />
'''2. Create an .mdp file for the grompp call'''<br />
<br />
We use the same .mdp from the solvent minimization step<br />
<br />
<br />
'''3. Calling grompp to prepare the tpr file for our mdrun step'''<br />
<br />
<code><br />
grompp -v -f minimize_solvent.mdp -c 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2_min.pdb -p 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.top -o 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2_min2.tpr | tee grompp_minimize_system.out<br />
</code><br />
<br />
<br />
'''4. Calling mdrun to minimize the system'''<br />
<br />
<code><br />
mdrun -v -deffnm 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2_min2 -c 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2_min2.pdb | tee mdrun_minimize_system.out <br />
</code><br />
<br />
<br />
== Equilibration of system ==<br />
<br />
=== Heating up the System with constant volume and particle number (NVT) ===<br />
<br />
Here we will heat up our system to a certain temperature. We will keep our particle number and volume constant, the pressure is changing. <br />
<br />
1. creating a .mdp configuration file for grompp<br />
<br />
We created a mdp file with the name nvt.mdp :<br />
<br />
<code><br />
vim nvt.mdp<br />
</code><br />
<br />
The content is as follows:<br />
<br />
<code><br />
define = -DPOSRES<br />
integrator = md<br />
dt = 0.005 <br />
nsteps = 10000<br />
nstxout = 0<br />
nstvout = 0<br />
nstfout = 0<br />
nstlog = 1000<br />
nstxtcout = 0<br />
nstenergy = 5<br />
energygrps = Protein Non-Protein<br />
nstcalcenergy = 5<br />
nstlist = 10<br />
ns-type = Grid<br />
pbc = xyz<br />
rlist = 0.9<br />
coulombtype = PME<br />
rcoulomb = 0.9<br />
rvdw = 0.9<br />
fourierspacing = 0.12<br />
pme_order = 4<br />
ewald_rtol = 1e-5 <br />
gen_vel = yes <br />
gen_temp = 200.0 <br />
gen_seed = 9999 <br />
constraints = all-bonds<br />
tcoupl = V-rescale <br />
tc-grps = Protein Non-Protein <br />
tau_t = 0.1 0.1 <br />
ref_t = 298 298 <br />
pcoupl = no<br />
</code><br />
<br />
<br />
'''2. Creating .tpr file for our MD run by employing grompp'''<br />
<br />
<code><br />
grompp -v -f nvt.mdp -c 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_solv2_min2.pdb -p 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.top -o 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_nvt.tpr | tee grompp_equi_system_nvt.out<br />
</code><br />
<br />
<br />
'''3. MD run'''<br />
<br />
<code><br />
mdrun -v -deffnm 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_nvt | tee mdrun_equi_system_nvt.out<br />
</code><br />
<br />
<br />
=== Heating up the System with constant pressure and particle number (NPT) ===<br />
<br />
Here we will heat up our system to a certain temperature. We will keep our particle number and pressure constant, the volume is changing. <br />
<br />
<br />
'''1. Creating an .mdp file for pressure coupling NPT'''<br />
<br />
We create a .mdp file with the following command:<br />
<br />
<code><br />
vim npt.mdp<br />
</code><br />
<br />
And the content:<br />
<br />
<code><br />
define = -DPOSRES<br />
integrator = md<br />
dt = 0.005 <br />
nsteps = 10000<br />
nstxout = 0<br />
nstvout = 0<br />
nstfout = 0<br />
nstlog = 1000<br />
nstxtcout = 0<br />
nstenergy = 5<br />
xtc_precision = 1000<br />
xtc-grps = System<br />
energygrps = Protein Non-Protein<br />
nstcalcenergy = 5<br />
nstlist = 5<br />
ns-type = Grid<br />
pbc = xyz<br />
rlist = 0.9<br />
coulombtype = PME<br />
rcoulomb = 0.9<br />
rvdw = 0.9<br />
fourierspacing = 0.12<br />
pme_order = 4<br />
ewald_rtol = 1e-5<br />
tcoupl = V-rescale <br />
tc-grps = Protein Non-Protein <br />
tau_t = 0.1 0.1 <br />
ref_t = 298 298 <br />
pcoupl = Berendsen<br />
Pcoupltype = Isotropic<br />
tau_p = 1.0<br />
compressibility = 4.5e-5<br />
ref_p = 1.0<br />
gen_vel = no<br />
constraints = all-bonds<br />
</code><br />
<br />
<br />
'''2. Creating .tpr file for our MD run by employing grompp'''<br />
<br />
<code><br />
grompp -v -f npt.mdp -c 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_nvt.gro -p 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.top -o 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_npt.tpr | tee grompp_equi_system_npt.out<br />
</code><br />
<br />
<br />
'''3. MD run'''<br />
<br />
<code><br />
mdrun -v -deffnm 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_npt | tee mdrun_equi_system_npt.out<br />
</code><br />
<br />
== Production run ==<br />
<br />
In this part we prepare the files for our final production run on the LRZ cluster.<br />
<br />
<br />
'''1. Creating a .mdp parameter file for grompp'''<br />
<br />
We created a file with the following command:<br />
<br />
<code><br />
vim md.mdp<br />
</code><br />
<br />
And this content: <br />
<br />
<code><br />
integrator = md<br />
tinit = 0<br />
dt = 0.005<br />
nsteps = 2000000<br />
nstxout = 50000<br />
nstvout = 50000<br />
nstfout = 0<br />
nstlog = 1000<br />
nstxtcout = 1000<br />
nstenergy = 1000<br />
energygrps = Protein Non-Protein <br />
nstcalcenergy = 5<br />
nstlist = 5<br />
ns-type = Grid<br />
pbc = xyz<br />
rlist = 0.9<br />
coulombtype = PME<br />
rcoulomb = 0.9<br />
rvdw = 0.9<br />
fourierspacing = 0.12<br />
pme_order = 4<br />
ewald_rtol = 1e-5<br />
tcoupl = V-rescale <br />
tc-grps = Protein Non-Protein <br />
tau_t = 0.1 0.1 <br />
ref_t = 298 298 <br />
pcoupl = Berendsen<br />
Pcoupltype = Isotropic<br />
tau_p = 2.0<br />
compressibility = 4.5e-5<br />
ref_p = 1.0<br />
gen_vel = no<br />
constraints = all-bonds<br />
constraint-algorithm = Lincs<br />
unconstrained-start = yes<br />
lincs-order = 4<br />
lincs-iter = 1<br />
lincs-warnangle = 30<br />
comm_mode = linear<br />
</code><br />
<br />
<br />
'''2. Creating a .tpr file for our MD run'''<br />
<br />
<code><br />
grompp -v -f md.mdp -c 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_npt.gro -p 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.top -o 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_md.tpr | tee grompp_prduction_run_md.out<br />
</code><br />
<br />
<br />
== LRZ ==<br />
<br />
=== Locally ===<br />
<br />
'''1. Copy files to LRZ'''<br />
<br />
We copied the following files to the LRZ:<br />
<br />
* md.mdp <br />
* 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water.top <br />
* 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_npt.gro<br />
* 1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_md.tpr <br />
<br />
By using the following command: <br />
<br />
<code><br />
scp -r . di69duj@lx64ia3.lrz.de:/home/cluster/pr58ni/di69duj<br />
</code><br />
<br />
<br />
=== On the LRZ ===<br />
<br />
'''1. Connecting to the LRZ'''<br />
<br />
<code><br />
ssh di69duj@lx64ia3.lrz.de<br />
</code><br />
<br />
<br />
'''2. Creating a working dictionary''' <br />
<br />
<code><br />
mkdir md_wt<br />
</code><br />
<br />
<br />
'''3. Copy all files to this working directory'''<br />
<br />
<code><br />
cp *.* md_wt<br />
</code><br />
<br />
And finally change to that directory:<br />
<br />
<code><br />
cd md_wt<br />
</code><br />
<br />
<br />
'''4. Creating a runfile with the name md_run.cmd'''<br />
<br />
Creating the file with vim:<br />
<br />
<code><br />
vim md_run.cmd<br />
</code><br />
<br />
With the following content:<br />
<br />
<code><br />
#!/bin/bash<br />
#$-o $HOME/md_wt/mdrun.WT_1.out -j y<br />
#$-N MD_WT<br />
#$-S /bin/bash<br />
#$-M erik.pfeiffenberger@gmail.com <br />
#$-l h_rt=32:00:00<br />
#$-l march=x86_64<br />
#$-pe mpi_32 32<br />
. /etc/profile<br />
cd $HOME/md_wt<br />
$HOME/md_wt/mpirun -np 32 mdrun_mpi -v -deffnm $HOME/md_wt/1J8U_nosol_after_SCWRL_no_h_merged_crystal_water_md<br />
</code><br />
<br />
<br />
'''5. Submitting job'''<br />
<br />
<code><br />
qsub md_run.cmd<br />
</code><br />
<br />
<br />
'''6. Waiting...'''<br />
<br />
While waiting, we recommend a cup of earl grey tea with one quarter milk and two tea spoons of sugar. Enjoy :)</div>
Offman