Re: CCL:Amber parameters and charges
Bernard Pirard of Zeneca inquired about the availability of AMBER
parameters for treating generic organic molecules. The primary
reference for the "official" AMBER force field is:
CORNELL WD; CIEPLAK P; BAYLY CI; GOULD IR; and others.
A SECOND GENERATION FORCE FIELD FOR THE SIMULATION OF PROTEINS, NUCLEIC
ACIDS, AND ORGANIC MOLECULES.
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1995 MAY 17, V117 N19:5179-5197.
The set of parameters described therein represents a complete set
for modelling proteins and nucleic acids containing standard residues.
The force field is also described as being appropriate for organic
molecules, however, the set of parameters available for modelling
such molecules is far from complete.
When studying the interactions between small molecules and proteins
or nucleic acids, bonded parameters can often be borrowed from other
force fields. Charges for a new molecule, however, need to be calculated
from an ab initio generated HF/6-31G* wave function using the RESP protocol
as described in these references:
BAYLY CI; CIEPLAK P; CORNELL WD; KOLLMAN PA.
A WELL-BEHAVED ELECTROSTATIC POTENTIAL BASED METHOD USING CHARGE
RESTRAINTS FOR DERIVING ATOMIC CHARGES - THE RESP MODEL.
JOURNAL OF PHYSICAL CHEMISTRY, 1993 OCT 7, V97 N40:10269-10280.
CORNELL WD; CIEPLAK P; BAYLY CI; KOLLMAN PA.
APPLICATION OF RESP CHARGES TO CALCULATE CONFORMATIONAL ENERGIES, HYDROGEN
BOND ENERGIES, AND FREE ENERGIES OF SOLVATION.
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1993 OCT 20, V115 N21:9620-9631.
Variations of the AMBER force field have been implemented in various
modelling packages and typically include the bonded and VDW parameters
but use charges derived using a simpler and less costly approach.
These implementations may be quite useful but it should be noted that
they differ from the original force field. The "bona fide" AMBER
force field was optimized to perform well at calculating relative energies
of interaction in solution rather than, for example, screening large
sets of compounds in a DOCKing procedure.
For those who want to carry out MD studies in solution involving
functional groups for which parameters are not available in the
JACS reference, members of the Kollman research group as well as
other members of the AMBER community often have derived such parameters.
We are currently working on developing additional parameters, however,
it will still fall to the user to calculate charges for a new molecule.
Cheers,
=Wendy Cornell
_________________________________________________________________________________
Wendy D. Cornell, Ph.D.
Senior Scientist
Biomolecular Structure and Drug Design Chemistry
Parke-Davis Research
Ann Arbor, MI 48105 USA