Editing AMBER

{{Short description|Molecular dynamics software package}}
{{other uses|Amber (disambiguation)}}
{{Infobox software
| name = Assisted Model Building with Energy Refinement (AMBER)
| logo = 
| screenshot = 
| caption = 
| collapsible = 
| author = [[Peter Kollman]], Version1: +Paul Weiner. Version 2: + U Chandra Singh; V3. + David Pearlman, James Caldwell, William Ross, Thomas Cheatham, Stephen Debolt, David Ferguson, George Seibel; Later versions: David Case, Tom Cheatham, Ken Merz, Adrian Roitberg, Carlos Simmerling, Ray Luo, Junmei Wang, Ross Walker
| developer = [[University of California, San Francisco]]
| released = {{Start date and age|1981}}
| latest release version = Amber23, AmberTools23<ref name="Amber23Manual">[https://ambermd.org/doc12/Amber23.pdf Amber 2023 Reference Manual]</ref>
| latest release date = {{Start date and age|2023|04|21}}
| latest preview version = 
| latest preview date = 
| programming language = [[C (programming language)|C]], [[C++]], [[Fortran]]
| operating system = [[Microsoft Windows|Windows]], [[OS X]], [[Linux]], [[Unix]], [[CNK operating system|CNK]]
| platform = [[x86]], [[List of Nvidia graphics processing units|Nvidia GPUs]], [[Blue Gene]]
| size = Varies
| language = English
| genre = [[Molecular dynamics]]
| license = Amber: [[Proprietary software|Proprietary]]<br />AmberTools: [[GNU General Public License|GPL]], [[public domain]], other [[open-source software|open-source]]
| website = {{URL|https://ambermd.org}}
}}

[[File:Bond stretching energy.png|thumb|right|AMBER is used to minimize the bond stretching energy of this [[ethane]] molecule.]]

'''Assisted Model Building with Energy Refinement''' ('''AMBER''') is the name of a widely-used molecular dynamics software package originally developed by [[Peter Kollman]]'s group at the [[University of California, San Francisco]]. It has also, subsequently, come to designate a family of [[Force field (chemistry)|force fields]] for [[molecular dynamics]] of [[biomolecule]]s that can be used both within the AMBER software suite and with many modern computational platforms. 

The original version of the '''AMBER''' [[Software package (installation)|software package]] was written by Paul Weiner as a post-doc in Peter Kollman's laboratory, and was released in 1981.<ref>{{Cite journal |last=Weiner |first=Paul K. |last2=Kollman |first2=Peter A. |date=1981 |title=AMBER : Assisted model building with energy refinement. A general program for modeling molecules and their interactions |url=https://onlinelibrary.wiley.com/doi/10.1002/jcc.540020311 |journal=Journal of Computational Chemistry |language=en |volume=2 |issue=3 |pages=287–303 |doi=10.1002/jcc.540020311 |issn=0192-8651|url-access=subscription }}</ref> 

Subsequently, U Chandra Singh expanded AMBER as a post-doc in Kollman's laboratory, adding molecular dynamics and free energy capabilities. 

The next iteration of AMBER was started around 1987 by a group of developers in (and associated with) the Kollman lab, including David Pearlman, David Case, James Caldwell, William Ross, Thomas Cheatham, Stephen DeBolt, David Ferguson, and George Seibel.<ref>{{Cite journal |last=Pearlman |first=David A. |last2=Case |first2=David A. |last3=Caldwell |first3=James W. |last4=Ross |first4=Wilson S. |last5=Cheatham |first5=Thomas E. |last6=DeBolt |first6=Steve |last7=Ferguson |first7=David |last8=Seibel |first8=George |last9=Kollman |first9=Peter |date=1995 |title=AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules |url=https://doi.org/10.1016/0010-4655(95)00041-D |journal=Computer Physics Communications |volume=91 |issue=1-3 |pages=1–41 |doi=10.1016/0010-4655(95)00041-d |issn=0010-4655|url-access=subscription }}</ref> This team headed development for more than a decade and introduced a variety of improvements, including significant expansion of the free energy capabilities, accommodation for modern parallel and array processing hardware platforms (Cray, Star, etc.), restructuring of the code and revision control for greater maintainability, PME Ewald summations, tools for NMR refinement, and many others. 

Currently, AMBER is maintained by an active collaboration between David Case at [[Rutgers University]], Tom Cheatham at the [[University of Utah]], Adrian Roitberg at [[University of Florida]], Ken Merz at Michigan State University, [[Carlos Simmerling]] at [[Stony Brook University]], Ray Luo at [[UC Irvine]], and Junmei Wang at [[University of Pittsburgh]].

== Force field ==
The term ''AMBER [[Force field (chemistry)|force field]]'' generally refers to the functional form used by the family of AMBER force fields. This form includes several parameters; each member of the family of AMBER force fields provides values for these parameters and has its own name.

=== Functional form ===
The functional form of the AMBER force field is<ref name="Cornell1995">{{cite journal |vauthors=Cornell WD, Cieplak P, Bayly CI, Gould IR, ((Merz KM Jr)), Ferguson DM, Spellmeyer DC, Fox T, Caldwell JW, Kollman PA |title=A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules |journal=J. Am. Chem. Soc. |volume=117 |issue=19 |pages=5179–5197 |year=1995 |doi=10.1021/ja00124a002|citeseerx=10.1.1.323.4450 }}</ref>
:<math>
V(r^N)=\sum_{i \in \text{bonds}} {k_b}_i (l_i-l_i^0)^2 + \sum_{i \in \text{angles}} {k_a}_i (\theta_i - \theta_i^0)^2</math>
<blockquote>
<math>+ \sum_{i \in \text{torsions}} \sum_n \frac{1}{2} V_i^n [1+\cos(n \omega_i - \gamma_i)]</math> 
<math>+\sum_{j=1} ^{N-1} \sum_{i=j+1} ^N f_{ij}\biggl\{\epsilon_{ij}\biggl[\left(\frac{r^{0}_{ij}}{r_{ij}} \right)^{12} - 2\left(\frac{r^{0}_{ij}}{r_{ij}} \right)^{6} \biggr]+ \frac{q_iq_j}{4\pi \epsilon_0 r_{ij}}\biggr\}
</math>
</blockquote>
Despite the term ''force field'', this equation defines the potential energy of the system; the force is the derivative of this potential relative to position.

The meanings of right hand side [[term (logic)|term]]s are:
* First term ([[summation|summing]] over bonds): represents the energy between covalently bonded atoms. This harmonic (ideal spring) force is a good approximation near the equilibrium bond length, but becomes increasingly poor as atoms separate.
* Second term (summing over angles): represents the energy due to the geometry of electron orbitals involved in covalent bonding.
* Third term (summing over torsions): represents the energy for twisting a bond due to bond order (e.g., double bonds) and neighboring bonds or lone pairs of electrons. One bond may have more than one of these terms, such that the total torsional energy is expressed as a [[Fourier series]].
* Fourth term (double summation over <math>i</math> and <math>j</math>): represents the non-bonded energy between all atom pairs, which can be decomposed into [[van der Waals force|van der Waals]] (first term of summation) and [[electrostatics|electrostatic]] (second term of summation) energies.

The form of the van der Waals energy is calculated using the equilibrium distance (<math> r^{0}_{ij} </math>) and well depth (<math> \epsilon </math>). The factor of <math>2</math> ensures that the equilibrium distance is <math> r^{0}_{ij} </math>. The energy is sometimes reformulated in terms of <math>\sigma</math>, where <math> r^{0}_{ij} = 2^{1/6}(\sigma)</math>, as used e.g. in the implementation of the softcore potentials.

The form of the electrostatic energy used here assumes that the charges due to the protons and electrons in an atom can be represented by a single point charge (or in the case of parameter sets that employ lone pairs, a small number of point charges.)

=== Parameter sets ===
To use the AMBER force field, it is necessary to have values for the parameters of the force field (e.g. force constants, equilibrium bond lengths and angles, charges). A fairly large number of these parameter sets exist, and are described in detail in the AMBER software user manual. Each parameter set has a name, and provides parameters for certain types of molecules.
*[[Peptide]], [[protein]], and [[nucleic acid]] parameters are provided by parameter sets with names starting with "ff" and containing a two digit year number, for instance "ff99". As of 2018 the primary protein model used by the AMBER suit is the ff14SB<ref>{{Cite journal |doi = 10.1021/acs.jctc.5b00255|pmid = 26574453|pmc = 4821407|title = Ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB|journal = Journal of Chemical Theory and Computation|volume = 11|issue = 8|pages = 3696–3713|year = 2015|last1 = Maier|first1 = James A|last2 = Martinez|first2 = Carmenza|last3 = Kasavajhala|first3 = Koushik|last4 = Wickstrom|first4 = Lauren|last5 = Hauser|first5 = Kevin E|last6 = Simmerling|first6 = Carlos}}</ref><ref>{{Cite web|url=http://ambermd.org/AmberModels.php|title = The Amber Force Fields}}</ref> force field.
* ''General AMBER force field'' (GAFF) provides parameters for small organic molecules to facilitate simulations of drugs and small molecule ligands in conjunction with biomolecules.
*The GLYCAM force fields have been developed by Rob Woods for simulating carbohydrates.
*The primary force field used in the AMBER suit for lipids is lipid14.<ref>{{Cite journal |doi = 10.1021/ct4010307|pmid = 24803855|pmc = 3985482|title = Lipid14: The Amber Lipid Force Field|journal = Journal of Chemical Theory and Computation|volume = 10|issue = 2|pages = 865–879|year = 2014|last1 = Dickson|first1 = Callum J|last2 = Madej|first2 = Benjamin D|last3 = Skjevik|first3 = Åge A|last4 = Betz|first4 = Robin M|last5 = Teigen|first5 = Knut|last6 = Gould|first6 = Ian R|last7 = Walker|first7 = Ross C}}</ref>

== Software ==
The AMBER software suite provides a set of programs to apply the AMBER forcefields to simulations of biomolecules. It is written in the programming languages [[Fortran 90]] and [[C (programming language)|C]], with support for most major [[Unix-like]] operating systems and [[compiler]]s. Development is conducted by a loose association of mostly academic labs. New versions are released usually in the spring of even numbered years; AMBER 10 was released in April 2008. The software is available under a [[site license]] agreement, which includes full source, currently priced at US$500 for non-commercial and US$20,000 for commercial organizations.

=== Programs ===
* ''LEaP'' prepares input files for the simulation programs.
* ''Antechamber'' automates the process of parameterizing small organic molecules using GAFF.
* ''Simulated Annealing with NMR-Derived Energy Restraints'' (SANDER) is the central simulation program and provides facilities for energy minimizing and molecular dynamics with a wide variety of options.
* ''pmemd'' is a somewhat more feature-limited reimplementation of SANDER by Bob Duke. It was designed for [[parallel computing]], and performs significantly better than SANDER when running on more than 8–16 processors.
** ''pmemd.cuda'' runs simulations on machines with [[graphics processing unit]]s (GPUs).
** ''pmemd.amoeba'' handles the extra parameters in the polarizable AMOEBA force field.
* ''nmode'' calculates normal modes.
* ''ptraj'' numerically analyzes simulation results. AMBER includes no visualizing abilities, which is commonly performed with [[Visual Molecular Dynamics]] (VMD). Ptraj is now unsupported as of AmberTools 13. 
* ''cpptraj'' is a rewritten version of ptraj made in [[C++]] to give faster analysis of simulation results. Several actions have been made parallelizable with OpenMP and MPI.
* ''MM-PBSA'' allows implicit solvent calculations on snap shots from molecular dynamics simulations.
* ''NAB'' is a built-in nucleic acid building environment made to aid in the process of manipulating proteins and nucleic acids where an atomic level of description will aid computing.

== See also ==
{{columns-list|colwidth=30em|
* [[Comparison of software for molecular mechanics modeling]]
* [[Comparison of force field implementations]]
* [[Molecular dynamics]]
* [[Molecular geometry]]
* [[Molecular design software]]
* [[Molecular mechanics]]
* [[MDynaMix]]
* [[Ascalaph Designer]]
* [[BOSS (molecular mechanics)]]
* [[CHARMM]]
* [[GROMACS]]
*[[Reverse Pharmacology]]
* [[OPLS]]
* [[Yasara]]
* [[Folding@home]]
}}

== References ==
{{reflist}}
===Related reading===<!-- These should be moved inline. -->
1. {{cite journal |last1=Duan |first1=Yong |last2=Wu |first2=Chun |last3=Chowdhury |first3=Shibasish |last4=Lee |first4=Mathew C. |last5=Xiong |first5=Guoming |last6=Zhang |first6=Wei |last7=Yang |first7=Rong |last8=Cieplak |first8=Piotr |last9=Luo |first9=Ray |display-authors=8|title=A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations |journal=Journal of Computational Chemistry |volume=24 |issue=16 |pages=1999–2012 |year=2003 |doi=10.1002/jcc.10349|pmid=14531054 |s2cid=283317 }}

== External links ==
* {{Official website|https://ambermd.org}}
*[http://archive.ambermd.org/ AMBER mailing list archive]
*[https://wiki.bwhpc.de/e/Amber Amber on the German HPC-C5 Cluster-Systems]

{{Chemistry software}}

[[Category:Fortran software]]
[[Category:Molecular dynamics software]]
[[Category:Force fields (chemistry)]]