CHARMM36m: an improved force field for folded and intrinsically disordered proteins

Jing Huang,Sarah Rauscher,Grzegorz Nawrocki,Ting Ran,Michael Feig,Bert L de Groot,Helmut Grubmüller,Alexander D MacKerell

DOI: https://doi.org/10.1038/nmeth.4067

IF: 48

2016-11-07

Nature Methods

Abstract:An all-atom protein force field, CHARMM36m, offers improved accuracy for simulating intrinsically disordered peptides and proteins.

biochemical research methods

What problem does this paper attempt to address?

Comparison and Evaluation of Force Fields for Intrinsically Disordered Proteins

Mueed Ur Rahman,Ashfaq Ur Rehman,Hao Liu,Hai-Feng Chen

DOI: https://doi.org/10.1021/acs.jcim.0c00762

IF: 6.162

2020-08-20

Journal of Chemical Information and Modeling

Abstract:Molecular dynamics (MD) simulations of six upgraded empirical force fields were compared and evaluated with short peptides, intrinsically disordered proteins, and folded proteins using trajectories of 1, 1.5, 5, or 10 μs (five replicates of 200 ns, 300 ns, 1 μs, or 2 μs) for each system. Previous studies have shown that different force fields, water models, simulation methods, and parameters can affect simulation outcomes. Here, the MD simulations were done in an explicit solvent with RS-peptide, HEWL19, HIV-rev, β amyloid (Aβ)-40, Aβ-42, phosphodiesterase-γ, CspTm, and ubiquitin using ff99IDPs, ff14IDPs, ff14IDPSFF, ff03w, CHARMM36m, and CHARMM22* force fields. The IDP ensembles generated by six all-atom empirical force fields were compared against NMR data. Despite using identical starting structures and simulation parameters, ensembles obtained with different force fields exhibit significant differences in NMR RMDs, secondary structure contents, and global properties such as the radius of gyration. The intrinsically disordered protein (IDP)-specific force fields could substantially reproduce the experimental observables in force field comparison: they have the lowest error in chemical shifts and J-couplings for short peptides/proteins, reasonably well for large IDPs and reasonably well with the radius of gyration. A high population of disorderness was observed in the IDP-specific force field for the IDP ensemble with a fraction of β sheets for β-amyloids. CHARMM22* performs better for many observables; however, it still has a preference toward the helicity for short peptides. The results of β-amyloid 42 starting from two different initial structures (Aβ421Z0Q and Aβ42model) were also compared with DSSP and NMR data. The results obtained with IDP-specific force fields within 2 μs simulation time are similar, even though starting from different structures. The current force fields perform equally well for folded proteins. The results of currently developed or modified force fields for IDPs are capable of enlightening the overall performance of the force field for disordered as well as folded proteins, thereby contributing to force field development.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jcim.0c00762?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jcim.0c00762</a>.Predicted measurement comparison (Table S1); secondary chemical shift conversion (Figure S1); time evolution of the cumulative Rg and the end-to-end distance (Figure S2); comparison of secondary chemical shifts between MD simulations and NMR experiments (Figure S3); cluster analyses (Figure S4); calculated NMR observables in simulations of Ab40 starting from the PDB structure and the extended structure (Figure S5); secondary structure propensity analysis of the trajectories (Figure S6); RMSD analyses of all of the trajectories under the tested force fields (Figure S7) (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.0c00762/suppl_file/ci0c00762_si_001.pdf">PDF</a>)This article has not yet been cited by other publications.

chemistry, multidisciplinary, medicinal,computer science, interdisciplinary applications, information systems
Significantly Improved Protein Folding Thermodynamics Using a Dispersion-Corrected Water Model and a New Residue-Specific Force Field

Hao-Nan Wu,Fan Jiang,Yun-Dong Wu

DOI: https://doi.org/10.1021/acs.jpclett.7b01213

2017-01-01

Abstract:An accurate potential energy model is crucial for biomolecular simulations. Despite many recent improvements of classical protein force fields, there are remaining key issues: much weaker temperature dependence of folding/unfolding equilibrium and overly collapsed unfolded or disordered states. For the latter problem, a new water model (TIP4P-D) has been proposed to correct the significantly underestimated water dispersion interactions. Here, using TIP4P-D, we reveal problems in current force fields through failures in folding model systems (a polyalanine peptide, Trp-cage, and the GB1 hairpin). By using residue-specific parameters to achieve better match between amino acid sequences and native structures and adding a small H-bond correction to partially compensate the missing many-body effects in a-helix formation, the new RSFF2+ force field with the TIP4P-D water model can excellently reproduce experimental melting curves of both alpha-helical and beta-hairpin systems. The RSFF2+/TIP4P-D method also gives less collapsed unfolded structures and describes well folded proteins simultaneously.
Polarizable Force Field of Intrinsically Disordered Proteins with CMAP and Reweighting Optimization

Xiaochen Cui,Hao Liu,Hai-Feng Chen

DOI: https://doi.org/10.1021/acs.jcim.2c00835

IF: 6.162

2022-10-01

Journal of Chemical Information and Modeling

Abstract:Intrinsically disordered proteins (IDPs) are highly structurally heterogeneous without a specific tertiary structure under physiology conditions and play key roles in the development of human diseases. Due to the characteristics of diverse conformations, as one of the important methods, molecular dynamics simulation can complement information for experimental methods. Because of the enrichment for charged amino acids for IDPs, polarizable force fields should be a good choice for the simulation...

chemistry, multidisciplinary, medicinal,computer science, interdisciplinary applications, information systems
A Test of AMBER Force Fields in Predicting the Secondary Structure of Α-Helical and Β-Hairpin Peptides

Ya Gao,Chaomin Zhang,Xianwei Wang,Tong Zhu

DOI: https://doi.org/10.1016/j.cplett.2017.04.074

IF: 2.719

2017-01-01

Chemical Physics Letters

Abstract:We tested the ability of some current AMBER force fields, namely, AMBER03, AMBER99SB, AMBER99SB-ildn, AMBER99SB-nmr, AMBER12SB, AMBER14SB, and AMBER14ipq, with implicit solvent model in reproducing the folding behavior of two peptides by REMD simulations. AMBER99SB-nmr force field provides the most reliable performance. After a novel polarized hydrogen bond charge model is considered, the α-helix successfully folded to its native state, while the further folding of the β-hairpin is not observed. This study strongly suggests that polarization effect and correct torsional term are important to investigate dynamic and conformational properties of peptides with different secondary structures.
Toward a general neural network force field for protein simulations: Refining the intramolecular interaction in protein

Pan Zhang,Weitao Yang

DOI: https://doi.org/10.1063/5.0142280

2023-07-14

Abstract:Molecular dynamics (MD) is an extremely powerful, highly effective, and widely used approach to understanding the nature of chemical processes in atomic details for proteins. The accuracy of results from MD simulations is highly dependent on force fields. Currently, molecular mechanical (MM) force fields are mainly utilized in MD simulations because of their low computational cost. Quantum mechanical (QM) calculation has high accuracy, but it is exceedingly time consuming for protein simulations. Machine learning (ML) provides the capability for generating accurate potential at the QM level without increasing much computational effort for specific systems that can be studied at the QM level. However, the construction of general machine learned force fields, needed for broad applications and large and complex systems, is still challenging. Here, general and transferable neural network (NN) force fields based on CHARMM force fields, named CHARMM-NN, are constructed for proteins by training NN models on 27 fragments partitioned from the residue-based systematic molecular fragmentation (rSMF) method. The NN for each fragment is based on atom types and uses new input features that are similar to MM inputs, including bonds, angles, dihedrals, and non-bonded terms, which enhance the compatibility of CHARMM-NN to MM MD and enable the implementation of CHARMM-NN force fields in different MD programs. While the main part of the energy of the protein is based on rSMF and NN, the nonbonded interactions between the fragments and with water are taken from the CHARMM force field through mechanical embedding. The validations of the method for dipeptides on geometric data, relative potential energies, and structural reorganization energies demonstrate that the CHARMM-NN local minima on the potential energy surface are very accurate approximations to QM, showing the success of CHARMM-NN for bonded interactions. However, the MD simulations on peptides and proteins indicate that more accurate methods to represent protein-water interactions in fragments and non-bonded interactions between fragments should be considered in the future improvement of CHARMM-NN, which can increase the accuracy of approximation beyond the current mechanical embedding QM/MM level.
High quality NMR structures: a new force field with implicit water and membrane solvation for Xplor-NIH

Ye Tian,Charles D. Schwieters,Stanley J. Opella,Francesca M. Marassi

DOI: https://doi.org/10.1007/s10858-016-0082-5

2016-12-29

Journal of Biomolecular NMR

Abstract:Structure determination of proteins by NMR is unique in its ability to measure restraints, very accurately, in environments and under conditions that closely mimic those encountered in vivo. For example, advances in solid-state NMR methods enable structure determination of membrane proteins in detergent-free lipid bilayers, and of large soluble proteins prepared by sedimentation, while parallel advances in solution NMR methods and optimization of detergent-free lipid nanodiscs are rapidly pushing the envelope of the size limit for both soluble and membrane proteins. These experimental advantages, however, are partially squandered during structure calculation, because the commonly used force fields are purely repulsive and neglect solvation, Van der Waals forces and electrostatic energy. Here we describe a new force field, and updated energy functions, for protein structure calculations with EEFx implicit solvation, electrostatics, and Van der Waals Lennard-Jones forces, in the widely used program Xplor-NIH. The new force field is based primarily on CHARMM22, facilitating calculations with a wider range of biomolecules. The new EEFx energy function has been rewritten to enable OpenMP parallelism, and optimized to enhance computation efficiency. It implements solvation, electrostatics, and Van der Waals energy terms together, thus ensuring more consistent and efficient computation of the complete nonbonded energy lists. Updates in the related python module allow detailed analysis of the interaction energies and associated parameters. The new force field and energy function work with both soluble proteins and membrane proteins, including those with cofactors or engineered tags, and are very effective in situations where there are sparse experimental restraints. Results obtained for NMR-restrained calculations with a set of five soluble proteins and five membrane proteins show that structures calculated with EEFx have significant improvements in accuracy, precision, and conformation, and that structure refinement can be obtained by short relaxation with EEFx to obtain improvements in these key metrics. These developments broaden the range of biomolecular structures that can be calculated with high fidelity from NMR restraints.

biochemistry & molecular biology,spectroscopy
Environment-Specific Force Field for Intrinsically Disordered and Ordered Proteins

Dong Song,Hao Liu,Ray Luo,Hai-Feng Chen

DOI: https://doi.org/10.1021/acs.jcim.0c00059

IF: 6.162

2020-03-31

Journal of Chemical Information and Modeling

Abstract:The need for accurate and efficient force fields for modeling 3D structures of macrobiomolecules and in particular intrinsically disordered proteins (IDPs) has increased with recent findings to associate IDPs and human diseases. However, most conventional protein force fields and recent IDP-specific force fields are limited in reproducing accurate structural features of IDPs. Here, we present an environmental specific precise force field (ESFF1) based on CMAP corrections of 71 different sequence environments to improve the accuracy and efficiency of MD simulation for both IDPs and folded proteins. MD simulations of 84 different short peptides, IDPs, and structured proteins show that ESFF1 can accurately reproduce spectroscopic properties for different peptides and proteins whether they are disordered or ordered. The successful ab initio folding of five fast-folding proteins further supports the reliability of ESFF1. The extensive analysis documented here shows that ESFF1 is able to achieve a reasonable balance between ordered and disordered states in protein simulations.

chemistry, multidisciplinary, medicinal,computer science, interdisciplinary applications, information systems
Likely Overstabilization of Charge-Charge Interactions in CHARMM36m(w): A Case for a99SB-disp Water

Xiping Gong,Yumeng Zhang,Jianhan Chen

DOI: https://doi.org/10.1021/acs.jpcb.4c04777

IF: 3.466

2024-11-18

The Journal of Physical Chemistry B

Abstract:Recent years have witnessed drastic improvements in general-purpose explicit solvent protein force fields, partially driven by the need to study intrinsically disordered proteins (IDPs), and yet the state-of-the-art force fields such as CHARMM36m (c36m) and a99SB-disp still provide different performances in simulating disordered protein states, where c36m has a bias toward overcompaction for large IDPs. Here, we examine the performance of c36m and a99SB-disp in describing the stabilities of a...

chemistry, physical
Improved side-chain torsion potentials for the Amber ff99SB protein force field

Kresten Lindorff‐Larsen,Stefano Piana,Kim Palmo,Paul Maragakis,John L. Klepeis,Ron O. Dror,David E. Shaw,Kresten Lindorff-Larsen

DOI: https://doi.org/10.1002/prot.22711

2010-03-09

Abstract:Recent advances in hardware and software have enabled increasingly long molecular dynamics (MD) simulations of biomolecules, exposing certain limitations in the accuracy of the force fields used for such simulations and spurring efforts to refine these force fields. Recent modifications to the Amber and CHARMM protein force fields, for example, have improved the backbone torsion potentials, remedying deficiencies in earlier versions. Here, we further advance simulation accuracy by improving the amino acid side-chain torsion potentials of the Amber ff99SB force field. First, we used simulations of model alpha-helical systems to identify the four residue types whose rotamer distribution differed the most from expectations based on Protein Data Bank statistics. Second, we optimized the side-chain torsion potentials of these residues to match new, high-level quantum-mechanical calculations. Finally, we used microsecond-timescale MD simulations in explicit solvent to validate the resulting force field against a large set of experimental NMR measurements that directly probe side-chain conformations. The new force field, which we have termed Amber ff99SB-ILDN, exhibits considerably better agreement with the NMR data.

biochemistry & molecular biology,biophysics
Further optimization of a hybrid united-atom and coarse-grained force field for folding simulations: Improved backbone hydration and interactions between charged side chains.

Wei Han,Klaus Schulten

DOI: https://doi.org/10.1021/ct300696c

2012-01-01

Journal of Chemical Theory and Computation

Abstract:PACE, a hybrid force field which couples united-atom protein models with coarse-grained (CG) solvent, has been further optimized, aiming to improve itse ciency for folding simulations. Backbone hydration parameters have been re-optimized based on hydration free energies of polyalanyl peptides through atomistic simulations. Also, atomistic partial charges from all-atom force fields were combined with PACE in order to provide a more realistic description of interactions between charged groups. Using replica exchange molecular dynamics (REMD), ab initio folding using the new PACE has been achieved for seven small proteins (16 - 23 residues) with different structural motifs. Experimental data about folded states, such as their stability at room temperature, melting point and NMR NOE constraints, were also well reproduced. Moreover, a systematic comparison of folding kinetics at room temperature has been made with experiments, through standard MD simulations, showing that the new PACE may speed up the actual folding kinetics 5-10 times. Together with the computational speedup benefited from coarse-graining, the force field provides opportunities to study folding mechanisms. In particular, we used the new PACE to fold a 73-residue protein, 3D, in multiple 10 - 30 μs simulations, to its native states (C(α) RMSD ~ 0.34 nm). Our results suggest the potential applicability of the new PACE for the study of folding and dynamics of proteins.
Developing a molecular dynamics force field for both folded and disordered protein states

Paul Robustelli,Stefano Piana,David E. Shaw

DOI: https://doi.org/10.1073/pnas.1800690115

IF: 11.1

2018-05-07

Proceedings of the National Academy of Sciences

Abstract:Significance Many proteins that perform important biological functions are completely or partially disordered under physiological conditions. Molecular dynamics simulations could be a powerful tool for the structural characterization of such proteins, but it has been unclear whether the physical models (force fields) used in simulations are sufficiently accurate. Here, we systematically compare the accuracy of a number of different force fields in simulations of both ordered and disordered proteins, finding that each force field has strengths and limitations. We then describe a force field that substantially improves on the state-of-the-art accuracy for simulations of disordered proteins without sacrificing accuracy for folded proteins, thus broadening the range of biological systems amenable to molecular dynamics simulations.
Residue-Speci fi c Force Field Based on the Protein Coil Library

Fan Jiang,Chen-Yang Zhou,Yun-Dong Wu

2014-01-01

Abstract:Traditional protein force fields use one set of parameters for most of the 20 amino acids (AAs), allowing transferability of the parameters. However, a significant shortcoming is the difficulty to fit the Ramachandran plots of all AA residues simultaneously, affecting the accuracy of the force field. In this Feature Article, we report a new strategy for protein force field parametrization. Backbone and side-chain conformational distributions of all 20 AA residues obtained from protein coil library were used as the target data. The dihedral angle (torsion) potentials and some local nonbonded (1-4/1-5/1-6) interactions in OPLS-AA/L force field were modified such that the target data can be excellently reproduced by molecular dynamics simulations of dipeptides (blocked AAs) in explicit water, resulting in a new force field with AA-specific parameters, RSFF1. An efficient free energy decomposition approach was developed to separate the corrections on φ and ψ from the two-dimensional Ramachandran plots. RSFF1 is shown to reproduce the experimental NMR J-coupling constants of AA dipeptides better than other force fields. It has a good balance between α-helical and β-sheet secondary structures. It can successfully fold a set of α-helix proteins (Trp-cage and Homeodomain) and β-hairpins (Trpzip-2, GB1 hairpin), which cannot be consistently stabilized by other state-of-the-art force fields. Interestingly, the RSFF1 force field systematically overestimates the melting temperature (and the stability of native state) of these peptides/proteins. It has a potential application in the simulation of protein folding and protein structure refinement.
A Force Balanced Fragmentation Method for Ab Initio Molecular Dynamic Simulation of Protein

Mingyuan Xu,Tong Zhu,John Z. H. Zhang

DOI: https://doi.org/10.3389/fchem.2018.00189

IF: 5.5

2018-01-01

Frontiers in Chemistry

Abstract:A force balanced generalized molecular fractionation with conjugate caps (FB-GMFCC) method is proposed for ab initio molecular dynamic simulation of proteins. In this approach, the energy of the protein is computed by a linear combination of the QM energies of individual residues and molecular fragments that account for the two-body interaction of hydrogen bond between backbone peptides. The atomic forces on the caped H atoms were corrected to conserve the total force of the protein. Using this approach, ab initio molecular dynamic simulation of an Ace-(ALA)9-NME linear peptide showed the conservation of the total energy of the system throughout the simulation. Further a more robust 110 ps ab initio molecular dynamic simulation was performed for a protein with 56 residues and 862 atoms in explicit water. Compared with the classical force field, the ab initio molecular dynamic simulations gave better description of the geometry of peptide bonds. Although further development is still needed, the current approach is highly efficient, trivially parallel, and can be applied to ab initio molecular dynamic simulation study of large proteins.
New-Generation Amber United-Atom Force Field

Lijiang Yang,Chun-hu Tan,Meng-Juei Hsieh,Junmei Wang,Yong Duan,Piotr Cieplak,James Caldwell,Peter A. Kollman,Ray Luo

DOI: https://doi.org/10.1021/jp060163v

2006-01-01

Abstract:We have developed a new-generation Amber united-atom force field for simulations involving highly demanding conformational sampling such as protein folding and protein-protein binding. In the new united-atom force field, all hydrogens on aliphatic carbons in all amino acids are united with carbons except those on C alpha. Our choice of explicit representation of all protein backbone atoms aims at minimizing perturbation to protein backbone conformational distributions and to simplify development of backbone torsion terms. Tests with dipeptides and solvated proteins show that our goal is achieved quite successfully. The new united-atom force field uses the same new RESP charging scheme based on B3LYP/cc-pVTZ//HF/6-31g** quantum mechanical calculations in the PCM continuum solvent as that in the Duan et al. force field. van der Waals parameters are empirically refitted starting from published values with respect to experimental solvation free energies of amino acid side-chain analogues. The suitability of mixing new point charges and van der Waals parameters with existing Amber covalent terms is tested on alanine dipeptide and is found to be reasonable. Parameters for all new torsion terms are refitted based on the new point charges and the van der Waals parameters. Molecular dynamics simulations of three small globular proteins in the explicit TIP3P solvent are performed to test the overall stability and accuracy of the new united-atom force field. Good agreements between the united-atom force field and the Duan et al. all-atom force field for both backbone and side-chain conformations are observed. In addition, the per-step efficiency of the new united-atom force field is demonstrated for simulations in the implicit generalized Born solvent. A speedup around two is observed over the Duan et al. all-atom force field for the three tested small proteins. Finally, the efficiency gain of the new united-atom force field in conformational sampling is further demonstrated with a well-known toy protein folding system, an 18 residue polyalanine in distance-dependent dielectric. The new united-atom force field is at least a factor of 200 more efficient than the Duan et al. all-atom force field for ab initio folding of the tested peptide.
CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

Jumin Lee,Xi Cheng,Jason M. Swails,Min Sun Yeom,Peter K. Eastman,Justin A. Lemkul,Shuai Wei,Joshua Buckner,Jong Cheol Jeong,Yifei Qi,Sunhwan Jo,Vijay S. Pande,David A. Case,Charles L. Brooks,Alexander D. MacKerell,Jeffery B. Klauda,Wonpil Im

DOI: https://doi.org/10.1021/acs.jctc.5b00935

2015-12-03

Journal of Chemical Theory and Computation

Abstract:Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.

chemistry, physical,physics, atomic, molecular & chemical
Universal Implementation of a Residue-Specific Force Field Based on CMAP Potentials and Free Energy Decomposition

Wei Kang,Fan Jiang,Yun-Dong Wu

DOI: https://doi.org/10.1021/acs.jctc.8b00285

2018-01-01

Journal of Chemical Theory and Computation

Abstract:The coupling between neighboring backbone ϕ and ψ dihedral angles (torsions) has been well appreciated in protein force field development, as in correction map (CMAP) potentials. However, although preferences of backbone torsions are significantly affected by side-chain conformation, there has been no easy way to optimize this coupling. Herein, we prove that the three-dimensional (3D) free energy hypersurface of joint (ϕ, ψ, χ1) torsions can be decomposed into three separated 2D surfaces. Thus, each of the 2D torsional surfaces can be efficiently and automatically optimized using a CMAP potential. This strategy is then used to reparameterize an AMBER force field such that the resulting χ1-dependent backbone conformational preference can agree excellently with the reference protein coil library statistics. In various validation simulations (including the folding of seven peptides/proteins, backbone dynamics of three folded proteins, and two intrinsically disordered peptides), the new RSFF2C (residue-specific force field with CMAP potentials) force field gives similar or better performance compared with RSFF2. This strategy can be used to implement our RSFF force fields into a variety of molecular dynamics packages easily.
Direct Folding Simulation of Helical Proteins Using an Effective Polarizable Bond Force Field

Lili Duan,Tong Zhu,Changge Ji,Qinggang Zhang,John Z. H. Zhang

DOI: https://doi.org/10.1039/c7cp02101a

2017-01-01

Abstract:Snapshots of the intermediate conformation of Trp-cage at various simulation times using AMBER03, EPB03, AMBER12SB, and EPB12SB. Here, the N terminal is always on the top.
Improved Parameterization of Amine-Carboxyate and Amine-Phosphate Interactions for Molecular Dynamics Simulations Using the CHARMM and AMBER Force Fields

Jejoong Yoo,Aleksei Aksimentiev

DOI: https://doi.org/10.1016/j.bpj.2015.11.3457

2015-10-08

Abstract:Over the past decades, molecular dynamics (MD) simulations of biomolecules have become a mainstream biophysics technique. As the length and time scales amenable to the MD method increase, shortcomings of the empirical force fields---which have been developed and validated using relatively short simulations of small molecules---become apparent. One common artifact is aggregation of water-soluble biomolecules driven by artificially strong charge--charge interactions. Here, we report a systematic refinement of Lennard-Jones parameters (NBFIX) describing amine--carboxylate and amine--phosphate interactions, which brings MD simulations of basic peptide-mediated nucleic acids assemblies and lipid bilayer membranes in better agreement with experiment. As our refinement neither affects the existing parameterization of bonded interaction nor does it alter the solvation free energies, it improves realism of an MD simulation without introducing additional artifacts.

Biological Physics,Chemical Physics,Biomolecules
Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

JL Banks,GA Kaminski,RH Zhou,DT Mainz,BJ Berne,RA Friesner

DOI: https://doi.org/10.1063/1.478043

1999-01-01

Abstract:We have developed a polarizable force field for peptides, using all-atom OPLS (OPLS-AA) nonelectrostatic terms and electrostatics based on a fluctuating charge model and fit to ab initio calculations of polarization responses. We discuss the fitting procedure, and specific techniques we have developed that are necessary in order to obtain an accurate, stable model. Our model is comparable to the best existing molecular mechanics force fields in reproducing quantum-chemical peptide energetics. It also accurately reproduces many-body effects in many cases. We believe that straightforward extensions of our linear-response electrostatic model will significantly improve the accuracy for those cases that the present model does not adequately address.
A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations

Yong Duan,Chun Wu,Shibasish Chowdhury,Mathew C. Lee,Guoming Xiong,Wei Zhang,Rong Yang,Piotr Cieplak,Ray Luo,Taisung Lee,James Caldwell,Junmei Wang,Peter Kollman

DOI: https://doi.org/10.1002/jcc.10349

2003-10-01

Journal of Computational Chemistry

Abstract:Molecular mechanics models have been applied extensively to study the dynamics of proteins and nucleic acids. Here we report the development of a third-generation point-charge all-atom force field for proteins. Following the earlier approach of Cornell et al., the charge set was obtained by fitting to the electrostatic potentials of dipeptides calculated using B3LYP/cc-pVTZ//HF/6-31G** quantum mechanical methods. The main-chain torsion parameters were obtained by fitting to the energy profiles of Ace-Ala-Nme and Ace-Gly-Nme di-peptides calculated using MP2/cc-pVTZ//HF/6-31G** quantum mechanical methods. All other parameters were taken from the existing AMBER data base. The major departure from previous force fields is that all quantum mechanical calculations were done in the condensed phase with continuum solvent models and an effective dielectric constant of epsilon = 4. We anticipate that this force field parameter set will address certain critical short comings of previous force fields in condensed-phase simulations of proteins. Initial tests on peptides demonstrated a high-degree of similarity between the calculated and the statistically measured Ramanchandran maps for both Ace-Gly-Nme and Ace-Ala-Nme di-peptides. Some highlights of our results include (1) well-preserved balance between the extended and helical region distributions, and (2) favorable type-II poly-proline helical region in agreement with recent experiments. Backward compatibility between the new and Cornell et al. charge sets, as judged by overall agreement between dipole moments, allows a smooth transition to the new force field in the area of ligand-binding calculations. Test simulations on a large set of proteins are also discussed.

chemistry, multidisciplinary

CHARMM36m: an improved force field for folded and intrinsically disordered proteins

Comparison and Evaluation of Force Fields for Intrinsically Disordered Proteins

Significantly Improved Protein Folding Thermodynamics Using a Dispersion-Corrected Water Model and a New Residue-Specific Force Field

Polarizable Force Field of Intrinsically Disordered Proteins with CMAP and Reweighting Optimization

A Test of AMBER Force Fields in Predicting the Secondary Structure of Α-Helical and Β-Hairpin Peptides

Toward a general neural network force field for protein simulations: Refining the intramolecular interaction in protein

High quality NMR structures: a new force field with implicit water and membrane solvation for Xplor-NIH

Environment-Specific Force Field for Intrinsically Disordered and Ordered Proteins

Likely Overstabilization of Charge-Charge Interactions in CHARMM36m(w): A Case for a99SB-disp Water

Improved side-chain torsion potentials for the Amber ff99SB protein force field

Further optimization of a hybrid united-atom and coarse-grained force field for folding simulations: Improved backbone hydration and interactions between charged side chains.

Developing a molecular dynamics force field for both folded and disordered protein states

Residue-Speci fi c Force Field Based on the Protein Coil Library

A Force Balanced Fragmentation Method for Ab Initio Molecular Dynamic Simulation of Protein

New-Generation Amber United-Atom Force Field

CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

Universal Implementation of a Residue-Specific Force Field Based on CMAP Potentials and Free Energy Decomposition

Direct Folding Simulation of Helical Proteins Using an Effective Polarizable Bond Force Field

Improved Parameterization of Amine-Carboxyate and Amine-Phosphate Interactions for Molecular Dynamics Simulations Using the CHARMM and AMBER Force Fields

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations