Dear Fabricio,<div><br></div><div>It&#39;s not only you who is confused here...</div><div><br></div><div>So, I will try to help the best I can and sorry if I ask you redundant questions you&#39;ve likely answer before. It happens that I read GMX mailing list in digest mode.</div>
<div><br></div><div>BTW, when replying gmx list, cc to my email (alanwilter _at gmail . com) and so you can have a prompt reply.<br><br></div><div>AFAIU, you have a protetin (receptor.itp?) and a ligand, right? And your ligand is not covalently bond to the protein, isn&#39;t it?</div>
<div><br></div><div>If so, have you see this?:</div><div><br></div><div><a href="http://code.google.com/p/acpypi/wiki/TutorialAcpypi4Gromacs">http://code.google.com/p/acpypi/wiki/TutorialAcpypi4Gromacs</a><br></div><div><br>
</div><div>Have you seen this note?:</div><div><br></div><div><strong>NB(1):</strong>&nbsp;<tt>#include &quot;Ligand.itp&quot;</tt>&nbsp;has to be inserted right after&nbsp;<tt>ffamber**.itp</tt>&nbsp;line and before&nbsp;<tt>Protein_*.itp</tt>&nbsp;line in&nbsp;<i>Complex.top</i>.<br>
</div><div><br></div><div>It seems that you missing the difference between *.top and *.itp files. So, as Justin said, read chapter 5 of GMX manual.</div><div><br></div><div>Another alternative, if you pretty familiar to AMBER and tleap is to build your system there and once you have your prmtop and inpcrd files, use acpypi to convert them to GMX files and put run EM and MD.</div>
<div><br></div><div>If all I am saying here is still not helping, consider to send me your input PDB file (protein and complex) and let me run it here.</div><div><br></div><div>Just one more thing, would you mind where are you from and what&#39;s your institute?</div>
<div><br></div><div>Many thanks in advance,</div><div><br></div><div>Alan</div><div><br><div class="gmail_quote">On Thu, Dec 4, 2008 at 04:01,  <span dir="ltr">&lt;<a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><br>
I am working with a system using ffamber99 forcefield and trying to<br>
set up a ligand protein complex system. I have encountered some<br>
problems with grompp now.<br>
In my past email, the problem was that : &gt; Fatal error: Found a second<br>
defaults directive, file &quot;ligand.itp&quot;, line 5. Then the following<br>
solution was given:<br>
<br>
&quot; You have to combine all atomtypes into one file something like<br>
<br>
filename ff_atomtypes.itp<br>
<br>
#define FF_ATOMTYPES<br>
<br>
[ defaults ]<br>
; nbfunc &nbsp; &nbsp; &nbsp; &nbsp;comb-rule &nbsp; &nbsp; &nbsp; gen-pairs &nbsp; &nbsp; &nbsp; fudgeLJ fudgeQQ<br>
1 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 2 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; yes &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 0.5 &nbsp; &nbsp; 0.8333<br>
<br>
<br>
[ atomtypes ]<br>
; name &nbsp; &nbsp; &nbsp; &nbsp;mass &nbsp; &nbsp; &nbsp;charge &nbsp; &nbsp;ptype &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; sigma<br>
epsilon<br>
<br>
 &nbsp; &nbsp; S &nbsp; &nbsp; &nbsp; 0.0000 &nbsp;0.0000 &nbsp;A &nbsp; 3.56359e-01 &nbsp;1.04600e+00<br>
 &nbsp; CC &nbsp; &nbsp; &nbsp;0.0000 &nbsp;0.0000 &nbsp;A &nbsp; 3.39967e-01 &nbsp;3.59824e-01<br>
 &nbsp; HO &nbsp; &nbsp; &nbsp;0.0000 &nbsp;0.0000 &nbsp;A &nbsp; 0.00000e+00 &nbsp;0.00000e+00<br>
 &nbsp; N3 &nbsp; &nbsp; &nbsp;0.0000 &nbsp;0.0000 &nbsp;A &nbsp; 3.25000e-01 &nbsp;7.11280e-01<br>
all atom types from amber force field 99 and also from your ligand (check if<br>
parameters are ok).<br>
<br>
<br>
delete atomtypes and defaults entry in your ligand.itp<br>
<br>
load in your .top<br>
<br>
#include &quot;ff_atomtypes.itp&quot;<br>
<br>
#include &quot;receptor.itp&quot;<br>
#include &quot;ligand.itp&quot;<br>
&quot;<br>
I tried to do this and ended up with the following .top file:<br>
<br>
; 15d-pgj2_GMX.top created by acpypi on Wed Dec &nbsp;3 14:10:55 2008<br>
<br>
[ defaults ]<br>
; nbfunc &nbsp; &nbsp; &nbsp; &nbsp;comb-rule &nbsp; &nbsp; &nbsp; gen-pairs &nbsp; &nbsp; &nbsp; fudgeLJ fudgeQQ<br>
1 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 2 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; yes &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 0.5 &nbsp; &nbsp; 0.8333<br>
<br>
; Include 15d-pgj2_GMX.itp topology<br>
#include &quot;15d-pgj2_GMX.itp&quot;<br>
<br>
[ system ]<br>
&nbsp;15d-pgj2<br>
<br>
[ molecules ]<br>
; Compound &nbsp; &nbsp; &nbsp; &nbsp;nmols<br>
&nbsp;15d-pgj2 &nbsp; &nbsp; &nbsp; &nbsp; 1<br>
<br>
And all the atomtypes are in the 15d-pgj2_GMX.itp, which is too big to<br>
put here. But is present in past posts.<br>
Ok, so no grompp tells me the same error it did before. Fatal error:<br>
Found a second defaults directive, file &quot;15d-pgj2_GMX.top&quot;<br>
<br>
When I substitute the 15d-pgj2_GMX.top for the 15d-pgj2_GMX.itp in my<br>
protein.top grompp gives out the following error:<br>
<br>
Fatal error:No such moleculetype PJ2<br>
<br>
Now I am quite confused and would like some advice on the matter.<br>
Thank you<br>
Fabrício Bracht<br>
<br></blockquote></div><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/~awd28">http://www.bio.cam.ac.uk/~awd28</a>&lt;&lt;<br>

</div>