<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'>
Hi,<br><br>I fixed it.<br>Thanks for the fast test and the complete instructions,<br><br>Berk<br><br><hr id="stopSpelling">From: alanwilter@gmail.com<br>Date: Tue, 24 Aug 2010 15:22:34 +0100<br>To: gmx-users@gromacs.org<br>Subject: [gmx-users] pdb2gmx e82fc gmx 4.5 is failing with DNA now<br><br>Hi there, in special Berk.<div><br></div><div>So pdb2gmx may be working with HIS and variants for oplsaa but now, something that was working before is failing:</div><div><br></div><div><div>wget -c "<a href="http://www.pdbe.org/download/1BNA" target="_blank">http://www.pdbe.org/download/1BNA</a>" -O 1BNA.pdb</div>

<div>grep 'ATOM &nbsp;' 1BNA.pdb &gt;| DNA.pdb</div><div><br></div><div><div>cat &lt;&lt; EOF &gt;| SPE.mdp</div><div>define = -DFLEXIBLE</div><div>integrator &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = md</div><div>nsteps &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 0</div>

<div>dt &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 0.001</div><div>constraints &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= none</div><div>emtol &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 10.0</div><div>emstep &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 0.01</div><div>nstcomm &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1</div><div>ns_type &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= simple</div>

<div>nstlist &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 0</div><div>rlist &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 0</div><div>rcoulomb &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 0</div><div>rvdw &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 0</div><div>Tcoupl &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = no</div><div>Pcoupl &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = no</div>

<div>gen_vel &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= no</div><div>nstxout &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1</div><div>pbc &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= no</div><div>nstlog = 1</div><div>nstenergy = 1</div><div>nstvout = 1</div><div>nstfout = 1</div><div>nstxtcout = 1</div>

<div>comm_mode = ANGULAR</div><div>continuation = yes</div><div>EOF</div></div><div><br></div><div><div>pdb2gmx -f DNA.pdb -o DnaAmberSBGMX45.pdb -ff amber99sb -water none -p DnaAmberSBGMX45</div><div>[snip]</div></div><div>

<div>8 out of 8 lines of specbond.dat converted successfully</div><div>[1] &nbsp; &nbsp;42209 segmentation fault &nbsp;pdb2gmx -f DNA.pdb -o DnaAmberSBGMX45.pdb -ff amber99sb -water none -p&nbsp;</div><div><br></div></div>with&nbsp;5e347 it worded fine, i.e, it opens files</div>

<div><br></div><div><div>Opening force field file /Volumes/CloneAmadeus/usr/local/share/gromacs/top/amber99sb.ff/aminoacids.arn</div><div>Opening force field file /Volumes/CloneAmadeus/usr/local/share/gromacs/top/amber99sb.ff/dna.arn</div>

<div>Opening force field file /Volumes/CloneAmadeus/usr/local/share/gromacs/top/amber99sb.ff/rna.arn</div></div><div><br></div><div>and proceed.</div><div><br></div><div>Alan<br>-- <br>Alan Wilter S. da Silva, D.Sc. - CCPN Research Associate<br>

Department of Biochemistry, University of Cambridge. <br>80 Tennis Court Road, Cambridge CB2 1GA, UK.<br>&gt;&gt;<a href="http://www.bio.cam.ac.uk/%7Eawd28" target="_blank">http://www.bio.cam.ac.uk/~awd28</a>&lt;&lt;<br>
</div>
<br>-- 
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-request@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php                                               </body>
</html>