1
|
Chung MKJ, Wang Z, Rackers JA, Ponder JW. Classical Exchange Polarization: An Anisotropic Variable Polarizability Model. J Phys Chem B 2022; 126:7579-7594. [PMID: 36166814 PMCID: PMC10868659 DOI: 10.1021/acs.jpcb.2c04237] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Polarizability, or the tendency of the electron distribution to distort under an electric field, often depends on the local chemical environment. For example, the polarizability of a chloride ion is larger in gas phase compared to a chloride ion solvated in water. This effect is due to the restriction the Pauli exclusion principle places on the allowed electron states. Because no two electrons can occupy the same state, when a highly polarizable atom comes in close contact with other atoms or molecules, the space of allowed states can dramatically decrease. This constraint suggests that an accurate molecular mechanics polarizability model should depend on the radial distance between neighboring atoms. This paper introduces a variable polarizability model within the framework of the HIPPO (Hydrogen-like Intermolecular Polarizable Potential) force field, by damping the polarizability as a function of the orbital overlap of two atoms. This effectively captures the quantum mechanical exchange polarization effects, without explicit utilization of antisymmetrized wave functions. We show that the variable polarizability model remarkably improves the two-body polarization energies and three-body energies of ion-ion and ion-water systems. Under this model, no manual tuning of atomic polarizabilities for monatomic ions is required; the gas-phase polarizability can be used because an appropriate damping function is able to correct the polarizability at short range.
Collapse
Affiliation(s)
- Moses K. J. Chung
- Medical Scientist Training Program, Washington University School of Medicine, Saint Louis, MO 63110, USA
- Department of Physics, Washington University in St. Louis, Saint Louis, MO 63130, USA
| | - Zhi Wang
- Department of Chemistry, Washington University in St. Louis, Saint Louis, MO 63130, USA
| | - Joshua A. Rackers
- Center for Computing Research, Sandia National Laboratories, Albuquerque, NM 87185, USA
| | - Jay W. Ponder
- Department of Chemistry, Washington University in St. Louis, Saint Louis, MO 63130, USA
- Department of Biochemistry & Molecular Biophysics, Washington University School of Medicine, Saint Louis, MO 63110, USA
| |
Collapse
|
2
|
Visscher KM, Geerke DP. Deriving a Polarizable Force Field for Biomolecular Building Blocks with Minimal Empirical Calibration. J Phys Chem B 2020; 124:1628-1636. [PMID: 32073849 PMCID: PMC7061328 DOI: 10.1021/acs.jpcb.9b10903] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2019] [Revised: 01/23/2020] [Indexed: 12/31/2022]
Abstract
Force field parametrization involves a complex set of linked optimization problems, with the goal of describing complex molecular interactions by using simple classical potential-energy functions that model Coulomb interactions, dispersion, and exchange repulsion. These functions comprise a set of atomic (and molecular) parameters and together with the bonded terms they constitute the molecular mechanics force field. Traditionally, many of these parameters have been fitted in a calibration approach in which experimental measures for thermodynamic and other relevant properties of small-molecule compounds are used for fitting and validation. As these approaches are laborious and time-consuming and represent an underdetermined optimization problem, we study methods to fit and derive an increasing number of parameters directly from electronic structure calculations, in order to greatly reduce possible parameter space for the remaining free parameters. In the current work we investigate a polarizable model with a higher order dispersion term for use in biomolecular simulation. Results for 49 biochemically relevant molecules are presented including updated parameters for hydrocarbon side chains. We show that our recently presented set of QM/MM derived atomic polarizabilities can be used in direct conjunction with partial charges and a higher order dispersion model that are quantum-mechanically determined, to freeze nearly all (i.e., 132 out of 138) nonbonded parameters to their quantum determined values.
Collapse
Affiliation(s)
- Koen M. Visscher
- AIMMS Division of Molecular
Toxicology, De Boelelaan
1108, 1081 HV Amsterdam, The Netherlands
| | - Daan P. Geerke
- AIMMS Division of Molecular
Toxicology, De Boelelaan
1108, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
3
|
Waibel C, Gross J. Polarizable Transferable Anisotropic United-Atom Force Field Based on the Mie Potential for Phase Equilibria: Ethers, n-Alkanes, and Nitrogen. J Chem Theory Comput 2019; 15:2561-2573. [DOI: 10.1021/acs.jctc.8b01238] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Christian Waibel
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Joachim Gross
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| |
Collapse
|
4
|
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: 167] [Impact Index Per Article: 33.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
|
5
|
Visscher KM, Geerke DP. Deriving Force-Field Parameters from First Principles Using a Polarizable and Higher Order Dispersion Model. J Chem Theory Comput 2019; 15:1875-1883. [PMID: 30763086 PMCID: PMC6581419 DOI: 10.1021/acs.jctc.8b01105] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2018] [Indexed: 11/30/2022]
Abstract
In this work we propose a strategy based on quantum mechanical (QM) calculations to parametrize a polarizable force field for use in molecular dynamics (MD) simulations. We investigate the use of multiple atoms-in-molecules (AIM) strategies to partition QM determined molecular electron densities into atomic subregions. The partitioned atomic densities are subsequently used to compute atomic dispersion coefficients from effective exchange-hole-dipole moment (XDM) calculations. In order to derive values for the repulsive van der Waals parameters from first principles, we use a simple volume relation to scale effective atomic radii. Explicit inclusion of higher order dispersion coefficients was tested for a series of alkanes, and we show that combining C6 and C8 attractive terms together with a C11 repulsive potential yields satisfying models when used in combination with our van der Waals parameters and electrostatic and bonded parameters as directly obtained from quantum calculations as well. This result highlights that explicit inclusion of higher order dispersion terms could be viable in simulation, and it suggests that currently available QM analysis methods allow for first-principles parametrization of molecular mechanics models.
Collapse
Affiliation(s)
- Koen M. Visscher
- AIMMS Division of Molecular Toxicology,
Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, the Netherlands
| | - Daan P. Geerke
- AIMMS Division of Molecular Toxicology,
Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, the Netherlands
| |
Collapse
|
6
|
Visscher KM, Swope WC, Geerke DP. A QM/MM Derived Polarizable Water Model for Molecular Simulation. Molecules 2018; 23:E3131. [PMID: 30501058 PMCID: PMC6321318 DOI: 10.3390/molecules23123131] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2018] [Revised: 11/19/2018] [Accepted: 11/21/2018] [Indexed: 11/26/2022] Open
Abstract
In this work, we propose an improved QM/MM-based strategy to determine condensed-phase polarizabilities and we use this approach to optimize a new and simple polarizable four-site water model for classical molecular simulation. For the determination of the model value for the polarizability from QM/MM, we show that our proposed consensus-fitting strategy significantly reduces the uncertainty in calculated polarizabilities in cases where the size of the local external electric field is small. By fitting electrostatic, polarization and dispersion properties of our water model based on quantum and/or combined QM/MM calculations, only a single model parameter (describing exchange repulsion) is left for empirical calibration. The resulting model performs well in describing relevant pure-liquid thermodynamic and transport properties, which illustrates the merit of our approach to minimize the number of free variables in our model.
Collapse
Affiliation(s)
- Koen M Visscher
- AIMMS Division of Molecular Toxicology, Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, De Boelelaan 1108, 1081 HV Amsterdam, The Netherlands.
| | - William C Swope
- IBM Almaden Research Center, 650 Harry Road, San Jose, CA 95120, USA.
| | - Daan P Geerke
- AIMMS Division of Molecular Toxicology, Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, De Boelelaan 1108, 1081 HV Amsterdam, The Netherlands.
| |
Collapse
|
7
|
Visscher KM, Vosmeer CR, Luirink RA, Geerke DP. A systematic approach to calibrate a transferable polarizable force field parameter set for primary alcohols. J Comput Chem 2018; 38:508-517. [PMID: 28133840 DOI: 10.1002/jcc.24702] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2016] [Revised: 10/24/2016] [Accepted: 11/14/2016] [Indexed: 11/05/2022]
Abstract
In this work, parameters are optimized for a charge-on-spring based polarizable force field for linear alcohols. We show that parameter transferability can be obtained using a systematic approach in which the effects of parameter changes on physico-chemical properties calculated from simulation are predicted. Our previously described QM/MM calculations are used to attribute condensed-phase polarizabilities, and starting from the non-polarizable GROMOS 53A5/53A6 parameter set, van der Waals and Coulomb interaction parameters are optimized to reproduce pure-liquid (thermodynamic, dielectric, and transport) properties, as well as hydration free energies. For a large set of models, which were obtained by combining small perturbations of 10 distinct parameters, values for pure-liquid properties of the series methanol to butanol were close to experiment. From this large set of models, we selected 34 models without special repulsive van der Waals parameters to distinguish between hydrogen-bonding and non-hydrogen-bonding atom pairs, to make the force field simple and transparent. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Koen M Visscher
- AIMMS Division of Molecular Toxicology, Department of Chemistry and Pharmaceutical Sciences, Faculty of Sciences, VU University Amsterdam, De Boelelaan 1108, HZ, Amsterdam, the Netherlands
| | - C Ruben Vosmeer
- AIMMS Division of Molecular Toxicology, Department of Chemistry and Pharmaceutical Sciences, Faculty of Sciences, VU University Amsterdam, De Boelelaan 1108, HZ, Amsterdam, the Netherlands
| | - Rosa A Luirink
- AIMMS Division of Molecular Toxicology, Department of Chemistry and Pharmaceutical Sciences, Faculty of Sciences, VU University Amsterdam, De Boelelaan 1108, HZ, Amsterdam, the Netherlands
| | - Daan P Geerke
- AIMMS Division of Molecular Toxicology, Department of Chemistry and Pharmaceutical Sciences, Faculty of Sciences, VU University Amsterdam, De Boelelaan 1108, HZ, Amsterdam, the Netherlands
| |
Collapse
|
8
|
Riquelme M, Lara A, Mobley DL, Verstraelen T, Matamala AR, Vöhringer-Martinez E. Hydration Free Energies in the FreeSolv Database Calculated with Polarized Iterative Hirshfeld Charges. J Chem Inf Model 2018; 58:1779-1797. [PMID: 30125107 PMCID: PMC6195221 DOI: 10.1021/acs.jcim.8b00180] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Computer simulations of biomolecular systems often use force fields, which are combinations of simple empirical atom-based functions to describe the molecular interactions. Even though polarizable force fields give a more detailed description of intermolecular interactions, nonpolarizable force fields, developed several decades ago, are often still preferred because of their reduced computation cost. Electrostatic interactions play a major role in biomolecular systems and are therein described by atomic point charges. In this work, we address the performance of different atomic charges to reproduce experimental hydration free energies in the FreeSolv database in combination with the GAFF force field. Atomic charges were calculated by two atoms-in-molecules approaches, Hirshfeld-I and Minimal Basis Iterative Stockholder (MBIS). To account for polarization effects, the charges were derived from the solute's electron density computed with an implicit solvent model, and the energy required to polarize the solute was added to the free energy cycle. The calculated hydration free energies were analyzed with an error model, revealing systematic errors associated with specific functional groups or chemical elements. The best agreement with the experimental data is observed for the AM1-BCC and the MBIS atomic charge methods. The latter includes the solvent polarization and presents a root-mean-square error of 2.0 kcal mol-1 for the 613 organic molecules studied. The largest deviation was observed for phosphorus-containing molecules and the molecules with amide, ester and amine functional groups.
Collapse
Affiliation(s)
- Maximiliano Riquelme
- Departamento de Físico-Química, Facultad de Ciencias Químicas , Universidad de Concepción , 4070386 Concepción , Chile
| | - Alejandro Lara
- Departamento de Físico-Química, Facultad de Ciencias Químicas , Universidad de Concepción , 4070386 Concepción , Chile
| | - David L Mobley
- Departments of Pharmaceutical Sciences and Chemistry, 147 Bison Modular , University of California, Irvine , Irvine , California 92617 , United States
| | - Toon Verstraelen
- Center for Molecular Modeling (CMM) , Ghent University , Technologiepark 903 , B-9052 Ghent , Belgium
| | - Adelio R Matamala
- Departamento de Físico-Química, Facultad de Ciencias Químicas , Universidad de Concepción , 4070386 Concepción , Chile
| | - Esteban Vöhringer-Martinez
- Departamento de Físico-Química, Facultad de Ciencias Químicas , Universidad de Concepción , 4070386 Concepción , Chile
| |
Collapse
|
9
|
Becker TM, Heinen J, Dubbeldam D, Lin LC, Vlugt TJH. Polarizable Force Fields for CO 2 and CH 4 Adsorption in M-MOF-74. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2017; 121:4659-4673. [PMID: 28286598 PMCID: PMC5338003 DOI: 10.1021/acs.jpcc.6b12052] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2016] [Revised: 01/25/2017] [Indexed: 05/28/2023]
Abstract
The family of M-MOF-74, with M = Co, Cr, Cu, Fe, Mg, Mn, Ni, Ti, V, and Zn, provides opportunities for numerous energy related gas separation applications. The pore structure of M-MOF-74 exhibits a high internal surface area and an exceptionally large adsorption capacity. The chemical environment of the adsorbate molecule in M-MOF-74 can be tuned by exchanging the metal ion incorporated in the structure. To optimize materials for a given separation process, insights into how the choice of the metal ion affects the interaction strength with adsorbate molecules and how to model these interactions are essential. Here, we quantitatively highlight the importance of polarization by comparing the proposed polarizable force field to orbital interaction energies from DFT calculations. Adsorption isotherms and heats of adsorption are computed for CO2, CH4, and their mixtures in M-MOF-74 with all 10 metal ions. The results are compared to experimental data, and to previous simulation results using nonpolarizable force fields derived from quantum mechanics. To the best of our knowledge, the developed polarizable force field is the only one so far trying to cover such a large set of possible metal ions. For the majority of metal ions, our simulations are in good agreement with experiments, demonstrating the effectiveness of our polarizable potential and the transferability of the adopted approach.
Collapse
Affiliation(s)
- Tim M Becker
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology , Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Jurn Heinen
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam , Science Park 904, 1098XH Amsterdam, The Netherlands
| | - David Dubbeldam
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands; Van't Hoff Institute for Molecular Sciences, University of Amsterdam, Science Park 904, 1098XH Amsterdam, The Netherlands
| | - Li-Chiang Lin
- William G. Lowrie Department of Chemical and Biomolecular Engineering, The Ohio State University , 151 West Woodruff Avenue, Columbus, Ohio 43210, United States
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology , Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| |
Collapse
|
10
|
Ploetz EA, Rustenburg AS, Geerke DP, Smith PE. To Polarize or Not to Polarize? Charge-on-Spring versus KBFF Models for Water and Methanol Bulk and Vapor–Liquid Interfacial Mixtures. J Chem Theory Comput 2016; 12:2373-87. [DOI: 10.1021/acs.jctc.5b01115] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Elizabeth A. Ploetz
- Department
of Chemistry, Kansas State University, Manhattan, Kansas 66506, United States
| | - Ariën S. Rustenburg
- AIMMS
Division of Molecular Toxicology, Department of Chemistry and Pharmaceutical
Sciences, Faculty of Sciences, VU University Amsterdam, 1081 HV Amsterdam, The Netherlands
| | - Daan P. Geerke
- AIMMS
Division of Molecular Toxicology, Department of Chemistry and Pharmaceutical
Sciences, Faculty of Sciences, VU University Amsterdam, 1081 HV Amsterdam, The Netherlands
| | - Paul E. Smith
- Department
of Chemistry, Kansas State University, Manhattan, Kansas 66506, United States
| |
Collapse
|
11
|
Kurnikov IV, Kurnikova M. Modeling Electronic Polarizability Changes in the Course of a Magnesium Ion Water Ligand Exchange Process. J Phys Chem B 2015; 119:10275-86. [PMID: 26109375 DOI: 10.1021/acs.jpcb.5b01295] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
This paper introduces explicit dependence of atomic polarizabilities on intermolecular interactions within the framework of a polarizable force field AMOEBA. Polarizable models used in biomolecular simulations often poorly describe molecular electrostatic induction in condensed phase, in part, due to neglect of a strong dependency of molecular electronic polarizability on intermolecular interactions at short distances. Our variable polarizability model parameters are derived from quantum chemical calculations of small clusters of atoms and molecules, and can be applied in simulations in condensed phase without additional scaling factors. The variable polarizability model is applied to simulate a ligand exchange reaction for a Mg(2+) ion solvated in water. Explicit dependence of water polarizability on a distance between a water oxygen and Mg(2+) is derived from in vacuum MP2 calculations of Mg(2+)-water dimer. The simulations yield a consistent description of the energetics of the Mg(2+)-water clusters of different size. Simulations also reproduce thermodynamics of ion solvation as well as kinetics of a water ligand exchange reaction. In contrast, simulations that used the additive force field or that used the constant polarizability models were not able to consistently and quantitatively describe the properties of the solvated Mg(2+) ion.
Collapse
Affiliation(s)
- Igor V Kurnikov
- Chemistry Department, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Maria Kurnikova
- Chemistry Department, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
12
|
Mei Y, Simmonett AC, Pickard FC, DiStasio RA, Brooks BR, Shao Y. Numerical study on the partitioning of the molecular polarizability into fluctuating charge and induced atomic dipole contributions. J Phys Chem A 2015; 119:5865-82. [PMID: 25945749 DOI: 10.1021/acs.jpca.5b03159] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In order to carry out a detailed analysis of the molecular static polarizability, which is the response of the molecule to a uniform external electric field, the molecular polarizability was computed using the finite-difference method for 21 small molecules, using density functional theory. Within nine charge population schemes (Löwdin, Mulliken, Becke, Hirshfeld, CM5, Hirshfeld-I, NPA, CHELPG, MK-ESP) in common use, the charge fluctuation contribution is found to dominate the molecular polarizability, with its ratio ranging from 59.9% with the Hirshfeld or CM5 scheme to 96.2% with the Mulliken scheme. The Hirshfeld-I scheme is also used to compute the other contribution to the molecular polarizability coming from the induced atomic dipoles, and the atomic polarizabilities in eight small molecules and water pentamer are found to be highly anisotropic for most atoms. Overall, the results suggest that (a) more emphasis probably should be placed on the charge fluctuation terms in future polarizable force field development and (b) an anisotropic polarizability might be more suitable than an isotropic one in polarizable force fields based entirely or partially on the induced atomic dipoles.
Collapse
Affiliation(s)
- Ye Mei
- †State Key Laboratory of Precision Spectroscopy, Department of Physics and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China.,‡NYU-ECNU Center for Computational Chemistry, NYU Shanghai, Shanghai 200062, China.,⊥Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Andrew C Simmonett
- ⊥Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Frank C Pickard
- ⊥Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Robert A DiStasio
- §Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Bernard R Brooks
- ⊥Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Yihan Shao
- ∥Q-Chem Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, United States
| |
Collapse
|
13
|
Vosmeer CR, Kiewisch K, Keijzer K, Visscher L, Geerke DP. A comparison between QM/MM and QM/QM based fitting of condensed-phase atomic polarizabilities. Phys Chem Chem Phys 2015; 16:17857-62. [PMID: 25042275 DOI: 10.1039/c4cp02401j] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Recently we reported a combined QM/MM approach to estimate condensed-phase values of atomic polarizabilities for use in (bio)molecular simulation. The setup relies on a MM treatment of the solvent when determining atomic polarizabilities to describe the response of a QM described solute to its external electric field. In this work, we study the effect of using alternative descriptions of the solvent molecules when evaluating atomic polarizabilities of a methanol solute. In a first step, we show that solute polarizabilities are not significantly affected upon substantially increasing the MM dipole moments towards values that are typically reported in literature for water solvent molecules. Subsequently, solute polarization is evaluated in the presence of a QM described solvent (using the frozen-density embedding method). In the latter case, lower oxygen polarizabilities were obtained than when using MM point charges to describe the solvent, due to introduction of Pauli-repulsion effects.
Collapse
Affiliation(s)
- C Ruben Vosmeer
- AIMMS Division of Molecular Toxicology, Department of Chemistry and Pharmaceutical Sciences, Faculty of Sciences, VU University Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, the Netherlands.
| | | | | | | | | |
Collapse
|
14
|
Graen T, Hoefling M, Grubmüller H. AMBER-DYES: Characterization of Charge Fluctuations and Force Field Parameterization of Fluorescent Dyes for Molecular Dynamics Simulations. J Chem Theory Comput 2014; 10:5505-12. [PMID: 26583233 DOI: 10.1021/ct500869p] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Recent advances in single molecule fluorescence experiments and theory allow a direct comparison and improved interpretation of experiment and simulation. To this end, force fields for a larger number of dyes are required which are compatible with and can be integrated into existing biomolecular force fields. Here, we developed, characterized, and implemented AMBER-DYES, a modular fluorescent label force field, for a set of 22 fluorescent dyes and their linkers from the Alexa, Atto, and Cy families, which are in common use for single molecule spectroscopy experiments. The force field is compatible with the AMBER protein force fields and the GROMACS molecular dynamics simulation program. The high electronic polarizability of the delocalized π-electron orbitals, as found in many fluorescent dyes, poses a particular challenge to point charge based force fields such as AMBER. To quantify the charge fluctuations due to the electronic polarizability, we simulated the 22 dyes in explicit solvent and sampled the charge fluctuations using QM/MM simulations at the B3LYP/6-31G*//TIP3P level of theory. The analysis of the simulations enabled us to derive ensemble fitted RESP charges from the solvated charge distributions of multiple trajectories. We observed broad, single peaked charge distributions for the conjugated ring atoms with well-defined mean values. The charge fitting procedure was validated against published charges of the dyelike amino acid tryptophan, which showed good agreement with existing tryptophan parameters from the AMBER, CHARMM, and OPLS force field families. A principal component analysis of the charge fluctuations revealed that a small number of collective coordinates suffices to describe most of the in-plane dye polarizability. The AMBER-DYES force field allows the rapid preparation of all atom molecular dynamics simulations of fluorescent systems for state of the art multi microsecond trajectories.
Collapse
Affiliation(s)
- Timo Graen
- Max Planck Institute for Biophysical Chemistry , Am Faβberg 11, 37077 Göttingen, Germany
| | - Martin Hoefling
- Max Planck Institute for Biophysical Chemistry , Am Faβberg 11, 37077 Göttingen, Germany
| | - Helmut Grubmüller
- Max Planck Institute for Biophysical Chemistry , Am Faβberg 11, 37077 Göttingen, Germany
| |
Collapse
|
15
|
Head-Gordon T, Rice J. Tribute to William C. Swope. J Phys Chem B 2014; 118:6357-9. [DOI: 10.1021/jp502812z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
16
|
Boulanger E, Thiel W. Solvent Boundary Potentials for Hybrid QM/MM Computations Using Classical Drude Oscillators: A Fully Polarizable Model. J Chem Theory Comput 2012; 8:4527-38. [DOI: 10.1021/ct300722e] [Citation(s) in RCA: 99] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Eliot Boulanger
- Max-Planck-Institut für Kohlenforschung,
Kaiser-Wilhelm-Platz
1, 45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung,
Kaiser-Wilhelm-Platz
1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|