Re: parm99

From: Junmei Wang <jwang.tbc.com>
Date: Wed, 21 May 2003 13:30:18 -0500

Hi, Dave,
I have tested the antechamber code for 40 molecules used by Christ Bayly.
The new scheme works fine. The programs are much faster and the algorithms
are more clear. I will load the revised code to chagall.scripps.edu after I
finish testing another 80 molecules. I will also add some force field
parameters to gaff.dat and then revise the papers ... Hopefully, I can
finish all the staff in two weeks.

As to the atom type of C-OH in tyrosine, I prefer to change its atom type
from C to CA in parm94 and parm96 based on the following reasons:
(1) first of all, we do not need make large change on parm94.dat and
parm96.dat (only change one bond length from 1.409 to 1.4)
(2) consistent with ff02 and parm99
(3) consistent with the atom type definition

The following is the parameters involved in parm94

C -CA 469.0 1.409 JCC,7,(1986),230; TYR
CA-CA 469.0 1.400 JCC,7,(1986),230; BENZENE,PHE,TRP,TYR
C -OH 450.0 1.364 JCC,7,(1986),230; TYR
CA-C -CA 63.0 120.00 changed from 85.0 bsd on C6H6 nmodes; AA
tyr
CA-C -OH 70.0 120.00 AA tyr
C -CA-CA 63.0 120.00 changed from 85.0 bsd on C6H6 nmodes
CA-CA-CA 63.0 120.00 changed from 85.0 bsd on C6H6 nmodes
X -C -CA-X 4 14.50 180.0 2. intrpol.bsd.on
C6H6
X -CA-CA-X 4 14.50 180.0 2. intrpol.bsd.on
C6H6
X -C -OH-X 2 1.80 180.0 2.
JCC,7,(1986),230

If we change C to CA, only bond length CA-C(CA) changed slightly from 1.409
to 1.400. Other parameters are exactly same.

If we want to do this, we need add three additional parameters

CA-OH 450.0 1.364 JCC,7,(1986),230; TYR
CA-CA-OH 70.0 120.00 AA tyr
X -CA-OH-X 2 1.80 180.0 2.
JCC,7,(1986),230

Best

Junmei


===============================================================
Dr. Junmei Wang
Chemistry & Biophysics
Encysive Pharmaceuticals
7000 Fannin, Houston TX 77030
Tel: 713-5786649
Email: jwang.tbc.com
Homepage: Http://sigyn.compchem.ucsf.edu/members/jmwang/index.html
===============================================================



                                                                           
             "David A. Case"
             <case.scripps.edu
> To
                                       cailliez <Fabien.Cailliez.ibpc.fr>
             05/20/2003 05:37 cc
             PM amber.heimdal.compchem.ucsf.edu,
                                       jwang.tbc.com
                                                                   Subject
             Please respond to Re: parm99
             amber.heimdal.com
              pchem.ucsf.edu
                                                                           
                                                                           
                                                                           
                                                                           




On Tue, May 20, 2003, cailliez wrote:

> There is still something I do not understand with parm99. In parm99.dat,
> there are lots of new parameters, including some new parameters for
tyrosine
> concerning valence angles or torsionnal angles. Reading parm99.dat, I
> understand that these changes for tyrosine must be done by changing the
type
> of CZ atom from C (in parm98) to CA (in parm99).

This looks like a bug to me: there doesn't seem to be any easy way to
use the parm99 tyrosine parameters for non-polarizable simulations, since
all_amino94.in still has the old atom types. [The polarizable libraries,
like all_amino02.in and all_amino02EP.in *do* have this change.]

Junmei: do you have a recommended solution for this?

The only workaround I see is to manually modify the TYR unit to change the
atom type of CZ, before loading in your pdb file.

Thanks for the report ... we'll try to figure out the best fix (so as not
to break older simulations that use ff94). This change is probably pretty
minor, but definitely looks to me like it is wrong.

...regards...dac

--
==================================================================
David A. Case                     |  e-mail:      case.scripps.edu
Dept. of Molecular Biology, TPC15 |  fax:          +1-858-784-8896
The Scripps Research Institute    |  phone:        +1-858-784-9768
10550 N. Torrey Pines Rd.         |  home page:
La Jolla CA 92037  USA            |    http://www.scripps.edu/case
==================================================================
Received on Wed May 21 2003 - 19:53:00 PDT
Custom Search