re: Mark's comment, pdb2gmx doesn't seem to be able to do the atom
typing, etc., necessary to be able to deal with general small
molecules. Or if it can in some cases, at least not the cases in which
I've tried it. :) I think it's geared primarily towards proteins. But
Mark is right -- the output from PDB in .itp format is not all-atom.
You can get nominally all-hydrogen mol2 files out but then you still
have to have some way of converting these to gromacs topologies. And it
sometimes leaves out hydrogens, too, in my experience.<br><br><div><span class="gmail_quote">On 11/17/05, <b class="gmail_sendername">Mark Abraham</b> &lt;<a href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>
&gt; wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">TJ Piggot wrote:<br>&gt; Hi I sent this email a few weeks ago and have had no replies yet I hope
<br>&gt; that someone will be able to help this time.<br>&gt;<br>&gt; I am trying to construct ATP and ADP topology files for use with the opls<br>&gt; force field. I have previously been able to do this for the ffG43a1 force
<br>&gt; field using prodrg.<br>&gt;<br>&gt; However prodrg does not produce output files for use with the opls force<br>&gt; field. I have tried to change the prodrg output files (.itp) for ADP and<br>&gt; ADP but the problem is that there is not enough atoms as the individual
<br>&gt; hydrogens are not defined.<br><br>OPLS/AA is an all-atom force field, whereas the GROMOS force fields are<br>not. Perhaps you should review why you are using whatever force field<br>you think you should be using. I have no idea what prodrg does, but if
<br>it writes OPLS/AA-suitable topologies and you want to go along this path<br>then you should go ahead and make prodrg do so. Otherwise you need to<br>look at alternatives, like generating a PDB structure (with or without
<br>hydrogens) and, if that structure is without hydrogens, then using some<br>tool to generate hydrogens (e.g. gromacs' pdb2gmx?)<br><br>Mark<br>_______________________________________________<br>gmx-users mailing list<br>
<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a href="http://www.gromacs.org/mailman/listinfo/gmx-users">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please don't post (un)subscribe requests to the list. Use the
<br>www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br></blockquote></div><br>