AMBER Archive (2007)Subject: Re: AMBER: prepin file format description; pointers on parameterization
From: FyD (fyd_at_q4md-forcefieldtools.org)
Date: Wed Oct 24 2007 - 14:09:09 CDT
Quoting David Mobley <dmobley_at_gmail.com>:
> 1) Can anybody point me to a description of the "prepin" file format
> (i.e., the format used by all_amino94.in, or heme_all.in, for
> example)? I can't find this on the file format description page, nor
> in the manual, and I'm trying to figure out what the columns
> represent.
Did you look in plep; see "Old Prep/Link/Edit/Parm Also Spasms, a
molecular dynamics program, and Resp, for charge-fitting. Resp Q&A" on
the AMBER web site.
Even reading what you report below, I still do not understand very
well why you do want to use the old prepin file format since it should
be possible to convert this file format into either mol2 file format
or OFF force field libraries.
For mol2 file format you could see the following pages:
http://q4md-forcefieldtools.org/Tutorial/Tutorial-1.php
and in particular:
-I.3.3- The Tripos mol2 file format @
http://q4md-forcefieldtools.org/Tutorial/Tutorial-1.php#15
Please, see also the new "savemol2" command incorporated in LEaP
-I.3.3.5- Using Tripos mol2 files in LEaP: Implementation of new
features in LEaP @
http://q4md-forcefieldtools.org/Tutorial/leap.php
> At a more detailed level, the basic problem I'm trying to solve is
> that I (to my chagrin) need to simulate an artificial mutant of
> cytochrome C peroxidase, which has a heme (heme B, particularly)
> bordering an artificial binding site I'm interested in. The chemistry
> of the heme isn't particularly important here, but I do at least need
> "sane" parameters for it so it stays stable, etc.
yes...
> There are a set of heme parameters distributed with AMBER for heme,
> apparently, and I hope to be able to modify them to describe my
> particular case. They're originally for myoglobin/hemoglobin, if I
> understand correctly, and are apparently for heme with the iron
> covalently bound additionally to histidine and O2. In my particular
> case, the iron should be covalently bound to histidine and a water
> molecule. I'm trying to figure out how I need to modify the existing
> parameters to achieve this (hence question #1).
>
> This brings me to my follow-up question:
>
> 2) My rough idea of the procedure I need to use follows; is this
> roughly right? (and any tips on how to accomplish the specific steps,
> especially step (a), will be appreciated):
> (a) Modify the existing heme parameters to get it to connect to water
> (b) separately run QM at the appropriate level on the iron+heme
> (without histidine/water, probably) and perform a RESP fit to get the
> partial charges
> (c) incorporate the partial charges into the parameter files
> (d) Load the whole thing together with my pdb into leap somehow and
> write out prmtop and crd files
Looking for heme from the J. Comput. Chem. web site gives:
http://www3.interscience.wiley.com/search/allsearch?WISsearch2=heme&WISindexid2=WISall&mode=quicksearch&products=journal&contextLink=%3Ca+href%3D%22%2Findex.html%22+target%3D%22_top%22%3EHome%3C%2Fa%3E+%2F+%3Ca+href%3D%22%2Fbrowse%2F%3Fsubject%3DCHEM%22+target%3D%22_top%22%3EChemistry%3C%2Fa%3E+%2F+%3Ca+href%3D%22%2Fbrowse%2F%3Fsubject%3DCH65%26titles%3Dtrue%22+target%3D%22_top%22%3EComputational+Chemistry+and+Molecular+Modeling%3C%2Fa%3E&contentOID=33822&contentTitle=Journal+of+Computational+Chemistry&WISindexid1=issn&WISsearch1=1096-987X&articleGo.x=0&articleGo.y=0
In particular, you could look at:
Classical force field parameters for the heme prosthetic group of cytochrome c
Journal of Computational Chemistry
Volume 25, Issue 13, Date: October 2004, Pages: 1613-1622
Felix Autenrieth, Emad Tajkhorshid, Jerome Baudry, Zaida Luthey-Schulten
Read the last paragraph about RESP charge derivation, & you will see
it is tricky on heme.
In this paper, the authors compute RESP charges for the heme + _all_
the ligands. The tricky point is to get a correct minimum with a
correct QM packages and to get SCF convergence for the correct spin
multiplicty/electronic configuration... You might look at the Jaguar
package to save you a lot of Gaussian/GAMESS cpu time ;-)
Finally, I think heme electronic configuration is highly depend on the
ligands. Thus, you should think about deciding or not deciding to
remove ligands from the whole heme system in your QM computation...
> Any general pointers on this will be greatly appreciated. I've used
> leap/antechamber before for setting up normal proteins and
> parameterizing small molecules that are non-covalently attached (i.e.
> that can be done in the "normal" way through antechamber) so I'm
> familiar with those "standard" steps, but am a bit at a loss when it
> comes to how to set up and tweak nonstandard residues.
My understanding is that antechamber does not work with metals yet.
regards, Francois
PS We will release a general heme force field in R.E.DD.B.
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|