<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Zoe Hall skrev 2011-01-11 16.29:
<blockquote
cite="mid:C1F6430AA5E06845A7E4CCB00CBCEDF64C337DAFB6@EXMBX05.ad.oak.ox.ac.uk"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=ISO-8859-1">
<meta name="Generator" content="Microsoft Word 12 (filtered
medium)">
<style>
<!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
span.EmailStyle17
        {mso-style-type:personal;
        font-family:"Calibri","sans-serif";
        color:windowtext;}
span.EmailStyle18
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page Section1
        {size:612.0pt 792.0pt;
        margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.Section1
        {page:Section1;}
-->
</style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div class="Section1">
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Hi gmx-users,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I am trying to set up a simulation of a
large protein <i>in
vacuo</i> using the OPLS-AA forcefield, with conditions
based on Patriksson et.
al (Biochemistry 2007, 46 p933). Basically after energy
minimisation, the
protein is run for 10ps in vacuum with constant temperature at
300K. This is
followed by the full production run of 10ns with temperature
and pressure
coupling turned off, H bonds are constrained using SHAKE with
tolerance 0.001.
For vacuum conditions, the periodic boundary conditions are
turned off, and no
cut-offs are used whatsoever. <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I am not sure, however, what is the
appropriate choice for
“coulombtype” and other electrostatic parameters for vacuum
simulations. I have
set “epsilon_r” to 2, as this seems to be an accepted value
for proteins but if
anyone has comments on this, I would appreciate it.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks very much,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Zoe <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Zoe Hall <o:p></o:p></p>
<p class="MsoNormal">Department of Chemistry<o:p></o:p></p>
<p class="MsoNormal">Oxford University<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><u><span style="color: rgb(77, 84, 219);"><a class="moz-txt-link-abbreviated" href="mailto:zoe.hall@chem.ox.ac.uk">zoe.hall@chem.ox.ac.uk</a><o:p></o:p></span></u></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</blockquote>
I should point out that there is probably an error in the methods
section. If I recall correctly we used a 0.5 fs timestep. We
couldn't get good conservation of total energy with 1 fs steps.
Remember to use double precision, or you will again loose
conservation of total erergy. Remember to monitor the total energy,
because if it's drifting you need to have a more conservative setup,
e.g., use a smaller timestep.<br>
<br>
Take a look at our "follow up" study (Marklund et. al., PCCP (2009)
11 p. 8069–8078) where we used LINCS and managed to conserve the
total energy with a 1 fs time step.<br>
<br>
Good luck!<br>
<br>
<pre class="moz-signature" cols="72">--
-----------------------------------------------
Erik Marklund, PhD student
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596, 75124 Uppsala, Sweden
phone: +46 18 471 4537 fax: +46 18 511 755
<a class="moz-txt-link-abbreviated" href="mailto:erikm@xray.bmc.uu.se">erikm@xray.bmc.uu.se</a> <a class="moz-txt-link-freetext" href="http://folding.bmc.uu.se/">http://folding.bmc.uu.se/</a>
</pre>
</body>
</html>