1
|
Hwang W, Austin SL, Blondel A, Boittier ED, Boresch S, Buck M, Buckner J, Caflisch A, Chang HT, Cheng X, Choi YK, Chu JW, Crowley MF, Cui Q, Damjanovic A, Deng Y, Devereux M, Ding X, Feig MF, Gao J, Glowacki DR, Gonzales JE, Hamaneh MB, Harder ED, Hayes RL, Huang J, Huang Y, Hudson PS, Im W, Islam SM, Jiang W, Jones MR, Käser S, Kearns FL, Kern NR, Klauda JB, Lazaridis T, Lee J, Lemkul JA, Liu X, Luo Y, MacKerell AD, Major DT, Meuwly M, Nam K, Nilsson L, Ovchinnikov V, Paci E, Park S, Pastor RW, Pittman AR, Post CB, Prasad S, Pu J, Qi Y, Rathinavelan T, Roe DR, Roux B, Rowley CN, Shen J, Simmonett AC, Sodt AJ, Töpfer K, Upadhyay M, van der Vaart A, Vazquez-Salazar LI, Venable RM, Warrensford LC, Woodcock HL, Wu Y, Brooks CL, Brooks BR, Karplus M. CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed. J Phys Chem B 2024. [PMID: 39303207 DOI: 10.1021/acs.jpcb.4c04100] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/22/2024]
Abstract
Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.
Collapse
Affiliation(s)
- Wonmuk Hwang
- Department of Biomedical Engineering, Texas A&M University, College Station, Texas 77843, United States
- Department of Materials Science and Engineering, Texas A&M University, College Station, Texas 77843, United States
- Department of Physics and Astronomy, Texas A&M University, College Station, Texas 77843, United States
- Center for AI and Natural Sciences, Korea Institute for Advanced Study, Seoul 02455, Republic of Korea
| | - Steven L Austin
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - Arnaud Blondel
- Institut Pasteur, Université Paris Cité, CNRS UMR3825, Structural Bioinformatics Unit, 28 rue du Dr. Roux F-75015 Paris, France
| | - Eric D Boittier
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Stefan Boresch
- Faculty of Chemistry, Department of Computational Biological Chemistry, University of Vienna, Wahringerstrasse 17, 1090 Vienna, Austria
| | - Matthias Buck
- Department of Physiology and Biophysics, Case Western Reserve University, School of Medicine, Cleveland, Ohio 44106, United States
| | - Joshua Buckner
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Amedeo Caflisch
- Department of Biochemistry, University of Zürich, CH-8057 Zürich, Switzerland
| | - Hao-Ting Chang
- Institute of Bioinformatics and Systems Biology, National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan, ROC
| | - Xi Cheng
- Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Yeol Kyo Choi
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Jhih-Wei Chu
- Institute of Bioinformatics and Systems Biology, Department of Biological Science and Technology, Institute of Molecular Medicine and Bioengineering, and Center for Intelligent Drug Systems and Smart Bio-devices (IDS2B), National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan, ROC
| | - Michael F Crowley
- Renewable Resources and Enabling Sciences Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Biomedical Engineering, Boston University, 44 Cummington Mall, Boston, Massachusetts 02215, United States
| | - Ana Damjanovic
- Department of Biophysics, Johns Hopkins University, Baltimore, Maryland 21218, United States
- Department of Physics and Astronomy, Johns Hopkins University, Baltimore, Maryland 21218, United States
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Yuqing Deng
- Shanghai R&D Center, DP Technology, Ltd., Shanghai 201210, China
| | - Mike Devereux
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Xinqiang Ding
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, United States
| | - Michael F Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Jiali Gao
- School of Chemical Biology & Biotechnology, Peking University Shenzhen Graduate School, Shenzhen, Guangdong 518055, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong 518055, China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - David R Glowacki
- CiTIUS Centro Singular de Investigación en Tecnoloxías Intelixentes da USC, 15705 Santiago de Compostela, Spain
| | - James E Gonzales
- Department of Biomedical Engineering, Texas A&M University, College Station, Texas 77843, United States
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Mehdi Bagerhi Hamaneh
- Department of Physiology and Biophysics, Case Western Reserve University, School of Medicine, Cleveland, Ohio 44106, United States
| | | | - Ryan L Hayes
- Department of Chemical and Biomolecular Engineering, University of California, Irvine, Irvine, California 92697, United States
- Department of Pharmaceutical Sciences, University of California, Irvine, Irvine, California 92697, United States
| | - Jing Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, Hangzhou, Zhejiang 310024, China
| | - Yandong Huang
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Phillip S Hudson
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
- Medicine Design, Pfizer Inc., Cambridge, Massachusetts 02139, United States
| | - Wonpil Im
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Shahidul M Islam
- Department of Chemistry, Delaware State University, Dover, Delaware 19901, United States
| | - Wei Jiang
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Michael R Jones
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Silvan Käser
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Fiona L Kearns
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - Nathan R Kern
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering, Institute for Physical Science and Technology, Biophysics Program, University of Maryland, College Park, Maryland 20742, United States
| | - Themis Lazaridis
- Department of Chemistry, City College of New York, New York, New York 10031, United States
| | - Jinhyuk Lee
- Disease Target Structure Research Center, Korea Research Institute of Bioscience and Biotechnology, Daejeon 34141, Republic of Korea
- Department of Bioinformatics, KRIBB School of Bioscience, University of Science and Technology, Daejeon 34141, Republic of Korea
| | - Justin A Lemkul
- Department of Biochemistry, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, United States
| | - Xiaorong Liu
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Yun Luo
- Department of Biotechnology and Pharmaceutical Sciences, College of Pharmacy, Western University of Health Sciences, Pomona, California 91766, United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, United States
| | - Dan T Major
- Department of Chemistry and Institute for Nanotechnology & Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, United States
| | - Lennart Nilsson
- Karolinska Institutet, Department of Biosciences and Nutrition, SE-14183 Huddinge, Sweden
| | - Victor Ovchinnikov
- Harvard University, Department of Chemistry and Chemical Biology, Cambridge, Massachusetts 02138, United States
| | - Emanuele Paci
- Dipartimento di Fisica e Astronomia, Universitá di Bologna, Bologna 40127, Italy
| | - Soohyung Park
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Richard W Pastor
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Amanda R Pittman
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - Carol Beth Post
- Borch Department of Medicinal Chemistry and Molecular Pharmacology, Purdue University, West Lafayette, Indiana 47907, United States
| | - Samarjeet Prasad
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Jingzhi Pu
- Department of Chemistry and Chemical Biology, Indiana University Indianapolis, Indianapolis, Indiana 46202, United States
| | - Yifei Qi
- School of Pharmacy, Fudan University, Shanghai 201203, China
| | | | - Daniel R Roe
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Benoit Roux
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | | | - Jana Shen
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, United States
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Alexander J Sodt
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Kai Töpfer
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Meenu Upadhyay
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Arjan van der Vaart
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | | | - Richard M Venable
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Luke C Warrensford
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - H Lee Woodcock
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - Yujin Wu
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Charles L Brooks
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Martin Karplus
- Harvard University, Department of Chemistry and Chemical Biology, Cambridge, Massachusetts 02138, United States
- Laboratoire de Chimie Biophysique, ISIS, Université de Strasbourg, 67000 Strasbourg, France
| |
Collapse
|
2
|
Musiał M, Riccardi D, Suiter CL, Sontarp EJ, Miller SL, Lirette RL, Rehmeier KC, Mahata A, Muzny CD, Stelson AC, Schwarz KA, Widegren JA. NMR Spectroscopy and Multiscale Modeling Shed Light on Ion-Solvent Interactions and Ion Pairing in Aqueous NaF Solutions. J Phys Chem B 2024; 128:8974-8983. [PMID: 39253766 DOI: 10.1021/acs.jpcb.4c03521] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/11/2024]
Abstract
The balance between ion solvation and ion pairing in aqueous solutions modulates chemical and physical processes from catalysis to protein folding. Yet, despite more than a century of investigation, experimental determination of the distribution of ion-solvation and ion-pairing states remains elusive, even for archetypal systems like aqueous alkali halides. Here, we combine nuclear magnetic resonance (NMR) spectroscopy and multiscale modeling to disentangle ion-solvent interactions from ion pairing in aqueous sodium fluoride solutions. We have developed a high-accuracy method to collect experimental NMR resonance frequencies for both ions as functions of temperature and concentration. Comparison of these data with resonance frequencies for nonassociating salts allows us to differentiate the influence of solvation and ion pairing on NMR spectra. These high-quality experimental NMR data are used to validate our modeling framework comprising polarizable force field molecular dynamics (MD) simulations and quantum chemical calculations of NMR resonance frequencies. Our experimental and theoretical resonance frequency shifts agree over a wide range of temperatures and concentrations. Structural analysis reveals how both trends are dominated by interactions with water molecules. For the more sensitive 19F nucleus, the NMR resonance frequency decreases as hydrogen bonds between fluoride and water molecules are reduced in number with increased temperature and molality. Through a detailed analysis of the theoretical NMR resonance frequencies for both ions, we show that NMR spectroscopy can distinguish both contact ion pairs and single-solvent-separated ion pairs from free ions. This quantitative framework can be applied directly to other systems.
Collapse
Affiliation(s)
- Małgorzata Musiał
- Department of Physics, University of Colorado, Boulder, Colorado 80309, United States
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, United States
| | - Demian Riccardi
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, United States
| | - Christopher L Suiter
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, United States
| | - Ethan J Sontarp
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, United States
- Department of Geosciences, Princeton University, Princeton, New Jersey 08544, United States
| | - Samantha L Miller
- Department of Physics, University of Colorado, Boulder, Colorado 80309, United States
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, United States
| | - Robert L Lirette
- Department of Physics, University of Colorado, Boulder, Colorado 80309, United States
- RF Technology Division, National Institute of Standards and Technology, Boulder, Colorado 80305, United States
| | - Kyle Covington Rehmeier
- Department of Physics, University of Colorado, Boulder, Colorado 80309, United States
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, United States
| | - Avik Mahata
- School of Engineering, Brown University, Providence, Rhode Island 02912, United States
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Chris D Muzny
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, United States
| | - Angela C Stelson
- RF Technology Division, National Institute of Standards and Technology, Boulder, Colorado 80305, United States
| | - Kathleen A Schwarz
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Jason A Widegren
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, United States
| |
Collapse
|
3
|
Yang Y, Jin Q, Yin S. Development of an anisotropic polarizable model for the all-atom AMOEBA force field. Phys Chem Chem Phys 2024; 26:22900-22911. [PMID: 39169824 DOI: 10.1039/d4cp01568a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2024]
Abstract
For planar and rigid π-conjugated molecular systems, electrostatic and inductive interactions are pivotal in governing molecular packing structures and electron polarization energies. These electrostatic interactions typically exhibit an anisotropic nature within π-conjugated systems. In this study, we utilize the atoms in molecules (AIM) theory in conjunction with linear response theory to decompose molecular polarizability into distributed atomic polarizability tensors. On the basis of atomic polarizability tensors, we extended an anisotropic polarizable model into the AMOEBA polarizable force field. Both anisotropic and isotropic polarizable models in combination with various density functional theory (DFT)-derived atomic multipoles were applied to optimize the experimental crystals of naphthalene and anthracene. Furthermore, these two types of electrostatic models, coupled with the evolutionary algorithm USPEX program, are utilized to predict the crystal structures of oligoacenes. Our findings demonstrate that the anisotropic polarizable model exhibits superior performance in crystal refinement and crystal structure prediction. This enriched anisotropic polarizable model is seamlessly integrated into the AMOEBA polarizable force field and readily applicable within our modified Tinker program.
Collapse
Affiliation(s)
- Yanyan Yang
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an City 710119, People's Republic of China.
| | - Qianqian Jin
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an City 710119, People's Republic of China.
| | - Shiwei Yin
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an City 710119, People's Republic of China.
| |
Collapse
|
4
|
Liao J, Wu M, Gao J, Chen C. Calculation of solvation force in molecular dynamics simulation by deep-learning method. Biophys J 2024; 123:2830-2838. [PMID: 38444159 PMCID: PMC11393703 DOI: 10.1016/j.bpj.2024.02.029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Revised: 01/31/2024] [Accepted: 02/29/2024] [Indexed: 03/07/2024] Open
Abstract
Electrostatic calculations are generally used in studying the thermodynamics and kinetics of biomolecules in solvent. Generally, this is performed by solving the Poisson-Boltzmann equation on a large grid system, a process known to be time consuming. In this study, we developed a deep neural network to predict the decomposed solvation free energies and forces of all atoms in a molecule. To train the network, the internal coordinates of the molecule were used as the input data, and the solvation free energies along with transformed atomic forces from the Poisson-Boltzmann equation were used as labels. Both the training and prediction tasks were accelerated on GPU. Formal tests demonstrated that our method can provide reasonable predictions for small molecules when the network is well-trained with its simulation data. This method is suitable for processing lots of snapshots of molecules in a long trajectory. Moreover, we applied this method in the molecular dynamics simulation with enhanced sampling. The calculated free energy landscape closely resembled that obtained from explicit solvent simulations.
Collapse
Affiliation(s)
- Jun Liao
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Mincong Wu
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Junyong Gao
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Changjun Chen
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China.
| |
Collapse
|
5
|
Li Y, Liu Y, Yang B, Li G, Chu H. Polarizable atomic multipole-based force field for cholesterol. J Biomol Struct Dyn 2024; 42:7747-7757. [PMID: 37565356 DOI: 10.1080/07391102.2023.2245045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2023] [Accepted: 07/21/2023] [Indexed: 08/12/2023]
Abstract
Cholesterol is one of the essential component of lipid in membrane. We present a polarizable atomic multipole force field (FF) for the molecular dynamic simulation of cholesterol. The FF building process follows the computational framework as the atomic multipole optimized energetics for biomolecular applications (AMOEBA) model. In this framework, the electronics parameters, including atomic monopole moments, dipole moments, and quadrupole moments calculated from ab initio calculations in the gas phase, are applied to represent the charge distribution. Furthermore, the many-body polarization is modeled by following the same pattern of distributed atomic polarizabilities. Then, the bilayers composed of two typical phospholipid molecules, 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), in a range of different cholesterol concentrations are built and implemented by molecular dynamics (MD) simulations based on the proposed polarizable FF. The simulation results are statistically analyzed to validate the feasibility of the proposed FF. The structural properties of the bilayers are calculated to compare with the related experimental values. The MD values show the same trend of experimental values changes.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Yan Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Dalian, Liaoning, China
| | - Ye Liu
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Dalian, Liaoning, China
| | - Boya Yang
- Dalian Municipal Central Hospital, Liaoning, China
| | - Guohui Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Dalian, Liaoning, China
| | - Huiying Chu
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Dalian, Liaoning, China
| |
Collapse
|
6
|
Blazhynska M, Lagardère L, Liu C, Adjoua O, Ren P, Piquemal JP. Water-glycan interactions drive the SARS-CoV-2 spike dynamics: insights into glycan-gate control and camouflage mechanisms. Chem Sci 2024:d4sc04364b. [PMID: 39220162 PMCID: PMC11359970 DOI: 10.1039/d4sc04364b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2024] [Accepted: 08/22/2024] [Indexed: 09/04/2024] Open
Abstract
To develop therapeutic strategies against COVID-19, we introduce a high-resolution all-atom polarizable model capturing many-body effects of protein, glycan, solvent, and membrane components in SARS-CoV-2 spike protein open and closed states. Employing μs-long molecular dynamics simulations powered by high-performance cloud-computing and unsupervised density-driven adaptive sampling, we investigated the differences in bulk-solvent-glycan and protein-solvent-glycan interfaces between these states. We unraveled a sophisticated solvent-glycan polarization interaction network involving the N165/N343 glycan-gate patterns that provide structural support for the open state and identified key water molecules that could potentially be targeted to destabilize this configuration. In the closed state, the reduced solvent polarization diminishes the overall N165/N343 dipoles, yet internal interactions and a reorganized sugar coat stabilize this state. Despite variations, our glycan-solvent accessibility analysis reveals the glycan shield capability to conserve constant interactions with the solvent, effectively camouflaging the virus from immune detection in both states. The presented insights advance our comprehension of viral pathogenesis at an atomic level, offering potential to combat COVID-19.
Collapse
Affiliation(s)
- Marharyta Blazhynska
- Laboratoire de Chimie Théorique, Sorbonne Université, UMR 7616 CNRS 75005 Paris France
| | - Louis Lagardère
- Laboratoire de Chimie Théorique, Sorbonne Université, UMR 7616 CNRS 75005 Paris France
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin Texas 78712 USA
- Qubit Pharmaceuticals 75014 Paris France
| | - Olivier Adjoua
- Laboratoire de Chimie Théorique, Sorbonne Université, UMR 7616 CNRS 75005 Paris France
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin Texas 78712 USA
| | - Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, Sorbonne Université, UMR 7616 CNRS 75005 Paris France
| |
Collapse
|
7
|
Bone RA, Chung MKJ, Ponder JW, Riccardi D, Muzny C, Sundararaman R, Schwarz K. A new method to calculate broadband dielectric spectra of solvents from molecular dynamics simulations demonstrated with polarizable force fields. J Chem Phys 2024; 161:064306. [PMID: 39132799 PMCID: PMC11324330 DOI: 10.1063/5.0217883] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2024] [Accepted: 07/26/2024] [Indexed: 08/13/2024] Open
Abstract
Simulating the dielectric spectra of solvents requires the nuanced definition of inter- and intra-molecular forces. Non-polarizable force fields, while thoroughly benchmarked for dielectric applications, do not capture all the spectral features of solvents, such as water. Conversely, polarizable force fields have been largely untested in the context of dielectric spectroscopy but include charge and dipole fluctuations that contribute to intermolecular interactions. We benchmark non-polarizable force fields and the polarizable force fields AMOEBA03 and HIPPO for liquid water and find that the polarizable force fields can capture all the experimentally observed spectral features with varying degrees of accuracy. However, the non-polarizable force fields miss at least one peak. To diagnose this deficiency, we decompose the liquid water spectra from polarizable force fields at multiple temperatures into static and induced dipole contributions and find that the peak originates from induced dipole contributions. Broadening our inquiry to other solvents parameterized with the AMOEBA09 force field, we demonstrate good agreement between the experimental and simulated dielectric spectra of methanol and formamide. To produce these spectra, we develop a new computational approach to calculate the dielectric spectrum via the fluctuation dissipation theorem. This method minimizes the error in both the low and high frequency portions of the spectrum, improving the overall accuracy of the simulated spectrum and broadening the computed frequency range.
Collapse
Affiliation(s)
| | - Moses K. J. Chung
- Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri 63130, USA
| | - Jay W. Ponder
- Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri 63130, USA
| | - Demian Riccardi
- Material Measurement Laboratory, National Institute of Standards and Technology, 325 Broadway, Boulder, Colorado 80305, USA
| | - Chris Muzny
- Material Measurement Laboratory, National Institute of Standards and Technology, 325 Broadway, Boulder, Colorado 80305, USA
| | - Ravishankar Sundararaman
- Department of Materials Science and Engineering, Rensselaer Polytechnic Institute, 110 8th St., Troy, New York 12180, USA
| | - Kathleen Schwarz
- Material Measurement Laboratory, National Institute of Standards and Technology, 100 Bureau Dr., Gaithersburg, Maryland 20899, USA
| |
Collapse
|
8
|
Mareš J, Mayorga Delgado P. Getting the intermolecular forces correct: introducing the ASTA strategy for a water model. RSC Adv 2024; 14:25712-25727. [PMID: 39148757 PMCID: PMC11325342 DOI: 10.1039/d4ra02685c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2024] [Accepted: 08/01/2024] [Indexed: 08/17/2024] Open
Abstract
Having a force field for water providing good bulk properties is paramount for modern studies of most biological systems. Some of the most common three-site force fields are TIP3, SPC/ε or OPC3, providing a decent range of bulk properties. That does not mean though, that they have realistic inter-atomic forces. These force fields have been parameterized with a top-down approach, meaning, by fitting the force field parameters to the experimental bulk properties. This approach has been the governing strategy also for many variants of four- and more-site models. We test a bottom-up approach, in which the force field is parameterized by optimizing the non-bonded inter-atomic forces. Our philosophy is that correct inter-atomic forces lead to correct geometrical and dynamical properties. The first system we try to optimize with the accurately system tailored atomic (ASTA) approach is water, but we aim to eventually probe other systems in the future as well. We applied our ASTA strategy to find a good set of parameters providing accurate bulk properties for the simple three-site force field forms, and also for AMOEBA, a more detailed and polarizable force field. Even though our bottom-up approach did not provide satisfactory results for the simple three-site force fields (with fixed charges), for the case of the AMOEBA force field it led to a modification of the original strategy, giving very good intra- and inter-molecular forces, as compared to accurate quantum chemically calculated reference forces. At the same time, important bulk properties, in this study restricted to the density and diffusion, were accurately reproduced with respect to the experimental values.
Collapse
Affiliation(s)
- Jiří Mareš
- Department of Physics, University of Oulu Finland
| | | |
Collapse
|
9
|
Elliott GR, Wanless EJ, Webber GB, Andersson GG, Craig VSJ, Page AJ. Dynamic Ion Correlations and Ion-Pair Lifetimes in Aqueous Alkali Metal Chloride Electrolytes. J Phys Chem B 2024; 128:7438-7444. [PMID: 39037039 DOI: 10.1021/acs.jpcb.4c01992] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/23/2024]
Abstract
Electrolytes are central to many technological applications, as well as life itself. The behavior and properties of electrolytes are often described in terms of ion pairs, whereby ions associate as either contact ion pairs (in which ions are "touching") solvent-separated ion pairs (in which ions' solvent shells overlap) or solvent-solvent-separated ion pairs (in which ions' solvent shells are distinct). However, this paradigm is generally restricted to statistically averaged descriptions of solution structure and ignores temporal behavior. Here we elucidate the time-resolved dynamics of these ion-ion interactions in aqueous metal chloride electrolytes using the partial van Hove correlation function, based on polarizable molecular dynamics simulations. Our results show that the existence and persistence of ion pairs in aqueous metal chloride electrolytes should not be assumed a priori, but in fact are ion specific features of the solution with lifetimes on subpicosecond time scales.
Collapse
Affiliation(s)
- Gareth R Elliott
- Discipline of Chemistry, The University of Newcastle, Callaghan, New South Wales 2308, Australia
| | - Erica J Wanless
- Discipline of Chemistry, The University of Newcastle, Callaghan, New South Wales 2308, Australia
| | - Grant B Webber
- Discipline of Chemical Engineering, The University of Newcastle, Callaghan, New South Wales 2308, Australia
| | - Gunther G Andersson
- Flinders Institute of Nanoscale Science and Technology, College of Science and Engineering, Flinders University, Adelaide, South Australia 5042, Australia
| | - Vincent S J Craig
- Department of Material Physics, Research School of Physics, Australian National University, Canberra, Australian Capital Territory 0200, Australia
| | - Alister J Page
- Discipline of Chemistry, The University of Newcastle, Callaghan, New South Wales 2308, Australia
| |
Collapse
|
10
|
Song C, Wang LP. A Polarizable QM/MM Model That Combines the State-Averaged CASSCF and AMOEBA Force Field for Photoreactions in Proteins. J Chem Theory Comput 2024. [PMID: 39088696 DOI: 10.1021/acs.jctc.4c00623] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/03/2024]
Abstract
This study presents the polarizable quantum mechanics/molecular mechanics (QM/MM) embedding of the state-averaged complete active space self-consistent field (SA-CASSCF) in the atomic multipole optimized energetics for biomolecular applications (AMOEBA) force field for the purpose of studying photoreactions in protein environments. We describe two extensions of our previous work that combine SA-CASSCF with AMOEBA water models, allowing it to be generalized to AMOEBA models for proteins and other macromolecules. First, we discuss how our QM/MM model accounts for the discrepancy between the direct and polarization electric fields that arises in the AMOEBA description of intramolecular polarization. A second improvement is the incorporation of link atom schemes to treat instances in which the QM/MM boundary goes through covalent bonds. A single-link atom scheme and double-link atom scheme are considered in this work, and we will discuss how electrostatic interaction, van der Waals interaction, and various kinds of valence terms are treated across the boundary. To test the accuracy of the link atom scheme, we will compare QM/MM with full QM calculations and study how the errors in ground state properties, excited state properties, and excitation energies change when tuning the parameters in the link atom scheme. We will also test the new SA-CASSCF/AMOEBA method on an elementary reaction step in NanoLuc, an artificial bioluminescence luciferase. We will show how the reaction mechanism is different when calculated in the gas phase, in polarizable continuum medium (PCM), versus in protein AMOEBA models.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Chemistry, University of California, Davis, 1 Shields Avenue, Davis, California 95616, United States
| | - Lee-Ping Wang
- Department of Chemistry, University of California, Davis, 1 Shields Avenue, Davis, California 95616, United States
| |
Collapse
|
11
|
Plé T, Adjoua O, Lagardère L, Piquemal JP. FeNNol: An efficient and flexible library for building force-field-enhanced neural network potentials. J Chem Phys 2024; 161:042502. [PMID: 39051830 DOI: 10.1063/5.0217688] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2024] [Accepted: 06/28/2024] [Indexed: 07/27/2024] Open
Abstract
Neural network interatomic potentials (NNPs) have recently proven to be powerful tools to accurately model complex molecular systems while bypassing the high numerical cost of ab initio molecular dynamics simulations. In recent years, numerous advances in model architectures as well as the development of hybrid models combining machine-learning (ML) with more traditional, physically motivated, force-field interactions have considerably increased the design space of ML potentials. In this paper, we present FeNNol, a new library for building, training, and running force-field-enhanced neural network potentials. It provides a flexible and modular system for building hybrid models, allowing us to easily combine state-of-the-art embeddings with ML-parameterized physical interaction terms without the need for explicit programming. Furthermore, FeNNol leverages the automatic differentiation and just-in-time compilation features of the Jax Python library to enable fast evaluation of NNPs, shrinking the performance gap between ML potentials and standard force-fields. This is demonstrated with the popular ANI-2x model reaching simulation speeds nearly on par with the AMOEBA polarizable force-field on commodity GPUs (graphics processing units). We hope that FeNNol will facilitate the development and application of new hybrid NNP architectures for a wide range of molecular simulation problems.
Collapse
Affiliation(s)
- Thomas Plé
- Sorbonne Université, LCT, UMR 7616 CNRS, 75005 Paris, France
| | - Olivier Adjoua
- Sorbonne Université, LCT, UMR 7616 CNRS, 75005 Paris, France
| | - Louis Lagardère
- Sorbonne Université, LCT, UMR 7616 CNRS, 75005 Paris, France
| | | |
Collapse
|
12
|
Sedano LF, Blazquez S, Vega C. Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant? J Chem Phys 2024; 161:044505. [PMID: 39046346 DOI: 10.1063/5.0211871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2024] [Accepted: 06/30/2024] [Indexed: 07/25/2024] Open
Abstract
The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.
Collapse
Affiliation(s)
- L F Sedano
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - S Blazquez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
13
|
Gogal RA, Nessler AJ, Thiel AC, Bernabe HV, Corrigan Grove RA, Cousineau LM, Litman JM, Miller JM, Qi G, Speranza MJ, Tollefson MR, Fenn TD, Michaelson JJ, Okada O, Piquemal JP, Ponder JW, Shen J, Smith RJH, Yang W, Ren P, Schnieders MJ. Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment. J Chem Phys 2024; 161:012501. [PMID: 38958156 PMCID: PMC11223778 DOI: 10.1063/5.0214652] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2024] [Accepted: 06/17/2024] [Indexed: 07/04/2024] Open
Abstract
Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid-base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.
Collapse
Affiliation(s)
- Rose A. Gogal
- Roy J. Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, USA
| | - Aaron J. Nessler
- Roy J. Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, USA
| | - Andrew C. Thiel
- Roy J. Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, USA
| | - Hernan V. Bernabe
- Roy J. Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, USA
| | - Rae A. Corrigan Grove
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Leah M. Cousineau
- Department of Biochemistry and Molecular Biology, University of Iowa, Iowa City, Iowa 52242, USA
| | - Jacob M. Litman
- Department of Biochemistry and Molecular Biology, University of Iowa, Iowa City, Iowa 52242, USA
| | - Jacob M. Miller
- Roy J. Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, USA
| | - Guowei Qi
- Department of Biochemistry and Molecular Biology, University of Iowa, Iowa City, Iowa 52242, USA
| | - Matthew J. Speranza
- Roy J. Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, USA
| | - Mallory R. Tollefson
- Roy J. Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, USA
| | - Timothy D. Fenn
- Analytical Development, LEXEO Therapeutics, New York, New York 10010, USA
| | - Jacob J. Michaelson
- Department of Psychiatry, University of Iowa Hospitals and Clinics, Iowa City, Iowa 52242, USA
| | - Okimasa Okada
- Sohyaku Innovative Research Division, Mitsubishi Tanabe Pharma Corporation, 1000 Kamoshida-cho, Aoba-ku, Yokohama, Kanagawa 227-0033, Japan
| | | | - Jay W. Ponder
- Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri 63130, USA
| | - Jana Shen
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, USA
| | - Richard J. H. Smith
- Molecular Otolaryngology and Renal Research Laboratories, Department of Otolaryngology, University of Iowa Hospitals and Clinics, Iowa City, Iowa 52242, USA
| | | | - Pengyu Ren
- Department of Biomedical Engineering, University of Texas, Austin, Texas 78712, USA
| | | |
Collapse
|
14
|
Kirov MV. Edge-sharing water prisms. Phys Chem Chem Phys 2024; 26:17777-17784. [PMID: 38873970 DOI: 10.1039/d4cp00745j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2024]
Abstract
Regular-shaped water clusters and nanostructures can be of particular interest for the self-assembly of complex structures and functional materials involving water molecules. Polyhedral water clusters are the most well studied. The low-energy structures of water hexamer, octamer and decamer are shaped like right prisms. The cavities of gas hydrates also have a polyhedral shape. Ice nanotubes are of no less interest both as configurations of global energy minima and in terms of possible applications. This article presents the results of the first systematic study of the structure and properties of a special class of water clusters. It reveals different combinations of three right prisms sharing an edge. According to modern calculations, the configurations of global energy minima of water clusters with the number of molecules being eighteen and twenty have this structure precisely. The topological characteristics of the edge-sharing water prisms are studied and the formulas for estimating the residual entropy are obtained. For these clusters (3-prisms), the usefulness of using a discrete model of intermolecular interaction [strong-weak-effective-bond model], previously developed for polyhedral water clusters, is shown. Calculations of the binding energy were performed using the non-additive Amoeba potential. To understand the relative stability of 3-prisms among other cluster forms, we use the structures recently reported in a published database containing over 3 × 106 unique water cluster networks (H2O)N of size N = 3-25 [Rakshit et al., J. Chem. Phys., 2019, 151(21), 214307, DOI: 10.1063/1.5128378].
Collapse
Affiliation(s)
- Mikhail V Kirov
- Earth Cryosphere Institute, Tyumen Scientific Centre SB RAS, Tyumen 625026, Malygina 86, Russian Federation.
| |
Collapse
|
15
|
Nikitin A, Wang F. Simulation of Linear and Cyclic Alkanes with Second-Order Møller-Plesset Perturbation Theory through Adaptive Force Matching. J Chem Theory Comput 2024; 20:5241-5249. [PMID: 38848512 PMCID: PMC11209940 DOI: 10.1021/acs.jctc.4c00509] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2024] [Revised: 05/28/2024] [Accepted: 05/29/2024] [Indexed: 06/09/2024]
Abstract
Predicting ensemble properties, such as density and heat of vaporization, of small hydrocarbons is challenging due to the dispersion-dominated weak interactions between these molecules. With the adaptive force matching (AFM) method, the bonded and short-range nonbonded interactions are fitted to second-order Møller-Plesset perturbation theory (MP2) references computed with the def2-TZVP basis set. The dispersion is modeled using symmetry adapted perturbation theory (SAPT) at MP4 accuracy using the def2-TZVPD basis set. A new charge matrix decomposition technique is described to obtain partial charges in AFM. Although the models developed do not have any empirical parameters, several properties of the resulting models are compared with experiments as validations. The density, heat of vaporization, pressure dependence of density, diffusion constants, and surface tensions all show quantitative agreement with experiments. Although the density shows a very small systematic error, which could be due to missing three-body dispersion, the heat of vaporization agrees with experiments of within 0.5%. The paper shows that AFM can be used as a reliable tool to enable simulations at post-Hartree-Fock quality at the cost of molecular mechanics force fields.
Collapse
Affiliation(s)
- Alexei Nikitin
- Department of Chemistry and
Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, United States
| | - Feng Wang
- Department of Chemistry and
Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, United States
| |
Collapse
|
16
|
Lagardère L, Maurin L, Adjoua O, El Hage K, Monmarché P, Piquemal JP, Hénin J. Lambda-ABF: Simplified, Portable, Accurate, and Cost-Effective Alchemical Free-Energy Computation. J Chem Theory Comput 2024; 20:4481-4498. [PMID: 38805379 DOI: 10.1021/acs.jctc.3c01249] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2024]
Abstract
We introduce the lambda-Adaptive Biasing Force (lambda-ABF) method for the computation of alchemical free-energy differences. We propose a software implementation and showcase it on biomolecular systems. The method arises from coupling multiple-walker adaptive biasing force with λ-dynamics. The sampling of the alchemical variable is continuous and converges toward a uniform distribution, making manual optimization of the λ schedule unnecessary. Contrary to most other approaches, alchemical free-energy estimates are obtained immediately without any postprocessing. Free diffusion of λ improves orthogonal relaxation compared to fixed-λ thermodynamic integration or free-energy perturbation. Furthermore, multiple walkers provide generic orthogonal space coverage with minimal user input and negligible computational overhead. We show that our high-performance implementations coupling the Colvars library with NAMD and Tinker-HP can address real-world cases including ligand-receptor binding with both fixed-charge and polarizable models, with a demonstrably richer sampling than fixed-λ methods. The implementation is fully open-source, publicly available, and readily usable by practitioners of current alchemical methods. Thanks to the portable Colvars library, lambda-ABF presents a unified user interface regardless of the back-end (NAMD, Tinker-HP, or any software to be interfaced in the future), sparing users the effort of learning multiple interfaces. Finally, the Colvars Dashboard extension of the visual molecular dynamics (VMD) software provides an interactive monitoring and diagnostic tool for lambda-ABF simulations.
Collapse
Affiliation(s)
- Louis Lagardère
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
- Sorbonne Université, Institut Parisien de Chimie Physique et Théorique, FR2622 CNRS, 75005 Paris, France
- Qubit Pharmaceuticals, 29 rue du Faubourg Saint Jacques, 75014 Paris, France
| | - Lise Maurin
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
- Sorbonne Université, Laboratoire Jacques-Louis Lions, UMR 7589 CNRS, 75005 Paris, France
| | - Olivier Adjoua
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
| | - Krystel El Hage
- Qubit Pharmaceuticals, 29 rue du Faubourg Saint Jacques, 75014 Paris, France
| | - Pierre Monmarché
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
- Sorbonne Université, Laboratoire Jacques-Louis Lions, UMR 7589 CNRS, 75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
- Qubit Pharmaceuticals, 29 rue du Faubourg Saint Jacques, 75014 Paris, France
| | - Jérôme Hénin
- Laboratoire de Biochimie Théorique, Université Paris Cité, CNRS, UPR 9080, 75005 Paris, France
| |
Collapse
|
17
|
Yang X, Liu C, Ren P. Exploring Biomolecular Conformational Dynamics with Polarizable Force Field AMOEBA and Enhanced Sampling Method Milestoning. J Chem Theory Comput 2024; 20:4065-4075. [PMID: 38742922 PMCID: PMC11187603 DOI: 10.1021/acs.jctc.4c00053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2024]
Abstract
Conformational dynamics play a crucial role in determining the behavior of the biomolecules. Polarizable force fields, such as AMOEBA, can accurately capture electrostatic interactions underlying the conformational space. However, applying a polarizable force field in molecular dynamics (MD) simulations can be computationally expensive, especially in studying long-time-scale dynamics. To overcome this challenge, we incorporated the AMOEBA potential with Milestoning, an enhanced sampling method in this work. This integration allows us to efficiently sample the rare and important conformational states of a biomolecule by using many short and independent molecular dynamics trajectories with the AMOEBA force field. We applied this method to investigate the conformational dynamics of alanine dipeptide, DNA, and RNA A-B form conversion. Well-converged thermodynamic and kinetic properties were obtained, including the free energy difference, mean first passage time, and critical transitions between states. Our results demonstrate the power of integrating polarizable force fields with enhanced sampling methods in quantifying the thermodynamic and kinetic properties of biomolecules at the atomic level.
Collapse
Affiliation(s)
- Xudong Yang
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712, USA
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712, USA
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712, USA
| |
Collapse
|
18
|
Dodia M, Rouxel JR, Cho D, Zhang Y, Keefer D, Bonn M, Nagata Y, Mukamel S. Water Solvent Reorganization upon Ultrafast Resonant Stimulated X-ray Raman Excitation of a Metalloporphyrin Dimer. J Chem Theory Comput 2024; 20:4254-4264. [PMID: 38727197 DOI: 10.1021/acs.jctc.4c00040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2024]
Abstract
We propose an X-ray Raman pump-X-ray diffraction probe scheme to follow solvation dynamics upon charge migration in a solute molecule. The X-ray Raman pump selectively prepares a valence electronic wavepacket in the solute, while the probe provides information about the entire molecular ensemble. A combination of molecular dynamics and ab initio quantum chemistry simulations is applied to a Zn-Ni porphyrin dimer in water. Using time-resolved X-ray diffraction and pair distribution functions, we extracted solvation shell dynamics.
Collapse
Affiliation(s)
- Mayank Dodia
- Max Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | - Jérémy R Rouxel
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Daeheum Cho
- Department of Chemistry and Green-Nano Materials Research Center, Kyungpook National University, Daegu 41566, Republic of Korea
| | - Yu Zhang
- Ames National Laboratory, Iowa State University, Ames, Iowa 50011, United States
| | - Daniel Keefer
- Max Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | - Mischa Bonn
- Max Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | - Yuki Nagata
- Max Planck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | - Shaul Mukamel
- Department of Chemistry, University of California, Irvine, California 92697-2025, United States
| |
Collapse
|
19
|
Chatterjee S, Nochebuena J, Cisneros GA. Impact of an Ionic Liquid Solution on Horseradish Peroxidase Activity. J Am Chem Soc 2024; 146:13247-13257. [PMID: 38701006 DOI: 10.1021/jacs.4c01100] [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: 05/05/2024]
Abstract
Horseradish peroxidase (HRP) is an enzyme that oxidizes pollutants from wastewater. A previous report indicated that peroxidases can have an enhancement in initial enzymatic activity in an aqueous solution of 0.26 M 1-ethyl-3-methylimidazolium ethyl sulfate ([EMIm][EtSO4]) at neutral pH. However, the atomistic details remain elusive. In the enzymatic landscape of HRP, compound II (Cpd II) plays a key role and involves a histidine (H42) residue. Cpd II exists as oxoferryl (2a) or hydroxoferryl (2b(FeIV)) forms, where 2a is the predominantly observed form in experimental studies. Intriguingly, the ferric 2b(FeIII) form seen in synthetic complexes has not been observed in HRP. Here, we have investigated the structure and dynamics of HRP in pure water and aqueous [EMIm][EtSO4] (0.26 M), as well as the reaction mechanism of 2a to 2b conversion using polarizable molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations. When HRP is solvated in aq [EMIm][EtSO4], the catalytic water displaces, and H42 directly orients over the ferryl moiety, allowing a direct proton transfer (PT) with a significant energy barrier reduction. Conversely, in neat water, the reaction of 2a to 2b follows the previously reported mechanism. We further investigated the deprotonated form of H42. Analysis of the electric fields at the active site indicates that the aq [EMIm][EtSO4] medium facilitates the reaction by providing a more favorable environment compared with the system solvated in neat water. Overall, the atomic level supports the previous experimental observations and underscores the importance of favorable electric fields in the active site to promote catalysis.
Collapse
Affiliation(s)
- Shubham Chatterjee
- Department of Chemistry and Biochemistry, University of Texas at Dallas, Richardson, Texas 75080, United States
| | - Jorge Nochebuena
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75080, United States
| | - G Andrés Cisneros
- Department of Chemistry and Biochemistry, University of Texas at Dallas, Richardson, Texas 75080, United States
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75080, United States
| |
Collapse
|
20
|
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
|
21
|
Heindel JP, LaCour RA, Head-Gordon T. The role of charge in microdroplet redox chemistry. Nat Commun 2024; 15:3670. [PMID: 38693110 DOI: 10.1038/s41467-024-47879-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2023] [Accepted: 04/12/2024] [Indexed: 05/03/2024] Open
Abstract
In charged water microdroplets, which occur in nature or in the lab upon ultrasonication or in electrospray processes, the thermodynamics for reactive chemistry can be dramatically altered relative to the bulk phase. Here, we provide a theoretical basis for the observation of accelerated chemistry by simulating water droplets of increasing charge imbalance to create redox agents such as hydroxyl and hydrogen radicals and solvated electrons. We compute the hydration enthalpy of OH- and H+ that controls the electron transfer process, and the corresponding changes in vertical ionization energy and vertical electron affinity of the ions, to create OH• and H• reactive species. We find that at ~ 20 - 50% of the Rayleigh limit of droplet charge the hydration enthalpy of both OH- and H+ have decreased by >50 kcal/mol such that electron transfer becomes thermodynamically favorable, in correspondence with the more favorable vertical electron affinity of H+ and the lowered vertical ionization energy of OH-. We provide scaling arguments that show that the nanoscale calculations and conclusions extend to the experimental microdroplet length scale. The relevance of the droplet charge for chemical reactivity is illustrated for the formation of H2O2, and has clear implications for other redox reactions observed to occur with enhanced rates in microdroplets.
Collapse
Affiliation(s)
- Joseph P Heindel
- Kenneth S. Pitzer Theory Center and Department of Chemistry, Berkeley, CA, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - R Allen LaCour
- Kenneth S. Pitzer Theory Center and Department of Chemistry, Berkeley, CA, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Theory Center and Department of Chemistry, Berkeley, CA, USA.
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA.
- Departments of Bioengineering and Chemical and Biomolecular Engineering University of CAlifornia, Berkeley, CA, USA.
| |
Collapse
|
22
|
Muthachikavil AV, Sun G, Peng B, Tanaka H, Kontogeorgis GM, Liang X. Unraveling thermodynamic anomalies of water: A molecular simulation approach to probe the two-state theory with atomistic and coarse-grained water models. J Chem Phys 2024; 160:154505. [PMID: 38624123 DOI: 10.1063/5.0194036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2023] [Accepted: 03/19/2024] [Indexed: 04/17/2024] Open
Abstract
Thermodynamic and dynamic anomalies of water play a crucial role in supporting life on our planet. The two-state theory attributes these anomalies to a dynamic equilibrium between locally favored tetrahedral structures (LFTSs) and disordered normal liquid structures. This theory provides a straightforward, phenomenological explanation for water's unique thermodynamic and dynamic characteristics. To validate this two-state feature, it is critical to unequivocally identify these structural motifs in a dynamically fluctuating disordered liquid. In this study, we employ a recently introduced structural parameter (θavg) that characterizes the local angular order within the first coordination shell to identify these LFTSs through molecular dynamics simulations. We employ both realistic water models with a liquid-liquid critical point (LLCP) and a coarse-grained water model without an LLCP to study water's anomalies in low-pressure regions below 2 kbar. The two-state theory consistently describes water's thermodynamic anomalies in these models, both with and without an LLCP. This suggests that the anomalies predominantly result from the two-state features rather than criticality, particularly within experimentally accessible temperature-pressure regions.
Collapse
Affiliation(s)
- Aswin V Muthachikavil
- Department of Chemical and Biochemical Engineering, Center for Energy Resources Engineering, Technical University of Denmark, Building 229, Lyngby DK-2800, Denmark
| | - Gang Sun
- Department of Physics, Center for Advanced Quantum Studies, Beijing Normal University, Beijing 100875, China
| | - Baoliang Peng
- Research Institute of Petroleum Exploration & Development (RIPED), PetroChina, Beijing 100083, China
| | - Hajime Tanaka
- Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8904, Japan
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Georgios M Kontogeorgis
- Department of Chemical and Biochemical Engineering, Center for Energy Resources Engineering, Technical University of Denmark, Building 229, Lyngby DK-2800, Denmark
| | - Xiaodong Liang
- Department of Chemical and Biochemical Engineering, Center for Energy Resources Engineering, Technical University of Denmark, Building 229, Lyngby DK-2800, Denmark
| |
Collapse
|
23
|
McDonogh DP, Gale JD, Raiteri P, Gebauer D. Redefined ion association constants have consequences for calcium phosphate nucleation and biomineralization. Nat Commun 2024; 15:3359. [PMID: 38637527 PMCID: PMC11026415 DOI: 10.1038/s41467-024-47721-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2023] [Accepted: 04/10/2024] [Indexed: 04/20/2024] Open
Abstract
Calcium orthophosphates (CaPs), as hydroxyapatite (HAP) in bones and teeth are the most important biomineral for humankind. While clusters in CaP nucleation have long been known, their speciation and mechanistic pathways to HAP remain debated. Evidently, mineral nucleation begins with two ions interacting in solution, fundamentally underlying solute clustering. Here, we explore CaP ion association using potentiometric methods and computer simulations. Our results agree with literature association constants for Ca2+ and H2PO4-, and Ca2+ and HPO42-, but not for Ca2+ and PO43- ions, which previously has been strongly overestimated by two orders of magnitude. Our data suggests that the discrepancy is due to a subtle, premature phase separation that can occur at low ion activity products, especially at higher pH. We provide an important revision of long used literature constants, where association of Ca2+ and PO43- actually becomes negligible below pH 9.0, in contrast to previous values. Instead, [CaHPO4]0 dominates the aqueous CaP speciation between pH ~6-10. Consequently, calcium hydrogen phosphate association is critical in cluster-based precipitation in the near-neutral pH regime, e.g., in biomineralization. The revised thermodynamics reveal significant and thus far unexplored multi-anion association in computer simulations, constituting a kinetic trap that further complicates aqueous calcium phosphate speciation.
Collapse
Affiliation(s)
- David P McDonogh
- Institute of Inorganic Chemistry, Leibniz University Hannover, Callinstr. 9, 30167, Hannover, Germany
| | - Julian D Gale
- Curtin Institute for Computation and School of Molecular and Life Sciences, Curtin University, P.O. Box U1987, Perth, WA, 6845, Australia
| | - Paolo Raiteri
- Curtin Institute for Computation and School of Molecular and Life Sciences, Curtin University, P.O. Box U1987, Perth, WA, 6845, Australia.
| | - Denis Gebauer
- Institute of Inorganic Chemistry, Leibniz University Hannover, Callinstr. 9, 30167, Hannover, Germany.
| |
Collapse
|
24
|
Lambros E, Fetherolf JH, Hammes-Schiffer S, Li X. A Many-Body Perspective of Nuclear Quantum Effects in Aqueous Clusters. J Phys Chem Lett 2024; 15:4070-4075. [PMID: 38587257 DOI: 10.1021/acs.jpclett.4c00439] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/09/2024]
Abstract
Nuclear quantum effects play an important role in the structure and thermodynamics of aqueous systems. By performing a many-body expansion with nuclear-electronic orbital (NEO) theory, we show that proton quantization can give rise to significant energetic contributions for many-body interactions spanning several molecules in single-point energy calculations of water clusters. Although zero-point motion produces a large increase in energy at the one-body level, nuclear quantum effects serve to stabilize higher-order molecular interactions. These results are significant because they demonstrate that nuclear quantum effects play a nontrivial role in many-body interactions of aqueous systems. Our approach also provides a pathway for incorporating nuclear quantum effects into water potential energy surfaces. The NEO approach is advantageous for many-body expansion analyses because it includes nuclear quantum effects directly in the energies.
Collapse
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Jonathan H Fetherolf
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
25
|
Thiel A, Speranza MJ, Jadhav S, Stevens LL, Unruh DK, Ren P, Ponder JW, Shen J, Schnieders MJ. Constant-pH Simulations with the Polarizable Atomic Multipole AMOEBA Force Field. J Chem Theory Comput 2024; 20:2921-2933. [PMID: 38507252 PMCID: PMC11008096 DOI: 10.1021/acs.jctc.3c01180] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2023] [Revised: 03/05/2024] [Accepted: 03/05/2024] [Indexed: 03/22/2024]
Abstract
Accurately predicting protein behavior across diverse pH environments remains a significant challenge in biomolecular simulations. Existing constant-pH molecular dynamics (CpHMD) algorithms are limited to fixed-charge force fields, hindering their application to biomolecular systems described by permanent atomic multipoles or induced dipoles. This work overcomes these limitations by introducing the first polarizable CpHMD algorithm in the context of the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field. Additionally, our implementation in the open-source Force Field X (FFX) software has the unique ability to handle titration state changes for crystalline systems including flexible support for all 230 space groups. The evaluation of constant-pH molecular dynamics (CpHMD) with the AMOEBA force field was performed on 11 crystalline peptide systems that span the titrating amino acids (Asp, Glu, His, Lys, and Cys). Titration states were correctly predicted for 15 out of the 16 amino acids present in the 11 systems, including for the coordination of Zn2+ by cysteines. The lone exception was for a HIS-ALA peptide where CpHMD predicted both neutral histidine tautomers to be equally populated, whereas the experimental model did not consider multiple conformers and diffraction data are unavailable for rerefinement. This work demonstrates the promise polarizable CpHMD simulations for pKa predictions, the study of biochemical mechanisms such as the catalytic triad of proteases, and for improved protein-ligand binding affinity accuracy in the context of pharmaceutical lead optimization.
Collapse
Affiliation(s)
- Andrew
C. Thiel
- Department
of Biomedical Engineering, University of
Iowa, Iowa City, Iowa 52242, United States
| | - Matthew J. Speranza
- Department
of Biomedical Engineering, University of
Iowa, Iowa City, Iowa 52242, United States
| | - Sanika Jadhav
- Department
of Pharmaceutical Sciences and Experimental Therapeutics, University of Iowa, Iowa City, Iowa 52242, United States
| | - Lewis L. Stevens
- Department
of Pharmaceutical Sciences and Experimental Therapeutics, University of Iowa, Iowa City, Iowa 52242, United States
| | - Daniel K. Unruh
- Office
of the Vice President for Research, University
of Iowa, Iowa City, Iowa 52242, United
States
| | - Pengyu Ren
- Department
of Biomedical Engineering, University of
Texas, Austin, Texas 78712, United States
| | - Jay W. Ponder
- Department
of Chemistry, Washington University in St.
Louis, St. Louis, Missouri 63130, United
States
| | - Jana Shen
- Department
of Pharmaceutical Sciences, University of
Maryland School of Pharmacy, Baltimore, Maryland 21201, United States
| | - Michael J. Schnieders
- Department
of Biomedical Engineering, University of
Iowa, Iowa City, Iowa 52242, United States
- Department
of Biochemistry, University of Iowa, Iowa City, Iowa 52242, United States
| |
Collapse
|
26
|
Bondanza M, Nottoli T, Nottoli M, Cupellini L, Lipparini F, Mennucci B. The OpenMMPol library for polarizable QM/MM calculations of properties and dynamics. J Chem Phys 2024; 160:134106. [PMID: 38557842 DOI: 10.1063/5.0198251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Accepted: 03/13/2024] [Indexed: 04/04/2024] Open
Abstract
We present a new library designed to provide a simple and straightforward way to implement QM/AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications) and other polarizable QM/MM (Molecular Mechanics) methods based on induced point dipoles. The library, herein referred to as OpenMMPol, is free and open-sourced and is engineered to address the increasing demand for accurate and efficient QM/MM simulations. OpenMMPol is specifically designed to allow polarizable QM/MM calculations of ground state energies and gradients and excitation properties. Key features of OpenMMPol include a modular architecture facilitating extensibility, parallel computing capabilities for enhanced performance on modern cluster architectures, a user-friendly interface for intuitive implementation, and a simple and flexible structure for providing input data. To show the capabilities offered by the library, we present an interface with PySCF to perform QM/AMOEBA molecular dynamics, geometry optimization, and excited-state calculation based on (time-dependent) density functional theory.
Collapse
Affiliation(s)
- Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Tommaso Nottoli
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Michele Nottoli
- Institute of Applied Analysis and Numerical Simulation, Universität Stuttgart, Pfaffenwaldring 57, D-70569 Stuttgart, Germany
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
27
|
Nochebuena J, Liu S, Cisneros GA. Relative cooperativity in neutral and charged molecular clusters using QM/MM calculations. J Chem Phys 2024; 160:134301. [PMID: 38557841 DOI: 10.1063/5.0203020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2024] [Accepted: 03/13/2024] [Indexed: 04/04/2024] Open
Abstract
QM/MM methods have been used to study electronic structure properties and chemical reactivity in complex molecular systems where direct electronic structure calculations are not feasible. In our previous work, we showed that non-polarizable force fields, by design, describe intermolecular interactions through pairwise interactions, overlooking many-body interactions involving three or more particles. In contrast, polarizable force fields account partially for many-body effects through polarization, but still handle van der Waals and permanent electrostatic interactions pairwise. We showed that despite those limitations, polarizable and non-polarizable force fields can reproduce relative cooperativity achieved using density functional theory due to error compensation mechanisms. In this contribution, we assess the performance of QM/MM methods in reproducing these phenomena. Our study highlights the significance of the QM region size and force field choice in QM/MM calculations, emphasizing the importance of parameter validation to obtain accurate interaction energy predictions.
Collapse
Affiliation(s)
- Jorge Nochebuena
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75080, USA
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599, USA
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| | - G Andrés Cisneros
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75080, USA
- Department of Chemistry and Biochemistry, University of Texas at Dallas, Richardson, Texas 75080, USA
| |
Collapse
|
28
|
Wang Y, Inizan TJ, Liu C, Piquemal JP, Ren P. Incorporating Neural Networks into the AMOEBA Polarizable Force Field. J Phys Chem B 2024; 128:2381-2388. [PMID: 38445577 PMCID: PMC10985787 DOI: 10.1021/acs.jpcb.3c08166] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
Neural network potentials (NNPs) offer significant promise to bridge the gap between the accuracy of quantum mechanics and the efficiency of molecular mechanics in molecular simulation. Most NNPs rely on the locality assumption that ensures the model's transferability and scalability and thus lack the treatment of long-range interactions, which are essential for molecular systems in the condensed phase. Here we present an integrated hybrid model, AMOEBA+NN, which combines the AMOEBA potential for the short- and long-range noncovalent atomic interactions and an NNP to capture the remaining local covalent contributions. The AMOEBA+NN model was trained on the conformational energy of the ANI-1x data set and tested on several external data sets ranging from small molecules to tetrapeptides. The hybrid model demonstrated substantial improvements over the baseline models in term of accuracy as the molecule size increased, suggesting its potential as a next-generation approach for chemically accurate molecular simulations.
Collapse
Affiliation(s)
- Yanxing Wang
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Théo Jaffrelot Inizan
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Jean-Philip Piquemal
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, Paris 75005, France
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
29
|
Quoika PK, Zacharias M. Liquid-Vapor Coexistence and Spontaneous Evaporation at Atmospheric Pressure of Common Rigid Three-Point Water Models in Molecular Simulations. J Phys Chem B 2024; 128:2457-2468. [PMID: 38427971 PMCID: PMC10945489 DOI: 10.1021/acs.jpcb.3c08183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2023] [Revised: 02/21/2024] [Accepted: 02/22/2024] [Indexed: 03/03/2024]
Abstract
Molecular dynamics (MD) simulations are widely used to investigate molecular systems at atomic resolution including biomolecular structures, drug-receptor interactions, and novel materials. Frequently, MD simulations are performed in an aqueous solution with explicit models of water molecules. Commonly, such models are parameterized to reproduce the liquid phase of water under ambient conditions. However, often, simulations at significantly higher temperatures are also of interest. Hence, it is important to investigate the equilibrium of the liquid and vapor phases of molecular models of water at elevated temperatures. Here, we evaluate the behavior of 11 common rigid three-point water models over a wide range of temperatures. From liquid-vapor coexistence simulations, we estimated the critical points and studied the spontaneous evaporation of these water models. Moreover, we investigated the influence of the system size, choice of the pressure-coupling algorithm, and rate of heating on the process and compared them with the experimental data. We found that modern rigid three-point water models reproduce the critical point surprisingly well. Furthermore, we discovered that the critical temperature correlates with the quadrupole moment of the respective water model. This indicates that the spatial arrangement of the partial charges is important for reproducing the liquid-vapor phase transition. Our findings may guide the selection of water models for simulations conducted at high temperatures.
Collapse
Affiliation(s)
- Patrick K. Quoika
- Center for Functional Protein
Assemblies, Technical University of Munich, Ernst-Otto-Fischer-Str. 8, Garching 85748, Germany
| | - Martin Zacharias
- Center for Functional Protein
Assemblies, Technical University of Munich, Ernst-Otto-Fischer-Str. 8, Garching 85748, Germany
| |
Collapse
|
30
|
Zhao S, Cieplak P, Duan Y, Luo R. Assessment of Amino Acid Electrostatic Parametrizations of the Polarizable Gaussian Multipole Model. J Chem Theory Comput 2024; 20:2098-2110. [PMID: 38394331 PMCID: PMC11060985 DOI: 10.1021/acs.jctc.3c01347] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2024]
Abstract
Accurate parametrization of amino acids is pivotal for the development of reliable force fields for molecular modeling of biomolecules such as proteins. This study aims to assess amino acid electrostatic parametrizations with the polarizable Gaussian Multipole (pGM) model by evaluating the performance of the pGM-perm (with atomic permanent dipoles) and pGM-ind (without atomic permanent dipoles) variants compared to the traditional RESP model. The 100-conf-combterm fitting strategy on tetrapeptides was adopted, in which (1) all peptide bond atoms (-CO-NH-) share identical set of parameters and (2) the total charges of the two terminal N-acetyl (ACE) and N-methylamide (NME) groups were set to neutral. The accuracy and transferability of electrostatic parameters across peptides with varying lengths and real-world examples were examined. The results demonstrate the enhanced performance of the pGM-perm model in accurately representing the electrostatic properties of amino acids. This insight underscores the potential of the pGM-perm model and the 100-conf-combterm strategy for the future development of the pGM force field.
Collapse
Affiliation(s)
- Shiji Zhao
- Nurix Therapeutics, Inc., 1700 Owens St. Suite 205, San Francisco, CA 94158, USA
| | - Piotr Cieplak
- SBP Medical Discovery Institute, 10901 North Torrey Pines Road, La Jolla, California 92037, United States
| | - Yong Duan
- UC Davis Genome Center and Department of Biomedical Engineering, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - Ray Luo
- Departments of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine. Irvine, California 92697, United States
| |
Collapse
|
31
|
Debnath T, Cisneros GA. Investigation of the stability of D5SIC-DNAM-incorporated DNA duplex in Taq polymerase binary system: a systematic classical MD approach. Phys Chem Chem Phys 2024; 26:7287-7295. [PMID: 38353000 PMCID: PMC11078294 DOI: 10.1039/d3cp05571j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/29/2024]
Abstract
DNA polymerases are fundamental enzymes that play a crucial role in processing DNA with high fidelity and accuracy ensuring the faithful transmission of genetic information. The recognition of unnatural base pairs (UBPs) by polymerases, enabling their replication, represents a significant and groundbreaking discovery with profound implications for genetic expansion. Romesberg et al. examined the impact of DNA containing 2,6-dimethyl-2H-isoquiniline-1-thione: D5SIC (DS) and 2-methoxy-3-methylnaphthalene: DNAM (DN) UBPs bound to T. aquaticus DNA polymerase (Taq) through crystal structure analysis. Here, we have used polarizable and nonpolarizable classical molecular dynamics (MD) simulations to investigate the structural aspects and stability of Taq in complex with a DNA duplex including a DS-DN pair in the terminal 3' and 5' positions. Our results suggest that the flexibility of UBP-incorporated DNA in the terminal position is arrested by the polymerase, thus preventing fraying and mispairing. Our investigation also reveals that the UBP remains in an intercalated conformation inside the active site, exhibiting two distinct orientations in agreement with experimental findings. Our analysis pinpoints particular residues responsible for favorable interactions with the UBP, with some relying on van der Waals interactions while other on Coulombic forces.
Collapse
Affiliation(s)
- Tanay Debnath
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75080, Dallas, USA.
| | - G Andrés Cisneros
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75080, Dallas, USA.
- Department of Chemistry and Biochemistry, University of Texas at Dallas, Richardson, Texas 75080, Dallas, USA
| |
Collapse
|
32
|
Debnath T, Cisneros GA. Investigation of dynamical flexibility of D5SIC-DNAM inside DNA duplex in aqueous solution: a systematic classical MD approach. Phys Chem Chem Phys 2024; 26:7435-7445. [PMID: 38353005 PMCID: PMC11080001 DOI: 10.1039/d3cp05572h] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/29/2024]
Abstract
Incorporation of artificial 3rd base pairs (unnatural base pairs, UBPs) has emerged as a fundamental technique in pursuit of expanding the genetic alphabet. 2,6-Dimethyl-2H-isoquiniline-1-thione: D5SIC (DS) and 2-methoxy-3-methylnaphthalene: DNAM (DN), a potential unnatural base pair (UBP) developed by Romesberg and colleagues, has been shown to have remarkable capability for replication within DNA. Crystal structures of a Taq polymerase/double-stranded DNA (ds-DNA) complex containing a DS-DN pair in the 3' terminus showed a parallelly stacked geometry for the pre-insertion, and an intercalated geometry for the post-insertion structure. Unconventional orientations of DS-DN inside a DNA duplex have inspired scientists to investigate the conformational orientations and structural properties of UBP-incorporated DNA. In recent years, computational simulations have been used to investigate the geometry of DS-DN within the DNA duplex; nevertheless, unresolved questions persist owing to inconclusive findings. In this work, we investigate the structural and dynamical properties of DS and DN inside a ds-DNA strand in aqueous solution considering both short and long DNA templates using polarizable, and non-polarizable classical MD simulations. Flexible conformational change of UBP with major populations of Watson-Crick-Franklin (WCF) and three distinct non-Watson-Crick-Franklin (nWCFP1, nWCFP2, nWCFO) conformations through intra and inter-strand flipping have been observed. Our results suggest that a dynamical conformational change leads to the production of diffierent conformational distribution for the systems. Simulations with a short ds-DNA duplex suggest nWCF (P1 and O) as the predominant structures, whereas long ds-DNA duplex simulations indicate almost equal populations of WCF, nWCFP1, nWCFO. DS-DN in the terminal position is found to be more flexible with occasional mispairing and fraying. Overall, these results suggest flexibility and dynamical conformational change of the UBP as well as indicate varied conformational distribution irrespective of starting orientation of the UBP and length og DNA strand.
Collapse
Affiliation(s)
- Tanay Debnath
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75080, Dallas, USA.
| | - G Andrés Cisneros
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75080, Dallas, USA.
- Department of Chemistry and Biochemistry, University of Texas at Dallas, Richardson, Texas 75080, Dallas, USA
| |
Collapse
|
33
|
Hu X, Amin KS, Schneider M, Lim C, Salahub D, Baldauf C. System-Specific Parameter Optimization for Nonpolarizable and Polarizable Force Fields. J Chem Theory Comput 2024; 20:1448-1464. [PMID: 38279917 PMCID: PMC10867808 DOI: 10.1021/acs.jctc.3c01141] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2023] [Revised: 12/04/2023] [Accepted: 12/05/2023] [Indexed: 01/29/2024]
Abstract
The accuracy of classical force fields (FFs) has been shown to be limited for the simulation of cation-protein systems despite their importance in understanding the processes of life. Improvements can result from optimizing the parameters of classical FFs or by extending the FF formulation by terms describing charge transfer (CT) and polarization (POL) effects. In this work, we introduce our implementation of the CTPOL model in OpenMM, which extends the classical additive FF formula by adding CT and POL. Furthermore, we present an open-source parametrization tool, called FFAFFURR, that enables the (system-specific) parametrization of OPLS-AA and CTPOL models. The performance of our workflow was evaluated by its ability to reproduce quantum chemistry energies and by molecular dynamics simulations of a zinc-finger protein.
Collapse
Affiliation(s)
- Xiaojuan Hu
- Fritz-Haber-Institut
der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Kazi S. Amin
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Markus Schneider
- Fritz-Haber-Institut
der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Carmay Lim
- Institute
of Biomedical Sciences, Academia Sinica, Taipei 115, Taiwan
- Department
of Chemistry, National Tsing Hua University, Hsinchu 300, Taiwan
| | - Dennis Salahub
- Centre
for Molecular Simulation and Department of Chemistry, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Carsten Baldauf
- Fritz-Haber-Institut
der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| |
Collapse
|
34
|
Delgado JM, Nagy PR, Varma S. Polarizable AMOEBA Model for Simulating Mg 2+·Protein·Nucleotide Complexes. J Chem Inf Model 2024; 64:378-392. [PMID: 38051630 PMCID: PMC11345861 DOI: 10.1021/acs.jcim.3c01513] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/07/2023]
Abstract
Molecular mechanics (MM) simulations have the potential to provide detailed insights into the mechanisms of enzymes that utilize nucleotides as cofactors. In most cases, the activities of these enzymes also require the binding of divalent cations to catalytic sites. However, modeling divalent cations in MM simulations has been challenging. The inclusion of explicit polarization was considered promising, but despite improvements over nonpolarizable force fields and despite the inclusion of "Nonbonded-fix (NB-fix)" corrections, errors in interaction energies of divalent cations with proteins remain large. Importantly, the application of these models fails to reproduce the experimental structural data on Mg2+·Protein·ATP complexes. Focusing on these complexes, here we provide a systematic assessment of the polarizable AMOEBA model and recommend critical changes that substantially improve its predictive performance. Our key results are as follows. We first show that our recent revision of the AMOEBA protein model (AMOEBABIO18-HFC), which contains high field corrections (HFCs) to induced dipoles, dramatically improves Mg2+-protein interaction energies, reducing the mean absolute error (MAE) from 17 to 10 kcal/mol. This further supports the general applicability of AMOEBABIO18-HFC. The inclusion of many-body NB-fix corrections further reduces MAE to 6 kcal/mol, which amounts to less than 2% error. The errors are estimated with respect to vdW-inclusive density functional theory that we benchmark against CCSD(T) calculations and experiments. We also present a new model of ATP with revised polarization parameters to better capture its high field response, as well as new vdW and dihedral parameters. The ATP model accurately predicts experimental Mg2+-ATP binding free energy in the aqueous phase and provides new insights into how Mg2+ associates with ATP. Finally, we show that molecular dynamics (MD) simulations of Mg2+·Kinase·ATP complexes carried out with these improvements lead to a better agreement in global and local catalytic site structures between MD and X-ray crystallography.
Collapse
Affiliation(s)
- Julian M Delgado
- Department of Molecular Biosciences, University of South Florida, 4202 E. Fowler Avenue, Tampa, Florida 33620, United States
| | - Péter R Nagy
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., Budapest H-1111, Hungary
- HUN-REN-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., Budapest H-1111, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., Budapest H-1111, Hungary
| | - Sameer Varma
- Department of Molecular Biosciences, University of South Florida, 4202 E. Fowler Avenue, Tampa, Florida 33620, United States
- Department of Physics, University of South Florida, 4202 E. Fowler Avenue, Tampa, Florida 33620, United States
| |
Collapse
|
35
|
Plé T, Lagardère L, Piquemal JP. Force-field-enhanced neural network interactions: from local equivariant embedding to atom-in-molecule properties and long-range effects. Chem Sci 2023; 14:12554-12569. [PMID: 38020379 PMCID: PMC10646944 DOI: 10.1039/d3sc02581k] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2023] [Accepted: 10/03/2023] [Indexed: 12/01/2023] Open
Abstract
We introduce FENNIX (Force-Field-Enhanced Neural Network InteraXions), a hybrid approach between machine-learning and force-fields. We leverage state-of-the-art equivariant neural networks to predict local energy contributions and multiple atom-in-molecule properties that are then used as geometry-dependent parameters for physically-motivated energy terms which account for long-range electrostatics and dispersion. Using high-accuracy ab initio data (small organic molecules/dimers), we trained a first version of the model. Exhibiting accurate gas-phase energy predictions, FENNIX is transferable to the condensed phase. It is able to produce stable Molecular Dynamics simulations, including nuclear quantum effects, for water predicting accurate liquid properties. The extrapolating power of the hybrid physically-driven machine learning FENNIX approach is exemplified by computing: (i) the solvated alanine dipeptide free energy landscape; (ii) the reactive dissociation of small molecules.
Collapse
Affiliation(s)
- Thomas Plé
- Sorbonne Université, LCT, UMR 7616 CNRS F-75005 Paris France thomas.ple@sorbonne-université louis.lagardere@sorbonne-université jean-philip.piquemal@sorbonne-université
| | - Louis Lagardère
- Sorbonne Université, LCT, UMR 7616 CNRS F-75005 Paris France thomas.ple@sorbonne-université louis.lagardere@sorbonne-université jean-philip.piquemal@sorbonne-université
| | - Jean-Philip Piquemal
- Sorbonne Université, LCT, UMR 7616 CNRS F-75005 Paris France thomas.ple@sorbonne-université louis.lagardere@sorbonne-université jean-philip.piquemal@sorbonne-université
| |
Collapse
|
36
|
Nochebuena J, Piquemal JP, Liu S, Cisneros GA. Cooperativity and Frustration Effects (or Lack Thereof) in Polarizable and Non-polarizable Force Fields. J Chem Theory Comput 2023; 19:7715-7730. [PMID: 37888874 DOI: 10.1021/acs.jctc.3c00762] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2023]
Abstract
Understanding cooperativity and frustration is crucial for studying biological processes such as molecular recognition and protein aggregation. Force fields have been extensively utilized to explore cooperativity in the formation of protein secondary structures and self-assembled systems. Multiple studies have demonstrated that polarizable force fields provide more accurate descriptions of this phenomenon compared to fixed-charge pairwise nonpolarizable force fields, thanks to the incorporation of polarization effects. In this study, we assess the performance of the AMOEBA polarizable force field and the AMBER and OPLS nonpolarizable pairwise force fields in capturing positive and negative cooperativity recently explored in neutral and charged molecular clusters using density functional theory. Our findings show that polarizable and nonpolarizable force fields qualitatively reproduce the relative cooperativity observed in electron structure calculations. However, AMBER and OPLS fail to describe absolute cooperativity. In contrast, AMOEBA accounts for the absolute cooperativity by considering interactions beyond pairwise interactions. According to the energy decomposition analysis, it is observed that the electrostatic interactions calculated with the AMBER and OPLS force fields seem to play an important and counterintuitive role in reproducing the adiabatic interaction energies calculated with density functional theory. However, it is important to note that these force fields, due to their nature, do not explicitly incorporate many-body effects, which limits their ability to accurately describe cooperativity. On the other hand, frustration in polarizable and nonpolarizable force fields is caused by changes in bond stretching and angle bending terms of the building blocks when they are forming a complex.
Collapse
Affiliation(s)
- Jorge Nochebuena
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75080, United States
| | - Jean-Philip Piquemal
- Laboratoire de Chimie théorique, Sorbonne Université, UMR 7616 CNRS, Paris 75005, France
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599, United States
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599, United States
| | - G Andrés Cisneros
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75080, United States
- Department of Chemistry and Biochemistry, University of Texas at Dallas, Richardson, Texas 75080, United States
| |
Collapse
|
37
|
Gudkovskikh SV, Kirov MV. Ice structures assembled from cubic water clusters of D 2d and S 4 symmetry. Acta Crystallogr A Found Adv 2023; 79:527-535. [PMID: 37706368 DOI: 10.1107/s2053273323007428] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Accepted: 08/23/2023] [Indexed: 09/15/2023] Open
Abstract
The study of self-assembly processes is of key importance for fundamental science and modern technologies. Cubic water clusters of D2d and S4 symmetry show great potential as building blocks for self-assembly. The objective of this paper is to construct possible ice structures formed by hydrogen bonding of these very stable water clusters. A number of such structures are herein presented, including quasi-2D and 3D ices as well as spatial layered and tubular ices. The energetics and structure of many configurations differing in the arrangement of hydrogen atoms in hydrogen bonds have been studied. It was established that the proton disorder of all such ices is of island type. The residual entropy of these ices is equal to ln(3)/4 in dimensionless form. For layered structures formed by the stacking of multiple bilayers, the determining role of the van der Waals interactions is shown. Note that, for all considered ices, the lowest-energy configurations are formed only by clusters of D2d symmetry.
Collapse
Affiliation(s)
- Sergey V Gudkovskikh
- Earth Cryosphere Institute, Tyumen Scientific Center Cryosphere SB RAS, Tyumen, 625000, Russian Federation
| | - Mikhail V Kirov
- Earth Cryosphere Institute, Tyumen Scientific Center Cryosphere SB RAS, Tyumen, 625000, Russian Federation
| |
Collapse
|
38
|
Ogrin P, Urbic T. Thermodynamics perturbation theory for solvation of nonpolar solutes in rose model. Phys Rev E 2023; 108:054135. [PMID: 38115497 DOI: 10.1103/physreve.108.054135] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2023] [Accepted: 11/07/2023] [Indexed: 12/21/2023]
Abstract
A simple model of water, called the rose model, is used in this work. The rose model is a very simple model that can provide insight into the anomalous properties of water. In the rose water model, the molecules are represented as two-dimensional Lennard-Jones disks with potentials for orientation-dependent pairwise interactions mimicking formations of hydrogen bonds. We have recently applied a Wertheim integral equation theory (IET) and a thermodynamic perturbation theory (TPT) to the rose model in bulk. These analytical theories offer the advantage of being computationally less intensive than computer simulations by orders of magnitudes. Here we have applied the TPT to study the transfer of a nonpolar solute into rose water, the so-called hydrophobic effect. Similarly as in our previous work for bulk water, we have found that the theory reproduces the computer simulation results quite well at higher temperatures, while the theories predict the qualitative trends at low temperatures.
Collapse
Affiliation(s)
- Peter Ogrin
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| | - Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
39
|
Muniz MC, Car R, Panagiotopoulos AZ. Neural Network Water Model Based on the MB-Pol Many-Body Potential. J Phys Chem B 2023; 127:9165-9171. [PMID: 37824703 DOI: 10.1021/acs.jpcb.3c04629] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2023]
Abstract
The MB-pol many-body potential accurately predicts many properties of water, including cluster, liquid phase, and vapor-liquid equilibrium properties, but its high computational cost can make applying it in large-scale simulations quite challenging. In order to address this limitation, we developed a "deep potential" neural network (DPMD) model based on the MB-pol potential for water. We find that a DPMD model trained on mostly liquid configurations yields a good description of the bulk liquid phase but severely underpredicts vapor-liquid coexistence densities. By contrast, adding cluster configurations to the neural network training set leads to a good agreement for the vapor coexistence densities. Liquid phase densities under supercooled conditions are also represented well, even though they were not included in the training set. These results confirm that neural network models can combine accuracy and transferability if sufficient attention is given to the construction of a representative training set for the target system.
Collapse
Affiliation(s)
- Maria Carolina Muniz
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Roberto Car
- Department of Chemistry, Department of Physics, Program in Applied and Computational Mathematics, and Princeton Materials Institute, Princeton University, Princeton, New Jersey 08544, United States
| | | |
Collapse
|
40
|
Herman KM, Stone AJ, Xantheas SS. Accurate Calculation of Many-Body Energies in Water Clusters Using a Classical Geometry-Dependent Induction Model. J Chem Theory Comput 2023; 19:6805-6815. [PMID: 37703063 DOI: 10.1021/acs.jctc.3c00575] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/14/2023]
Abstract
We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The moment expansion is carried out up to the hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA classical interaction potentials for water. The classical induction model with no adjustable parameters reproduces the ab initio 3-/4-body terms contained in these two databases with a root-mean-square error (RMSE) of 0.104/0.058 and a mean-absolute error (MAE) of 0.054/0.026 kcal/mol, respectively. These results are on par with the ones obtained by fitting the same data using over 14,000 (for the 3-body) and 200 (for the 4-body) parameters via Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy, included in the calculation of the 3-body energy, was found to be small but not quite negligible. The model represents a practical, efficient, and transferable approach for obtaining accurate nonadditive interactions for multicomponent systems without the need to perform tens of thousands of high level electronic structure calculations and fitting them with PIPs.
Collapse
Affiliation(s)
- Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, Washington 98185, United States
| | - Anthony J Stone
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, U.K
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98185, United States
- Advanced Computing Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MSIN J7-10, Richland, Washington 99352, United States
| |
Collapse
|
41
|
Ogrin P, Urbic T. Angle-dependent integral equation theory improves results of thermodynamics and structure of rose water model. J Chem Phys 2023; 159:114505. [PMID: 37732557 PMCID: PMC10908565 DOI: 10.1063/5.0159438] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2023] [Accepted: 08/18/2023] [Indexed: 09/22/2023] Open
Abstract
Orientation-dependent integral equation theory (ODIET) was applied to the rose water model. Structural and thermodynamic properties of water modeled with the rose model were calculated using ODIET and compared to results from orientation-averaged integral equation theory (IET) and Monte Carlo simulations. Rose water model is a simple two-dimensional water model where molecules of water are represented as Lennard-Jones disks with explicit hydrogen bonding potential in form of rose functions. Orientational dependency significantly improves IET, as the thermodynamic results obtained using ODIET are significantly more in agreement with results calculated using MC than in the case of the orientationally averaged version. At high temperatures, the agreement between the simulation and theory is quantitative; however, when temperatures lower, a slight deviation between results obtained with different methods appear. ODIET correctly predicts the radial distribution function; moreover, ODIet also enables the calculation of angular distributions. While the angular distributions obtained with ODIET are in qualitative agreement with distributions from MC simulations, the height of the peaks in angular distributions differs between methods. Using results from ODIET, the spatial distribution of water molecules was constructed, which aids in the interpretation of other structural properties. ODIET was also used to calculate fractions of molecules with different number of hydrogen bonds, which is in the agreement with the simulations. Overall, use of ODIET significantly improves the obtained results in comparison to standard IET.
Collapse
Affiliation(s)
- Peter Ogrin
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| | - Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
42
|
Sameera WMC, Jayaweera AP, Ishibashi A, Hidaka H, Oba Y, Watanabe N. A systematic mechanistic survey on the reactions between OH radical and CH 3OH on ice. Faraday Discuss 2023; 245:508-518. [PMID: 37335277 DOI: 10.1039/d3fd00033h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/21/2023]
Abstract
A systematic mechanistic survey was performed for the CH3OH + OH reaction on ice. ONIOM(ωB97X-D/Def2-TZVP:AMOEBA09) calculations suggested a range of binding energies for the CH2OH radical (0.29-0.69 eV) and CH3OH (0.15-0.72 eV) molecule on hexagonal water ice (Ih) and amorphous solid water (ASW). Computed average binding energies of CH2OH radical (0.49 eV) and CH3OH (0.41 eV) are relatively stronger compared to the CH3O radical binding energies (0.32 eV, Sameera et al., J. Phy. Chem. A, 2021, 125, 387-393). Thus, the CH3OH molecule, CH2OH and CH3O radicals can adsorb on ice, where the binding energies follow the order CH2OH > CH3OH > CH3O. The multi-component artificial force-induced reaction (MC-AFIR) method systematically determined the reaction mechanisms for the CH3OH + OH reaction on ice, where two reaction paths, giving rise to CH2OH and CH3O radicals, were confirmed. A range of reaction barriers, employing the ωB97X-D/Def2-TZVP level of theory, was found for each reaction (0.03-0.11 eV for CH2OH radical formation, and 0.03-0.44 eV for CH3O radical formation). Based on the lowest energy reaction paths, we suspect that both reactions operate on ice. The computed data in this study evidence that the nature of the binding site or the reaction site has a significant effect on the computed binding energies or reaction barriers. Thus, the outcomes of the present study will be very useful for the computational astrochemistry community to determine reliable binding energies and reaction barriers on ice.
Collapse
Affiliation(s)
- W M C Sameera
- Institute of Low Temperature Science, Hokkaido University, Sapporo, Hokkaido 060-0819, Japan.
- Department of Chemistry, University of Colombo, Colombo 00300, Sri Lanka
| | - Avon P Jayaweera
- Department of Chemistry, University of Colombo, Colombo 00300, Sri Lanka
| | - Atsuki Ishibashi
- Institute of Low Temperature Science, Hokkaido University, Sapporo, Hokkaido 060-0819, Japan.
| | - Hiroshi Hidaka
- Institute of Low Temperature Science, Hokkaido University, Sapporo, Hokkaido 060-0819, Japan.
| | - Yasuhiro Oba
- Institute of Low Temperature Science, Hokkaido University, Sapporo, Hokkaido 060-0819, Japan.
| | - Naoki Watanabe
- Institute of Low Temperature Science, Hokkaido University, Sapporo, Hokkaido 060-0819, Japan.
| |
Collapse
|
43
|
Cinq N, Simon A, Louisnard F, Cuny J. Accurate SCC-DFTB Parametrization of Liquid Water with Improved Atomic Charges and Iterative Boltzmann Inversion. J Phys Chem B 2023; 127:7590-7601. [PMID: 37603798 DOI: 10.1021/acs.jpcb.3c03479] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
This work presents improvements of the description of liquid water within the self-consistent-charge density-functional based tight-binding scheme combining the use of Weighted Mulliken (WMull) charges and optimized O-H repulsive potential through the iterative Boltzmann inversion (IBI) process. The quality of the newly developed models is validated considering pair radial distribution functions (RDFs), as well as other structural, energetic, thermodynamic, and dynamic properties. The use of WMull charges certainly improves the agreement with experimental data, however leading to over-structured RDFs at short distance, that can be further improved by considering an optimized O-H repulsive potential obtained by the IBI process. Three different schemes were used to optimize this potential: (i) optimization including short O-H distances. This led to accurate RDFs as well as improved self-diffusion coefficient and heat of vaporization, while the proton transfer energy barrier is severely deteriorated; (ii) optimization starting at long distance. The proton transfer energy barrier is recovered while the heat of vaporization is deteriorated and the O-H RDF is less accurate at short distance; (iii) optimization within the path-integral molecular dynamics scheme which allows us to exclude nuclear quantum effects from the repulsive potential. The latter potential, in conjunction with the WMull improved atomic charges, provides similar results as (i) for structural, dynamic, and thermodynamic properties while recovering a large part of the proton transfer energy barrier. It therefore offers a good compromise to study both dynamic properties and chemistry within liquid water at a quantum chemical level.
Collapse
Affiliation(s)
- Nicolas Cinq
- Laboratoire de Chimie et Physique Quantiques (LCPQ), FeRMI Institute, Université de Toulouse [UT3] and CNRS, Toulouse F-31062, France
| | - Aude Simon
- Laboratoire de Chimie et Physique Quantiques (LCPQ), FeRMI Institute, Université de Toulouse [UT3] and CNRS, Toulouse F-31062, France
| | - Fernand Louisnard
- Laboratoire de Chimie et Physique Quantiques (LCPQ), FeRMI Institute, Université de Toulouse [UT3] and CNRS, Toulouse F-31062, France
| | - Jérôme Cuny
- Laboratoire de Chimie et Physique Quantiques (LCPQ), FeRMI Institute, Université de Toulouse [UT3] and CNRS, Toulouse F-31062, France
| |
Collapse
|
44
|
Mulvey D, Jordan KD. Application of a Fluctuating Charge Polarization Model to Large Polyaromatic Hydrocarbons and Graphene Nanoflakes. J Phys Chem Lett 2023; 14:7869-7875. [PMID: 37639228 PMCID: PMC10494230 DOI: 10.1021/acs.jpclett.3c02013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 08/21/2023] [Indexed: 08/29/2023]
Abstract
We present a polarization model incorporating coupled fluctuating charges and point inducible dipoles that is able to accurately describe the dipole polarizabilities of small hydrocarbons and, for sufficiently large graphene nanoflakes, reproduce the classical image potential of an infinite conducting sheet. When our fluctuating charge model is applied to the hexagonal carbon nanoflake C60000 we attain excellent agreement with the image potential and induced charge distribution of a conducting sheet. With the inclusion of inducible dipole terms, the model predicts an image plane of zim = 1.3334 a0, which falls in line with prior estimates for graphene. We consider the case of two charges placed on opposite sides of C60000 and find that the fluctuating charge model reproduces classical electrostatics once again. By testing opposing and similar signs of the external charges, we conclude that an atomically thin molecule or extended system does not fully screen their interaction.
Collapse
Affiliation(s)
- Devin
M. Mulvey
- Department of Chemistry, University
of Pittsburgh, 219 Parkman Avenue, Pittsburgh, Pennsylvania 15260, United States
| | - Kenneth D. Jordan
- Department of Chemistry, University
of Pittsburgh, 219 Parkman Avenue, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
45
|
Huang Z, Zhao S, Cieplak P, Duan Y, Luo R, Wei H. Optimal Scheme to Achieve Energy Conservation in Induced Dipole Models. J Chem Theory Comput 2023; 19:5047-5057. [PMID: 37441805 PMCID: PMC10434752 DOI: 10.1021/acs.jctc.3c00226] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/15/2023]
Abstract
Induced dipole models have proven to be effective tools for simulating electronic polarization effects in biochemical processes, yet their potential has been constrained by energy conservation issue, particularly when historical data is utilized for dipole prediction. This study identifies error outliers as the primary factor causing this failure of energy conservation and proposes a comprehensive scheme to overcome this limitation. Leveraging maximum relative errors as a convergence metric, our data demonstrates that energy conservation can be upheld even when using historical information for dipole predictions. Our study introduces the multi-order extrapolation method to quicken induction iteration and optimize the use of historical data, while also developing the preconditioned conjugate gradient with local iterations to refine the iteration process and effectively remove error outliers. This scheme further incorporates a "peek" step via Jacobi under-relaxation for optimal performance. Simulation evidence suggests that our proposed scheme can achieve energy convergence akin to that of point-charge models within a limited number of iterations, thus promising significant improvements in efficiency and accuracy.
Collapse
Affiliation(s)
- Zhen Huang
- Chemical and Materials Physics Graduate Program, University of California, Irvine. Irvine, California 92697, United States
- Departments of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine. Irvine, California 92697, United States
| | - Shiji Zhao
- Departments of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine. Irvine, California 92697, United States
- Nurix Therapeutics, Inc., 1700 Owens St. Suite 205, San Francisco, CA 94158, United States
| | - Piotr Cieplak
- SBP Medical Discovery Institute, 10901 North Torrey Pines Road, La Jolla, California 92037, United States
| | - Yong Duan
- UC Davis Genome Center and Department of Biomedical Engineering, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - Ray Luo
- Chemical and Materials Physics Graduate Program, University of California, Irvine. Irvine, California 92697, United States
- Departments of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine. Irvine, California 92697, United States
| | - Haixin Wei
- Departments of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine. Irvine, California 92697, United States
- Department of Chemistry & Biochemistry, University of California, San Diego. 9500 Gilman Drive, La Jolla, California 92093, United States
| |
Collapse
|
46
|
Jaffrelot Inizan T, Plé T, Adjoua O, Ren P, Gökcan H, Isayev O, Lagardère L, Piquemal JP. Scalable hybrid deep neural networks/polarizable potentials biomolecular simulations including long-range effects. Chem Sci 2023; 14:5438-5452. [PMID: 37234902 PMCID: PMC10208042 DOI: 10.1039/d2sc04815a] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2022] [Accepted: 04/03/2023] [Indexed: 07/28/2023] Open
Abstract
Deep-HP is a scalable extension of the Tinker-HP multi-GPU molecular dynamics (MD) package enabling the use of Pytorch/TensorFlow Deep Neural Network (DNN) models. Deep-HP increases DNNs' MD capabilities by orders of magnitude offering access to ns simulations for 100k-atom biosystems while offering the possibility of coupling DNNs to any classical (FFs) and many-body polarizable (PFFs) force fields. It allows therefore the introduction of the ANI-2X/AMOEBA hybrid polarizable potential designed for ligand binding studies where solvent-solvent and solvent-solute interactions are computed with the AMOEBA PFF while solute-solute ones are computed by the ANI-2X DNN. ANI-2X/AMOEBA explicitly includes AMOEBA's physical long-range interactions via an efficient Particle Mesh Ewald implementation while preserving ANI-2X's solute short-range quantum mechanical accuracy. The DNN/PFF partition can be user-defined allowing for hybrid simulations to include key ingredients of biosimulation such as polarizable solvents, polarizable counter ions, etc.… ANI-2X/AMOEBA is accelerated using a multiple-timestep strategy focusing on the model's contributions to low-frequency modes of nuclear forces. It primarily evaluates AMOEBA forces while including ANI-2X ones only via correction-steps resulting in an order of magnitude acceleration over standard Velocity Verlet integration. Simulating more than 10 μs, we compute charged/uncharged ligand solvation free energies in 4 solvents, and absolute binding free energies of host-guest complexes from SAMPL challenges. ANI-2X/AMOEBA average errors are discussed in terms of statistical uncertainty and appear in the range of chemical accuracy compared to experiment. The availability of the Deep-HP computational platform opens the path towards large-scale hybrid DNN simulations, at force-field cost, in biophysics and drug discovery.
Collapse
Affiliation(s)
- Théo Jaffrelot Inizan
- Sorbonne Université, Laboratoire de Chimie Théorique UMR 7616 CNRS Paris 75005 France
| | - Thomas Plé
- Sorbonne Université, Laboratoire de Chimie Théorique UMR 7616 CNRS Paris 75005 France
| | - Olivier Adjoua
- Sorbonne Université, Laboratoire de Chimie Théorique UMR 7616 CNRS Paris 75005 France
| | - Pengyu Ren
- Department of Biomedical Engineering, University of Texas at Austin Austin Texas USA
| | - Hatice Gökcan
- Department of Chemistry, Carnegie Mellon University Pittsburgh Pennsylvania USA
| | - Olexandr Isayev
- Department of Chemistry, Carnegie Mellon University Pittsburgh Pennsylvania USA
| | - Louis Lagardère
- Sorbonne Université, Laboratoire de Chimie Théorique UMR 7616 CNRS Paris 75005 France
- Sorbonne Université, Institut Parisien de Chimie Physique et Théorique FR 2622 CNRS Paris France
| | - Jean-Philip Piquemal
- Sorbonne Université, Laboratoire de Chimie Théorique UMR 7616 CNRS Paris 75005 France
- Department of Biomedical Engineering, University of Texas at Austin Austin Texas USA
| |
Collapse
|
47
|
Rick SW, Thompson WH. Effects of polarizability and charge transfer on water dynamics and the underlying activation energies. J Chem Phys 2023; 158:2890774. [PMID: 37191215 DOI: 10.1063/5.0151253] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2023] [Accepted: 05/02/2023] [Indexed: 05/17/2023] Open
Abstract
A large number of force fields have been proposed for describing the behavior of liquid water within classical atomistic simulations, particularly molecular dynamics. In the past two decades, models that incorporate molecular polarizability and even charge transfer have become more prevalent, in attempts to develop more accurate descriptions. These are frequently parameterized to reproduce the measured thermodynamics, phase behavior, and structure of water. On the other hand, the dynamics of water is rarely considered in the construction of these models, despite its importance in their ultimate applications. In this paper, we explore the structure and dynamics of polarizable and charge-transfer water models, with a focus on timescales that directly or indirectly relate to hydrogen bond (H-bond) making and breaking. Moreover, we use the recently developed fluctuation theory for dynamics to determine the temperature dependence of these properties to shed light on the driving forces. This approach provides key insight into the timescale activation energies through a rigorous decomposition into contributions from the different interactions, including polarization and charge transfer. The results show that charge transfer effects have a negligible effect on the activation energies. Furthermore, the same tension between electrostatic and van der Waals interactions that is found in fixed-charge water models also governs the behavior of polarizable models. The models are found to involve significant energy-entropy compensation, pointing to the importance of developing water models that accurately describe the temperature dependence of water structure and dynamics.
Collapse
Affiliation(s)
- Steven W Rick
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, USA
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, USA
| |
Collapse
|
48
|
González-Pérez R, Adams S, Dowling AW, Phillip WA, Whitmer JK. Thermodynamics of Li +-Crown Ether Interactions in Aqueous Solvent. J Phys Chem A 2023. [PMID: 37196205 DOI: 10.1021/acs.jpca.3c00094] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/19/2023]
Abstract
Lithium ion-based batteries are ubiquitous in modern technology due to applications in personal electronics and high-capacity storage for electric vehicles. Concerns about lithium supply and battery waste have prompted interest in lithium recycling methods. The crown ether 12-crown-4 has been studied for its abilities to form stable complexes with lithium ions (Li+). In this paper, molecular dynamics simulations are applied to examine the binding properties of a 12-crown-4-Li+ system in aqueous solution. It was found that 12-crown-4 did not form stable complexes with Li+ in aqueous solution due to the binding geometry which was prone to interference by surrounding water molecules. In addition, the binding properties of sodium ions (Na+) to 12-crown-4 are examined for comparison. Subsequently, calculations were performed with the crown ethers 15-crown-5 and 18-crown-6 to study their complexation with Li+ as well as Na+. It was determined that binding was unfavorable for both types of ions for all three crown ethers tested, though 15-crown-5 and 18-crown-6 showed a marginally greater affinity for Li+ than 12-crown-4. Metastable minima present in the potential of mean force for Na+ render binding marginally more likely there. We discuss these results in the context of membrane-based applications of crown ethers for Li+ separations.
Collapse
Affiliation(s)
- Ramón González-Pérez
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Stephen Adams
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Alexander W Dowling
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - William A Phillip
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Jonathan K Whitmer
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, United States
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| |
Collapse
|
49
|
Chung MKJ, Miller RJ, Novak B, Wang Z, Ponder JW. Accurate Host-Guest Binding Free Energies Using the AMOEBA Polarizable Force Field. J Chem Inf Model 2023; 63:2769-2782. [PMID: 37075788 PMCID: PMC10878370 DOI: 10.1021/acs.jcim.3c00155] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/21/2023]
Abstract
A grand challenge of computational biophysics is accurate prediction of interactions between molecules. Molecular dynamics (MD) simulations have recently gained much interest as a tool to directly compute rigorous intermolecular binding affinities. The choice of a fixed point-charge or polarizable multipole force field used in MD is a topic of ongoing discussion. To compare alternative methods, we participated in the SAMPL7 and SAMPL8 Gibb octaacid host-guest challenges to assess the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) polarizable multipole force field. Advantages of AMOEBA over fixed charge models include improved representation of molecular electrostatic potentials and better description of water occupying the unligated host cavity. Prospective predictions for 26 host-guest systems exhibit a mean unsigned error vs experiment of 0.848 kcal/mol across all absolute binding free energies, demonstrating excellent agreement between computational and experimental results. In addition, we explore two topics related to the inclusion of ions in MD simulations: use of a neutral co-alchemical protocol and the effect of salt concentration on binding affinity. Use of the co-alchemical method minimally affects computed energies, but salt concentration significantly perturbs our binding results. Higher salt concentration strengthens binding through classical charge screening. In particular, added Na+ ions screen negatively charged carboxylate groups near the binding cavity, thereby diminishing repulsive coulomb interactions with negatively charged guests. Overall, the AMOEBA results demonstrate the accuracy available through a force field providing a detailed energetic description of the four octaacid hosts and 13 charged organic guests. Use of the AMOEBA polarizable atomic multipole force field in conjunction with an alchemical free energy protocol can achieve chemical accuracy in application to realistic molecular systems.
Collapse
Affiliation(s)
- Moses K. J. Chung
- Medical Scientist Training Program, Washington University School of Medicine, Saint Louis, MO 63110, USA
- Department of Physics, Washington University in St. Louis, Saint Louis, MO 63130, USA
| | - Ryan J. Miller
- Department of Chemistry, Washington University in St. Louis, Saint Louis, MO 63130, USA
| | - Borna Novak
- Medical Scientist Training Program, Washington University School of Medicine, Saint Louis, MO 63110, USA
| | - Zhi Wang
- Department of Chemistry, Washington University in St. Louis, Saint Louis, MO 63130, USA
| | - Jay W. Ponder
- Department of Chemistry, Washington University in St. Louis, Saint Louis, MO 63130, USA
- Department of Biochemistry & Molecular Biophysics, Washington University School of Medicine, Saint Louis, MO 63110, USA
| |
Collapse
|
50
|
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
|