1
|
Ramasami P, Murray JS. Radial Behavior of Electrostatic Potentials of Atoms and Ions Revisited: Isotropy and Anisotropy. Chemphyschem 2024; 25:e202400450. [PMID: 38775267 DOI: 10.1002/cphc.202400450] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2024] [Revised: 05/17/2024] [Indexed: 07/12/2024]
Abstract
In this paper we revisit earlier work relating to monoatomic atoms and ions published by pioneers in the area of electrostatic potentials. We include plots of the radial distributions of the electrostatic potentials for spherically symmetric atoms and cations, and for singly, doubly and triply negative anions. For atoms with anisotropy in their densities and electrostatic potentials, such as the halonium cations, it is shown how the molecular surface approach for plotting electrostatic potentials complements that achieved by directional radial distributions.
Collapse
Affiliation(s)
- Ponnadurai Ramasami
- Computational Chemistry Group, Department of Chemistry Faculty of Science, University of Mauritius, Reduit, 80837, Mauritius
- Centre of Natural Product Research, Department of Chemical Sciences University of Johannesburg, Doornfontein Campus, Johannesburg, 2028, South Africa
| | - Jane S Murray
- Department of Chemistry, University of New Orleans, New Orleans, LA, 70148, USA
| |
Collapse
|
2
|
Oliveira MP, Gonçalves YMH, Ol Gheta SK, Rieder SR, Horta BAC, Hünenberger PH. Comparison of the United- and All-Atom Representations of (Halo)alkanes Based on Two Condensed-Phase Force Fields Optimized against the Same Experimental Data Set. J Chem Theory Comput 2022; 18:6757-6778. [PMID: 36190354 DOI: 10.1021/acs.jctc.2c00524] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The level of accuracy that can be achieved by a force field is influenced by choices made in the interaction-function representation and in the relevant simulation parameters. These choices, referred to here as functional-form variants (FFVs), include for example the model resolution, the charge-derivation procedure, the van der Waals combination rules, the cutoff distance, and the treatment of the long-range interactions. Ideally, assessing the effect of a given FFV on the intrinsic accuracy of the force-field representation requires that only the specific FFV is changed and that this change is performed at an optimal level of parametrization, a requirement that may prove extremely challenging to achieve in practice. Here, we present a first attempt at such a comparison for one specific FFV, namely the choice of a united-atom (UA) versus an all-atom (AA) resolution in a force field for saturated acyclic (halo)alkanes. Two force-field versions (UA vs AA) are optimized in an automated way using the CombiFF approach against 961 experimental values for the pure-liquid densities ρliq and vaporization enthalpies ΔHvap of 591 compounds. For the AA force field, the torsional and third-neighbor Lennard-Jones parameters are also refined based on quantum-mechanical rotational-energy profiles. The comparison between the UA and AA resolutions is also extended to properties that have not been included as parameterization targets, namely the surface-tension coefficient γ, the isothermal compressibility κT, the isobaric thermal-expansion coefficient αP, the isobaric heat capacity cP, the static relative dielectric permittivity ϵ, the self-diffusion coefficient D, the shear viscosity η, the hydration free energy ΔGwat, and the free energy of solvation ΔGche in cyclohexane. For the target properties ρliq and ΔHvap, the UA and AA resolutions reach very similar levels of accuracy after optimization. For the nine other properties, the AA representation leads to more accurate results in terms of η; comparably accurate results in terms of γ, κT, αP, ϵ, D, and ΔGche; and less accurate results in terms of cP and ΔGwat. This work also represents a first step toward the calibration of a GROMOS-compatible force field at the AA resolution.
Collapse
Affiliation(s)
- Marina P Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Yan M H Gonçalves
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - S Kashef Ol Gheta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Salomé R Rieder
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Bruno A C Horta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Philippe H Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| |
Collapse
|
3
|
Rackers JA, Silva RR, Wang Z, Ponder JW. Polarizable Water Potential Derived from a Model Electron Density. J Chem Theory Comput 2021; 17:7056-7084. [PMID: 34699197 DOI: 10.1021/acs.jctc.1c00628] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A new empirical potential for efficient, large scale molecular dynamics simulation of water is presented. The HIPPO (Hydrogen-like Intermolecular Polarizable POtential) force field is based upon the model electron density of a hydrogen-like atom. This framework is used to derive and parametrize individual terms describing charge penetration damped permanent electrostatics, damped polarization, charge transfer, anisotropic Pauli repulsion, and damped dispersion interactions. Initial parameter values were fit to Symmetry Adapted Perturbation Theory (SAPT) energy components for ten water dimer configurations, as well as the radial and angular dependence of the canonical dimer. The SAPT-based parameters were then systematically refined to extend the treatment to water bulk phases. The final HIPPO water model provides a balanced representation of a wide variety of properties of gas phase clusters, liquid water, and ice polymorphs, across a range of temperatures and pressures. This water potential yields a rationalization of water structure, dynamics, and thermodynamics explicitly correlated with an ab initio energy decomposition, while providing a level of accuracy comparable or superior to previous polarizable atomic multipole force fields. The HIPPO water model serves as a cornerstone around which similarly detailed physics-based models can be developed for additional molecular species.
Collapse
Affiliation(s)
- Joshua A Rackers
- Program in Computational & Molecular Biophysics, Washington University School of Medicine, Saint Louis, Missouri 63110, United States.,Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87123, United States
| | - Roseane R Silva
- 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
| | - Jay W Ponder
- Department of Chemistry, Washington University in Saint Louis, Saint Louis, Missouri 63130, United States.,Department of Biochemistry & Molecular Biophysics, Washington University School of Medicine, Saint Louis, Missouri 63110, United States
| |
Collapse
|
4
|
Jiménez-Grávalos F, Suárez D. A Quantum Chemical Topology Picture of Intermolecular Electrostatic Interactions and Charge Penetration Energy. J Chem Theory Comput 2021; 17:4981-4995. [PMID: 34279923 PMCID: PMC8901103 DOI: 10.1021/acs.jctc.1c00263] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Based on the Interacting Quantum Atoms approach, we present herein a conceptual and theoretical framework of short-range electrostatic interactions, whose accurate description is still a challenging problem in molecular modeling. For all the noncovalent complexes in the S66 database, the fragment-based and atomic decomposition of the electrostatic binding energies is performed using both the charge density of the dimers and the unrelaxed densities of the monomers. This energy decomposition together with dispersion corrections gives rise to a pairwise approximation to the total binding energy. It also provides energetic descriptors at varying distance that directly address the atomic and molecular electrostatic interactions as described by point-charge or multipole-based potentials. Additionally, we propose a consistent definition of the charge penetration energy within quantum chemical topology, which is mainly characterized in terms of the intramolecular electrostatic energy. Finally, we discuss some practical implications of our results for the design and validation of electrostatic potentials.
Collapse
Affiliation(s)
| | - Dimas Suárez
- Departamento de Química Física y Analítica, Universidad de Oviedo, E-33006 Oviedo, Spain
| |
Collapse
|
5
|
Combining Molecular Dynamic Information and an Aspherical-Atom Data Bank in the Evaluation of the Electrostatic Interaction Energy in Multimeric Protein-Ligand Complex: A Case Study for HIV-1 Protease. Molecules 2021; 26:molecules26133872. [PMID: 34202892 PMCID: PMC8270314 DOI: 10.3390/molecules26133872] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2021] [Revised: 06/19/2021] [Accepted: 06/22/2021] [Indexed: 11/23/2022] Open
Abstract
Computational analysis of protein–ligand interactions is of crucial importance for drug discovery. Assessment of ligand binding energy allows us to have a glimpse of the potential of a small organic molecule to be a ligand to the binding site of a protein target. Available scoring functions, such as in docking programs, all rely on equations that sum each type of protein–ligand interactions in order to predict the binding affinity. Most of the scoring functions consider electrostatic interactions involving the protein and the ligand. Electrostatic interactions constitute one of the most important part of total interactions between macromolecules. Unlike dispersion forces, they are highly directional and therefore dominate the nature of molecular packing in crystals and in biological complexes and contribute significantly to differences in inhibition strength among related enzyme inhibitors. In this study, complexes of HIV-1 protease with inhibitor molecules (JE-2147 and darunavir) were analyzed by using charge densities from the transferable aspherical-atom University at Buffalo Databank (UBDB). Moreover, we analyzed the electrostatic interaction energy for an ensemble of structures, using molecular dynamic simulations to highlight the main features of electrostatic interactions important for binding affinity.
Collapse
|
6
|
Wang C, You ZX, Xing YH, Bai FY, Shi Z. A three-dimensional supramolecular network structure through hydrogen bonding and π–π interaction: synthesis, structure, and the fluorescence detection of balsalazide disodium. CrystEngComm 2021. [DOI: 10.1039/d1ce00373a] [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/31/2022]
Abstract
Three co-crystal compounds with different dimensions were synthesized. The fluorescence properties of co-crystal 1 were studied, especially detecting balsalazide disodium.
Collapse
Affiliation(s)
- Chen Wang
- College of Chemistry and Chemical Engineering
- Liaoning Normal University
- Dalian 116029
- P. R. China
| | - Zi Xin You
- College of Chemistry and Chemical Engineering
- Liaoning Normal University
- Dalian 116029
- P. R. China
| | - Yong-Heng Xing
- College of Chemistry and Chemical Engineering
- Liaoning Normal University
- Dalian 116029
- P. R. China
| | - Feng-Ying Bai
- College of Chemistry and Chemical Engineering
- Liaoning Normal University
- Dalian 116029
- P. R. China
| | - Zhan Shi
- State Key Laboratory of Inorganic Synthesis and Preparative Chemistry
- College of Chemistry
- Jilin University
- Changchun 130012
- P. R. China
| |
Collapse
|
7
|
Li F, Zheng Z, Xia S, Yu L. Synthesis, co-crystal structure, and DFT calculations of a multicomponent co-crystal constructed from 1H-benzotriazole and tetrafluoroterephthalic acid. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2020.128480] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
|
8
|
Peels M, Knizia G. Fast Evaluation of Two-Center Integrals over Gaussian Charge Distributions and Gaussian Orbitals with General Interaction Kernels. J Chem Theory Comput 2020; 16:2570-2583. [PMID: 32040326 DOI: 10.1021/acs.jctc.9b01296] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
We present efficient algorithms for computing two-center integrals and integral derivatives, with general interaction kernels K(r12), over Gaussian charge distributions of general angular momenta l. While formulated in terms of traditional ab initio integration techniques, full derivations and required secondary information, as well as a reference implementation, are provided to make the content accessible to other fields. Concretely, the presented algorithms are based on an adaption of the McMurchie-Davidson Recurrence Relation (MDRR) combined with analytical properties of the solid harmonic transformation; this obviates all intermediate recurrences except the adapted MDRR itself, and allows it to be applied to fully contracted auxiliary kernel integrals. The technique is particularly well-suited for semiempirical molecular orbital methods, where it can serve as a more general and efficient replacement of Slater-Koster tables, and for first-principles quantum chemistry methods employing density fitting. But the formalism's high efficiency and ability of handling general interaction kernels K(r12) and multipolar Gaussian charge distributions may also be of interest for modeling electrostatic interactions and short-range exchange and charge penetration effects in classical force fields and model potentials. With the presented technique, a 4894 × 4894 univ-JKFIT Coulomb matrix JAB = (A|1/r12|B) (183 MiB) can be computed in 50 ms on a Q2'2018 notebook CPU, without any screening or approximations.
Collapse
Affiliation(s)
- Mieke Peels
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Gerald Knizia
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| |
Collapse
|
9
|
Das AK, Urban L, Leven I, Loipersberger M, Aldossary A, Head-Gordon M, Head-Gordon T. Development of an Advanced Force Field for Water Using Variational Energy Decomposition Analysis. J Chem Theory Comput 2019; 15:5001-5013. [DOI: 10.1021/acs.jctc.9b00478] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|
10
|
Bojarowski SA, Kumar P, Wandtke CM, Dittrich B, Dominiak PM. Universal Method for Electrostatic Interaction Energies Estimation with Charge Penetration and Easily Attainable Point Charges. J Chem Theory Comput 2018; 14:6336-6345. [PMID: 30359528 DOI: 10.1021/acs.jctc.8b00781] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Our new model of electron density augmented by point charges (aug-PROmol) provides an estimation of electrostatic interaction energies including penetration effects ( ChemPhysChem 2016, 17, 2455-2460). In this paper we prove that it can be applied using sources of point charges other than those from direct restrained fitting to electrostatic potential (RESP). We used a newly established databank of tabulated invariom point charges and a widely known semiempirical method. Both sources perform equivalently to the basic aug-PROmol method as well as to reference energies at the DFT-SAPT/aug-cc-pVTZ level of theory. This is possible due to the universal character of the penetration model included in the aug-PROmol. Aug-PROmol may become a basis for development of new nonbonded terms in force fields or a high success rate scoring function.
Collapse
Affiliation(s)
- Sławomir A Bojarowski
- Biological and Chemical Research Centre, Department of Chemistry , University of Warsaw , ul. Żwirki i Wigury 101 , 02-089 Warszawa , Poland
| | - Prashant Kumar
- Biological and Chemical Research Centre, Department of Chemistry , University of Warsaw , ul. Żwirki i Wigury 101 , 02-089 Warszawa , Poland
| | - Claudia M Wandtke
- Institut fur Anorganische Chemie , Georg-August-Universitat , Tammannstrasse 4 , 37077 Göttingen , Germany
| | - Birger Dittrich
- Anorganische Chemie und Strukturchemie , Heinrich-Heine Universitat Düsseldorf , Universitatsstrasse 1, Gëbaude 26.42.01.21 , 40225 Düsseldorf , Germany
| | - Paulina M Dominiak
- Biological and Chemical Research Centre, Department of Chemistry , University of Warsaw , ul. Żwirki i Wigury 101 , 02-089 Warszawa , Poland
| |
Collapse
|
11
|
Ghahremanpour MM, van Maaren PJ, Caleman C, Hutchison GR, van der Spoel D. Polarizable Drude Model with s-Type Gaussian or Slater Charge Density for General Molecular Mechanics Force Fields. J Chem Theory Comput 2018; 14:5553-5566. [PMID: 30281307 DOI: 10.1021/acs.jctc.8b00430] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Gas-phase electric properties of molecules can be computed routinely using wave function methods or density functional theory (DFT). However, these methods remain computationally expensive for high-throughput screening of the vast chemical space of virtual compounds. Therefore, empirical force fields are a more practical choice in many cases, particularly since force field methods allow one to routinely predict the physicochemical properties in the condensed phases. This work presents Drude polarizable models, to increase the physical realism in empirical force fields, where the core particle is treated as a point charge and the Drude particle is treated either as a 1 s-Gaussian or a ns-Slater ( n = 1, 2, 3) charge density. Systematic parametrization to large high-quality quantum chemistry data obtained from the open access Alexandria Library ( https://doi.org/10.5281/zenodo.1004711 ) ensures the transferability of these parameters. The dipole moments and isotropic polarizabilities of the isolated molecules predicted by the proposed Drude models are in agreement with experiment with accuracy similar to DFT calculations at the B3LYP/aug-cc-pVTZ level of theory. The results show that the inclusion of explicit polarization into the models reduces the root-mean-square deviation with respect to DFT calculations of the predicted dipole moments of 152 dimers and clusters by more than 50%. Finally, we show that the accuracy of the electrostatic interaction energy of the water dimers can be improved systematically by the introduction of polarizable smeared charges as a model for charge penetration.
Collapse
Affiliation(s)
- Mohammad Mehdi Ghahremanpour
- Uppsala Center for Computational Chemistry, Department of Cell and Molecular Biology , Uppsala University , Husargatan 3 , Box 596, SE-75124 Uppsala , Sweden
| | - Paul J van Maaren
- Uppsala Center for Computational Chemistry, Department of Cell and Molecular Biology , Uppsala University , Husargatan 3 , Box 596, SE-75124 Uppsala , Sweden
| | - Carl Caleman
- Department of Physics and Astronomy , Uppsala University , Box 516, SE-75120 Uppsala , Sweden.,Center for Free-Electron Laser Science , Deutsches Elektronen-Synchrotron , DE-22607 Hamburg , Germany
| | - Geoffrey R Hutchison
- Department of Chemistry , University of Pittsburgh , Pittsburgh , Pennsylvania 15260 , United States
| | - David van der Spoel
- Uppsala Center for Computational Chemistry, Department of Cell and Molecular Biology , Uppsala University , Husargatan 3 , Box 596, SE-75124 Uppsala , Sweden
| |
Collapse
|
12
|
Walz MM, Ghahremanpour MM, van Maaren PJ, van der Spoel D. Phase-Transferable Force Field for Alkali Halides. J Chem Theory Comput 2018; 14:5933-5948. [DOI: 10.1021/acs.jctc.8b00507] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Marie-Madeleine Walz
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - Mohammad M. Ghahremanpour
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - Paul J. van Maaren
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - David van der Spoel
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| |
Collapse
|
13
|
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
|
14
|
A series of cocrystals formed by 2,3-dimethylpyrazine bridging various aromatic acids through hydrogen bonds: Synthesis, structural characterization and synthon discussion. J Mol Struct 2018. [DOI: 10.1016/j.molstruc.2018.03.105] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
|
15
|
Gryn’ova G, Corminboeuf C. Steric "attraction": not by dispersion alone. Beilstein J Org Chem 2018; 14:1482-1490. [PMID: 30013675 PMCID: PMC6037011 DOI: 10.3762/bjoc.14.125] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2018] [Accepted: 06/02/2018] [Indexed: 12/23/2022] Open
Abstract
Non-covalent interactions between neutral, sterically hindered organic molecules generally involve a strong stabilizing contribution from dispersion forces that in many systems turns the 'steric repulsion' into a 'steric attraction'. In addition to London dispersion, such systems benefit from electrostatic stabilization, which arises from a short-range effect of charge penetration and gets bigger with increasing steric bulk. In the present work, we quantify this contribution for a diverse set of molecular cores, ranging from unsubstituted benzene and cyclohexane to their derivatives carrying tert-butyl, phenyl, cyclohexyl and adamantyl substituents. While the importance of electrostatic interactions in the dimers of sp2-rich (e.g., π-conjugated) cores is well appreciated, less polarizable assemblies of sp3-rich systems with multiple short-range CH···HC contacts between the bulky cyclohexyl and adamantyl moieties are also significantly influenced by electrostatics. Charge penetration is drastically larger in absolute terms for the sp2-rich cores, but still has a non-negligible effect on the sp3-rich dimers, investigated herein, both in terms of their energetics and equilibrium interaction distances. These results emphasize the importance of this electrostatic effect, which has so far been less recognized in aliphatic systems compared to London dispersion, and are therefore likely to have implications for the development of force fields and methods for crystal structure prediction.
Collapse
Affiliation(s)
- Ganna Gryn’ova
- Institut des Sciences et Ingénierie Chimiques, École polytechnique fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Clémence Corminboeuf
- Institut des Sciences et Ingénierie Chimiques, École polytechnique fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| |
Collapse
|
16
|
Esser A, Belsare S, Marx D, Head-Gordon T. Mode specific THz spectra of solvated amino acids using the AMOEBA polarizable force field. Phys Chem Chem Phys 2018; 19:5579-5590. [PMID: 28165073 DOI: 10.1039/c6cp07388c] [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/21/2022]
Abstract
We have used the AMOEBA model to simulate the THz spectra of two zwitterionic amino acids in aqueous solution, which is compared to the results on these same systems using ab initio molecular dynamics (AIMD) simulations. Overall we find that the polarizable force field shows promising agreement with AIMD data for both glycine and valine in water. This includes the THz spectral assignments and the mode-specific spectral decomposition into intramolecular solute motions as well as distinct solute-water cross-correlation modes some of which cannot be captured by non-polarizable force fields that rely on fixed partial charges. This bodes well for future studies for simulating and decomposing the THz spectra for larger solutes such as proteins or polymers for which AIMD studies are presently intractable. Furthermore, we believe that the current study on rather simple aqueous solutions offers a way to systematically investigate the importance of charge transfer, nuclear quantum effects, and the validity of computationally practical density functionals, all of which are needed to fully quantitatively capture complex dynamical motions in the condensed phase.
Collapse
Affiliation(s)
- Alexander Esser
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany.
| | - Saurabh Belsare
- The UC Berkeley - UCSF Graduate Program in Bioengineering, University of California, Berkeley, CA, 94720, USA
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany.
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, Department of Bioengineering, Department of Chemical and Biomolecular Engineering, University of California Berkeley, CA 94720, USA.
| |
Collapse
|
17
|
Rackers JA, Wang Q, Liu C, Piquemal JP, Ren P, Ponder JW. An optimized charge penetration model for use with the AMOEBA force field. Phys Chem Chem Phys 2018; 19:276-291. [PMID: 27901142 DOI: 10.1039/c6cp06017j] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The principal challenge of using classical physics to model biomolecular interactions is capturing the nature of short-range interactions that drive biological processes from nucleic acid base stacking to protein-ligand binding. In particular most classical force fields suffer from an error in their electrostatic models that arises from an ability to account for the overlap between charge distributions occurring when molecules get close to each other, known as charge penetration. In this work we present a simple, physically motivated model for including charge penetration in the AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications) force field. With a function derived from the charge distribution of a hydrogen-like atom and a limited number of parameters, our charge penetration model dramatically improves the description of electrostatics at short range. On a database of 101 biomolecular dimers, the charge penetration model brings the error in the electrostatic interaction energy relative to the ab initio SAPT electrostatic interaction energy from 13.4 kcal mol-1 to 1.3 kcal mol-1. The model is shown not only to be robust and transferable for the AMOEBA model, but also physically meaningful as it universally improves the description of the electrostatic potential around a given molecule.
Collapse
Affiliation(s)
- Joshua A Rackers
- Program in Computational & Molecular Biophysics, Washington University, School of Medicine, Saint Louis, Missouri 63110, USA
| | - Qiantao Wang
- 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
| | - 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, USA
| | - Jay W Ponder
- Department of Chemistry, Washington University in Saint Louis, Saint Louis, Missouri 63130, USA.
| |
Collapse
|
18
|
Bojarowski SA, Kumar P, Dominiak PM. Interplay of point multipole moments and charge penetration for intermolecular electrostatic interaction energies from the University at Buffalo pseudoatom databank model of electron density. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2017; 73:598-609. [PMID: 28762970 DOI: 10.1107/s2052520617005510] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/07/2017] [Accepted: 04/12/2017] [Indexed: 06/07/2023]
Abstract
The strength of the University at Buffalo DataBank (UBDB) in Ees estimation is mainly due to charge overlap effects because the UBDB offers continuous representation of charge density which allows for a direct account of charge penetration in the derivation of electrostatic energies. In the UBDB model, these effects begin to play an important role at distances below twice the equilibrium distance and significantly increase as distances decrease. At equilibrium distances they are responsible for 30-50% of Ees for polar molecules and around 90% of Ees for nonpolar molecules. When the energy estimation from the UBDB is reduced to point multipoles, the results are comparable to point charges fitted to electrostatic potentials. On the other hand, particular components of energy from point multipole moments from the UBDB model are sensitive to the type of interaction and might be helpful in the characterization of interactions.
Collapse
Affiliation(s)
- Sławomir A Bojarowski
- Biological and Chemical Research Center, Department of Chemistry, University of Warsaw, ul. Zwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Prashant Kumar
- Biological and Chemical Research Center, Department of Chemistry, University of Warsaw, ul. Zwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Paulina M Dominiak
- Biological and Chemical Research Center, Department of Chemistry, University of Warsaw, ul. Zwirki i Wigury 101, 02-089 Warszawa, Poland
| |
Collapse
|
19
|
Mao Y, Shao Y, Dziedzic J, Skylaris CK, Head-Gordon T, Head-Gordon M. Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis. J Chem Theory Comput 2017; 13:1963-1979. [PMID: 28430427 DOI: 10.1021/acs.jctc.7b00089] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The importance of incorporating solvent polarization effects into the modeling of solvation processes has been well-recognized, and therefore a new generation of hybrid quantum mechanics/molecular mechanics (QM/MM) approaches that accounts for this effect is desirable. We present a fully self-consistent, mutually polarizable QM/MM scheme using the AMOEBA force field, in which the total energy of the system is variationally minimized with respect to both the QM electronic density and the MM induced dipoles. This QM/AMOEBA model is implemented through the Q-Chem/LibEFP code interface and then applied to the evaluation of solute-solvent interaction energies for various systems ranging from the water dimer to neutral and ionic solutes (NH3, NH4+, CN-) surrounded by increasing numbers of water molecules (up to 100). In order to analyze the resulting interaction energies, we also utilize an energy decomposition analysis (EDA) scheme which identifies contributions from permanent electrostatics, polarization, and van der Waals (vdW) interaction for the interaction between the QM solute and the solvent molecules described by AMOEBA. This facilitates a component-wise comparison against full QM calculations where the corresponding energy components are obtained via a modified version of the absolutely localized molecular orbitals (ALMO)-EDA. The results show that the present QM/AMOEBA model can yield reasonable solute-solvent interaction energies for neutral and cationic species, while further scrutiny reveals that this accuracy highly relies on the delicate balance between insufficiently favorable permanent electrostatics and softened vdW interaction. For anionic solutes where the charge penetration effect becomes more pronounced, the QM/MM interface turns out to be unbalanced. These results are consistent with and further elucidate our findings in a previous study using a slightly different QM/AMOEBA model ( Dziedzic et al. J. Chem. Phys. 2016 , 145 , 124106 ). The implications of these results for further refinement of this model are also discussed.
Collapse
Affiliation(s)
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma , Norman, Oklahoma 73019, United States
| | - Jacek Dziedzic
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, U.K.,Faculty of Applied Physics and Mathematics, Gdańsk University of Technology , Gdańsk 80-233, Poland
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, U.K
| | | | - Martin Head-Gordon
- Chemical Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
20
|
Thornton AW, Simon CM, Kim J, Kwon O, Deeg KS, Konstas K, Pas SJ, Hill MR, Winkler DA, Haranczyk M, Smit B. Materials Genome in Action: Identifying the Performance Limits of Physical Hydrogen Storage. CHEMISTRY OF MATERIALS : A PUBLICATION OF THE AMERICAN CHEMICAL SOCIETY 2017; 29:2844-2854. [PMID: 28413259 PMCID: PMC5390509 DOI: 10.1021/acs.chemmater.6b04933] [Citation(s) in RCA: 77] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/21/2016] [Revised: 03/06/2017] [Indexed: 05/29/2023]
Abstract
The Materials Genome is in action: the molecular codes for millions of materials have been sequenced, predictive models have been developed, and now the challenge of hydrogen storage is targeted. Renewably generated hydrogen is an attractive transportation fuel with zero carbon emissions, but its storage remains a significant challenge. Nanoporous adsorbents have shown promising physical adsorption of hydrogen approaching targeted capacities, but the scope of studies has remained limited. Here the Nanoporous Materials Genome, containing over 850 000 materials, is analyzed with a variety of computational tools to explore the limits of hydrogen storage. Optimal features that maximize net capacity at room temperature include pore sizes of around 6 Å and void fractions of 0.1, while at cryogenic temperatures pore sizes of 10 Å and void fractions of 0.5 are optimal. Our top candidates are found to be commercially attractive as "cryo-adsorbents", with promising storage capacities at 77 K and 100 bar with 30% enhancement to 40 g/L, a promising alternative to liquefaction at 20 K and compression at 700 bar.
Collapse
Affiliation(s)
- Aaron W. Thornton
- Future Industries, Commonwealth Scientific and Industrial Research Organisation, Private Bag 10, Clayton Soutth MDC, Victoria 3169, Australia
| | - Cory M. Simon
- Department of Chemical and Biomolecular Engineering and Department of Chemistry, University of California, Berkeley, California 94720-1462, United States
| | - Jihan Kim
- Department
of Chemical and Biomolecular Engineering, Korea Advanced Institute of Science and Technology, 291 Daehak-ro Yuseong-gu, Daejeon, 305-701, Korea
| | - Ohmin Kwon
- Department
of Chemical and Biomolecular Engineering, Korea Advanced Institute of Science and Technology, 291 Daehak-ro Yuseong-gu, Daejeon, 305-701, Korea
| | - Kathryn S. Deeg
- Department of Chemical and Biomolecular Engineering and Department of Chemistry, University of California, Berkeley, California 94720-1462, United States
| | - Kristina Konstas
- Future Industries, Commonwealth Scientific and Industrial Research Organisation, Private Bag 10, Clayton Soutth MDC, Victoria 3169, Australia
| | - Steven J. Pas
- Power & Energy Systems, Maritime Division, Defence Science and
Technology Group, Department of Defence, 506 Lorimer Street, Fishermans Bend, Victoria 3207, Australia
- School of Chemistry and Department of Chemical
Engineering, Monash University, Clayton, Victoria 3800, Australia
| | - Matthew R. Hill
- Future Industries, Commonwealth Scientific and Industrial Research Organisation, Private Bag 10, Clayton Soutth MDC, Victoria 3169, Australia
- School of Chemistry and Department of Chemical
Engineering, Monash University, Clayton, Victoria 3800, Australia
| | - David A. Winkler
- Future Industries, Commonwealth Scientific and Industrial Research Organisation, Private Bag 10, Clayton Soutth MDC, Victoria 3169, Australia
- Monash
Institute of Pharmaceutical Sciences, 381 Royal Parade, Parkville, Victoria 3052, Australia
- Latrobe Institute for Molecular Science, Bundoora, Victoria 3046, Australia
- School of Chemical
and Physical Sciences, Flinders University, Bedford Park, South Australia 5042, Australia
| | - Maciej Haranczyk
- Computational
Research Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720-8139, United States
| | - Berend Smit
- Department of Chemical and Biomolecular Engineering and Department of Chemistry, University of California, Berkeley, California 94720-1462, United States
- Laboratory of Molecular Simulation, Institut des Sciences et Ingénierie Chimiques, Valais, Rue de l’Industrie
17, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1950 Sion, Switzerland
| |
Collapse
|
21
|
Verma P, Wang B, Fernandez LE, Truhlar DG. Physical Molecular Mechanics Method for Damped Dispersion. J Phys Chem A 2017; 121:2855-2862. [PMID: 28328203 DOI: 10.1021/acs.jpca.7b02384] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Damped dispersion can be a significant component of the interaction energy in many physical and chemical processes, for example, physisorption and noncovalent complexation. For physically interpreting and modeling such processes, it is convenient to have an analytic method to calculate damped dispersion that is readily applicable across the entire periodic table. Of the available methods to calculate damped dispersion energy for interacting systems with overlapping charge distributions, we select symmetry-adapted perturbation theory (SAPT) as providing a reasonable definition, and of the possible analytic forms, we choose the D3(BJ) method. However, the available parametrizations of D3(BJ) include not only damped dispersion energy but also corrections for errors in specific exchange-correlation functionals. Here we present a parametrization that provides a physical measure of damped dispersion without such density functional corrections. The method generalizes an earlier method of Pernal and co-workers to all elements from hydrogen to plutonium.
Collapse
Affiliation(s)
- Pragya Verma
- Department of Chemistry, Nanoporous Materials Genome Center, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota , 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455-0431, United States
| | - Bo Wang
- Department of Chemistry, Nanoporous Materials Genome Center, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota , 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455-0431, United States
| | - Laura E Fernandez
- Department of Chemistry, Nanoporous Materials Genome Center, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota , 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Nanoporous Materials Genome Center, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota , 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
22
|
Diller DJ. Development of a Simple Electron Transfer and Polarization Model and Its Application to Biological Systems. J Chem Theory Comput 2017; 13:329-339. [PMID: 27958737 DOI: 10.1021/acs.jctc.6b00852] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Here we present a new method for point charge calculation which we call QET (charges by electron transfer). The intent of this work is to develop a method that can be useful for studying charge transfer in large biological systems. It is based on the intuitive framework of the QEQ method with the key difference being that the QET method tracks all pairwise electron transfers by augmenting the QEQ pseudoenergy function with a distance dependent cost function for each electron transfer. This approach solves the key limitation of the QEQ method which is its handling of formally charged groups. First, we parametrize the QET method by fitting to electrostatic potentials calculated using ab initio quantum mechanics on over 11,000 small molecules. On an external test set of over 2500 small molecules the QET method achieves a mean absolute error of 1.37 kcal/mol/electron when compared to the ab initio electrostatic potentials. Second, we examine the conformational dependence of the charges on over 2700 tripeptides. With the tripeptide data set, we show that the conformational effects account for approximately 0.4 kcal/mol/electron on the electrostatic potentials. Third, we test the QET method for its ability to reproduce the effects of polarization and electron transfer on 1000 water clusters. For the water clusters, we show that the QET method captures about 50% of the polarization and electron transfer effects. Finally, we examine the effects of electron transfer and polarizability on the electrostatic interaction between p38 and 94 small molecule ligands. When used in conjunction with the Generalized-Born continuum solvent model, polarization and electron transfer with the QET model lead to an average change of 17 kcal/mol on the calculated electrostatic component of ΔG.
Collapse
Affiliation(s)
- David J Diller
- Data2Discovery Consulting, East Windsor, New Jersey 08520, United States
| |
Collapse
|
23
|
Liu W, Jia X, Wang M, Li P, Wang X, Hu W, Zheng J, Mei Y. Calculations of the absolute binding free energies for Ralstonia solanacearum lectins bound with methyl-α-l-fucoside at molecular mechanical and quantum mechanical/molecular mechanical levels. RSC Adv 2017. [DOI: 10.1039/c7ra06215j] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
In this work, both a molecular mechanical (MM) method and a hybrid quantum mechanical/molecular mechanical (QM/MM) method have been applied in the study of the binding affinities of methyl-α-l-fucoside to Ralstonia solanacearum lectins.
Collapse
Affiliation(s)
- Wei Liu
- State Key Laboratory of Precision Spectroscopy
- School of Physics and Materials Science
- East China Normal University
- Shanghai 200062
- China
| | - Xiangyu Jia
- State Key Laboratory of Precision Spectroscopy
- School of Physics and Materials Science
- East China Normal University
- Shanghai 200062
- China
| | - Meiting Wang
- State Key Laboratory of Precision Spectroscopy
- School of Physics and Materials Science
- East China Normal University
- Shanghai 200062
- China
| | - Pengfei Li
- State Key Laboratory of Precision Spectroscopy
- School of Physics and Materials Science
- East China Normal University
- Shanghai 200062
- China
| | - Xiaohui Wang
- State Key Laboratory of Precision Spectroscopy
- School of Physics and Materials Science
- East China Normal University
- Shanghai 200062
- China
| | - Wenxin Hu
- The Computer Center
- School of Computer Science and Software Engineering
- East China Normal University
- Shanghai 200062
- China
| | - Jun Zheng
- The Computer Center
- School of Computer Science and Software Engineering
- East China Normal University
- Shanghai 200062
- China
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy
- School of Physics and Materials Science
- East China Normal University
- Shanghai 200062
- China
| |
Collapse
|
24
|
Mao Y, Demerdash O, Head-Gordon M, Head-Gordon T. Assessing Ion-Water Interactions in the AMOEBA Force Field Using Energy Decomposition Analysis of Electronic Structure Calculations. J Chem Theory Comput 2016; 12:5422-5437. [PMID: 27709939 DOI: 10.1021/acs.jctc.6b00764] [Citation(s) in RCA: 65] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
AMOEBA is a molecular mechanics force field that addresses some of the shortcomings of a fixed partial charge model, by including permanent atomic point multipoles through quadrupoles, as well as many-body polarization through the use of point inducible dipoles. In this work, we investigate how well AMOEBA formulates its non-bonded interactions, and how it implicitly incorporates quantum mechanical effects such as charge penetration (CP) and charge transfer (CT), for water-water and water-ion interactions. We find that AMOEBA's total interaction energies, as a function of distance and over angular scans for the water dimer and for a range of water-monovalent cations, agree well with an advanced density functional theory (DFT) model, whereas the water-halides and water-divalent cations show significant disagreement with the DFT result, especially in the compressed region when the two fragments overlap. We use a second-generation energy decomposition analysis (EDA) scheme based on absolutely localized molecular orbitals (ALMOs) to show that in the best cases AMOEBA relies on cancellation of errors by softening of the van der Waals (vdW) wall to balance permanent electrostatics that are too unfavorable, thereby compensating for the missing CP effect. CT, as another important stabilizing effect not explicitly taken into account in AMOEBA, is also found to be incorporated by the softened vdW interaction. For the water-halides and water-divalent cations, this compensatory approach is not as well executed by AMOEBA over all distances and angles, wherein permanent electrostatics remains too unfavorable and polarization is overdamped in the former while overestimated in the latter. We conclude that the DFT-based EDA approach can help refine a next-generation AMOEBA model that either realizes a better cancellation of errors for problematic cases like those illustrated here, or serves to guide the parametrization of explicit functional forms for short-range contributions from CP and/or CT.
Collapse
Affiliation(s)
- Yuezhi Mao
- Kenneth S. Pitzer Center for Theoretical Chemistry, and Department of Chemistry, University of California at Berkeley , Berkeley, California 94720, United States
| | - Omar Demerdash
- Kenneth S. Pitzer Center for Theoretical Chemistry, and Department of Chemistry, University of California at Berkeley , Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, and Department of Chemistry, University of California at Berkeley , Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, and Department of Chemistry, University of California at Berkeley , Berkeley, California 94720, United States.,Department of Bioengineering, and Department of Chemical and Biomolecular Engineering, University of California at Berkeley , Berkeley, California 94720, United States
| |
Collapse
|
25
|
Verstraelen T, Vandenbrande S, Heidar-Zadeh F, Vanduyfhuys L, Van Speybroeck V, Waroquier M, Ayers PW. Minimal Basis Iterative Stockholder: Atoms in Molecules for Force-Field Development. J Chem Theory Comput 2016; 12:3894-912. [DOI: 10.1021/acs.jctc.6b00456] [Citation(s) in RCA: 83] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Toon Verstraelen
- Center
for Molecular Modeling (CMM), Member of the QCMM Ghent−Brussels
Alliance, Ghent University, Technologiepark 903, B9000 Ghent, Belgium
| | - Steven Vandenbrande
- Center
for Molecular Modeling (CMM), Member of the QCMM Ghent−Brussels
Alliance, Ghent University, Technologiepark 903, B9000 Ghent, Belgium
| | - Farnaz Heidar-Zadeh
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 West
Main Street, Hamilton, Ontario L8S 4M1, Canada
| | - Louis Vanduyfhuys
- Center
for Molecular Modeling (CMM), Member of the QCMM Ghent−Brussels
Alliance, Ghent University, Technologiepark 903, B9000 Ghent, Belgium
| | - Veronique Van Speybroeck
- Center
for Molecular Modeling (CMM), Member of the QCMM Ghent−Brussels
Alliance, Ghent University, Technologiepark 903, B9000 Ghent, Belgium
| | - Michel Waroquier
- Center
for Molecular Modeling (CMM), Member of the QCMM Ghent−Brussels
Alliance, Ghent University, Technologiepark 903, B9000 Ghent, Belgium
| | - Paul W. Ayers
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 West
Main Street, Hamilton, Ontario L8S 4M1, Canada
| |
Collapse
|
26
|
Bojarowski SA, Kumar P, Dominiak PM. A Universal and Straightforward Approach to Include Penetration Effects in Electrostatic Interaction Energy Estimation. Chemphyschem 2016; 17:2455-60. [PMID: 27166026 DOI: 10.1002/cphc.201600390] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2016] [Indexed: 12/23/2022]
Abstract
To compensate for the lack of the explicit treatment of charge penetration in classical force fields, we propose a new charge-distribution model based on a promolecule augmented with point charges (aug-PROmol). It relies on a superposition of spherical atomic electron densities obtained for each chemical element from SCF energy optimized atomic orbitals. Atomic densities are further rescaled by partial point charges computed from fits to the molecular electrostatic potential. Aug-PROmol was tested on the S66 benchmark dataset extended to nonequilibrium geometries (J. Chem. Theory Comput., 2011, 7, 3466). The model does not need any additional parametrization other than point charges. Despite its simplicity, aug-PROmol approximates the electrostatic energy with good agreement (RMSE=0.76 kcal mol(-1) to DFT-SAPT with B3LYP/aug-cc-pVTZ).
Collapse
Affiliation(s)
- Sławomir A Bojarowski
- Department of Chemistry, Biological and Chemical Research Centre, University of Warsaw, ul. Żwirki i Wigury 101, 02-089, Warszawa, Poland
| | - Prashant Kumar
- Department of Chemistry, Biological and Chemical Research Centre, University of Warsaw, ul. Żwirki i Wigury 101, 02-089, Warszawa, Poland
| | - Paulina M Dominiak
- Department of Chemistry, Biological and Chemical Research Centre, University of Warsaw, ul. Żwirki i Wigury 101, 02-089, Warszawa, Poland.
| |
Collapse
|
27
|
Öhrn A, Hermida-Ramon JM, Karlström G. Method for Slater-Type Density Fitting for Intermolecular Electrostatic Interactions with Charge Overlap. I. The Model. J Chem Theory Comput 2016; 12:2298-311. [DOI: 10.1021/acs.jctc.5b01155] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Anders Öhrn
- Department
of Theoretical Chemistry, Chemical Centre, P.O. Box 124, S-221 00 Lund, Sweden
| | - Jose M. Hermida-Ramon
- Departamento
de Química Física, Facultade de Química, Universidade de Vigo, 36310 Vigo, Spain
| | - Gunnar Karlström
- Department
of Theoretical Chemistry, Chemical Centre, P.O. Box 124, S-221 00 Lund, Sweden
| |
Collapse
|
28
|
Narth C, Lagardère L, Polack É, Gresh N, Wang Q, Bell DR, Rackers JA, Ponder JW, Ren PY, Piquemal JP. Scalable improvement of SPME multipolar electrostatics in anisotropic polarizable molecular mechanics using a general short-range penetration correction up to quadrupoles. J Comput Chem 2016; 37:494-506. [PMID: 26814845 PMCID: PMC4730919 DOI: 10.1002/jcc.24257] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2015] [Revised: 08/02/2015] [Accepted: 10/21/2015] [Indexed: 12/25/2022]
Abstract
We propose a general coupling of the Smooth Particle Mesh Ewald SPME approach for distributed multipoles to a short-range charge penetration correction modifying the charge-charge, charge-dipole and charge-quadrupole energies. Such an approach significantly improves electrostatics when compared to ab initio values and has been calibrated on Symmetry-Adapted Perturbation Theory reference data. Various neutral molecular dimers have been tested and results on the complexes of mono- and divalent cations with a water ligand are also provided. Transferability of the correction is adressed in the context of the implementation of the AMOEBA and SIBFA polarizable force fields in the TINKER-HP software. As the choices of the multipolar distribution are discussed, conclusions are drawn for the future penetration-corrected polarizable force fields highlighting the mandatory need of non-spurious procedures for the obtention of well balanced and physically meaningful distributed moments. Finally, scalability and parallelism of the short-range corrected SPME approach are addressed, demonstrating that the damping function is computationally affordable and accurate for molecular dynamics simulations of complex bio- or bioinorganic systems in periodic boundary conditions.
Collapse
Affiliation(s)
- Christophe Narth
- UPMC Univ. Paris 06, UMR 7616, Laboratoire de Chimie Théorique, F-75005, Paris, France
| | - Louis Lagardère
- UPMC Univ. Paris 06, Institut du Calcul et de la Simulation, F-75005, Paris, France
| | - Étienne Polack
- UPMC Univ. Paris 06, UMR 7616, Laboratoire de Chimie Théorique, F-75005, Paris, France
- UPMC Univ. Paris 06, UMR 7598, Laboratoire Jacques-Louis Lions, F-75005, Paris, France
| | - Nohad Gresh
- UPMC Univ. Paris 06, UMR 7616, Laboratoire de Chimie Théorique, F-75005, Paris, France
- Chemistry and Biology Nucleo(s)tides and immunology for Therapy (CBNIT), UMR 8601 CNRS, UFR Biomédicale, Paris 75006, France
| | - Qiantao Wang
- Department of Biomedical Engineering, The University of Texas at Austin, Texas 78712
| | - David R. Bell
- Department of Biomedical Engineering, The University of Texas at Austin, Texas 78712
| | - Joshua A. Rackers
- Department of Biochemistry and Molecular Biophysics, Washington University, St. Louis, Missouri 63110
| | - Jay W. Ponder
- Department of Biochemistry and Molecular Biophysics, Washington University, St. Louis, Missouri 63110
| | - Pengyu Y. Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Texas 78712
| | - Jean-Philip Piquemal
- UPMC Univ. Paris 06, UMR 7616, Laboratoire de Chimie Théorique, F-75005, Paris, France
| |
Collapse
|
29
|
Tafipolsky M, Ansorg K. Toward a Physically Motivated Force Field: Hydrogen Bond Directionality from a Symmetry-Adapted Perturbation Theory Perspective. J Chem Theory Comput 2016; 12:1267-79. [DOI: 10.1021/acs.jctc.5b01057] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Maxim Tafipolsky
- Institut
für Physikalische
und Theoretische Chemie, Universität Würzburg, Campus
Hubland Nord, Emil-Fischer-Strasse 42, D-97074 Würzburg, Germany
| | - Kay Ansorg
- Institut
für Physikalische
und Theoretische Chemie, Universität Würzburg, Campus
Hubland Nord, Emil-Fischer-Strasse 42, D-97074 Würzburg, Germany
| |
Collapse
|
30
|
Wallace AJ, Jayasinghe CD, Polson MIJ, Curnow OJ, Crittenden DL. Cyclopropenium Cations Break the Rules of Attraction to Form Closely Bound Dimers. J Am Chem Soc 2015; 137:15528-32. [DOI: 10.1021/jacs.5b10388] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Andrew J. Wallace
- Department of Chemistry, University of Canterbury, Private Bag 4800, Christchurch 8041, New Zealand
| | - Chaminda D. Jayasinghe
- Department of Chemistry, University of Canterbury, Private Bag 4800, Christchurch 8041, New Zealand
| | - Matthew I. J. Polson
- Department of Chemistry, University of Canterbury, Private Bag 4800, Christchurch 8041, New Zealand
| | - Owen J. Curnow
- Department of Chemistry, University of Canterbury, Private Bag 4800, Christchurch 8041, New Zealand
| | - Deborah L. Crittenden
- Department of Chemistry, University of Canterbury, Private Bag 4800, Christchurch 8041, New Zealand
| |
Collapse
|
31
|
Energetics of proton release on the first oxidation step in the water-oxidizing enzyme. Nat Commun 2015; 6:8488. [PMID: 26442814 PMCID: PMC4617610 DOI: 10.1038/ncomms9488] [Citation(s) in RCA: 101] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2015] [Accepted: 08/26/2015] [Indexed: 12/31/2022] Open
Abstract
In photosystem II (PSII), the Mn4CaO5 cluster catalyses the water splitting reaction. The crystal structure of PSII shows the presence of a hydrogen-bonded water molecule directly linked to O4. Here we show the detailed properties of the H-bonds associated with the Mn4CaO5 cluster using a quantum mechanical/molecular mechanical approach. When O4 is taken as a μ-hydroxo bridge acting as a hydrogen-bond donor to water539 (W539), the S0 redox state best describes the unusually short O4–OW539 distance (2.5 Å) seen in the crystal structure. We find that in S1, O4 easily releases the proton into a chain of eight strongly hydrogen-bonded water molecules. The corresponding hydrogen-bond network is absent for O5 in S1. The present study suggests that the O4-water chain could facilitate the initial deprotonation event in PSII. This unexpected insight is likely to be of real relevance to mechanistic models for water oxidation. The availability of crystal structures of photosystem II opens up the possibility of gaining insights into its mechanism. Here, the authors use a computational approach and propose a deprotonation event at O4 followed by long-range proton-transfer along a chain of strongly bonded water molecules.
Collapse
|
32
|
Wang Q, Rackers JA, He C, Qi R, Narth C, Lagardere L, Gresh N, Ponder JW, Piquemal JP, Ren P. General Model for Treating Short-Range Electrostatic Penetration in a Molecular Mechanics Force Field. J Chem Theory Comput 2015; 11:2609-2618. [PMID: 26413036 PMCID: PMC4570253 DOI: 10.1021/acs.jctc.5b00267] [Citation(s) in RCA: 83] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2015] [Indexed: 11/30/2022]
Abstract
Classical molecular mechanics force fields typically model interatomic electrostatic interactions with point charges or multipole expansions, which can fail for atoms in close contact due to the lack of a description of penetration effects between their electron clouds. These short-range penetration effects can be significant and are essential for accurate modeling of intermolecular interactions. In this work we report parametrization of an empirical charge-charge function previously reported (Piquemal J.-P.; J. Phys. Chem. A2003, 107, 10353) to correct for the missing penetration term in standard molecular mechanics force fields. For this purpose, we have developed a database (S101×7) of 101 unique molecular dimers, each at 7 different intermolecular distances. Electrostatic, induction/polarization, repulsion, and dispersion energies, as well as the total interaction energy for each complex in the database are calculated using the SAPT2+ method (Parker T. M.; J. Chem. Phys.2014, 140, 094106). This empirical penetration model significantly improves agreement between point multipole and quantum mechanical electrostatic energies across the set of dimers and distances, while using only a limited set of parameters for each chemical element. Given the simplicity and effectiveness of the model, we expect the electrostatic penetration correction will become a standard component of future molecular mechanics force fields.
Collapse
Affiliation(s)
- Qiantao Wang
- Department of Biomedical Engineering and Division of Medicinal Chemistry, College of Pharmacy, The University of Texas at Austin , Austin, Texas 78712, United States ; Department of Biomedical Engineering and Division of Medicinal Chemistry, College of Pharmacy, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Joshua A Rackers
- Computational and Molecular Biophysics Program, Division of Biology & Biomedical Sciences, Washington University in St. Louis , St. Louis, Missouri 63110, United States
| | - Chenfeng He
- Department of Biomedical Engineering and Division of Medicinal Chemistry, College of Pharmacy, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Rui Qi
- Department of Biomedical Engineering and Division of Medicinal Chemistry, College of Pharmacy, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Christophe Narth
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC Paris 06, UMR 7616 , Case Courrier 137, 4 Place Jussieu, F-75005 Paris, France
| | - Louis Lagardere
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC Paris 06, UMR 7616 , Case Courrier 137, 4 Place Jussieu, F-75005 Paris, France
| | - Nohad Gresh
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC Paris 06, UMR 7616 , Case Courrier 137, 4 Place Jussieu, F-75005 Paris, France
| | - Jay W Ponder
- Department of Chemistry, Washington University in St. Louis , St. Louis, Missouri 63130, 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 and Division of Medicinal Chemistry, College of Pharmacy, The University of Texas at Austin , Austin, Texas 78712, United States
| |
Collapse
|
33
|
Odoh SO, Cramer CJ, Truhlar DG, Gagliardi L. Quantum-Chemical Characterization of the Properties and Reactivities of Metal–Organic Frameworks. Chem Rev 2015; 115:6051-111. [DOI: 10.1021/cr500551h] [Citation(s) in RCA: 206] [Impact Index Per Article: 22.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Samuel O. Odoh
- Department of Chemistry,
Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Christopher J. Cramer
- Department of Chemistry,
Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G. Truhlar
- Department of Chemistry,
Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Laura Gagliardi
- Department of Chemistry,
Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
34
|
Laury ML, Wang LP, Pande VS, Head-Gordon T, Ponder JW. Revised Parameters for the AMOEBA Polarizable Atomic Multipole Water Model. J Phys Chem B 2015; 119:9423-9437. [PMID: 25683601 DOI: 10.1021/jp510896n] [Citation(s) in RCA: 154] [Impact Index Per Article: 17.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
A set of improved parameters for the AMOEBA polarizable atomic multipole water model is developed. An automated procedure, ForceBalance, is used to adjust model parameters to enforce agreement with ab initio-derived results for water clusters and experimental data for a variety of liquid phase properties across a broad temperature range. The values reported here for the new AMOEBA14 water model represent a substantial improvement over the previous AMOEBA03 model. The AMOEBA14 model accurately predicts the temperature of maximum density and qualitatively matches the experimental density curve across temperatures from 249 to 373 K. Excellent agreement is observed for the AMOEBA14 model in comparison to experimental properties as a function of temperature, including the second virial coefficient, enthalpy of vaporization, isothermal compressibility, thermal expansion coefficient, and dielectric constant. The viscosity, self-diffusion constant, and surface tension are also well reproduced. In comparison to high-level ab initio results for clusters of 2-20 water molecules, the AMOEBA14 model yields results similar to AMOEBA03 and the direct polarization iAMOEBA models. With advances in computing power, calibration data, and optimization techniques, we recommend the use of the AMOEBA14 water model for future studies employing a polarizable water model.
Collapse
Affiliation(s)
- Marie L Laury
- Department of Chemistry, Washington University in St. Louis, St. Louis, MO 63105
| | - Lee-Ping Wang
- Department of Chemistry, Stanford University, Stanford, CA 94305
| | - Vijay S Pande
- Department of Chemistry, Stanford University, Stanford, CA 94305
| | | | - Jay W Ponder
- Department of Chemistry, Washington University in St. Louis, St. Louis, MO 63105
| |
Collapse
|