1
|
Maruyama Y, Yoshida N. RISMiCal: A software package to perform fast RISM/3D-RISM calculations. J Comput Chem 2024; 45:1470-1482. [PMID: 38472097 DOI: 10.1002/jcc.27340] [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: 12/27/2023] [Revised: 02/23/2024] [Accepted: 02/28/2024] [Indexed: 03/14/2024]
Abstract
Solvent plays an essential role in a variety of chemical, physical, and biological processes that occur in the solution phase. The reference interaction site model (RISM) and its three-dimensional extension (3D-RISM) serve as powerful computational tools for modeling solvation effects in chemical reactions, biological functions, and structure formations. We present the RISM integrated calculator (RISMiCal) program package, which is based on RISM and 3D-RISM theories with fast GPU code. RISMiCal has been developed as an integrated RISM/3D-RISM program that has interfaces with external programs such as Gaussian16, GAMESS, and Tinker. Fast 3D-RISM programs for single- and multi-GPU codes written in CUDA would enhance the availability of these hybrid methods because they require the performance of many computationally expensive 3D-RISM calculations. We expect that our package can be widely applied for chemical and biological processes in solvent. The RISMiCal package is available at https://rismical-dev.github.io.
Collapse
Affiliation(s)
- Yutaka Maruyama
- Data Science Center for Creative Design and Manufacturing, The Institute of Statistical Mathematics, Tachikawa, Tokyo, Japan
- Department of Physics, School of Science and Technology, Meiji University, Kawasaki-shi, Kanagawa, Japan
| | - Norio Yoshida
- Graduate School of Informatics, Nagoya University, Chikusa, Nagoya, Japan
| |
Collapse
|
2
|
Fedotova MV, Chuev GN. The Three-Dimensional Reference Interaction Site Model Approach as a Promising Tool for Studying Hydrated Viruses and Their Complexes with Ligands. Int J Mol Sci 2024; 25:3697. [PMID: 38612508 PMCID: PMC11011341 DOI: 10.3390/ijms25073697] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2024] [Revised: 03/22/2024] [Accepted: 03/25/2024] [Indexed: 04/14/2024] Open
Abstract
Viruses are the most numerous biological form living in any ecosystem. Viral diseases affect not only people but also representatives of fauna and flora. The latest pandemic has shown how important it is for the scientific community to respond quickly to the challenge, including critically assessing the viral threat and developing appropriate measures to counter this threat. Scientists around the world are making enormous efforts to solve these problems. In silico methods, which allow quite rapid obtention of, in many cases, accurate information in this field, are effective tools for the description of various aspects of virus activity, including virus-host cell interactions, and, thus, can provide a molecular insight into the mechanism of virus functioning. The three-dimensional reference interaction site model (3D-RISM) seems to be one of the most effective and inexpensive methods to compute hydrated viruses, since the method allows us to provide efficient calculations of hydrated viruses, remaining all molecular details of the liquid environment and virus structure. The pandemic challenge has resulted in a fast increase in the number of 3D-RISM calculations devoted to hydrated viruses. To provide readers with a summary of this literature, we present a systematic overview of the 3D-RISM calculations, covering the period since 2010. We discuss various biophysical aspects of the 3D-RISM results and demonstrate capabilities, limitations, achievements, and prospects of the method using examples of viruses such as influenza, hepatitis, and SARS-CoV-2 viruses.
Collapse
Affiliation(s)
- Marina V. Fedotova
- G.A. Krestov Institute of Solution Chemistry, The Russian Academy of Sciences, Akademicheskaya St., 1, 153045 Ivanovo, Russia
| | - Gennady N. Chuev
- Institute of Theoretical and Experimental Biophysics, The Russian Academy of Sciences, Institutskaya St., 142290 Pushchino, Russia
| |
Collapse
|
3
|
Bruch N, Binninger T, Huang J, Eikerling M. Incorporating Electrolyte Correlation Effects into Variational Models of Electrochemical Interfaces. J Phys Chem Lett 2024; 15:2015-2022. [PMID: 38349906 PMCID: PMC10895655 DOI: 10.1021/acs.jpclett.3c03295] [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/23/2023] [Revised: 01/24/2024] [Accepted: 01/26/2024] [Indexed: 02/15/2024]
Abstract
We propose a way for obtaining a classical free energy density functional for electrolytes based on a first-principle many-body partition function. Via a one-loop expansion, we include coulombic correlations beyond the conventional mean-field approximation. To examine electrochemical interfaces, we integrate the electrolyte free energy functional into a hybrid quantum-classical model. This scheme self-consistently couples electronic, ionic, and solvent degrees of freedom and incorporates electrolyte correlation effects. The derived free energy functional causes a correlation-induced enhancement in interfacial counterion density and leads to an overall increase in capacitance. This effect is partially compensated by a reduction of the dielectric permittivity of interfacial water. At larger surface charge densities, ion crowding at the interface stifles these correlation effects. While scientifically intriguing already at planar interfaces, we anticipate these correlation effects to play an essential role for electrolytes in nanoconfinement.
Collapse
Affiliation(s)
- Nils Bruch
- Theory
and Computation of Energy Materials (IEK-13), Institute of Energy
and Climate Research, Forschungszentrum
Jülich GmbH, 52425, Jülich, Germany
- Chair
of Theory and Computation of Energy Materials, Faculty of Georesources
and Materials Engineering, RWTH Aachen University, 52062, Aachen Germany
| | - Tobias Binninger
- Theory
and Computation of Energy Materials (IEK-13), Institute of Energy
and Climate Research, Forschungszentrum
Jülich GmbH, 52425, Jülich, Germany
| | - Jun Huang
- Theory
and Computation of Energy Materials (IEK-13), Institute of Energy
and Climate Research, Forschungszentrum
Jülich GmbH, 52425, Jülich, Germany
- Chair
of Theory and Computation of Energy Materials, Faculty of Georesources
and Materials Engineering, RWTH Aachen University, 52062, Aachen Germany
| | - Michael Eikerling
- Theory
and Computation of Energy Materials (IEK-13), Institute of Energy
and Climate Research, Forschungszentrum
Jülich GmbH, 52425, Jülich, Germany
- Chair
of Theory and Computation of Energy Materials, Faculty of Georesources
and Materials Engineering, RWTH Aachen University, 52062, Aachen Germany
| |
Collapse
|
4
|
Imamura K, Yokogawa D, Sato H. Recent developments and applications of reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED): A hybrid model of quantum chemistry and integral equation theory of molecular liquids. J Chem Phys 2024; 160:050901. [PMID: 38341702 DOI: 10.1063/5.0190116] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2023] [Accepted: 01/04/2024] [Indexed: 02/13/2024] Open
Abstract
The significance of solvent effects in electronic structure calculations has long been noted, and various methods have been developed to consider this effect. The reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED) is a hybrid model that combines the integral equation theory of molecular liquids with quantum chemistry. This method can consider the statistically convergent solvent distribution at a significantly lower cost than molecular dynamics simulations. Because the RISM theory explicitly considers the solvent structure, it performs well for systems where hydrogen bonds are formed between the solute and solvent molecules, which is a challenge for continuum solvent models. Taking advantage of being founded on the variational principle, theoretical developments have been made in calculating various properties and incorporating electron correlation effects. In this review, we organize the theoretical aspects of RISM-SCF-cSED and its distinctions from other hybrid methods involving integral equation theories. Furthermore, we carefully present its progress in terms of theoretical developments and recent applications.
Collapse
Affiliation(s)
- Kosuke Imamura
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
| | - Daisuke Yokogawa
- Graduate School of Arts and Science, The University of Tokyo, Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| |
Collapse
|
5
|
Cao S, Kalin ML, Huang X. EPISOL: A software package with expanded functions to perform 3D-RISM calculations for the solvation of chemical and biological molecules. J Comput Chem 2023; 44:1536-1549. [PMID: 36856731 DOI: 10.1002/jcc.27088] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2022] [Revised: 12/24/2022] [Accepted: 01/29/2023] [Indexed: 03/02/2023]
Abstract
Integral equation theory (IET) provides an effective solvation model for chemical and biological systems that balances computational efficiency and accuracy. We present a new software package, the expanded package for IET-based solvation (EPISOL), that performs 3D-reference interaction site model (3D-RISM) calculations to obtain the solvation structure and free energies of solute molecules in different solvents. In EPISOL, we have implemented 22 different closures, multiple free energy functionals, and new variations of 3D-RISM theory, including the recent hydrophobicity-induced density inhomogeneity (HI) theory for hydrophobic solutes and ion-dipole correction (IDC) theory for negatively charged solutes. To speed up the convergence and enhance the stability of the self-consistent iterations, we have introduced several numerical schemes in EPISOL, including a newly developed dynamic mixing approach. We show that these schemes have significantly reduced the failure rate of 3D-RISM calculations compared to AMBER-RISM software. EPISOL consists of both a user-friendly graphic interface and a kernel library that allows users to call its routines and adapt them to other programs. EPISOL is compatible with the force-field and coordinate files from both AMBER and GROMACS simulation packages. Moreover, EPISOL is equipped with an internal memory control to efficiently manage the use of physical memory, making it suitable for performing calculations on large biomolecules. We demonstrate that EPISOL can efficiently and accurately calculate solvation density distributions around various solute molecules (including a protein chaperone consisting of 120,715 atoms) and obtain solvent free energy for a wide range of organic compounds. We expect that EPISOL can be widely applied as a solvation model for chemical and biological systems. EPISOL is available at https://github.com/EPISOLrelease/EPISOL.
Collapse
Affiliation(s)
- Siqin Cao
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin, USA
| | - Michael L Kalin
- Biophysics Graduate Program, University of Wisconsin-Madison, Madison, Wisconsin, USA
| | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin, USA
| |
Collapse
|
6
|
Raposo DJ. Effect of Conformational Equilibrium on Solvation Properties of 1,2-DCE in Water: A Solvation Thermodynamics and 3D-RISM Study. J Phys Chem B 2023; 127:757-765. [PMID: 36626710 DOI: 10.1021/acs.jpcb.2c07836] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
The contributions of the enthalpy and entropy of solvation for the study of chemical and biological systems are important in the prediction, interpretation, and manipulation of these processes. The relation between solvation Gibbs energies, enthalpies, and entropies of solvation, and their rigorous relation with the conformational equilibrium, are derived for the first time and applied with a computational method, in accordance with the Solvation Thermodynamics previous results, to 1,2-dichloroethane solvation in water. The rigid conformer calculations in solution were performed by using PC+/3D-RISM approach, with the conformational averaged results for enthalpy and solvation Gibbs energy reproducing the experimental results quite successfully. A qualitative agreement in the entropy of solvation predictions was also observed.
Collapse
Affiliation(s)
- Diego J Raposo
- Departamento de Química Fundamental, Universidade Federal de Pernambuco, Cidade Universitária, Recife, Pernambuco50740-560, Brazil
| |
Collapse
|
7
|
Roy D, Patel C. Revisiting the Use of Quantum Chemical Calculations in LogP octanol-water Prediction. Molecules 2023; 28:801. [PMID: 36677858 PMCID: PMC9866719 DOI: 10.3390/molecules28020801] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Revised: 01/06/2023] [Accepted: 01/10/2023] [Indexed: 01/15/2023] Open
Abstract
The partition coefficients of drug and drug-like molecules between an aqueous and organic phase are an important property for developing new therapeutics. The predictive power of computational methods is used extensively to predict partition coefficients of molecules. The application of quantum chemical calculations is used to develop methods to develop structure-activity relationship models for such prediction, either based on molecular fragment methods, or via direct calculation of solvation free energy in solvent continuum. The applicability, merits, and shortcomings of these developments are revisited here.
Collapse
Affiliation(s)
- Dipankar Roy
- Department of Biological Sciences, University of Alberta, Edmonton, AB T6G 2E9, Canada
| | - Chandan Patel
- Department of Applied Sciences, COEP Technological University, Wellesely Road, Shivajinagar, Pune 411005, Maharashtra, India
| |
Collapse
|
8
|
Bokhovchuk F, Mesrouze Y, Meyerhofer M, Fontana P, Zimmermann C, Villard F, Erdmann D, Kallen J, Clemens S, Velez‐Vega C, Chène P. N-terminal β-strand in YAP is critical for stronger binding to scalloped relative to TEAD transcription factor. Protein Sci 2023; 32:e4545. [PMID: 36522189 PMCID: PMC9798255 DOI: 10.1002/pro.4545] [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: 07/02/2022] [Revised: 11/30/2022] [Accepted: 12/06/2022] [Indexed: 12/23/2022]
Abstract
The yes-associated protein (YAP) regulates the transcriptional activity of the TEAD transcription factors that are key in the control of organ morphogenesis. YAP interacts with TEAD via three secondary structure elements: a β-strand, an α-helix, and an Ω-loop. Earlier results have shown that the β-strand has only a marginal contribution in the YAP:TEAD interaction, but we show here that it significantly enhances the affinity of YAP for the Drosophila homolog of TEAD, scalloped (Sd). Nuclear magnetic resonance shows that the β-strand adopts a more rigid conformation once bound to Sd; pre-steady state kinetic measurements show that the YAP:Sd complex is more stable. Although the crystal structures of the YAP:TEAD and YAP:Sd complexes reveal no differences at the binding interface that could explain these results. Molecular Dynamics simulations are in line with our experimental findings regarding β-strand stability and overall binding affinity of YAP to TEAD and Sd. In particular, RMSF, correlated motion and MMGBSA analyses suggest that β-sheet fluctuations play a relevant role in YAP53-57 β-strand dissociation from TEAD4 and contribute to the lower affinity of YAP for TEAD4. Identifying a clear mechanism leading to the difference in YAP's β-strand stability proved to be challenging, pointing to the potential relevance of multiple modest structural changes or fluctuations for regulation of binding affinity.
Collapse
Affiliation(s)
- Fedir Bokhovchuk
- Disease Area OncologyNovartis Institutes for Biomedical ResearchBaselSwitzerland
| | - Yannick Mesrouze
- Disease Area OncologyNovartis Institutes for Biomedical ResearchBaselSwitzerland
| | - Marco Meyerhofer
- Disease Area OncologyNovartis Institutes for Biomedical ResearchBaselSwitzerland
| | - Patrizia Fontana
- Disease Area OncologyNovartis Institutes for Biomedical ResearchBaselSwitzerland
| | - Catherine Zimmermann
- Disease Area OncologyNovartis Institutes for Biomedical ResearchBaselSwitzerland
| | - Frédéric Villard
- Chemical Biology and TherapeuticsNovartis Institutes for Biomedical ResearchBaselSwitzerland
| | - Dirk Erdmann
- Disease Area OncologyNovartis Institutes for Biomedical ResearchBaselSwitzerland
| | - Joerg Kallen
- Chemical Biology and TherapeuticsNovartis Institutes for Biomedical ResearchBaselSwitzerland
| | - Scheufler Clemens
- Chemical Biology and TherapeuticsNovartis Institutes for Biomedical ResearchBaselSwitzerland
| | - Camilo Velez‐Vega
- Global Discovery ChemistryNovartis Institutes for Biomedical ResearchCambridgeMassachusettsUSA
| | - Patrick Chène
- Disease Area OncologyNovartis Institutes for Biomedical ResearchBaselSwitzerland
| |
Collapse
|
9
|
Schauperl M, Denny RA. AI-Based Protein Structure Prediction in Drug Discovery: Impacts and Challenges. J Chem Inf Model 2022; 62:3142-3156. [PMID: 35727311 DOI: 10.1021/acs.jcim.2c00026] [Citation(s) in RCA: 36] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Proteins are the molecular machinery of the human body, and their malfunctioning is often responsible for diseases, making them crucial targets for drug discovery. The three-dimensional structure of a protein determines its biological function, its conformational state determines substrates, cofactors, and protein binding. Rational drug discovery employs engineered small molecules to selectively interact with proteins to modulate their function. To selectively target a protein and to design small molecules, knowing the protein structure with all its specific conformation is critical. Unfortunately, for a large number of proteins relevant for drug discovery, the three-dimensional structure has not yet been experimentally solved. Therefore, accurately predicting their structure based on their amino acid sequence is one of the grant challenges in biology. Recently, AlphaFold2, a machine learning application based on a deep neural network, was able to predict unknown structures of proteins with an unprecedented accuracy. Despite the impressive progress made by AlphaFold2, nature still challenges the field of structure prediction. In this Perspective, we explore how AlphaFold2 and related methods help make drug design more efficient. Furthermore, we discuss the roles of predicting domain-domain orientations, all relevant conformational states, the influence of posttranslational modifications, and conformational changes due to protein binding partners. We highlight where further improvements are needed for advanced machine learning methods to be successfully and frequently used in the pharmaceutical industry.
Collapse
Affiliation(s)
- Michael Schauperl
- Department of Computational Sciences HotSpot Therapeutics 50 Milk Street, Boston, Massachusetts 02110, United States
| | - Rajiah Aldrin Denny
- Department of Computational Sciences HotSpot Therapeutics 50 Milk Street, Boston, Massachusetts 02110, United States
| |
Collapse
|
10
|
Reimann M, Kaupp M. Reaction Entropies in Solution from Analytical Three-Dimensional Reference Interaction Site Model Derivatives with Application to Redox and Spin-Crossover Processes. J Phys Chem A 2022; 126:3708-3716. [PMID: 35652546 DOI: 10.1021/acs.jpca.2c02317] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An analytical approach to compute the excess entropy of solvation at constant pressure in three-dimensional reference interaction site model (3D-RISM) calculations is presented. It includes the changes in the macroscopic dielectric constant of the solvent upon variation of temperature and density. The approach is exact within the framework of force-field descriptions of the solute and gives reasonable results for self-consistently determined electrostatics as used in the 3D-RISM-self-consistent field approach, particularly for entropy differences. The new method is applied to simple examples of reaction entropies of iron complexes in aqueous solution, for which simple gas-phase calculations and many other approaches give unreliable estimates. For both redox half-reactions and spin-crossover processes, (semi)quantitative agreement with experimental reaction entropies can be achieved out of the box.
Collapse
Affiliation(s)
- Marc Reimann
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17.Juni 135, Berlin D-10623, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17.Juni 135, Berlin D-10623, Germany
| |
Collapse
|
11
|
Wilson L, Krasny R, Luchko T. Accelerating the 3D reference interaction site model theory of molecular solvation with treecode summation and cut-offs. J Comput Chem 2022; 43:1251-1270. [PMID: 35567580 DOI: 10.1002/jcc.26889] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Revised: 02/25/2022] [Accepted: 04/15/2022] [Indexed: 11/07/2022]
Abstract
The 3D reference interaction site model (3D-RISM) of molecular solvation is a powerful tool for computing the equilibrium thermodynamics and density distributions of solvents, such as water and co-ions, around solute molecules. However, 3D-RISM solutions can be expensive to calculate, especially for proteins and other large molecules where calculating the potential energy between solute and solvent requires more than half the computation time. To address this problem, we have developed and implemented treecode summation for long-range interactions and analytically corrected cut-offs for short-range interactions to accelerate the potential energy and long-range asymptotics calculations in non-periodic 3D-RISM in the AmberTools molecular modeling suite. For the largest single protein considered in this work, tubulin, the total computation time was reduced by a factor of 4. In addition, parallel calculations with these new methods scale almost linearly and the iterative solver remains the largest impediment to parallel scaling. To demonstrate the utility of our approach for large systems, we used 3D-RISM to calculate the solvation thermodynamics and density distribution of 7-ring microtubule, consisting of 910 tubulin dimers, over 1.2 million atoms.
Collapse
Affiliation(s)
- Leighton Wilson
- Department of Mathematics, University of Michigan, Ann Arbor, Michigan, USA
| | - Robert Krasny
- Department of Mathematics, University of Michigan, Ann Arbor, Michigan, USA
| | - Tyler Luchko
- Department of Physics and Astronomy, California State University, Los Angeles, California, USA
| |
Collapse
|
12
|
Ganyecz Á, Kállay M. Implementation and Optimization of the Embedded Cluster Reference Interaction Site Model with Atomic Charges. J Phys Chem A 2022; 126:2417-2429. [PMID: 35394778 PMCID: PMC9036516 DOI: 10.1021/acs.jpca.1c07904] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
In this work, we
implemented the embedded cluster reference interaction
site model (EC-RISM) originally developed by Kloss, Heil, and Kast
(J. Phys. Chem. B2008, 112, 4337–4343).
This method combines quantum mechanical calculations with the 3D reference
interaction site model (3D-RISM). Numerous options, such as buffer,
grid space, basis set, charge model, water model, closure relation,
and so forth, were investigated to find the best settings. Additionally,
the small point charges, which are derived from the solvent distribution
from the 3D-RISM solution to represent the solvent in the QM calculation,
were neglected to reduce the overhead without the loss of accuracy.
On the MNSOL[a], MNSOL, and FreeSolv databases, our implemented and
optimized method provides solvation free energies in water with 5.70,
6.32, and 6.44 kJ/mol root-mean-square deviations, respectively, but
with different settings, 5.22, 6.08, and 6.63 kJ/mol can also be achieved.
Only solvent models containing fitting parameters, like COSMO-RS and
EC-RISM with universal correction and directly used electrostatic
potential, perform better than our EC-RISM implementation with atomic
charges.
Collapse
Affiliation(s)
- Ádám Ganyecz
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest P.O. Box 91, H-1521 Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest P.O. Box 91, H-1521 Hungary
| |
Collapse
|
13
|
Salahub DR. Multiscale molecular modelling: from electronic structure to dynamics of nanosystems and beyond. Phys Chem Chem Phys 2022; 24:9051-9081. [PMID: 35389399 DOI: 10.1039/d1cp05928a] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Important contemporary biological and materials problems often depend on interactions that span orders of magnitude differences in spatial and temporal dimensions. This Tutorial Review attempts to provide an introduction to such fascinating problems through a series of case studies, aimed at beginning researchers, graduate students, postdocs and more senior colleagues who are changing direction to focus on multiscale aspects of their research. The choice of specific examples is highly personal, with examples either chosen from our own work or outstanding multiscale efforts from the literature. I start with various embedding schemes, as exemplified by polarizable continuum models, 3-D RISM, molecular DFT and frozen-density embedding. Next, QM/MM (quantum mechanical/molecular mechanical) techniques are the workhorse of pm-to-nm/ps-to-ns simulations; examples are drawn from enzymes and from nanocatalysis for oil-sands upgrading. Using polarizable force-fields in the QM/MM framework represents a burgeoning subfield; with examples from ion channels and electron dynamics in molecules subject to strong external fields, probing the atto-second dynamics of the electrons with RT-TDDFT (real-time - time-dependent density functional theory) eventually coupled with nuclear motion through the Ehrenfest approximation. This is followed by a section on coarse graining, bridging dimensions from atoms to cells. The penultimate chapter gives a quick overview of multiscale approaches that extend into the meso- and macro-scales, building on atomistic and coarse-grained techniques to enter the world of materials engineering, on the one hand, and cell biology, on the other. A final chapter gives just a glimpse of the burgeoning impact of machine learning on the structure-dynamics front. I aim to capture the excitement of contemporary leading-edge breakthroughs in the description of physico-chemical systems and processes in complex environments, with only enough historical content to provide context and aid the next generation of methodological development. While I aim also for a clear description of the essence of methodological breakthroughs, equations are kept to a minimum and detailed formalism and implementation details are left to the references. My approach is very selective (case studies) rather than exhaustive. I think that these case studies should provide fodder to build as complete a reference tree on multiscale modelling as the reader may wish, through forward and backward citation analysis. I hope that my choices of cases will excite interest in newcomers and help to fuel the growth of multiscale modelling in general.
Collapse
Affiliation(s)
- Dennis R Salahub
- Department of Chemistry, Department of Physics and Astronomy, CMS-Centre for Molecular Simulation, IQST-Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary, Calgary, Alberta, T2N 1N4, Canada.
| |
Collapse
|
14
|
Ringe S, Hörmann NG, Oberhofer H, Reuter K. Implicit Solvation Methods for Catalysis at Electrified Interfaces. Chem Rev 2021; 122:10777-10820. [PMID: 34928131 PMCID: PMC9227731 DOI: 10.1021/acs.chemrev.1c00675] [Citation(s) in RCA: 51] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
Implicit solvation
is an effective, highly coarse-grained approach
in atomic-scale simulations to account for a surrounding liquid electrolyte
on the level of a continuous polarizable medium. Originating in molecular
chemistry with finite solutes, implicit solvation techniques are now
increasingly used in the context of first-principles modeling of electrochemistry
and electrocatalysis at extended (often metallic) electrodes. The
prevalent ansatz to model the latter electrodes and the reactive surface
chemistry at them through slabs in periodic boundary condition supercells
brings its specific challenges. Foremost this concerns the difficulty
of describing the entire double layer forming at the electrified solid–liquid
interface (SLI) within supercell sizes tractable by commonly employed
density functional theory (DFT). We review liquid solvation methodology
from this specific application angle, highlighting in particular its
use in the widespread ab initio thermodynamics approach
to surface catalysis. Notably, implicit solvation can be employed
to mimic a polarization of the electrode’s electronic density
under the applied potential and the concomitant capacitive charging
of the entire double layer beyond the limitations of the employed
DFT supercell. Most critical for continuing advances of this effective
methodology for the SLI context is the lack of pertinent (experimental
or high-level theoretical) reference data needed for parametrization.
Collapse
Affiliation(s)
- Stefan Ringe
- Department of Energy Science and Engineering, Daegu Institute of Science and Technology (DGIST), Daegu 42988, Republic of Korea.,Energy Science & Engineering Research Center, Daegu Institute of Science and Technology (DGIST), Daegu 42988, Republic of Korea
| | - Nicolas G Hörmann
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany.,Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstraße 4, D-85747 Garching, Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstraße 4, D-85747 Garching, Germany.,Chair for Theoretical Physics VII and Bavarian Center for Battery Technology (BayBatt), University of Bayreuth, Universitätsstraße 30, 95447 Bayreuth, Germany
| | - Karsten Reuter
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| |
Collapse
|
15
|
Schuitemaker A, Aufort J, Koziara KB, Demichelis R, Raiteri P, Gale JD. Simulating the binding of key organic functional groups to aqueous calcium carbonate species. Phys Chem Chem Phys 2021; 23:27253-27265. [PMID: 34870292 DOI: 10.1039/d1cp04226b] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The interaction of organic molecules with mineral systems is relevant to a wide variety of scientific problems both in the environment and minerals processing. In this study, the coordination of small organics that contain the two most relevant functional groups for biomineralisation of calcium carbonate, namely carboxylate and ammonium, with the corresponding mineral ions are examined in aqueous solution. Specifically, two force fields have been examined based on rigid-ion or polarisable models, with the latter being within the AMOEBA formalism. Here the parameters for the rigid-ion model are determined to target the accurate reproduction of the hydration structure and solvation thermodynamics, while both force fields are designed to be compatible with the corresponding recently published models for aqueous calcium carbonate. The application of these force fields to ion pairing in aqueous solution is studied in order to quantitatively determine the extent of association.
Collapse
Affiliation(s)
- Alicia Schuitemaker
- Curtin Institute for Computation, The Institute for Geoscience Research (TIGeR), School of Molecular and Life Sciences, Curtin University, GPO Box U1987, 6845 Perth, Western Australia, Australia.
| | - Julie Aufort
- Curtin Institute for Computation, The Institute for Geoscience Research (TIGeR), School of Molecular and Life Sciences, Curtin University, GPO Box U1987, 6845 Perth, Western Australia, Australia.
| | - Katarzyna B Koziara
- Curtin Institute for Computation, The Institute for Geoscience Research (TIGeR), School of Molecular and Life Sciences, Curtin University, GPO Box U1987, 6845 Perth, Western Australia, Australia.
| | - Raffaella Demichelis
- Curtin Institute for Computation, The Institute for Geoscience Research (TIGeR), School of Molecular and Life Sciences, Curtin University, GPO Box U1987, 6845 Perth, Western Australia, Australia.
| | - Paolo Raiteri
- Curtin Institute for Computation, The Institute for Geoscience Research (TIGeR), School of Molecular and Life Sciences, Curtin University, GPO Box U1987, 6845 Perth, Western Australia, Australia.
| | - Julian D Gale
- Curtin Institute for Computation, The Institute for Geoscience Research (TIGeR), School of Molecular and Life Sciences, Curtin University, GPO Box U1987, 6845 Perth, Western Australia, Australia.
| |
Collapse
|
16
|
Garay-Ruiz D, Bo C. Rationalizing the Mechanism of Peroxyformate Decomposition: Computational Insights To Understand Solvent Influence. Chemistry 2021; 27:11618-11626. [PMID: 34076322 PMCID: PMC8457178 DOI: 10.1002/chem.202100755] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2021] [Indexed: 11/29/2022]
Abstract
The heterolytic decomposition of tert‐butyl peroxyformate to tert‐butanol and carbon dioxide, catalyzed by pyridine, is a long‐known example of a reaction whose kinetics are strongly affected by solvent polarity. From DFT and ab initio methods together with the SMD implicit solvation model, an extension on the formerly accepted mechanism is proposed. This novel proposal involves the formation of a carbonic acid ester intermediate and its further decomposition, through an unreported pyridine‐mediated stepwise route. Computed barriers for this mechanism at DLPNO/CCSD(T)‐def2‐TZVP are in excellent agreement with experimental kinetic data across different solvents. Furthermore, the strong relationships between activation energies, geometric parameters in the transition state and the characteristics of the different solvents are also analyzed in depth.
Collapse
Affiliation(s)
- Diego Garay-Ruiz
- Institute of Chemical Research of Catalonia (ICIQ), Barcelona Institute of Science & Technology (BIST), Av. Països Catalans, 16, 43007, Tarragona, Spain
| | - Carles Bo
- Institute of Chemical Research of Catalonia (ICIQ), Barcelona Institute of Science & Technology (BIST), Av. Països Catalans, 16, 43007, Tarragona, Spain.,Departament de Química Física i Inorgànica, Universitat Rovira i Virgili (URV) C/ Marcel⋅lí Domingo s/n, 43007, Tarragona, Spain
| |
Collapse
|
17
|
Biomolecular Simulations with the Three-Dimensional Reference Interaction Site Model with the Kovalenko-Hirata Closure Molecular Solvation Theory. Int J Mol Sci 2021; 22:ijms22105061. [PMID: 34064655 PMCID: PMC8151972 DOI: 10.3390/ijms22105061] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2021] [Revised: 05/08/2021] [Accepted: 05/10/2021] [Indexed: 11/17/2022] Open
Abstract
The statistical mechanics-based 3-dimensional reference interaction site model with the Kovalenko-Hirata closure (3D-RISM-KH) molecular solvation theory has proven to be an essential part of a multiscale modeling framework, covering a vast region of molecular simulation techniques. The successful application ranges from the small molecule solvation energy to the bulk phase behavior of polymers, macromolecules, etc. The 3D-RISM-KH successfully predicts and explains the molecular mechanisms of self-assembly and aggregation of proteins and peptides related to neurodegeneration, protein-ligand binding, and structure-function related solvation properties. Upon coupling the 3D-RISM-KH theory with a novel multiple time-step molecular dynamic (MD) of the solute biomolecule stabilized by the optimized isokinetic Nosé-Hoover chain thermostat driven by effective solvation forces obtained from 3D-RISM-KH and extrapolated forward by generalized solvation force extrapolation (GSFE), gigantic outer time-steps up to picoseconds to accurately calculate equilibrium properties were obtained in this new quasidynamics protocol. The multiscale OIN/GSFE/3D-RISM-KH algorithm was implemented in the Amber package and well documented for fully flexible model of alanine dipeptide, miniprotein 1L2Y, and protein G in aqueous solution, with a solvent sampling rate ~150 times faster than a standard MD simulation in explicit water. Further acceleration in computation can be achieved by modifying the extent of solvation layers considered in the calculation, as well as by modifying existing closure relations. This enhanced simulation technique has proven applications in protein-ligand binding energy calculations, ligand/solvent binding site prediction, molecular solvation energy calculations, etc. Applications of the RISM-KH theory in molecular simulation are discussed in this work.
Collapse
|
18
|
Tanimoto S, Tamura K, Hayashi S, Yoshida N, Nakano H. A computational method to simulate global conformational changes of proteins induced by cosolvent. J Comput Chem 2021; 42:552-563. [PMID: 33433010 DOI: 10.1002/jcc.26481] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2020] [Revised: 12/09/2020] [Accepted: 12/29/2020] [Indexed: 12/14/2022]
Abstract
A computational method to investigate the global conformational change of a protein is proposed by combining the linear response path following (LRPF) method and three-dimensional reference interaction site model (3D-RISM) theory, which is referred to as the LRPF/3D-RISM method. The proposed method makes it possible to efficiently simulate protein conformational changes caused by either solutions of varying concentrations or the presence of cosolvent species by taking advantage of the LRPF and 3D-RISM. The proposed method is applied to the urea-induced denaturation of ubiquitin. The LRPF/3D-RISM trajectories successfully simulate the early stage of the denaturation process within the simulation time of 300 ns, whereas no significant structural change is observed even in the 1 μs standard MD simulation. The obtained LRPF/3D-RISM trajectories reproduce the mechanism of the urea denaturation of ubiquitin reported in previous studies, and demonstrate the high efficiency of the method.
Collapse
Affiliation(s)
- Shoichi Tanimoto
- Department of Chemistry, Graduate School of Science, Kyushu University, Fukuoka, Japan
| | - Koichi Tamura
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Shigehiko Hayashi
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Norio Yoshida
- Department of Chemistry, Graduate School of Science, Kyushu University, Fukuoka, Japan
| | - Haruyuki Nakano
- Department of Chemistry, Graduate School of Science, Kyushu University, Fukuoka, Japan
| |
Collapse
|
19
|
Kobryn AE, Maruyama Y, Velázquez-Martínez CA, Yoshida N, Gusarov S. Modeling the interaction of SARS-CoV-2 binding to the ACE2 receptor via molecular theory of solvation. NEW J CHEM 2021. [DOI: 10.1039/d1nj02015c] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The angiotensin-converting enzyme 2 (ACE2) protein is a cell gate receptor for the SARS-CoV-2 virus, responsible for the development of symptoms associated with the Covid-19 disease.
Collapse
Affiliation(s)
- Alexander E. Kobryn
- Nanotechnology Research Centre, National Research Council Canada, 11421 Saskatchewan Drive NW, Edmonton, Alberta, T6G 2M9, Canada
| | - Yutaka Maruyama
- Architecture Development Team, FLAGSHIP 2020 Project, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Carlos A. Velázquez-Martínez
- 2142-L Katz Group Centre for Research, University of Alberta, 11315-87 Avenue NW, Edmonton, Alberta, T6G 2H5, Canada
| | - Norio Yoshida
- Department of Chemistry, Graduate School of Science, Kyushu University, 744, Motooka, Nishi-ku, Fukuoka, 819-0395, Japan
| | - Sergey Gusarov
- Nanotechnology Research Centre, National Research Council Canada, 11421 Saskatchewan Drive NW, Edmonton, Alberta, T6G 2M9, Canada
| |
Collapse
|
20
|
Jeanmairet G, Levesque M, Borgis D. Tackling Solvent Effects by Coupling Electronic and Molecular Density Functional Theory. J Chem Theory Comput 2020; 16:7123-7134. [PMID: 32894674 DOI: 10.1021/acs.jctc.0c00729] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Solvation effects can have a tremendous influence on chemical reactions. However, precise quantum chemistry calculations are most often done either in vacuum neglecting the role of the solvent or using continuum solvent model ignoring its molecular nature. We propose a new method coupling a quantum description of the solute using electronic density functional theory with a classical grand-canonical treatment of the solvent using molecular density functional theory. Unlike a previous work, both densities are minimized self-consistently, accounting for mutual polarization of the molecular solvent and the solute. The electrostatic interaction is accounted using the full electron density of the solute rather than fitted point charges. The introduced methodology represents a good compromise between the two main strategies to tackle solvation effects in quantum calculation. It is computationally more effective than a direct quantum mechanics/molecular mechanics coupling, requiring the exploration of many solvent configurations. Compared to continuum methods, it retains the full molecular-level description of the solvent. We validate this new framework onto two usual benchmark systems: a water solvated in water and the symmetrical nucleophilic substitution between chloromethane and chloride in water. The prediction for the free energy profiles are not yet fully quantitative compared to experimental data, but the most important features are qualitatively recovered. The method provides a detailed molecular picture of the evolution of the solvent structure along the reaction pathway.
Collapse
Affiliation(s)
- Guillaume Jeanmairet
- Sorbonne Université, CNRS, Physico-Chimie des Électrolytes et Nanosystèmes, Interfaciaux, PHENIX, F-75005 Paris, France.,Réseau sur le Stockage Électrochimique de l'Énergie (RS2E), FR CNRS 3459, 80039 Amiens Cedex, France
| | - Maximilien Levesque
- PASTEUR, Département de chimie, École normale supérieure, PSL University, Sorbonne, Université, CNRS, 75005 Paris, France.,Aqemia, 75006 Paris, France
| | - Daniel Borgis
- PASTEUR, Département de chimie, École normale supérieure, PSL University, Sorbonne, Université, CNRS, 75005 Paris, France.,Maison de la Simulation, CEA, CNRS, Université Paris-Sud, UVSQ, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| |
Collapse
|
21
|
Reimann M, Kaupp M. Evaluation of an Efficient 3D-RISM-SCF Implementation as a Tool for Computational Spectroscopy in Solution. J Phys Chem A 2020; 124:7439-7452. [PMID: 32838530 DOI: 10.1021/acs.jpca.0c06322] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
The 3D-RISM-SCF solvent-model implementation of Gusarov et al. [ J. Phys. Chem. A 2006, 110, 6083-6090] in the Amsterdam density functional program has been improved and extended. In particular, an accurate yet efficient representation of the solute electrostatic potential is provided. The Coulomb-potential fitting of many DFT codes can be used advantageously in this context. The extra effort compared to a point-charge representation is small for a given SCF cycle and compensated by faster SCF convergence. This allows applications to large solutes, as demonstrated by evaluation of the solvatochromism of Reichardt's dye. In general, TDDFT applications to excitation energies in solution stand out and are highlighted. Applications to the 17O NMR chemical shifts of N-methylformamide in different solvents also demonstrate the distinct advantages of 3D-RISM over continuum solvents. Limitations are observed in this case for water solvent, where the solvent shielding is overestimated. This shortcoming applies also to the 17O gas-to-liquid shift of water, where we used localized molecular orbital analyses for a deeper understanding. For such cases of extremely strong solute-solvent interactions, couplings between solute and solvent orbitals induced by the magnetic perturbation are relevant. These clearly require a quantum-mechanical treatment of the most closely bound solvent molecules. Except for such extreme cases, 3D-RISM-SCF is very well suited to treat solvent effects on NMR parameters. More serious limitations pertain to the treatment of vibrational spectra, where the absence of the coupling between solute and solvent vibrational modes limits the accuracy of applications of 3D-RISM-SCF. The reported extended, efficient, and numerically accurate 3D-RISM-SCF implementation should provide a useful tool to study chemical and spectroscopic properties of molecules of appreciable size in a realistic solvent environment.
Collapse
Affiliation(s)
- Marc Reimann
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623, Berlin, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623, Berlin, Germany
| |
Collapse
|
22
|
Schwarz K, Sundararaman R. The electrochemical interface in first-principles calculations. SURFACE SCIENCE REPORTS 2020; 75:10.1016/j.surfrep.2020.100492. [PMID: 34194128 PMCID: PMC8240516 DOI: 10.1016/j.surfrep.2020.100492] [Citation(s) in RCA: 49] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
First-principles predictions play an important role in understanding chemistry at the electrochemical interface. Electronic structure calculations are straightforward for vacuum interfaces, but do not easily account for the interfacial fields and solvation that fundamentally change the nature of electrochemical reactions. Prevalent techniques for first-principles prediction of electrochemical processes range from expensive explicit solvation using ab initio molecular dynamics, through a hierarchy of continuum solvation techniques, to neglecting solvation and interfacial field effects entirely. Currently, no single approach reliably captures all relevant effects of the electrochemical double layer in first-principles calculations. This review systematically lays out the relation between all major approaches to first-principles electrochemistry, including the key approximations and their consequences for accuracy and computational cost. Focusing on ab initio methods for thermodynamic properties of aqueous interfaces, we first outline general considerations for modeling electrochemical interfaces, including solvent and electrolyte dynamics and electrification. We then present the specifics of various explicit and implicit models of the solvent and electrolyte. Finally, we discuss the compromise between computational efficiency and accuracy, and identify key outstanding challenges and future opportunities in the wide range of techniques for first-principles electrochemistry.
Collapse
Affiliation(s)
- Kathleen Schwarz
- Material Measurement Laboratory, National Institute of Standards and Technology, 100 Bureau Dr., Gaithersburg, Maryland 20899, USA
| | - Ravishankar Sundararaman
- Department of Materials Science and Engineering, Rensselaer Polytechnic Institute, 110 8th St., Troy, New York 12180, USA
| |
Collapse
|
23
|
Maldonado AM, Basdogan Y, Berryman JT, Rempe SB, Keith JA. First-principles modeling of chemistry in mixed solvents: Where to go from here? J Chem Phys 2020; 152:130902. [PMID: 32268733 DOI: 10.1063/1.5143207] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
Mixed solvents (i.e., binary or higher order mixtures of ionic or nonionic liquids) play crucial roles in chemical syntheses, separations, and electrochemical devices because they can be tuned for specific reactions and applications. Apart from fully explicit solvation treatments that can be difficult to parameterize or computationally expensive, there is currently no well-established first-principles regimen for reliably modeling atomic-scale chemistry in mixed solvent environments. We offer our perspective on how this process could be achieved in the near future as mixed solvent systems become more explored using theoretical and computational chemistry. We first outline what makes mixed solvent systems far more complex compared to single-component solvents. An overview of current and promising techniques for modeling mixed solvent environments is provided. We focus on so-called hybrid solvation treatments such as the conductor-like screening model for real solvents and the reference interaction site model, which are far less computationally demanding than explicit simulations. We also propose that cluster-continuum approaches rooted in physically rigorous quasi-chemical theory provide a robust, yet practical, route for studying chemical processes in mixed solvents.
Collapse
Affiliation(s)
- Alex M Maldonado
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Yasemin Basdogan
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Joshua T Berryman
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Susan B Rempe
- Center for Computational Biology and Biophysics, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - John A Keith
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| |
Collapse
|
24
|
Okamoto D, Watanabe Y, Yoshida N, Nakano H. Implementation of state-averaged MCSCF method to RISM- and 3D-RISM-SCF schemes. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.05.051] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
25
|
Roy D, Kovalenko A. Performance of 3D-RISM-KH in Predicting Hydration Free Energy: Effect of Solute Parameters. J Phys Chem A 2019; 123:4087-4093. [PMID: 30993994 DOI: 10.1021/acs.jpca.9b01623] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The three-dimensional reference interaction site model molecular solvation theory with the Kovalenko-Hirata closure relation has been shown to produce excellent solvation characteristics for a large class of (bio)chemical systems in solution. Correct calculation of hydration free energy is central to successful application of any solvation model. In order to find out the best possible force-field parameters to be used for hydration free energy calculation with the aforementioned theory, we have developed an extended database containing a large number of experimental solvation free energies available in the current literature and used a plethora of theoretical models for assessment. The general Amber force field was found to perform satisfactorily, whereas special care should be taken in solute charge assignment with the universal force field.
Collapse
Affiliation(s)
- Dipankar Roy
- Department of Mechanical Engineering , University of Alberta , 10-203 Donadeo Innovation Centre for Engineering, 9211-116 Street NW , Edmonton , Alberta T6G 1H9 , Canada
| | - Andriy Kovalenko
- Department of Mechanical Engineering , University of Alberta , 10-203 Donadeo Innovation Centre for Engineering, 9211-116 Street NW , Edmonton , Alberta T6G 1H9 , Canada.,Nanotechnology Research Centre , 11421 Saskatchewan Drive , Edmonton , Alberta T6G 2M9 , Canada
| |
Collapse
|
26
|
Fujiki R, Kasai Y, Seno Y, Matsui T, Shigeta Y, Yoshida N, Nakano H. A computational scheme of pK a values based on the three-dimensional reference interaction site model self-consistent field theory coupled with the linear fitting correction scheme. Phys Chem Chem Phys 2018; 20:27272-27279. [PMID: 30167611 DOI: 10.1039/c8cp04354j] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
A scheme for quantitatively computing the acid dissociation constant, pKa, of hydrated molecules is proposed. It is based on the three-dimensional reference interaction site model self-consistent field (3D-RISM-SCF) theory coupled with the linear fitting correction (LFC) scheme. In LFC/3D-RISM-SCF, pKa values of target molecules are evaluated using the Gibbs energy difference between the protonated and unprotonated states calculated by 3D-RISM-SCF and the parameters fitted by the LFC scheme to the experimental values of training set systems. The pKa values computed by LFC/3D-RISM-SCF show quantitative agreement with the experimental data.
Collapse
Affiliation(s)
- Ryo Fujiki
- Department of Chemistry, Graduate School of Science, Kyushu University, 744, Motooka, Nishiku, Fukuoka, 819-0395, Japan.
| | | | | | | | | | | | | |
Collapse
|
27
|
Kovalenko A, Gusarov S. Multiscale methods framework: self-consistent coupling of molecular theory of solvation with quantum chemistry, molecular simulations, and dissipative particle dynamics. Phys Chem Chem Phys 2018; 20:2947-2969. [DOI: 10.1039/c7cp05585d] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
In this work, we will address different aspects of self-consistent field coupling of computational chemistry methods at different time and length scales in modern materials and biomolecular science.
Collapse
Affiliation(s)
- Andriy Kovalenko
- National Institute for Nanotechnology
- National Research Council of Canada
- Edmonton
- Canada
- Department of Mechanical Engineering
| | - Sergey Gusarov
- National Institute for Nanotechnology
- National Research Council of Canada
- Edmonton
- Canada
| |
Collapse
|
28
|
Roy D, Blinov N, Kovalenko A. Predicting Accurate Solvation Free Energy in n-Octanol Using 3D-RISM-KH Molecular Theory of Solvation: Making Right Choices. J Phys Chem B 2017; 121:9268-9273. [PMID: 28880087 DOI: 10.1021/acs.jpcb.7b06375] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular theory of solvation, a.k.a., three-dimensional reference interaction site model theory of solvation with Kovalenko-Hirata closure relation (3D-RISM-KH), is an accurate and fast theory predicting solvation free energy and structure. Here we report a benchmark study of n-octanol solvation free energy calculations using this theory. The choice of correct force field parameters is quintessential for the success of 3D-RISM theory, and we present a guideline to obtain them for n-octanol solvent. Our best prediction of the solvation free energy on a set of 205 small organic molecules supplemented with the so-called "universal correction" scheme yields relative mean unsigned error of 0.94 kcal/mol against the reported database. The best agreement is obtained with the united atom (UA) type force field parametrization of n-octanol with the van der Waals parameters of hydroxyl hydrogen reported by Kobryn et al. [ Kobryn , A. E. ; Kovalenko , A. J. Chem. Phys. 2008 , 129 , 134701 ].
Collapse
Affiliation(s)
- Dipankar Roy
- Department of Mechanical Engineering, University of Alberta 10-203 Donadeo Innovation Centre for Engineering , 9211-116 Street NW, Edmonton, Alberta T6G 1H9, Canada.,National Institute for Nanotechnology , 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| | - Nikolay Blinov
- Department of Mechanical Engineering, University of Alberta 10-203 Donadeo Innovation Centre for Engineering , 9211-116 Street NW, Edmonton, Alberta T6G 1H9, Canada.,National Institute for Nanotechnology , 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| | - Andriy Kovalenko
- Department of Mechanical Engineering, University of Alberta 10-203 Donadeo Innovation Centre for Engineering , 9211-116 Street NW, Edmonton, Alberta T6G 1H9, Canada.,National Institute for Nanotechnology , 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| |
Collapse
|
29
|
Nocito D, Beran GJO. Averaged Condensed Phase Model for Simulating Molecules in Complex Environments. J Chem Theory Comput 2017; 13:1117-1129. [PMID: 28170251 DOI: 10.1021/acs.jctc.6b00890] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The need for configurational sampling dramatically increases the cost of combined quantum mechanics/molecular mechanics (QM/MM) simulations of chemical processes in solution. We developed an averaged condensed phase environment (ACPE) model that constructs an effective polarizable environment directly from explicitly sampled molecular dynamics configurations via the K-means++ algorithm and a mathematically rigorous translation of the molecular mechanics parameters. The model captures detailed heterogeneous features in the environment that may be difficult to describe using a conventional polarizable continuum model. Instead of performing repeated QM/MM calculations for each new configuration of the environment, the ACPE approach allows one to perform a single QM calculation on an averaged configuration. Here, we demonstrate the model by computing electronic excitation energies for several small molecules in solution. The ACPE model predicts the excitation energies in excellent agreement with conventional configurational averaging yet with orders of magnitude of reduction in the computational cost.
Collapse
Affiliation(s)
- Dominique Nocito
- Department of Chemistry, University of California , Riverside, California 92521, United States
| | - Gregory J O Beran
- Department of Chemistry, University of California , Riverside, California 92521, United States
| |
Collapse
|
30
|
Aono S, Nakagaki M, Sakaki S. Theoretical study of one-electron-oxidized salen complexes of group 7 (Mn(iii), Tc(iii), and Re(iii)) and group 10 metals (Ni(ii), Pd(ii), and Pt(ii)) with the 3D-RISM-GMC-QDPT method: localized vs. delocalized ground and excited states in solution. Phys Chem Chem Phys 2017. [DOI: 10.1039/c7cp02992f] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Mixed-valence natures of one-electron oxidized salen complexes of group 7 and 10 metals are theoretically investigated by GMC-QDPT, incorporating solvation effects.
Collapse
Affiliation(s)
- Shinji Aono
- Fukui Institute for Fundamental Chemistry
- Kyoto University
- Sakyo-ku
- Japan
| | - Masayuki Nakagaki
- Fukui Institute for Fundamental Chemistry
- Kyoto University
- Sakyo-ku
- Japan
| | - Shigeyoshi Sakaki
- Fukui Institute for Fundamental Chemistry
- Kyoto University
- Sakyo-ku
- Japan
| |
Collapse
|
31
|
Kobryn AE, Gusarov S, Kovalenko A. A closure relation to molecular theory of solvation for macromolecules. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:404003. [PMID: 27549008 DOI: 10.1088/0953-8984/28/40/404003] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We propose a closure to the integral equations of molecular theory of solvation, particularly suitable for polar and charged macromolecules in electrolyte solution. This includes such systems as oligomeric polyelectrolytes at a finite concentration in aqueous and various non-aqueous solutions, as well as drug-like compounds in solution. The new closure by Kobryn, Gusarov, and Kovalenko (KGK closure) imposes the mean spherical approximation (MSA) almost everywhere in the solvation shell but levels out the density distribution function to zero (with the continuity at joint boundaries) inside the repulsive core and in the spatial regions of strong density depletion emerging due to molecular associative interactions. Similarly to MSA, the KGK closure reduces the problem to a linear equation for the direct correlation function which is predefined analytically on most of the solvation shells and has to be determined numerically on a relatively small (three-dimensional) domain of strong depletion, typically within the repulsive core. The KGK closure leads to the solvation free energy in the form of the Gaussian fluctuation (GF) functional. We first test the performance of the KGK closure coupled to the reference interaction site model (RISM) integral equations on the examples of Lennard-Jones liquids, polar and nonpolar molecular solvents, including water, and aqueous solutions of simple ions. The solvation structure, solvation chemical potential, and compressibility obtained from RISM with the KGK closure favorably compare to the results of the hypernetted chain (HNC) and Kovalenko-Hirata (KH) closures, including their combination with the GF solvation free energy. We then use the KGK closure coupled to RISM to obtain the solvation structure and thermodynamics of oligomeric polyelectrolytes and drug-like compounds at a finite concentration in electrolyte solution, for which no convergence is obtained with other closures. For comparison, we calculate their solvation structure from molecular dynamics (MD) simulations. We further couple the 3D-RISM integral equation with the 3D-version of the KGK closure, and solve it for molecular mixtures as well as oligomeric polyelectrolytes and drug-like molecules in electrolyte solutions.
Collapse
Affiliation(s)
- Alexander E Kobryn
- National Institute for Nanotechnology, National Research Council Canada, 11421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada
| | | | | |
Collapse
|
32
|
Omelyan I, Kovalenko A. MTS-MD of Biomolecules Steered with 3D-RISM-KH Mean Solvation Forces Accelerated with Generalized Solvation Force Extrapolation. J Chem Theory Comput 2016; 11:1875-95. [PMID: 26574393 DOI: 10.1021/ct5010438] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
We developed a generalized solvation force extrapolation (GSFE) approach to speed up multiple time step molecular dynamics (MTS-MD) of biomolecules steered with mean solvation forces obtained from the 3D-RISM-KH molecular theory of solvation (three-dimensional reference interaction site model with the Kovalenko-Hirata closure). GSFE is based on a set of techniques including the non-Eckart-like transformation of coordinate space separately for each solute atom, extension of the force-coordinate pair basis set followed by selection of the best subset, balancing the normal equations by modified least-squares minimization of deviations, and incremental increase of outer time step in motion integration. Mean solvation forces acting on the biomolecule atoms in conformations at successive inner time steps are extrapolated using a relatively small number of best (closest) solute atomic coordinates and corresponding mean solvation forces obtained at previous outer time steps by converging the 3D-RISM-KH integral equations. The MTS-MD evolution steered with GSFE of 3D-RISM-KH mean solvation forces is efficiently stabilized with our optimized isokinetic Nosé-Hoover chain (OIN) thermostat. We validated the hybrid MTS-MD/OIN/GSFE/3D-RISM-KH integrator on solvated organic and biomolecules of different stiffness and complexity: asphaltene dimer in toluene solvent, hydrated alanine dipeptide, miniprotein 1L2Y, and protein G. The GSFE accuracy and the OIN efficiency allowed us to enlarge outer time steps up to huge values of 1-4 ps while accurately reproducing conformational properties. Quasidynamics steered with 3D-RISM-KH mean solvation forces achieves time scale compression of conformational changes coupled with solvent exchange, resulting in further significant acceleration of protein conformational sampling with respect to real time dynamics. Overall, this provided a 50- to 1000-fold effective speedup of conformational sampling for these systems, compared to conventional MD with explicit solvent. We have been able to fold the miniprotein from a fully denatured, extended state in about 60 ns of quasidynamics steered with 3D-RISM-KH mean solvation forces, compared to the average physical folding time of 4-9 μs observed in experiment.
Collapse
Affiliation(s)
- Igor Omelyan
- Department of Mechanical Engineering, University of Alberta , Mechanical Engineering Building 4-9, Edmonton, Alberta T6G 2G8, Canada.,National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada.,Institute for Condensed Matter Physics, National Academy of Sciences of Ukraine , 1 Svientsitskii Street, Lviv 79011, Ukraine
| | - Andriy Kovalenko
- Department of Mechanical Engineering, University of Alberta , Mechanical Engineering Building 4-9, Edmonton, Alberta T6G 2G8, Canada.,National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| |
Collapse
|
33
|
Peeples CA, Schreckenbach G. Implementation of the SM12 Solvation Model into ADF and Comparison with COSMO. J Chem Theory Comput 2016; 12:4033-41. [PMID: 27322858 DOI: 10.1021/acs.jctc.6b00410] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this article, an implementation of the newest iteration of the Minnesota solvation model, SM12, into the Amsterdam density functional (ADF) computational package is presented. ADF makes exclusive use of Slater-type orbitals (STO), which correctly represent the true atomic orbitals for atoms, whereas SM12 and the underlying charge model 5 (CM5) have previously only been tested on Gaussian-type orbitals (GTO). This new implementation is used to prove the basis set independence of both CM5 and SM12. A detailed comparison of the SM12 and COSMO solvation models, as implemented in ADF, is also presented. We show that this new implementation of SM12 has a mean unsigned error (MUE) of 0.68 kcal/mol for 272 molecules in water solvent, 4.10 kcal/mol MUE for 112 charged ions in water, and 0.92 kcal/mol MUE for 197 solvent calculations of various molecules. SM12 outperforms COSMO for all neutral molecules and performs as well as COSMO for cationic molecules, only falling short when anionic molecules are taken into consideration, likely due to CM5's use of Hirshfeld charges and their poor description of anionic molecules, though CM5 seems to improve upon this discrepancy.
Collapse
Affiliation(s)
- Craig A Peeples
- Department of Chemistry, University of Manitoba Winnipeg , Winnipeg, MB R3T 2N2, Canada
| | - Georg Schreckenbach
- Department of Chemistry, University of Manitoba Winnipeg , Winnipeg, MB R3T 2N2, Canada
| |
Collapse
|
34
|
|
35
|
Hölzl C, Kibies P, Imoto S, Frach R, Suladze S, Winter R, Marx D, Horinek D, Kast SM. Design principles for high–pressure force fields: Aqueous TMAO solutions from ambient to kilobar pressures. J Chem Phys 2016; 144:144104. [DOI: 10.1063/1.4944991] [Citation(s) in RCA: 67] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Christoph Hölzl
- Institut für Physikalische und Theoretische Chemie, Universität Regensburg, 93040 Regensburg, Germany
| | - Patrick Kibies
- Physikalische Chemie III, Technische Universität Dortmund, 44227 Dortmund, Germany
| | - Sho Imoto
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Roland Frach
- Physikalische Chemie III, Technische Universität Dortmund, 44227 Dortmund, Germany
| | - Saba Suladze
- Physikalische Chemie I, Technische Universität Dortmund, 44227 Dortmund, Germany
| | - Roland Winter
- Physikalische Chemie I, Technische Universität Dortmund, 44227 Dortmund, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Dominik Horinek
- Institut für Physikalische und Theoretische Chemie, Universität Regensburg, 93040 Regensburg, Germany
| | - Stefan M. Kast
- Physikalische Chemie III, Technische Universität Dortmund, 44227 Dortmund, Germany
| |
Collapse
|
36
|
Aono S, Mori T, Sakaki S. 3D-RISM-MP2 Approach to Hydration Structure of Pt(II) and Pd(II) Complexes: Unusual H-Ahead Mode vs Usual O-Ahead One. J Chem Theory Comput 2016; 12:1189-206. [PMID: 26863511 DOI: 10.1021/acs.jctc.5b01137] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Solvation of transition metal complexes with water has been one of the fundamental topics in physical and coordination chemistry. In particular, Pt(II) complexes have recently attracted considerable interest for their relation to anticancer activity in cisplatin and its analogues, yet the interaction of the water molecule and the metal center has been obscured. The challenge from a theoretical perspective remains that both the microscopic solvation effect and the dynamical electron correlation (DEC) effect have to be treated simultaneously in a reasonable manner. In this work we derive the analytical gradient for the three-dimensional reference interaction site model Møller-Plesset second order (3D-RISM-MP2) free energy. On the basis of the three-regions 3D-RISM self-consistent field (SCF) method recently proposed by us, we apply a new layer of the Z-vector method to the CP-RISM equation as well as point-charge approximation to the derivatives with respect to the density matrix elements in the RISM-CPHF equation to remarkably reduce the computational cost. This method is applied to study the interaction of H2O with the d(8) square planar transition metal complexes in aqueous solution, trans-[Pt(II)Cl2(NH3)(glycine)] (1a), [Pt(II)(NH3)4](2+) (1b), [Pt(II)(CN)4](2-) (1c), and their Pd(II) analogues 2a, 2b, and 2c, respectively, to elucidate whether the usual H2O interaction through O atom (O-ahead mode) or unusual one through H atom (H-ahead mode) is stable in these complexes. We find that the interaction energy of the coordinating water and the transition metal complex changes little when switching from gas to aqueous phase, but the solvation free energy differs remarkably between the two interaction modes, thereby affecting the relative stability of the H-ahead and O-ahead modes. Particularly, in contrast to the expectation that the O-ahead mode is preferred due to the presence of positive charges in 1b, the H-ahead mode is also found to be more stable. The O-ahead mode is found to be more stable than the H-ahead one only in 2b. The energy decomposition analysis (EDA) at the 3D-RISM-MP2 level revealed that the O-ahead mode is stabilized by the electrostatic (ES) interaction, whereas the H-ahead one is mainly stabilized by the DEC effect. The ES interaction is also responsible for the difference between the Pd(II) and Pt(II) complexes; because the electrostatic potential is more negative along the z-axis in the Pt(II) complex than in the Pd(II) one, the O-ahead mode prefers the Pd(II) complexes, whereas the H-ahead becomes predominant in the Pt(II) complexes.
Collapse
Affiliation(s)
- Shinji Aono
- Fukui Institute for Fundamental Chemistry, Kyoto University , Nishihiraki-cho, Takano, Sakyo-ku, Kyoto 606-8103, Japan
| | - Toshifumi Mori
- Institute for Molecular Science , Okazaki, Aichi 444-8585, Japan.,School of Physical Sciences, The Graduate University for Advanced Studies , Okazaki, Aichi 444-8585, Japan
| | - Shigeyoshi Sakaki
- Fukui Institute for Fundamental Chemistry, Kyoto University , Nishihiraki-cho, Takano, Sakyo-ku, Kyoto 606-8103, Japan
| |
Collapse
|
37
|
Casanova D, Gusarov S, Kovalenko A, Ziegler T. Evaluation of the SCF Combination of KS-DFT and 3D-RISM-KH; Solvation Effect on Conformational Equilibria, Tautomerization Energies, and Activation Barriers. J Chem Theory Comput 2015; 3:458-76. [PMID: 26637028 DOI: 10.1021/ct6001785] [Citation(s) in RCA: 54] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The effect of solvation on conformational equilibria, tautomerization energies, and activation barriers in simple SN2 reactions is reproduced by using the self-consistent field coupling of the Kohn-Sham density functional theory (KS-DFT) for electronic structure and the three-dimensional reference interaction site model with the closure approximation of Kovalenko and Hirata (3D-RISM-KH) for molecular solvation structure. These examples are used in order to validate the implementation of the 3D-RISM-KH method in the Amsterdam Density Functional package. The computations of the free energy difference in the trans/gauche conformational equilibrium for 1,2-dichloroethane in different solvents; the relative tautomerization free energy for cytosine, isocytosine, and guanine; and the free energy activation barrier for a CH3X-type (X = F, Cl, Br) SN2 reaction exhibit agreement with the experimental data. The method is also applied to the electronic and hydration structure of carbon single-wall nanotubes of different diameters, including the effect of water located in the inner space of the nanotubes. A comparison with continuum models of solvation (including COSMO) as well as with other more precise and computationally more expensive calculations is made to demonstrate the accuracy and predictive capability of the new KS-DFT/3D-RISM-KH method.
Collapse
Affiliation(s)
- David Casanova
- Departament de Química Inorgànica and Centre de Recerca en Química Teòrica, Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain, National Institute for Nanotechnology, National Research Council of Canada, 421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada, and Department of Chemistry, University of Calgary, Calgary, Alberta, T2N 1N4, Canada
| | - Sergey Gusarov
- Departament de Química Inorgànica and Centre de Recerca en Química Teòrica, Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain, National Institute for Nanotechnology, National Research Council of Canada, 421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada, and Department of Chemistry, University of Calgary, Calgary, Alberta, T2N 1N4, Canada
| | - Andriy Kovalenko
- Departament de Química Inorgànica and Centre de Recerca en Química Teòrica, Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain, National Institute for Nanotechnology, National Research Council of Canada, 421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada, and Department of Chemistry, University of Calgary, Calgary, Alberta, T2N 1N4, Canada
| | - Tom Ziegler
- Departament de Química Inorgànica and Centre de Recerca en Química Teòrica, Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain, National Institute for Nanotechnology, National Research Council of Canada, 421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada, and Department of Chemistry, University of Calgary, Calgary, Alberta, T2N 1N4, Canada
| |
Collapse
|
38
|
Three-dimensional reference interaction site model solvent combined with a quantum mechanical treatment of the solute. COMPUT THEOR CHEM 2015. [DOI: 10.1016/j.comptc.2015.08.003] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
39
|
Matveev A, Li B, Rösch N. Uranyl Solvation by a Three-Dimensional Reference Interaction Site Model. J Phys Chem A 2015; 119:8702-13. [PMID: 26167741 DOI: 10.1021/acs.jpca.5b03712] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We report an implementation of the three-dimensional reference interaction site model (3D RISM) that in particular addresses the treatment of the long-range Coulomb field of charged species, represented by point charges and/or a distributed charge density. A comparison of 1D and 3D results for atomic ions demonstrates a reasonable accuracy, even for a moderate size of the unit cell and a moderate grid resolution. In an application to uranyl complexes with 4-6 explicit aqua ligands and an implicit bulk solvent modeled by RISM, we show that the 3D technique is not susceptible to the deficiencies of the 1D technique exposed in our previous work [Li, Matveev, Krüger, Rösch, Comp. Theor. Chem. 2015, 1051, 151]. The 3D method eliminates the artificial superposition of explicit aqua ligands and the RISM medium and predicts essentially the same values for uranyl and uranyl-water bond lengths as a state-of-the-art polarizable continuum model. With the first solvation shell treated explicitly, the observables are nearly independent of the order of the closure relationship used when solving the set of integral equations for the various distribution functions. Furthermore, we calculated the activation barrier of water exchange with a hybrid approach that combines the 3D RISM model for the bulk aqueous solvent and a quantum mechanical description (at the level of electronic density functional theory) of uranyl interacting with explicitly represented water molecules. The calculated result agrees very well with experiment and the best theoretical estimates.
Collapse
Affiliation(s)
| | | | - Notker Rösch
- §Institute of High Performance Computing, Agency for Science, Technology and Research, 1 Fusionopolis Way, Connexis #16-16, Singapore 138632, Singapore
| |
Collapse
|
40
|
Ratkova EL, Palmer DS, Fedorov MV. Solvation thermodynamics of organic molecules by the molecular integral equation theory: approaching chemical accuracy. Chem Rev 2015; 115:6312-56. [PMID: 26073187 DOI: 10.1021/cr5000283] [Citation(s) in RCA: 135] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Ekaterina L Ratkova
- †G. A. Krestov Institute of Solution Chemistry of the Russian Academy of Sciences, Akademicheskaya Street 1, Ivanovo 153045, Russia.,‡The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig 04103, Germany
| | - David S Palmer
- ‡The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig 04103, Germany.,§Department of Chemistry, University of Strathclyde, Thomas Graham Building, 295 Cathedral Street, Glasgow, Scotland G1 1XL, United Kingdom
| | - Maxim V Fedorov
- ‡The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig 04103, Germany.,∥Department of Physics, Scottish Universities Physics Alliance (SUPA), University of Strathclyde, John Anderson Building, 107 Rottenrow East, Glasgow G4 0NG, United Kingdom
| |
Collapse
|
41
|
Huang W, Blinov N, Kovalenko A. Octanol-Water Partition Coefficient from 3D-RISM-KH Molecular Theory of Solvation with Partial Molar Volume Correction. J Phys Chem B 2015; 119:5588-97. [PMID: 25844645 DOI: 10.1021/acs.jpcb.5b01291] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The octanol-water partition coefficient is an important physical-chemical characteristic widely used to describe hydrophobic/hydrophilic properties of chemical compounds. The partition coefficient is related to the transfer free energy of a compound from water to octanol. Here, we introduce a new protocol for prediction of the partition coefficient based on the statistical-mechanical, 3D-RISM-KH molecular theory of solvation. It was shown recently that with the compound-solvent correlation functions obtained from the 3D-RISM-KH molecular theory of solvation, the free energy functional supplemented with the correction linearly related to the partial molar volume obtained from the Kirkwood-Buff/3D-RISM theory, also called the "universal correction" (UC), provides accurate prediction of the hydration free energy of small compounds, compared to explicit solvent molecular dynamics [ Palmer , D. S. ; J. Phys.: Condens. Matter 2010 , 22 , 492101 ]. Here we report that with the UC reparametrized accordingly this theory also provides an excellent agreement with the experimental data for the solvation free energy in nonpolar solvent (1-octanol) and so accurately predicts the octanol-water partition coefficient. The performance of the Kovalenko-Hirata (KH) and Gaussian fluctuation (GF) functionals of the solvation free energy, with and without UC, is tested on a large library of small compounds with diverse functional groups. The best agreement with the experimental data for octanol-water partition coefficients is obtained with the KH-UC solvation free energy functional.
Collapse
Affiliation(s)
- WenJuan Huang
- †Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta Canada.,‡National Institute for Nanotechnology, National Research Council of Canada, Edmonton, Alberta Canada
| | - Nikolay Blinov
- †Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta Canada.,‡National Institute for Nanotechnology, National Research Council of Canada, Edmonton, Alberta Canada
| | - Andriy Kovalenko
- †Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta Canada.,‡National Institute for Nanotechnology, National Research Council of Canada, Edmonton, Alberta Canada
| |
Collapse
|
42
|
Shedge SV, Wesolowski TA. Nonuniform Continuum Model for Solvatochromism Based on Frozen-Density Embedding Theory. Chemphyschem 2014; 15:3291-300. [DOI: 10.1002/cphc.201402351] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2014] [Revised: 07/09/2014] [Indexed: 11/07/2022]
|
43
|
Kovalenko A. Predictive Multiscale Modeling of Nanocellulose Based Materials and Systems. ACTA ACUST UNITED AC 2014. [DOI: 10.1088/1757-899x/64/1/012040] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
44
|
Omelyan I, Kovalenko A. Multiple time step molecular dynamics in the optimized isokinetic ensemble steered with the molecular theory of solvation: accelerating with advanced extrapolation of effective solvation forces. J Chem Phys 2014; 139:244106. [PMID: 24387356 DOI: 10.1063/1.4848716] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
We develop efficient handling of solvation forces in the multiscale method of multiple time step molecular dynamics (MTS-MD) of a biomolecule steered by the solvation free energy (effective solvation forces) obtained from the 3D-RISM-KH molecular theory of solvation (three-dimensional reference interaction site model complemented with the Kovalenko-Hirata closure approximation). To reduce the computational expenses, we calculate the effective solvation forces acting on the biomolecule by using advanced solvation force extrapolation (ASFE) at inner time steps while converging the 3D-RISM-KH integral equations only at large outer time steps. The idea of ASFE consists in developing a discrete non-Eckart rotational transformation of atomic coordinates that minimizes the distances between the atomic positions of the biomolecule at different time moments. The effective solvation forces for the biomolecule in a current conformation at an inner time step are then extrapolated in the transformed subspace of those at outer time steps by using a modified least square fit approach applied to a relatively small number of the best force-coordinate pairs. The latter are selected from an extended set collecting the effective solvation forces obtained from 3D-RISM-KH at outer time steps over a broad time interval. The MTS-MD integration with effective solvation forces obtained by converging 3D-RISM-KH at outer time steps and applying ASFE at inner time steps is stabilized by employing the optimized isokinetic Nosé-Hoover chain (OIN) ensemble. Compared to the previous extrapolation schemes used in combination with the Langevin thermostat, the ASFE approach substantially improves the accuracy of evaluation of effective solvation forces and in combination with the OIN thermostat enables a dramatic increase of outer time steps. We demonstrate on a fully flexible model of alanine dipeptide in aqueous solution that the MTS-MD/OIN/ASFE/3D-RISM-KH multiscale method of molecular dynamics steered by effective solvation forces allows huge outer time steps up to tens of picoseconds without affecting the equilibrium and conformational properties, and thus provides a 100- to 500-fold effective speedup in comparison to conventional MD with explicit solvent. With the statistical-mechanical 3D-RISM-KH account for effective solvation forces, the method provides efficient sampling of biomolecular processes with slow and/or rare solvation events such as conformational transitions of hydrated alanine dipeptide with the mean life times ranging from 30 ps up to 10 ns for "flip-flop" conformations, and is particularly beneficial for biomolecular systems with exchange and localization of solvent and ions, ligand binding, and molecular recognition.
Collapse
Affiliation(s)
- Igor Omelyan
- National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| | - Andriy Kovalenko
- National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| |
Collapse
|
45
|
State-specific solvation effect on the intramolecular charge transfer reaction in solution: A linear-response free energy TDDFT method. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.05.104] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
46
|
Yoshida N. Efficient implementation of the three-dimensional reference interaction site model method in the fragment molecular orbital method. J Chem Phys 2014; 140:214118. [DOI: 10.1063/1.4879795] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
|
47
|
Maruyama Y, Yoshida N, Tadano H, Takahashi D, Sato M, Hirata F. Massively parallel implementation of 3D-RISM calculation with volumetric 3D-FFT. J Comput Chem 2014; 35:1347-55. [PMID: 24771232 DOI: 10.1002/jcc.23619] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2014] [Revised: 03/26/2014] [Accepted: 04/06/2014] [Indexed: 01/19/2023]
Abstract
A new three-dimensional reference interaction site model (3D-RISM) program for massively parallel machines combined with the volumetric 3D fast Fourier transform (3D-FFT) was developed, and tested on the RIKEN K supercomputer. The ordinary parallel 3D-RISM program has a limitation on the number of parallelizations because of the limitations of the slab-type 3D-FFT. The volumetric 3D-FFT relieves this limitation drastically. We tested the 3D-RISM calculation on the large and fine calculation cell (2048(3) grid points) on 16,384 nodes, each having eight CPU cores. The new 3D-RISM program achieved excellent scalability to the parallelization, running on the RIKEN K supercomputer. As a benchmark application, we employed the program, combined with molecular dynamics simulation, to analyze the oligomerization process of chymotrypsin Inhibitor 2 mutant. The results demonstrate that the massive parallel 3D-RISM program is effective to analyze the hydration properties of the large biomolecular systems.
Collapse
Affiliation(s)
- Yutaka Maruyama
- Department of Physics, Faculty of Science and Technology, Keio University, Yokohama, Kanagawa, 223-8522, Japan
| | | | | | | | | | | |
Collapse
|
48
|
Minezawa N. Excited-state free energy surfaces in solution: time-dependent density functional theory∕reference interaction site model self-consistent field method. J Chem Phys 2014; 138:244101. [PMID: 23822221 DOI: 10.1063/1.4811201] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Constructing free energy surfaces for electronically excited states is a first step toward the understanding of photochemical processes in solution. For that purpose, the analytic free energy gradient is derived and implemented for the linear-response time-dependent density functional theory combined with the reference interaction site model self-consistent field method. The proposed method is applied to study (1) the fluorescence spectra of aqueous acetone and (2) the excited-state intramolecular proton transfer reaction of ortho-hydroxybenzaldehyde in an acetonitrile solution.
Collapse
Affiliation(s)
- Noriyuki Minezawa
- Fukui Institute for Fundamental Chemistry, Kyoto University, Sakyo-ku, Kyoto 606-8103, Japan
| |
Collapse
|
49
|
Fujishige S, Kawashima Y, Yoshida N, Nakano H. Three-Dimensional Reference Interaction Site Model Self-Consistent Field Study of the Electronic Structure of [Cr(H2O)6]3+ in Aqueous Solution. J Phys Chem A 2013; 117:8314-22. [DOI: 10.1021/jp405876g] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Shinya Fujishige
- Department of Chemistry, Graduate
School of Sciences, Kyushu University,
Fukuoka 812-8581, Japan
| | - Yukio Kawashima
- Department of Chemistry, Graduate
School of Sciences, Kyushu University,
Fukuoka 812-8581, Japan
| | - Norio Yoshida
- Department of Chemistry, Graduate
School of Sciences, Kyushu University,
Fukuoka 812-8581, Japan
| | - Haruyuki Nakano
- Department of Chemistry, Graduate
School of Sciences, Kyushu University,
Fukuoka 812-8581, Japan
| |
Collapse
|
50
|
Multiscale modeling of solvation in chemical and biological nanosystems and in nanoporous materials. PURE APPL CHEM 2013. [DOI: 10.1351/pac-con-12-06-03] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Statistical–mechanical, 3D-RISM-KH molecular theory of solvation (3D reference interaction site model with the Kovalenko–Hirata closure) is promising as an essential part of multiscale methodology for chemical and biomolecular nanosystems in solution. 3D-RISM-KH explains the molecular mechanisms of self-assembly and conformational stability of synthetic organic rosette nanotubes (RNTs), aggregation of prion proteins and β-sheet amyloid oligomers, protein-ligand binding, and function-related solvation properties of complexes as large as the Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC) and GroEL/ES chaperone. Molecular mechanics/Poisson–Boltzmann (generalized Born) surface area [MM/PB(GB)SA] post-processing of molecular dynamics (MD) trajectories involving SA empirical nonpolar terms is replaced with MM/3D-RISM-KH statistical–mechanical evaluation of the solvation thermodynamics. 3D-RISM-KH has been coupled with multiple time-step (MTS) MD of the solute biomolecule driven by effective solvation forces, which are obtained analytically by converging the 3D-RISM-KH integral equations at outer time-steps and are calculated in between by using solvation force coordinate extrapolation (SFCE) in the subspace of previous solutions to 3D-RISM-KH. The procedure is stabilized by the optimized isokinetic Nosé–Hoover (OIN) chain thermostatting, which enables gigantic outer time-steps up to picoseconds to accurately calculate equilibrium properties. The multiscale OIN/SFCE/3D-RISM-KH algorithm is implemented in the Amber package and illustrated on a fully flexible model of alanine dipeptide in aqueous solution, exhibiting the computational rate of solvent sampling 20 times faster than standard MD with explicit solvent. Further substantial acceleration can be achieved with 3D-RISM-KH efficiently sampling essential events with rare statistics such as exchange and localization of solvent, ions, and ligands at binding sites and pockets of the biomolecule. 3D-RISM-KH was coupled with ab initio complete active space self-consistent field (CASSCF) and orbital-free embedding (OFE) Kohn–Sham (KS) density functional theory (DFT) quantum chemistry methods in an SCF description of electronic structure, optimized geometry, and chemical reactions in solution. The (OFE)KS-DFT/3D-RISM-KH multi-scale method is implemented in the Amsterdam Density Functional (ADF) package and extensively validated against experiment for solvation thermochemistry, photochemistry, conformational equilibria, and activation barriers of various nanosystems in solvents and ionic liquids (ILs). Finally, the replica RISM-KH-VM molecular theory for the solvation structure, thermodynamics, and electrochemistry of electrolyte solutions sorbed in nanoporous materials reveals the molecular mechanisms of sorption and supercapacitance in nanoporous carbon electrodes, which is drastically different from a planar electrical double layer.
Collapse
|