1
|
Ferraz-Caetano J, Teixeira F, Cordeiro MNDS. Explainable Supervised Machine Learning Model To Predict Solvation Gibbs Energy. J Chem Inf Model 2024; 64:2250-2262. [PMID: 37603608 PMCID: PMC11005042 DOI: 10.1021/acs.jcim.3c00544] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2023] [Indexed: 08/23/2023]
Abstract
Many challenges persist in developing accurate computational models for predicting solvation free energy (ΔGsol). Despite recent developments in Machine Learning (ML) methodologies that outperformed traditional quantum mechanical models, several issues remain concerning explanatory insights for broad chemical predictions with an acceptable speed-accuracy trade-off. To overcome this, we present a novel supervised ML model to predict the ΔGsol for an array of solvent-solute pairs. Using two different ensemble regressor algorithms, we made fast and accurate property predictions using open-source chemical features, encoding complex electronic, structural, and surface area descriptors for every solvent and solute. By integrating molecular properties and chemical interaction features, we have analyzed individual descriptor importance and optimized our model though explanatory information form feature groups. On aqueous and organic solvent databases, ML models revealed the predictive relevance of solutes with increasing polar surface area and decreasing polarizability, yielding better results than state-of-the-art benchmark Neural Network methods (without complex quantum mechanical or molecular dynamic simulations). Both algorithms successfully outperformed previous ΔGsol predictions methods, with a maximum absolute error of 0.22 ± 0.02 kcal mol-1, further validated in an external benchmark database and with solvent hold-out tests. With these explanatory and statistical insights, they allow a thoughtful application of this method for predicting other thermodynamic properties, stressing the relevance of ML modeling for further complex computational chemistry problems.
Collapse
Affiliation(s)
- José Ferraz-Caetano
- Department
of Chemistry and Biochemistry − Faculty of Sciences, University of Porto - Rua do Campo Alegre, S/N, 4169-007 Porto, Portugal
| | - Filipe Teixeira
- Centre
of Chemistry, University of Minho, Campus
de Gualtar, 4710-057 Braga, Portugal
| | - M. Natália D. S. Cordeiro
- Department
of Chemistry and Biochemistry − Faculty of Sciences, University of Porto - Rua do Campo Alegre, S/N, 4169-007 Porto, Portugal
| |
Collapse
|
2
|
Xu M, Wu Z. Preparation of small peptoid-modified zwitterion-exchange/reversed-phase mixed-mode chromatography stationary phases for the separation of arenes and antibiotics. Anal Chim Acta 2024; 1296:342335. [PMID: 38401942 DOI: 10.1016/j.aca.2024.342335] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Revised: 01/19/2024] [Accepted: 02/04/2024] [Indexed: 02/26/2024]
Abstract
In this study, three small peptoids with different structures, named Sil-peptoids, were developed to improve the separation selectivity of zwitterion-exchange/reversed-phase mixed-mode chromatography stationary phases for multi-component complex drugs. Nonpolar, amphoteric, and alkaline drugs were used as test samples to demonstrate their retention behaviors in reversed-phase, ionic, and mixed-mode interactions. It was observed that different carboxyl anions in the small peptoids of the Sil-peptoids had vast differences in their stereo-selectivity. The stereo-selectivity and the influence of Sil-peptoids on the retention behavior of complex drugs and their interaction mechanism for the drug molecules were effectively evaluated through the combination of chromatographic analysis and molecular modeling. Finally, a mixture of drugs consisting of two polar and six non-polar drugs was used to obtain a separation effect with a resolution >1.5. Two other groups of polar antibiotics were used to verify the separation ability of the Sil-peptoids. The results indicated that the Sil-peptoids could separate multiple substances simultaneously. These novel stationary phases can be applied to the analysis of complex multi-component drugs.
Collapse
Affiliation(s)
- Meng Xu
- School of Pharmaceutical Science & Technology, Tianjin University, 92 Weijin Road, Tianjin, 300072, China.
| | - Zhenwei Wu
- School of Pharmaceutical Science & Technology, Tianjin University, 92 Weijin Road, Tianjin, 300072, China
| |
Collapse
|
3
|
Ghosh S, Hassan SH, Das A. Role of Explicit Solvation in Computational Modeling of Chemical Reactions: Mechanism of Cu(I) Transfer Between Thiolate-Based Chelators in Water. J Phys Chem B 2024. [PMID: 38503566 DOI: 10.1021/acs.jpcb.3c07327] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/21/2024]
Abstract
Solvation plays important roles in controlling the thermodynamic and kinetic aspects of chemical reactions. The conventional approaches to treat solvation via electronic structure methods are likely to become inadequate, when the reacting solutes have strong electrostatic and hydrogen bonding interactions with the solvent and undergo significant structural changes during the course of the reaction. In this article, we present evidence of such solvent and structural effects in the computational study of the Cu(I) transfer reaction between thiolate-based chelators dithiobutylamine (DTBA) and dithiotheritol (DTT) in water, inspired from biological copper trafficking phenomena. We propose a general solution to the problem by combining classical molecular dynamics (MD) simulations of the bulk system and static quantum chemistry calculations. The fluctuating solvation shell was estimated from MD, and energetics was assessed by averaging QM energies of a series of molecular clusters constructed from the MD snapshots. Applying this approach, we propose a reaction pathway with estimates of relative intermediate stabilities and barriers, which suggest the overall reaction to be reversible in nature and likely to go through both two and three coordinated intermediates, confirming previous studies of similar protein analogues. An interesting fact that emerged from our study was the strong indication that the rate-determining step is the deprotonation of initial thiol bound Cu(I) complex, without involving any Cu(I)-S bonds. The proposed method will lead to a better treatment of solvations, and these mechanistic insights will aid our understanding of biological copper(I) trafficking.
Collapse
Affiliation(s)
- Soumak Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Kolkata 700032, India
| | - Sk Hasibo Hassan
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Kolkata 700032, India
| | - Avisek Das
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Kolkata 700032, India
| |
Collapse
|
4
|
Šebesta F, Sovová Ž, Burda JV. Determination of Amino Acids' p Ka: Importance of Cavity Scaling within Implicit Solvation Models and Choice of DFT Functionals. J Phys Chem B 2024; 128:1627-1637. [PMID: 38345944 PMCID: PMC10895671 DOI: 10.1021/acs.jpcb.3c07007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Revised: 01/25/2024] [Accepted: 01/26/2024] [Indexed: 02/23/2024]
Abstract
Protonation states of molecules significantly influence the thermodynamics and kinetics of chemical reactions. This is especially important in biochemical processes, where appropriate protonation states of amino acids control the exo/endoergicity of practically all biochemical cycles. This paper is focused on appraisal of the impact of DFT functionals and PCM solvation models on the accuracy of pKa evaluations for all proteinogenic amino acids. Eight functionals (B3LYP, PBE0, revPBE0, M06-2X, M11, M11-L, TPSSh, and ωB97X-D) and four basis sets are considered, together with four kinds of implicit solvation models when additional attention is paid to a cavity construction. An influence of nonelectrostatic contributions and Wertz's corrections on Gibbs free energy is investigated together with accuracy of provided proton solvation energy. The best model is based on the M06-2X/6-311++G**/D-PCM/UAKS computational level. The fitting procedure is utilized to improve the accuracy of the evaluated models. All of these results are also compared with values obtained from the COSMOtherm program and CCSD(T) calculations. Results for cysteine and histidine are discussed individually, as they can be found in different protonation states at neutral pH.
Collapse
Affiliation(s)
- Filip Šebesta
- Department of Chemical Physics
and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16 Prague 2, Czech Republic
| | - Žofie Sovová
- Department of Chemical Physics
and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16 Prague 2, Czech Republic
| | - Jaroslav V. Burda
- Department of Chemical Physics
and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16 Prague 2, Czech Republic
| |
Collapse
|
5
|
Boukar O, Fifen JJ, Conradie J, Conradie MM. Solvation energies of the ferrous ion in water and in ammonia at various temperatures. J Mol Model 2024; 30:52. [PMID: 38285315 DOI: 10.1007/s00894-024-05839-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Accepted: 01/09/2024] [Indexed: 01/30/2024]
Abstract
CONTEXT The solvation of metal ions is crucial to understanding relevant properties in physics, chemistry, or biology. Therefore, we present solvation enthalpies and solvation free energies of the ferrous ion in water and ammonia. Our results agree well with the experimental reports for the hydration free energy and hydration enthalpy. We obtained [Formula: see text] kJ mol[Formula: see text] for the hydration free energy and [Formula: see text] kJ mol[Formula: see text] for the hydration enthalpy of ferrous ion in water at room temperature. At ambient temperature, we obtained [Formula: see text] kJ mol[Formula: see text] as the [Formula: see text] ammoniation free energy and [Formula: see text] kJ mol[Formula: see text] for the ammoniation enthalpy. In addition, the free energy of solvation is deeply affected when the temperature increases. This pattern can be attributed to the rise of entropy when the temperature rises. Besides, the temperature does not affect the ammoniation enthalpies and the hydration enthalpy of the [Formula: see text] ion. METHOD All the geometry optimizations are performed at the MP2 methods associated with the 6-31++g(d,p) basis set of Pople. solvated phase structures of [Formula: see text] ion in water or in ammonia are performed using the PCM model. The [Formula: see text] program suite was used to perform all the calculations. The program TEMPO was also used to evaluate the temperature sensitivity of the different obtained geometries.
Collapse
Affiliation(s)
- Ousman Boukar
- Department of Physics, Faculty of Science, University of Maroua, P.O. Box 46, Maroua, Cameroon.
| | - Jean Jules Fifen
- Quantum Theory and Aplications Unit, Department of Physics, Faculty of Science, The University of Ngaoundere, P.O. Box 454, Ngaoundere, Cameroon.
| | - Jeanet Conradie
- Department of Chemistry, University of the Free State, PO Box 339, Bloemfontein, 9300, South Africa
| | | |
Collapse
|
6
|
Liu J, Wan J, Ren Y, Shao X, Xu X, Rao L. DOX_BDW: Incorporating Solvation and Desolvation Effects of Cavity Water into Nonfitting Protein-Ligand Binding Affinity Prediction. J Chem Inf Model 2023; 63:4850-4863. [PMID: 37539963 DOI: 10.1021/acs.jcim.3c00776] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/05/2023]
Abstract
Accurate prediction of the protein-ligand binding affinity (PLBA) with an affordable cost is one of the ultimate goals in the field of structure-based drug design (SBDD), as well as a great challenge in the computational and theoretical chemistry. Herein, we have systematically addressed the complicated solvation and desolvation effects on the PLBA brought by the difference of the explicit water in the protein cavity before and after ligands bind to the protein-binding site. Based on the new solvation model, a nonfitting method at the first-principles level for the PLBA prediction was developed by taking the bridging and displaced water (BDW) molecules into account simultaneously. The newly developed method, DOX_BDW, was validated against a total of 765 noncovalent and covalent protein-ligand binding pairs, including the CASF2016 core set, Cov_2022 covalent binding testing set, and six testing sets for the hit and lead compound optimization (HLO) simulation. In all of the testing sets, the DOX_BDW method was able to produce PLBA predictions that were strongly correlated with the corresponding experimental data (R = 0.66-0.85). The overall performance of DOX_BDW is better than the current empirical scoring functions that are heavily parameterized. DOX_BDW is particularly outstanding for the covalent binding situation, implying the need for considering an electronic structure in covalent drug design. Furthermore, the method is especially recommended to be used in the HLO scenario of SBDD, where hundreds of similar derivatives need to be screened and refined. The computational cost of DOX_BDW is affordable, and its accuracy is remarkable.
Collapse
Affiliation(s)
- Jiaqi Liu
- National Key Laboratory of Green Pesticide, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, College of Chemistry, Central China Normal University, Wuhan 43009, People's Republic of China
| | - Jian Wan
- National Key Laboratory of Green Pesticide, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, College of Chemistry, Central China Normal University, Wuhan 43009, People's Republic of China
| | - Yanliang Ren
- National Key Laboratory of Green Pesticide, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, College of Chemistry, Central China Normal University, Wuhan 43009, People's Republic of China
| | - Xubo Shao
- National Key Laboratory of Green Pesticide, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, College of Chemistry, Central China Normal University, Wuhan 43009, People's Republic of China
| | - Xin Xu
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education (MOE) Laboratory for Computational Physical Science, Department of Chemistry, Fudan University, Shanghai 200433, People's Republic of China
| | - Li Rao
- National Key Laboratory of Green Pesticide, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, College of Chemistry, Central China Normal University, Wuhan 43009, People's Republic of China
| |
Collapse
|
7
|
Zhu Q, Costentin C, Stubbe J, Nocera DG. Disulfide radical anion as a super-reductant in biology and photoredox chemistry. Chem Sci 2023; 14:6876-6881. [PMID: 37389245 PMCID: PMC10306091 DOI: 10.1039/d3sc01867a] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2023] [Accepted: 05/17/2023] [Indexed: 07/01/2023] Open
Abstract
Disulfides are involved in a broad range of radical-based synthetic organic and biochemical transformations. In particular, the reduction of a disulfide to the corresponding radical anion, followed by S-S bond cleavage to yield a thiyl radical and a thiolate anion plays critical roles in radical-based photoredox transformations and the disulfide radical anion in conjunction with a proton donor, mediates the enzymatic synthesis of deoxynucleotides from nucleotides within the active site of the enzyme, ribonucleotide reductase (RNR). To gain fundamental thermodynamic insight into these reactions, we have performed experimental measurements to furnish the transfer coefficient from which the standard E0(RSSR/RSSR˙-) reduction potential has been determined for a homologous series of disulfides. The electrochemical potentials are found to be strongly dependent on the structures and electronic properties of the substituents of the disulfides. In the case of cysteine, a standard potential of E0(RSSR/RSSR˙-) = -1.38 V vs. NHE is determined, making the disulfide radical anion of cysteine one of the most reducing cofactors in biology.
Collapse
Affiliation(s)
- Qilei Zhu
- Department of Chemistry and Chemical Biology, Harvard University 12 Oxford Street Cambridge Massachusetts 02138 USA
- Department of Chemistry, University of Utah Salt Lake City Utah 84112 USA
| | | | - JoAnne Stubbe
- Department of Chemistry and Chemical Biology, Harvard University 12 Oxford Street Cambridge Massachusetts 02138 USA
- Departments of Chemistry and Department of Biology, Massachusetts Institute of Technology Cambridge Massachusetts 02139 USA
| | - Daniel G Nocera
- Department of Chemistry and Chemical Biology, Harvard University 12 Oxford Street Cambridge Massachusetts 02138 USA
| |
Collapse
|
8
|
Elder J, Broome JA, Bushnell EAC. Computational Insights into the Regeneration of Ovothiol and Ergothioneine and Their Selenium Analogues by Glutathione. ACS OMEGA 2022; 7:31813-31821. [PMID: 36120043 PMCID: PMC9476190 DOI: 10.1021/acsomega.2c02506] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/21/2022] [Accepted: 07/18/2022] [Indexed: 06/15/2023]
Abstract
Ovothiol and ergothioneine are powerful antioxidants that readily react with oxidants by forming their respective disulfides. In fact, ovothiol is widely considered one of the most powerful natural antioxidants. However, for these antioxidants to be again involved in reacting with oxidants, they must be regenerated via the reduction of the disulfide bonds. In the present work, the regeneration of the antioxidants ovothiol and ergothioneine and their selenium analogues, by the closed-shell nucleophilic attack of glutathione, was investigated using density functional theory. From the calculated thermodynamic data, the attack of glutathione on OSSO and EYYE (where Y = S and/or Se) will readily occur in solution. Moreover, in comparison to the reference reaction GSH + GSSG → GSSG + GSH, all reactions are expected to be faster. Overall, the results presented herein show that the key antioxidant GSH should readily recycle ovothiol, ovoselenol, ergothioneine, and ergoseloneine from OYYO and EYYE (where Y = S and/or Se).
Collapse
|
9
|
Gupta M, da Silva EF, Svendsen HF. Modeling Differential Enthalpy of Absorption of CO 2 with Piperazine as a Function of Temperature. J Phys Chem B 2022; 126:1980-1991. [PMID: 35226495 PMCID: PMC8919260 DOI: 10.1021/acs.jpcb.1c10755] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Temperature-dependent correlations for equilibrium constants (ln K) and heat of absorption (ΔHabs) of different reactions (i.e., deprotonation, double deprotonation, carbamate formation, protonated carbamate formation, dicarbamate formation) involved in the piperazine (PZ)/CO2/H2O system have been calculated using computational chemistry based ln K values input to the Gibbs-Helmholtz equation. This work also presents an extensive study of gaseous phase free energy and enthalpy for different reactions using composite (G3MP2B3, G3MP2, CBS-QB3, and G4MP2) and density functional theory [B3LYP/6-311++G(d,p)] methods. The explicit solvation shell (ESS) model and SM8T solvation free energy coupled with gaseous phase density functional theory calculations give temperature-dependent reaction equilibrium constants for different reactions. Calculated individual and overall reaction equilibrium constants and enthalpies of different reactions involved in CO2 absorption in piperazine solution are compared against experimental data, where available, in the temperature range 273.15-373 K. Postcombustion CO2 capture (PCC) is a temperature swing absorption-desorption process. The enthalpy of the solution directly correlates with the steam requirement of the amine regeneration step. Temperature-dependent correlations for ln K and ΔHabs calculated using computational chemistry tools can help evaluate potential PCC solvents' thermodynamics and cost-efficiency. These correlations can also be employed in thermodynamic models (e.g., e-UNIQUAC, e-NRTL) to better understand postcombustion CO2 capture solvent chemistry.
Collapse
Affiliation(s)
- Mayuri Gupta
- Department of Chemical Engineering, Norwegian University of Science and Technology, Sem Sælands vei 4, Trondheim 7491, Norway
| | | | - Hallvard F Svendsen
- Department of Chemical Engineering, Norwegian University of Science and Technology, Sem Sælands vei 4, Trondheim 7491, Norway
| |
Collapse
|
10
|
Gupta M, Weaver DF. Microsecond molecular dynamics studies of cholesterol-mediated myelin sheath degeneration in early Alzheimer's disease. Phys Chem Chem Phys 2021; 24:222-239. [PMID: 34878462 DOI: 10.1039/d1cp03844c] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Cholesterol-mediated perturbations of membrane structural integrity are key early events in the molecular pathogenesis of Alzheimer's disease (AD). In AD, protein misfolding (proteopathy) and pro-inflammatory conditions (immunopathy) culminate in neuronal death, a process enabled by altered membrane biophysical properties which render neurons more susceptible to proteopathic and immunopathic cytotoxicities. Since cholesterol is a principal neuronal membrane lipid, normal cholesterol homeostasis is central to membrane health; also, since increased cholesterol composition is especially present in neuronal myelin sheath (i.e. brain "white matter"), recent studies have not surprisingly revealed that white matter atrophy precedes the conventional biomarkers of AD (amyloid plaques, tau tangles). Employing extensive microsecond all-atom molecular dynamics simulations, we investigated biophysical and mechanical properties of myelin sheath membrane as a function of cholesterol mole fraction (χCHL). Impaired χCHL modulates multiple bilayer properties, including surface area per lipid (APL), chain order, number and mass density profiles, area compressibility and bending moduli, bilayer thickness, lipid tilt angles, H-bonding interactions and tail interdigitation. The increased orientational ordering of both palmitoyl and oleoyl chains in model healthy myelin sheath (HMS) membranes illustrates the condensing effect of cholesterol. With an increase in χCHL, number density profiles of water tend to attain bulk water number density more quickly, indicating shrinkage in the interfacial region with increasing χCHL. The average tilt value is 11.5° for the C10-C13 angle in cholesterol and 64.2° for the P-N angle in POPC lipids in HMS. These calculations provide a molecular-level understanding of myelin sheath susceptibility to pathology as an early event in the pathogenesis of AD.
Collapse
Affiliation(s)
- Mayuri Gupta
- Krembil Research Institute, University Health Network, 60 Leonard Avenue, Toronto, M5T 0S8, Canada.
| | - Donald F Weaver
- Krembil Research Institute, University Health Network, 60 Leonard Avenue, Toronto, M5T 0S8, Canada. .,Departments of Chemistry, Medicine and Pharmaceutical Sciences, University of Toronto, Toronto, ON, M55 3H6, Canada
| |
Collapse
|
11
|
Gupta M, Weaver DF. Axonal plasma membrane-mediated toxicity of cholesterol in Alzheimer's disease: A microsecond molecular dynamics study. Biophys Chem 2021; 281:106718. [PMID: 34808480 DOI: 10.1016/j.bpc.2021.106718] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2021] [Revised: 11/03/2021] [Accepted: 11/10/2021] [Indexed: 12/18/2022]
Abstract
Alzheimer's disease is increasingly being recognized as an immune-mediated disease of brain. Since physiological brain health and brain immune function is dependent upon homeostatic neuronal membrane structure and function, alterations in membrane lipid biochemistry may predispose to disease. Brain is rich in cholesterol, and cholesterol metabolism dysfunction is a known risk factor for AD. Employing extensive microsecond all-atom molecular dynamics simulations, we investigated the properties of model neuronal membranes as a function of cholesterol concentration; phospholipid and phospholipid/cholesterol bilayers were also simulated to compare against available experimental data. Increased cholesterol concentrations compact and stiffen the lipid membrane, reducing permeability while modulating local water densities in the peri-membranous environment. Conversely, lower cholesterol mole fraction yields membranes with increased molecular disorder, enhanced fluidity, higher molecular tilting, and augmented interdigitation between bilayer leaflet lipids. Our findings provide a molecular insight on effect of cholesterol composition on various biochemical processes occurring at neuronal axon plasma membrane. These calculations also endeavor to establish a membrane-based link between cholesterol as an AD risk factor and possible AD pathology.
Collapse
Affiliation(s)
- Mayuri Gupta
- Krembil Research Institute, University Health Network, 60 Leonard Avenue, Toronto M5T 0S8, Canada
| | - Donald F Weaver
- Krembil Research Institute, University Health Network, 60 Leonard Avenue, Toronto M5T 0S8, Canada; Department of Chemistry, University of Toronto, Toronto M55 3H6, Canada; Department of Medicine, University of Toronto, Toronto M5G 2C4, Canada; Department of Pharmaceutical Sciences, University of Toronto, Toronto M5S 3M2, Canada.
| |
Collapse
|
12
|
Gupta M, Svendsen HF. Understanding Carbamate Formation Reaction Thermochemistry of Amino Acids as Solvents for Postcombustion CO 2 Capture. J Phys Chem B 2019; 123:8433-8447. [PMID: 31513397 DOI: 10.1021/acs.jpcb.9b06447] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The carbamate stability constant for a data set of 10 amino acids, having potential for being postcombustion CO2 capture (PCC) solvents, has been calculated using various implicit and explicit solvation shell models. This work also includes an extensive study of gas-phase free energy and enthalpy for the amino acid carbamate formation reaction with the Hartree Fock method, density functional methods [B3LYP/6-311++G(d,p)], and composite methods (G3MP2B3, G3MP2, CBS-QB3, and G4MP2). Ideal PCC solvent properties require finding a profitable tradeoff between various thermodynamic and system optimization parameters. Benchmark gaseous-phase and solution-phase thermodynamic properties given in this work can help in making informed decisions when choosing promising PCC solvents. The temperature dependency of the carbamate stability constant of amino acids is predicted using PCM and SM8T implicit solvation models. PCC is a temperature swing absorption-desorption process, and the high-temperature sensitivity of the ln KcAmCOO- value is of vital importance in attaining cost-efficient processes.
Collapse
Affiliation(s)
- Mayuri Gupta
- Department of Chemical Engineering , Norwegian University of Science and Technology , 7491 Trondheim , Norway
| | - Hallvard F Svendsen
- Department of Chemical Engineering , Norwegian University of Science and Technology , 7491 Trondheim , Norway
| |
Collapse
|
13
|
Baetzold RC. Density Functional Calculations for Aqueous Silver Clusters Containing Water and Nitrate Ligands. J Phys Chem B 2019; 123:8300-8312. [PMID: 31478667 DOI: 10.1021/acs.jpcb.9b06634] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The incorporation of a nitrate ion into silver-water aqueous clusters has been examined using PBE0 density functional theory with the solvent model density (SMD) solvation model. The Gibbs free energy of solvation and other thermodynamic variables are calculated using the harmonic/rigid rotor/ideal gas model at 298.15 K for aqueous solutes including the effects of solute relaxation in water and with London dispersive forces at the D3 level. Free energies of solvation for Ag+ and NO3- were found to agree well with experimental values of -118.2 and -60.1 kcal/mol, respectively, calculated using cluster-continuum models with six to eight water molecules and including solute relaxation and London D3 dispersive interactions. An analysis of data of varying cluster size upon calculated free energy is presented. A direct procedure is applied to aqueous clusters such as Agnz(NO3-)(H2O)5, Agnz(H2O)5, and (NO3-)(H2O)6 n = 1-4; z = 0, +1 in the SMD solvent representation to calculate equilibrium constants for nitrate association with silver clusters in solution that includes fully relaxed solutes. The equilibrium structures of the nitrate-containing clusters involve one or more bonds from nitrate oxygen to positive silver clusters. Water molecules interact with nitrate through H atoms, and overall, the structure represents a silver nitrate cluster with water ligands having similarity to a close ion pair in many aspects. The neutral silver atom is attached to nitrate through H-bonded water molecules. The ratio of nitrate-containing silver clusters to nitrate-free clusters using a calculated equilibrium constant of 0.51 L/mol for Ag+ is small in the range of many experiments. Similar values are found for positive silver clusters up to four atoms in size. The resulting procedures were applied to aqueous clusters of Agn(NO3)m+(n-m) that have been previously experimentally studied for silver reduction in aqueous solution. A chain-like structure with collinear and bidentate oxygen bonds to silver was found, and the equilibrium constants for clustering were determined. A simplified model calculation for the reduction of Ag(H2O)6+ clusters in the presence of silver clusters in aqueous media was studied to understand catalytic effects observed in these systems. The reduction potentials vary with silver cluster size indicating a more favorable reduction caused by the presence of larger silver clusters.
Collapse
|
14
|
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
|
15
|
Mirzaei S, Ivanov MV, Timerghazin QK. Improving Performance of the SMD Solvation Model: Bondi Radii Improve Predicted Aqueous Solvation Free Energies of Ions and pKa Values of Thiols. J Phys Chem A 2019; 123:9498-9504. [DOI: 10.1021/acs.jpca.9b02340] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Saber Mirzaei
- Department of Chemistry, Marquette University, Milwaukee, Wisconsin 53201-1414, United States
| | - Maxim V. Ivanov
- Department of Chemistry, Marquette University, Milwaukee, Wisconsin 53201-1414, United States
| | - Qadir K. Timerghazin
- Department of Chemistry, Marquette University, Milwaukee, Wisconsin 53201-1414, United States
| |
Collapse
|
16
|
Yang X, Zou J, Wang Y, Xue Y, Yang S. Role of Water in the Reaction Mechanism and endo/exo Selectivity of 1,3-Dipolar Cycloadditions Elucidated by Quantum Chemistry and Machine Learning. Chemistry 2019; 25:8289-8303. [PMID: 30887586 DOI: 10.1002/chem.201900617] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2019] [Indexed: 02/05/2023]
Abstract
Asymmetric 1,3-dipolar cycloadditions of azomethine ylides with activated olefins are among the most important and versatile methods for the synthesis of enantioenriched pyrroline and pyrrolidine derivatives. Despite both theoretical and practical importance, the role of water molecules in the reactivity and endo/exo selectivity remains unclear. To explore how water accelerates the reactions and improves the endo/exo selectivity of the cycloadditions of 1,3-dipole phthalazinium-2-dicyanomethanide (1) and two dipolarophiles, an ab initio-quality neural network potential that overcomes the computational bottleneck of explicitly considering water molecules was used. It is demonstrated that not only the nature of both the dipolarophile and the 1,3-dipole, but also the solvent medium, can perturb or even alter the reaction mechanism. An extreme case was found for the reaction of 1,3-dipole 1 with methyl vinyl ketone, in which the reaction mechanism changes from a concerted to a stepwise mode on going from MeCN to H2 O as solvent, with formation of a zwitterionic intermediate that is a very shallow minimum on the energy surface. Thus, high stereocontrol can still be expected despite the stepwise nature of the mechanism. The results indicate that water can induce global polarization along the reaction coordinate and highlight the role of microsolvation effects and bulk-phase effects in reproducing the experimentally observed aqueous acceleration and enhanced endo/exo selectivity.
Collapse
Affiliation(s)
- Xin Yang
- State Key Laboratory of Biotherapy and Cancer Center, West China Hospital, Sichuan University, Chengdu, Sichuan, 610041, P.R. China
| | - Jun Zou
- State Key Laboratory of Biotherapy and Cancer Center, West China Hospital, Sichuan University, Chengdu, Sichuan, 610041, P.R. China
| | - Yifei Wang
- State Key Laboratory of Biotherapy and Cancer Center, West China Hospital, Sichuan University, Chengdu, Sichuan, 610041, P.R. China
| | - Ying Xue
- College of Chemistry, Key Lab of Green Chemistry and Technology in Ministry of Education, Sichuan University, Chengdu, Sichuan, 610041, P.R. China
| | - Shengyong Yang
- State Key Laboratory of Biotherapy and Cancer Center, West China Hospital, Sichuan University, Chengdu, Sichuan, 610041, P.R. China
| |
Collapse
|
17
|
Arunachalam V, Tummanapelli AK, Vasudevan S. The multiple dissociation constants of glutathione disulfide: interpreting experimental pH-titration curves with ab initio MD simulations. Phys Chem Chem Phys 2019; 21:9212-9217. [PMID: 30993274 DOI: 10.1039/c9cp00761j] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The hexapeptide glutathione disulfide (GSSG) has six ionizable groups with six associated dissociation constants. The experimentally measured pH-titration curve, however, does not exhibit the six corresponding equivalence points and bears little resemblance to standard textbook examples of acid-base pH-titration curves. The curve highlights the difficulties in determining multiple pKa values of polyprotic acids - typically proteins and peptides - from experiment. The six pKa values of GSSG can, however, be estimated using Car-Parrinello molecular dynamics (CPMD) simulations in conjunction with metadynamics sampling of the underlying free energy landscape of the dissociation reactions. Ab initio MD simulations were performed on a GSSG molecule solvated by 200 water molecules. Using the estimated pKa values the theoretical titration curve was calculated and found to be in good agreement with experiment. The results clearly highlight how dissociation constants estimated from ab initio MD simulations can facilitate the interpretation of the pH-titration curves of complex chemical and biological systems.
Collapse
|
18
|
Roy D, Kovalenko A. Performance of 3D-RISM-KH in Predicting Hydration Free Energy: Effect of Solute Parameters. J Phys Chem A 2019; 123:4087-4093. [PMID: 30993994 DOI: 10.1021/acs.jpca.9b01623] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The three-dimensional reference interaction site model molecular solvation theory with the Kovalenko-Hirata closure relation has been shown to produce excellent solvation characteristics for a large class of (bio)chemical systems in solution. Correct calculation of hydration free energy is central to successful application of any solvation model. In order to find out the best possible force-field parameters to be used for hydration free energy calculation with the aforementioned theory, we have developed an extended database containing a large number of experimental solvation free energies available in the current literature and used a plethora of theoretical models for assessment. The general Amber force field was found to perform satisfactorily, whereas special care should be taken in solute charge assignment with the universal force field.
Collapse
Affiliation(s)
- Dipankar Roy
- Department of Mechanical Engineering , University of Alberta , 10-203 Donadeo Innovation Centre for Engineering, 9211-116 Street NW , Edmonton , Alberta T6G 1H9 , Canada
| | - Andriy Kovalenko
- Department of Mechanical Engineering , University of Alberta , 10-203 Donadeo Innovation Centre for Engineering, 9211-116 Street NW , Edmonton , Alberta T6G 1H9 , Canada.,Nanotechnology Research Centre , 11421 Saskatchewan Drive , Edmonton , Alberta T6G 2M9 , Canada
| |
Collapse
|
19
|
Zheng Y, Zheng W, Zhu D, Chang H. Theoretical modeling of pKa's of thiol compounds in aqueous solution. NEW J CHEM 2019. [DOI: 10.1039/c8nj06259e] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The pKa's of different kinds of thiols (R-SH) were investigated by using the M06-2X method with a SMDsSAS model.
Collapse
Affiliation(s)
- Yuanyuan Zheng
- College of Chemistry and Chemical Engineering
- Shanghai University of Engineering Science
- Shanghai 201620
- China
| | - Wenrui Zheng
- College of Chemistry and Chemical Engineering
- Shanghai University of Engineering Science
- Shanghai 201620
- China
| | - Danfeng Zhu
- College of Chemistry and Chemical Engineering
- Shanghai University of Engineering Science
- Shanghai 201620
- China
| | - Huifang Chang
- College of Chemistry and Chemical Engineering
- Shanghai University of Engineering Science
- Shanghai 201620
- China
| |
Collapse
|
20
|
Yang X, Rees RJ, Conway W, Puxty G, Yang Q, Winkler DA. Computational Modeling and Simulation of CO2 Capture by Aqueous Amines. Chem Rev 2017; 117:9524-9593. [PMID: 28517929 DOI: 10.1021/acs.chemrev.6b00662] [Citation(s) in RCA: 68] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Xin Yang
- CSIRO Manufacturing, Bayview Avenue, Clayton 3169, Australia
- College
of Chemistry, Key Lab of Green Chemistry and Technology in Ministry
of Education, Sichuan University, Chengdu 610064, People’s Republic of China
| | - Robert J. Rees
- Data61
- CSIRO, Door 34 Goods
Shed, Village Street, Docklands VIC 3008, Australia
| | | | | | - Qi Yang
- CSIRO Manufacturing, Bayview Avenue, Clayton 3169, Australia
| | - David A. Winkler
- CSIRO Manufacturing, Bayview Avenue, Clayton 3169, Australia
- Monash Institute of Pharmaceutical Sciences, 392 Royal Parade, Parkville 3052, Australia
- Latrobe Institute for Molecular Science, Bundoora 3046, Australia
- School
of
Chemical and Physical Science, Flinders University, Bedford Park 5042, Australia
| |
Collapse
|
21
|
Sastre S, Casasnovas R, Muñoz F, Frau J. Isodesmic reaction for accurate theoretical pKa calculations of amino acids and peptides. Phys Chem Chem Phys 2017; 18:11202-12. [PMID: 27052591 DOI: 10.1039/c5cp07053h] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Theoretical and quantitative prediction of pKa values at low computational cost is a current challenge in computational chemistry. We report that the isodesmic reaction scheme provides semi-quantitative predictions (i.e. mean absolute errors of 0.5-1.0 pKa unit) for the pKa1 (α-carboxyl), pKa2 (α-amino) and pKa3 (sidechain groups) of a broad set of amino acids and peptides. This method fills the gaps of thermodynamic cycles for the computational pKa calculation of molecules that are unstable in the gas phase or undergo proton transfer reactions or large conformational changes from solution to the gas phase. We also report the key criteria to choose a reference species to make accurate predictions. This method is computationally inexpensive and makes use of standard density functional theory (DFT) and continuum solvent models. It is also conceptually simple and easy to use for researchers not specialized in theoretical chemistry methods.
Collapse
Affiliation(s)
- S Sastre
- Department de Química, Universitat de les Illes Balears, Palma de Mallorca 07122, Spain.
| | - R Casasnovas
- Institute of Neuroscience and Medicine INM-9, Institute for Advanced Simulations IAS-5. Computational Biomedicine, Forschungszentrum Jülich 52428, Germany
| | - F Muñoz
- Department de Química, Universitat de les Illes Balears, Palma de Mallorca 07122, Spain. and Instituto de Investigación Sanitaria de Palma (IdISPa), 07010, Palma, Spain
| | - J Frau
- Department de Química, Universitat de les Illes Balears, Palma de Mallorca 07122, Spain. and Instituto de Investigación Sanitaria de Palma (IdISPa), 07010, Palma, Spain
| |
Collapse
|
22
|
Kannan R, Muthuvijayan V, Prasad E. In vitro study of a glucose attached poly(aryl ether) dendron based gel as a drug carrier for a local anaesthetic. NEW J CHEM 2017. [DOI: 10.1039/c7nj01420a] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A poly(aryl ether) dendron based gelator as an efficient transdermal drug carrier for the controlled release of prilocaine hydrochloride.
Collapse
Affiliation(s)
- Ramya Kannan
- Department of Chemistry
- Indian Institute of Technology
- Chennai
- India
- Department of Biotechnology
| | - Vignesh Muthuvijayan
- Department of Biotechnology
- Bhupat and Jyoti Mehta School of Biosciences
- Indian Institute of Technology
- Chennai
- India
| | - Edamana Prasad
- Department of Chemistry
- Indian Institute of Technology
- Chennai
- India
| |
Collapse
|
23
|
Satapathy S, Prasad E. Charge Transfer Modulated Self-Assembly in Poly(aryl ether) Dendron Derivatives with Improved Stability and Transport Characteristics. ACS APPLIED MATERIALS & INTERFACES 2016; 8:26176-26189. [PMID: 27603727 DOI: 10.1021/acsami.6b09175] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Alteration of native gelation properties of anthracene and pyrene cored first generation poly(aryl ether) dendrons, G1-An and G1-Py, by introducing a common acceptor, 2,4,7-trinitro-9H-fluoren-9-one (TNF), results in forming charge transfer gels in long chain alcoholic solvents. This strategy leads to significant perturbation of optical and electronic properties within the gel matrix. Consequently, a noticeable increase of their electrical conductivities is observed, making these poly(aryl ether) dendron based gels potential candidates for organic electronics. While the dc-conductivity (σ) value for the native gel from G1-An is 2.8 × 10-4 S m-1, the value increased 3 times (σ = 8.7 × 10-4 S m-1) for its corresponding charge transfer gel. Further, the dc-conductivity for the native gel self-assembled from G1-Py dramatically enhanced by approximately an order of magnitude from 4.9 × 10-4 to 1.3 × 10-3 S m-1, under the influence of an acceptor. Apart from H-bonding and π···π interactions, charge transfer results in the formation of a robust 3D network of fibers, with improved aspect ratio, providing high thermo-mechanical stability to the gels compared to the native ones. The charge transfer gels self-assembled from G1-An/TNF (1:1) and G1-Py/TNF exhibit a 7.3- and 2.5-fold increase in their yield stress, respectively, compared to their native assemblies. A similar trend follows in the case of their thermal stabilities. This is attributed to the typical bilayer self-assembly of the former which is not present in the case of G1-Py/TNF charge transfer gel. Density functional calculations provide deeper insights accounting for the role of charge transfer interactions in the mode of self-assembly. The 1D potential energy surface for the G1-An/TNF dimer and G1-Py/TNF dimer is found to be 11.8 and 1.9 kcal mol-1 more stable than their corresponding native gel dimers, G1-An/G1-An and G1-Py/G1-Py, respectively.
Collapse
Affiliation(s)
- Sitakanta Satapathy
- Department of Chemistry, Indian Institute of Technology Madras (IIT M) , Chennai 600 036, India
| | - Edamana Prasad
- Department of Chemistry, Indian Institute of Technology Madras (IIT M) , Chennai 600 036, India
| |
Collapse
|
24
|
Gupta M, da Silva EF, Svendsen HF. Postcombustion CO2 Capture Solvent Characterization Employing the Explicit Solvation Shell Model and Continuum Solvation Models. J Phys Chem B 2016; 120:9034-50. [DOI: 10.1021/acs.jpcb.6b04049] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Mayuri Gupta
- Department
of Chemical Engineering, Norwegian University of Science and Technology, Trondheim 7491, Norway
| | - Eirik F. da Silva
- Department
of Process Technology, SINTEF Materials and Chemistry, Trondheim 7034, Norway
| | - Hallvard F. Svendsen
- Department
of Chemical Engineering, Norwegian University of Science and Technology, Trondheim 7491, Norway
| |
Collapse
|
25
|
Thapa B, Schlegel HB. Density Functional Theory Calculation of pKa’s of Thiols in Aqueous Solution Using Explicit Water Molecules and the Polarizable Continuum Model. J Phys Chem A 2016; 120:5726-35. [DOI: 10.1021/acs.jpca.6b05040] [Citation(s) in RCA: 120] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Bishnu Thapa
- Chemistry Department, Wayne State University, Detroit, Michigan 48202, United States
| | - H. Bernhard Schlegel
- Chemistry Department, Wayne State University, Detroit, Michigan 48202, United States
| |
Collapse
|
26
|
Abstract
We recently showed for a large data set of pKas and reduction potentials that free energies calculated directly within the SMD continuum model compares very well with corresponding thermodynamic cycle calculations in both aqueous and organic solvents [ Phys. Chem. Chem. Phys. 2015 , 17 , 2859 ]. In this paper, we significantly expand the scope of our study to examine the suitability of this approach for calculating general solution phase kinetics and thermodynamics, in conjunction with several commonly used solvation models (SMD-M062X, SMD-HF, CPCM-UAKS, and CPCM-UAHF) for a broad range of systems. This includes cluster-continuum schemes for pKa calculations as well as various neutral, radical, and ionic reactions such as enolization, cycloaddition, hydrogen and chlorine atom transfer, and SN2 and E2 reactions. On the basis of this benchmarking study, we conclude that the accuracies of both approaches are generally very similar-the mean errors for Gibbs free energy changes of neutral and ionic reactions are approximately 5 and 25 kJ mol(-1), respectively. In systems where there are significant structural changes due to solvation, as is the case for certain ionic transition states and amino acids, the direct approach generally afford free energy changes that are in better agreement with experiment.
Collapse
Affiliation(s)
- Junming Ho
- Agency for Science, Technology and Research, Institute of High Performance Computing , 1 Fusionopolis Way, #16-16 Connexis, Singapore 138632.,Department of Chemistry, Yale University , P.O. Box 208107, New Haven, Connecticut 06520, United States
| | - Mehmed Z Ertem
- Chemistry Department, Brookhaven National Laboratory , Upton, New York 11973, United States
| |
Collapse
|
27
|
Tummanapelli AK, Vasudevan S. Ab Initio MD Simulations of the Brønsted Acidity of Glutathione in Aqueous Solutions: Predicting pKa Shifts of the Cysteine Residue. J Phys Chem B 2015; 119:15353-8. [PMID: 26550841 DOI: 10.1021/acs.jpcb.5b10093] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The tripeptide glutathione (GSH) is one of the most abundant peptides and the major repository for nonprotein sulfur in both animal and plant cells. It plays a critical role in intracellular oxidative stress management by the reversible formation of glutathione disulfide with the thiol-disulfide pair acting as a redox buffer. The state of charge of the ionizable groups of GSH can influence the redox couple, and hence the pKa value of the cysteine residue of GSH is critical to its functioning. Here we report ab initio Car-Parrinello molecular dynamics simulations of glutathione solvated by 200 water molecules, all of which are considered in the simulation. We show that the free-energy landscape for the protonation-deprotonation reaction of the cysteine residue of GSH computed using metadynamics sampling provides accurate estimates of the pKa and correctly predicts the shift in the dissociation constant values as compared with the isolated cysteine amino acid.
Collapse
Affiliation(s)
- Anil Kumar Tummanapelli
- Department of Inorganic and Physical Chemistry, Indian Institute of Science , Bangalore 560012, India
| | - Sukumaran Vasudevan
- Department of Inorganic and Physical Chemistry, Indian Institute of Science , Bangalore 560012, India
| |
Collapse
|
28
|
Tummanapelli AK, Vasudevan S. Ab Initio Molecular Dynamics Simulations of Amino Acids in Aqueous Solutions: Estimating pKa Values from Metadynamics Sampling. J Phys Chem B 2015; 119:12249-55. [PMID: 26331783 DOI: 10.1021/acs.jpcb.5b05211] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Changes in the protonation and deprotonation of amino acid residues in proteins play a key role in many biological processes and pathways. Here, we report calculations of the free-energy profile for the protonation-deprotonation reaction of the 20 canonical α amino acids in aqueous solutions using ab initio Car-Parrinello molecular dynamics simulations coupled with metadynamics sampling. We show here that the calculated change in free energy of the dissociation reaction provides estimates of the multiple pKa values of the amino acids that are in good agreement with experiment. We use the bond-length-dependent number of the protons coordinated to the hydroxyl oxygen of the carboxylic and the amine groups as the collective variables to explore the free-energy profiles of the Bronsted acid-base chemistry of amino acids in aqueous solutions. We ensure that the amino acid undergoing dissociation is solvated by at least three hydrations shells with all water molecules included in the simulations. The method works equally well for amino acids with neutral, acidic and basic side chains and provides estimates of the multiple pKa values with a mean relative error, with respect to experimental results, of 0.2 pKa units.
Collapse
Affiliation(s)
- Anil Kumar Tummanapelli
- Department of Inorganic and Physical Chemistry, Indian Institute of Science , Bangalore 560012, India
| | - Sukumaran Vasudevan
- Department of Inorganic and Physical Chemistry, Indian Institute of Science , Bangalore 560012, India
| |
Collapse
|
29
|
Thapa B, Schlegel HB. Calculations of pKa's and redox potentials of nucleobases with explicit waters and polarizable continuum solvation. J Phys Chem A 2014; 119:5134-44. [PMID: 25291241 DOI: 10.1021/jp5088866] [Citation(s) in RCA: 99] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The SMD implicit solvation model augmented with one and four explicit water molecules was used to calculate pKa's and redox potentials of N-methyl-substituted nucleic acid bases guanine, adenine, cytosine, thymine, and uracil. Calculations were carried out with the B3LYP/6-31+G(d,p) level of theory. The same numbers of water molecules were hydrogen bonded to the neutral, protonated, and deprotonated nucleobases in their unoxidized and oxidized forms. The improvement in pKa1 involving neutrals and cations was modest. By contrast, the improvement in pKa2 involving neutrals and anions was quite significant, reducing the mean absolute error from 4.6 pKa units with no waters, to 2.6 with one water and 1.7 with four waters. For the oxidation of nucleobases, adding explicit waters did little to improve E(X(•),H(+)/XH), possibly because both species in the redox couple are neutral molecules at pH 7.
Collapse
Affiliation(s)
- Bishnu Thapa
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - H Bernhard Schlegel
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| |
Collapse
|
30
|
Yang C, Xue XS, Li X, Cheng JP. Computational Study on the Acidic Constants of Chiral Brønsted Acids in Dimethyl Sulfoxide. J Org Chem 2014; 79:4340-51. [DOI: 10.1021/jo500158e] [Citation(s) in RCA: 63] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Chen Yang
- State Key
Laboratory of Elemento-Organic
Chemistry, Department of Chemistry, and Collaborative Innovation Center
of Chemical Science and Engineering, Nankai University, Tianjin 300071, China
| | - Xiao-Song Xue
- State Key
Laboratory of Elemento-Organic
Chemistry, Department of Chemistry, and Collaborative Innovation Center
of Chemical Science and Engineering, Nankai University, Tianjin 300071, China
| | - Xin Li
- State Key
Laboratory of Elemento-Organic
Chemistry, Department of Chemistry, and Collaborative Innovation Center
of Chemical Science and Engineering, Nankai University, Tianjin 300071, China
| | - Jin-Pei Cheng
- State Key
Laboratory of Elemento-Organic
Chemistry, Department of Chemistry, and Collaborative Innovation Center
of Chemical Science and Engineering, Nankai University, Tianjin 300071, China
| |
Collapse
|
31
|
Xue XS, Yang C, Li X, Cheng JP. Computational Study on the pKa Shifts in Proline Induced by Hydrogen-Bond-Donating Cocatalysts. J Org Chem 2014; 79:1166-73. [DOI: 10.1021/jo402605n] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Xiao-Song Xue
- State Key Laboratory of Elemento-Organic
Chemistry, Collaborative Innovation Center of Chemical Science and
Engineering, Department of Chemistry, Nankai University, Tianjin 300071, China
| | - Chen Yang
- State Key Laboratory of Elemento-Organic
Chemistry, Collaborative Innovation Center of Chemical Science and
Engineering, Department of Chemistry, Nankai University, Tianjin 300071, China
| | - Xin Li
- State Key Laboratory of Elemento-Organic
Chemistry, Collaborative Innovation Center of Chemical Science and
Engineering, Department of Chemistry, Nankai University, Tianjin 300071, China
| | - Jin-Pei Cheng
- State Key Laboratory of Elemento-Organic
Chemistry, Collaborative Innovation Center of Chemical Science and
Engineering, Department of Chemistry, Nankai University, Tianjin 300071, China
| |
Collapse
|
32
|
Abstract
Computational prediction of condensed phase acidity is a topic of much interest in the field today. We introduce the methods available for predicting gas phase acidity and pKas in aqueous and non-aqueous solvents including high-level electronic structure methods, empirical linear free energy relationships (LFERs), implicit solvent methods, explicit solvent statistical free energy methods, and hybrid implicit–explicit approaches. The focus of this paper is on implicit solvent methods, and we review recent developments including new electronic structure methods, cluster-continuum schemes for calculating ionic solvation free energies, as well as address issues relating to the choice of proton solvation free energy to use with implicit solvation models, and whether thermodynamic cycles are necessary for the computation of pKas. A comparison of the scope and accuracy of implicit solvent methods with ab initio molecular dynamics free energy methods is also presented. The present status of the theory and future directions are outlined.
Collapse
|
33
|
Xue XS, Yang C, Li X, Cheng JP. Computation of standard equilibrium acidity of C–H acids in ionic media: shedding light on predicting changes of chemical behavior by switching solvent system from molecular to ionic. Org Chem Front 2014. [DOI: 10.1039/c3qo00070b] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|