RE: AMBER: Density of Cyclohexane box

From: chaiann ng <andy810915.hotmail.com>
Date: Fri, 04 Nov 2005 11:10:50 +0800

Dear Ross:

Thanks for your reply. I did use constant pressure where (ntb=2, ntp=1).
I somehow found the right way to doing it but it is rather strange.
Initially I am using antechamber to create the prepin file from the PDB file
that I download from www.molecularmodels.ca
It is these prmtop and inpcrd files create from this prepin file are not
working as the density is decreasing from the very beginning.
Colleague fo mine suggest that delete the H from the prepin file and use
Xleap to build the H back into the structure and save as PDB file, then
create a new prepin file from this new PDB file. This is time the density is
increasing to 0.74 after 150ps of MD.

So I do not really know what is going wrong but it seems to work this way.

<BLOCKQUOTE style='PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #A0C6E5
2px solid; MARGIN-RIGHT: 0px'><font
style='FONT-SIZE:11px;FONT-FAMILY:tahoma,sans-serif'><hr color=#A0C6E5
size=1>
From: <i>&quot;Ross Walker&quot;
&lt;ross.rosswalker.co.uk&gt;</i><br>Reply-To:
<i>amber.scripps.edu</i><br>To:
<i>&lt;amber.scripps.edu&gt;</i><br>Subject: <i>RE: AMBER: Density of
Cyclohexane box</i><br>Date: <i>Wed, 2 Nov 2005 20:52:32 -0800</i><br>Dear
Chaiann,<br><br> &gt; I am trying to create a cyclohexane box for my
simulation. I<br> &gt; am using the<br> &gt; input for under the folder
&quot;solvent&quot; in the amber directory<br> &gt; and also the<br> &gt;
input file from the web. I am getting density of less than<br> &gt; 0.4 at
300K after<br> &gt; 50ps MD, which is the initial density when i solvate
my<br> &gt; cyclohexane (chair<br> &gt; conformation) with a box of
cyclohexane. Can anyone help?<br> &gt; Attached is the<br> &gt; PDB file of
the cyclohexane and prepin file. Thanks.<br><br>Are you running constant
pressure (ntb=2, ntp=1)? The default for periodic<br>boundaries is to run
constant volume (ntb=1). In this case the box size will<br>be fixed at the
value setup by Leap (a not very good guess at the density)<br>and thus the
density will not change.<br><br>With constant pressure you should be able to
equilibrate to a density of<br>around 0.78 g/cm^3 or so for cyclohexane.
Although you may need more than<br>50ps of equilibration to reach an
equilibrate density.<br><br>If you still have problems please post your
prmtop, inpcrd, mdout and mdin<br>files.<br><br>All the
best<br>Ross<br><br>/\<br>\/<br>|\oss Walker<br><br>| Department of
Molecular Biology TPC15 |<br>| The Scripps Research Institute |<br>| Tel: +1
858 784 8889 | EMail:- ross.rosswalker.co.uk |<br>|
http://www.rosswalker.co.uk | PGP Key available on request |<br><br>Note:
Electronic Mail is not secure, has no guarantee of delivery, may not<br>be
read every day, and should not be used for urgent or sensitive
issues.<br><br><br>-----------------------------------------------------------------------<br>The
AMBER Mail Reflector<br>To post, send mail to amber.scripps.edu<br>To
unsubscribe, send &quot;unsubscribe amber&quot; to
majordomo.scripps.edu<br></font></BLOCKQUOTE>

_________________________________________________________________
Get an advanced look at the new version of MSN Messenger.
http://messenger.msn.com.my/Beta/Default.aspx

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Fri Nov 04 2005 - 03:53:00 PST
Custom Search