1
|
Cheng Z, Bi H, Liu S, Chen J, Misquitta AJ, Yu K. Developing a Differentiable Long-Range Force Field for Proteins with E(3) Neural Network-Predicted Asymptotic Parameters. J Chem Theory Comput 2024; 20:5598-5608. [PMID: 38888427 DOI: 10.1021/acs.jctc.4c00337] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/20/2024]
Abstract
Accurately describing long-range interactions is a significant challenge in molecular dynamics (MD) simulations of proteins. High-quality long-range potential is also an important component of the range-separated machine learning force field. This study introduces a comprehensive asymptotic parameter database encompassing atomic multipole moments, polarizabilities, and dispersion coefficients. Leveraging active learning, our database comprehensively represents protein fragments with up to 8 heavy atoms, capturing their conformational diversity with merely 78,000 data points. Additionally, the E(3) neural network (E3NN) is employed to predict the asymptotic parameters directly from the local geometry. The E3NN models demonstrate exceptional accuracy and transferability across all asymptotic parameters, achieving an R2 of 0.999 for both protein fragments and 20 amino acid dipeptide test sets. The long-range electrostatic and dispersion energies can be obtained using the E3NN-predicted parameters, with an error of 0.07 and 0.02 kcal/mol, respectively, when compared to symmetry-adapted perturbation theory (SAPT). Therefore, our force fields demonstrate the capability to accurately describe long-range interactions in proteins, paving the way for next-generation protein force fields.
Collapse
Affiliation(s)
- Zheng Cheng
- School of Mathematical Sciences, Peking University, Beijing 100871, China
- AI for Science Institute, Beijing 100084, P. R. China
| | - Hangrui Bi
- School of Mathematical Sciences, Peking University, Beijing 100871, China
- DP Technology, Beijing 100080, P. R. China
| | - Siyuan Liu
- DP Technology, Beijing 100080, P. R. China
| | - Junmin Chen
- Tsinghua-Berkeley Shenzhen Institute, Shenzhen 518055, Guangdong, P. R. China
- Tsinghua Shenzhen International Graduate School, Shenzhen 518055, Guangdong, P. R. China
| | - Alston J Misquitta
- School of Physics and Astronomy, Queen Mary, University of London, London E1 4NS, U.K
| | - Kuang Yu
- Tsinghua-Berkeley Shenzhen Institute, Shenzhen 518055, Guangdong, P. R. China
- Tsinghua Shenzhen International Graduate School, Shenzhen 518055, Guangdong, P. R. China
| |
Collapse
|
2
|
Chung MKJ, Ponder JW. Beyond isotropic repulsion: Classical anisotropic repulsion by inclusion of p orbitals. J Chem Phys 2024; 160:174118. [PMID: 38748037 PMCID: PMC11078554 DOI: 10.1063/5.0203678] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Accepted: 04/15/2024] [Indexed: 05/19/2024] Open
Abstract
Accurate modeling of intermolecular repulsion is an integral component in force field development. Although repulsion can be explicitly calculated by applying the Pauli exclusion principle, this approach is computationally viable only for systems of limited sizes. Instead, it has previously been shown that repulsion can be reformulated in a "classical" picture: the Pauli exclusion principle prohibits electrons from occupying the same state, leading to a depletion of electronic charge between atoms, giving rise to an enhanced nuclear-nuclear electrostatic repulsion. This classical picture is called the isotropic S2/R approximation, where S is the overlap and R is the interatomic distance. This approximation accurately captures the repulsion of isotropic atoms such as noble gas dimers; however, a key deficiency is that it fails to capture the angular dependence of the repulsion of anisotropic molecules. To include directionality, the wave function must at least be a linear combination of s and p orbitals. We derive a new anisotropic S2/R repulsion model through the inclusion of the anisotropic p orbital term in the total wave function. Because repulsion is pairwise and decays rapidly, it can be truncated at a short range, making it amenable for efficient calculation of energy and forces in complex biomolecular systems. We present a parameterization of the S101 dimer database against the ab initio benchmark symmetry-adapted perturbation theory, which yields an rms error of only 0.9 kcal/mol. The importance of the anisotropic term is demonstrated through angular scans of water-water dimers and dimers involving halobenzene. Simulation of liquid water shows that the model can be computed efficiently for realistic system sizes.
Collapse
Affiliation(s)
| | - Jay W. Ponder
- Author to whom correspondence should be addressed: . Tel.: 314-935-4275
| |
Collapse
|
3
|
Chen J, Yu K. PhyNEO: A Neural-Network-Enhanced Physics-Driven Force Field Development Workflow for Bulk Organic Molecule and Polymer Simulations. J Chem Theory Comput 2024; 20:253-265. [PMID: 38118076 DOI: 10.1021/acs.jctc.3c01045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2023]
Abstract
An accurate, generalizable, and transferable force field plays a crucial role in the molecular dynamics simulations of organic polymers and biomolecules. Conventional empirical force fields often fail to capture precise intermolecular interactions due to their negligence of important physics, such as polarization, charge penetration, many-body dispersion, etc. Moreover, the parameterization of these force fields relies heavily on top-down fittings, limiting their transferabilities to new systems where the experimental data are often unavailable. To address these challenges, we introduce a general and fully ab initio force field construction strategy, named PhyNEO. It features a hybrid approach that combines both the physics-driven and the data-driven methods and is able to generate a bulk potential with chemical accuracy using only quantum chemistry data of very small clusters. Careful separations of long-/short-range interactions and nonbonding/bonding interactions are the key to the success of PhyNEO. By such a strategy, we mitigate the limitations of pure data-driven methods in long-range interactions, thus largely increasing the data efficiency and the scalability of machine learning models. The new approach is thoroughly tested on poly(ethylene oxide) and polyethylene glycol systems, giving superior accuracies in both microscopic and bulk properties compared to conventional force fields. This work thus offers a promising framework for the development of advanced force fields in a wide range of organic molecular systems.
Collapse
Affiliation(s)
- Junmin Chen
- Tsinghua-Berkeley Shenzhen Institute, Tsinghua University, Shenzhen, Guangdong 518055, P. R. China
| | - Kuang Yu
- Tsinghua-Berkeley Shenzhen Institute, Tsinghua University, Shenzhen, Guangdong 518055, P. R. China
- Institute of Materials Research (iMR), Tsinghua Shenzhen International Graduate School, Tsinghua University, Shenzhen, Guangdong 518055, P. R. China
| |
Collapse
|
4
|
Tillotson MJ, Diamantonis NI, Buda C, Bolton LW, Müller EA. Molecular modelling of the thermophysical properties of fluids: expectations, limitations, gaps and opportunities. Phys Chem Chem Phys 2023; 25:12607-12628. [PMID: 37114325 DOI: 10.1039/d2cp05423j] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
This manuscript provides an overview of the current state of the art in terms of the molecular modelling of the thermophysical properties of fluids. It is intended to manage the expectations and serve as guidance to practising physical chemists, chemical physicists and engineers in terms of the scope and accuracy of the more commonly available intermolecular potentials along with the peculiarities of the software and methods employed in molecular simulations while providing insights on the gaps and opportunities available in this field. The discussion is focused around case studies which showcase both the precision and the limitations of frequently used workflows.
Collapse
Affiliation(s)
- Marcus J Tillotson
- Department of Chemical Engineering, Imperial College London, London, UK.
| | | | | | | | - Erich A Müller
- Department of Chemical Engineering, Imperial College London, London, UK.
| |
Collapse
|
5
|
Janicki TD, Van Vleet MJ, Schmidt JR. Development and Implementation of Atomically Anisotropic First-Principles Force Fields: A Benzene Case Study. J Phys Chem A 2023; 127:1736-1749. [PMID: 36780209 DOI: 10.1021/acs.jpca.2c07244] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/14/2023]
Abstract
π-interactions are an important motif in chemical and biochemical systems. However, due to their anisotropic electron densities and complex balance of intermolecular interactions, aromatic molecules represent an ongoing challenge for accurate and transferable force field development. Historically, ab initio force fields for aromatics have not exhibited good accuracy with respect to bulk properties or have only been used to study gas-phase dimers. Using benzene as a proof of concept, herein we show how our own ab initio MASTIFF force field incorporates an atomically anisotropic description of intermolecular interactions to yield an accurate and robust model for aromatic interactions irrespective of phase. Compared to existing models, the MASTIFF benzene force field not only is accurate for liquid phase properties but also offers transferability to the gas and solid phases. Additionally, we introduce a computationally efficient OpenMM plugin which enables customizable anisotropic intermolecular functional forms and which can be generically used in any MD simulation where a model for nonspherical atomic features is required. Overall, our results demonstrate the importance of atomic-level anisotropy in enabling next-generation ab initio force field development.
Collapse
Affiliation(s)
- Tesia D Janicki
- Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Mary J Van Vleet
- Department of Chemistry and Biochemistry, Spelman College, 350 Spelman Ln SW, Atlanta, Georgia 30314, United States
| | - J R Schmidt
- Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
6
|
Brown M, Skelton JM, Popelier PLA. Construction of a Gaussian Process Regression Model of Formamide for Use in Molecular Simulations. J Phys Chem A 2023; 127:1702-1714. [PMID: 36756842 PMCID: PMC9969515 DOI: 10.1021/acs.jpca.2c06566] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
FFLUX, a novel force field based on quantum chemical topology, can perform molecular dynamics simulations with flexible multipole moments that change with geometry. This is enabled by Gaussian process regression machine learning models, which accurately predict atomic energies and multipole moments up to the hexadecapole. We have constructed a model of the formamide monomer at the B3LYP/aug-cc-pVTZ level of theory capable of sub-kJ mol-1 accuracy, with the maximum prediction error for the molecule being 0.8 kJ mol-1. This model was used in FFLUX simulations along with Lennard-Jones parameters to successfully optimize the geometry of formamide dimers with errors smaller than 0.1 Å compared to those obtained with D3-corrected B3LYP/aug-cc-pVTZ. Comparisons were also made to a force field constructed with static multipole moments and Lennard-Jones parameters. FFLUX recovers the expected energy ranking of dimers compared to the literature, and changes in C═O and C-N bond lengths associated with hydrogen bonding were found to be consistent with density functional theory.
Collapse
|
7
|
James A, Swathi RS. Modeling the Adsorption of Polycyclic Aromatic Hydrocarbons on Graphynes: An Improved Lennard-Jones Formulation. J Phys Chem A 2022; 126:3472-3485. [PMID: 35609299 DOI: 10.1021/acs.jpca.2c01777] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Research on the development of theoretical methodologies for modeling noncovalent interactions governing the adsorption of polycyclic aromatic hydrocarbons (PAHs) on graphene and other two-dimensional materials is being intensely pursued in recent times. Highly accurate empirical potentials have emerged as a viable alternative to first-principles calculations for performing large-scale simulations. Herein, we report exploration of the potential energy surfaces for the adsorption of cata-condensed and peri-condensed PAHs on graphynes (GYs) using the improved Lennard-Jones (ILJ) potential. Initially, the ILJ potential is parametrized against benchmark electronic structure calculations performed on a selected set of PAH-GY complexes using dispersion-corrected density functional theory. The accuracy of the parametrization scheme is then assessed by a comparison of the adsorption features predicted from the ILJ potential with those computed using electronic structure calculations. The potential energy profiles as well as the single point energy calculations and geometry reoptimizations performed on the minimum-energy configurations predicted by the ILJ potential for a broader range of PAH-GY complexes provided a validation of the parametrization scheme. Finally, by an extrapolation of the PAH adsorption energies on various GYs, we estimated the interlayer cohesion energies for the van der Waals bilayer heterostructures of GYs with graphene to be in the range of 25-50 meV/atom.
Collapse
Affiliation(s)
- Anto James
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram (IISER TVM), Vithura, Thiruvananthapuram 695551, India
| | - Rotti Srinivasamurthy Swathi
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram (IISER TVM), Vithura, Thiruvananthapuram 695551, India
| |
Collapse
|
8
|
Das AK, Liu M, Head-Gordon T. Development of a Many-Body Force Field for Aqueous Alkali Metal and Halogen Ions: An Energy Decomposition Analysis Guided Approach. J Chem Theory Comput 2022; 18:953-967. [DOI: 10.1021/acs.jctc.1c00537] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Akshaya Kumar Das
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
| | - Meili Liu
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
- Department of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Teresa Head-Gordon
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
- Departments of Bioengineering and Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
9
|
Glick ZL, Koutsoukas A, Cheney DL, Sherrill CD. Cartesian message passing neural networks for directional properties: Fast and transferable atomic multipoles. J Chem Phys 2021; 154:224103. [PMID: 34241239 DOI: 10.1063/5.0050444] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
The message passing neural network (MPNN) framework is a promising tool for modeling atomic properties but is, until recently, incompatible with directional properties, such as Cartesian tensors. We propose a modified Cartesian MPNN (CMPNN) suitable for predicting atom-centered multipoles, an essential component of ab initio force fields. The efficacy of this model is demonstrated on a newly developed dataset consisting of 46 623 chemical structures and corresponding high-quality atomic multipoles, which was deposited into the publicly available Molecular Sciences Software Institute QCArchive server. We show that the CMPNN accurately predicts atom-centered charges, dipoles, and quadrupoles and that errors in the predicted atomic multipoles have a negligible effect on multipole-multipole electrostatic energies. The CMPNN is accurate enough to model conformational dependencies of a molecule's electronic structure. This opens up the possibility of recomputing atomic multipoles on the fly throughout a simulation in which they might exhibit strong conformational dependence.
Collapse
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
| | - Alexios Koutsoukas
- 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
|
10
|
Schriber JB, Nascimento DR, Koutsoukas A, Spronk SA, Cheney DL, Sherrill CD. CLIFF: A component-based, machine-learned, intermolecular force field. J Chem Phys 2021; 154:184110. [PMID: 34241025 DOI: 10.1063/5.0042989] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Computation of intermolecular interactions is a challenge in drug discovery because accurate ab initio techniques are too computationally expensive to be routinely applied to drug-protein models. Classical force fields are more computationally feasible, and force fields designed to match symmetry adapted perturbation theory (SAPT) interaction energies can remain accurate in this context. Unfortunately, the application of such force fields is complicated by the laborious parameterization required for computations on new molecules. Here, we introduce the component-based machine-learned intermolecular force field (CLIFF), which combines accurate, physics-based equations for intermolecular interaction energies with machine-learning models to enable automatic parameterization. The CLIFF uses functional forms corresponding to electrostatic, exchange-repulsion, induction/polarization, and London dispersion components in SAPT. Molecule-independent parameters are fit with respect to SAPT2+(3)δMP2/aug-cc-pVTZ, and molecule-dependent atomic parameters (atomic widths, atomic multipoles, and Hirshfeld ratios) are obtained from machine learning models developed for C, N, O, H, S, F, Cl, and Br. The CLIFF achieves mean absolute errors (MAEs) no worse than 0.70 kcal mol-1 in both total and component energies across a diverse dimer test set. For the side chain-side chain interaction database derived from protein fragments, the CLIFF produces total interaction energies with an MAE of 0.27 kcal mol-1 with respect to reference data, outperforming similar and even more expensive methods. In applications to a set of model drug-protein interactions, the CLIFF is able to accurately rank-order ligand binding strengths and achieves less than 10% error with respect to SAPT reference values for most complexes.
Collapse
Affiliation(s)
- Jeffrey B Schriber
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| | - Daniel R Nascimento
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30318, 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 30318, USA
| |
Collapse
|
11
|
Aina AA, Misquitta AJ, Price SL. A non-empirical intermolecular force-field for trinitrobenzene and its application in crystal structure prediction. J Chem Phys 2021; 154:094123. [PMID: 33685142 DOI: 10.1063/5.0043746] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An anisotropic atom-atom distributed intermolecular force-field (DIFF) for rigid trinitrobenzene (TNB) is developed using distributed multipole moments, dipolar polarizabilities, and dispersion coefficients derived from the charge density of the isolated molecule. The short-range parameters of the force-field are fitted to first- and second-order symmetry-adapted perturbation theory dimer interaction energy calculations using the distributed density-overlap model to guide the parameterization of the short-range anisotropy. The second-order calculations are used for fitting the damping coefficients of the long-range dispersion and polarization and also for relaxing the isotropic short-range coefficients in the final model, DIFF-srL2(rel). We assess the accuracy of the unrelaxed model, DIFF-srL2(norel), and its equivalent without short-range anisotropy, DIFF-srL0(norel), as these models are easier to derive. The model potentials are contrasted with empirical models for the repulsion-dispersion fitted to organic crystal structures with multipoles of iterated stockholder atoms (ISAs), FIT(ISA,L4), and with Gaussian Distributed Analysis (GDMA) multipoles, FIT(GDMA,L4), commonly used in modeling organic crystals. The potentials are tested for their ability to model the solid state of TNB. The non-empirical models provide more reasonable relative lattice energies of the three polymorphs of TNB and propose more sensible hypothetical structures than the empirical force-field (FIT). The DIFF-srL2(rel) model successfully has the most stable structure as one of the many structures that match the coordination sphere of form III. The neglect of the conformational flexibility of the nitro-groups is a significant approximation. This methodology provides a step toward force-fields capable of representing all phases of a molecule in molecular dynamics simulations.
Collapse
Affiliation(s)
- Alex A Aina
- Department of Chemistry, University College London, 20 Gordon St., London WC1H 0AJ, United Kingdom
| | - Alston J Misquitta
- School of Physics and Astronomy and The Thomas Young Centre for Theory and Simulation of Materials at Queen Mary, University of London, London E1 4NS, United Kingdom
| | - Sarah L Price
- Department of Chemistry, University College London, 20 Gordon St., London WC1H 0AJ, United Kingdom
| |
Collapse
|
12
|
Riera M, Hirales A, Ghosh R, Paesani F. Data-Driven Many-Body Models with Chemical Accuracy for CH4/H2O Mixtures. J Phys Chem B 2020; 124:11207-11221. [DOI: 10.1021/acs.jpcb.0c08728] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Marc Riera
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Alan Hirales
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Raja Ghosh
- 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
|
13
|
Abstract
Solid-state polymer electrolytes and high-concentration liquid electrolytes, such as water-in-salt electrolytes and ionic liquids, are emerging materials to replace the flammable organic electrolytes widely used in industrial lithium-ion batteries. Extensive efforts have been made to understand the ion transport mechanisms and optimize the ion transport properties. This perspective reviews the current understanding of the ion transport and polymer dynamics in liquid and polymer electrolytes, comparing the similarities and differences in the two types of electrolytes. Combining recent experimental and theoretical findings, we attempt to connect and explain ion transport mechanisms in different types of small-molecule and polymer electrolytes from a theoretical perspective, linking the macroscopic transport coefficients to the microscopic, molecular properties such as the solvation environment of the ions, salt concentration, solvent/polymer molecular weight, ion pairing, and correlated ion motion. We emphasize universal features in the ion transport and polymer dynamics by highlighting the relevant time and length scales. Several outstanding questions and anticipated developments for electrolyte design are discussed, including the negative transference number, control of ion transport through precision synthesis, and development of predictive multiscale modeling approaches.
Collapse
Affiliation(s)
- Chang Yun Son
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Zhen-Gang Wang
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
14
|
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
|
15
|
Garcia J, Podeszwa R, Szalewicz K. SAPT codes for calculations of intermolecular interaction energies. J Chem Phys 2020; 152:184109. [PMID: 32414261 DOI: 10.1063/5.0005093] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Symmetry-adapted perturbation theory (SAPT) is a method for calculations of intermolecular (noncovalent) interaction energies. The set of SAPT codes that is described here, the current version named SAPT2020, includes virtually all variants of SAPT developed so far, among them two-body SAPT based on perturbative, coupled cluster, and density functional theory descriptions of monomers, three-body SAPT, and two-body SAPT for some classes of open-shell monomers. The properties of systems governed by noncovalent interactions can be predicted only if potential energy surfaces (force fields) are available. SAPT is the preferred approach for generating such surfaces since it is seamlessly connected to the asymptotic expansion of interaction energy. SAPT2020 includes codes for automatic development of such surfaces, enabling generation of complete dimer surfaces with a rigid monomer approximation for dimers containing about one hundred atoms. These codes can also be used to obtain surfaces including internal degrees of freedom of monomers.
Collapse
Affiliation(s)
- Javier Garcia
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Rafał Podeszwa
- Institute of Chemistry, University of Silesia at Katowice, Szkolna 9, Katowice, Poland
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
16
|
Kantonen SM, Muddana HS, Schauperl M, Henriksen NM, Wang LP, Gilson MK. Data-Driven Mapping of Gas-Phase Quantum Calculations to General Force Field Lennard-Jones Parameters. J Chem Theory Comput 2020; 16:1115-1127. [PMID: 31917572 PMCID: PMC7101068 DOI: 10.1021/acs.jctc.9b00713] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular dynamics simulations are helpful tools for a range of applications, ranging from drug discovery to protein structure determination. The successful use of this technology largely depends on the potential function, or force field, used to determine the potential energy at each configuration of the system. Most force fields encode all of the relevant parameters to be used in distinct atom types, each associated with parameters for all parts of the force field, typically bond stretches, angle bends, torsions, and nonbonded terms accounting for van der Waals and electrostatic interactions. Much attention has been paid to the nonbonded parameters and their derivation, which are important in particular due to their governance of noncovalent interactions, such as protein-ligand binding. Parametrization involves adjusting the nonbonded parameters to minimize the error between simulation results and experimental properties, such as heats of vaporization and densities of neat liquids. In this setting, determining the best set of atom types is far from trivial, and the large number of parameters to be fit for the atom types in a typical force field can make it difficult to approach a true optimum. Here, we utilize a previously described Minimal Basis Iterative Stockholder (MBIS) method to carry out an atoms-in-molecules partitioning of electron densities. Information from these atomic densities is then mapped to Lennard-Jones parameters using a set of mapping parameters much smaller than the typical number of atom types in a force field. This approach is advantageous in two ways: it eliminates atom types by allowing each atom to have unique Lennard-Jones parameters, and it greatly reduces the number of parameters to be optimized. We show that this approach yields results comparable to those obtained with the typed GAFF 1.7 force field, even when trained on a relatively small amount of experimental data.
Collapse
Affiliation(s)
- Sophie M Kantonen
- Skaggs School of Pharmacy and Pharmaceutical Sciences , University of California San Diego , 9500 Gilman Drive , La Jolla , California 92093-0736 , United States
| | - Hari S Muddana
- Skaggs School of Pharmacy and Pharmaceutical Sciences , University of California San Diego , 9500 Gilman Drive , La Jolla , California 92093-0736 , United States
- OpenEye Scientific Software, Inc. , 9 Bisbee Court, Suite D , Santa Fe , New Mexico 87508 , United States
| | - Michael Schauperl
- Skaggs School of Pharmacy and Pharmaceutical Sciences , University of California San Diego , 9500 Gilman Drive , La Jolla , California 92093-0736 , United States
| | - Niel M Henriksen
- Skaggs School of Pharmacy and Pharmaceutical Sciences , University of California San Diego , 9500 Gilman Drive , La Jolla , California 92093-0736 , United States
- AtomWise, Inc. , 717 Market Street, Suite 800 , San Francisco , California 94103 , United States
| | - Lee-Ping Wang
- Department of Chemistry , University of California Davis , One Shields Avenue , Davis , California 95616 , United States
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences , University of California San Diego , 9500 Gilman Drive , La Jolla , California 92093-0736 , United States
| |
Collapse
|
17
|
Gilmore RAJ, Dove MT, Misquitta AJ. First-Principles Many-Body Nonadditive Polarization Energies from Monomer and Dimer Calculations Only: A Case Study on Water. J Chem Theory Comput 2020; 16:224-242. [PMID: 31769980 DOI: 10.1021/acs.jctc.9b00819] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The many-body polarization energy is the major source of nonadditivity in strongly polar systems such as water. This nonadditivity is often considerable and must be included, if only in an average manner, to correctly describe the physical properties of the system. Models for the polarization energy are usually parametrized using experimental data, or theoretical estimates of the many-body effects. Here we show how many-body polarization models can be developed for water complexes using data for the monomer and dimer only using ideas recently developed in the field of intermolecular perturbation theory and state-of-the-art approaches for calculating distributed molecular properties based on the iterated stockholder atoms (ISA) algorithm. We show how these models can be calculated, and we validate their accuracy in describing the many-body nonadditive energies of a range of water clusters. We further investigate their sensitivity to the details of the polarization damping models used. We show how our very best polarization models yield many-body energies that agree with those computed with coupled-cluster methods, but at a fraction of the computational cost.
Collapse
Affiliation(s)
- Rory A J Gilmore
- School of Physics and Astronomy and the Thomas Young Centre for Theory and Simulation of Materials at Queen Mary University of London , London E1 4NS , U.K
| | - Martin T Dove
- School of Physics and Astronomy and the Thomas Young Centre for Theory and Simulation of Materials at Queen Mary University of London , London E1 4NS , U.K
| | - Alston J Misquitta
- School of Physics and Astronomy and the Thomas Young Centre for Theory and Simulation of Materials at Queen Mary University of London , London E1 4NS , U.K
| |
Collapse
|
18
|
Patkowski K. Recent developments in symmetry‐adapted perturbation theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1452] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Konrad Patkowski
- Department of Chemistry and Biochemistry Auburn University Auburn Alabama
| |
Collapse
|
19
|
Jing Z, Liu C, Cheng SY, Qi R, Walker BD, Piquemal JP, Ren P. Polarizable Force Fields for Biomolecular Simulations: Recent Advances and Applications. Annu Rev Biophys 2019; 48:371-394. [PMID: 30916997 DOI: 10.1146/annurev-biophys-070317-033349] [Citation(s) in RCA: 217] [Impact Index Per Article: 43.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Realistic modeling of biomolecular systems requires an accurate treatment of electrostatics, including electronic polarization. Due to recent advances in physical models, simulation algorithms, and computing hardware, biomolecular simulations with advanced force fields at biologically relevant timescales are becoming increasingly promising. These advancements have not only led to new biophysical insights but also afforded opportunities to advance our understanding of fundamental intermolecular forces. This article describes the recent advances and applications, as well as future directions, of polarizable force fields in biomolecular simulations.
Collapse
Affiliation(s)
- Zhifeng Jing
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Sara Y Cheng
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Rui Qi
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Brandon D Walker
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Jean-Philip Piquemal
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA; .,Sorbonne Université, CNRS, Laboratoire de Chimie Theórique, 75252 Paris CEDEX 05, France.,Institut Universitaire de France, 75005 Paris, France
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| |
Collapse
|
20
|
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
|
21
|
Veit M, Jain SK, Bonakala S, Rudra I, Hohl D, Csányi G. Equation of State of Fluid Methane from First Principles with Machine Learning Potentials. J Chem Theory Comput 2019; 15:2574-2586. [DOI: 10.1021/acs.jctc.8b01242] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Max Veit
- Engineering Laboratory, University of Cambridge, Trumpington Street, Cambridge CB2 1PZ, United Kingdom
| | | | | | - Indranil Rudra
- Shell India Markets
Pvt. Ltd., Bengaluru 562149, Karnataka, India
| | - Detlef Hohl
- Shell Global Solutions
International BV, Grasweg 31, 1031 HW Amsterdam, The Netherlands
| | - Gábor Csányi
- Engineering Laboratory, University of Cambridge, Trumpington Street, Cambridge CB2 1PZ, United Kingdom
| |
Collapse
|
22
|
Misquitta AJ, Stone AJ. ISA-Pol: distributed polarizabilities and dispersion models from a basis-space implementation of the iterated stockholder atoms procedure. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2371-4] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
23
|
Bereau T, DiStasio RA, Tkatchenko A, von Lilienfeld OA. Non-covalent interactions across organic and biological subsets of chemical space: Physics-based potentials parametrized from machine learning. J Chem Phys 2018; 148:241706. [DOI: 10.1063/1.5009502] [Citation(s) in RCA: 109] [Impact Index Per Article: 18.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tristan Bereau
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Robert A. DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg,
Luxembourg
| | - O. Anatole von Lilienfeld
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel,
Switzerland
| |
Collapse
|