From tinkergpu
Jump to navigation Jump to search

AMOEBA library

Official AMOEBA parameters are distributed with Tinker: Tinkertools Github

  • Small molecules, ions (2009): amoeba09.prm;
  • Proteins (2013): amoebapro13.prm;
  • Nucleic acids DNA/RNA (2017): amoebanuc17.prm
  • Most current, combined biomolecule force field (2018): amoebabio18.prm
  • Water: please use amoeba03 parameters within the prm files above.

Improvements made after the Mg/Ca studies

(Jing et al PNAS 2018, JCP 2017)

These changes are incoporated into the latest amoebabio18.prm and amoebapro13.prm as of 2020

  • Refit the Glu gamma C and H multipoles. The quadrupoles on H (Qxz=0.97) were too large and caused instability if two Glu are close to each other. The refit to EPS was done with Qxz forced to 0.
  • Refit Ca and Mg parameters based on Mg/Ca-water MP2/CBS interaction and HFE. The original parameter development used MP2/6-311++G(3df,3pd) or aTZ.
  • Change COO- O vdW and polarizability based on interactions between COO- and water/Ca/Mg and acetate HFE. Polarizability 1.200 is between 0.837 (amoeba09) and 1.7 (amoeba13) and better reproduces the many-body interactions between multiple COO- and Mg/Ca.
  • Set the quadupoles on COO- (C and O) to 0 to reproduce the angular dependence of COO- Ca/Mg interactions. Before the changes, the bidentate binding was too strong and monodentate binding too weak.
  • Use vdwpr for O=C and Mg/Ca to better reproduce the interaction energy, but the equilibrium distance is still slightly shorter compared to QM
  • Use vdwpr for COO- O and amide H. The interactions were too weak (probably due to lack of charge transfer)
  • Thr O and H quardrupoles were scaled by 60%. But we found it is better to reverse this, based on Thr interaction with COO-
  • Refit Asp chi1 and chi2 and Glu chi1 and chi3 torsion parameters based on wB97xD/PCM. The old torsion parameters could not produce the correct pocket geometries for EF hands.

For new small molecules

Recent updates & next-gen development

  • Halegen compounds (2014): Mu et al, JPCB 2014, 118 (24), pp 6456. Not yet in official release of small molecule params.
  • Mg, Ca ions (2017): minor update included in ameobapro13.prm in summer 2017
  • SAPT database between organic molecules, protein side chain analogs
      http://biomol.bme.utexas.edu/~ch38988/s101x7 [1]
  • General Model for Treating Short-Range Electrostatic Penetration in a Molecular Mechanics Force Field, Qiantao Wang, Joshua A. Rackers, Chenfeng He, Rui Qi, Christophe Narth, Louis Lagardere, Nohad Gresh, Jay W. Ponder, Jean-Philip Piquemal , and Pengyu Ren J. Chem. Theory Comput., 2015, 11 (6), pp 2609–2618
  • General van der Waals potential for common organic molecules, Rui Qi, Qiantao Wang, Pengyu Ren (2016) Bioorganic & Medicinal Chemistry, 24, 4911–4919
  • Capturing Many-body Interactions with Classical Dipole Induction Models, C Liu, R Qi, Q Wang, JP Piquemal, P Ren, Journal of Chemical Theory and Computation 2017, 13 (6), 2751–2761





  • ForceBalance can be installed easily, only if your local python settings are proper. Otherwise, I recommend you to install a fresh python 2.7  
  • Install Python 2.7.14
    • Download the 2.7.14 compressed file from python website; release the file then do the following 
      • cd ...../Python-2.7.14/
      • ./configure --prefix=/your/local/directory
      • make
      • make install
      • After installation, add the prefix to your $PATH environmental variable
    • Install the required python modules: numpy v1.5, scipy v0.9 (scipy depends on numpy, so install numpy first), lxml (if you want to use OpenMM interface, otherwise you can ignore this)
      • Exactly the same commands are used:  {python setup.py build} then {python setup.py install}
      • If other modules are also needed, I recommend git clone the modules then use the above two commands to install
  • Install ForceBalance
    • git clone https://github.com/leeping/forcebalance (this will clone the whole FB files to your $currentworkingdirectory)
    • cd ...../forcebalance
    • python setup.py build (note that if --prefix is not specified, forcebalance module will appear in your python site-packages) 
    • python setup.py install
    • Note that python does not always recognize the modules in the custom location. You may need to explicitly set $PYTHONPATH to the location where the modules are.
  • Run a case to test your installation (in forcebalance/studies/ choose one to run.  Note that if you want to run with TINKER, you need to set the tinker_path)

References and citations

Ren, P. and J.W. Ponder, Polarizable atomic multipole water model for molecular mechanics simulation. The Journal of Physical Chemistry B, 2003. 107(24): p. 5933-5947.

Ponder, J.W., C. Wu, P. Ren, V.S. Pande, J.D. Chodera, M.J. Schnieders, I. Haque, D.L. Mobley, D.S. Lambrecht, and R.A. DiStasio Jr, Current status of the AMOEBA polarizable force field. The journal of physical chemistry B, 2010. 114(8): p. 2549-2564.

Ren, P., C. Wu, and J.W. Ponder, Polarizable atomic multipole-based molecular mechanics for organic molecules. Journal of chemical theory and computation, 2011. 7(10): p. 3143-3161.

Wu, J.C., G. Chattree, and P. Ren, Automation of AMOEBA polarizable force field parameterization for small molecules. Theoretical chemistry accounts, 2012. 131(3): p. 1138.

Shi, Y., Z. Xia, J. Zhang, R. Best, C. Wu, J.W. Ponder, and P. Ren, Polarizable atomic multipole-based AMOEBA force field for proteins. Journal of chemical theory and computation, 2013. 9(9): p. 4046-4063.