Orbifold notation: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
en>Giftlite
 
en>Muale
→‎Wallpaper groups: Correction the wallpaper groups of some pictures
Line 1: Line 1:
== Ray Ban Wayfarer ICH m jeg gerne ville. ==
{{Other uses|Force field (disambiguation)}}


Hvis din krop forbliver i konstant alarmberedskab, kan dine cortisol niveauer aldrig falde natten, forstyrrer din søvn. Under 30 CEO har også en liste over de bedste byer for unge iværksættere .. For det meste, er det harmløst, siger Dowdell bortset fra den langsigtede effekt på vores søvn sundhed . <br><br>Det kunne være meget værre, og sandsynligvis vil være inden længe .. Køb billetter enten på billetlugen på modtagelser niveau eller på de gule maskiner (med euromønter eller med et betalingskort, der bruger en "PIN" koden). Hinke, hvis administrerende direktør Troy Lum er Mei partner, vil frigive filmen i Australien og New Zealand .. <br><br>At skabe noget godt billede, skal du kende dit emne at være i stand til at fange sin essens. De tager deres publikum til "ubehageligt, men nødvendigt," steder at udsætte dem for realiteterne på arbejdspladsen, uden at "corporate [http://www.viking-stamps.dk/Email/emotion.asp Ray Ban Wayfarer] sukker belægning." . <br><br>Århundrede fvt, og kun 32 af de mere end tusind, der blev udført i løbet af det 5. Wilcox væsentligt forbedret den codec retningsbestemte forbedrede svartider og diskrethed CBS SQ System .. Jeg lovede at glatte tingene ud. I 1970'erne, ornerne var igen til stede i Danmark og Sverige som følge af en befolkning, der er undsluppet fra fangenskab. <br><br>Jeg bare don vide, at endnu. Tesla lavede overskrifter i 2013 for succes i sin Model S, skyhøje aktiekursen, bil brand kontroverser, og administrerende direktør Elon Musk tager på New York Times i en ugunstig gennemgang. Han er i Gambia. Du kan forblive afslappet og fokuseret ved at tillade dig selv at være åben over for den mulighed, [http://www.flex-godning.dk/menumachine/core/banner.asp Nike Dk] at dine intentioner vil komme til dig på måder, du aldrig havde drømt muligt. <br><br>Det findes i mange almindelige udtryk (ICH m jeg gerne ville.) Og anvendes til at udtrykke tvivl eller høflighed. Bedste replik i bogen er Narcissisme avler en ny narcissist generation såvel som folk, der synes at være magneter for denne personlighed type.As en introduktion entry level til narcissisme og dens virkninger dette er et godt valg for dem, der måtte have lidt i hænderne af grandiose tænkning og / eller skamløse personer med dårlig boundaries.I fundet strukturen af ​​"Seven Deadly Sins" for at [http://www.marinedesign.dk/application/header.php Nike Air Force 1] være en nyttig og kreativ måde at forklare læserne, hvad der kan være i alvorlige tilfælde, en mental disorder.In Dr. <br><br>Texas ventede blot to timer efter Shelby afgørelse at meddele, at staten ville begynde at håndhæve sin vælgertilslutning ID lov, som tidligere var blevet blokeret i henhold til § 5, og at dens diskriminerende redistricting plan ville [http://www.flex-godning.dk/menumachine/flexmenu/preview/cache.asp Vibram Fivefingers Kso] også gå i kraft..<ul>
[[File:Bond stretching energy.png|thumb|right|A force field is used to minimize the bond stretching energy of this ethane molecule.]]
 
  <li>[http://verdamilio.net/tonio/spip.php?article227/ http://verdamilio.net/tonio/spip.php?article227/]</li>
 
  <li>[http://cgi.www5c.biglobe.ne.jp/~kk_aoi/bbs/apeboard_plus.cgi/ http://cgi.www5c.biglobe.ne.jp/~kk_aoi/bbs/apeboard_plus.cgi/]</li>
 
  <li>[http://metransparent.nfrance.com/~k1001/spip.php?article8359&lang=ar&id_forum=8701/ http://metransparent.nfrance.com/~k1001/spip.php?article8359&lang=ar&id_forum=8701/]</li>
 
  <li>[http://www.tssfoundation.com/tssdemo/index.php/forum/15-upcoming-events/121124-vibram-fivefingers-kso-hvis-du-onsker-det#121124 http://www.tssfoundation.com/tssdemo/index.php/forum/15-upcoming-events/121124-vibram-fivefingers-kso-hvis-du-onsker-det#121124]</li>
 
  <li>[http://www.dysh.net/bbs?mod=spacecp&ac=blog&blogid= http://www.dysh.net/bbs?mod=spacecp&ac=blog&blogid=]</li>
 
</ul>


== Ray Ban Clubmaster  Brokaw stillede det første spørgsmål ==
In the context of [[molecular modeling]], a '''force field''' refers to the [[Mathematical formula|form]] and [[parameter]]s of [[Function (mathematics)|mathematical functions]] used to describe the [[potential energy]] of a system of particles (typically [[molecule]]s and [[atom]]s). Force field functions and parameter sets are derived from both experimental work and high-level [[quantum mechanics|quantum mechanical]] calculations. "All-atom" force fields provide parameters for every type of atom in a system, including [[hydrogen]], while "united-atom" force fields treat the hydrogen and [[carbon]] atoms in each [[methyl group|terminal methyl]] and each [[methylene bridge]] as a single interaction center. "Coarse-grained" force fields, which are frequently used in long-time simulations of [[protein]]s, provide even more crude representations for increased computational efficiency.


Hidtil kritikere siger, har strategien været ineffektiv. De sagde [http://www.viking-stamps.dk/Email/emotion.asp Ray Ban Clubmaster] nej. (I hans fravær, Det Hvide Hus torsdag lancerede en hjemmeside dedikeret til at gøre tilfældet for en Syrien strejke.) I et brev til Obama, opfordrede præsidenten til at "formulere den trussel mod vores nationale sikkerhed, til [http://www.flex-godning.dk/menumachine/flexmenu/preview/cache.asp Vibram Fivefingers Kso] at formulere din plan om at eliminere, at trussel og søge folkets godkendelse for din handling. <br><br>Min kone er på F1 og im høvling at ansøge om cos til F2 i januar 18,2013. Som du kan se, Ryan elskede at være en ring ihændehaver. Hvad en blast vi havde! Hvad en mini saga! Vi først forelagt Apothecary [http://www.marinedesign.dk/application/header.php Nike Air Force 1 Dame] tilbage i sensommeren 2011, overbeviste om at vi havde undgået alle større rifter og omhyggelig med at holde app frække bits især et mildt tung vejrtrækning afsnittet, rammende, MOAN sikkert inden for (hvad vi troede var) parametrene for Apples 17 + rating.. <br><br>Jeg har skrevet om det før behovet for forældre at få tid til hinanden, mens hæve en familie. Mitchell fra University of Pennsylvania i et papir for National Bureau of Economic Research. Washington og Lee pensum blander traditionelle liberale kunst og videnskab med præ professionelle programmer i erhvervslivet, journalistik og jura, der giver eleverne et nutidigt perspektiv er nødvendige for at trives i en kompleks verden.. <br><br>Jeg stræber efter at have den "nye" hjem rumme mine klienters liv uden at bekæmpe dem ved hvert knudepunkt. Derfor vil jeg prale desto mere glæde af mine svagheder, således at Kristi magt kan være over mig .. "Jeg var ked af, at de ikke rigtig vise den side af os, fordi jeg følte, at en masse af de andre racers .. <br><br>Jeg ser ikke meget tv, og har ikke set en fuld episode af The Dog Whisperer. Men en masse tid og penge kunne være brugt, før det er overstået. Jeg kan ikke tale for designerne sandheden skal frem, ville jeg være bange for at være alene i et rum med dem. <br><br>Ugen af ​​6. Juni begyndte han at sætte farten ned på at spise. Jeg troede ham ikke, men jeg stadig kiggede tilbage over mit spørgsmål sørge for jeg ikke ville se dum! (Uanset om jeg gjorde eller ikke at for dig at beslutte) Brokaw stillede det første spørgsmål, så kom til mig. <br><br>Den øverste ned, kontrolleret kommunikation model, skinner er død og borte .. "Hver eneste morgen jeg går mine tvillinger i skole, det er noget, vi gør sammen, og jeg elsker at sende dem ud på deres dag.". Ring til politiet på dem hver chance du får. <br><br>Men jeg var sikker på om mig, om alt, sikrere end han nogensinde kunne være, sikker på mit liv og sikker død, jeg havde ventet på mig. Gud altid tænker med barmhjertighed: Glem ikke dette. The Singles: To af de [http://www.viking-stamps.dk/Email/emotion.asp Ray Ban Aviator] kandidater til registrering og sang af året deler en fællesnævner: Pharrell Williams..<ul>
The usage of the term "force field" in chemistry and computational biology differs from the standard usage in [[physics]]. In chemistry it is a system of potential energy functions rather than the [[gradient]] of a [[scalar potential]], as defined [[force field (physics)|in physics]].
 
  <li>[http://trendhr.egloos.com/565451/ http://trendhr.egloos.com/565451/]</li>
 
  <li>[http://www.pierre-le-jeune.info/spip.php?article3 http://www.pierre-le-jeune.info/spip.php?article3]</li>
 
  <li>[http://bbs.ewche.com/forum.php?mod=viewthread&tid=612332 http://bbs.ewche.com/forum.php?mod=viewthread&tid=612332]</li>
 
  <li>[http://enseignement-lsf.com/spip.php?article64#forum24189074 http://enseignement-lsf.com/spip.php?article64#forum24189074]</li>
 
  <li>[http://bbs.zzabc.cn/forum.php?mod=viewthread&tid=1352807&fromuid=16592 http://bbs.zzabc.cn/forum.php?mod=viewthread&tid=1352807&fromuid=16592]</li>
 
</ul>


== Nike Air Force 1 Danmark Dag ==
==Functional form==
{{further2|[[Molecular mechanics]]}}
[[File:MM PEF.png|thumb|right|Molecular mechanics potential energy function with continuum solvent.]]
The basic functional form of a force field encapsulates both [[chemical bond|bonded]] terms relating to atoms that are linked by [[covalent bond]]s, and nonbonded (also called "noncovalent") terms describing the long-range [[electrostatic]] and [[van der Waals force|van der Waals]] forces. The specific decomposition of the terms depends on the force field, but a general form for the total energy in an additive force field can be written as
<math>\ E_{total} = E_{bonded} + E_{nonbonded} </math>
where the components of the covalent and noncovalent contributions are given by the following summations:


Der er mere: De vil også nødt til at regne ud, hvad man [http://www.marinedesign.dk/application/header.php Nike Air Force 1 Danmark] skal gøre ved den automatiske over hele linjen nedskæringer, som formodes at være prisen for kongressen 'fiasko sidste år at skære underskuddet med 1200 milliarder dollar. Og Kongressen bliver nødt til at stemme endnu en gang på en gæld grænse forlængelse. Alt dette vil sandsynligvis ske i en lame duck kongres, der skal forhandle med Obama, der enten vil [http://www.flex-godning.dk/menumachine/core/banner.asp Nike Store] være blevet genvalgt, eller vil selv være på vej ud af døren.. <br><br>Må ikke blive hjemme hele tiden. Sørg for at du går ud, og lad ham nu, at du komme videre med dit liv. Hvis han mener, du bevæger dig på, så han er mere tilbøjelige til at beslutte, at han vil have dig tilbage i sit liv, og du kan opleve, at han kontakter dig for at se, hvad du synes om den mulighed. <br><br>'I ud pause af mytteri, Maharaja Gulab Singh kaldte sit Rådet og sendt et kontingent på 2.000 infanteri, 200 Cavalry og 6 kanoner under kommando af hans søn og arving, Mian Ranbir Singh Ji. Oprørerne der kom ind på Jammu område blev leveret til den britiske regering. Den Maharaja døde før Mytteriet var blevet nedkæmpet. <br><br>Bestemt, de steder (artiklen mapper, blogs, osv) du udgiver eller du planlægger at udgive dine artikler er specielt vigtige. Selvfølgelig vil du forsøge at udgive dine artikler på disse hjemmesider og blogs, der har en høj PR-niveau på de vigtigste søgemaskiner og betydeligt antal af besøgende (for eksempel på Alexa rating data). Dagens Internettet tilbyder en meget bred vifte af muligheder for artikler publikation. <br><br>Dag, selvom det er sandsynligvis sikrere, end det var dengang, er ofte bange for at lade børnene strejfe, som endeløse alarmerende [http://www.flex-godning.dk/menumachine/flexmenu/preview/cache.asp Vibram Fivefingers Kso] tv-programmer overbevise dem om, at børn, der forlader voksent opsyn for selv et øjeblik uundgåeligt vil blive snuppet forældre. Sikkert bare at lade dem grønt ud foran fjernsynet. Det er næsten som om hele fearmongering ting var virkelig om markedsføring.. <br><br>Interval træning er gode træning for at få hurtigere, tabe sig, og øge din kondition hurtigere end end at gøre en stabil, vedvarende cardio indsats. Ved at inkorporere anfald af højere intensitet i dit fitness træningsprogram vil du ikke kun få styrke, men øge dit stofskifte, som til gengæld vil brænde flere kalorier. Der er mange måder at tilføje intervaller i din træning program og en sådan metode er i hvad der kaldes en pyramide træning.. <br><br>Undersøgelsen har været meget udbredt fejlrapporteres, hvilket er virkelig en skam, da det er en meget interessant undersøgelse. Hvis der findes sammenfiltret dækning i dunhammer, you got at få din flippin pind ud og flip måtterne. Udstationering på hende ved 1:16, skriver Palin, at "på intet tidspunkt ikke under valgkampen, og ikke i de seneste interviews har jeg bedt formanden om at producere hans fødselsattest eller antydet, at han ikke var født i USA." Hun forsøger at afskedige hende birther flirt som blot en erkendelse af vælgernes ret til at vide: "Vælgerne har al mulig ret til at bede kandidater til oplysninger, hvis de ønsker [http://www.flex-godning.dk/menumachine/core/banner.asp Nike København] det..<ul>
<math>\ E_{bonded} = E_{bond} + E_{angle} + E_{dihedral}</math>
 
 
  <li>[http://jewel.egloos.com/3850007/ http://jewel.egloos.com/3850007/]</li>
<math>\ E_{nonbonded} = E_{electrostatic} + E_{van der Waals} </math>
 
 
  <li>[http://ilivewebsolutions.com/awus/DEV/index.php/forum/7-health-a-beauty/268783-ray-ban-wayfarer-som-selskabet-arbejder-altid-hen-imod-hvert-mal-satte-de-bor-arbejde-hen-imod-at-gore-deres-vision-til-virkelighed#268783 http://ilivewebsolutions.com/awus/DEV/index.php/forum/7-health-a-beauty/268783-ray-ban-wayfarer-som-selskabet-arbejder-altid-hen-imod-hvert-mal-satte-de-bor-arbejde-hen-imod-at-gore-deres-vision-til-virkelighed#268783]</li>
The bond and angle terms are usually modeled as [[harmonic oscillator]]s in force fields that do not allow bond breaking. A more realistic description of a covalent bond at higher stretching is provided by the more expensive [[Morse potential]].  Just as the potential energy can be written as a quadratic form in the internal coordinates, so it can also be written in terms of generalized forces. The resulting coefficients are termed compliance constants.
 
 
  <li>[http://old.feipai.net:8010/bbs/showtopic.aspx?topicid=3007666&forumpage=1 http://old.feipai.net:8010/bbs/showtopic.aspx?topicid=3007666&forumpage=1]</li>
The functional form for the rest of the bonded terms is highly variable. Proper dihedral potentials are usually included. Additionally, "improper torsional" terms may be added to enforce the planarity of [[aromatic]] rings and other [[conjugated system]]s, and "cross-terms" that describe coupling of different internal variables, such as angles and bond lengths. Some force fields also include explicit terms for [[hydrogen bond]]s.
 
 
  <li>[http://www.sofapaka.com/index.php?option=com_community&view=profile&userid=9206 http://www.sofapaka.com/index.php?option=com_community&view=profile&userid=9206]</li>
The nonbonded terms are most computationally intensive because they include many more interactions per atom. A popular choice is to limit interactions to pairwise energies. The van der Waals term is usually computed with a [[Lennard-Jones potential]] and the electrostatic term with [[Coulomb's law]], although both can be buffered or scaled by a constant factor to account for electronic [[polarizability]] and produce better agreement with experimental observations.
 
 
  <li>[http://www.harrietoneilllaw.com/node/1834829 http://www.harrietoneilllaw.com/node/1834829]</li>
==Parameterization==
 
In addition to the functional form of the potentials, a force field defines a set of parameters for each type of atom. For example, a force field would include distinct parameters for an [[oxygen]] atom in a [[carbonyl]] [[functional group]] and in a [[hydroxyl]] group. The typical parameter set includes values for [[atomic mass]], [[van der Waals radius]], and [[partial charge]] for individual atoms, and equilibrium values of [[bond length]]s, bond angles, and [[dihedral angle]]s for pairs, triplets, and quadruplets of bonded atoms, and values corresponding to the effective [[spring constant]] for each potential. Most current force fields use a "fixed-charge" model by which each atom is assigned a single value for the atomic [[electrostatic charge|charge]] that is not affected by the local [[electrostatic]] environment; proposed developments in next-generation force fields incorporate models for [[polarizability]], in which a particle's charge is influenced by electrostatic interactions with its neighbors. For example, polarizability can be approximated by the introduction of induced dipoles; it can also be represented by [[Drude particle]]s, or massless, charge-carrying virtual sites attached by a springlike [[harmonic oscillator|harmonic potential]] to each polarizable atom. The introduction of polarizability into force fields in common use has been inhibited by the high computational expense associated with calculating the local electrostatic field.
</ul>
 
Although many molecular simulations involve biological [[macromolecule]]s such as [[protein]]s, [[DNA]], and [[RNA]], the parameters for given atom types are generally derived from observations on small [[organic compound|organic]] molecules that are more tractable for experimental studies and quantum calculations. Different force fields can be derived from dissimilar types of experimental data, such as [[enthalpy]] of [[vaporization]] ([[OPLS]]), [[enthalpy]] of [[Sublimation (chemistry)|sublimation]], [[Molecular dipole moment|dipole moments]], or various spectroscopic parameters.
 
Parameter sets and functional forms are defined by force field developers to be [[self-consistent]]. Because the functional forms of the potential terms vary extensively between even closely related force fields (or successive versions of the same force field), the parameters from one force field should never be used in conjunction with the potential from another.
 
==Deficiencies==
All force fields are based on numerous approximations and derived from different types of experimental data.  Therefore they are called ''empirical''. Some existing force fields do not account for electronic [[Dielectric polarization|polarization]] of the environment, an effect that can significantly reduce electrostatic interactions of partial atomic charges.  This problem was addressed by developing “polarizable force fields” <ref name="Ponder">Ponder JW and Case DA. (2003). "Force fields for protein simulations". ''Adv. Prot. Chem.'' '''66''' 27-85.</ref><ref name="warshel">Warshel A, Sharma PK, Kato M and Parson WW (2006). "Modeling Electrostatic Effects in Proteins." ''Biochim. Biophys. Acta'' '''1764''' 1647-1676.</ref> or using macroscopic [[dielectric constant]]. However, application of a single value of [[dielectric constant]] is questionable in the highly heterogeneous environments of  proteins or biological membranes, and the nature of the dielectric depends on the model used.<ref name="Shultz">Schutz CN. and Warshel A. (2001). "What are the dielectric "constants" of proteins and how to validate electrostatic models?". ''Proteins'' '''44''' 400-417.</ref>
 
All types of [[van der Waals force]]s are also strongly environment-dependent, because these forces originate from interactions of induced and “instantaneous” dipoles (see [[Intermolecular force]]).  The original [[Fritz London]] theory of these forces can only be applied in vacuum.  A more general theory of [[van der Waals force]]s in condensed media was developed by A. D. McLachlan in 1963 (this theory includes the original London’s approach as a special case).<ref name="Israelachvili">Israelachvili, J.N. (1992). ''Intermolecular and surface forces.'' Academic Press, San Diego.</ref>  The McLachlan theory predicts that van der Waals attractions in media are weaker than in vacuum and follow the "like dissolves like" rule, which means that different types of atoms interact more weakly than identical types of atoms.<ref name="Leckband">Leckband, D. and Israelachvili, J. (2001). "Intermolecular forces in biology". ''Quart. Rev. Biophys.'' '''34''' 105-267.</ref>  This is in contrast to “combinatorial rules” or Slater-Kirkwood equation applied for development of the classical force fields.  The “combinatorial rules” state that interaction energy of two dissimilar atoms (e.g. C…N) is an average of the interaction energies of corresponding identical atom pairs (i.e. C…C and N…N).  According to McLachlan theory, the interactions of particles in a media can even be completely repulsive, as  observed for liquid [[helium]].<ref name="Israelachvili" />  The conclusions of McLachlan theory are supported by direct measurements of attraction forces between different materials ([[Hamaker constant]]), as explained by [[Jacob Israelachvili]] in his book "Intermolecular and surface forces". It was concluded that "''the interaction between hydrocarbons across water is about 10% of that across vacuum''".<ref name="Israelachvili" /> Such effects are unaccounted in the standard molecular mechanics.
 
Another round of criticism came from practical applications, such as protein structure refinement.  It was noted that [[CASP]] participants did not try to refine their models to avoid "''a central embarrassment of molecular mechanics, namely that energy minimization or molecular dynamics generally leads to a model that is less like the experimental structure''".<ref name="Koehl">Koehl P. and Levitt M. (1999). "A brighter future for protein structure prediction".  ''Nature Struct. Biol.'' '''6''' 108-111.</ref>  Actually, the force fields have been successfully applied for protein structure refinement in different [[X-ray crystallography]] and [[NMR spectroscopy]] applications, especially using program XPLOR.<ref name="Brunger">Brunger AT and Adams PD. (2002). "Molecular dynamics applied to X-ray structure refinement". ''Acc. Chem. Res.'' '''35''' 404-412.</ref>  However, such refinement is driven primarily by a set of experimental constraints, whereas the force fields serve merely to remove interatomic hindrances. The results of calculations are practically the same with rigid sphere potentials implemented in program DYANA <ref name="Guntert">Guntert P. (1998). "Structure calculation of biological macromolecules from NMR data". ''Quart. Rev. Biophys.'' '''31''' 145-237.</ref> (calculations from NMR data), or with programs for crystallographic refinement that do not use any energy functions. The deficiencies of the force fields remain a major bottleneck in [[homology modeling]] of proteins.<ref name="Tramontano">Tramontano A. and Morea V. (2003). "Assessment of homology-based predictions in CASP5". ''Proteins.'' '''53''' 352-368.</ref> Such situation gave rise to development of alternative empirical scoring functions specifically for [[ligand docking]],<ref name="Gohlke">Gohlke H. and Klebe G. (2002). "Approaches to the description and prediction of the binding affinity of small-molecule ligands to macromolecular receptors". ''Angew. Chem. Internat. Ed.'' '''41''' 2644-2676.</ref> [[protein folding]],<ref name="Edgcomb">Edgcomb S.P. and Murphy K.P. (2000). "Structural energetics of protein folding and binding". ''Current Op. Biotechnol.'' '''11''' 62-66.</ref><ref name="Lazaridis">Lazaridis T. and Karplus M. (2000). "Effective energy functions for protein structure prediction". ''Curr. Op. Struct. Biol.'' '''10''' 139-145.</ref><ref name="awml">Levitt M. and Warshel A. (1975). "Computer Simulations of Protein Folding". ''Nature'' '''253''' 694-698.</ref> homology model refinement,<ref name="krieger">Krieger E., Joo K., Lee J., Lee J., Raman S., Thompson J., Tyka M., Baker D. and Karplus K. (2009). "Improving physical realism, stereochemistry, and side-chain accuracy in homology modeling: Four approaches that performed well in CASP8". ''Proteins'' '''77 Suppl 9''' 114-122.</ref> computational [[protein design]],<ref name="Gordon">Gordon DB, Marshall SA, and Mayo SL (1999). "Energy functions for protein design". ''Curr. Op. Struct. Biol.'' '''9''' 509-513.</ref><ref name="Mendes">Mendes J., Guerois R, and Serrano L (2002). "Energy estimation in protein design". ''Curr. Op. Struct. Biol.'' '''12''' 441-446.</ref><ref name="Rohl">Rohl CA, Strauss CEM, Misura KMS, and Baker D. (2004). "Protein structure prediction using Rosetta". ''Meth. Enz.'' '''383''' 66-93.</ref> and modeling of proteins in membranes.<ref name="Lomize1">Lomize AL, Pogozheva ID, Lomize MA, Mosberg HI (2006). "Positioning of proteins in membranes: A computational approach". ''Protein Sci.'' '''15''' 1318-1333.</ref>
 
There is also an opinion that molecular mechanics may operate with energy which is irrelevant to protein folding or ligand binding.<ref name="Lomize2">Lomize A.L., Reibarkh M.Y. and Pogozheva I.D. (2002). "Interatomic potentials and solvation parameters from protein engineering data for buried residues". ''Protein Sci.'' '''11''' 1984-2000.</ref> The parameters of typical force fields reproduce [[enthalpy]] of [[sublimation (physics)|sublimation]], i.e. energy of evaporation of molecular crystals.  However, it was recognized that [[protein folding]] and [[ligand binding]] are thermodynamically very similar to [[crystallization]], or liquid-solid transitions, because all these processes represent “freezing” of mobile molecules in condensed media.<ref name="Murphy">Murphy K.P. and Gill S.J. (1991). "Solid model compounds and the thermodynamics of protein unfolding". ''J. Mol. Biol.'' '''222''' 699-709.</ref><ref name="Shakhnovich">Shakhnovich, E.I. and Finkelstein, A.V. (1989). "Theory of cooperative transitions in protein molecules.  I. Why denaturation of globular proteins is a first-order phase transition". ''Biopolymers'' '''28''' 1667-1680.</ref><ref name="Graziano">Graziano, G., Catanzano, F., Del Vecchio, P., Giancola, C., and Barone, G. (1996). "Thermodynamic stability of globular proteins: a reliable model from small molecule studies". ''Gazetta Chim. Italiana'' '''126''' 559-567.</ref>  Therefore, free energy changes during protein folding or ligand binding are expected to represent a combination of an energy similar to [[heat of fusion]] (energy absorbed during melting of molecular crystals), a [[conformational entropy]] contribution, and [[solvation]] free energy.  The [[heat of fusion]] is significantly smaller than enthalpy of sublimation.<ref name="Israelachvili" /> Hence, the potentials  describing protein folding or ligand binding must be weaker than potentials in molecular mechanics.  Indeed, the energies of [[H-bond]]s in proteins are ~ -1.5 kcal/mol when estimated from [[protein engineering]] or [[alpha helix]] to [[Random coil|coil]] transition data,<ref name="Myers">Myers J.K. and Pace C.N. (1996). "Hydrogen bonding stabilizes globular proteins". ''Biophys. J.'' '''71''' 2033-2039.</ref><ref name="Scholtz">Scholtz J.M., Marqusee S., Baldwin R.L., York E.J., Stewart J.M., Santoro M., and Bolen D.W. (1991). "Calorimetric determination of the enthalpy change for the alpha-helix to coil transition of an alanine peptide in water". ''Proc. Natl. Acad. Sci. USA'' '''88''' 2854-2858.</ref> but the same energies estimated from [[sublimation (physics)|sublimation]] [[enthalpy]] of molecular [[crystal]]s were -4 to -6 kcal/mol.<ref name="Gavezotti">Gavezotti A. and Filippini G. (1994). "Geometry of intermolecular X-H...Y (X,Y=N,O) hydrogen bond and the calibration of empirical hydrogen-bond potentials". ''J. Phys. Chem.'' '''98''' 4831-4837.</ref>  The depths of modified [[Lennard-Jones potential]]s derived from protein engineering data were also smaller than in typical force fields and followed the “like dissolves like” rule, as predicted by McLachlan theory.<ref name="Lomize2" />
 
==Future perspectives==
Molecular mechanics or force field was first introduced apparently independently by Hill and by Westheimer in 1949,{{Citation needed|date=August 2012}} primarily applied to organic chemistry to estimate properties such as strain energies among others. The functional form of the force field, focused in this article applied to biological systems, was established by Lifson in the 1960s.{{Citation needed|date=August 2012}}  For over a half century, force fields have served us well, providing useful insights into and interpretation of biomolecular structure and function.  Undoubtedly, it will continue to be widely used, thanks to its computational efficiency, while its reliability will continue to be improved.  Yet, there are many well-known deficiencies as noted above.  In addition, the number of energy terms used in a given force field cannot be uniquely determined and a highly redundant number of degrees of freedom are typically used.  Consequently, the "parameters" in different force fields can be vastly different.  Of course, the emphasis to incorporate polarization into the standard pair-wise potentials can be very useful; however, there is no unique way of treating polarization in molecular mechanics because it is of quantum mechanical origin.<ref name="XPOL0">Gao, J. (1997). "Toward a Molecular Orbital Derived Empirical Potential for Liquid Simulations". ''J. Phys. Chem. B'' '''101''' 657-663.</ref><ref name="XPOL">Xie, W., and Gao, J. (2007). "Design of a Next Generation Force Field: The X-POL Potential". ''J. Chem. Theory Comput.'' '''3''' 1890-1900.</ref> Furthermore, often we are more interested in the properties derived from the dynamic dependence of the force field itself on molecular fluctuations.
 
One possibility is that the future development of force field ought to move beyond the current molecular mechanics approach, by using quantum mechanics explicitly to construct the force field.  A number of the "polarizable force fields" listed below, such as density fitting and bond-polarization, already included some of the key ingredients towards this goal. The explicit polarization (X-Pol) method appears to have established the fundamental theoretical framework for a quantal force field; the next step is to develop the necessary parameters to achieve more accurate results than classical mechanics can offer.<ref name="XPOL0"/><ref name="XPOL"/>
 
== Popular force fields ==
Different force fields are designed for different purposes.
 
MM2 was developed by [[Norman Allinger]] primarily for conformational analysis of [[hydrocarbon]]s and other small organic molecules. It is designed to reproduce the equilibrium covalent geometry of molecules as precisely as possible.  It implements a large set of parameters that is continuously refined and updated for many different classes of organic compounds (MM3 and MM4).<ref>Allinger, Norman L. (1977). "Conformational analysis. 130. MM2. A hydrocarbon force field utilizing V1 and V2 torsional terms". ''J. Am. Chem. Soc.'' '''99''' 8127-8134.</ref><ref>[http://europa.chem.uga.edu/allinger/mm2mm3.html MM2 and MM3 home page]</ref><ref>Allinger, N.L., Yuh, Y.H., & Lii, J-H. (1989). "Molecular Mechanics. The MM3 Force Field for Hydrocarbons. 1". ''J. Am. Chem. Soc.'' '''111''' 8551-8565.</ref><ref>
Lii, J-H., & Allinger, N.L. (1989). "Molecular Mechanics. The MM3 Force Field for Hydrocarbons. 2. Vibrational Frequencies and Thermodynamics". ''J. Am. Chem. Soc.'' '''111''' 8566-8575.</ref><ref>
Lii, J-H., & Allinger, N.L. (1989). "Molecular Mechanics. The MM3 Force Field for Hydrocarbons. 3. The van der Waals Potentials and Crystal data for Aliphatic and Aromatic Hydrocarbons". ''J. Am. Chem. Soc.'' '''111''' 8576-8582.</ref>
 
CFF was developed by [[Arieh Warshel]], Lifson and coworkers as a general method for unifying studies of energies, structures and vibration of general molecules and molecular crystals. The CFF program, developed by Levitt and Warshel, is based on the Cartesian representation of all the atoms, and it served as the basis for many subsequent simulation programs.
 
ECEPP was developed specifically for modeling of peptides and proteins. It uses fixed geometries of amino acid residues to simplify the potential energy surface.  Thus, the energy minimization is conducted in the space of protein torsion angles. Both MM2 and ECEPP include potentials for H-bonds and torsion potentials for describing rotations around single bonds. ECEPP/3 was implemented (with some modifications) in [[Internal Coordinate Mechanics]] and FANTOM.<ref name="Schaumann">Schaumann, T., Braun, W. and Wutrich, K. (1990). "The program FANTOM for energy refinement of polypeptides and proteins using a Newton-Raphson minimizer in torsion angle space". ''Biopolymers'' '''29''' 679-694.</ref>
 
[[AMBER]], [[CHARMM]] and [[GROMOS]] have been developed primarily for [[molecular dynamics]] of macromolecules, although they are also commonly applied for energy minimization. Therefore, the coordinates of all atoms are considered as free variables.
 
===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]].
*[[CVFF]] - also broadly used for small molecules and macromolecules.{{Citation needed|date=May 2010}}
*[[COSMOS-NMR]] - hybrid QM/MM forcefield adapted to a variety of inorganic compounds, organic compounds and biological macromolecules, including semi-empirical calculation of atomic charges and NMR properties. COSMOS-NMR is optimized for NMR based structure elucidation and implemented in COSMOS molecular modelling package.<ref>Möllhoff, Margit and Sternberg, Ulrich (2001). "Molecular Mechanics with fluctuating atomic charges - a new force field with semi-empirical charge calculation". ''J. Mol. Model.'' '''7''' 90-102.</ref>
*[[GROMOS]] - a force field that comes as part of the [http://www.igc.ethz.ch/gromos/ GROMOS] (GROningen MOlecular Simulation package), a general-purpose molecular dynamics computer simulation package for the study of biomolecular systems. GROMOS force field (A-version) has been developed for application to aqueous or apolar solutions of proteins, nucleotides and sugars. However, a gas phase version (B-version) for simulation of isolated molecules is also available.
*[[OPLS]] (Optimized Potential for Liquid Simulations) (variations include OPLS-AA, OPLS-UA, OPLS-2001, OPLS-2005) - developed by [http://zarbi.chem.yale.edu/ William L. Jorgensen] at the Yale University Department of Chemistry.
*[http://laetro.usc.edu/programs/#molaris ENZYMIX] – a general polarizable force field for modeling chemical reactions in biological molecules. This force field is implemented with the empirical valence bond (EVB) method and is also combined with the semimacroscopic PDLD approach in the program in the [http://futura.usc.edu/programs/index.html#molaris MOLARIS] package.
*[http://msdlocal.ebi.ac.uk/docs/mmrefs.html#ecepp ECEPP/2] - first force field for polypeptide molecules - developed by F.A. Momany, H.A. Scheraga and colleagues.{{Citation needed|date=May 2010}}
*QCFF/PI – A general force field for conjugated molecules.<ref name="warshel2">Warshel, A. (1973). "Quantum Mechanical Consistent Force Field (QCFF/PI) Method: Calculations of Energies, Conformations and Vibronic Interactions of Ground and Excited States of Conjugated Molecules" ''Israel J. Chem.'' '''11''' 709.</ref><ref name="warshel-levitt">Warshel, A. and Levitt, M. (1974). "QCFF/PI: A Program for the Consistent Force Field Evaluation of Equilibrium Geometries and Vibrational Frequencies of Molecules". QCPE 247, Quantum Chemistry Program Exchange, Indiana University.</ref>
*[[Universal Force Field|UFF]] - A general force field with parameters for the full periodic table up to and including the actinoids - developed at Colorado State University.<ref>{{cite doi|10.1021/ja00051a040}}</ref>
'''''Note'''''
* While it is common to refer to the "GROMACS force field", such a force field has never existed. For a period of time, the ffgmx variant of the standard GROMOS forcefield was referred to as the GROMACS forcefield, but this was only to distinguish it from the main GROMOS forcefield.<ref>http://www.gromacs.org/Documentation/Terminology/Force_Fields/GROMOS</ref> This has been confirmed by David van der Spoel, one of GROMACS' chief authors.<ref>http://lists.gromacs.org/pipermail/gmx-users/2008-March/032661.html</ref>
 
===Second-generation force fields===
*[[Consistent Force Field|CFF]] (Consistent Force Field) - a family of forcefields adapted to a broad variety of organic compounds, includes force fields for polymers, metals, etc.
*[[Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies|COMPASS]] (Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies) - developed by H. Sun at Molecular Simulations Inc., parameterized for a variety of molecules in the condensed phase, now available through [[Accelrys]].<ref>{{cite doi|10.1021/jp980939v}}</ref>
*[[Merck Molecular Force Field|MMFF]] (Merck Molecular Force Field)- developed at Merck, for a broad range of molecules.
*MM2, MM3, MM4 - developed by [[Norman Allinger]], parametrized for a broad range of molecules.
*[http://www.exorga.com/ QVBMM] - developed by Vernon G. S. Box, parameterized for all biomolecules and a broad range of organic molecules, and implemented in StruMM3D (STR3DI32).
 
===Polarizable force fields based on electronic structural theory===
* [[X-Pol: the Explicit Polarization Theory]][http://www.x-pol.org]<ref name="XPOL0"/><ref name="XPOL"/> - a fragment-based electronic structure method introduced by Jiali Gao [http://jialigao.org] at the University of Minnesota, which can be used at any level of theory—ab initio Hartree-Fock (HF), semiempirical molecular orbital theory, correlated wave function theory, or Kohn-Sham (KS) density functional theory (DFT). It is capable of performing more than 3200 steps (3.2 ps) of MD simulations of a fully solvated protein in water with periodic boundary conditions, consisting of about 15000 atoms and 30000 basis functions on a single processor in 24 hours in 2008, with a full quantum mechanical representation of the entire system.<ref name="XPOLbpti">Xie, W., Orozco, M., Truhlar, D. G., and Gao, J. (2009) "X-Pol Potential: An Electronic Structure-Based Force Field for Molecular Dynamics Simulation of a Solvated Protein in Water". ''J. Chem. Theory Comput.'' '''5''' 459-467.</ref> Note that the first MD simulation of a protein by McCammon, Gelin, and Karplus in 1977 lasted 8.8 ps using a united-atom force field without solvent.<ref name="MDProtein">McCammon, J. A., Gelin, B. R., and Karplus, M. (1977) "Dynamics of folded proteins". ''Nature'' '''267''' 585-590.</ref>
 
===Polarizable force field based on induced dipole===
* CFF/ind and ENZYMIX – The first polarizable force field <ref name="awml3">Warshel A. and Levitt M. (1976) "Theoretical Studies of Enzymatic Reactions: Dielectric Electrostatic and Steric Stabilization of the Carbonium Ion in the Reaction of Lysozyme". ''J. Mol. Biol.'' '''103''' 227-249.</ref> which has subsequently been used in many applications to biological systems.<ref name="warshel"/>
* [http://iqc.udg.es/~marcel/manuals/drf90/ DRF90] developed by P. Th. van Duijnen and coworkers.
* PIPF – The polarizable intermolecular potential for fluids is an induced point-dipole force field for organic liquids and biopolymers. The molecular polarization is based on Thole's interacting dipole (TID) model and was developed by Jiali Gao [http://jialigao.org] at the University of Minnesota.<ref name="pipf">Gao, J., Habibollahzadeh, D., and Shao, L. (1995) "A Polarizable Intermolecular Potential Functions for Simulations of Liquid Alcohols". ''J. Phys. Chem.'' '''99''' 16460-16467.</ref><ref name="pipf2">Xie, W., Pu, J., MacKerell, A. D., Jr., and Gao, J. (2007) "Development of a Polarizable Intermolecular Potential Function (PIPF) for Liquid Amides and Alkanes". ''J. Chem. Theory Comput.'' '''3''' 1878-1889.</ref>
 
===Polarizable force fields based on point charges===
* PFF (Polarizable Force Field) developed by Richard A. Friesner and coworkers.{{Citation needed|date=May 2010}}
* SP-basis [[Chemical Potential Equalization]] (CPE) approach developed by R. Chelli and P. Procacci.{{Citation needed|date=May 2010}}
* CHARMM polarizable force field developed by S. Patel (University of Delaware) and C. L. Brooks III (University of Michigan).<ref name="CPFF0">Patel, S., Brooks, C.L., III (2004) "CHARMM fluctuating charge force field for proteins: I parameterization and application to bulk organic liquid simulations." ''J. Comput. Chem.'' '''25''' 1-16.</ref><ref name="CPFF1">Patel, S., MacKerell, A.D., Jr., Brooks, C.L., III (2004) "CHARMM fluctuating charge force field for proteins: II Protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model." ''J. Comput. Chem.'' '''25''' 1504-1514.</ref>
* AMBER polarizable force field developed by Jim Caldwell and coworkers.{{Citation needed|date=May 2010}}
* CHARMM polarizable force field based on the classical Drude oscillator developed by A. MacKerell (University of Maryland, Baltimore) and B. Roux (University of Chicago).<ref name="drude1">Anisimov, V.M., Lamoureux, G., Vorobyov, I.V., Huang, N., Roux, B. and MacKerell, A.D., Jr. (2005) "Determination of Electrostatic Parameters for a Polarizable Force Field Based on the Classical Drude Oscillator". ''J. Chem. Theory Comput.'' '''1''' 153-168.</ref><ref name="drude2">Yu, H., Whitfield, T.W., Harder, E., Lamoureux, G., Vorobyov, I., Anisimov, V. M., MacKerell, A.D., Jr., and Roux, B. (2010) "Simulating Monovalent and Divalent Ions in Aqueous Solution Using a Drude Polarizable Force Field". ''J. Chem. Theory Comput.'' '''6''' 774–786.</ref>
 
===Polarizable force fields based on distributed multipoles===
* The [http://www.lct.jussieu.fr/pagesperso/jpp/ SIBFA (Sum of Interactions Between Fragments Ab initio computed)] force field <ref name="SIBGEM">Gresh, N., Cisneros, G. A., Darden, T. A. and Piquemal, J-P. (2007) "Anisotropic, polarizable molecular mechanics studies of inter-, intra-molecular interactions, and ligand-macromolecule complexes. A bottom-up strategy." ''J. Chem. Theory. Comput.'' '''3''' 1960.</ref> for small molecules and flexible proteins, developed by Nohad Gresh (Paris V, René Descartes University) and Jean-Philip Piquemal (Paris VI, Pierre & Marie Curie University). SIBFA is a molecular mechanics procedure formulated and calibrated on the basis of ab initio supermolecule computations. Its purpose is to enable the simultaneous and reliable computations of both intermolecular and conformational energies governing the binding specificities of biologically and pharmacologically relevant molecules. This procedure enables an accurate treatment of transition metals. The inclusion of a ligand field contribution allows computations on "open-shell" metalloproteins.
* [http://dasher.wustl.edu/tinker/ AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications) ]  force field developed by Pengyu Ren (University of Texas at Austin) and Jay W. Ponder (Washington University).
* [http://www-stone.ch.cam.ac.uk/programs.html ORIENT procedure] developed by Anthony J. Stone (Cambridge University) and coworkers.
* Non-Empirical Molecular Orbital (NEMO) procedure developed by Gunnar Karlström and coworkers at Lund University (Sweden)<ref name="NEMO">O. Engkvist, P.-O. Åstrand, and G. Karlström. "Accurate intermolecular potentials obtained from molecular wave functions: Bridging the gap between quantum chemistry and molecular simulations". Chem. Rev., 100:4087-4108, 2000.</ref>
 
===Polarizable force fields based on density===
* Gaussian Electrostatic Model (GEM)<ref name="SIBGEM"/><ref name="GEM0">Piquemal, J-P, Cisneros, G.A., Reinhardt, P., Gresh, N. and Darden, T.A. (2006). "Towards a Force Field based on Density Fitting".
''J. Chem. Phys.'' '''124''' 104101.</ref><ref name="GEM2">Cisneros, G.A., Piquemal, J-P and Darden, T.A. (2006). "Generalization of the Gaussian Electrostatic Model: extension to arbitrary angular momentum, distributed multipoles and speedup with reciprocal space methods". ''J. Chem. Phys.'' '''125''' 184101.</ref> - a polarizable force field based on Density Fitting developed by Thomas A. Darden and G. Andrés Cisneros at NIEHS; and Jean-Philip Piquemal (Paris VI University).
* Polarizable procedure based on the Kim-Gordon approach developed by Jürg Hutter and coworkers (University of Zürich){{Citation needed|date=May 2010}}
 
===Polarizable force fields based on Bond Polarization Theory (BPT)===
*[[COSMOS-NMR]] (Computer Simulation of Molecular Structure) - developed by Ulrich Sternberg and coworkers. Hybrid QM/MM force field enables explicit quantum-mechanical calculation of electrostatic properties using localized bond orbitals with fast BPT formalism.<ref>Sternberg, Ulrich,  Koch, Frank-Thomas and Möllhoff, Margit (1994). "New Approach to the Semiempirical Calculation of Atomic Charges for Polypeptides and Large Molecular Systems" ''J. Comp. Chem.'' '''15''' 524-531.</ref> Atomic charge fluctuation is possible in each molecular dynamics step.
 
===Reactive force fields===
*[[ReaxFF]] - reactive force field developed by Adri van Duin, [[William Andrew Goddard III|William Goddard]] and coworkers. It is fast, transferable and is the computational method of choice for atomistic-scale dynamical simulations of chemical reactions.<ref>{{cite doi|10.1021/jp004368u}}</ref> Parallelized ReaxFF allows reactive simulations on >>1000,000 atoms.
*EVB (empirical valence bond) - this reactive force field, introduced by Warshel and coworkers, is probably the most reliable and physically consistent way of using force fields in modeling chemical reactions in different environments.{{According to whom|date=May 2010}} The EVB facilitates calculations of actual activation free energies in condensed phases and in enzymes.
*[[RWFF]] - reactive force field for water developed by Detlef W. M. Hofmann, Liudmila N. Kuleshova and Bruno D'Aguanno. It is very fast{{Quantify|date=May 2010}}, reproduces the experimental data of neutron scattering accurately, and allows the simulation of bond formation/breaking of water and acids.<ref>{{cite doi|10.1016/j.cplett.2007.09.063}}</ref>
 
===Coarse-grained force fields===
*[[VAMM]] (Virtual atom molecular mechanics) - a coarse-grained force field developed by Korkut and Hendrickson for molecular mechanics calculations such as large scale conformational transitions based on the virtual interactions of C-alpha atoms. It is a knowledge based force field and formulated to capture features dependent on secondary structure and on residue-specific contact information in proteins.<ref>A. Korkut, W.A. Hendrickson (2009). "A force field for virtual atom molecular mechanics of proteins." ''Proc. Natl. Acad. Sci. USA'' '''106''' 15667–15672.</ref>
 
*[[MARTINI]] - a coarse-grained force field developed by Marrink and coworkers at the [[University of Groningen]], initially developed for molecular dynamics simulations of lipids,<ref>S.J. Marrink, H.J. Risselada, S. Yefimov,D.P. Tieleman,A.H. de Vries (2007). "The MARTINI force field: coarse grained model for biomolecular simulations." ''J. Phys. Chem. B.'' '''111''' 7812-7824.</ref> later extended to various other molecules. The force field applies a mapping of four heavy atoms to one CG interaction site and is parametrized with the aim of reproducing thermodynamic properties.
 
===Water models===
{{Main|water model}}
 
The set of parameters used to model water or aqueous solutions (basically a force field for water) is called a [[water model]]. Water has attracted a great deal of attention due to its unusual properties and its importance as a solvent. Many water models have been proposed; some examples are TIP3P, TIP4P, SPC, [[Flexible SPC water model|Flexible SPC]], and ST2.
 
===Post-translational Modifications and Unnatural Amino Acids===
*[http://selene.princeton.edu/FFPTM/ Forcefield_PTM] - An AMBER-based forcefield and webtool for modeling common post-translational modifications of amino acids in proteins developed by Chris Floudas and coworkers. It utilizes the ff03 charge model.<ref>{{cite doi|10.1021/ct400556v}}</ref>
*[http://selene.princeton.edu/FFNCAA/ Forcefield_NCAA] - An AMBER-based forcefield and webtool for modeling common unnatural amino acids in proteins in condensed-phase simulations using the ff03 charge model.<ref>{{cite doi|10.1021/sb400168u}}</ref>
 
===Other===
*[[VALBOND]] - a function for angle bending that is based on [[valence bond theory]] and works for large angular distortions, [[hypervalent molecule]]s, and [[complex (chemistry)|transition metal complexes]]. It can be incorporated into other force fields such as CHARMM and UFF.
 
== See also ==
{{columns-list|2|
* [[Force field implementation]]
* [[Molecular dynamics]]
* [[Molecular design software]]
* [[Molecular modelling]]
* [[List of software for molecular mechanics modeling|Software for molecular mechanics modeling]]
* [[Statistical potential]]
}}
 
==References==
{{Reflist|30em}}
 
==Further reading==
*{{cite book |last=Israelachvili |first=J. N. |year=1992 |title=Intermolecular and surface forces |publisher=Academic Press |location=San Diego |isbn=0-12-375181-0 |ref=harv }}
*{{cite book |last=Schlick |first=T. |year=2002 |title=Molecular Modeling and Simulation: An Interdisciplinary Guide |series=Interdisciplinary Applied Mathematics: Mathematical Biology |publisher=Springer-Verlag |location=New York |isbn=0-387-95404-X |ref=harv }}
*{{cite book |last=Warshel |first=A. |year=1991 |title=Computer Modeling of Chemical Reactions in Enzymes and Solutions |publisher=John Wiley & Sons |location=New York |isbn=0-471-53395-5 |ref=harv }}
 
{{DEFAULTSORT:Force Field (Chemistry)}}
[[Category:Force fields| ]]
[[Category:Intermolecular forces]]
[[Category:Molecular physics]]
[[Category:Molecular modelling]]

Revision as of 12:01, 30 July 2013

I'm Fernando (21) from Seltjarnarnes, Iceland.
I'm learning Norwegian literature at a local college and I'm just about to graduate.
I have a part time job in a the office.

my site; wellness [continue reading this..]

A force field is used to minimize the bond stretching energy of this ethane molecule.

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). Force field functions and parameter sets are derived from both experimental work and high-level quantum mechanical calculations. "All-atom" force fields provide parameters for every type of atom in a system, including hydrogen, while "united-atom" force fields treat the hydrogen and carbon atoms in each terminal methyl and each methylene bridge as a single interaction center. "Coarse-grained" force fields, which are frequently used in long-time simulations of proteins, provide even more crude representations for increased computational efficiency.

The usage of the term "force field" in chemistry and computational biology differs from the standard usage in physics. In chemistry it is a system of potential energy functions rather than the gradient of a scalar potential, as defined in physics.

Functional form

Template:Further2

Molecular mechanics potential energy function with continuum solvent.

The basic functional form of a force field encapsulates both bonded terms relating to atoms that are linked by covalent bonds, and nonbonded (also called "noncovalent") terms describing the long-range electrostatic and van der Waals forces. The specific decomposition of the terms depends on the force field, but a general form for the total energy in an additive force field can be written as where the components of the covalent and noncovalent contributions are given by the following summations:

The bond and angle terms are usually modeled as harmonic oscillators in force fields that do not allow bond breaking. A more realistic description of a covalent bond at higher stretching is provided by the more expensive Morse potential. Just as the potential energy can be written as a quadratic form in the internal coordinates, so it can also be written in terms of generalized forces. The resulting coefficients are termed compliance constants.

The functional form for the rest of the bonded terms is highly variable. Proper dihedral potentials are usually included. Additionally, "improper torsional" terms may be added to enforce the planarity of aromatic rings and other conjugated systems, and "cross-terms" that describe coupling of different internal variables, such as angles and bond lengths. Some force fields also include explicit terms for hydrogen bonds.

The nonbonded terms are most computationally intensive because they include many more interactions per atom. A popular choice is to limit interactions to pairwise energies. The van der Waals term is usually computed with a Lennard-Jones potential and the electrostatic term with Coulomb's law, although both can be buffered or scaled by a constant factor to account for electronic polarizability and produce better agreement with experimental observations.

Parameterization

In addition to the functional form of the potentials, a force field defines a set of parameters for each type of atom. For example, a force field would include distinct parameters for an oxygen atom in a carbonyl functional group and in a hydroxyl group. The typical parameter set includes values for atomic mass, van der Waals radius, and partial charge for individual atoms, and equilibrium values of bond 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. Most current force fields use a "fixed-charge" model by which each atom is assigned a single value for the atomic charge that is not affected by the local electrostatic environment; proposed developments in next-generation force fields incorporate models for polarizability, in which a particle's charge is influenced by electrostatic interactions with its neighbors. For example, polarizability can be approximated by the introduction of induced dipoles; it can also be represented by Drude particles, or massless, charge-carrying virtual sites attached by a springlike harmonic potential to each polarizable atom. The introduction of polarizability into force fields in common use has been inhibited by the high computational expense associated with calculating the local electrostatic field.

Although many molecular simulations involve biological macromolecules such as proteins, DNA, and RNA, the parameters for given atom types are generally derived from observations on small organic molecules that are more tractable for experimental studies and quantum calculations. Different force fields can be derived from dissimilar types of experimental data, such as enthalpy of vaporization (OPLS), enthalpy of sublimation, dipole moments, or various spectroscopic parameters.

Parameter sets and functional forms are defined by force field developers to be self-consistent. Because the functional forms of the potential terms vary extensively between even closely related force fields (or successive versions of the same force field), the parameters from one force field should never be used in conjunction with the potential from another.

Deficiencies

All force fields are based on numerous approximations and derived from different types of experimental data. Therefore they are called empirical. Some existing force fields do not account for electronic polarization of the environment, an effect that can significantly reduce electrostatic interactions of partial atomic charges. This problem was addressed by developing “polarizable force fields” [1][2] or using macroscopic dielectric constant. However, application of a single value of dielectric constant is questionable in the highly heterogeneous environments of proteins or biological membranes, and the nature of the dielectric depends on the model used.[3]

All types of van der Waals forces are also strongly environment-dependent, because these forces originate from interactions of induced and “instantaneous” dipoles (see Intermolecular force). The original Fritz London theory of these forces can only be applied in vacuum. A more general theory of van der Waals forces in condensed media was developed by A. D. McLachlan in 1963 (this theory includes the original London’s approach as a special case).[4] The McLachlan theory predicts that van der Waals attractions in media are weaker than in vacuum and follow the "like dissolves like" rule, which means that different types of atoms interact more weakly than identical types of atoms.[5] This is in contrast to “combinatorial rules” or Slater-Kirkwood equation applied for development of the classical force fields. The “combinatorial rules” state that interaction energy of two dissimilar atoms (e.g. C…N) is an average of the interaction energies of corresponding identical atom pairs (i.e. C…C and N…N). According to McLachlan theory, the interactions of particles in a media can even be completely repulsive, as observed for liquid helium.[4] The conclusions of McLachlan theory are supported by direct measurements of attraction forces between different materials (Hamaker constant), as explained by Jacob Israelachvili in his book "Intermolecular and surface forces". It was concluded that "the interaction between hydrocarbons across water is about 10% of that across vacuum".[4] Such effects are unaccounted in the standard molecular mechanics.

Another round of criticism came from practical applications, such as protein structure refinement. It was noted that CASP participants did not try to refine their models to avoid "a central embarrassment of molecular mechanics, namely that energy minimization or molecular dynamics generally leads to a model that is less like the experimental structure".[6] Actually, the force fields have been successfully applied for protein structure refinement in different X-ray crystallography and NMR spectroscopy applications, especially using program XPLOR.[7] However, such refinement is driven primarily by a set of experimental constraints, whereas the force fields serve merely to remove interatomic hindrances. The results of calculations are practically the same with rigid sphere potentials implemented in program DYANA [8] (calculations from NMR data), or with programs for crystallographic refinement that do not use any energy functions. The deficiencies of the force fields remain a major bottleneck in homology modeling of proteins.[9] Such situation gave rise to development of alternative empirical scoring functions specifically for ligand docking,[10] protein folding,[11][12][13] homology model refinement,[14] computational protein design,[15][16][17] and modeling of proteins in membranes.[18]

There is also an opinion that molecular mechanics may operate with energy which is irrelevant to protein folding or ligand binding.[19] The parameters of typical force fields reproduce enthalpy of sublimation, i.e. energy of evaporation of molecular crystals. However, it was recognized that protein folding and ligand binding are thermodynamically very similar to crystallization, or liquid-solid transitions, because all these processes represent “freezing” of mobile molecules in condensed media.[20][21][22] Therefore, free energy changes during protein folding or ligand binding are expected to represent a combination of an energy similar to heat of fusion (energy absorbed during melting of molecular crystals), a conformational entropy contribution, and solvation free energy. The heat of fusion is significantly smaller than enthalpy of sublimation.[4] Hence, the potentials describing protein folding or ligand binding must be weaker than potentials in molecular mechanics. Indeed, the energies of H-bonds in proteins are ~ -1.5 kcal/mol when estimated from protein engineering or alpha helix to coil transition data,[23][24] but the same energies estimated from sublimation enthalpy of molecular crystals were -4 to -6 kcal/mol.[25] The depths of modified Lennard-Jones potentials derived from protein engineering data were also smaller than in typical force fields and followed the “like dissolves like” rule, as predicted by McLachlan theory.[19]

Future perspectives

Molecular mechanics or force field was first introduced apparently independently by Hill and by Westheimer in 1949,Potter or Ceramic Artist Truman Bedell from Rexton, has interests which include ceramics, best property developers in singapore developers in singapore and scrabble. Was especially enthused after visiting Alejandro de Humboldt National Park. primarily applied to organic chemistry to estimate properties such as strain energies among others. The functional form of the force field, focused in this article applied to biological systems, was established by Lifson in the 1960s.Potter or Ceramic Artist Truman Bedell from Rexton, has interests which include ceramics, best property developers in singapore developers in singapore and scrabble. Was especially enthused after visiting Alejandro de Humboldt National Park. For over a half century, force fields have served us well, providing useful insights into and interpretation of biomolecular structure and function. Undoubtedly, it will continue to be widely used, thanks to its computational efficiency, while its reliability will continue to be improved. Yet, there are many well-known deficiencies as noted above. In addition, the number of energy terms used in a given force field cannot be uniquely determined and a highly redundant number of degrees of freedom are typically used. Consequently, the "parameters" in different force fields can be vastly different. Of course, the emphasis to incorporate polarization into the standard pair-wise potentials can be very useful; however, there is no unique way of treating polarization in molecular mechanics because it is of quantum mechanical origin.[26][27] Furthermore, often we are more interested in the properties derived from the dynamic dependence of the force field itself on molecular fluctuations.

One possibility is that the future development of force field ought to move beyond the current molecular mechanics approach, by using quantum mechanics explicitly to construct the force field. A number of the "polarizable force fields" listed below, such as density fitting and bond-polarization, already included some of the key ingredients towards this goal. The explicit polarization (X-Pol) method appears to have established the fundamental theoretical framework for a quantal force field; the next step is to develop the necessary parameters to achieve more accurate results than classical mechanics can offer.[26][27]

Popular force fields

Different force fields are designed for different purposes.

MM2 was developed by Norman Allinger primarily for conformational analysis of hydrocarbons and other small organic molecules. It is designed to reproduce the equilibrium covalent geometry of molecules as precisely as possible. It implements a large set of parameters that is continuously refined and updated for many different classes of organic compounds (MM3 and MM4).[28][29][30][31][32]

CFF was developed by Arieh Warshel, Lifson and coworkers as a general method for unifying studies of energies, structures and vibration of general molecules and molecular crystals. The CFF program, developed by Levitt and Warshel, is based on the Cartesian representation of all the atoms, and it served as the basis for many subsequent simulation programs.

ECEPP was developed specifically for modeling of peptides and proteins. It uses fixed geometries of amino acid residues to simplify the potential energy surface. Thus, the energy minimization is conducted in the space of protein torsion angles. Both MM2 and ECEPP include potentials for H-bonds and torsion potentials for describing rotations around single bonds. ECEPP/3 was implemented (with some modifications) in Internal Coordinate Mechanics and FANTOM.[33]

AMBER, CHARMM and GROMOS have been developed primarily for molecular dynamics of macromolecules, although they are also commonly applied for energy minimization. Therefore, the coordinates of all atoms are considered as free variables.

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.
  • CVFF - also broadly used for small molecules and macromolecules.Potter or Ceramic Artist Truman Bedell from Rexton, has interests which include ceramics, best property developers in singapore developers in singapore and scrabble. Was especially enthused after visiting Alejandro de Humboldt National Park.
  • COSMOS-NMR - hybrid QM/MM forcefield adapted to a variety of inorganic compounds, organic compounds and biological macromolecules, including semi-empirical calculation of atomic charges and NMR properties. COSMOS-NMR is optimized for NMR based structure elucidation and implemented in COSMOS molecular modelling package.[34]
  • GROMOS - a force field that comes as part of the GROMOS (GROningen MOlecular Simulation package), a general-purpose molecular dynamics computer simulation package for the study of biomolecular systems. GROMOS force field (A-version) has been developed for application to aqueous or apolar solutions of proteins, nucleotides and sugars. However, a gas phase version (B-version) for simulation of isolated molecules is also available.
  • OPLS (Optimized Potential for Liquid Simulations) (variations include OPLS-AA, OPLS-UA, OPLS-2001, OPLS-2005) - developed by William L. Jorgensen at the Yale University Department of Chemistry.
  • ENZYMIX – a general polarizable force field for modeling chemical reactions in biological molecules. This force field is implemented with the empirical valence bond (EVB) method and is also combined with the semimacroscopic PDLD approach in the program in the MOLARIS package.
  • ECEPP/2 - first force field for polypeptide molecules - developed by F.A. Momany, H.A. Scheraga and colleagues.Potter or Ceramic Artist Truman Bedell from Rexton, has interests which include ceramics, best property developers in singapore developers in singapore and scrabble. Was especially enthused after visiting Alejandro de Humboldt National Park.
  • QCFF/PI – A general force field for conjugated molecules.[35][36]
  • UFF - A general force field with parameters for the full periodic table up to and including the actinoids - developed at Colorado State University.[37]

Note

  • While it is common to refer to the "GROMACS force field", such a force field has never existed. For a period of time, the ffgmx variant of the standard GROMOS forcefield was referred to as the GROMACS forcefield, but this was only to distinguish it from the main GROMOS forcefield.[38] This has been confirmed by David van der Spoel, one of GROMACS' chief authors.[39]

Second-generation force fields

  • CFF (Consistent Force Field) - a family of forcefields adapted to a broad variety of organic compounds, includes force fields for polymers, metals, etc.
  • COMPASS (Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies) - developed by H. Sun at Molecular Simulations Inc., parameterized for a variety of molecules in the condensed phase, now available through Accelrys.[40]
  • MMFF (Merck Molecular Force Field)- developed at Merck, for a broad range of molecules.
  • MM2, MM3, MM4 - developed by Norman Allinger, parametrized for a broad range of molecules.
  • QVBMM - developed by Vernon G. S. Box, parameterized for all biomolecules and a broad range of organic molecules, and implemented in StruMM3D (STR3DI32).

Polarizable force fields based on electronic structural theory

  • X-Pol: the Explicit Polarization Theory[1][26][27] - a fragment-based electronic structure method introduced by Jiali Gao [2] at the University of Minnesota, which can be used at any level of theory—ab initio Hartree-Fock (HF), semiempirical molecular orbital theory, correlated wave function theory, or Kohn-Sham (KS) density functional theory (DFT). It is capable of performing more than 3200 steps (3.2 ps) of MD simulations of a fully solvated protein in water with periodic boundary conditions, consisting of about 15000 atoms and 30000 basis functions on a single processor in 24 hours in 2008, with a full quantum mechanical representation of the entire system.[41] Note that the first MD simulation of a protein by McCammon, Gelin, and Karplus in 1977 lasted 8.8 ps using a united-atom force field without solvent.[42]

Polarizable force field based on induced dipole

  • CFF/ind and ENZYMIX – The first polarizable force field [43] which has subsequently been used in many applications to biological systems.[2]
  • DRF90 developed by P. Th. van Duijnen and coworkers.
  • PIPF – The polarizable intermolecular potential for fluids is an induced point-dipole force field for organic liquids and biopolymers. The molecular polarization is based on Thole's interacting dipole (TID) model and was developed by Jiali Gao [3] at the University of Minnesota.[44][45]

Polarizable force fields based on point charges

  • PFF (Polarizable Force Field) developed by Richard A. Friesner and coworkers.Potter or Ceramic Artist Truman Bedell from Rexton, has interests which include ceramics, best property developers in singapore developers in singapore and scrabble. Was especially enthused after visiting Alejandro de Humboldt National Park.
  • SP-basis Chemical Potential Equalization (CPE) approach developed by R. Chelli and P. Procacci.Potter or Ceramic Artist Truman Bedell from Rexton, has interests which include ceramics, best property developers in singapore developers in singapore and scrabble. Was especially enthused after visiting Alejandro de Humboldt National Park.
  • CHARMM polarizable force field developed by S. Patel (University of Delaware) and C. L. Brooks III (University of Michigan).[46][47]
  • AMBER polarizable force field developed by Jim Caldwell and coworkers.Potter or Ceramic Artist Truman Bedell from Rexton, has interests which include ceramics, best property developers in singapore developers in singapore and scrabble. Was especially enthused after visiting Alejandro de Humboldt National Park.
  • CHARMM polarizable force field based on the classical Drude oscillator developed by A. MacKerell (University of Maryland, Baltimore) and B. Roux (University of Chicago).[48][49]

Polarizable force fields based on distributed multipoles

  • The SIBFA (Sum of Interactions Between Fragments Ab initio computed) force field [50] for small molecules and flexible proteins, developed by Nohad Gresh (Paris V, René Descartes University) and Jean-Philip Piquemal (Paris VI, Pierre & Marie Curie University). SIBFA is a molecular mechanics procedure formulated and calibrated on the basis of ab initio supermolecule computations. Its purpose is to enable the simultaneous and reliable computations of both intermolecular and conformational energies governing the binding specificities of biologically and pharmacologically relevant molecules. This procedure enables an accurate treatment of transition metals. The inclusion of a ligand field contribution allows computations on "open-shell" metalloproteins.
  • AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications) force field developed by Pengyu Ren (University of Texas at Austin) and Jay W. Ponder (Washington University).
  • ORIENT procedure developed by Anthony J. Stone (Cambridge University) and coworkers.
  • Non-Empirical Molecular Orbital (NEMO) procedure developed by Gunnar Karlström and coworkers at Lund University (Sweden)[51]

Polarizable force fields based on density

  • Gaussian Electrostatic Model (GEM)[50][52][53] - a polarizable force field based on Density Fitting developed by Thomas A. Darden and G. Andrés Cisneros at NIEHS; and Jean-Philip Piquemal (Paris VI University).
  • Polarizable procedure based on the Kim-Gordon approach developed by Jürg Hutter and coworkers (University of Zürich)Potter or Ceramic Artist Truman Bedell from Rexton, has interests which include ceramics, best property developers in singapore developers in singapore and scrabble. Was especially enthused after visiting Alejandro de Humboldt National Park.

Polarizable force fields based on Bond Polarization Theory (BPT)

  • COSMOS-NMR (Computer Simulation of Molecular Structure) - developed by Ulrich Sternberg and coworkers. Hybrid QM/MM force field enables explicit quantum-mechanical calculation of electrostatic properties using localized bond orbitals with fast BPT formalism.[54] Atomic charge fluctuation is possible in each molecular dynamics step.

Reactive force fields

  • ReaxFF - reactive force field developed by Adri van Duin, William Goddard and coworkers. It is fast, transferable and is the computational method of choice for atomistic-scale dynamical simulations of chemical reactions.[55] Parallelized ReaxFF allows reactive simulations on >>1000,000 atoms.
  • EVB (empirical valence bond) - this reactive force field, introduced by Warshel and coworkers, is probably the most reliable and physically consistent way of using force fields in modeling chemical reactions in different environments.Template:According to whom The EVB facilitates calculations of actual activation free energies in condensed phases and in enzymes.
  • RWFF - reactive force field for water developed by Detlef W. M. Hofmann, Liudmila N. Kuleshova and Bruno D'Aguanno. It is very fastTemplate:Quantify, reproduces the experimental data of neutron scattering accurately, and allows the simulation of bond formation/breaking of water and acids.[56]

Coarse-grained force fields

  • VAMM (Virtual atom molecular mechanics) - a coarse-grained force field developed by Korkut and Hendrickson for molecular mechanics calculations such as large scale conformational transitions based on the virtual interactions of C-alpha atoms. It is a knowledge based force field and formulated to capture features dependent on secondary structure and on residue-specific contact information in proteins.[57]
  • MARTINI - a coarse-grained force field developed by Marrink and coworkers at the University of Groningen, initially developed for molecular dynamics simulations of lipids,[58] later extended to various other molecules. The force field applies a mapping of four heavy atoms to one CG interaction site and is parametrized with the aim of reproducing thermodynamic properties.

Water models

Mining Engineer (Excluding Oil ) Truman from Alma, loves to spend time knotting, largest property developers in singapore developers in singapore and stamp collecting. Recently had a family visit to Urnes Stave Church.

The set of parameters used to model water or aqueous solutions (basically a force field for water) is called a water model. Water has attracted a great deal of attention due to its unusual properties and its importance as a solvent. Many water models have been proposed; some examples are TIP3P, TIP4P, SPC, Flexible SPC, and ST2.

Post-translational Modifications and Unnatural Amino Acids

  • Forcefield_PTM - An AMBER-based forcefield and webtool for modeling common post-translational modifications of amino acids in proteins developed by Chris Floudas and coworkers. It utilizes the ff03 charge model.[59]
  • Forcefield_NCAA - An AMBER-based forcefield and webtool for modeling common unnatural amino acids in proteins in condensed-phase simulations using the ff03 charge model.[60]

Other

See also

Template:Columns-list

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.

Further reading

  • 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.

    My blog: http://www.primaboinca.com/view_profile.php?userid=5889534
  • 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.

    My blog: http://www.primaboinca.com/view_profile.php?userid=5889534
  • 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.

    My blog: http://www.primaboinca.com/view_profile.php?userid=5889534
  1. Ponder JW and Case DA. (2003). "Force fields for protein simulations". Adv. Prot. Chem. 66 27-85.
  2. 2.0 2.1 Warshel A, Sharma PK, Kato M and Parson WW (2006). "Modeling Electrostatic Effects in Proteins." Biochim. Biophys. Acta 1764 1647-1676.
  3. Schutz CN. and Warshel A. (2001). "What are the dielectric "constants" of proteins and how to validate electrostatic models?". Proteins 44 400-417.
  4. 4.0 4.1 4.2 4.3 Israelachvili, J.N. (1992). Intermolecular and surface forces. Academic Press, San Diego.
  5. Leckband, D. and Israelachvili, J. (2001). "Intermolecular forces in biology". Quart. Rev. Biophys. 34 105-267.
  6. Koehl P. and Levitt M. (1999). "A brighter future for protein structure prediction". Nature Struct. Biol. 6 108-111.
  7. Brunger AT and Adams PD. (2002). "Molecular dynamics applied to X-ray structure refinement". Acc. Chem. Res. 35 404-412.
  8. Guntert P. (1998). "Structure calculation of biological macromolecules from NMR data". Quart. Rev. Biophys. 31 145-237.
  9. Tramontano A. and Morea V. (2003). "Assessment of homology-based predictions in CASP5". Proteins. 53 352-368.
  10. Gohlke H. and Klebe G. (2002). "Approaches to the description and prediction of the binding affinity of small-molecule ligands to macromolecular receptors". Angew. Chem. Internat. Ed. 41 2644-2676.
  11. Edgcomb S.P. and Murphy K.P. (2000). "Structural energetics of protein folding and binding". Current Op. Biotechnol. 11 62-66.
  12. Lazaridis T. and Karplus M. (2000). "Effective energy functions for protein structure prediction". Curr. Op. Struct. Biol. 10 139-145.
  13. Levitt M. and Warshel A. (1975). "Computer Simulations of Protein Folding". Nature 253 694-698.
  14. Krieger E., Joo K., Lee J., Lee J., Raman S., Thompson J., Tyka M., Baker D. and Karplus K. (2009). "Improving physical realism, stereochemistry, and side-chain accuracy in homology modeling: Four approaches that performed well in CASP8". Proteins 77 Suppl 9 114-122.
  15. Gordon DB, Marshall SA, and Mayo SL (1999). "Energy functions for protein design". Curr. Op. Struct. Biol. 9 509-513.
  16. Mendes J., Guerois R, and Serrano L (2002). "Energy estimation in protein design". Curr. Op. Struct. Biol. 12 441-446.
  17. Rohl CA, Strauss CEM, Misura KMS, and Baker D. (2004). "Protein structure prediction using Rosetta". Meth. Enz. 383 66-93.
  18. Lomize AL, Pogozheva ID, Lomize MA, Mosberg HI (2006). "Positioning of proteins in membranes: A computational approach". Protein Sci. 15 1318-1333.
  19. 19.0 19.1 Lomize A.L., Reibarkh M.Y. and Pogozheva I.D. (2002). "Interatomic potentials and solvation parameters from protein engineering data for buried residues". Protein Sci. 11 1984-2000.
  20. Murphy K.P. and Gill S.J. (1991). "Solid model compounds and the thermodynamics of protein unfolding". J. Mol. Biol. 222 699-709.
  21. Shakhnovich, E.I. and Finkelstein, A.V. (1989). "Theory of cooperative transitions in protein molecules. I. Why denaturation of globular proteins is a first-order phase transition". Biopolymers 28 1667-1680.
  22. Graziano, G., Catanzano, F., Del Vecchio, P., Giancola, C., and Barone, G. (1996). "Thermodynamic stability of globular proteins: a reliable model from small molecule studies". Gazetta Chim. Italiana 126 559-567.
  23. Myers J.K. and Pace C.N. (1996). "Hydrogen bonding stabilizes globular proteins". Biophys. J. 71 2033-2039.
  24. Scholtz J.M., Marqusee S., Baldwin R.L., York E.J., Stewart J.M., Santoro M., and Bolen D.W. (1991). "Calorimetric determination of the enthalpy change for the alpha-helix to coil transition of an alanine peptide in water". Proc. Natl. Acad. Sci. USA 88 2854-2858.
  25. Gavezotti A. and Filippini G. (1994). "Geometry of intermolecular X-H...Y (X,Y=N,O) hydrogen bond and the calibration of empirical hydrogen-bond potentials". J. Phys. Chem. 98 4831-4837.
  26. 26.0 26.1 26.2 Gao, J. (1997). "Toward a Molecular Orbital Derived Empirical Potential for Liquid Simulations". J. Phys. Chem. B 101 657-663.
  27. 27.0 27.1 27.2 Xie, W., and Gao, J. (2007). "Design of a Next Generation Force Field: The X-POL Potential". J. Chem. Theory Comput. 3 1890-1900.
  28. Allinger, Norman L. (1977). "Conformational analysis. 130. MM2. A hydrocarbon force field utilizing V1 and V2 torsional terms". J. Am. Chem. Soc. 99 8127-8134.
  29. MM2 and MM3 home page
  30. Allinger, N.L., Yuh, Y.H., & Lii, J-H. (1989). "Molecular Mechanics. The MM3 Force Field for Hydrocarbons. 1". J. Am. Chem. Soc. 111 8551-8565.
  31. Lii, J-H., & Allinger, N.L. (1989). "Molecular Mechanics. The MM3 Force Field for Hydrocarbons. 2. Vibrational Frequencies and Thermodynamics". J. Am. Chem. Soc. 111 8566-8575.
  32. Lii, J-H., & Allinger, N.L. (1989). "Molecular Mechanics. The MM3 Force Field for Hydrocarbons. 3. The van der Waals Potentials and Crystal data for Aliphatic and Aromatic Hydrocarbons". J. Am. Chem. Soc. 111 8576-8582.
  33. Schaumann, T., Braun, W. and Wutrich, K. (1990). "The program FANTOM for energy refinement of polypeptides and proteins using a Newton-Raphson minimizer in torsion angle space". Biopolymers 29 679-694.
  34. Möllhoff, Margit and Sternberg, Ulrich (2001). "Molecular Mechanics with fluctuating atomic charges - a new force field with semi-empirical charge calculation". J. Mol. Model. 7 90-102.
  35. Warshel, A. (1973). "Quantum Mechanical Consistent Force Field (QCFF/PI) Method: Calculations of Energies, Conformations and Vibronic Interactions of Ground and Excited States of Conjugated Molecules" Israel J. Chem. 11 709.
  36. Warshel, A. and Levitt, M. (1974). "QCFF/PI: A Program for the Consistent Force Field Evaluation of Equilibrium Geometries and Vibrational Frequencies of Molecules". QCPE 247, Quantum Chemistry Program Exchange, Indiana University.
  37. Template:Cite doi
  38. http://www.gromacs.org/Documentation/Terminology/Force_Fields/GROMOS
  39. http://lists.gromacs.org/pipermail/gmx-users/2008-March/032661.html
  40. Template:Cite doi
  41. Xie, W., Orozco, M., Truhlar, D. G., and Gao, J. (2009) "X-Pol Potential: An Electronic Structure-Based Force Field for Molecular Dynamics Simulation of a Solvated Protein in Water". J. Chem. Theory Comput. 5 459-467.
  42. McCammon, J. A., Gelin, B. R., and Karplus, M. (1977) "Dynamics of folded proteins". Nature 267 585-590.
  43. Warshel A. and Levitt M. (1976) "Theoretical Studies of Enzymatic Reactions: Dielectric Electrostatic and Steric Stabilization of the Carbonium Ion in the Reaction of Lysozyme". J. Mol. Biol. 103 227-249.
  44. Gao, J., Habibollahzadeh, D., and Shao, L. (1995) "A Polarizable Intermolecular Potential Functions for Simulations of Liquid Alcohols". J. Phys. Chem. 99 16460-16467.
  45. Xie, W., Pu, J., MacKerell, A. D., Jr., and Gao, J. (2007) "Development of a Polarizable Intermolecular Potential Function (PIPF) for Liquid Amides and Alkanes". J. Chem. Theory Comput. 3 1878-1889.
  46. Patel, S., Brooks, C.L., III (2004) "CHARMM fluctuating charge force field for proteins: I parameterization and application to bulk organic liquid simulations." J. Comput. Chem. 25 1-16.
  47. Patel, S., MacKerell, A.D., Jr., Brooks, C.L., III (2004) "CHARMM fluctuating charge force field for proteins: II Protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model." J. Comput. Chem. 25 1504-1514.
  48. Anisimov, V.M., Lamoureux, G., Vorobyov, I.V., Huang, N., Roux, B. and MacKerell, A.D., Jr. (2005) "Determination of Electrostatic Parameters for a Polarizable Force Field Based on the Classical Drude Oscillator". J. Chem. Theory Comput. 1 153-168.
  49. Yu, H., Whitfield, T.W., Harder, E., Lamoureux, G., Vorobyov, I., Anisimov, V. M., MacKerell, A.D., Jr., and Roux, B. (2010) "Simulating Monovalent and Divalent Ions in Aqueous Solution Using a Drude Polarizable Force Field". J. Chem. Theory Comput. 6 774–786.
  50. 50.0 50.1 Gresh, N., Cisneros, G. A., Darden, T. A. and Piquemal, J-P. (2007) "Anisotropic, polarizable molecular mechanics studies of inter-, intra-molecular interactions, and ligand-macromolecule complexes. A bottom-up strategy." J. Chem. Theory. Comput. 3 1960.
  51. O. Engkvist, P.-O. Åstrand, and G. Karlström. "Accurate intermolecular potentials obtained from molecular wave functions: Bridging the gap between quantum chemistry and molecular simulations". Chem. Rev., 100:4087-4108, 2000.
  52. Piquemal, J-P, Cisneros, G.A., Reinhardt, P., Gresh, N. and Darden, T.A. (2006). "Towards a Force Field based on Density Fitting". J. Chem. Phys. 124 104101.
  53. Cisneros, G.A., Piquemal, J-P and Darden, T.A. (2006). "Generalization of the Gaussian Electrostatic Model: extension to arbitrary angular momentum, distributed multipoles and speedup with reciprocal space methods". J. Chem. Phys. 125 184101.
  54. Sternberg, Ulrich, Koch, Frank-Thomas and Möllhoff, Margit (1994). "New Approach to the Semiempirical Calculation of Atomic Charges for Polypeptides and Large Molecular Systems" J. Comp. Chem. 15 524-531.
  55. Template:Cite doi
  56. Template:Cite doi
  57. A. Korkut, W.A. Hendrickson (2009). "A force field for virtual atom molecular mechanics of proteins." Proc. Natl. Acad. Sci. USA 106 15667–15672.
  58. S.J. Marrink, H.J. Risselada, S. Yefimov,D.P. Tieleman,A.H. de Vries (2007). "The MARTINI force field: coarse grained model for biomolecular simulations." J. Phys. Chem. B. 111 7812-7824.
  59. Template:Cite doi
  60. Template:Cite doi