<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV id=yiv1869022878>Dear <STRONG>Yanmei Song</STRONG>:</DIV>
<DIV>I used gromacs version 3.3.1. I update ffG45a3.rtp to include my molecule. When I run pdb2gmx using my molecules with the command:</DIV>
<DIV>&nbsp;</DIV>
<DIV><STRONG>"pdb2gmx -f min_C16.pdb -o C16.pdb -p topol.top -ff G45a3"</STRONG></DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;I got the error message: 
<DIV><I>"&nbsp; Opening library file ffG45a3.rtp<BR></I><I>&nbsp; Opening library file&nbsp; /user/localgromacs/share/gromacs/top/aminoacid.dat</I></DIV>
<DIV><EM>&nbsp; Reading min_C16.pdb...&nbsp;<BR>&nbsp;', 81 atomsin in water<BR>&nbsp; Opening library file&nbsp; /user/localgromacs/share/gromacs/top/xlateat.dat<BR></EM>&nbsp; <EM>26 out of 26 lines of xlateat.dat converted succesfully<BR></EM>&nbsp; <EM>Analyzing pdb file<BR></EM>&nbsp; <EM>There are 1 chains and 0 blocks of water and 1 residues with&nbsp;81 atoms</EM></DIV>
<DIV><BR>&nbsp;&nbsp; <EM>&nbsp;chain &nbsp;#res #atoms<BR></EM>&nbsp;&nbsp;&nbsp;&nbsp; <EM>1 '&nbsp; ' &nbsp; &nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; 81<BR></EM>&nbsp; </DIV>
<DIV>&nbsp; All occupancies are one<BR>&nbsp; <EM>Opening library file ffG45a3.atp<BR></EM>&nbsp; <EM>Atomtype 49<BR></EM>&nbsp; <EM>Reading residue database... (ffG45a3)<BR></EM>&nbsp;&nbsp; <EM>Opening library file ffG45a3.rtp<BR></EM><I>&nbsp;&nbsp;Using default: not generating all posible dihedrals</I></DIV>
<DIV><EM>&nbsp; Using default:excluding 3 bonded neighbors</EM></DIV>
<DIV><EM>&nbsp; Using default: generating 1,4 H--H interactions</EM></DIV>
<DIV><EM>&nbsp;Using default: removing impropers on same bond as a proper</EM></DIV>
<DIV><I>&nbsp; -------------------------------------------------------<BR></I>&gt;&gt;&gt;&gt;<I> Program pdb2gmx, VERSION 3.3.1<BR></I>&gt;&gt;&gt;&gt;<I> Source code file: resall.c, line: 331<BR></I>&gt;&gt;&gt;&gt;<I><BR></I>&gt;&gt;&gt;&gt;<I> Fatal error:<BR></I>&gt;&gt;&gt;&gt;<I> in .rtp file at line: "</I></DIV><I></I>&nbsp;</DIV>
<DIV>The molecule I added in the ffG45a3.rtp is in attached file</DIV>
<DIV><BR><EM></EM>&nbsp;</DIV>
<DIV>Thank you very much indeed.</DIV>
<DIV><EM></EM>&nbsp;</DIV>
<DIV><EM>Angel</EM></DIV>
<DIV><EM></EM>&nbsp;</DIV></td></tr></table><br>






      &nbsp;