<HTML dir=ltr><HEAD><TITLE>Re: [gmx-users] modified atomtypes.atp and grompp</TITLE>
<META http-equiv=Content-Type content="text/html; charset=unicode">
<META content="MSHTML 6.00.2800.1649" name=GENERATOR></HEAD>
<BODY>
<DIV id=idOWAReplyText98030 dir=ltr>
<DIV dir=ltr><FONT face=Arial color=#000000 size=2>Dear Justin,</FONT></DIV>
<DIV dir=ltr><FONT face=Arial size=2></FONT>&nbsp;</DIV>
<DIV dir=ltr><FONT face=Arial size=2>I checked my working directory (including ff-directory) again and found some entries (i added myself) with the old uppercase spelling in the nonbonded.itp. After editing these lines everything works well. Conclusion: I have to learn using the command grep properly :)</FONT></DIV>
<DIV dir=ltr><FONT face=Arial size=2></FONT>&nbsp;</DIV>
<DIV dir=ltr><FONT face=Arial size=2>Thank you very much for your help.</FONT></DIV>
<DIV dir=ltr><FONT face=Arial size=2>Felix</FONT></DIV></DIV>
<DIV dir=ltr><BR>
<HR tabIndex=-1>
<FONT face=Tahoma size=2><B>Von:</B> gmx-users-bounces@gromacs.org im Auftrag von Justin A. Lemkul<BR><B>Gesendet:</B> Fr 11.03.2011 13:07<BR><B>An:</B> Discussion list for GROMACS users<BR><B>Betreff:</B> Re: [gmx-users] modified atomtypes.atp and grompp<BR></FONT><BR></DIV>
<DIV><BR><BR>
<P><FONT size=2>Rausch, Felix wrote:<BR>&gt;&nbsp;<BR>&gt; Dear Gromacs users,<BR>&gt;&nbsp;<BR>&gt; I would like so simulate a protein with phosphorylated amino acids with<BR>&gt; Gromacs 4.5.3. To do so, I got the "ffG43a1p" forcefield available from<BR>&gt; the user contributions section (subsection force fields) and adapted the<BR>&gt; parameters given there to Gromos53a6.<BR>&gt; The problem now is that they introduced a new atom type called "OP",<BR>&gt; which I have to add to the atomtypes.atp of my modified force field. I<BR>&gt; already read about problems of a modified atomtypes.atp here in the<BR>&gt; mailing list, but I assumed the problems were sorted out with version 4.5.3.<BR><BR>Not entirely, but probably the post you saw related to an LP atomtype, which<BR>appears to be a special case.<BR><BR>&gt; But if I add OP to the atom list, pdb2gmx works fine, but grompp is<BR>&gt; raising the error<BR>&gt;&nbsp;<BR>&gt; Program grompp, VERSION 4.5.3.<BR>&gt; Source code file: toppush.c, line: 629<BR>&gt;&nbsp;<BR>&gt; Fatal error:<BR>&gt; Unknown atomtype CMET<BR>&gt;&nbsp;<BR>&gt; which is quite strange because&nbsp; "CMET" isn't appearing in any file I'm<BR>&gt; using. When I remove OP from the atom list, pdb2gmx is of course raising<BR>&gt; an error, so I'm a bit stuck here...<BR>&gt;&nbsp;<BR><BR>grompp reads all the force field files and is finding something it doesn't like.<BR>&nbsp; In 4.5.3, the proper capitalization should be "CMet," not "CMET," as it was in<BR>old versions.&nbsp; Capitalization and other nomenclature has been standardized as of<BR>version 4.5.3, so old files with outdated names cause errors.<BR><BR>-Justin<BR><BR>&gt; To cut a long story short: Is there any way to introduce a new atom type<BR>&gt; to atomtypes.atp or do I have to parameterize the phosphorylated<BR>&gt; residues with already presend atom types?&nbsp;<BR>&gt;&nbsp;<BR>&gt; Thank you in advance,<BR>&gt; Felix Rausch<BR>&gt;<BR>&gt;&nbsp;<BR>&gt;<BR><BR>--<BR>========================================<BR><BR>Justin A. Lemkul<BR>Ph.D. Candidate<BR>ICTAS Doctoral Scholar<BR>MILES-IGERT Trainee<BR>Department of Biochemistry<BR>Virginia Tech<BR>Blacksburg, VA<BR>jalemkul[at]vt.edu | (540) 231-9080<BR><A href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</A><BR><BR>========================================<BR>--<BR>gmx-users mailing list&nbsp;&nbsp;&nbsp; gmx-users@gromacs.org<BR><A href="http://lists.gromacs.org/mailman/listinfo/gmx-users">http://lists.gromacs.org/mailman/listinfo/gmx-users</A><BR>Please search the archive at <A href="http://www.gromacs.org/Support/Mailing_Lists/Search">http://www.gromacs.org/Support/Mailing_Lists/Search</A> before posting!<BR>Please don't post (un)subscribe requests to the list. Use the<BR>www interface or send it to gmx-users-request@gromacs.org.<BR>Can't post? Read <A href="http://www.gromacs.org/Support/Mailing_Lists">http://www.gromacs.org/Support/Mailing_Lists</A><BR></FONT></P></DIV></BODY></HTML>