<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:12pt"><br>
<div style="font-family: arial, helvetica, sans-serif; font-size: 12pt;"><div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"><div id="yiv1060757216"><div><div style="color:#000;background-color:#fff;font-family:arial, helvetica, sans-serif;font-size:12pt;"><div><span>Thanks for your suggestions.<br><br></span></div><div>The total charge is zero in real but what I got through swissparam is as below:</div><div>[ atoms ] <br>; nr type resnr resid atom cgnr charge mass <br>&nbsp;&nbsp; 1 C=O&nbsp; 1&nbsp; LIG C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp; 0.4500&nbsp; 12.0110 <br>&nbsp;&nbsp; 2 O=C&nbsp; 1&nbsp; LIG O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 -0.5700&nbsp; 15.9994 <br>&nbsp;&nbsp; 3 HCMM 1&nbsp; LIG H1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3&nbsp; 0.0600&nbsp;&nbsp; 1.0079 <br>&nbsp;&nbsp; 4 HCMM 1&nbsp; LIG H2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4&nbsp; 0.0600&nbsp;&nbsp; 1.0079 <br></div><div>So I entered as I showed
 you.</div><div><br></div><div>In addition how can I add the proton on it? Editing the C to CT1 is right? or I have to add any charge here?</div><div><br></div><div>Formyl is connected to valine in the peptide
 sequence in this .pdb file. Now, what about adding&nbsp; this line ?:&nbsp; <br></div><div>[ dihedrals ]<br>&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp; C&nbsp;&nbsp; +N&nbsp;&nbsp; CA <br></div><div><br></div><div>Cheers,</div><div>Shima<br></div>  <div style="font-family:arial, helvetica, sans-serif;font-size:12pt;"> <div style="font-family:times new roman, new york, times, serif;font-size:12pt;"> <div dir="ltr"> <font face="Arial" size="2"> <hr size="1">  <b><span style="font-weight:bold;">From:</span></b> Justin A. Lemkul &lt;jalemkul@vt.edu&gt;<br> <b><span style="font-weight:bold;">To:</span></b> Shima Arasteh &lt;shima_arasteh2001@yahoo.com&gt;; Discussion list for GROMACS users &lt;gmx-users@gromacs.org&gt; <br> <b><span style="font-weight:bold;">Sent:</span></b> Saturday, May 12, 2012 10:07 PM<br> <b><span style="font-weight:bold;">Subject:</span></b> Re: [gmx-users] Fatal error: No atoms found in .rtp file in residue pairs<br> </font> </div> <br>
<br><br>On 5/12/12 1:32 PM, Shima Arasteh wrote:<br>&gt; Dear gmx users,<br>&gt;<br>&gt; I want to simulate a peptide in water. The peptide has a formyl residue as the<br>&gt; N-terminus. I got the parameters of it and then add it to the .rtp file of<br>&gt; charmm36.ff as below:<br>&gt;<br>&gt; [ For ]<br>&gt; [ atoms ]<br>&gt; ; name type charge chargegroup<br>&gt; C C 0.4500 0<br>&gt; O O -0.5700 0<br>&gt; [ bonds ]<br>&gt; ; ai aj fu b0 kb, b0 kb<br>&gt; 1 2 1 0.12220 779866.6 0.12220 779866.6<br>&gt;<br>&gt; [ pairs ]<br>&gt; ; ai aj fu<br>&gt;<br>&gt; [ angles ]<br>&gt; ; ai aj ak fu th0 kth ub0 kub th0 kth ub0 kub<br>&gt; 2 1 3 1 123.4390 403.48 123.4390 403.48<br>&gt; 2 1 4 1 123.4390 403.48 123.4390 403.48<br>&gt;<br>&gt; [ dihedrals ]<br>&gt; ; ai aj ak al fu phi0 kphi mult phi0 kphi mult<br>&gt;<br>&gt; But when I do what I described, I face this fatal error:<br>&gt; Fatal error:<br>&gt; No atoms found in .rtp file in residue
 pairs<br>&gt;<br>&gt; Would you help me with this problem? Did I add the formyl in a wrong way?<br>&gt;<br><br>Your [pairs] directive is empty, hence the error.&nbsp; Your [dihedrals] directive <br>will also produce the same error.&nbsp; There are no possible pairs in a unit <br>containing so few atoms, so you don't even need this.&nbsp; You may need to define <br>dihedrals, however, but they also depend upon the next residue.<br><br>You also have several other mistakes:<br><br>1. A formyl group has a proton on it; yours has none<br>2. Your net charge on the formyl group is not zero, though this may or may not <br>be a consequence of point #1<br>3. There is no connectivity information for linking to the next residue, which <br>will mean the formyl group will not be chemically bonded.&nbsp; This also affects <br>your dihedrals.<br><br>-Justin<br><br>-- <br>========================================<br><br>Justin A. Lemkul, Ph.D.<br>Department of
 Biochemistry<br>Virginia Tech<br>Blacksburg, VA<br>jalemkul[at]<a rel="nofollow" target="_blank" href="http://vt.edu">vt.edu</a> | (540) 231-9080<br>http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin<br><br>========================================<br><br><br> </div> </div>  </div></div></div><br><br> </div> </div>  </div></body></html>