MOLECULAR MECHANICS and Force Field Calculations
MOLECULAR MECHANICS and Force Field Calculations
MOLECULAR MECHANICS and Force Field Calculations
MOLECULAR MECHANICS
AND
FORCE FIELD CALCULATIONS
BY
SUSNEHA PATRA
1
CONTENTS
Molecular mechanics
Force field
Functional form
Parametrization
Types of force fields
An example
Areas of application
Advantages and disadventages
Refernces
2
MOLECULAR MECHANICS
3
All-atomistic molecular mechanics methods have
the following properties:
Each atom is simulated as a single particle
4
The potential energy of all systems in
molecular mechanics is calculated using
force fields.
5
FORCE FIELD
In the context of molecular modeling, a force field refers
to the form and parameters of mathematical
functions used to describe the potential energy of a
system of particles (typically molecules and atoms).
2 types:
"All-atom" force fields
"united-atom" force fields
"Coarse-grained" force fields, which are frequently
used in long-time simulations of proteins, provide even
more crude representations for increased computational
efficiency.
6
A force field is used
To minimize the bond
streching energy of this
Ethane molecule
7
FUNCTIONAL FORM
Molecular system's potential energy (E) in a given
conformation is expressed as a sum of individual
energy terms.
E = E(covelent) + E(noncovelent)
where ,
8
Some of the energy terms that need to be
taken into account are:
(1) Bond stretching (l)
9
The exact functional form of the potential
function, or force field, depends on the particular
simulation program being used.
bond and angle terms are modeled
as harmonic potentials
For accurate reproduction of vibrational
spectra, the Morse potential can be used
instead
the van der Waals term falls off rapidly – it is
typically modeled using a "6–12 Lennard-Jones
potential“
Electrostatic term with coulomb potential
10
Molecular mechanics
Potential energy
function
with continuum
Solvent
11
PARAMETRIZATION
A force field defines a set of parameters for
each type of atom.
The typical parameter set includes
values for atomic mass, van der Waals radius,
and partial charge for individual atoms, and
equilibrium values ofbond lengths, bond
angles, and dihedral angles for pairs, triplets,
and quadruplets of bonded atoms, and values
corresponding to the effective spring
constant for each potential.
12
TYPES OF FORCE FIELDS
Classical force fields
AMBER (Assisted Model Building and Energy
Refinement) - widely used for proteins and DNA.
CHARMM (Chemistry at HARvard Molecular
Mechanics) - originally developed at Harvard,
widely used for both small molecules and
macromolecules
CHARMm - commercial version of CHARMM,
available through Accelrys.
13
Second-generation force fields
CFF - a family of forcefields adapted to a broad
variety of organic compounds, includes force
fields for polymers, metals, etc.
MMFF (Merck Molecular Force Field)-
developed at Merck, for a broad range of
molecules.
14
Coarse-grained Force fields
VAMM (Virtual atom molecular mechanics)
Polarizable force fields based on electronic
structural theory
X-Pol: the Explicit Polarization Theory[1]
Other
VALBOND - a function for angle bending that
is based on valence bond theory and works
for large angular distortions, hypervalent
molecules, and transition metal complexes
15
EXAMPLE
TABLE V. Computed and experimental fundamentals ~cm21! for various isomers of C2H4
16
Calculated
n1 3031.4a 3014.8c 2294.2 2283.0 2277.4j 2260.3
n2 1602.1 1582.3 1567.9 1568.5 1544.5 1513.8l
n3 1286.0 1029.1 1215.8 1283.0g 1041.8 984.1
n4 999.6 890.4 975.0 988.0 757.4 729.2
n5 3050.6 2328.0 3050.6d 3040.3h 2217.5 2308.2
n6 1121.6 1137.5 1035.7 1000.0 996.0 999.6
n7 803.0 746.7 842.4 724.7 724.0 720.0
n8 937.5 937.3 754.2 853.0 915.1 768.8
n9 3089.9 3087.6 3053.9e 3059.3 3041.5k 2337.1
n10 728.2 680.9 659.5 669.4 625.6 590.0
n11 2267.5b 2227.9 2255.4f 2270.4j 2327.0 2197.3
n12 1397.9 1379.8 1338.1 1296.0 1285.6 1075.8
17
AREAS OF APPLICATION
18
Global energy optimization can be
accomplished using simulated annealing,
the Metropolis algorithm and other Monte
Carlo methods
The main aim of optimization methods is
finding the lowest energy conformation of a
molecule or identifying a set of low-energy
conformers that are in equilibrium with each
other.
19
ADVANTAGES
20
DISADVANTAGES
21
REFERENCES
http://en.wikipedia.org/wiki/Force_field_(ch
emistry)
http://chemistry.ncssm.edu/book/Chap6Mo
lMech.pdf
evgenii.rudnyi.ru/doc/teaching/md/md7.pdf
http://theochem.weizmann.ac.il/AIPreprints
/69.pdf
linkinghub.elsevier.com/retrieve/pii/002228
5271900087
22
THANK
YOU
23