1
|
Maste S, Sharma B, Pongratz T, Grabe B, Hiller W, Erlach MB, Kremer W, Kalbitzer HR, Marx D, Kast SM. The accuracy limit of chemical shift predictions for species in aqueous solution. Phys Chem Chem Phys 2024; 26:6386-6395. [PMID: 38315169 DOI: 10.1039/d3cp05471c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2024]
Abstract
Interpreting NMR experiments benefits from first-principles predictions of chemical shifts. Reaching the accuracy limit of theory is relevant for unambiguous structural analysis and dissecting theoretical approximations. Since accurate chemical shift measurements are based on using internal reference compounds such as trimethylsilylpropanesulfonate (DSS), a detailed comparison of experimental with theoretical data requires simultaneous consideration of both target and reference species ensembles in the same solvent environment. Here we show that ab initio molecular dynamics simulations to generate liquid-state ensembles of target and reference compounds, including explicitly their short-range solvation environments and combined with quantum-mechanical solvation models, allows for predicting highly accurate 1H (∼0.1-0.5 ppm) and aliphatic 13C (∼1.5 ppm) chemical shifts for aqueous solutions of the model compounds trimethylamine N-oxide (TMAO) and N-methylacetamide (NMA), referenced to DSS without any system-specific adjustments. This encompasses the two peptide bond conformations of NMA identified by NMR. The results are used to derive a general-purpose guideline set for predictive NMR chemical shift calculations of NMA in the liquid state and to identify artifacts of force field models. Accurate predictions are only obtained if a sufficient number of explicit water molecules is included in the quantum-mechanical calculations, disproving a purely electrostatic model of the solvent effect on chemical shifts.
Collapse
Affiliation(s)
- Stefan Maste
- Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44227 Dortmund, Germany.
| | - Bikramjit Sharma
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany.
| | - Tim Pongratz
- Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44227 Dortmund, Germany.
| | - Bastian Grabe
- Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44227 Dortmund, Germany.
| | - Wolf Hiller
- Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44227 Dortmund, Germany.
| | - Markus Beck Erlach
- Fakultät für Biologie und Vorklinische Medizin, Universität Regensburg, 93040 Regensburg, Germany
| | - Werner Kremer
- Fakultät für Biologie und Vorklinische Medizin, Universität Regensburg, 93040 Regensburg, Germany
| | - Hans Robert Kalbitzer
- Fakultät für Biologie und Vorklinische Medizin, Universität Regensburg, 93040 Regensburg, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany.
| | - Stefan M Kast
- Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44227 Dortmund, Germany.
| |
Collapse
|
2
|
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
|
3
|
Reimann M, Kaupp M. Spin-State Splittings in 3d Transition-Metal Complexes Revisited: Toward a Reliable Theory Benchmark. J Chem Theory Comput 2023; 19:97-108. [PMID: 36576816 DOI: 10.1021/acs.jctc.2c00925] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
A new composite method for the calculation of spin-crossover energies in 3d transition-metal complexes based on multireference methods is presented. The method reduces to MRCISD+Q at the complete-basis-set (CBS) level for atomic ions, for which it gives excitation energies with a mean absolute error of only ca. 0.01 eV. For molecular complexes, the CASPT2+δMRCI composite approach corresponds to a CASPT2/CBS calculation augmented by a high-level MRCISD+Q-CASPT2 correction with a smaller ligand basis set. For a set of [Fe(He)6]n+ test complexes, the approach reproduces full MRCISD+Q/CBS results to within better than 0.04 eV, without depending on any arbitrary IPEA shifts. The high-quality CASPT2+δMRCI method has then been applied to a series of 3d transition-metal hexaqua complexes in aqueous solution, augmented by an elaborate 3D-RISM-SCF solvent treatment of the underlying structures. It provides unprecedented agreement with experiment for the lowest-lying vertical spin-flip excitation energies, except for the Fe3+ system. Closer examination of the latter case provides strong evidence that the observed lowest-energy excitation at 1.56 eV, which has been used frequently for evaluating quantum-chemical methods, does not arise from the iron(III) hexaqua complex in solution, but from its singly deprotonated counterpart, [Fe(H2O)5OH]2+.
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
|
4
|
Nakatani K, Teshigawara S, Tanahashi Y, Kasahara K, Higashi M, Sato H. Solvation in nitration of benzene and the valence electronic structure of the Wheland intermediate. Phys Chem Chem Phys 2022; 24:16453-16461. [PMID: 35647764 DOI: 10.1039/d2cp01699k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Nitration of benzene is a representative aromatic substitution reaction related to the σ-complex (arenium ion or "Wheland" intermediate) concept. This reaction is typically carried out in a mixed acid solution to generate nitronium ions, and how solvent molecules play roles in the reaction has been of great interest. Here we will shed new light on the reaction, namely the electronic structure and the microscopic insights of the solvation, which have been rarely discussed so far. We studied this process using the reference interaction site model-self consistent field with constrained spatial electron density distribution (RISM-SCF-cSED) method, considering sulfuric acid or water molecules as a solvent. In this method, the electronic structure of the solute and the solvation structure are self-consistently determined based on quantum chemistry and statistical mechanics of molecular liquids. The solvation free energy surfaces in solution and solvation structures were verified. In the bond formation process of benzene and nitronium ions, the solvation structure by sulfuric acid molecules drastically changes and the solvation effect on the free energy is quite large. We revealed largely contributing resonance structures in the π-electron system of the σ-complex in gas and solution phases by analysing the valence electronic structures.
Collapse
Affiliation(s)
- Kaho Nakatani
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Nishikyo-ku, Kyoto 615-8510, Japan.
| | - Sho Teshigawara
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Nishikyo-ku, Kyoto 615-8510, Japan.
| | - Yuta Tanahashi
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Nishikyo-ku, Kyoto 615-8510, Japan.
| | - Kento Kasahara
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Nishikyo-ku, Kyoto 615-8510, Japan.
| | - Masahiro Higashi
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Nishikyo-ku, Kyoto 615-8510, Japan. .,Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Nishikyo-ku, Kyoto 615-8520, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Nishikyo-ku, Kyoto 615-8510, Japan. .,Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Nishikyo-ku, Kyoto 615-8520, Japan.,Fukui Institute for Fundamental Chemistry, Kyoto University, Takano Nishihiraki-cho 34-4, Sakyo-ku, Kyoto 606-8103, Japan
| |
Collapse
|
5
|
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
|
6
|
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
|
7
|
Sharma B, Tran VA, Pongratz T, Galazzo L, Zhurko I, Bordignon E, Kast SM, Neese F, Marx D. A Joint Venture of Ab Initio Molecular Dynamics, Coupled Cluster Electronic Structure Methods, and Liquid-State Theory to Compute Accurate Isotropic Hyperfine Constants of Nitroxide Probes in Water. J Chem Theory Comput 2021; 17:6366-6386. [PMID: 34516119 PMCID: PMC8515807 DOI: 10.1021/acs.jctc.1c00582] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2021] [Indexed: 01/11/2023]
Abstract
The isotropic hyperfine coupling constant (HFCC, Aiso) of a pH-sensitive spin probe in a solution, HMI (2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water, is computed using an ensemble of state-of-the-art computational techniques and is gauged against X-band continuous wave electron paramagnetic resonance (EPR) measurement spectra at room temperature. Fundamentally, the investigation aims to delineate the cutting edge of current first-principles-based calculations of EPR parameters in aqueous solutions based on using rigorous statistical mechanics combined with correlated electronic structure techniques. In particular, the impact of solvation is described by exploiting fully atomistic, RISM integral equation, and implicit solvation approaches as offered by ab initio molecular dynamics (AIMD) of the periodic bulk solution (using the spin-polarized revPBE0-D3 hybrid functional), embedded cluster reference interaction site model integral equation theory (EC-RISM), and polarizable continuum embedding (using CPCM) of microsolvated complexes, respectively. HFCCs are obtained from efficient coupled cluster calculations (using open-shell DLPNO-CCSD theory) as well as from hybrid density functional theory (using revPBE0-D3). Re-solvation of "vertically desolvated" spin probe configuration snapshots by EC-RISM embedding is shown to provide significantly improved results compared to CPCM since only the former captures the inherent structural heterogeneity of the solvent close to the spin probe. The average values of the Aiso parameter obtained based on configurational statistics using explicit water within AIMD and from EC-RISM solvation are found to be satisfactorily close. Using either such explicit or RISM solvation in conjunction with DLPNO-CCSD calculations of the HFCCs provides an average Aiso parameter for HMI in aqueous solution at 300 K and 1 bar that is in good agreement with the experimentally determined one. The developed computational strategy is general in the sense that it can be readily applied to other spin probes of similar molecular complexity, to aqueous solutions beyond ambient conditions, as well as to other solvents in the longer run.
Collapse
Affiliation(s)
- Bikramjit Sharma
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44780 Bochum, Germany
| | - Van Anh Tran
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Tim Pongratz
- Physikalische
Chemie III, Technische Universität
Dortmund, Otto-Hahn-Str. 4a, 44227 Dortmund, Germany
| | - Laura Galazzo
- Faculty
of Chemistry and Biochemistry, Ruhr University
Bochum, 44780 Bochum, Germany
| | - Irina Zhurko
- Laboratory
of Nitrogen Compounds, N.N. Vorozhtsov Novosibirsk Institute of Organic
Chemistry, NIOCH SB RAS, 9 Lavrentiev Avenue, 630090, Novosibirsk, Russia
| | - Enrica Bordignon
- Faculty
of Chemistry and Biochemistry, Ruhr University
Bochum, 44780 Bochum, Germany
| | - Stefan M. Kast
- Physikalische
Chemie III, Technische Universität
Dortmund, Otto-Hahn-Str. 4a, 44227 Dortmund, Germany
| | - Frank Neese
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Dominik Marx
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44780 Bochum, Germany
| |
Collapse
|
8
|
Tielker N, Güssregen S, Kast SM. SAMPL7 physical property prediction from EC-RISM theory. J Comput Aided Mol Des 2021; 35:933-941. [PMID: 34278539 PMCID: PMC8367877 DOI: 10.1007/s10822-021-00410-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2021] [Accepted: 07/05/2021] [Indexed: 01/08/2023]
Abstract
Inspired by the successful application of the embedded cluster reference interaction site model (EC-RISM), a combination of quantum-mechanical calculations with three-dimensional RISM theory to predict Gibbs energies of species in solution within the SAMPL6.1 (acidity constants, pKa) and SAMPL6.2 (octanol-water partition coefficients, log P) the methodology was applied to the recent SAMPL7 physical property challenge on aqueous pKa and octanol-water log P values. Not part of the challenge but provided by the organizers, we also computed distribution coefficients log D7.4 from predicted pKa and log P data. While macroscopic pKa predictions compared very favorably with experimental data (root mean square error, RMSE 0.72 pK units), the performance of the log P model (RMSE 1.84) fell behind expectations from the SAMPL6.2 challenge, leading to reasonable log D7.4 predictions (RMSE 1.69) from combining the independent calculations. In the post-submission phase, conformations generated by different methodology yielded results that did not significantly improve the original predictions. While overall satisfactory compared to previous log D challenges, the predicted data suggest that further effort is needed for optimizing the robustness of the partition coefficient model within EC-RISM calculations and for shaping the agreement between experimental conditions and the corresponding model description.
Collapse
Affiliation(s)
- Nicolas Tielker
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Stefan Güssregen
- Sanofi-Aventis Deutschland GmbH, R&D Integrated Drug Discovery, 65926, Frankfurt am Main, Germany
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany.
| |
Collapse
|
9
|
Borgis D, Luukkonen S, Belloni L, Jeanmairet G. Accurate prediction of hydration free energies and solvation structures using molecular density functional theory with a simple bridge functional. J Chem Phys 2021; 155:024117. [PMID: 34266282 DOI: 10.1063/5.0057506] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper assesses the ability of molecular density functional theory to predict efficiently and accurately the hydration free energies of molecular solutes and the surrounding microscopic water structure. A wide range of solutes were investigated, including hydrophobes, water as a solute, and the FreeSolv database containing 642 drug-like molecules having a variety of shapes and sizes. The usual second-order approximation of the theory is corrected by a third-order, angular-independent bridge functional. The overall functional is parameter-free in the sense that the only inputs are bulk water properties, independent of the solutes considered. These inputs are the direct correlation function, compressibility, liquid-gas surface tension, and excess chemical potential of the solvent. Compared to molecular simulations with the same force field and the same fixed solute geometries, the present theory is shown to describe accurately the solvation free energy and structure of both hydrophobic and hydrophilic solutes. Overall, the method yields a precision of order 0.5 kBT for the hydration free energies of the FreeSolv database, with a computer speedup of 3 orders of magnitude. The theory remains to be improved for a better description of the H-bonding structure and the hydration free energy of charged solutes.
Collapse
Affiliation(s)
- Daniel Borgis
- Maison de la Simulation, USR 3441 CNRS-CEA-Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Sohvi Luukkonen
- Maison de la Simulation, USR 3441 CNRS-CEA-Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Luc Belloni
- Universié Paris-Saclay, CEA, CNRS, NIMBE, 91191 Gif-sur-Yvette, France
| | - Guillaume Jeanmairet
- Sorbonne Université, CNRS, Physico-Chimie des Électrolytes et Nanosystèmes Interfaciaux, PHENIX, F-75005 Paris, France
| |
Collapse
|
10
|
Viayna A, Pinheiro S, Curutchet C, Luque FJ, Zamora WJ. Prediction of n-octanol/water partition coefficients and acidity constants (pK a) in the SAMPL7 blind challenge with the IEFPCM-MST model. J Comput Aided Mol Des 2021; 35:803-811. [PMID: 34244905 PMCID: PMC8295120 DOI: 10.1007/s10822-021-00394-6] [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: 02/23/2021] [Accepted: 06/04/2021] [Indexed: 12/17/2022]
Abstract
Within the scope of SAMPL7 challenge for predicting physical properties, the Integral Equation Formalism of the Miertus-Scrocco-Tomasi (IEFPCM/MST) continuum solvation model has been used for the blind prediction of n-octanol/water partition coefficients and acidity constants of a set of 22 and 20 sulfonamide-containing compounds, respectively. The log P and pKa were computed using the B3LPYP/6-31G(d) parametrized version of the IEFPCM/MST model. The performance of our method for partition coefficients yielded a root-mean square error of 1.03 (log P units), placing this method among the most accurate theoretical approaches in the comparison with both globally (rank 8th) and physical (rank 2nd) methods. On the other hand, the deviation between predicted and experimental pKa values was 1.32 log units, obtaining the second best-ranked submission. Though this highlights the reliability of the IEFPCM/MST model for predicting the partitioning and the acid dissociation constant of drug-like compounds compound, the results are discussed to identify potential weaknesses and improve the performance of the method.
Collapse
Affiliation(s)
- Antonio Viayna
- Department of Nutrition, Food Sciences and Gastronomy, Faculty of Pharmacy and Food Sciences, Institute of Biomedicine (IBUB), and Institute of Theoretical and Computational Chemistry (IQTC-UB), University of Barcelona (UB), Avda. Prat de La Riba, 171, 08921, Santa Coloma de Gramenet, Spain.
| | - Silvana Pinheiro
- Institute of Exact and Natural Sciences, Federal University of Pará, Belém, Pará, 66075-110, Brazil
| | - Carles Curutchet
- Department of Pharmacy and Pharmaceutical Technology and Physical Chemistry, Faculty of Pharmacy and Food Sciences, and Institute of Theoretical and Computational Chemistry (IQTC-UB), University of Barcelona, Av. de Joan XXIII, 27-31, 08028, Barcelona, Spain
| | - F Javier Luque
- Department of Nutrition, Food Sciences and Gastronomy, Faculty of Pharmacy and Food Sciences, Institute of Biomedicine (IBUB), and Institute of Theoretical and Computational Chemistry (IQTC-UB), University of Barcelona (UB), Avda. Prat de La Riba, 171, 08921, Santa Coloma de Gramenet, Spain
| | - William J Zamora
- School of Chemistry and Faculty of Pharmacy, University of Costa Rica, San Pedro, San José, Costa Rica.,Advanced Computing Lab (CNCA), National High Technology Center (CeNAT), Pavas, San José, Costa Rica
| |
Collapse
|
11
|
Tielker N, Eberlein L, Hessler G, Schmidt KF, Güssregen S, Kast SM. Quantum-mechanical property prediction of solvated drug molecules: what have we learned from a decade of SAMPL blind prediction challenges? J Comput Aided Mol Des 2021; 35:453-472. [PMID: 33079358 PMCID: PMC8018924 DOI: 10.1007/s10822-020-00347-5] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2020] [Accepted: 09/26/2020] [Indexed: 01/26/2023]
Abstract
Joint academic-industrial projects supporting drug discovery are frequently pursued to deploy and benchmark cutting-edge methodical developments from academia in a real-world industrial environment at different scales. The dimensionality of tasks ranges from small molecule physicochemical property assessment over protein-ligand interaction up to statistical analyses of biological data. This way, method development and usability both benefit from insights gained at both ends, when predictiveness and readiness of novel approaches are confirmed, but the pharmaceutical drug makers get early access to novel tools for the quality of drug products and benefit of patients. Quantum-mechanical and simulation methods particularly fall into this group of methods, as they require skills and expense in their development but also significant resources in their application, thus are comparatively slowly dripping into the realm of industrial use. Nevertheless, these physics-based methods are becoming more and more useful. Starting with a general overview of these and in particular quantum-mechanical methods for drug discovery we review a decade-long and ongoing collaboration between Sanofi and the Kast group focused on the application of the embedded cluster reference interaction site model (EC-RISM), a solvation model for quantum chemistry, to study small molecule chemistry in the context of joint participation in several SAMPL (Statistical Assessment of Modeling of Proteins and Ligands) blind prediction challenges. Starting with early application to tautomer equilibria in water (SAMPL2) the methodology was further developed to allow for challenge contributions related to predictions of distribution coefficients (SAMPL5) and acidity constants (SAMPL6) over the years. Particular emphasis is put on a frequently overlooked aspect of measuring the quality of models, namely the retrospective analysis of earlier datasets and predictions in light of more recent and advanced developments. We therefore demonstrate the performance of the current methodical state of the art as developed and optimized for the SAMPL6 pKa and octanol-water log P challenges when re-applied to the earlier SAMPL5 cyclohexane-water log D and SAMPL2 tautomer equilibria datasets. Systematic improvement is not consistently found throughout despite the similarity of the problem class, i.e. protonation reactions and phase distribution. Hence, it is possible to learn about hidden bias in model assessment, as results derived from more elaborate methods do not necessarily improve quantitative agreement. This indicates the role of chance or coincidence for model development on the one hand which allows for the identification of systematic error and opportunities toward improvement and reveals possible sources of experimental uncertainty on the other. These insights are particularly useful for further academia-industry collaborations, as both partners are then enabled to optimize both the computational and experimental settings for data generation.
Collapse
Affiliation(s)
- Nicolas Tielker
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Lukas Eberlein
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Gerhard Hessler
- R&D Integrated Drug Discovery, Sanofi-Aventis Deutschland GmbH, 65926, Frankfurt am Main, Germany
| | - K Friedemann Schmidt
- R&D Preclinical Safety, Sanofi-Aventis Deutschland GmbH, 65926, Frankfurt am Main, Germany
| | - Stefan Güssregen
- R&D Integrated Drug Discovery, Sanofi-Aventis Deutschland GmbH, 65926, Frankfurt am Main, Germany.
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany.
| |
Collapse
|
12
|
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
|
13
|
Karbalaei Khani S, Marefat Khah A, Hättig C. Comparison of Reaction Field Schemes for Coupling Continuum Solvation Models with Wave Function Methods for Excitation Energies. J Chem Theory Comput 2020; 16:4554-4564. [DOI: 10.1021/acs.jctc.0c00316] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
| | | | - Christof Hättig
- Arbeitsgruppe Quantenchemie, Ruhr-Universität, Bochum 44780, Germany
| |
Collapse
|
14
|
Işık M, Bergazin TD, Fox T, Rizzi A, Chodera JD, Mobley DL. Assessing the accuracy of octanol-water partition coefficient predictions in the SAMPL6 Part II log P Challenge. J Comput Aided Mol Des 2020; 34:335-370. [PMID: 32107702 PMCID: PMC7138020 DOI: 10.1007/s10822-020-00295-0] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Accepted: 01/24/2020] [Indexed: 12/12/2022]
Abstract
The SAMPL Challenges aim to focus the biomolecular and physical modeling community on issues that limit the accuracy of predictive modeling of protein-ligand binding for rational drug design. In the SAMPL5 log D Challenge, designed to benchmark the accuracy of methods for predicting drug-like small molecule transfer free energies from aqueous to nonpolar phases, participants found it difficult to make accurate predictions due to the complexity of protonation state issues. In the SAMPL6 log P Challenge, we asked participants to make blind predictions of the octanol-water partition coefficients of neutral species of 11 compounds and assessed how well these methods performed absent the complication of protonation state effects. This challenge builds on the SAMPL6 p[Formula: see text] Challenge, which asked participants to predict p[Formula: see text] values of a superset of the compounds considered in this log P challenge. Blind prediction sets of 91 prediction methods were collected from 27 research groups, spanning a variety of quantum mechanics (QM) or molecular mechanics (MM)-based physical methods, knowledge-based empirical methods, and mixed approaches. There was a 50% increase in the number of participating groups and a 20% increase in the number of submissions compared to the SAMPL5 log D Challenge. Overall, the accuracy of octanol-water log P predictions in SAMPL6 Challenge was higher than cyclohexane-water log D predictions in SAMPL5, likely because modeling only the neutral species was necessary for log P and several categories of method benefited from the vast amounts of experimental octanol-water log P data. There were many highly accurate methods: 10 diverse methods achieved RMSE less than 0.5 log P units. These included QM-based methods, empirical methods, and mixed methods with physical modeling supported with empirical corrections. A comparison of physical modeling methods showed that QM-based methods outperformed MM-based methods. The average RMSE of the most accurate five MM-based, QM-based, empirical, and mixed approach methods based on RMSE were 0.92 ± 0.13, 0.48 ± 0.06, 0.47 ± 0.05, and 0.50 ± 0.06, respectively.
Collapse
Affiliation(s)
- Mehtap Işık
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA.
- Tri-Institutional PhD Program in Chemical Biology, Weill Cornell Graduate School of Medical Sciences, Cornell University, New York, NY, 10065, USA.
| | | | - Thomas Fox
- Computational Chemistry, Medicinal Chemistry, Boehringer Ingelheim Pharma GmbH & Co KG, 88397, Biberach, Germany
| | - Andrea Rizzi
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA
- Tri-Institutional Training Program in Computational Biology and Medicine, New York, NY, 10065, USA
| | - John D Chodera
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA
| | - David L Mobley
- Department of Pharmaceutical Sciences, University of California, Irvine, CA, 92697, USA
- Department of Chemistry, University of California, Irvine, CA, 92697, USA
| |
Collapse
|
15
|
Eberlein L, Beierlein FR, van Eikema Hommes NJR, Radadiya A, Heil J, Benner SA, Clark T, Kast SM, Richards NGJ. Tautomeric Equilibria of Nucleobases in the Hachimoji Expanded Genetic Alphabet. J Chem Theory Comput 2020; 16:2766-2777. [PMID: 32125859 DOI: 10.1021/acs.jctc.9b01079] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Evolution has yielded biopolymers that are constructed from exactly four building blocks and are able to support Darwinian evolution. Synthetic biology aims to extend this alphabet, and we recently showed that 8-letter (hachimoji) DNA can support rule-based information encoding. One source of replicative error in non-natural DNA-like systems, however, is the occurrence of alternative tautomeric forms, which pair differently. Unfortunately, little is known about how structural modifications impact free-energy differences between tautomers of the non-natural nucleobases used in the hachimoji expanded genetic alphabet. Determining experimental tautomer ratios is technically difficult, and so, strategies for improving hachimoji DNA replication efficiency will benefit from accurate computational predictions of equilibrium tautomeric ratios. We now report that high-level quantum-chemical calculations in aqueous solution by the embedded cluster reference interaction site model, benchmarked against free-energy molecular simulations for solvation thermodynamics, provide useful quantitative information on the tautomer ratios of both Watson-Crick and hachimoji nucleobases. In agreement with previous computational studies, all four Watson-Crick nucleobases adopt essentially only one tautomer in water. This is not the case, however, for non-natural nucleobases and their analogues. For example, although the enols of isoguanine and a series of related purines are not populated in water, these heterocycles possess N1-H and N3-H keto tautomers that are similar in energy, thereby adversely impacting accurate nucleobase pairing. These robust computational strategies offer a firm basis for improving experimental measurements of tautomeric ratios, which are currently limited to studying molecules that exist only as two tautomers in solution.
Collapse
Affiliation(s)
- Lukas Eberlein
- Physikalische Chemie III, Technische Universität Dortmund, Dortmund 44227, Germany
| | - Frank R Beierlein
- Computer-Chemistry-Centre and Interdisciplinary Centre for Molecular Materials, Department of Chemistry & Pharmacy, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen 91054, Germany
| | - Nico J R van Eikema Hommes
- Computer-Chemistry-Centre and Interdisciplinary Centre for Molecular Materials, Department of Chemistry & Pharmacy, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen 91054, Germany
| | - Ashish Radadiya
- School of Chemistry, Cardiff University, Cardiff CF10 3AT, U.K
| | - Jochen Heil
- Physikalische Chemie III, Technische Universität Dortmund, Dortmund 44227, Germany
| | - Steven A Benner
- Foundation for Applied Molecular Evolution, Alachua, Florida 32615, United States
| | - Timothy Clark
- Computer-Chemistry-Centre and Interdisciplinary Centre for Molecular Materials, Department of Chemistry & Pharmacy, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen 91054, Germany
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, Dortmund 44227, Germany
| | - Nigel G J Richards
- School of Chemistry, Cardiff University, Cardiff CF10 3AT, U.K.,Foundation for Applied Molecular Evolution, Alachua, Florida 32615, United States
| |
Collapse
|
16
|
Pongratz T, Kibies P, Eberlein L, Tielker N, Hölzl C, Imoto S, Beck Erlach M, Kurrmann S, Schummel PH, Hofmann M, Reiser O, Winter R, Kremer W, Kalbitzer HR, Marx D, Horinek D, Kast SM. Pressure-dependent electronic structure calculations using integral equation-based solvation models. Biophys Chem 2020; 257:106258. [DOI: 10.1016/j.bpc.2019.106258] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2019] [Accepted: 08/25/2019] [Indexed: 12/18/2022]
|
17
|
Chuev GN, Fedotova MV, Valiev M. Chemical bond effects in classical site density functional theory of inhomogeneous molecular liquids. J Chem Phys 2020; 152:041101. [PMID: 32007044 DOI: 10.1063/1.5139619] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
Intra-molecular interactions or chemical bonds represent one of the main distinguishing characteristics of molecular fluids. Development of accurate and practical methods to treat these effects is one of the long standing problems in classical site density functional theory (SDFT). One particular instance when these issues become particularly severe is the case of classical interaction potentials with auxiliary sites or dummy atoms. In this situation, current SDFT implementations, such as the three-dimensional reference interaction site model, lead to nonphysical results. We re-examine this issue in this work using our recent reformulation of SDFT (Valiev and Chuev, J. Stat. Mech.: Theory Exp. 2018, 093201). We put forward a simple practical solution to this problem and illustrate its utility for the case of spherical solutes in diatomic liquids.
Collapse
Affiliation(s)
- Gennady N Chuev
- Institute of Theoretical and Experimental Biophysics, Russian Academy of Science, Pushchino, Moscow Region 142290, Russia
| | - Marina V Fedotova
- G.A. Krestov Institute of Solution Chemistry, Russian Academy of Sciences, Akademicheskaya St., 1, 153045 Ivanovo, Russia
| | - Marat Valiev
- Molecular Sciences Software Group, Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| |
Collapse
|
18
|
Tielker N, Tomazic D, Eberlein L, Güssregen S, Kast SM. The SAMPL6 challenge on predicting octanol-water partition coefficients from EC-RISM theory. J Comput Aided Mol Des 2020; 34:453-461. [PMID: 31981015 PMCID: PMC7125249 DOI: 10.1007/s10822-020-00283-4] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2019] [Accepted: 01/08/2020] [Indexed: 12/14/2022]
Abstract
Results are reported for octanol–water partition coefficients (log P) of the neutral states of drug-like molecules provided during the SAMPL6 (Statistical Assessment of Modeling of Proteins and Ligands) blind prediction challenge from applying the “embedded cluster reference interaction site model” (EC-RISM) as a solvation model for quantum-chemical calculations. Following the strategy outlined during earlier SAMPL challenges we first train 1- and 2-parameter water-free (“dry”) and water-saturated (“wet”) models for n-octanol solvation Gibbs energies with respect to experimental values from the “Minnesota Solvation Database” (MNSOL), yielding a root mean square error (RMSE) of 1.5 kcal mol−1 for the best-performing 2-parameter wet model, while the optimal water model developed for the pKa part of the SAMPL6 challenge is kept unchanged (RMSE 1.6 kcal mol−1 for neutral compounds from a model trained on both neutral and ionic species). Applying these models to the blind prediction set yields a log P RMSE of less than 0.5 for our best model (2-parameters, wet). Further analysis of our results reveals that a single compound is responsible for most of the error, SM15, without which the RMSE drops to 0.2. Since this is the only compound in the challenge dataset with a hydroxyl group we investigate other alcohols for which Gibbs energy of solvation data for both water and n-octanol are available in the MNSOL database to demonstrate a systematic cause of error and to discuss strategies for improvement.
Collapse
Affiliation(s)
- Nicolas Tielker
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Daniel Tomazic
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Lukas Eberlein
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Stefan Güssregen
- Sanofi-Aventis Deutschland GmbH, R&D Integrated Drug Discovery, 65926, Frankfurt am Main, Germany
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany.
| |
Collapse
|
19
|
Zamora WJ, Pinheiro S, German K, Ràfols C, Curutchet C, Luque FJ. Prediction of the n-octanol/water partition coefficients in the SAMPL6 blind challenge from MST continuum solvation calculations. J Comput Aided Mol Des 2019; 34:443-451. [PMID: 31776809 DOI: 10.1007/s10822-019-00262-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2019] [Accepted: 11/21/2019] [Indexed: 12/20/2022]
Abstract
The IEFPCM/MST continuum solvation model is used for the blind prediction of n-octanol/water partition of a set of 11 fragment-like small molecules within the SAMPL6 Part II Partition Coefficient Challenge. The partition coefficient of the neutral species (log P) was determined using an extended parametrization of the B3LYP/6-31G(d) version of the Miertus-Scrocco-Tomasi continuum solvation model in n-octanol. Comparison with the experimental data provided for partition coefficients yielded a root-mean square error (rmse) of 0.78 (log P units), which agrees with the accuracy reported for our method (rmse = 0.80) for nitrogen-containing heterocyclic compounds. Out of the 91 sets of log P values submitted by the participants, our submission is within those with an rmse < 1 and among the four best ranked physical methods. The largest errors involve three compounds: two with the largest positive deviations (SM13 and SM08), and one with the largest negative deviations (SM15). Here we report the potentiometric determination of the log P for SM13, leading to a value of 3.62 ± 0.02, which is in better agreement with most empirical predictions than the experimental value reported in SAMPL6. In addition, further inclusion of several conformations for SM08 significantly improved our results. Inclusion of these refinements led to an overall error of 0.51 (log P units), which supports the reliability of the IEFPCM/MST model for predicting the partitioning of neutral compounds.
Collapse
Affiliation(s)
- William J Zamora
- Department of Nutrition, Food Science and Gastronomy, Faculty of Pharmacy and Food Sciences, Institute of Biomedicine (IBUB), Institute of Theoretical and Computational Chemistry (IQTCUB) and Campus Torribera, University of Barcelona, 08921, Santa Coloma De Gramenet, Spain
- Department of Pharmacy and Pharmaceutical Technology and Physical Chemistry and Institute of Theoretical and Computational Chemistry (IQTCUB), Faculty of Pharmacy and Food Sciences, University of Barcelona, Av. Joan XXIII 27-31, 08028, Barcelona, Spain
| | - Silvana Pinheiro
- Department of Pharmacy and Pharmaceutical Technology and Physical Chemistry and Institute of Theoretical and Computational Chemistry (IQTCUB), Faculty of Pharmacy and Food Sciences, University of Barcelona, Av. Joan XXIII 27-31, 08028, Barcelona, Spain
| | - Kilian German
- Department of Chemical Engineering and Analytical Chemistry and Institute of Biomedicine (IBUB), Universitat de Barcelona, Martí i Franquès 1-11, 08028, Barcelona, Spain
| | - Clara Ràfols
- Department of Chemical Engineering and Analytical Chemistry and Institute of Biomedicine (IBUB), Universitat de Barcelona, Martí i Franquès 1-11, 08028, Barcelona, Spain
| | - Carles Curutchet
- Department of Pharmacy and Pharmaceutical Technology and Physical Chemistry and Institute of Theoretical and Computational Chemistry (IQTCUB), Faculty of Pharmacy and Food Sciences, University of Barcelona, Av. Joan XXIII 27-31, 08028, Barcelona, Spain
| | - F Javier Luque
- Department of Nutrition, Food Science and Gastronomy, Faculty of Pharmacy and Food Sciences, Institute of Biomedicine (IBUB), Institute of Theoretical and Computational Chemistry (IQTCUB) and Campus Torribera, University of Barcelona, 08921, Santa Coloma De Gramenet, Spain.
| |
Collapse
|
20
|
Hölzl C, Kibies P, Imoto S, Noetzel J, Knierbein M, Salmen P, Paulus M, Nase J, Held C, Sadowski G, Marx D, Kast SM, Horinek D. Structure and thermodynamics of aqueous urea solutions from ambient to kilobar pressures: From thermodynamic modeling, experiments, and first principles simulations to an accurate force field description. Biophys Chem 2019; 254:106260. [PMID: 31522071 DOI: 10.1016/j.bpc.2019.106260] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2019] [Revised: 08/29/2019] [Accepted: 08/29/2019] [Indexed: 12/17/2022]
Abstract
Molecular simulations based on classical force fields are a powerful method for shedding light on the complex behavior of biomolecules in solution. When cosolutes are present in addition to water and biomolecules, subtle balances of weak intermolecular forces have to be accounted for. This imposes high demands on the quality of the underlying force fields, and therefore force field development for small cosolutes is still an active field. Here, we present the development of a new urea force field from studies of urea solutions at ambient and elevated hydrostatic pressures based on a combination of experimental and theoretical approaches. Experimental densities and solvation shell properties from ab initio molecular dynamics simulations at ambient conditions served as the target properties for the force field optimization. Since urea is present in many marine life forms, elevated hydrostatic pressure was rigorously addressed: densities at high pressure were measured by vibrating tube densitometry up to 500 bar and by X-ray absorption up to 5 kbar. Densities were determined by the perturbed-chain statistical associating fluid theory equation of state. Solvation properties were determined by embedded cluster integral equation theory and ab initio molecular dynamics. Our new force field is able to capture the properties of urea solutions at high pressures without further high-pressure adaption, unlike trimethylamine-N-oxide, for which a high-pressure adaption is necessary.
Collapse
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
| | - Jan Noetzel
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Michael Knierbein
- Laboratory of Thermodynamics, Technische Universität Dortmund, 44227 Dortmund, Germany
| | - Paul Salmen
- Fakultät Physik/DELTA, Technische Universität Dortmund, 44221 Dortmund, Germany
| | - Michael Paulus
- Fakultät Physik/DELTA, Technische Universität Dortmund, 44221 Dortmund, Germany
| | - Julia Nase
- Fakultät Physik/DELTA, Technische Universität Dortmund, 44221 Dortmund, Germany
| | - Christoph Held
- Laboratory of Thermodynamics, Technische Universität Dortmund, 44227 Dortmund, Germany
| | - Gabriele Sadowski
- Laboratory of Thermodynamics, Technische Universität Dortmund, 44227 Dortmund, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, 44227 Dortmund, Germany.
| | - Dominik Horinek
- Institut für Physikalische und Theoretische Chemie, Universität Regensburg, 93040 Regensburg, Germany.
| |
Collapse
|
21
|
Munte CE, Karl M, Kauter W, Eberlein L, Pham TV, Erlach MB, Kast SM, Kremer W, Kalbitzer HR. High pressure response of 1H NMR chemical shifts of purine nucleotides. Biophys Chem 2019; 254:106261. [PMID: 31522070 DOI: 10.1016/j.bpc.2019.106261] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2019] [Revised: 09/01/2019] [Accepted: 09/01/2019] [Indexed: 11/25/2022]
Abstract
The study of the pressure response by NMR spectroscopy provides information on the thermodynamics of conformational equilibria in proteins and nucleic acids. For obtaining a database for expected pressure effects on free nucleotides and nucleotides bound in macromolecular complexes, the pressure response of 1H chemical shifts and J-coupling constants of the purine 5'-ribonucleotides AMP, ADP, ATP, GMP, GDP, and GTP were studied in the absence and presence of Mg2+-ions. Experiments are supported by quantum-chemical calculations of populations and chemical shift differences in order to corroborate structural interpretations and to estimate missing data for AMP. The preference of the ribose S puckering obtained from the analysis of the experimental J-couplings is also confirmed by the calculations. In addition, the pressure response of the non-hydrolysable GTP analogues GppNHp, GppCH2p, and GTPγS was examined within a pressure range up to 200 MPa. As observed earlier for 31P NMR chemical shifts of these nucleotides the pressure dependence of chemical shifts is clearly non-linear in most cases. In di- and tri-phospho nucleosides, the resonances of the two protons bound to the ribose 5' carbon are non-equivalent and can be observed separately. The gg-rotamer at C4'- C5' bond is strongly preferred and the downfield shifted resonance can be assigned to the H5″ proton in the nucleotides. In contrast, in adenosine itself the frequencies of the two resonances are interchanged.
Collapse
Affiliation(s)
- Claudia E Munte
- University of Regensburg, Institute of Biophysics and Physical Biochemistry, Center of Magnetic Resonance in Chemistry and Biomedicine, Universitätsstraße 31, 93053 Regensburg, Germany
| | - Matthias Karl
- University of Regensburg, Institute of Biophysics and Physical Biochemistry, Center of Magnetic Resonance in Chemistry and Biomedicine, Universitätsstraße 31, 93053 Regensburg, Germany
| | - Waldemar Kauter
- University of Regensburg, Institute of Biophysics and Physical Biochemistry, Center of Magnetic Resonance in Chemistry and Biomedicine, Universitätsstraße 31, 93053 Regensburg, Germany
| | - Lukas Eberlein
- TU Dortmund University, Physical Chemistry III, Otto-Hahn-Straße 4a, 44227 Dortmund, Germany
| | - Thuy-Vy Pham
- University of Regensburg, Institute of Biophysics and Physical Biochemistry, Center of Magnetic Resonance in Chemistry and Biomedicine, Universitätsstraße 31, 93053 Regensburg, Germany
| | - Markus Beck Erlach
- University of Regensburg, Institute of Biophysics and Physical Biochemistry, Center of Magnetic Resonance in Chemistry and Biomedicine, Universitätsstraße 31, 93053 Regensburg, Germany
| | - Stefan M Kast
- TU Dortmund University, Physical Chemistry III, Otto-Hahn-Straße 4a, 44227 Dortmund, Germany
| | - Werner Kremer
- University of Regensburg, Institute of Biophysics and Physical Biochemistry, Center of Magnetic Resonance in Chemistry and Biomedicine, Universitätsstraße 31, 93053 Regensburg, Germany
| | - Hans Robert Kalbitzer
- University of Regensburg, Institute of Biophysics and Physical Biochemistry, Center of Magnetic Resonance in Chemistry and Biomedicine, Universitätsstraße 31, 93053 Regensburg, Germany.
| |
Collapse
|
22
|
Tielker N, Eberlein L, Chodun C, Güssregen S, Kast SM. pK a calculations for tautomerizable and conformationally flexible molecules: partition function vs. state transition approach. J Mol Model 2019; 25:139. [PMID: 31041535 DOI: 10.1007/s00894-019-4033-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2019] [Accepted: 04/07/2019] [Indexed: 11/26/2022]
Abstract
Calculations of acidities of molecules with multiple tautomeric and/or conformational states require adequate treatment of the relative energetics of accessible states accompanied by a statistical-mechanical formulation of their contribution to the macroscopic pKa value. Here, we demonstrate rigorously the formal equivalence of two such approaches: a partition function treatment and statistics over transitions between molecular tautomeric and conformational states in the limit of a theory that does not require adjustment by empirical parameters correcting energetic values. However, for a frequently employed correction scheme, linear scaling of (free) energies and regression with respect to reference data taking an additive constant into account, this equivalence breaks down if more than one acid or base state is involved. The consequences of the resulting inconsistency are discussed on our datasets developed for aqueous pKa predictions during the recent SAMPL6 challenge, where molecular state energetics were computed based on the "embedded cluster reference interaction site model" (EC-RISM). This method couples integral equation theory as a solvation model to quantum-chemical calculations and yielded a test set root mean square error of 1.1 pK units from a partition function ansatz. For all practical purposes, the present results indicate that a state transition approach yields comparable accuracy despite the formal theoretical inconsistency, and that an additive regression intercept, which is strictly constant in the limit of large compound mass only, is a valid approximation. Graphical abstract Embedded cluster reference interaction site model-derived vs. experimental pKa for the test set calculated with either the partition function (blue) or the state transition approach (red), using m as a free parameter.
Collapse
Affiliation(s)
- Nicolas Tielker
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Lukas Eberlein
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Christian Chodun
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Stefan Güssregen
- R&D Integrated Drug Discovery, Sanofi-Aventis Deutschland GmbH, 65926, Frankfurt am Main, Germany
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany.
| |
Collapse
|
23
|
Tielker N, Eberlein L, Güssregen S, Kast SM. The SAMPL6 challenge on predicting aqueous pKa values from EC-RISM theory. J Comput Aided Mol Des 2018; 32:1151-1163. [DOI: 10.1007/s10822-018-0140-z] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2018] [Accepted: 07/28/2018] [Indexed: 12/15/2022]
|
24
|
Matter H, Güssregen S. Characterizing hydration sites in protein-ligand complexes towards the design of novel ligands. Bioorg Med Chem Lett 2018; 28:2343-2352. [PMID: 29880400 DOI: 10.1016/j.bmcl.2018.05.061] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2018] [Revised: 05/28/2018] [Accepted: 05/30/2018] [Indexed: 11/18/2022]
Abstract
Water is an essential part of protein binding sites and mediates interactions to ligands. Its displacement by ligand parts affects the free binding energy of resulting protein-ligand complexes. Therefore the characterization of solvation properties is important for design. Of particular interest is the propensity of localized water to be favorably displaced by a ligand. This review discusses two popular computational approaches addressing these questions, namely WaterMap based on statistical mechanics analysis of MD simulations and 3D RISM based on integral equation theory of liquids. The theoretical background and recent applications in structure-based design will be presented.
Collapse
Affiliation(s)
- Hans Matter
- Sanofi-Aventis Deutschland GmbH, Integrated Drug Discovery (IDD), Synthetic Molecular Design, Building G838, Industriepark Höchst, 65926 Frankfurt am Main, Germany.
| | - Stefan Güssregen
- Sanofi-Aventis Deutschland GmbH, Integrated Drug Discovery (IDD), Synthetic Molecular Design, Building G838, Industriepark Höchst, 65926 Frankfurt am Main, Germany
| |
Collapse
|
25
|
Ding L, Levesque M, Borgis D, Belloni L. Efficient molecular density functional theory using generalized spherical harmonics expansions. J Chem Phys 2017; 147:094107. [DOI: 10.1063/1.4994281] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Lu Ding
- Maison de la Simulation, USR 3441 CNRS-CEA-Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Maximilien Levesque
- PASTEUR, Département de Chimie, École Normale Supérieure, UPMC Univ. Paris 06, CNRS, PSL Research University, 75005 Paris, France
- Sorbonne Universités, UPMC Univ. Paris 06, École Normale Supérieure, CNRS, Processus d’Activation Sélective par Transfert d’Énergie Uni-Électronique ou Radiatif (PASTEUR), 75005 Paris, France
| | - Daniel Borgis
- Maison de la Simulation, USR 3441 CNRS-CEA-Université Paris-Saclay, 91191 Gif-sur-Yvette, France
- PASTEUR, Département de Chimie, École Normale Supérieure, UPMC Univ. Paris 06, CNRS, PSL Research University, 75005 Paris, France
- Sorbonne Universités, UPMC Univ. Paris 06, École Normale Supérieure, CNRS, Processus d’Activation Sélective par Transfert d’Énergie Uni-Électronique ou Radiatif (PASTEUR), 75005 Paris, France
| | - Luc Belloni
- LIONS, NIMBE, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| |
Collapse
|
26
|
Thomas HB, Hennemann M, Kibies P, Hoffgaard F, Güssregen S, Hessler G, Kast SM, Clark T. The hpCADD NDDO Hamiltonian: Parametrization. J Chem Inf Model 2017; 57:1907-1922. [DOI: 10.1021/acs.jcim.7b00080] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Heike B. Thomas
- Computer-Chemie-Centrum,
Department Chemie und Pharmazie, Friedrich-Alexander-Universität Erlangen-Nürnberg, Nägelsbachstr. 25, 91052 Erlangen, Germany
| | - Matthias Hennemann
- Computer-Chemie-Centrum,
Department Chemie und Pharmazie, Friedrich-Alexander-Universität Erlangen-Nürnberg, Nägelsbachstr. 25, 91052 Erlangen, Germany
| | - Patrick Kibies
- Physikalische
Chemie III, Technische Universität Dortmund, Otto-Hahn-Str.
4a, 44227 Dortmund, Germany
| | - Franziska Hoffgaard
- Physikalische
Chemie III, Technische Universität Dortmund, Otto-Hahn-Str.
4a, 44227 Dortmund, Germany
| | - Stefan Güssregen
- R&D, IDD, Structural Design and Informatics, Sanofi-Aventis Deutschland GmbH, Industriepark Höchst, 65926 Frankfurt am Main, Germany
| | - Gerhard Hessler
- R&D, IDD, Structural Design and Informatics, Sanofi-Aventis Deutschland GmbH, Industriepark Höchst, 65926 Frankfurt am Main, Germany
| | - Stefan M. Kast
- Physikalische
Chemie III, Technische Universität Dortmund, Otto-Hahn-Str.
4a, 44227 Dortmund, Germany
| | - Timothy Clark
- Computer-Chemie-Centrum,
Department Chemie und Pharmazie, Friedrich-Alexander-Universität Erlangen-Nürnberg, Nägelsbachstr. 25, 91052 Erlangen, Germany
| |
Collapse
|
27
|
Sure R, Grimme S. Comprehensive Benchmark of Association (Free) Energies of Realistic Host-Guest Complexes. J Chem Theory Comput 2016; 11:3785-801. [PMID: 26574460 DOI: 10.1021/acs.jctc.5b00296] [Citation(s) in RCA: 153] [Impact Index Per Article: 19.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The S12L test set for supramolecular Gibbs free energies of association ΔGa (Grimme, S. Chem. Eur. J. 2012, 18, 9955-9964) is extended to 30 complexes (S30L), featuring more diverse interaction motifs, anions, and higher charges (-1 up to +4) as well as larger systems with up to 200 atoms. Various typical noncovalent interactions like hydrogen and halogen bonding, π-π stacking, nonpolar dispersion, and CH-π and cation-dipolar interactions are represented by "real" complexes. The experimental Gibbs free energies of association (ΔGa exp) cover a wide range from -0.7 to -24.7 kcal mol-1. In order to obtain a theoretical best estimate for ΔGa, we test various dispersion corrected density functionals in combination with quadruple-ζ basis sets for calculating the association energies in the gas phase. Further, modern semiempirical methods are employed to obtain the thermostatistical corrections from energy to Gibbs free energy, and the COSMO-RS model with several parametrizations as well as the SMD model are used to include solvation contributions. We investigate the effect of including counterions for the charged systems (S30L-CI), which is found to overall improve the results. Our best method combination consists of PW6B95-D3 (for neutral and charged systems) or ωB97X-D3 (for systems with counterions) energies, HF-3c thermostatistical corrections, and Gibbs free energies of solvation obtained with the COSMO-RS 2012 parameters for nonpolar solvents and 2013-fine for water. This combination gives a mean absolute deviation for ΔGa of only 2.4 kcal mol-1 (S30L) and 2.1 kcal mol-1 (S30L-CI), with a mean deviation of almost zero compared to experiment. Regarding the relative Gibbs free energies of association for the 13 pairs of complexes which share the same host, the correct trend in binding affinities could be reproduced except for two cases. The MAD compared to experiment amounts to 1.2 kcal mol-1, and the MD is almost zero. The best-estimate theoretical corrections are used to back-correct the experimental ΔGa values in order to get an empirical estimate for the "experimental", zero-point vibrational energy exclusive, gas phase binding energies. These are then utilized to benchmark the performance of various "lowcost" quantum chemical methods for noncovalent interactions in large systems. The performance of other common DFT methods as well as the use of semiempirical methods for structure optimizations is discussed.
Collapse
Affiliation(s)
- Rebecca Sure
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn Beringstr. 4, D-53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
28
|
Tielker N, Tomazic D, Heil J, Kloss T, Ehrhart S, Güssregen S, Schmidt KF, Kast SM. The SAMPL5 challenge for embedded-cluster integral equation theory: solvation free energies, aqueous pK a, and cyclohexane-water log D. J Comput Aided Mol Des 2016; 30:1035-1044. [PMID: 27554666 DOI: 10.1007/s10822-016-9939-7] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2016] [Accepted: 08/17/2016] [Indexed: 12/01/2022]
Abstract
We predict cyclohexane-water distribution coefficients (log D 7.4) for drug-like molecules taken from the SAMPL5 blind prediction challenge by the "embedded cluster reference interaction site model" (EC-RISM) integral equation theory. This task involves the coupled problem of predicting both partition coefficients (log P) of neutral species between the solvents and aqueous acidity constants (pK a) in order to account for a change of protonation states. The first issue is addressed by calibrating an EC-RISM-based model for solvation free energies derived from the "Minnesota Solvation Database" (MNSOL) for both water and cyclohexane utilizing a correction based on the partial molar volume, yielding a root mean square error (RMSE) of 2.4 kcal mol-1 for water and 0.8-0.9 kcal mol-1 for cyclohexane depending on the parametrization. The second one is treated by employing on one hand an empirical pK a model (MoKa) and, on the other hand, an EC-RISM-derived regression of published acidity constants (RMSE of 1.5 for a single model covering acids and bases). In total, at most 8 adjustable parameters are necessary (2-3 for each solvent and two for the pK a) for training solvation and acidity models. Applying the final models to the log D 7.4 dataset corresponds to evaluating an independent test set comprising other, composite observables, yielding, for different cyclohexane parametrizations, 2.0-2.1 for the RMSE with the first and 2.2-2.8 with the combined first and second SAMPL5 data set batches. Notably, a pure log P model (assuming neutral species only) performs statistically similarly for these particular compounds. The nature of the approximations and possible perspectives for future developments are discussed.
Collapse
Affiliation(s)
- Nicolas Tielker
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Daniel Tomazic
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Jochen Heil
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Thomas Kloss
- IPhT, L'Orme des Merisiers, CEA-Saclay, 91191, Gif-sur-Yvette, France
| | | | - Stefan Güssregen
- Sanofi-Aventis Deutschland GmbH, Industriepark Höchst, 65926, Frankfurt am Main, Germany
| | - K Friedemann Schmidt
- Sanofi-Aventis Deutschland GmbH, Industriepark Höchst, 65926, Frankfurt am Main, Germany
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany.
| |
Collapse
|
29
|
Imoto S, Kibies P, Rosin C, Winter R, Kast SM, Marx D. Toward Extreme Biophysics: Deciphering the Infrared Response of Biomolecular Solutions at High Pressures. Angew Chem Int Ed Engl 2016. [DOI: 10.1002/ange.201602757] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Sho Imoto
- Lehrstuhl für Theoretische Chemie; Ruhr-Universität Bochum; 44780 Bochum Germany
| | - Patrick Kibies
- Physikalische Chemie III; Technische Universität Dortmund; 44227 Dortmund Germany
| | - Christopher Rosin
- Physikalische Chemie I; Technische Universität Dortmund; 44227 Dortmund Germany
| | - Roland Winter
- Physikalische Chemie I; Technische Universität Dortmund; 44227 Dortmund Germany
| | - Stefan M. Kast
- Physikalische Chemie III; Technische Universität Dortmund; 44227 Dortmund Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie; Ruhr-Universität Bochum; 44780 Bochum Germany
| |
Collapse
|
30
|
Imoto S, Kibies P, Rosin C, Winter R, Kast SM, Marx D. Toward Extreme Biophysics: Deciphering the Infrared Response of Biomolecular Solutions at High Pressures. Angew Chem Int Ed Engl 2016; 55:9534-8. [DOI: 10.1002/anie.201602757] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2016] [Revised: 05/18/2016] [Indexed: 11/10/2022]
Affiliation(s)
- Sho Imoto
- Lehrstuhl für Theoretische Chemie; Ruhr-Universität Bochum; 44780 Bochum Germany
| | - Patrick Kibies
- Physikalische Chemie III; Technische Universität Dortmund; 44227 Dortmund Germany
| | - Christopher Rosin
- Physikalische Chemie I; Technische Universität Dortmund; 44227 Dortmund Germany
| | - Roland Winter
- Physikalische Chemie I; Technische Universität Dortmund; 44227 Dortmund Germany
| | - Stefan M. Kast
- Physikalische Chemie III; Technische Universität Dortmund; 44227 Dortmund Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie; Ruhr-Universität Bochum; 44780 Bochum Germany
| |
Collapse
|
31
|
Jeanmairet G, Levy N, Levesque M, Borgis D. Molecular density functional theory of water including density-polarization coupling. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:244005. [PMID: 27116250 DOI: 10.1088/0953-8984/28/24/244005] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
We present a three-dimensional molecular density functional theory of water derived from first-principles that relies on the particle's density and multipolar polarization density and includes the density-polarization coupling. This brings two main benefits: (i) scalar density and vectorial multipolar polarization density fields are much more tractable and give more physical insight than the full position and orientation densities, and (ii) it includes the full density-polarization coupling of water, that is known to be non-vanishing but has never been taken into account. Furthermore, the theory requires only the partial charge distribution of a water molecule and three measurable bulk properties, namely the structure factor and the Fourier components of the longitudinal and transverse dielectric susceptibilities.
Collapse
Affiliation(s)
- Guillaume Jeanmairet
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, Stuttgart 70569, Germany
| | | | | | | |
Collapse
|
32
|
Frach R, Kibies P, Böttcher S, Pongratz T, Strohfeldt S, Kurrmann S, Koehler J, Hofmann M, Kremer W, Kalbitzer HR, Reiser O, Horinek D, Kast SM. The Chemical Shift Baseline for High-Pressure NMR Spectra of Proteins. Angew Chem Int Ed Engl 2016; 55:8757-60. [DOI: 10.1002/anie.201602054] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2016] [Indexed: 11/09/2022]
Affiliation(s)
- Roland Frach
- Physikalische Chemie III, Fakultät für Chemie und Chemische Biologie; TU Dortmund; Otto-Hahn-Strasse 4a 44227 Dortmund Germany
| | - Patrick Kibies
- Physikalische Chemie III, Fakultät für Chemie und Chemische Biologie; TU Dortmund; Otto-Hahn-Strasse 4a 44227 Dortmund Germany
| | - Saraphina Böttcher
- Physikalische Chemie III, Fakultät für Chemie und Chemische Biologie; TU Dortmund; Otto-Hahn-Strasse 4a 44227 Dortmund Germany
| | - Tim Pongratz
- Physikalische Chemie III, Fakultät für Chemie und Chemische Biologie; TU Dortmund; Otto-Hahn-Strasse 4a 44227 Dortmund Germany
| | - Steven Strohfeldt
- Physikalische Chemie III, Fakultät für Chemie und Chemische Biologie; TU Dortmund; Otto-Hahn-Strasse 4a 44227 Dortmund Germany
| | - Simon Kurrmann
- Institut für Biophysik und physikalische Biochemie, Fakultät für Biologie und Vorklin. Medizin; Universität Regensburg; Germany
| | - Joerg Koehler
- Institut für Biophysik und physikalische Biochemie, Fakultät für Biologie und Vorklin. Medizin; Universität Regensburg; Germany
| | - Martin Hofmann
- Institut für Organische Chemie, Fakultät für Chemie und Pharmazie; Universität Regensburg; Germany
| | - Werner Kremer
- Institut für Biophysik und physikalische Biochemie, Fakultät für Biologie und Vorklin. Medizin; Universität Regensburg; Germany
| | - Hans Robert Kalbitzer
- Institut für Biophysik und physikalische Biochemie, Fakultät für Biologie und Vorklin. Medizin; Universität Regensburg; Germany
| | - Oliver Reiser
- Institut für Organische Chemie, Fakultät für Chemie und Pharmazie; Universität Regensburg; Germany
| | - Dominik Horinek
- Institut für Physikalische und Theoretische Chemie, Fakultät für Chemie und Pharmazie; Universität Regensburg; Germany
| | - Stefan M. Kast
- Physikalische Chemie III, Fakultät für Chemie und Chemische Biologie; TU Dortmund; Otto-Hahn-Strasse 4a 44227 Dortmund Germany
| |
Collapse
|
33
|
Frach R, Kibies P, Böttcher S, Pongratz T, Strohfeldt S, Kurrmann S, Koehler J, Hofmann M, Kremer W, Kalbitzer HR, Reiser O, Horinek D, Kast SM. Die Basislinie der chemischen Verschiebung in Hochdruck-NMR-Spektren von Proteinen. Angew Chem Int Ed Engl 2016. [DOI: 10.1002/ange.201602054] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Roland Frach
- Physikalische Chemie III, Fakultät für Chemie und Chemische Biologie; TU Dortmund; Otto-Hahn-Straße 4a 44227 Dortmund Deutschland
| | - Patrick Kibies
- Physikalische Chemie III, Fakultät für Chemie und Chemische Biologie; TU Dortmund; Otto-Hahn-Straße 4a 44227 Dortmund Deutschland
| | - Saraphina Böttcher
- Physikalische Chemie III, Fakultät für Chemie und Chemische Biologie; TU Dortmund; Otto-Hahn-Straße 4a 44227 Dortmund Deutschland
| | - Tim Pongratz
- Physikalische Chemie III, Fakultät für Chemie und Chemische Biologie; TU Dortmund; Otto-Hahn-Straße 4a 44227 Dortmund Deutschland
| | - Steven Strohfeldt
- Physikalische Chemie III, Fakultät für Chemie und Chemische Biologie; TU Dortmund; Otto-Hahn-Straße 4a 44227 Dortmund Deutschland
| | - Simon Kurrmann
- Institut für Biophysik und physikalische Biochemie, Fakultät für Biologie und Vorklin. Medizin; Universität Regensburg; Deutschland
| | - Joerg Koehler
- Institut für Biophysik und physikalische Biochemie, Fakultät für Biologie und Vorklin. Medizin; Universität Regensburg; Deutschland
| | - Martin Hofmann
- Institut für Organische Chemie, Fakultät für Chemie und Pharmazie; Universität Regensburg; Deutschland
| | - Werner Kremer
- Institut für Biophysik und physikalische Biochemie, Fakultät für Biologie und Vorklin. Medizin; Universität Regensburg; Deutschland
| | - Hans Robert Kalbitzer
- Institut für Biophysik und physikalische Biochemie, Fakultät für Biologie und Vorklin. Medizin; Universität Regensburg; Deutschland
| | - Oliver Reiser
- Institut für Organische Chemie, Fakultät für Chemie und Pharmazie; Universität Regensburg; Deutschland
| | - Dominik Horinek
- Institut für Physikalische und Theoretische Chemie, Fakultät für Chemie und Pharmazie; Universität Regensburg; Deutschland
| | - Stefan M. Kast
- Physikalische Chemie III, Fakultät für Chemie und Chemische Biologie; TU Dortmund; Otto-Hahn-Straße 4a 44227 Dortmund Deutschland
| |
Collapse
|
34
|
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
|
35
|
Frach R, Heil J, Kast SM. Structure and thermodynamics of nondipolar molecular liquids and solutions from integral equation theory. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1167266] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- Roland Frach
- Physikalische Chemie III, Technische Universität Dortmund, Dortmund, Germany
| | - Jochen Heil
- Physikalische Chemie III, Technische Universität Dortmund, Dortmund, Germany
| | - Stefan M. Kast
- Physikalische Chemie III, Technische Universität Dortmund, Dortmund, Germany
| |
Collapse
|
36
|
Misin M, Fedorov MV, Palmer DS. Hydration Free Energies of Molecular Ions from Theory and Simulation. J Phys Chem B 2016; 120:975-83. [PMID: 26756333 DOI: 10.1021/acs.jpcb.5b10809] [Citation(s) in RCA: 57] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a theoretical/computational framework for accurate calculation of hydration free energies of ionized molecular species. The method is based on a molecular theory, 3D-RISM, combined with a recently developed pressure correction (PC+). The 3D-RISM/PC+ model can provide ∼3 kcal/mol hydration free energy accuracy for a large variety of ionic compounds, provided that the Galvani potential of water is taken into account. The results are compared with direct atomistic simulations. Several methodological aspects of hydration free energy calculations for charged species are discussed.
Collapse
Affiliation(s)
| | | | - David S Palmer
- Department of Pure and Applied Chemistry, University of Strathclyde , 295 Cathedral Street, Glasgow, G1 1XL, United Kingdom
| |
Collapse
|
37
|
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]
|
38
|
Engel J, Richters A, Getlik M, Tomassi S, Keul M, Termathe M, Lategahn J, Becker C, Mayer-Wrangowski S, Grütter C, Uhlenbrock N, Krüll J, Schaumann N, Eppmann S, Kibies P, Hoffgaard F, Heil J, Menninger S, Ortiz-Cuaran S, Heuckmann JM, Tinnefeld V, Zahedi RP, Sos ML, Schultz-Fademrecht C, Thomas RK, Kast SM, Rauh D. Targeting Drug Resistance in EGFR with Covalent Inhibitors: A Structure-Based Design Approach. J Med Chem 2015; 58:6844-63. [PMID: 26275028 DOI: 10.1021/acs.jmedchem.5b01082] [Citation(s) in RCA: 76] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Receptor tyrosine kinases represent one of the prime targets in cancer therapy, as the dysregulation of these elementary transducers of extracellular signals, like the epidermal growth factor receptor (EGFR), contributes to the onset of cancer, such as non-small cell lung cancer (NSCLC). Strong efforts were directed to the development of irreversible inhibitors and led to compound CO-1686, which takes advantage of increased residence time at EGFR by alkylating Cys797 and thereby preventing toxic effects. Here, we present a structure-based approach, rationalized by subsequent computational analysis of conformational ligand ensembles in solution, to design novel and irreversible EGFR inhibitors based on a screening hit that was identified in a phenotype screen of 80 NSCLC cell lines against approximately 1500 compounds. Using protein X-ray crystallography, we deciphered the binding mode in engineered cSrc (T338M/S345C), a validated model system for EGFR-T790M, which constituted the basis for further rational design approaches. Chemical synthesis led to further compound collections that revealed increased biochemical potency and, in part, selectivity toward mutated (L858R and L858R/T790M) vs nonmutated EGFR. Further cell-based and kinetic studies were performed to substantiate our initial findings. Utilizing proteolytic digestion and nano-LC-MS/MS analysis, we confirmed the alkylation of Cys797.
Collapse
Affiliation(s)
- Julian Engel
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - André Richters
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Matthäus Getlik
- Chemical Genomics Centre of the Max-Planck Society , Otto-Hahn-Straße 15, D-44227 Dortmund, Germany
| | - Stefano Tomassi
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Marina Keul
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Martin Termathe
- Chemical Genomics Centre of the Max-Planck Society , Otto-Hahn-Straße 15, D-44227 Dortmund, Germany
| | - Jonas Lategahn
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Christian Becker
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Svenja Mayer-Wrangowski
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Christian Grütter
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Niklas Uhlenbrock
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Jasmin Krüll
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Niklas Schaumann
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Simone Eppmann
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Patrick Kibies
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Franziska Hoffgaard
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Jochen Heil
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Sascha Menninger
- Lead Discovery Center GmbH , Otto-Hahn-Straße 15, D-44227 Dortmund, Germany
| | - Sandra Ortiz-Cuaran
- Department of Translational Genomics, Medical Faculty, University of Cologne , Weyertal 115b, D-50931 Cologne, Germany
| | - Johannes M Heuckmann
- Department of Translational Genomics, Medical Faculty, University of Cologne , Weyertal 115b, D-50931 Cologne, Germany
| | - Verena Tinnefeld
- Leibniz-Institut für Analytische Wissenschaften - ISAS - e.V. , D-44139 Dortmund, Germany
| | - René P Zahedi
- Leibniz-Institut für Analytische Wissenschaften - ISAS - e.V. , D-44139 Dortmund, Germany
| | - Martin L Sos
- Department of Translational Genomics, Medical Faculty, University of Cologne , Weyertal 115b, D-50931 Cologne, Germany.,Molecular Pathology, University Hospital of Cologne , Kerpenerstraße 62, D-50937 Cologne, Germany
| | | | - Roman K Thomas
- Department of Translational Genomics, Medical Faculty, University of Cologne , Weyertal 115b, D-50931 Cologne, Germany.,Department of Pathology, University of Cologne , Joseph-Stelzmann Straße 9, D-50931 Cologne, Germany
| | - Stefan M Kast
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany
| | - Daniel Rauh
- Department of Chemistry and Chemical Biology, TU Dortmund University , Otto-Hahn-Straße 6, D-44227 Dortmund, Germany.,Chemical Genomics Centre of the Max-Planck Society , Otto-Hahn-Straße 15, D-44227 Dortmund, Germany
| |
Collapse
|
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
|
Abstract
The calculation of electrostatic solute-solvent interactions in 3D RISM ("three-dimensional reference interaction site model") integral equation theory is recast in a form that allows for a computational treatment analogous to the "particle-mesh Ewald" formalism as used for molecular simulations. In addition, relations that connect 3D RISM correlation functions and interaction potentials with thermodynamic quantities such as the chemical potential and average solute-solvent interaction energy are reformulated in a way that calculations of expensive real-space electrostatic terms on the 3D grid are completely avoided. These methodical enhancements allow for both, a significant speedup particularly for large solute systems and a smoother convergence of predicted thermodynamic quantities with respect to box size, as illustrated for several benchmark systems.
Collapse
Affiliation(s)
- Jochen Heil
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 6, 44227 Dortmund, Germany
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 6, 44227 Dortmund, Germany
| |
Collapse
|
42
|
Misin M, Fedorov MV, Palmer DS. Communication: Accurate hydration free energies at a wide range of temperatures from 3D-RISM. J Chem Phys 2015; 142:091105. [DOI: 10.1063/1.4914315] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Affiliation(s)
- Maksim Misin
- Department of Physics, SUPA, University of Strathclyde, 107 Rottenrow, Glasgow G4 0NG, United Kingdom
| | - Maxim V. Fedorov
- Department of Physics, SUPA, University of Strathclyde, 107 Rottenrow, Glasgow G4 0NG, United Kingdom
| | - David S. Palmer
- Department of Pure and Applied Chemistry, University of Strathclyde, Thomas Graham Building, 295 Cathedral Street, Glasgow G1 1XL, United Kingdom
| |
Collapse
|
43
|
Li B, Matveev AV, Krüger S, Rösch N. Uranyl solvation by a reference interaction site model. COMPUT THEOR CHEM 2015. [DOI: 10.1016/j.comptc.2014.10.032] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
44
|
Antony J, Sure R, Grimme S. Using dispersion-corrected density functional theory to understand supramolecular binding thermodynamics. Chem Commun (Camb) 2015; 51:1764-74. [DOI: 10.1039/c4cc06722c] [Citation(s) in RCA: 116] [Impact Index Per Article: 12.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
A recently published theoretical approach employing a nondynamic structure model using dispersion-corrected density functional theory (DFT-D3) to calculate equilibrium free energies of association (Chem. – Eur. J., 2012, 18, 9955–9964) is illustrated by its application to eight supramolecular complexes.
Collapse
Affiliation(s)
- Jens Antony
- Mulliken Center for Theoretical Chemistry
- Institut für Physikalische und Theoretische Chemie der Universität Bonn
- D-53115 Bonn
- Germany
| | - Rebecca Sure
- Mulliken Center for Theoretical Chemistry
- Institut für Physikalische und Theoretische Chemie der Universität Bonn
- D-53115 Bonn
- Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry
- Institut für Physikalische und Theoretische Chemie der Universität Bonn
- D-53115 Bonn
- Germany
| |
Collapse
|
45
|
Frach R, Kast SM. Solvation Effects on Chemical Shifts by Embedded Cluster Integral Equation Theory. J Phys Chem A 2014; 118:11620-8. [DOI: 10.1021/jp5084407] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Roland Frach
- Physikalische Chemie III, TU Dortmund, Otto-Hahn-Straße 6, 44227 Dortmund, Germany
| | - Stefan M. Kast
- Physikalische Chemie III, TU Dortmund, Otto-Hahn-Straße 6, 44227 Dortmund, Germany
| |
Collapse
|
46
|
Sergiievskyi VP, Jeanmairet G, Levesque M, Borgis D. Fast Computation of Solvation Free Energies with Molecular Density Functional Theory: Thermodynamic-Ensemble Partial Molar Volume Corrections. J Phys Chem Lett 2014; 5:1935-1942. [PMID: 26273876 DOI: 10.1021/jz500428s] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Molecular density functional theory (MDFT) offers an efficient implicit-solvent method to estimate molecule solvation free-energies, whereas conserving a fully molecular representation of the solvent. Even within a second-order approximation for the free-energy functional, the so-called homogeneous reference fluid approximation, we show that the hydration free-energies computed for a data set of 500 organic compounds are of similar quality as those obtained from molecular dynamics free-energy perturbation simulations, with a computer cost reduced by 2-3 orders of magnitude. This requires to introduce the proper partial volume correction to transform the results from the grand canonical to the isobaric-isotherm ensemble that is pertinent to experiments. We show that this correction can be extended to 3D-RISM calculations, giving a sound theoretical justification to empirical partial molar volume corrections that have been proposed recently.
Collapse
Affiliation(s)
- Volodymyr P Sergiievskyi
- †École Normale Supérieure - PSL Research University, Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Guillaume Jeanmairet
- †École Normale Supérieure - PSL Research University, Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Maximilien Levesque
- †École Normale Supérieure - PSL Research University, Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Daniel Borgis
- †École Normale Supérieure - PSL Research University, Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
- ‡Maison de la Simulation, CNRS USR 3441, CEA Saclay, 91191 Gif-sur-Yvette, France
| |
Collapse
|
47
|
Jeanmairet G, Marry V, Levesque M, Rotenberg B, Borgis D. Hydration of clays at the molecular scale: the promising perspective of classical density functional theory. Mol Phys 2014. [DOI: 10.1080/00268976.2014.899647] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
48
|
Heil J, Tomazic D, Egbers S, Kast SM. Acidity in DMSO from the embedded cluster integral equation quantum solvation model. J Mol Model 2014; 20:2161. [PMID: 24664119 DOI: 10.1007/s00894-014-2161-4] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2013] [Accepted: 01/27/2014] [Indexed: 10/25/2022]
Abstract
The embedded cluster reference interaction site model (EC-RISM) is applied to the prediction of acidity constants of organic molecules in dimethyl sulfoxide (DMSO) solution. EC-RISM is based on a self-consistent treatment of the solute's electronic structure and the solvent's structure by coupling quantum-chemical calculations with three-dimensional (3D) RISM integral equation theory. We compare available DMSO force fields with reference calculations obtained using the polarizable continuum model (PCM). The results are evaluated statistically using two different approaches to eliminating the proton contribution: a linear regression model and an analysis of pK(a) shifts for compound pairs. Suitable levels of theory for the integral equation methodology are benchmarked. The results are further analyzed and illustrated by visualizing solvent site distribution functions and comparing them with an aqueous environment.
Collapse
Affiliation(s)
- Jochen Heil
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Straße 6, 44227, Dortmund, Germany
| | | | | | | |
Collapse
|
49
|
Jeanmairet G, Levesque M, Borgis D. Molecular density functional theory of water describing hydrophobicity at short and long length scales. J Chem Phys 2013; 139:154101. [DOI: 10.1063/1.4824737] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
|
50
|
Hoffgaard F, Heil J, Kast SM. Three-Dimensional RISM Integral Equation Theory for Polarizable Solute Models. J Chem Theory Comput 2013; 9:4718-26. [PMID: 26583390 DOI: 10.1021/ct400699q] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Modeling solute polarizability is a key ingredient for improving the description of solvation phenomena. In recent years, polarizable molecular mechanics force fields have emerged that circumvent the limitations of classical fixed charge force fields by the ability to adapt their electrostatic potential distribution to a polarizing environment. Solvation phenomena are characterized by the solute's excess chemical potential, which can be computed by expensive fully atomistic free energy simulations. The alternative is to employ an implicit solvent model, which poses a challenge to the formulation of the solute-solvent interaction term within a polarizable framework. Here, we adapt the three-dimensional reference interaction site model (3D RISM) integral equation theory as a solvent model, which analytically yields the chemical potential, to the polarizable AMOEBA force field using an embedding cluster (EC-RISM) strategy. The methodology is analogous to our earlier approach to the coupling of a quantum-chemical solute description with a classical 3D RISM solvent. We describe the conceptual physical and algorithmic basis as well as the performance for several benchmark cases as a proof of principle. The results consistently show reasonable agreement between AMOEBA and quantum-chemical free energies in solution in general and allow for separate assessment of energetic and solvation-related contributions. We find that, depending on the parametrization, AMOEBA reproduces the chemical potential in better agreement with reference quantum-chemical calculations than the intramolecular energies, which suggests possible routes toward systematic improvement of polarizable force fields.
Collapse
Affiliation(s)
- Franziska Hoffgaard
- Physikalische Chemie III, TU Dortmund , Otto-Hahn-Str. 6, 44227 Dortmund, Germany
| | - Jochen Heil
- Physikalische Chemie III, TU Dortmund , Otto-Hahn-Str. 6, 44227 Dortmund, Germany
| | - Stefan M Kast
- Physikalische Chemie III, TU Dortmund , Otto-Hahn-Str. 6, 44227 Dortmund, Germany
| |
Collapse
|