1
|
Wang L, Schauperl M, Mobley DL, Bayly C, Gilson MK. A Fast, Convenient, Polarizable Electrostatic Model for Molecular Dynamics. J Chem Theory Comput 2024; 20:1293-1305. [PMID: 38240687 PMCID: PMC10867846 DOI: 10.1021/acs.jctc.3c01171] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2024]
Abstract
We present an efficient polarizable electrostatic model, utilizing typed, atom-centered polarizabilities and the fast direct approximation, designed for efficient use in molecular dynamics (MD) simulations. The model provides two convenient approaches for assigning partial charges in the context of atomic polarizabilities. One is a generalization of RESP, called RESP-dPol, and the other, AM1-BCC-dPol, is an adaptation of the widely used AM1-BCC method. Both are designed to accurately replicate gas-phase quantum mechanical electrostatic potentials. Benchmarks of this polarizable electrostatic model against gas-phase dipole moments, molecular polarizabilities, bulk liquid densities, and static dielectric constants of organic liquids show good agreement with the reference values. Of note, the model yields markedly more accurate dielectric constants of organic liquids, relative to a matched nonpolarizable force field. MD simulations with this method, which is currently parametrized for molecules containing elements C, N, O, and H, run only about 3.6-fold slower than fixed charge force fields, while simulations with the self-consistent mutual polarization average 4.5-fold slower. Our results suggest that RESP-dPol and AM1-BCC-dPol afford improved accuracy relative to fixed charge force fields and are good starting points for developing general, affordable, and transferable polarizable force fields. The software implementing these approaches has been designed to utilize the force field fitting frameworks developed and maintained by the Open Force Field Initiative, setting the stage for further exploration of this approach to polarizable force field development.
Collapse
Affiliation(s)
- Liangyue Wang
- Department
of Chemistry and Biochemistry, University
of California, San Diego, California 92093, United States
| | - Michael Schauperl
- HotSpot
Therapeutics, Inc., Boston, Massachusetts 02210, United States
| | - David L. Mobley
- Department
of Pharmaceutical Sciences, University of
California, Irvine, California 92697, United States
| | - Christopher Bayly
- OpenEye
Scientific, Cadence Molecular Sciences, Santa Fe, New Mexico 87508, United States
| | - Michael K. Gilson
- Skaggs
School of Pharmacy and Pharmaceutical Sciences, University of California, San
Diego, California 92093, United States
| |
Collapse
|
2
|
Deng J, Cui Q. Electronic Polarization Is Essential for the Stabilization and Dynamics of Buried Ion Pairs in Staphylococcal Nuclease Mutants. J Am Chem Soc 2022; 144:4594-4610. [PMID: 35239338 PMCID: PMC9616648 DOI: 10.1021/jacs.2c00312] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Buried charged residues play important roles in the modulation of protein stabilities and conformational dynamics and make crucial contributions to protein functions. Considering the generally nonpolar nature of protein interior, a key question concerns the contribution of electronic polarization to the stabilization and properties of buried charges. We answer this question by conducting free energy simulations using the latest polarizable CHARMM force field based on Drude oscillators for a series of Staphylococcal nuclease mutants that involve a buried Glu-Lys pair in different titration states and orientations. While a nonpolarizable model suggests that the ionized form of the buried Glu-Lys pair is more than 40 kcal/mol less stable than the charge-neutral form, the two titration states are comparable in stability when electronic polarization is included explicitly, a result better reconcilable with available experimental data. Analysis of free energy components suggests that additional stabilization of the ionized Glu-Lys pair has contributions from both the enhanced salt-bridge strength and stronger interaction between the ion-pair and surrounding protein residues and penetrated water. Despite the stronger direct interaction between Glu and Lys, the ion-pair exhibits considerably larger and faster structural fluctuations when polarization is included, due to compensation of interactions in the cavity. Collectively, observations from this work provide compelling evidence that electronic polarization is essential to the stability, hydration, dynamics, and therefore function of buried charges in proteins. Therefore, our study advocates for the explicit consideration of electronic polarization for mechanistic and engineering studies that implicate buried charged residues, such as enzymes and ion transporters.
Collapse
Affiliation(s)
- Jiahua Deng
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States.,Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States.,Department of Biomedical Engineering, Boston University, 44 Cummington Mall, Boston, Massachusetts 02215, United States
| |
Collapse
|
3
|
Pan Z, Huang J, Zhuang W. Protein-Ligand Binding Molecular Details Revealed by Terahertz Optical Kerr Spectroscopy: A Simulation Study. JACS AU 2021; 1:1788-1797. [PMID: 34723281 PMCID: PMC8549111 DOI: 10.1021/jacsau.1c00356] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/15/2021] [Indexed: 06/13/2023]
Abstract
Picosecond fast motions and their involvement in the biochemical processes such as protein-ligand binding has engaged significant attention. Terahertz optical Kerr spectroscopy (OKE) has the superior potential to probe these fast motions directly. Application of OKE in protein-ligand binding study is, however, limited by the difficulty of quantitative atomistic interpretation, and the calculation of Kerr spectrum for entire solvated protein complex was considered not yet feasible, due to the lack of one consistent polarizable model for both configuration sampling and polarizability calculation. Here, we analyzed the biochemical relevance of OKE to the lysozyme-triacetylchitotriose binding based on the first OKE simulation using one consistent Drude polarizable model. An analytical multipole and induced dipole scheme was employed to calculate the off-diagonal Drude polarizability more efficiently and accurately. Further theoretical analysis revealed how the subtle twisting and stiffening of aromatic protein residues' spatial arrangement as well as the confinement of small water clusters between ligand and protein cavity due to the ligand binding can be examined using Kerr spectroscopy. Comparison between the signals of bound complex and that of uncorrelated protein/ligand demonstrated that binding action alone has reflection in the OKE spectrum. Our study indicated OKE as a powerful terahertz probe for protein-ligand binding chemistry and dynamics.
Collapse
Affiliation(s)
- Zhijun Pan
- State
Key Laboratory of Structural Chemistry, Fujian Institute of Research
on the Structure of Matter, Chinese Academy
of Sciences, Fuzhou, Fujian 35000, China
- Shenzhen
Bay Laboratory, Guangming
District, Shenzhen 518107, China
| | - Jing Huang
- Westlake
Laboratory of Life Sciences and Biomedicine, Key Laboratory of Structural
Biology of Zhejiang Province, School of Life Sciences, Westlake University, Hangzhou, Zhejiang 310024, China
| | - Wei Zhuang
- State
Key Laboratory of Structural Chemistry, Fujian Institute of Research
on the Structure of Matter, Chinese Academy
of Sciences, Fuzhou, Fujian 35000, China
- Institute
of Urban Environment, Chinese Academy of
Sciences, XiaMen, Fujian 361021, China
| |
Collapse
|
4
|
Xu T, Cao K, Wang C, Yin S. The effect of asymmetric external reorganization energy on electron and hole transport in organic semiconductors. Phys Chem Chem Phys 2021; 23:15236-15244. [PMID: 34235520 DOI: 10.1039/d0cp06686a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Understanding the relationship between charge mobilities and the molecular stacking structures of π-conjugated organic semiconducting materials is essential for their development. In this study, a quantum mechanics (QM)-derived state-specific polarizable force field (SS-PFF) is applied to explicitly estimate the external reorganization energies (λext) during electron transfer (ET) or hole transfer (HT) processes using our recently proposed two-point model (J. Phys. Chem. A, 2018, 122, 8957-8964). Different from the Marcus two-sphere model, the application of the explicit two-point model produces a notably asymmetric λext for ET and HT processes in oligoacene crystals. For the same charge transfer channels, the λext of ET is 7-10 times higher than that of HT, which results in a larger intrinsic hole mobility. This successfully rationalizes why acenes are prone to be p-type conducting materials. Perfluorination can change the polarity of the molecular surface electrostatic potential (ESP). Thus, perfluorination is a possible approach to reduce the external reorganization energies for ET reactions, which can be used to tune the order of λext of ET and HT processes. The two-point model with the SS-PFF, therefore, opens the door to revisit the intrinsic mobilities of electrons and holes in organic semiconductors.
Collapse
Affiliation(s)
- Tao Xu
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an City 710119, People's Republic of China.
| | | | | | | |
Collapse
|
5
|
Corrigan RA, Qi G, Thiel AC, Lynn JR, Walker BD, Casavant TL, Lagardere L, Piquemal JP, Ponder JW, Ren P, Schnieders MJ. Implicit Solvents for the Polarizable Atomic Multipole AMOEBA Force Field. J Chem Theory Comput 2021; 17:2323-2341. [PMID: 33769814 DOI: 10.1021/acs.jctc.0c01286] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Computational protein design, ab initio protein/RNA folding, and protein-ligand screening can be too computationally demanding for explicit treatment of solvent. For these applications, implicit solvent offers a compelling alternative, which we describe here for the polarizable atomic multipole AMOEBA force field based on three treatments of continuum electrostatics: numerical solutions to the nonlinear and linearized versions of the Poisson-Boltzmann equation (PBE), the domain-decomposition conductor-like screening model (ddCOSMO) approximation to the PBE, and the analytic generalized Kirkwood (GK) approximation. The continuum electrostatics models are combined with a nonpolar estimator based on novel cavitation and dispersion terms. Electrostatic model parameters are numerically optimized using a least-squares style target function based on a library of 103 small-molecule solvation free energy differences. Mean signed errors for the adaptive Poisson-Boltzmann solver (APBS), ddCOSMO, and GK models are 0.05, 0.00, and 0.00 kcal/mol, respectively, while the mean unsigned errors are 0.70, 0.63, and 0.58 kcal/mol, respectively. Validation of the electrostatic response of the resulting implicit solvents, which are available in the Tinker (or Tinker-HP), OpenMM, and Force Field X software packages, is based on comparisons to explicit solvent simulations for a series of proteins and nucleic acids. Overall, the emergence of performative implicit solvent models for polarizable force fields opens the door to their use for folding and design applications.
Collapse
Affiliation(s)
- Rae A Corrigan
- Roy J Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, United States
| | - Guowei Qi
- Department of Biochemistry, University of Iowa, Iowa City, Iowa 52242, United States
| | - Andrew C Thiel
- Roy J Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, United States
| | - Jack R Lynn
- Roy J Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, United States
| | - Brandon D Walker
- Department of Biomedical Engineering, University of Texas in Austin, Austin, Texas 78712, United States
| | - Thomas L Casavant
- Roy J Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, United States
| | - Louis Lagardere
- Department of Chemistry, Sorbonne Université, F-75005 Paris, France
| | | | - Jay W Ponder
- Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri 63130, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, University of Texas in Austin, Austin, Texas 78712, United States
| | - Michael J Schnieders
- Roy J Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, United States.,Department of Biochemistry, University of Iowa, Iowa City, Iowa 52242, United States
| |
Collapse
|
6
|
Rupakheti C, Lamoureux G, MacKerell AD, Roux B. Statistical mechanics of polarizable force fields based on classical Drude oscillators with dynamical propagation by the dual-thermostat extended Lagrangian. J Chem Phys 2020; 153:114108. [PMID: 32962358 PMCID: PMC7656322 DOI: 10.1063/5.0019987] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2020] [Accepted: 08/13/2020] [Indexed: 12/11/2022] Open
Abstract
Polarizable force fields based on classical Drude oscillators offer a practical and computationally efficient avenue to carry out molecular dynamics (MD) simulations of large biomolecular systems. To treat the polarizable electronic degrees of freedom, the Drude model introduces a virtual charged particle that is attached to its parent nucleus via a harmonic spring. Traditionally, the need to relax the electronic degrees of freedom for each fixed set of nuclear coordinates is achieved by performing an iterative self-consistent field (SCF) calculation to satisfy a selected tolerance. This is a computationally demanding procedure that can increase the computational cost of MD simulations by nearly one order of magnitude. To avoid the costly SCF procedure, a small mass is assigned to the Drude particles, which are then propagated as dynamic variables during the simulations via a dual-thermostat extended Lagrangian algorithm. To help clarify the significance of the dual-thermostat extended Lagrangian propagation in the context of the polarizable force field based on classical Drude oscillators, the statistical mechanics of a dual-temperature canonical ensemble is formulated. The conditions for dynamically maintaining the dual-temperature properties in the case of the classical Drude oscillator are analyzed using the generalized Langevin equation.
Collapse
Affiliation(s)
- Chetan Rupakheti
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois 60637, USA
| | - Guillaume Lamoureux
- Department of Chemistry and Center for Computational and Integrative Biology, Rutgers University, Camden, New Jersey 08102, USA
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland 21201, USA
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
7
|
Marcellini M, Nguyen MH, Martin M, Hologne M, Walker O. Accurate Prediction of Protein NMR Spin Relaxation by Means of Polarizable Force Fields. Application to Strongly Anisotropic Rotational Diffusion. J Phys Chem B 2020; 124:5103-5112. [DOI: 10.1021/acs.jpcb.0c01922] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Moreno Marcellini
- Institut des Sciences Analytiques (ISA), Univ Lyon, CNRS, UMR5280, Université Claude Bernard Lyon1, Lyon, France
| | - Minh-Ha Nguyen
- Institut des Sciences Analytiques (ISA), Univ Lyon, CNRS, UMR5280, Université Claude Bernard Lyon1, Lyon, France
| | - Marie Martin
- Institut des Sciences Analytiques (ISA), Univ Lyon, CNRS, UMR5280, Université Claude Bernard Lyon1, Lyon, France
| | - Maggy Hologne
- Institut des Sciences Analytiques (ISA), Univ Lyon, CNRS, UMR5280, Université Claude Bernard Lyon1, Lyon, France
| | - Olivier Walker
- Institut des Sciences Analytiques (ISA), Univ Lyon, CNRS, UMR5280, Université Claude Bernard Lyon1, Lyon, France
| |
Collapse
|
8
|
Liu C, Piquemal JP, Ren P. Implementation of Geometry-Dependent Charge Flux into the Polarizable AMOEBA+ Potential. J Phys Chem Lett 2020; 11:419-426. [PMID: 31865706 PMCID: PMC7384396 DOI: 10.1021/acs.jpclett.9b03489] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Molecular dynamics (MD) simulations employing classical force fields (FFs) have been widely used to model molecular systems. The important ingredient of the current FFs, atomic charge, remains fixed during MD simulations despite the atomic environment or local geometry changes. This approximation hinders the transferability of the potential being used in multiple phases. Here we implement a geometry-dependent charge flux (GDCF) model into the multipole-based AMOEBA+ polarizable potential. The CF in the current work explicitly depends on the local geometry (bond and angle) of the molecule. To our knowledge, this is the first study that derives energy and force expressions due to GDCF in a multipole-based polarizable FF framework. Due to the inclusion of GDCF, the AMOEBA+ water model is noticeably improved in terms of describing the monomer properties, cluster binding/interaction energy, and a variety of liquid properties, including the infrared spectra that previous flexible water models were not able to capture.
Collapse
Affiliation(s)
- Chengwen Liu
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States
| | - Jean-Philip Piquemal
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States
- Laboratoire de Chimie Théorique , Sorbonne Université, UMR7616 CNRS , 75252 Paris , France
- Institut Universitaire de France , 75005 , Paris , France
| | - Pengyu Ren
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States
| |
Collapse
|
9
|
Pandey P, Aytenfisu AH, MacKerell AD, Mallajosyula SS. Drude Polarizable Force Field Parametrization of Carboxylate and N-Acetyl Amine Carbohydrate Derivatives. J Chem Theory Comput 2019; 15:4982-5000. [PMID: 31411469 PMCID: PMC6852669 DOI: 10.1021/acs.jctc.9b00327] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
In this work, we report the development of Drude polarizable force field parameters for the carboxylate and N-acetyl amine derivatives, extending the functionality of the existing Drude polarizable carbohydrate force field. The force field parameters have been developed in a hierarchical manner, reproducing the quantum mechanical gas-phase properties of small model compounds representing the key functional group in the carbohydrate derivatives, including optimization of the electrostatic and bonded parameters. The optimized parameters were then used to generate the models for carboxylate and N-acetyl amine carbohydrate derivatives. The transferred parameters were further tested and optimized to reproduce crystal geometries and J-coupling data from nuclear magnetic resonance experiments. The parameter development resulted in the incorporation of d-glucuronate, l-iduronate, N-acetyl-d-glucosamine (GlcNAc), and N-acetyl-d-galactosamine (GalNAc) sugars into the Drude polarizable force field. The parameters developed in this study were then applied to study the conformational properties of glycosaminoglycan polymer hyaluronan, composed of d-glucuronate and N-acetyl-d-glucosamine, in aqueous solution. Upon comparing the results from the additive and polarizable simulations, it was found that the inclusion of polarization improved the description of the electrostatic interactions observed in hyaluronan, resulting in enhanced conformational flexibility. The developed Drude polarizable force field parameters in conjunction with the remainder of the Drude polarizable force field parameters can be used for future studies involving carbohydrates and their conjugates in complex, heterogeneous systems.
Collapse
Affiliation(s)
| | - Asaminew H Aytenfisu
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , 20 Penn Street , Baltimore , Maryland 21201 , United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , 20 Penn Street , Baltimore , Maryland 21201 , United States
| | | |
Collapse
|
10
|
Bedrov D, Piquemal JP, Borodin O, MacKerell AD, Roux B, Schröder C. Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields. Chem Rev 2019; 119:7940-7995. [PMID: 31141351 PMCID: PMC6620131 DOI: 10.1021/acs.chemrev.8b00763] [Citation(s) in RCA: 284] [Impact Index Per Article: 56.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2018] [Indexed: 11/30/2022]
Abstract
Many applications in chemistry, biology, and energy storage/conversion research rely on molecular simulations to provide fundamental insight into structural and transport properties of materials with high ionic concentrations. Whether the system is comprised entirely of ions, like ionic liquids, or is a mixture of a polar solvent with a salt, e.g., liquid electrolytes for battery applications, the presence of ions in these materials results in strong local electric fields polarizing solvent molecules and large ions. To predict properties of such systems from molecular simulations often requires either explicit or mean-field inclusion of the influence of polarization on electrostatic interactions. In this manuscript, we review the pros and cons of different treatments of polarization ranging from the mean-field approaches to the most popular explicit polarization models in molecular dynamics simulations of ionic materials. For each method, we discuss their advantages and disadvantages and emphasize key assumptions as well as their adjustable parameters. Strategies for the development of polarizable models are presented with a specific focus on extracting atomic polarizabilities. Finally, we compare simulations using polarizable and nonpolarizable models for several classes of ionic systems, discussing the underlying physics that each approach includes or ignores, implications for implementation and computational efficiency, and the accuracy of properties predicted by these methods compared to experiments.
Collapse
Affiliation(s)
- Dmitry Bedrov
- Department
of Materials Science & Engineering, University of Utah, 122 South Central Campus Drive, Room 304, Salt Lake City, Utah 84112, United States
| | - Jean-Philip Piquemal
- Laboratoire
de Chimie Théorique, Sorbonne Université,
UMR 7616 CNRS, CC137, 4 Place Jussieu, Tour 12-13, 4ème étage, 75252 Paris Cedex 05, France
- Institut
Universitaire de France, 75005, Paris Cedex 05, France
- Department
of Biomedical Engineering, The University
of Texas at Austin, Austin, Texas 78712, United States
| | - Oleg Borodin
- Electrochemistry
Branch, Sensors and Electron Devices Directorate, Army Research Laboratory, 2800 Powder Mill Road, Adelphi, Maryland 20703, United
States
| | - Alexander D. MacKerell
- Department
of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland 21201, United
States
| | - Benoît Roux
- Department
of Biochemistry and Molecular Biology, Gordon Center for Integrative
Science, University of Chicago, 929 57th Street, Chicago, Illinois 60637, United States
| | - Christian Schröder
- Department
of Computational Biological Chemistry, University
of Vienna, Währinger Strasse 17, A-1090 Vienna, Austria
| |
Collapse
|
11
|
Nocito D, Beran GJO. Reduced computational cost of polarizable force fields by a modification of the always stable predictor-corrector. J Chem Phys 2019; 150:151103. [PMID: 31005093 DOI: 10.1063/1.5092133] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Classical polarizable force fields effectively incorporate the dynamic response of the electronic charge distributions into molecular dynamics simulations, but they do so at a significant increase in computational cost compared to simpler models. Here, we demonstrate how one can improve the stability of a polarizable force field molecular dynamics simulation or accelerate the evaluation of self-consistent polarization via a simple extension of the predictor in the always stable predictor-corrector method. Specifically, increasing the number of prior steps used in the predictor from 6 to 16 reduces the energy drift by an order of magnitude. Alternatively, for a given level of energy drift, the induced dipoles can be obtained ∼20% faster due to the reduced number of self-consistent field iterations required to maintain energetic stability. The extended-history predictor is straightforward to implement and involves minimal computational overhead.
Collapse
Affiliation(s)
- Dominique Nocito
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Gregory J O Beran
- Department of Chemistry, University of California, Riverside, California 92521, USA
| |
Collapse
|
12
|
Jing Z, Liu C, Cheng SY, Qi R, Walker BD, Piquemal JP, Ren P. Polarizable Force Fields for Biomolecular Simulations: Recent Advances and Applications. Annu Rev Biophys 2019; 48:371-394. [PMID: 30916997 DOI: 10.1146/annurev-biophys-070317-033349] [Citation(s) in RCA: 217] [Impact Index Per Article: 43.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Realistic modeling of biomolecular systems requires an accurate treatment of electrostatics, including electronic polarization. Due to recent advances in physical models, simulation algorithms, and computing hardware, biomolecular simulations with advanced force fields at biologically relevant timescales are becoming increasingly promising. These advancements have not only led to new biophysical insights but also afforded opportunities to advance our understanding of fundamental intermolecular forces. This article describes the recent advances and applications, as well as future directions, of polarizable force fields in biomolecular simulations.
Collapse
Affiliation(s)
- Zhifeng Jing
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Sara Y Cheng
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Rui Qi
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Brandon D Walker
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Jean-Philip Piquemal
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA; .,Sorbonne Université, CNRS, Laboratoire de Chimie Theórique, 75252 Paris CEDEX 05, France.,Institut Universitaire de France, 75005 Paris, France
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| |
Collapse
|
13
|
Abstract
Molecular dynamics (MD) simulations have been widely applied to computer-aided drug design (CADD). While MD has been used in a variety of applications such as free energy perturbation and long-time simulations, the accuracy of the results from those methods depends strongly on the force field used. Force fields for small molecules are crucial, as they not only serve as building blocks for developing force fields for larger biomolecules but also act as model compounds that will be transferred to ligands used in CADD. Currently, a wide range of small molecule force fields based on additive or nonpolarizable models have been developed. While these nonpolarizable force fields can produce reasonable estimations of physical properties and have shown success in a variety of systems, there is still room for improvements due to inherent limitations in these models including the lack of an electronic polarization response. For this reason, incorporating polarization effects into the energy function underlying a force field is believed to be an important step forward, giving rise to the development of polarizable force fields. Recent simulations of biological systems have indicated that polarizable force fields are able to provide a better physical representation of intermolecular interactions and, in many cases, better agreement with experimental properties than nonpolarizable, additive force fields. Therefore, this chapter focuses on the development of small molecule force fields with emphasis on polarizable models. It begins with a brief introduction on the importance of small molecule force fields and their evolution from additive to polarizable force fields. Emphasis is placed on the additive CHARMM General Force Field and the polarizable force field based on the classical Drude oscillator. The theory for the Drude polarizable force field and results for small molecules are presented showing their improvements over the additive model. The potential importance of polarization for their application in a wide range of biological systems including CADD is then discussed.
Collapse
Affiliation(s)
- Fang-Yu Lin
- Department of Pharmaceutical Sciences, Computer-Aided Drug Design Center, School of Pharmacy, University of Maryland, Baltimore, MD, USA
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, Computer-Aided Drug Design Center, School of Pharmacy, University of Maryland, Baltimore, MD, USA.
| |
Collapse
|
14
|
Xu T, Wang W, Yin S. Explicit Method To Evaluate the External Reorganization Energy of Charge-Transfer Reactions in Oligoacene Crystals Using the State-Specific Polarizable Force Field. J Phys Chem A 2018; 122:8957-8964. [DOI: 10.1021/acs.jpca.8b08998] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Tao Xu
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi’an City 710119, People’s Republic of China
| | - Wenliang Wang
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi’an City 710119, People’s Republic of China
| | - Shiwei Yin
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi’an City 710119, People’s Republic of China
| |
Collapse
|
15
|
Rackers JA, Wang Z, Lu C, Laury ML, Lagardère L, Schnieders MJ, Piquemal JP, Ren P, Ponder JW. Tinker 8: Software Tools for Molecular Design. J Chem Theory Comput 2018; 14:5273-5289. [PMID: 30176213 PMCID: PMC6335969 DOI: 10.1021/acs.jctc.8b00529] [Citation(s) in RCA: 290] [Impact Index Per Article: 48.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The Tinker software, currently released as version 8, is a modular molecular mechanics and dynamics package written primarily in a standard, easily portable dialect of Fortran 95 with OpenMP extensions. It supports a wide variety of force fields, including polarizable models such as the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field. The package runs on Linux, macOS, and Windows systems. In addition to canonical Tinker, there are branches, Tinker-HP and Tinker-OpenMM, designed for use on message passing interface (MPI) parallel distributed memory supercomputers and state-of-the-art graphical processing units (GPUs), respectively. The Tinker suite also includes a tightly integrated Java-based graphical user interface called Force Field Explorer (FFE), which provides molecular visualization capabilities as well as the ability to launch and control Tinker calculations.
Collapse
Affiliation(s)
- Joshua A. Rackers
- Program in Computational & Molecular Biophysics, Washington University School of Medicine, Saint Louis, Missouri 63110, United States
| | - Zhi Wang
- Department of Chemistry, Washington University in Saint Louis, Saint Louis, Missouri 63130, United States
| | - Chao Lu
- Department of Chemistry, Washington University in Saint Louis, Saint Louis, Missouri 63130, United States
| | - Marie L. Laury
- Department of Chemistry, Washington University in Saint Louis, Saint Louis, Missouri 63130, United States
| | - Louis Lagardère
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC Paris 06, UMR 7616, case courrier 137, 4 place Jussieu, F-75005, Paris, France
| | - Michael J. Schnieders
- Department of Biomedical Engineering, The University of Iowa, Iowa City, IA 52242, United States
| | - Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC Paris 06, UMR 7616, case courrier 137, 4 place Jussieu, F-75005, Paris, France
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Jay W. Ponder
- Program in Computational & Molecular Biophysics, Washington University School of Medicine, Saint Louis, Missouri 63110, United States
- Department of Chemistry, Washington University in Saint Louis, Saint Louis, Missouri 63130, United States
| |
Collapse
|
16
|
Huang J, Simmonett AC, Pickard FC, MacKerell AD, Brooks BR. Mapping the Drude polarizable force field onto a multipole and induced dipole model. J Chem Phys 2018; 147:161702. [PMID: 29096511 DOI: 10.1063/1.4984113] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
The induced dipole and the classical Drude oscillator represent two major approaches for the explicit inclusion of electronic polarizability into force field-based molecular modeling and simulations. In this work, we explore the equivalency of these two models by comparing condensed phase properties computed using the Drude force field and a multipole and induced dipole (MPID) model. Presented is an approach to map the electrostatic model optimized in the context of the Drude force field onto the MPID model. Condensed phase simulations on water and 15 small model compounds show that without any reparametrization, the MPID model yields properties similar to the Drude force field with both models yielding satisfactory reproduction of a range of experimental values and quantum mechanical data. Our results illustrate that the Drude oscillator model and the point induced dipole model are different representations of essentially the same physical model. However, results indicate the presence of small differences between the use of atomic multipoles and off-center charge sites. Additionally, results on the use of dispersion particle mesh Ewald further support its utility for treating long-range Lennard Jones dispersion contributions in the context of polarizable force fields. The main motivation in demonstrating the transferability of parameters between the Drude and MPID models is that the more than 15 years of development of the Drude polarizable force field can now be used with MPID formalism without the need for dual-thermostat integrators nor self-consistent iterations. This opens up a wide range of new methodological opportunities for polarizable models.
Collapse
Affiliation(s)
- Jing Huang
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, 20 Penn St., Baltimore, Maryland 21201, USA
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, 5635 Fishers Lane, Rockville, Maryland 20852, USA
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, 5635 Fishers Lane, Rockville, Maryland 20852, USA
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, 20 Penn St., Baltimore, Maryland 21201, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, 5635 Fishers Lane, Rockville, Maryland 20852, USA
| |
Collapse
|
17
|
Aviat F, Lagardère L, Piquemal JP. The truncated conjugate gradient (TCG), a non-iterative/fixed-cost strategy for computing polarization in molecular dynamics: Fast evaluation of analytical forces. J Chem Phys 2018; 147:161724. [PMID: 29096518 DOI: 10.1063/1.4985911] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In a recent paper [F. Aviat et al., J. Chem. Theory Comput. 13, 180-190 (2017)], we proposed the Truncated Conjugate Gradient (TCG) approach to compute the polarization energy and forces in polarizable molecular simulations. The method consists in truncating the conjugate gradient algorithm at a fixed predetermined order leading to a fixed computational cost and can thus be considered "non-iterative." This gives the possibility to derive analytical forces avoiding the usual energy conservation (i.e., drifts) issues occurring with iterative approaches. A key point concerns the evaluation of the analytical gradients, which is more complex than that with a usual solver. In this paper, after reviewing the present state of the art of polarization solvers, we detail a viable strategy for the efficient implementation of the TCG calculation. The complete cost of the approach is then measured as it is tested using a multi-time step scheme and compared to timings using usual iterative approaches. We show that the TCG methods are more efficient than traditional techniques, making it a method of choice for future long molecular dynamics simulations using polarizable force fields where energy conservation matters. We detail the various steps required for the implementation of the complete method by software developers.
Collapse
Affiliation(s)
- Félix Aviat
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC Université Paris 06, UMR 7616, F-75005 Paris, France
| | - Louis Lagardère
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC Université Paris 06, UMR 7616, F-75005 Paris, France
| | - Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC Université Paris 06, UMR 7616, F-75005 Paris, France
| |
Collapse
|
18
|
New developments in force fields for biomolecular simulations. Curr Opin Struct Biol 2018; 49:129-138. [DOI: 10.1016/j.sbi.2018.02.002] [Citation(s) in RCA: 136] [Impact Index Per Article: 22.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2017] [Revised: 01/30/2018] [Accepted: 02/04/2018] [Indexed: 11/18/2022]
|
19
|
Albaugh A, Head-Gordon T, Niklasson AMN. Higher-Order Extended Lagrangian Born–Oppenheimer Molecular Dynamics for Classical Polarizable Models. J Chem Theory Comput 2018; 14:499-511. [DOI: 10.1021/acs.jctc.7b01041] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | | | - Anders M. N. Niklasson
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
20
|
Demerdash O, Wang L, Head‐Gordon T. Advanced models for water simulations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1355] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Omar Demerdash
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley CA USA
- Department of Chemistry University of California Berkeley CA USA
| | - Lee‐Ping Wang
- Department of Chemistry University of California, Davis Davis CA USA
| | - Teresa Head‐Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley CA USA
- Department of Chemistry University of California Berkeley CA USA
- Department of Bioengineering University of California Berkeley CA USA
- Department of Chemical and Biomolecular Engineering University of California Berkeley CA USA
| |
Collapse
|
21
|
Albaugh A, Head-Gordon T. A New Method for Treating Drude Polarization in Classical Molecular Simulation. J Chem Theory Comput 2017; 13:5207-5216. [PMID: 28965397 DOI: 10.1021/acs.jctc.7b00838] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
With polarization becoming an increasingly common feature in classical molecular simulation, it is important to develop methods that can efficiently and accurately evaluate the many-body polarization solution. In this work, we expand the theoretical framework of our inertial extended Langrangian, self-consistent field iteration-free method (iEL/0-SCF), introduced for point induced dipoles, to the polarization model of a Drude oscillator. When applied to the polarizable simple point charge model (PSPC) for water, our iEL/0-SCF method for Drude polarization is as stable as a well-converged SCF solution and more stable than traditional extended Lagrangian (EL) approaches or EL formulations based on two temperature ensembles where Drude particles are kept "colder" than the real degrees of freedom. We show that the iEL/0-SCF method eliminates the need for mass repartitioning from parent atoms onto Drude particles, obeys system conservation of linear and angular momentum, and permits the extension of the integration time step of a basic molecular dynamics simulation to 6.0 fs for PSPC water.
Collapse
Affiliation(s)
- Alex Albaugh
- Departments of Chemical & Biomolecular Engineering, ‡Chemistry, and §Bioengineering, ∥Chemical Sciences Division, Lawrence Berkeley National Laboratory, University of California , Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Departments of Chemical & Biomolecular Engineering, ‡Chemistry, and §Bioengineering, ∥Chemical Sciences Division, Lawrence Berkeley National Laboratory, University of California , Berkeley, California 94720, United States
| |
Collapse
|
22
|
Albaugh A, Niklasson AMN, Head-Gordon T. Accurate Classical Polarization Solution with No Self-Consistent Field Iterations. J Phys Chem Lett 2017; 8:1714-1723. [PMID: 28350167 DOI: 10.1021/acs.jpclett.7b00450] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We present a new solution for classical polarization that does not require any self-consistent field iterations, the aspect of classical polarization that makes it computationally expensive. The new method builds upon our iEL/SCF Lagrangian scheme that defines a set of auxiliary induced dipoles whose original purpose was to serve as a time-reversible initial guess to the SCF solution of the set of real induced dipoles. In the new iEL/0-SCF approach the auxiliary dipoles now drive the time evolution of the real induced dipoles such that they stay close to the Born-Oppenheimer surface in order to achieve a truly SCF-less method. We show that the iEL/0-SCF exhibits no loss of simulation accuracy when analyzed across bulk water, low to high concentration salt solutions, and small solutes to large proteins in water. In addition, iEL/0-SCF offers significant computational savings over more expensive SCF calculations based on traditional 1 fs time step integration using symplectic integrators and is as fast as reversible reference system propagator algorithms with an outer 2 fs time step.
Collapse
Affiliation(s)
| | - Anders M N Niklasson
- Theoretical Division, Los Alamos National Laboratory , Los Alamos, New Mexico 87545, United States
| | - Teresa Head-Gordon
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, University of California , Berkeley, California 94720, United States
| |
Collapse
|
23
|
Nocito D, Beran GJO. Fast divide-and-conquer algorithm for evaluating polarization in classical force fields. J Chem Phys 2017; 146:114103. [DOI: 10.1063/1.4977981] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Dominique Nocito
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Gregory J. O. Beran
- Department of Chemistry, University of California, Riverside, California 92521, USA
| |
Collapse
|
24
|
Aviat F, Levitt A, Stamm B, Maday Y, Ren P, Ponder JW, Lagardère L, Piquemal JP. Truncated Conjugate Gradient: An Optimal Strategy for the Analytical Evaluation of the Many-Body Polarization Energy and Forces in Molecular Simulations. J Chem Theory Comput 2016; 13:180-190. [PMID: 28068773 PMCID: PMC5228058 DOI: 10.1021/acs.jctc.6b00981] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We
introduce a new class of methods, denoted as Truncated Conjugate
Gradient(TCG), to solve the many-body polarization energy and its
associated forces in molecular simulations (i.e. molecular dynamics
(MD) and Monte Carlo). The method consists in a fixed number of Conjugate
Gradient (CG) iterations. TCG approaches provide a scalable solution
to the polarization problem at a user-chosen cost and a corresponding
optimal accuracy. The optimality of the CG-method guarantees that
the number of the required matrix-vector products are reduced to a
minimum compared to other iterative methods. This family of methods
is non-empirical, fully adaptive, and provides analytical gradients,
avoiding therefore any energy drift in MD as compared to popular iterative
solvers. Besides speed, one great advantage of this class of approximate
methods is that their accuracy is systematically improvable. Indeed,
as the CG-method is a Krylov subspace method, the associated error
is monotonically reduced at each iteration. On top of that, two improvements
can be proposed at virtually no cost: (i) the use of preconditioners
can be employed, which leads to the Truncated Preconditioned Conjugate
Gradient (TPCG); (ii) since the residual of the final step of the
CG-method is available, one additional Picard fixed point iteration
(“peek”), equivalent to one step of Jacobi Over Relaxation
(JOR) with relaxation parameter ω, can be made at almost no
cost. This method is denoted by TCG-n(ω). Black-box adaptive
methods to find good choices of ω are provided and discussed.
Results show that TPCG-3(ω) is converged to high accuracy (a
few kcal/mol) for various types of systems including proteins and
highly charged systems at the fixed cost of four matrix-vector products:
three CG iterations plus the initial CG descent direction. Alternatively,
T(P)CG-2(ω) provides robust results at a reduced cost (three
matrix-vector products) and offers new perspectives for long polarizable
MD as a production algorithm. The T(P)CG-1(ω) level provides
less accurate solutions for inhomogeneous systems, but its applicability
to well-conditioned problems such as water is remarkable, with only
two matrix-vector product evaluations.
Collapse
Affiliation(s)
- Félix Aviat
- Laboratoire de Chimie Théorique, UPMC Université Paris 06, UMR 7617 , F-75005, Paris, France
| | - Antoine Levitt
- Inria Paris, F-75589 Paris Cedex 12, Université Paris-Est, CERMICS (ENPC) , Marne-la-Vallée, F-77455, France
| | - Benjamin Stamm
- MATHCCES, Department of Mathematics, RWTH Aachen University , Schinkelstraße 2, D-52062 Aachen, Germany.,Computational Biomedicine, Institute for Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum Jülich , Jülich, 52425, Germany
| | - Yvon Maday
- Laboratoire Jacques-Louis Lions, UPMC Université Paris 06, UMR 7598 , F-75005, Paris, France.,Institut Universitaire de France , Paris Cedex 05, 75231, France.,Division of Applied Maths, Brown University , Providence, Rhode Island 02912, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Jay W Ponder
- Department of Chemistry, Washington University in Saint Louis , Campus Box 1134, One Brookings Drive, Saint Louis, Missouri 63130, United States
| | - Louis Lagardère
- Laboratoire de Chimie Théorique, UPMC Université Paris 06, UMR 7617 , F-75005, Paris, France.,Institut du Calcul et de la Simulation, UPMC Université Paris 06 , F-75005, Paris, France
| | - Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, UPMC Université Paris 06, UMR 7617 , F-75005, Paris, France.,Institut Universitaire de France , Paris Cedex 05, 75231, France
| |
Collapse
|