A Computational Protocol Combining DFT and Cheminformatics for Prediction of pH-Dependent Redox Potentials.
Ontology highlight
ABSTRACT:
SUBMITTER: Fornari RP
PROVIDER: S-EPMC8271517 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:This study investigates four machine-learning (ML) models to predict the redox potentials of phenazine derivatives in dimethoxyethane using density functional theory (DFT). A small data set of 151 phenazine derivatives having only one type of functional group per molecule (20 unique groups) was used for the training. Prediction accuracy was improved by a combined strategy of feature selection and hyperparameter optimization, using the external validation set. Models were evaluated on the external test set containing new functional groups and diverse molecular structures. High prediction accuracies of R 2 > 0.74 were obtained on the external test set. Despite being trained on the molecules with a single type of functional group, models were able to predict the redox potentials of derivatives containing multiple and different types of functional groups with good accuracies (R 2 > 0.7). This type of performance for predicting redox potential from such a small and simple data set of phenazine derivatives has never been reported before. Redox flow batteries (RFBs) are emerging as promising candidates for energy storage systems. However, new green and efficient materials are required for their widespread usage. We believe that the hybrid DFT-ML approach demonstrated in this report would help in accelerating the virtual screening of phenazine derivatives, thus saving computational and experimental costs. Using this approach, we have identified promising phenazine derivatives for green energy storage systems such as RFBs.
| S-EPMC9017108 | biostudies-literature
Project description:The pH dependence of the reduction potential E° for a metalloprotein indicates that the protonation state of at least one residue near the redox site changes and may be important for its activity. The responsible residue is usually identified by site-specific mutagenesis, which may be time-consuming. Here, the titration of E° for Chromatium vinosum high-potential iron-sulfur protein is predicted to be in good agreement with experiment using density functional theory and Poisson-Boltzmann calculations if only the sole histidine undergoes changes in protonation. The implementation of this approach into CHARMMing, a user-friendly web-based portal, allows users to identify residues in other proteins causing similar pH dependence.
| S-EPMC3691860 | biostudies-literature
Project description:The data presented in this paper are related to the research article titled "Redox Behaviour of [Ru(?-diketonato)3] Compounds" [1]. This paper presents structural and energy data obtained from the density functional theory (DFT) computations. The energy data is related to experimentally obtained redox potential values. Various relationships are presented for the RuIII/II and RuIII/IV redox couples, involving both their experimental redox data as well as DFT calculated data, such as frontier orbital energies (E HOMO and E LUMO) and calculated Mulliken electronegativity values.
| S-EPMC7206207 | biostudies-literature
Project description:BackgroundOpen Source cheminformatics toolkits such as OpenBabel, the CDK and the RDKit share the same core functionality but support different sets of file formats and forcefields, and calculate different fingerprints and descriptors. Despite their complementary features, using these toolkits in the same program is difficult as they are implemented in different languages (C++ versus Java), have different underlying chemical models and have different application programming interfaces (APIs).ResultsWe describe Cinfony, a Python module that presents a common interface to all three of these toolkits, allowing the user to easily combine methods and results from any of the toolkits. In general, the run time of the Cinfony modules is almost as fast as accessing the underlying toolkits directly from C++ or Java, but Cinfony makes it much easier to carry out common tasks in cheminformatics such as reading file formats and calculating descriptors.ConclusionBy providing a simplified interface and improving interoperability, Cinfony makes it easy to combine complementary features of OpenBabel, the CDK and the RDKit.
| S-EPMC2646723 | biostudies-literature
Project description:This paper describes salient features of the C++ programming language and its programming ecosystem, with emphasis on how the language affects scientific software development. Brief history of C++ and its predecessor the C language is provided. Most important aspects of the language that define models of programming are described in greater detail and illustrated with code examples. Special attention is paid to the interoperability between C++ and other high-level languages commonly used in cheminformatics, machine learning, data processing and statistical computing.
| S-EPMC7006409 | biostudies-literature
Project description:The ability to rationally modify targeted physical and biological features of a protein of interest holds promise in numerous academic and industrial applications and paves the way towards de novo protein design. In particular, bioprocesses that utilize the remarkable properties of enzymes would often benefit from mutants that remain active at temperatures that are either higher or lower than the physiological temperature, while maintaining the biological activity. Many in silico methods have been developed in recent years for predicting the thermodynamic stability of mutant proteins, but very few have focused on thermostability. To bridge this gap, we developed an algorithm for predicting the best descriptor of thermostability, namely the melting temperature Tm, from the protein's sequence and structure. Our method is applicable when the Tm of proteins homologous to the target protein are known. It is based on the design of several temperature-dependent statistical potentials, derived from datasets consisting of either mesostable or thermostable proteins. Linear combinations of these potentials have been shown to yield an estimation of the protein folding free energies at low and high temperatures, and the difference of these energies, a prediction of the melting temperature. This particular construction, that distinguishes between the interactions that contribute more than others to the stability at high temperatures and those that are more stabilizing at low T, gives better performances compared to the standard approach based on T-independent potentials which predict the thermal resistance from the thermodynamic stability. Our method has been tested on 45 proteins of known Tm that belong to 11 homologous families. The standard deviation between experimental and predicted Tm's is equal to 13.6°C in cross validation, and decreases to 8.3°C if the 6 worst predicted proteins are excluded. Possible extensions of our approach are discussed.
| S-EPMC3960129 | biostudies-literature
Project description:The unraveling and control of protein stability at different temperatures is a fundamental problem in biophysics that is substantially far from being quantitatively and accurately solved, as it requires a precise knowledge of the temperature dependence of amino acid interactions. In this paper we attempt to gain insight into the thermal stability of proteins by designing a tool to predict the full stability curve as a function of the temperature for a set of 45 proteins belonging to 11 homologous families, given their sequence and structure, as well as the melting temperature (Tm) and the change in heat capacity (ΔCP) of proteins belonging to the same family. Stability curves constitute a fundamental instrument to analyze in detail the thermal stability and its relation to the thermodynamic stability, and to estimate the enthalpic and entropic contributions to the folding free energy. In summary, our approach for predicting the protein stability curves relies on temperature-dependent statistical potentials derived from three datasets of protein structures with targeted thermal stability properties. Using these potentials, the folding free energies (ΔG) at three different temperatures were computed for each protein. The Gibbs-Helmholtz equation was then used to predict the protein's stability curve as the curve that best fits these three points. The results are quite encouraging: the standard deviations between the experimental and predicted Tm's, ΔCP's and folding free energies at room temperature (ΔG25) are equal to 13° C, 1.3 kcal/(mol° C) and 4.1 kcal/mol, respectively, in cross-validation. The main sources of error and some further improvements and perspectives are briefly discussed.
| S-EPMC4102405 | biostudies-literature
Project description:Hydrogen bonding (H-bonding) is generally thought to play an important role in tuning the electronic structure and reactivity of metal-sulfur sites in proteins. To develop a quantitative understanding of this effect, S K-edge X-ray absorption spectroscopy (XAS) has been employed to directly probe ligand-metal bond covalency, where it has been found that protein active sites are significantly less covalent than their related model complexes. Sulfur K-edge XAS data are reported here on a series of P450 model complexes with increasing H-bonding to the ligated thiolate from its substituent. The XAS spectroscopic results show a dramatic decrease in preedge intensity. DFT calculations reproduce these effects and show that the observed changes are in fact solely due to H-bonding and not from the inductive effect of the substituent on the thiolate. These calculations also indicate that the H-bonding interaction in these systems is mainly dipolar in nature. The -2.5 kcal/mol energy of the H-bonding interaction was small relative to the large change in ligand-metal bond covalency (30%) observed in the data. A bond decomposition analysis of the total energy is developed to correlate the preedge intensity change to the change in Fe-S bonding interaction on H-bonding. This effect is greater for the reduced than the oxidized state, leading to a 260 mV increase in the redox potential. A simple model shows that E degrees should vary approximately linearly with the covalency of the Fe-S bond in the oxidized state, which can be determined directly from S K-edge XAS.
| S-EPMC2880190 | biostudies-literature
Project description:Ions are involved in multiple biological processes and may exist bound to biomolecules or may be associated with their surface. Although the presence of ions in nucleic acids has traditionally gained more interest, ion-protein interactions, often with a marked dependency on pH, are beginning to gather attention. Here we present a detailed analysis on the binding and distribution of ions around β-lactoglobulin using a constant-pH MD (CpHMD) method, at a pH range 3-8, and compare it with the more traditional Poisson-Boltzmann (PB) model and the existing experimental data. Most analyses used ion concentration maps built around the protein, obtained from either the CpHMD simulations or PB calculations. The requirements of approximate charge neutrality and ionic strength equal to bulk, imposed on the MD box, imply that the absolute value of the ion excess should be half the protein charge, which is in agreement with experimental observation on other proteins ( Proc. Natl. Acad. Sci. U.S.A. 2021, 118, e2015879118) and lends support to this protocol. In addition, the protein total charge (including territorially bound ions) estimated with MD is in excellent agreement with electrophoretic measurements. Overall, the CpHMD simulations show good agreement with the nonlinear form of the PB (NLPB) model but not with its linear form, which involves a theoretical inconsistency in the calculation of the concentration maps. In several analyses, the observed pH-dependent trends for the counterions and co-ions are those generally expected, and the ion concentration maps correctly converge to the bulk ionic strength as one moves away from the protein. Despite the overall similarity, the CpHMD and NLPB approaches show some discrepancies when analyzed in more detail, which may be related to an apparent overestimation of counterion excess and underestimation of co-ion exclusion by the NLPB model, particularly at short distances from the protein.
| S-EPMC9776516 | biostudies-literature
Project description:We present a computational scheme to compute the pH-dependence of binding free energy with explicit solvent. Despite the importance of pH, the effect of pH has been generally neglected in binding free energy calculations because of a lack of accurate methods to model it. To address this limitation, we use a constant-pH methodology to obtain a true ensemble of multiple protonation states of a titratable system at a given pH and analyze the ensemble using the Bennett acceptance ratio (BAR) method. The constant pH method is based on the combination of enveloping distribution sampling (EDS) with the Hamiltonian replica exchange method (HREM), which yields an accurate semi-grand canonical ensemble of a titratable system. By considering the free energy change of constraining multiple protonation states to a single state or releasing a single protonation state to multiple states, the pH dependent binding free energy profile can be obtained. We perform benchmark simulations of a host-guest system: cucurbit[7]uril (CB[7]) and benzimidazole (BZ). BZ experiences a large pKa shift upon complex formation. The pH-dependent binding free energy profiles of the benchmark system are obtained with three different long-range interaction calculation schemes: a cutoff, the particle mesh Ewald (PME), and the isotropic periodic sum (IPS) method. Our scheme captures the pH-dependent behavior of binding free energy successfully. Absolute binding free energy values obtained with the PME and IPS methods are consistent, while cutoff method results are off by 2 kcal mol(-1) . We also discuss the characteristics of three long-range interaction calculation methods for constant-pH simulations.
| S-EPMC4815317 | biostudies-literature