AW: [gmx-users] modified atomtypes.atp and grompp

[Rausch, Felix]

Dear Justin,
 
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 :)
 
Thank you very much for your help.
Felix

________________________________

Von: gmx-users-bounces at gromacs.org im Auftrag von Justin A. Lemkul
Gesendet: Fr 11.03.2011 13:07
An: Discussion list for GROMACS users
Betreff: Re: [gmx-users] modified atomtypes.atp and grompp





Rausch, Felix wrote:
> 
> Dear Gromacs users,
> 
> I would like so simulate a protein with phosphorylated amino acids with
> Gromacs 4.5.3. To do so, I got the "ffG43a1p" forcefield available from
> the user contributions section (subsection force fields) and adapted the
> parameters given there to Gromos53a6.
> The problem now is that they introduced a new atom type called "OP",
> which I have to add to the atomtypes.atp of my modified force field. I
> already read about problems of a modified atomtypes.atp here in the
> mailing list, but I assumed the problems were sorted out with version 4.5.3.

Not entirely, but probably the post you saw related to an LP atomtype, which
appears to be a special case.

> But if I add OP to the atom list, pdb2gmx works fine, but grompp is
> raising the error
> 
> Program grompp, VERSION 4.5.3.
> Source code file: toppush.c, line: 629
> 
> Fatal error:
> Unknown atomtype CMET
> 
> which is quite strange because  "CMET" isn't appearing in any file I'm
> using. When I remove OP from the atom list, pdb2gmx is of course raising
> an error, so I'm a bit stuck here...
> 

grompp reads all the force field files and is finding something it doesn't like.
  In 4.5.3, the proper capitalization should be "CMet," not "CMET," as it was in
old versions.  Capitalization and other nomenclature has been standardized as of
version 4.5.3, so old files with outdated names cause errors.

-Justin

> To cut a long story short: Is there any way to introduce a new atom type
> to atomtypes.atp or do I have to parameterize the phosphorylated
> residues with already presend atom types? 
> 
> Thank you in advance,
> Felix Rausch
>
> 
>

--
========================================

Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================
--
gmx-users mailing list    gmx-users at gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-request at gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20110311/9baad5d7/attachment.html>