Mechanisms and free energies of enzymatic reactions.
Ontology highlight
ABSTRACT:
SUBMITTER: Gao J
PROVIDER: S-EPMC4477011 | biostudies-literature | 2006 Aug
REPOSITORIES: biostudies-literature
Publications
Mechanisms and free energies of enzymatic reactions.
Gao Jiali J Ma Shuhua S Major Dan T DT Nam Kwangho K Pu Jingzhi J Truhlar Donald G DG
Chemical reviews 20060801 8
Similar Datasets
Project description:The temperature dependence of the kinetics and the mechanisms of ozone reactions with 19 olefins and 3 alkynes were investigated. The second-order rate constants (kO3) for ozone reactions with olefins were mostly in the range of 103-106 M-1 s-1, with activation energies of 17.4-37.7 kJ mol-1. In comparison, alkynes had lower kO3 (∼102 M-1 s-1) and higher activation energies (36.7-48.1 kJ mol-1). Reactivities of both olefins and alkynes are mainly influenced by inductive effects of substituents, with steric effects observed for cyclic olefins. 2-Buten-1,4-dial (BDA), synthesized with a novel method, is a toxic olefinic oxidation product from phenols. Its cis- and trans-isomers show distinct reactivities with ozone, with kO3 (20 °C) of 3.0 × 103 and 1.2 × 104 M-1 s-1, respectively. Two mols of glyoxal were formed per mol of ozonated BDA, with a slow release of the second mol from an α-hydroxyalkylhydroperoxide intermediate. 2-Ethynylbenzaldehyde reacts with ozone with a stoichiometry of 1:1 and kO3 (20 °C) = 1.6 × 102 M-1 s-1. Ozone attacks the ethynyl group, yielding a carboxyl product (2-carboxybenzaldehyde, 54%), an aldehyde product (phthaldialdehyde), and a dicarbonyl product with a stoichiometric release of H2O2 (21%). This study provides kinetic and mechanistic information for assessing the abatement of olefin- and alkyne-containing micropollutants by ozonation at various temperatures.
| S-EPMC11912337 | biostudies-literature
Project description:Transition path sampling (TPS) is widely used for the calculations of reaction rates, transition state structures, and reaction coordinates of condensed phase systems. Here we discuss a scheme for the calculation of free energies using the ensemble of TPS reactive trajectories in combination with a window-based sampling technique for enzyme-catalyzed reactions. We calculate the free energy profiles of the reactions catalyzed by the human methionine S-adenosyltransferase (MAT2A) enzyme and the Plasmodium vivax adenosine deaminase (pvADA) enzyme to assess the accuracy of this method. MAT2A catalyzes the formation of S-adenosine-l-methionine following a SN2 mechanism, and using our method, we estimate the free energy barrier for this reaction to be 16 kcal mol-1, which is in excellent agreement with the experimentally measured activation energy of 17.27 kcal mol-1. Furthermore, for the pvADA enzyme-catalyzed reaction we estimate a free energy barrier of 21 kcal mol-1, and the calculated free energy profile is similar to that predicted from experimental observations. Calculating free energies by employing our simple method within TPS provides significant advantages over methods such as umbrella sampling because it is free from any applied external bias, is accurate compared to experimental measurements, and has a reasonable computational cost.
| S-EPMC9444332 | biostudies-literature
Project description:Molecular dynamics is a powerful tool for studying the thermodynamics and kinetics of complex molecular events. However, these simulations can rarely sample the required time scales in practice. Transition path sampling overcomes this limitation by collecting unbiased trajectories and capturing the relevant events. Moreover, the integration of machine learning can boost the sampling while simultaneously learning a quantitative representation of the mechanism. Still, the resulting trajectories are by construction non-Boltzmann-distributed, preventing the calculation of free energies and rates. We developed an algorithm to approximate the equilibrium path ensemble from machine-learning-guided path sampling data. At the same time, our algorithm provides efficient sampling, mechanism, free energy, and rates of rare molecular events at a very moderate computational cost. We tested the method on the folding of the mini-protein chignolin. Our algorithm is straightforward and data-efficient, opening the door to applications in many challenging molecular systems.
| S-EPMC10753783 | biostudies-literature
Project description:In enzyme mechanistic studies and mutant design, it is highly desirable to know the individual residue contributions to the reaction free energy and barrier. In this work, we show that such free energy contributions from each residue can be readily obtained by postprocessing ab initio quantum mechanical molecular mechanical (ai-QM/MM) free energy simulation trajectories. Specifically, through a mean force integration along the minimum free energy pathway, one can obtain the electrostatic, polarization, and van der Waals contributions from each residue to the free energy barrier. Separately, a similar analysis procedure allows us to assess the contribution from different collective variables along the reaction coordinate. The chorismate mutase reaction is used to demonstrate the utilization of these two trajectory analysis tools.
| S-EPMC10835707 | biostudies-literature
Project description:Linear free energy relationships (LFERs) for substituent effects on reactions that proceed through similar transition states provide insight into transition state structures. A classical approach to the analysis of LFERs showed that differences in the slopes of Brønsted correlations for addition of substituted alkyl alcohols to ring-substituted 1-phenylethyl carbocations and to the β-galactopyranosyl carbocation intermediate of reactions catalyzed by β-galactosidase provide evidence that the enzyme catalyst modifies the curvature of the energy surface at the saddle point for the transition state for nucleophile addition. We have worked to generalize the use of LFERs in the determination of enzyme mechanisms. The defining property of enzyme catalysts is their specificity for binding the transition state with a much higher affinity than the substrate. Triosephosphate isomerase (TIM), orotidine 5'-monophosphate decarboxylase (OMPDC), and glycerol 3-phosphate dehydrogenase (GPDH) show effective catalysis of reactions of phosphorylated substrates and strong phosphite dianion activation of reactions of phosphodianion truncated substrates, with rate constants kcat/Km (M-1 s-1) and kcat/KdKHPi (M-2 s-1), respectively. Good linear logarithmic correlations, with a slope of 1.1, between these kinetic parameters determined for reactions catalyzed by five or more variant forms of each catalyst are observed, where the protein substitutions are mainly at side chains which function to stabilize the cage complex between the enzyme and substrate. This shows that the enzyme-catalyzed reactions of a whole substrate and substrate pieces proceed through transition states of similar structures. It provides support for the proposal that the dianion binding energy of whole phosphodianion substrates and of phosphite dianion is used to drive the conversion of these protein catalysts from flexible and entropically rich ground states to stiff and catalytically active Michaelis complexes that show the same activity toward catalysis of the reactions of whole and phosphodianion truncated substrates. There is a good linear correlation, with a slope of 0.73, between values of the dissociation constants log Ki for release of the transition state analog phosphoglycolate (PGA) trianion and log kcat/Km for isomerization of GAP for wild-type and variants of TIM. This correlation shows that the substituted amino acid side chains act to stabilize the complex between TIM and the PGA trianion and that ca. 70% of this stabilization is observed at the transition state for substrate deprotonation. The correlation provides evidence that these side chains function to enhance the basicity of the E165 side chain of TIM, which deprotonates the bound carbon acid substrate. There is a good linear correlation, with a slope of 0.74, between the values of ΔG‡ and ΔG° determined by electron valence bond (EVB) calculations to model deprotonation of dihydroxyacetone phosphate (DHAP) in water and when bound to wild-type and variant forms of TIM to form the enediolate reaction intermediate. This correlation provides evidence that the stabilizing interactions of the transition state for TIM-catalyzed deprotonation of DHAP are optimized by placement of amino acid side chains in positions that provide for the maximum stabilization of the charged reaction intermediate, relative to the neutral substrate.
| S-EPMC8157535 | biostudies-literature
Project description:SKBR3 breast cancer cell extracts were digested with trypsin (or LysC) on short time scales (7min, 15min, 30min, 1h, and 18h), and the results were compared to overnight digestion protocols.
2018-09-24 | PXD010306 | Pride
Project description:Hydride transfer reactions are very important chemical reactions in organic chemistry. It has been a chemist's dream to predict the rate constants of hydride transfer reactions by only using the physical parameters of the reactants. To realize this dream, we have developed a kinetic equation (Zhu equation) in our previous papers to predict the activation free energies of various chemical reactions using the activation free energies of the corresponding self-exchange reactions and the related bond dissociation energies or redox potentials of the reactants. Because the activation free energy of the hydride self-exchange reaction is difficult to measure using the experimental method, in this study, the activation free energies of 159 hydride self-exchange reactions in acetonitrile at 298 K were systematically computed using an accurately benchmarked density functional theory method with a precision of 1.1 kcal mol-1. The results show that the range of the activation free energies of the 159 hydride self-exchange reactions is from 16.1 to 46.6 kcal mol-1. The activation free energies of 25 122 hydride transfer reactions in acetonitrile at 298 K can be estimated using the activation free energies of the 159 hydride self-exchange reactions and the corresponding heterolytic bond dissociation free energies of the reactants. The effects of the heteroatom, substituent, and aromaticity on the activation free energies of hydride self-exchange reactions were examined. The results show that heteroatoms, substituents at the reaction center, and the aromaticity of reactants, all have remarkable effects on the activation free energy of hydride self-exchange reactions. All kinetic information provided in this work on the hydride self-exchange reactions in acetonitrile at 298 K should be very useful in chemical labs and chemical industry.
| S-EPMC6641257 | biostudies-literature
Project description:The calculation of absolute binding free energies (ABFEs) for protein-ligand systems has long been a challenge. Recently, refined force fields and algorithms have improved the quality of the ABFE calculations. However, achieving the level of accuracy required to inform drug discovery efforts remains difficult. Here, we present a transferable enhanced sampling strategy to accurately calculate absolute binding free energies using OneOPES with simple geometric collective variables. We tested the strategy on two protein targets, BRD4 and Hsp90, complexed with a total of 17 chemically diverse ligands, including both molecular fragments and drug-like molecules. Our results show that OneOPES accurately predicts protein-ligand binding affinities with a mean unsigned error within 1 kcal mol-1 of experimentally determined free energies, without the need to tailor the collective variables to each system. Furthermore, our strategy effectively samples different ligand binding modes and consistently matches the experimentally determined structures regardless of the initial protein-ligand configuration. Our results suggest that the proposed OneOPES strategy can be used to inform lead optimization campaigns in drug discovery and to study protein-ligand binding and unbinding mechanisms.
| S-EPMC11457222 | biostudies-literature
Project description:Understanding molecular mechanisms of enzymatic reactions is of vital importance in biochemistry and biophysics. Here, we introduce new functions of hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in the GENESIS program to compute the minimum-energy pathways (MEPs) and free-energy profiles of enzymatic reactions. For this purpose, an interface in GENESIS is developed to utilize a highly parallel electronic structure program, QSimulate-QM (https://qsimulate.com), calling it as a shared library from GENESIS. Second, algorithms to search the MEP are implemented, combining the string method (E et al. J. Chem. Phys. 2007, 126, 164103) with the energy minimization of the buffer MM region. The method implemented in GENESIS is applied to an enzyme, triosephosphate isomerase, which converts dihyroxyacetone phosphate to glyceraldehyde 3-phosphate in four proton-transfer processes. QM/MM-molecular dynamics simulations show performances of greater than 1 ns/day with the density functional tight binding (DFTB), and 10-30 ps/day with the hybrid density functional theory, B3LYP-D3. These performances allow us to compute not only MEP but also the potential of mean force (PMF) of the enzymatic reactions using the QM/MM calculations. The barrier height obtained as 13 kcal mol-1 with B3LYP-D3 in the QM/MM calculation is in agreement with the experimental results. The impact of conformational sampling in PMF calculations and the level of electronic structure calculations (DFTB vs B3LYP-D3) suggests reliable computational protocols for enzymatic reactions without high computational costs.
| S-EPMC10986901 | biostudies-literature
Project description:We present an energy function for predicting binding free energies of protein-protein complexes, using the three-dimensional structures of the complex and unbound proteins as input. Our function is a linear combination of nine terms and achieves a correlation coefficient of 0.63 with experimental measurements when tested on a benchmark of 144 complexes using leave-one-out cross validation. Although we systematically tested both atomic and residue-based scoring functions, the selected function is dominated by residue-based terms. Our function is stable for subsets of the benchmark stratified by experimental pH and extent of conformational change upon complex formation, with correlation coefficients ranging from 0.61 to 0.66.
| S-EPMC3375440 | biostudies-literature