1
|
Metcalf DP, Smith AJ, Glick ZL, Sherrill CD. Range-dependence of two-body intermolecular interactions and their energy components in molecular crystals. J Chem Phys 2022; 157:084503. [DOI: 10.1063/5.0103644] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Routinely assessing the stability of molecular crystals with high accuracy remains an open challenge in the computational sciences. The many-body expansion decomposes computation of the crystal lattice energy into an embarrassingly parallel collection of computations over molecular dimers, trimers, and so forth, making quantum chemistry techniques tractable for many crystals of small organic molecules. By examining the range-dependence of different types of energetic contributions to the crystal lattice energy, we can glean qualitative understanding of solid-state intermolecular interactions as well as practical, exploitable reductions in the number of computations required for accurate energies. Here, we assess the range-dependent character of two-body interactions of 24 small organic molecular crystals using the physically interpretable components from symmetry-adapted perturbation theory (electrostatics, exchange repulsion, induction/polarization, and London dispersion). We also examine correlations between the convergence rates of electrostatics and London dispersion terms with molecular dipole moments and polarizabilities, to provide guidance for estimating convergence rates in other molecular crystals.
Collapse
Affiliation(s)
- Derek P Metcalf
- Chemistry & Biochemistry, Georgia Institute of Technology, United States of America
| | | | - Zachary Lee Glick
- Chemistry and Biochemistry, Georgia Institute of Technology College of Sciences, United States of America
| | - C. David Sherrill
- School of Chemistry and Biochemistry, Georgia Institute of Technology College of Sciences, United States of America
| |
Collapse
|
2
|
Konrad M, Wenzel W. CONI-Net: Machine Learning of Separable Intermolecular Force Fields. J Chem Theory Comput 2021; 17:4996-5006. [PMID: 34247485 DOI: 10.1021/acs.jctc.1c00328] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Noncovalent interactions (NCIs) play an essential role in soft matter and biomolecular simulations. The ab initio method symmetry-adapted perturbation theory allows a precise quantitative analysis of NCIs and offers an inherent energy decomposition, enabling a deeper understanding of the nature of intermolecular interactions. However, this method is limited to small systems, for instance, dimers of molecules. Here, we present a scale-bridging approach to systematically derive an intermolecular force field from ab initio data while preserving the energy decomposition of the underlying method. We apply the model in molecular dynamics simulations of several solvents and compare two predicted thermodynamic observables-mass density and enthalpy of vaporization-to experiments and established force fields. For a data set limited to hydrocarbons, we investigate the extrapolation capabilities to molecules absent from the training set. Overall, despite the affordable moderate quality of the reference ab initio data, we find promising results. With the straightforward data set generation procedure and the lack of target data in the fitting process, we have developed a method that enables the rapid development of predictive force fields with an extra dimension of insights into the balance of NCIs.
Collapse
Affiliation(s)
- Manuel Konrad
- Institute of Nanotechnology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, Eggenstein-Leopoldshafen 76344, Germany
| | - Wolfgang Wenzel
- Institute of Nanotechnology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, Eggenstein-Leopoldshafen 76344, Germany
| |
Collapse
|
3
|
Vilhena JG, Greff da Silveira L, Livotto PR, Cacelli I, Prampolini G. Automated Parameterization of Quantum Mechanically Derived Force Fields for Soft Materials and Complex Fluids: Development and Validation. J Chem Theory Comput 2021; 17:4449-4464. [PMID: 34185536 DOI: 10.1021/acs.jctc.1c00213] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The reliability of molecular dynamics (MD) simulations in predicting macroscopic properties of complex fluids and soft materials, such as liquid crystals, colloidal suspensions, or polymers, relies on the accuracy of the adopted force field (FF). We present an automated protocol to derive specific and accurate FFs, fully based on ab initio quantum mechanical (QM) data. The integration of the Joyce and Picky procedures, recently proposed by our group to provide an accurate description of simple liquids, is here extended to larger molecules, capable of exhibiting more complex fluid phases. While the standard Joyce protocol is employed to parameterize the intramolecular FF term, a new automated procedure is here proposed to handle the computational cost of the QM calculations required for the parameterization of the intermolecular FF term. The latter is thus obtained by integrating the old Picky procedure with a fragmentation reconstruction method (FRM) that allows for a reliable, yet computationally feasible sampling of the intermolecular energy surface at the QM level. The whole FF parameterization protocol is tested on a benchmark liquid crystal, and the performances of the resulting quantum mechanically derived (QMD) FF were compared with those delivered by a general-purpose, transferable one, and by the third, "hybrid" FF, where only the bonded terms were refined against QM data. Lengthy atomistic MD simulations are carried out with each FF on extended 5CB systems in both isotropic and nematic phases, eventually validating the proposed protocol by comparing the resulting macroscopic properties with other computational models and with experiments. The QMD-FF yields the best performances, reproducing both phases in the correct range of temperatures and well describing their structure, dynamics, and thermodynamic properties, thus providing a clear protocol that may be explored to predict such properties on other complex fluids or soft materials.
Collapse
Affiliation(s)
- J G Vilhena
- Department of Physics, University of Basel, Klingelbergstrasse 82, 4056 Basel, Switzerland
| | - Leandro Greff da Silveira
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Paolo Roberto Livotto
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Ivo Cacelli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici, ICCOM-CNR, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
4
|
Oliveira MP, Andrey M, Rieder SR, Kern L, Hahn DF, Riniker S, Horta BAC, Hünenberger PH. Systematic Optimization of a Fragment-Based Force Field against Experimental Pure-Liquid Properties Considering Large Compound Families: Application to Saturated Haloalkanes. J Chem Theory Comput 2020; 16:7525-7555. [DOI: 10.1021/acs.jctc.0c00683] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Marina P. Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Maurice Andrey
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Salomé R. Rieder
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Leyla Kern
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - David F. Hahn
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Bruno A. C. Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Philippe H. Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| |
Collapse
|
5
|
Glick ZL, Metcalf DP, Koutsoukas A, Spronk SA, Cheney DL, Sherrill CD. AP-Net: An atomic-pairwise neural network for smooth and transferable interaction potentials. J Chem Phys 2020; 153:044112. [DOI: 10.1063/5.0011521] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Affiliation(s)
- Zachary L. Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Derek P. Metcalf
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Alexios Koutsoukas
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Steven A. Spronk
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Daniel L. Cheney
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - C. David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
6
|
Abstract
A hydrogen bond (HB) is an essential interaction in countless phenomena, regulating the chemistry of life. HBs are characterized by two features, strength and directionality, with a high degree of heterogeneity across different chemical groups. These characteristics are dependent on the electronic configuration of the atoms involved in the interaction, which, in turn, is influenced strongly by the local molecular environment. Studies based on the analysis of HB in the solid phase, such as X-ray crystallography, suffer from significant biases due to packing forces. These will tend to better describe strong HBs at the expenses of weak ones, which will be either distorted or under-represented. Using quantum mechanics (QM), we calculated interaction energies for about a hundred acceptors and donors in a rigorously defined set of geometries. We performed 180,000 independent QM calculations, covering all relevant angular components, mapping strength and directionality in a context free from external biases, with both single-site and cooperative HBs. By quantifying directionality, we show that there is no correlation with strength; therefore, these two components need to be addressed separately. Results demonstrate that there are very strong HB acceptors (e.g., dimethyl sulfoxide) with nearly isotropic interactions and weak ones (e.g., thioacetone) with a sharp directional profile. Similarly, groups can have comparable directional propensity but be very distant in the strength spectrum (e.g., thioacetone and pyridine). Results provide a new perspective on the way HB directionality is described, with implications for biophysics and molecular recognition that ultimately can influence chemical biology, protein engineering, and drug design.
Collapse
Affiliation(s)
- Diogo Santos-Martins
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, California 92037, United States
| | - Stefano Forli
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, California 92037, United States
| |
Collapse
|
7
|
Derricotte WD. Symmetry-Adapted Perturbation Theory Decomposition of the Reaction Force: Insights into Substituent Effects Involved in Hemiacetal Formation Mechanisms. J Phys Chem A 2019; 123:7881-7891. [PMID: 31429558 DOI: 10.1021/acs.jpca.9b06865] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The decomposition of the reaction force based on symmetry-adapted perturbation theory (SAPT) has been proposed. This approach was used to investigate the substituent effects along the reaction coordinate pathway for the hemiacetal formation mechanism between methanol and substituted aldehydes of the form CX3CHO (X = H, F, Cl, and Br), providing a quantitative evaluation of the reaction-driving and reaction-retarding force components. Our results highlight the importance of more favorable electrostatic and induction effects in the reactions involving halogenated aldehydes that leads to lower activation energy barriers. These substituent effects are further elucidated by applying the functional-group partition of symmetry-adapted perturbation theory (F-SAPT). The results show that the reaction is largely driven by favorable direct noncovalent interactions between the CX3 group on the aldehyde and the OH group on methanol.
Collapse
Affiliation(s)
- Wallace D Derricotte
- Department of Chemistry , Morehouse College , Atlanta , Georgia 30314 , United States
| |
Collapse
|
8
|
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]
|
9
|
Rackers JA, Ponder JW. Classical Pauli repulsion: An anisotropic, atomic multipole model. J Chem Phys 2019; 150:084104. [PMID: 30823770 DOI: 10.1063/1.5081060] [Citation(s) in RCA: 40] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Pauli repulsion is a key component of any theory of intermolecular interactions. Although Pauli or exchange repulsion has its origin in the quantum mechanical nature of electrons, it is possible to describe the resulting energetic effects via a classical model in terms of the overlap of electron densities. In fact, closed shell intermolecular repulsion can be explained as a diminution of election density in the internuclear region resulting in decreased screening of nuclear charges and increased nuclear-nuclear repulsion. We provide a concise anisotropic repulsion formulation using the atomic multipoles from the Atomic Multipole Optimized Energetics for Biomolecular Applications force field to describe the electron density at each atom in a larger system. Mathematically, the proposed model consists of damped pairwise exponential multipolar repulsion interactions truncated at short range, which are suitable for use in compute-intensive biomolecular force fields and molecular dynamics simulations. Parameters for 26 atom classes encompassing most organic molecules are derived from a fit to Symmetry Adapted Perturbation Theory exchange repulsion energies for the S101 dimer database. Several applications of the multipolar Pauli repulsion model are discussed, including noble gas interactions, analysis of stationary points on the water dimer potential surface, and the directionality of several halogen bonding interactions.
Collapse
Affiliation(s)
- Joshua A Rackers
- Program in Computational and Molecular Biophysics, Washington University, School of Medicine, Saint Louis, Missouri 63110, USA
| | - Jay W Ponder
- Department of Chemistry, Washington University in Saint Louis, Saint Louis, Missouri 63130, USA
| |
Collapse
|
10
|
Hagler AT. Force field development phase II: Relaxation of physics-based criteria… or inclusion of more rigorous physics into the representation of molecular energetics. J Comput Aided Mol Des 2018; 33:205-264. [DOI: 10.1007/s10822-018-0134-x] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2018] [Accepted: 07/18/2018] [Indexed: 01/04/2023]
|
11
|
Gryn'ova G, Lin KH, Corminboeuf C. Read between the Molecules: Computational Insights into Organic Semiconductors. J Am Chem Soc 2018; 140:16370-16386. [PMID: 30395466 PMCID: PMC6287891 DOI: 10.1021/jacs.8b07985] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
![]()
The
performance and key electronic properties of molecular organic
semiconductors are dictated by the interplay between the chemistry
of the molecular core and the intermolecular factors of which manipulation
has inspired both experimentalists and theorists. This Perspective
presents major computational challenges and modern methodological
strategies to advance the field. The discussion ranges from insights
and design principles at the quantum chemical level, in-depth atomistic
modeling based on multiscale protocols, morphological prediction and
characterization as well as energy-property maps involving data-driven
analysis. A personal overview of the past achievements and future
direction is also provided.
Collapse
Affiliation(s)
- Ganna Gryn'ova
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland
| | - Kun-Han Lin
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland.,Laboratory for Computational Molecular Design and National Center for Computational Design and Discovery of Novel Materials (MARVEL) , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland
| | - Clémence Corminboeuf
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland.,Laboratory for Computational Molecular Design and National Center for Computational Design and Discovery of Novel Materials (MARVEL) , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland
| |
Collapse
|
12
|
Rackers JA, Liu C, Ren P, Ponder JW. A physically grounded damped dispersion model with particle mesh Ewald summation. J Chem Phys 2018; 149:084115. [PMID: 30193468 DOI: 10.1063/1.5030434] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Accurate modeling of dispersion is critical to the goal of predictive biomolecular simulations. To achieve this accuracy, a model must be able to correctly capture both the short-range and asymptotic behavior of dispersion interactions. We present here a damped dispersion model based on the overlap of charge densities that correctly captures both regimes. The overlap damped dispersion model represents a classical physical interpretation of dispersion: the interaction between the instantaneous induced dipoles of two distinct charge distributions. This model is shown to be an excellent fit with symmetry adapted perturbation theory dispersion energy calculations, yielding an RMS error on the S101x7 database of 0.5 kcal/mol. Moreover, the damping function used in this model is wholly derived and parameterized from the electrostatic dipole-dipole interaction, making it not only physically grounded but transferable as well.
Collapse
Affiliation(s)
- Joshua A Rackers
- Program in Computational and Molecular Biophysics, Washington University School of Medicine, Saint Louis, Missouri 63110, USA
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Jay W Ponder
- Program in Computational and Molecular Biophysics, Washington University School of Medicine, Saint Louis, Missouri 63110, USA
| |
Collapse
|
13
|
Greff da Silveira L, Jacobs M, Prampolini G, Livotto PR, Cacelli I. Development and Validation of Quantum Mechanically Derived Force-Fields: Thermodynamic, Structural, and Vibrational Properties of Aromatic Heterocycles. J Chem Theory Comput 2018; 14:4884-4900. [PMID: 30040902 DOI: 10.1021/acs.jctc.8b00218] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
A selection of several aromatic molecules, representative of the important class of heterocyclic compounds, has been considered for testing and validating an automated Force Field (FF) parametrization protocol, based only on Quantum Mechanical data. The parametrization is carried out separately for the intra- and intermolecular contributions, employing respectively the Joyce and Picky software packages, previously implemented and refined in our research group. The whole approach is here automated and integrated with a computationally effective yet accurate method, devised very recently ( J. Chem. THEORY Comput., 2018, 14, 543-556) to evaluate a large number of dimer interaction energies. The resulting quantum mechanically derived FFs are then used in extensive molecular dynamics simulations, in order to evaluate a number of thermodynamic, structural, and dynamic properties of the heterocycle's gas and liquid phases. The comparison with the available experimental data is good and furnishes a validation of the presented approach, which can be confidently exploited for the design of novel and more complex materials.
Collapse
Affiliation(s)
- Leandro Greff da Silveira
- Instituto de Química , Universidade Federal do Rio Grande do Sul , Avenida Bento Gonçalves 9500 , CEP 91501-970 Porto , Alegre , Brazil
| | - Matheus Jacobs
- Instituto de Química , Universidade Federal do Rio Grande do Sul , Avenida Bento Gonçalves 9500 , CEP 91501-970 Porto , Alegre , Brazil.,Institut für Physik , Humboldt-Universität zu Berlin , Newtonstrasse 15 , 12489 , Berlin , Germany.,IRIS Adelrshof , Humboldt-Universität zu Berlin , Zum Großen Windkanal 6 , 12489 , Berlin , Germany
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR) , Area della Ricerca, via G. Moruzzi 1 , I-56124 Pisa , Italy
| | - Paolo Roberto Livotto
- Instituto de Química , Universidade Federal do Rio Grande do Sul , Avenida Bento Gonçalves 9500 , CEP 91501-970 Porto , Alegre , Brazil
| | - Ivo Cacelli
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR) , Area della Ricerca, via G. Moruzzi 1 , I-56124 Pisa , Italy.,Dipartimento di Chimica e Chimica Industriale , Università di Pisa , Via G. Moruzzi 13 , I-56124 Pisa , Italy
| |
Collapse
|
14
|
Egan CK, Paesani F. Assessing Many-Body Effects of Water Self-Ions. I: OH–(H2O)n Clusters. J Chem Theory Comput 2018. [DOI: 10.1021/acs.jctc.7b01273] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Colin K. Egan
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
15
|
Ang SJ, Mak AM, Sullivan MB, Wong MW. Site specificity of halogen bonding involving aromatic acceptors. Phys Chem Chem Phys 2018. [DOI: 10.1039/c7cp08343b] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Based on Cambridge structural database survey and quantum chemical studies, aromatic halogen bond (XB) acceptors are found to have unique pattern of XB binding sites and rim specificity.
Collapse
Affiliation(s)
- Shi Jun Ang
- NUS Graduate School for Integrative Sciences and Engineering
- Centre for Life Sciences
- Singapore 117456
- Singapore
- Institute of High Performance Computing
| | - Adrian M. Mak
- Institute of High Performance Computing
- Singapore 138632
- Singapore
| | - Michael B. Sullivan
- Institute of High Performance Computing
- Singapore 138632
- Singapore
- Department of Chemistry
- National University of Singapore
| | - Ming Wah Wong
- NUS Graduate School for Integrative Sciences and Engineering
- Centre for Life Sciences
- Singapore 117456
- Singapore
- Department of Chemistry
| |
Collapse
|
16
|
Demerdash O, Wang L, Head‐Gordon T. Advanced models for water simulations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1355] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Omar Demerdash
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley CA USA
- Department of Chemistry University of California Berkeley CA USA
| | - Lee‐Ping Wang
- Department of Chemistry University of California, Davis Davis CA USA
| | - Teresa Head‐Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley CA USA
- Department of Chemistry University of California Berkeley CA USA
- Department of Bioengineering University of California Berkeley CA USA
- Department of Chemical and Biomolecular Engineering University of California Berkeley CA USA
| |
Collapse
|
17
|
Grimme S, Bannwarth C, Caldeweyher E, Pisarek J, Hansen A. A general intermolecular force field based on tight-binding quantum chemical calculations. J Chem Phys 2017; 147:161708. [DOI: 10.1063/1.4991798] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| | - Christoph Bannwarth
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| | - Eike Caldeweyher
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| | - Jana Pisarek
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| |
Collapse
|
18
|
Aina AA, Misquitta AJ, Price SL. From dimers to the solid-state: Distributed intermolecular force-fields for pyridine. J Chem Phys 2017; 147:161722. [DOI: 10.1063/1.4999789] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alexander A. Aina
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Alston J. Misquitta
- School of Physics and Astronomy, Queen Mary, University of London, London E1 4NS, United Kingdom
| | - Sarah L. Price
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| |
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
|
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
|
21
|
Cui Q. Perspective: Quantum mechanical methods in biochemistry and biophysics. J Chem Phys 2017; 145:140901. [PMID: 27782516 DOI: 10.1063/1.4964410] [Citation(s) in RCA: 60] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
In this perspective article, I discuss several research topics relevant to quantum mechanical (QM) methods in biophysical and biochemical applications. Due to the immense complexity of biological problems, the key is to develop methods that are able to strike the proper balance of computational efficiency and accuracy for the problem of interest. Therefore, in addition to the development of novel ab initio and density functional theory based QM methods for the study of reactive events that involve complex motifs such as transition metal clusters in metalloenzymes, it is equally important to develop inexpensive QM methods and advanced classical or quantal force fields to describe different physicochemical properties of biomolecules and their behaviors in complex environments. Maintaining a solid connection of these more approximate methods with rigorous QM methods is essential to their transferability and robustness. Comparison to diverse experimental observables helps validate computational models and mechanistic hypotheses as well as driving further development of computational methodologies.
Collapse
Affiliation(s)
- Qiang Cui
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, USA
| |
Collapse
|
22
|
Vandenbrande S, Waroquier M, Speybroeck VV, Verstraelen T. The Monomer Electron Density Force Field (MEDFF): A Physically Inspired Model for Noncovalent Interactions. J Chem Theory Comput 2016; 13:161-179. [DOI: 10.1021/acs.jctc.6b00969] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Steven Vandenbrande
- Center for Molecular Modeling
(CMM), QCMM Ghent−Brussels Alliance, Ghent University, Technologiepark
903, B9000 Ghent, Belgium
| | - Michel Waroquier
- Center for Molecular Modeling
(CMM), QCMM Ghent−Brussels Alliance, Ghent University, Technologiepark
903, B9000 Ghent, Belgium
| | - Veronique Van Speybroeck
- Center for Molecular Modeling
(CMM), QCMM Ghent−Brussels Alliance, Ghent University, Technologiepark
903, B9000 Ghent, Belgium
| | - Toon Verstraelen
- Center for Molecular Modeling
(CMM), QCMM Ghent−Brussels Alliance, Ghent University, Technologiepark
903, B9000 Ghent, Belgium
| |
Collapse
|
23
|
Prampolini G, Campetella M, De Mitri N, Livotto PR, Cacelli I. Systematic and Automated Development of Quantum Mechanically Derived Force Fields: The Challenging Case of Halogenated Hydrocarbons. J Chem Theory Comput 2016; 12:5525-5540. [DOI: 10.1021/acs.jctc.6b00705] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Marco Campetella
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi
13, I-56124 Pisa, Italy
| | - Nicola De Mitri
- Department
of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Paolo Roberto Livotto
- Instituto
de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Ivo Cacelli
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi
13, I-56124 Pisa, Italy
| |
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
|
Abstract
Hydrogen bond directionality in the water dimer is explained on the basis of symmetry-adapted intermolecular perturbation theory which directly separates the intermolecular interaction energy into four physically interpretable components: electrostatics, exchange-repulsion, dispersion, and induction. Analysis of these four main contributions to the binding energy allows a deeper understanding of the dominant factors ruling the mutual arrangement of the two monomers. A preference for the linear configuration is shown to be due to a subtle interplay of all four energy components. While the first-order terms, electrostatic and exchange-repulsion, almost perfectly cancel each other near the equilibrium geometry of the dimer, the importance of the second- and higher-order terms, induction and dispersion, becomes evident.
Collapse
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
| |
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
|