1
|
da Silva Filho ASN, de Souza JIR, Carolina Roma A, da Silva JAB, Longo RL, Alves S, Belian MF. A combined experimental-molecular modeling study of crown ether europium complexes: Effects of the coordinated anion on structural and luminescence properties. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 308:123704. [PMID: 38070311 DOI: 10.1016/j.saa.2023.123704] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/03/2023] [Revised: 11/20/2023] [Accepted: 11/26/2023] [Indexed: 01/13/2024]
Abstract
It is reported the synthesis, characterization by elemental analysis, thermogravimetry; electronic absorption, infrared, excitation, and emission spectroscopies of the [Eu(12C4)(phen)2(X)n]X2 complexes, where 12C4 = 12-crown-4, phen = 1,10-phenanthroline, and X = F-, Cl-, Br-, SCN-, ClO4-, and NO3-. It is verified that the polarizability of the anion X- exerts remarkable effects on the emission process. As a general trend, lower wavenumbers for the 7F0→5L6, 7F0→5D2 and 7F0→5D1 transitions are associated with the anions with higher volumes and, consequently, higher polarizability. The molecular modeling results performed with quantum methods (RHF and DFT) suggest some relationships between the calculated structures, electronic, and luminescence properties with the presence of the LMCT (ligand-to-metal charge transfer) states, which explains the differences in the emission spectra of these complexes due to the coordinated anion.
Collapse
Affiliation(s)
| | - Jéssica I R de Souza
- Programa de Pós-Graduação em Ciência de Materiais, Universidade Federal de Pernambuco, 50670-901 Recife, PE, Brazil
| | - Ana Carolina Roma
- Departamento de Química Fundamental, Universidade Federal de Pernambuco, 50670-901, Recife, Pernambuco, Brazil
| | - Juliana A B da Silva
- Programa de Pós-Graduação em Química, Universidade Federal Rural de Pernambuco, Recife, PE 52171-900, Brazil; Núcleo Interdisciplinar de Ciências Exatas e da Natureza, CAA, Universidade Federal de Pernambuco, Nova Caruaru, Caruaru, PE 55.014-900, Brazil
| | - Ricardo L Longo
- Programa de Pós-Graduação em Ciência de Materiais, Universidade Federal de Pernambuco, 50670-901 Recife, PE, Brazil; Departamento de Química Fundamental, Universidade Federal de Pernambuco, 50670-901, Recife, Pernambuco, Brazil
| | - S Alves
- Programa de Pós-Graduação em Ciência de Materiais, Universidade Federal de Pernambuco, 50670-901 Recife, PE, Brazil; Departamento de Química Fundamental, Universidade Federal de Pernambuco, 50670-901, Recife, Pernambuco, Brazil
| | - Mônica F Belian
- Programa de Pós-Graduação em Química, Universidade Federal Rural de Pernambuco, Recife, PE 52171-900, Brazil; Departamento de Química, Universidade Federal Rural de Pernambuco, 52171-900 Recife, Pernambuco, Brazil.
| |
Collapse
|
2
|
Fan K, Zhang Y, Qiu Y, Zhang H. Impacts of targeting different hydration free energy references on the development of ion potentials. Phys Chem Chem Phys 2022; 24:16244-16262. [PMID: 35758314 DOI: 10.1039/d2cp01237e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Hydration free energy (HFE) as the most important solvation parameter is often targeted in ion model development, even though the reported values differ by dozens of kcal mol-1 mainly due to the experimentally undetermined HFE of the proton ΔG°(H+). The choice of ΔG°(H+) obviously affects the hydration of single ions and the relative HFE between the ions with different (magnitude or sign) charges, and the impacts of targeted HFEs on the ion solvation and ion-ion interactions are largely unrevealed. Here we designed point charge models of K+, Mg2+, Al3+, and Cl- ions targeting a variety of HFE references and then investigated the HFE influences on the simulations of dilute and concentrated ion solutions and of the salt ion pairs in gas, liquid, and solid phases. Targeting one more property of ion-water oxygen distances (IOD) leaves the ion-water binding distance invariant, while the binding strength increases with the decreasing (more negative) HFE of ions as a result of a decrease in ΔG°(H+) for the cation and an increase in ΔG°(H+) for the anion. The increase in ΔG°(H+) leads to strengthened cation-anion interactions and thus to close ion-ion contacts, low osmotic pressures, and small activity derivatives in concentrated ion solutions as well as too stable ion pairs of the salts in different phases. The ion diffusivity and water exchange rates around the ions are simply not HFE dependent but rather more complex. Targeting both the aqueous IOD and salt crystal properties of KCl was also attempted and the comparison between different models indicates the complexity and challenge in obtaining a balanced performance between different phases using classical force fields. Our results also support that a real ΔG°(H+) value of -259.8 kcal mol-1 recommended by Hünenberger and Reif guides ion models to reproduce ion-water and ion-ion interactions reasonably at relatively low salt concentrations. Simulations of a metalloprotein show that a relatively more positive ΔG°(H+) for Mg2+ model is better for a reasonable description of the metal binding network.
Collapse
Affiliation(s)
- Kun Fan
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Yejie Qiu
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| |
Collapse
|
3
|
Qiu Y, Jiang Y, Zhang Y, Zhang H. Rational Design of Nonbonded Point Charge Models for Monovalent Ions with Lennard-Jones 12-6 Potential. J Phys Chem B 2021; 125:13502-13518. [PMID: 34860517 DOI: 10.1021/acs.jpcb.1c09103] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Ions are of central importance in nature, and a variety of potential models was proposed to model ions in different phases for an in-depth exploration of ion-related systems. Here, we developed point charge models of 14 monovalent ions with the traditional 12-6 Lennard-Jones (LJ) potential for use in conjunction with 11 water models of TIP3P, OPC3, SPC/E, SPC/Eb, TIP3P-FB, a99SB-disp, TIP4P-Ew, OPC, TIP4P/2005, TIP4P-D, and TIP4P-FB. The designed models reproduced the real hydration free energy (HFE) of ions and the ion-oxygen distance (IOD) in the first hydration shell accurately and simultaneously, a performance similar to the previously reported 12-6-4 LJ-type ion models (12-6 LJ plus an attractive C4 term for cations or a repulsive one for anions). This work, along with our previous work on di-, tri-, and tetravalent metal cations (J. Chem. Inf. Model. 2021, 61, 4031-4044; J. Chem. Inf. Model. 2021, 61, 4613-4629), demonstrates the feasibility of the simple 12-6 LJ potential in ion modeling. In order for the 12-6 LJ potential to reproduce both the HFE and IOD, the LJ R parameters need to be close to Shannon's ionic radii for the highly charged cations and to the Stokes's van der Waals (vdW) radii for the monovalent ions. With an additional C4 term, the R parameters of 12-6-4 LJ ion models agree well with the Stokes's vdW radii and have a more physical meaning. It appears that the C4 term can be merged into the 12-6 LJ potential by a rational tuning of R and the LJ well depth. Simulations of the osmotic coefficients of alkali chloride solutions and the properties of gaseous and solid alkali halides indicate the necessity of further optimizing ion-ion interactions via, for instance, targeting more properties or using a more physical (polarizable) model.
Collapse
Affiliation(s)
- Yejie Qiu
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
4
|
Zhang Y, Jiang Y, Qiu Y, Zhang H. Rational Design of Nonbonded Point Charge Models for Highly Charged Metal Cations with Lennard-Jones 12-6 Potential. J Chem Inf Model 2021; 61:4613-4629. [PMID: 34467756 DOI: 10.1021/acs.jcim.1c00723] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Here, we developed nonbonded point charge models using a simple Lennard-Jones (LJ) 12-6 potential for highly charged metal cations (18 trivalent and 6 tetravalent ions) for use with 11 water models of TIP3P, OPC3, SPC/E, SPC/Eb, TIP3P-FB, a99SB-disp, TIP4P-Ew, OPC, TIP4P/2005, TIP4P-D, and TIP4P-FB. The designed models simultaneously reproduce the hydration free energy (HFE) and ion-oxygen distance (IOD) in the first hydration shell with an error within 1 kcal/mol and 0.01 Å on average, respectively, and yield reasonable coordination numbers for most cations. Such performance is equivalent to the previously reported point charge models using a more complex 12-6-4 LJ-type potential, while the LJ R parameters of our models are much close to Shannon's revised effective ion radii than that of the 12-6-4 models. Our designed models overestimate the diffusion constants of several trivalent ions by 5-68%. The performance in predicting osmotic coefficients of trivalent chlorides in aqueous solution depends on the salt type. A calibration of cation-anion interacting LJ parameters reproduces the experimental osmotic coefficients of an AlCl3 solution at 0.2-3.0 mol/L. The effectiveness of our new models is further demonstrated by simulating a metalloprotein system with four force field/water combinations. This work facilitates accurate modeling of metal-containing systems by a variety of force fields and water models in aqueous solution.
Collapse
Affiliation(s)
- Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Yejie Qiu
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
5
|
|
6
|
Zhang Y, Jiang Y, Peng J, Zhang H. Rational Design of Nonbonded Point Charge Models for Divalent Metal Cations with Lennard-Jones 12-6 Potential. J Chem Inf Model 2021; 61:4031-4044. [PMID: 34313132 DOI: 10.1021/acs.jcim.1c00580] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Exploring a metal-involved biochemical process at a molecular level often requires a reliable description of metal properties in aqueous solution by classical nonbonded models. An additional C4 term for considering ion-induced dipole interactions was previously proposed to supplement the widely used Lennard-Jones 12-6 potential (known as the 12-6-4 LJ-type model) with good accuracy. Here, we demonstrate an alternative to modeling divalent metal cations (M2+) with the traditional 12-6 LJ potential by developing nonbonded point charge models for use with 11 water models: TIP3P, SPC/E, SPC/Eb, TIP4P-Ew, TIP4P-D, and TIP4P/2005 and the more recent OPC3, TIP3P-FB, OPC, TIP4P-FB, and a99SB-disp. Our designed models simultaneously reproduce the experimental hydration free energy, ion-oxygen distance, and coordination number in the first hydration shell accurately for most of the metal cations, an accuracy equivalent to that of the complex 12-6-4 LJ-type and double exponential potential models. A systematic comparison with the existing M2+ models is presented as well in terms of effective ion radii, diffusion constants, water exchange rates, and ion-water interactions. Molecular dynamics simulations of metal substitution in Escherichia coli glyoxalase I variants show the great potential of our new models for metalloproteins.
Collapse
Affiliation(s)
- Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park 16802, Pennsylvania, United States
| | - Jiarong Peng
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| |
Collapse
|
7
|
Peng J, Zhang Y, Jiang Y, Zhang H. Developing and Assessing Nonbonded Dummy Models of Magnesium Ion with Different Hydration Free Energy References. J Chem Inf Model 2021; 61:2981-2997. [PMID: 34080414 DOI: 10.1021/acs.jcim.1c00281] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
A large diversity in the targeted hydration free energies (HFEs) during model parameterization of metal ions was reported in the literature with a difference by dozens of kcal/mol. Here, we developed a series of nonbonded dummy models of the Mg2+ ion targeting different HFE references in TIP3P water, followed by assessments of the designed models in the simulations of MgCl2 solution and biological systems. Together with the comparison of existing models, we conclude that the difference in the targeted HFEs has a limited influence on the model performance, while the usability of these models differs from case to case. The feasibility of reproducing more properties of Mg2+ such as diffusion constants and water exchange rates using a nonbonded dummy model is demonstrated. Underestimated activity derivative and osmotic coefficient of MgCl2 solutions in high concentration reveal a necessity for further optimization of ion-pair interactions. The developed dummy models are applicable to metal coordination with Asp, Glu, and His residues in metalloenzymes, and the performance in predicting monodentate or bidentate binding modes of Asp/Glu residues depends on the complexity of metal centers and the choice of protein force fields. When both the binding modes coexist, the nonbonded dummy models outperform point charge models, probably in need of considering polarization of metal-binding residues by, for instance, charge calibration in classical force fields. This work is valuable for the use and further development of magnesium ion models for simulations of metal-containing systems with good accuracy.
Collapse
Affiliation(s)
- Jiarong Peng
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
8
|
Herbert JM. Dielectric continuum methods for quantum chemistry. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1519] [Citation(s) in RCA: 37] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- John M. Herbert
- Department of Chemistry and Biochemistry The Ohio State University Columbus Ohio USA
| |
Collapse
|
9
|
Malloum A, Fifen JJ, Conradie J. Determination of the absolute solvation free energy and enthalpy of the proton in solutions. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2020.114919] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
|
10
|
Itkis D, Cavallo L, Yashina LV, Minenkov Y. Ambiguities in solvation free energies from cluster-continuum quasichemical theory: lithium cation in protic and aprotic solvents. Phys Chem Chem Phys 2021; 23:16077-16088. [PMID: 34291782 DOI: 10.1039/d1cp01454d] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Gibbs free energies for Li+ solvation in water, methanol, acetonitrile, DMSO, dimethylacetamide, dimethoxyethane, dimethylformamide, gamma-butyrolactone, pyridine, and sulfolane have been calculated using the cluster-continuum quasichemical theory. With n independent solvent molecules S initial state forming the "monomer" thermodynamic cycle, Li+ solvation free energies are found to be on average 14 kcal mol-1 more positive compared to those from the "cluster" thermodynamic cycle where the initial state is the cluster Sn. We ascribe the inconsistency between the "monomer" and "cluster" cycles mainly to the incorrectly predicted solvation free energies of solvent clusters Sn from the SMD and CPCM continuum solvation models, which is in line with the earlier study of Bryantsev et al., J. Phys. Chem. B, 2008, 112, 9709-9719. When experimental-based solvation free energies of individual solvent molecules and solvent clusters are employed, the "monomer" and "cluster" cycles result in identical numbers. The best overall agreement with experimental-based "bulk" scale lithium cation solvation free energies was obtained for the "monomer" scale, and we recommend this set of values. We expect that further progress in the field is possible if (i) consensus on the accuracy of experimental reference values is achieved; (ii) the most recent continuum solvation models are properly parameterized for all solute-solvent combinations and become widely accessible for testing.
Collapse
Affiliation(s)
- Daniil Itkis
- N.N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russia.
| | | | | | | |
Collapse
|
11
|
Voityuk AA, Vyboishchikov SF. Fast and accurate calculation of hydration energies of molecules and ions. Phys Chem Chem Phys 2020; 22:14591-14598. [DOI: 10.1039/d0cp02667k] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
An accurate and efficient method for calculation of hydration free energy of ions and neutral molecules is presented.
Collapse
Affiliation(s)
- Alexander A. Voityuk
- Institució Catalana de Recerca i Estudis Avançats (ICREA)
- 08010 Barcelona
- Spain
- Institut de Química Computacional i Catàlisi and Departament de Química
- Universitat de Girona
| | - Sergei F. Vyboishchikov
- Institut de Química Computacional i Catàlisi and Departament de Química
- Universitat de Girona
- 17003 Girona
- Spain
- Peoples’ Friendship University of Russia (RUDN University)
| |
Collapse
|
12
|
van der Spoel D, Manzetti S, Zhang H, Klamt A. Prediction of Partition Coefficients of Environmental Toxins Using Computational Chemistry Methods. ACS OMEGA 2019; 4:13772-13781. [PMID: 31497695 PMCID: PMC6713992 DOI: 10.1021/acsomega.9b01277] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/03/2019] [Accepted: 06/27/2019] [Indexed: 05/05/2023]
Abstract
The partitioning of compounds between aqueous and other phases is important for predicting toxicity. Although thousands of octanol-water partition coefficients have been measured, these represent only a small fraction of the anthropogenic compounds present in the environment. The octanol phase is often taken to be a mimic of the inner parts of phospholipid membranes. However, the core of such membranes is typically more hydrophobic than octanol, and other partition coefficients with other compounds may give complementary information. Although a number of (cheap) empirical methods exist to compute octanol-water (log k OW) and hexadecane-water (log k HW) partition coefficients, it would be interesting to know whether physics-based models can predict these crucial values more accurately. Here, we have computed log k OW and log k HW for 133 compounds from seven different pollutant categories as well as a control group using the solvation model based on electronic density (SMD) protocol based on Hartree-Fock (HF) or density functional theory (DFT) and the COSMO-RS method. For comparison, XlogP3 (log k OW) values were retrieved from the PubChem database, and KowWin log k OW values were determined as well. For 24 of these compounds, log k OW was computed using potential of mean force (PMF) calculations based on classical molecular dynamics simulations. A comparison of the accuracy of the methods shows that COSMO-RS, KowWin, and XlogP3 all have a root-mean-square deviation (rmsd) from the experimental data of ≈0.4 log units, whereas the SMD protocol has an rmsd of 1.0 log units using HF and 0.9 using DFT. PMF calculations yield the poorest accuracy (rmsd = 1.1 log units). Thirty-six out of 133 calculations are for compounds without known log k OW, and for these, we provide what we consider a robust prediction, in the sense that there are few outliers, by averaging over the methods. The results supplied may be instrumental when developing new methods in computational ecotoxicity. The log k HW values are found to be strongly correlated to log k OW for most compounds.
Collapse
Affiliation(s)
- David van der Spoel
- Uppsala Center for
Computational Chemistry, Science for Life Laboratory, Department of
Cell and Molecular Biology, Uppsala University, Husargatan 3, Box
596, SE-75124 Uppsala, Sweden
- E-mail: . Phone: +46 18 4714205
| | - Sergio Manzetti
- Uppsala Center for
Computational Chemistry, Science for Life Laboratory, Department of
Cell and Molecular Biology, Uppsala University, Husargatan 3, Box
596, SE-75124 Uppsala, Sweden
- Fjordforsk A.S., Institute
for Science and Technology, Midtun, 6894 Vangsnes, Norway
| | - Haiyang Zhang
- Department of Biological Science and Engineering,
School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Andreas Klamt
- COSMOlogic GmbH & Co. KG, Imbacher Weg 46, D-51379 Leverkusen, Germany
- Institute of Physical and Theoretical Chemistry, University of Regensburg, 93053 Regensburg, Germany
| |
Collapse
|
13
|
Pliego JR, Riveros JM. Hybrid discrete‐continuum solvation methods. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1440] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Josefredo R. Pliego
- Departamento de Ciências Naturais Universidade Federal de São João del‐Rei São João del‐Rei Brazil
| | - Jose M. Riveros
- Departamento de Química Fundamental Instituto de Química, Universidade de São Paulo São Paulo Brazil
| |
Collapse
|
14
|
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
|
15
|
Zhang H, Jiang Y, Cui Z, Yin C. Force Field Benchmark of Amino Acids. 2. Partition Coefficients between Water and Organic Solvents. J Chem Inf Model 2018; 58:1669-1681. [DOI: 10.1021/acs.jcim.8b00493] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Beijing Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Box 53, 100029 Beijing, China
| | - Ziheng Cui
- Beijing Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Box 53, 100029 Beijing, China
| | - Chunhua Yin
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
16
|
Zhang H, Yin C, Jiang Y, van der Spoel D. Force Field Benchmark of Amino Acids: I. Hydration and Diffusion in Different Water Models. J Chem Inf Model 2018; 58:1037-1052. [DOI: 10.1021/acs.jcim.8b00026] [Citation(s) in RCA: 61] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Chunhua Yin
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Yang Jiang
- Beijing Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Box 53, Beijing 100029, China
| | - David van der Spoel
- Uppsala Center for Computational Chemistry, Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box
596, SE-75124 Uppsala, Sweden
| |
Collapse
|