<br><br>----- Original Message -----<br>From: Sai Pooja &lt;saipooja@gmail.com&gt;<br>Date: Tuesday, October 26, 2010 7:26<br>Subject: [gmx-users] CT3 terminal<br>To: Discussion list for GROMACS users &lt;gmx-users@gromacs.org&gt;<br><br><div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>Hi,</div> <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>&nbsp;</div> <div><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>I use pdb2gmx command to convert a pdb file to a .gro and .top file(using CHARMM forcefield). For the C terminus, I use the CT3 option which is NCH3. When I use editconf to get a pdbfile from the .gro file generated, the C=0 and the NCH3&nbsp; seem to clash in space in pymol.<br><br></div>It looks like Par has borrowed the mechanism for generating two carboxylate oxygens on a C atom (one each cis and trans wrt the third generating atom) for generating one amide nitrogen, and the code doesn't check which of cis and trans is already occupied by the carboxyl oxygen. Please file a Bugzilla so that this gets fixed.<br><br>A real solution requires some coding, but you can get a crude work-around by replacing<br><br>"1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NT&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; CA&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; N" in share/top/charmm27.ff/aminoacids.c.tdb with<br>"1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NT&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; CA&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; O"<br><br>See chapter 5 for details.<br><br>This NT will be too close to carboxyl C, and EM will be needed to fix things.<br><br>Mark