Protein Conformational Landscapes and Catalysis. Influence of Active Site Conformations in the Reaction Catalyzed by L-Lactate Dehydrogenase

Constrained Nuclear-Electronic Orbital Transition State Theory Using Energy Surfaces with Nuclear Quantum Effects

Hydrogen-atom transfer is crucial in a myriad of chemical and biological processes, yet the accurate and efficient description of hydrogen-atom transfer reactions and kinetic isotope effects remains challenging due to significant quantum effects on hydrogenic motion, especially tunneling and zero-point energy. In this paper, we combine transition state theory (TST) with the recently developed constrained nuclear-electronic orbital (CNEO) theory to propose a new transition state theory denoted CNEO-TST. We use CNEO-TST with CNEO density functional theory (CNEO-DFT) to predict reaction rate constants for two prototypical gas-phase hydrogen-atom transfer reactions and their deuterated isotopologic reactions. CNEO-TST is similar to conventional TST except that it employs constrained minimized energy surfaces to include zero-point energy and shallow tunneling effects in the effective potential. We find that the new theory predicts reaction rates quite accurately at room temperature. The effective potential surface must be generated by CNEO theory rather than by ordinary electronic structure theory, but because of the favorable computational scaling of CNEO-DFT, the cost is economical even for large systems. Our results show that dynamics calculations with this approach achieve accuracy comparable to variational TST with a semiclassical multidimensional tunneling transmission coefficient at and above room temperature. Therefore, CNEO-TST can be a useful tool for rate prediction, even for reactions involving highly quantal motion, such as many chemical and biochemical reactions involving transfers of hydrogen atoms, protons, or hydride ions.

Electrostatic Preorganization in Three Distinct Heterogeneous Proteasome β-Subunits

The origin of the enzyme's powerful role in accelerating chemical reactions is one of the most critical and still widely discussed questions. It is already accepted that enzymes impose an electrostatic field onto their substrates by adopting complex three-dimensional structures; therefore, the preorganization of electric fields inside protein active sites has been proposed as a crucial contributor to catalytic mechanisms and rate constant enhancement. In this work, we focus on three catalytically active β-subunits of 20S proteasomes with low sequence identity (∼30%) whose active sites, although situated in an electrostatically miscellaneous environment, catalyze the same chemical reaction with similar catalytic efficiency. Our in silico experiments reproduce the experimentally observed equivalent reactivity of the three sites and show that obliteration of the electrostatic potential in all active sites would deprive the enzymes of their catalytic power by slowing down the chemical process by a factor of 1035. To regain enzymatic efficiency, besides catalytic Thr1 and Lys33 residues, the presence of aspartic acid in position 17 and an aqueous solvent is required, proving that the electrostatic potential generated by the remaining residues is insignificant for catalysis. Moreover, it was found that the gradual decay of atomic charges on Asp17 strongly correlates with the enzyme's catalytic rate deterioration as well as with a change in the charge distributions due to introduced mutations. The computational procedure used and described here may help identify key residues for catalysis in other biomolecular systems and consequently may contribute to the process of designing enzyme-like synthetic catalysts.

Cooperative Conformational Transitions Underpin the Activation Heat Capacity in the Temperature Dependence of Enzyme Catalysis

Many enzymes display non-Arrhenius behavior with curved Arrhenius plots in the absence of denaturation. There has been significant debate about the origin of this behavior and recently the role of the activation heat capacity (ΔCP⧧) has been widely discussed. If enzyme-catalyzed reactions occur with appreciable negative values of ΔCP⧧ (arising from narrowing of the conformational space along the reaction coordinate), then curved Arrhenius plots are a consequence. To investigate these phenomena in detail, we have collected high precision temperature-rate data over a wide temperature interval for a model glycosidase enzyme MalL, and a series of mutants that change the temperature-dependence of the enzyme-catalyzed rate. We use these data to test a range of models including macromolecular rate theory (MMRT) and an equilibrium model. In addition, we have performed extensive molecular dynamics (MD) simulations to characterize the conformational landscape traversed by MalL in the enzyme-substrate complex and an enzyme-transition state complex. We have crystallized the enzyme in a transition state-like conformation in the absence of a ligand and determined an X-ray crystal structure at very high resolution (1.10 Å). We show (using simulation) that this enzyme-transition state conformation has a more restricted conformational landscape than the wildtype enzyme. We coin the term "transition state-like conformation (TLC)" to apply to this state of the enzyme. Together, these results imply a cooperative conformational transition between an enzyme-substrate conformation (ES) and a transition-state-like conformation (TLC) that precedes the chemical step. We present a two-state model as an extension of MMRT (MMRT-2S) that describes the data along with a convenient approximation with linear temperature dependence of the activation heat capacity (MMRT-1L) that can be used where fewer data points are available. Our model rationalizes disparate behavior seen for MalL and previous results for a thermophilic alcohol dehydrogenase and is consistent with a raft of data for other enzymes. Our model can be used to characterize the conformational changes required for enzyme catalysis and provides insights into the role of cooperative conformational changes in transition state stabilization that are accompanied by changes in heat capacity for the system along the reaction coordinate. TLCs are likely to be of wide importance in understanding the temperature dependence of enzyme activity and other aspects of enzyme catalysis.

Oliveira AS, Mulholland AJ, Serapian SA, Colombo G. Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations

Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.

Structure, Mutagenesis, and QM:MM Modeling of 3-Ketosteroid Δ1-Dehydrogenase from Sterolibacterium denitrificans─The Role of a New Putative Membrane-Associated Domain and Proton-Relay System in Catalysis

3-Ketosteroid Δ1-dehydrogenases (KstD) are important microbial flavin enzymes that initiate the metabolism of steroid ring A and find application in the synthesis of steroid drugs. We present a structure of the KstD from Sterolibacterium denitrificans (AcmB), which contains a previously uncharacterized putative membrane-associated domain and extended proton-relay system. The experimental and theoretical studies show that the steroid Δ1-dehydrogenation proceeds according to the Ping-Pong bi-bi kinetics and a two-step base-assisted elimination (E2cB) mechanism. The mechanism is validated by evaluating the experimental and theoretical kinetic isotope effect for deuterium-substituted substrates. The role of the active-site residues is quantitatively assessed by point mutations, experimental activity assays, and QM/MM MD modeling of the reductive half-reaction (RHR). The pre-steady-state kinetics also reveals that the low pH (6.5) optimum of AcmB is dictated by the oxidative half-reaction (OHR), while the RHR exhibits a slight optimum at the pH usual for the KstD family of 8.5. The modeling confirms the origin of the enantioselectivity of C2-H activation and substrate specificity for Δ4-3-ketosteroids. Finally, the cholest-4-en-3-one turns out to be the best substrate of AcmB in terms of ΔG of binding and predicted rate of dehydrogenation.

Revisiting the Warburg Effect with Focus on Lactate

Rewired metabolism is acknowledged as one of the drivers of tumor growth. As a result, aerobic glycolysis, or the Warburg effect, is a feature of many cancers. Increased glucose uptake and glycolysis provide intermediates for anabolic reactions necessary for cancer cell proliferation while contributing sufficient energy. However, the accompanying increased lactate production, seemingly wasting glucose carbon, was originally explained only by the need to regenerate NAD⁺ for successive rounds of glycolysis by the lactate dehydrogenase (LDH) reaction in the cytosol. After the discovery of a mitochondrial LDH isoform, lactate oxidation entered the picture, and lactate was recognized as an important oxidative fuel. It has also been revealed that lactate serves a variety of signaling functions and helps cells adapt to the new environment. Here, we discuss recent findings on lactate metabolism and signaling in cancer while attempting to explain why the Warburg effect is adopted by cancer cells.

Natural products targeting human lactate dehydrogenases for cancer therapy: A mini review

Reprogramming cancer metabolism has become the hallmark of cancer progression. As the key enzyme catalyzing the conversion of pyruvate to lactate in aerobic glycolysis of cancer cells, human lactate dehydrogenase (LDH) has been a promising target in the discovery of anticancer agents. Natural products are important sources of new drugs. Up to now, some natural compounds have been reported with the activity to target LDH. To give more information on the development of LDH inhibitors and application of natural products, herein, we reviewed the natural compounds with inhibition of LDH from diverse structures and discussed the future direction of the discovery of natural LDH inhibitors for cancer therapy.

Covalent Inhibition of the Human 20S Proteasome with Homobelactosin C Inquired by QM/MM Studies

20S proteasome is a main player in the protein degradation pathway in the cytosol, thus intervening in multiple pivotal cellular processes. Over the years the proteasome has emerged as a crucial target for the treatment of many diseases such as neurodegenerative diseases, cancer, autoimmune diseases, developmental disorders, cystic fibrosis, diabetes, cardiac diseases, atherosclerosis, and aging. In this work, the mechanism of proteasome covalent inhibition with bisbenzyl-protected homobelactosin C (hBelC) was explored using quantum mechanics/molecular mechanics (QM/MM) methods. Molecular dynamic simulations were used to describe key interactions established between the hBelC and its unique binding mode in the primed site of the β5 subunit. The free energy surfaces were computed to characterize the kinetics and thermodynamics of the inhibition process. This study revealed that although the final inhibition product for hBelC is formed according to the same molecular mechanism as one described for hSalA, the free energy profile of the reaction pathway differs significantly from the one previously reported for γ-lactam-β-lactone containing inhibitors in terms of the height of the activation barrier as well as the stabilization of the final product. Moreover, it was proved that high stabilization of the covalent adduct formed between β5-subunit and hBelC, together with the presence of aminocarbonyl side chain in the structure of the inhibitor which prevents the hydrolysis of the ester bond from taking place, determines its irreversible character.

Protein dynamics implications of the low- and high-temperature denaturation of myoglobin

We reinvestigate a simple model used in the literature concerning the thermodynamic analysis of protein cold denaturation. We derive an exact thermodynamic expression for cold denaturation and give a better approximation than exists in the literature for predicting cold denaturation temperatures in the two-state model. We discuss the "dark-side" implications of this work for previous temperature-dependent protein dynamics experiments and discuss microfluidic experimental technologies, which could explore the thermal stability range of proteins below the bulk freezing point of water.

High Pretreatment LDH Predicts Poor Prognosis in Hypopharyngeal Cancer

Background

Elevated pretreatment lactate dehydrogenase (LDH) has been associated with poor prognosis in various malignancies; however, its prognostic role in hypopharyngeal cancer remains elusive. In this study, we aimed to assess the association between pretreatment LDH and clinical outcome of hypopharyngeal cancer.

Methods

We retrospectively collected 198 hypopharyngeal cancer patients treated with surgery in our institution between 2004 and 2018. The prognostic role of pretreatment LDH was explored by using univariate and multivariate analyses. Besides, subgroup analysis was performed based on T stage.

Results

Three-year and Five-year of disease-free survival (DFS, 67.0 vs. 57.4%, 65.8 vs. 39.8%, p = 0.007) and overall survival (OS, 74.8 vs. 68.9%, 66.8 vs. 50.8%, p = 0.006) exhibited significant differences between low LDH level and high LDH level groups. Univariate analysis showed that pretreatment elevated serum LDH served as an unfavorable determinant with regard to DFS and OS. Further multivariate analysis also confirmed that LDH was an independent predictor for DFS and OS. Additionally, N status and age were also found to be significantly associated with both DFS and OS.

Conclusion

Pretreatment elevated serum LDH is an inferior prognostic factor for patients with hypopharyngeal cancer. These results should be validated by more multicenter and prospective studies.

Transfer hydrogenations catalyzed by streptavidin-hosted secondary amine organocatalysts

Here, the streptavidin-biotin technology was applied to enable organocatalytic transfer hydrogenation. By introducing a biotin-tethered pyrrolidine (1) to the tetrameric streptavidin (T-Sav), the resulting hybrid catalyst was able to mediate hydride transfer from dihydro-benzylnicotinamide (BNAH) to α,β-unsaturated aldehydes. Hydrogenation of cinnamaldehyde and some of its aryl-substituted analogues was found to be nearly quantitative. Kinetic measurements revealed that the T-Sav:1 assembly possesses enzyme-like behavior, whereas isotope effect analysis, performed by QM/MM simulations, illustrated that the step of hydride transfer is at least partially rate-limiting. These results have proven the concept that T-Sav can be used to host secondary amine-catalyzed transfer hydrogenations.

Pressure tolerance of deep-sea enzymes can be evolved through increasing volume changes in protein transitions: a study with lactate dehydrogenases from abyssal and hadal fishes

We explore the principles of pressure tolerance in enzymes of deep-sea fishes using lactate dehydrogenases (LDH) as a case study. We compared the effects of pressure on the activities of LDH from hadal snailfishes Notoliparis kermadecensis and Pseudoliparis swirei with those from a shallow-adapted Liparis florae and an abyssal grenadier Coryphaenoides armatus. We then quantified the LDH content in muscle homogenates using mass-spectrometric determination of the LDH-specific conserved peptide LNLVQR. Existing theory suggests that adaptation to high pressure requires a decrease in volume changes in enzymatic catalysis. Accordingly, evolved pressure tolerance must be accompanied with an important reduction in the volume change associated with pressure-promoted alteration of enzymatic activity ( Δ V PP ∘ ). Our results suggest an important revision to this paradigm. Here, we describe an opposite effect of pressure adaptation-a substantial increase in the absolute value of Δ V PP ∘ in deep-living species compared to shallow-water counterparts. With this change, the enzyme activities in abyssal and hadal species do not substantially decrease their activity with pressure increasing up to 1-2 kbar, well beyond full-ocean depth pressures. In contrast, the activity of the enzyme from the tidepool snailfish, L. florae, decreases nearly linearly from 1 to 2500 bar. The increased tolerance of LDH activity to pressure comes at the expense of decreased catalytic efficiency, which is compensated with increased enzyme contents in high-pressure-adapted species. The newly discovered strategy is presumably used when the enzyme mechanism involves the formation of potentially unstable excited transient states associated with substantial changes in enzyme-solvent interactions.

Mechanism of inhibition of SARS-CoV-2 Mpro by peptidyl Michael acceptor explained by QM/MM simulations and design of new derivatives with tunable chemical reactivity

The SARS-CoV-2 main protease (Mpro) is essential for replication of the virus responsible for the COVID-19 pandemic, and one of the main targets for drug design. Here, we simulate the inhibition process of SARS-CoV-2 Mpro with a known Michael acceptor (peptidyl) inhibitor, N3. The free energy landscape for the mechanism of the formation of the covalent enzyme-inhibitor product is computed with QM/MM molecular dynamics methods. The simulations show a two-step mechanism, and give structures and calculated barriers in good agreement with experiment. Using these results and information from our previous investigation on the proteolysis reaction of SARS-CoV-2 Mpro, we design two new, synthetically accessible N3-analogues as potential inhibitors, in which the recognition and warhead motifs are modified. QM/MM modelling of the mechanism of inhibition of Mpro by these novel compounds indicates that both may be promising candidates as drug leads against COVID-19, one as an irreversible inhibitor and one as a potential reversible inhibitor.

Kramers’ Theory and the Dependence of Enzyme Dynamics on Trehalose-Mediated Viscosity

The disaccharide trehalose is accumulated in the cytoplasm of some organisms in response to harsh environmental conditions. Trehalose biosynthesis and accumulation are important for the survival of such organisms by protecting the structure and function of proteins and membranes. Trehalose affects the dynamics of proteins and water molecules in the bulk and the protein hydration shell. Enzyme catalysis and other processes dependent on protein dynamics are affected by the viscosity generated by trehalose, as described by the Kramers’ theory of rate reactions. Enzyme/protein stabilization by trehalose against thermal inactivation/unfolding is also explained by the viscosity mediated hindering of the thermally generated structural dynamics, as described by Kramers’ theory. The analysis of the relationship of viscosity–protein dynamics, and its effects on enzyme/protein function and other processes (thermal inactivation and unfolding/folding), is the focus of the present work regarding the disaccharide trehalose as the viscosity generating solute. Finally, trehalose is widely used (alone or in combination with other compounds) in the stabilization of enzymes in the laboratory and in biotechnological applications; hence, considering the effect of viscosity on catalysis and stability of enzymes may help to improve the results of trehalose in its diverse uses/applications.

Structural and functional regulation of lactate dehydrogenase-A in cancer

Dysregulated metabolism is one of the hallmarks of cancer. Under normal physiological conditions, ATP is primarily generated by oxidative phosphorylation. Cancers commonly undergo a dramatic shift toward glycolysis, despite the presence of oxygen. This phenomenon is known as the Warburg effect, and requires the activity of LDHA. LDHA converts pyruvate to lactate in the final step of glycolysis and is often upregulated in cancer. LDHA inhibitors present a promising therapeutic option, as LDHA blockade leads to apoptosis in cancer cells. Despite this, existing LDHA inhibitors have shown limited clinical efficacy. Here, we review recent progress in LDHA structure, function and regulation as well as strategies to target this critical enzyme.

A protocol to obtain multidimensional quantum tunneling corrections derived from QM(DFT)/MM calculations for an enzyme reaction

The multidimensional small-curvature tunneling (SCT) method with Electrostatic Embedding calculations is a compromise between an accessible computational cost and the attainment of an accurate enough estimation of tunneling for an enzyme reaction.

Loss of Hyperconjugative Effects Drives Hydride Transfer during Dihydrofolate Reductase Catalysis

Hydride transfer is widespread in nature and has an essential role in applied research. However, the mechanisms of how this transformation occurs in living organisms remain a matter of vigorous debate. Here, we examined dihydrofolate reductase (DHFR), an enzyme that catalyzes hydride from C4′ of NADPH to C6 of 7,8-dihydrofolate (H₂F). Despite many investigations of the mechanism of this reaction, the contribution of polarization of the π-bond of H₂F in driving hydride transfer remains unclear. H₂F was stereospecifically labeled with deuterium β to the reacting center, and β-deuterium kinetic isotope effects were measured. Our experimental results combined with analysis derived from QM/MM simulations reveal that hydride transfer is triggered by polarization at the C6 of H₂F. The σ C_β–H bonds contribute to the buildup of the cationic character during the chemical transformation, and hyperconjugation influences the formation of the transition state. Our findings provide key insights into the hydride transfer mechanism of the DHFR-catalyzed reaction, which is a target for antiproliferative drugs and a paradigmatic model in mechanistic enzymology.

Water Plays a Cocatalytic Role in Epoxide Ring Opening Reaction in Aspartate Proteases: A QM/MM Study

Aspartate proteases are potential targets for various diseases, and many of their inhibitors are FDA-approved drugs. However, these peptidomimetic and reversibly bound drugs become ineffective upon prolonged use. Attempts have been made to design and synthesize various nonpeptidic epoxide-based irreversible inhibitors to combat the drug-resistance enigma. Here, we study the mechanism of epoxide ring opening in two widely studied aspartate proteases, HIV-1 protease and pepsin. Our results from QM/MM molecular dynamics show that the epoxide ring opening in aspartate proteases follow a two-step mechanism with the formation of an oxyanion intermediate, stabilized by a set of water molecules in the protein active site. These water molecules by virtue of "low-barrier hydrogen bonds" with the epoxide ring reduce the intrinsic reaction barrier while remaining structurally unperturbed, thus playing a cocatalytic role. We validated our results by reproducing the experimentally observed protease/pepsin-epoxide covalent complexes as end products. The observed stability of our oxyanion intermediate in a four-water-coordinated state is also consistent with the reported stable state of the hydroxide ion in water as OH^-(H₂O)₄. Our study could pave the way for the design of new class "HIV protease irreversible inhibitors" from the acquired knowledge of the structures of intermediate and transition states traced during the explored reaction mechanism.

Osmolytes modify protein dynamics and function of tetrameric lactate dehydrogenase upon pressurization

We present a study of the combined effects of natural cosolvents (TMAO, glycine, urea) and pressure on the activity of the tetrameric enzyme lactate dehydrogenase (LDH). To this end, high-pressure stopped-flow methodology in concert with fast UV/Vis spectroscopic detection of product formation was applied. To reveal possible pressure effects on the stability and dynamics of the enzyme, FTIR spectroscopic and neutron scattering measurements were carried out. In neat buffer solution, the catalytic turnover number of the enzyme, k_cat, increases up to 1000 bar, the pressure range where dissociation of the tetrameric species to dimers sets in. Accordingly, we obtain a negative activation volume, ΔV^# = -45.3 mL mol^-1. Further, the enzyme substrate complex has a larger volume compared to the enzyme and substrate in the unbound state. The neutron scattering data show that changes in the fast internal dynamics of the enzyme are not responsible for the increase of k_cat upon compression. Whereas the magnitude of k_cat is similar in the presence of the osmolytes, the pressure of deactivation is modulated by the addition of cosolvents. TMAO and glycine increase the pressure of deactivation, and in accordance with the observed stabilizing effect both cosolvents exhibit against denaturation and/or dissociation of proteins. While urea does not markedly affect the magnitude of the Michaelis constant, K_M, both 1 M TMAO and 1 M glycine exhibit smaller K_M values of about 0.07 mM and 0.05 mM below about 1 kbar. Such positive effect on the substrate affinity could be rationalized by the effect the two cosolutes impose on the thermodynamic activities of the reactants, which reflect changes in water-mediated intermolecular interactions. Our data show that the intracellular milieu, i.e., the solution conditions that have evolved, may be sufficient to maintain enzymatic activity under extreme environmental conditions, including the whole pressure range encountered on Earth.

Convergence of theory and experiment on the role of preorganization, quantum tunneling and enzyme motions into flavoenzyme-catalyzed hydride transfer

Hydride transfer is one of the most common reactions catalyzed by enzymatic systems and it has become an object of study due to possible significant quantum tunneling effects. In the present work, we provide a combination of theoretical QM/MM simulations and experimental measurements of the rate constants and kinetic isotopic effects (KIEs) for the hydride transfer reaction catalyzed by morphinone reductase, MR. Quantum mechanical tunneling coefficients, computed in the framework of variational transition-state theory, play a significant role in this reaction, reaching values of 23.8 ± 5.5 for the lightest isotopologue; one of the largest values reported for enzymatic systems. This prediction is supported by the agreement between the theoretically predicted rate constants and the corresponding experimental values. Simulations indicate that the role of protein motions can be satisfactorily described as equilibrium fluctuations along the reaction coordinate, in line with a high degree of preorganization displayed by this enzyme.

Roles of Closed- and Open-Loop Conformations in Large-Scale Structural Transitions of l-Lactate Dehydrogenase

The mechanism of l-lactate generation from pyruvate by l-lactate dehydrogenase (LDH) from the rabbit muscle was studied theoretically by the multistructural microiteration (MSM) method combined with the quantum mechanics/molecular mechanics (QM/MM)-ONIOM method, where the MSM method describes the MM environment as a weighted average of multiple different structures that are fully relaxed during geometry optimization or a reaction path calculation for the QM part. The results showed that the substrate binding and product states were stabilized only in the open-loop conformation of LDH and the reaction occurred in the closed-loop conformation. In other words, before and after the chemical reaction, a large-scale structural transition from the open-loop conformation to the closed-loop conformation and vice versa occurred. The closed-loop conformation stabilized the transition state of the reaction. In contrast, the open-loop conformation stabilized the substrate binding and final states. In other words, the closed- to open-loop transition at the substrate binding state urges capture of the substrate molecule, the subsequent open- to closed-loop transition promotes the product generation, and the final closed- to open-loop transition at the final state prevents the reverse reaction going back to the substrate binding state. It is thus suggested that the exchange of stability between the closed- and open-loop conformations at different states promotes the catalytic cycle.

Dynamic Description of the Catalytic Cycle of Malate Enzyme: Stereoselective Recognition of Substrate, Chemical Reaction, and Ligand Release

In protein engineering, investigations of catalytic cycle facilitate rational design of enzymes. In the present work, deeper analysis on the catalytic cycle of malate enzyme (EC 1.1.1.40), an enzyme involved in cancer metabolic and fatty acid synthesis, was performed. In substrate binding, stereoselective recognition of a substrate originates from distance and angle difference between two chiral substrates and Mn²⁺ as well as monodentate or coplanar ion reaction with Arg165. In catalytic transformation, the activation barrier for the hydride transfer of d-malate is 20.28 kcal/mol higher than that for l-malate. The activation barrier for β-decarboxylation of oxaloacetate is about 4.59 kcal/mol higher than the activation barrier for the hydride transfer of l-malate. The effective activation barrier is 16.44 kcal/mol, which is in close agreement with the value derived from the application of transition-state theory and the Eyring equation to k_cat. In ligand release, l/d-malate needs to overcome a higher barrier than pyruvate to break all bonds in parallel and then to escape from the binding pocket. Leu167 and Asn421 comprise a swinging gate to control the product release. The more open gate is possibly required in the direction of pyruvate to l-malate. Our studies are focused on extending structural knowledge regarding the malate enzyme and provided a powerful strategy for future experimental investigations.

Experimental and Computational Studies Delineate the Role of Asparagine 177 in Hydride Transfer for E. coli Thymidylate Synthase

Thymidylate synthase (TSase), an enzyme responsible for the de novo biosynthesis of 2'-deoxythymidine 5'-monophosphate (thymidylate, dTMP) necessary for DNA synthesis, has been a drug target for decades. TSase is a highly conserved enzyme across species ranging from very primitive organisms to mammals. Among the many conserved active site residues, an asparagine (N177, using Escherichia coli residues numbering) appears to make direct hydrogen bonds with both the C4=O4 carbonyl of the 2'-deoxyuridine 5'-monophosphate (uridylate, dUMP) substrate and its pyrimidine ring's N3. Recent studies have reassessed the TSase catalytic mechanism, focusing on the degree of negative charge accumulation at the O4 carbonyl of the substrate during two critical H-transfers - a proton abstraction and a hydride transfer. To obtain insights into the role of this conserved N177 on the hydride transfer, we examined its aspartic acid (D) and serine (S) mutants - each of which is expected to alter hydrogen bonding and charge stabilization around the C4=O4 carbonyl of the 2'-deoxyuridine 5'-monophosphate (uridylate, dUMP) substrate. Steady-state kinetics, substrate binding order studies and temperature-dependency analysis of intrinsic KIEs for the hydride transfer step of the TSase catalytic cycle suggest the active site of N177D is not precisely organized for that step. A smaller disruption was observed for N177S, which could be rationalized by partial compensation by water molecules and rearrangement of other residues toward preparation of the system for the hydride transfer under study. These experimental findings are qualitatively mirrored by QM/MM computational simulations, thereby shedding light on the sequence and synchronicity of steps in the TSase-catalyzed reaction. This information could potentially inform the design of mechanism-based drugs targeting this enzyme.

Temperature dependence of dynamic, tunnelling and kinetic isotope effects in formate dehydrogenase

The origin of the catalytic power of enzymes has been a question of debate for a long time. In this regard, the possible contribution of protein dynamics in enzymatic catalysis has become one of the most controversial topics. In the present work, the hydride transfer step in the formate dehydrogenase (FDH EC 1.2.1.2) enzyme is studied by means of molecular dynamic (MD) simulations with quantum mechanics/molecular mechanics (QM/MM) potentials in order to explore any correlation between dynamics, tunnelling effects and the rate constant. The temperature dependence of the kinetic isotope effects (KIEs), which is one of the few tests that can be studied by experiments and simulations to shed light on this debate, has been computed and the results have been compared with previous experimental data. The classical mechanical free energy barrier and the number of recrossing trajectories seem to be temperature-independent while the quantum vibrational corrections and the tunnelling effects are slightly temperature-dependent over the interval of 5-45 °C. The computed primary KIEs are in very good agreement with previous experimental data, being almost temperature-independent within the standard deviations. The modest dependence on the temperature is due to just the quantum vibrational correction contribution. These results, together with the analysis of the evolution of the collective variables such as the electrostatic potential or the electric field created by the protein on the key atoms involved in the reaction, confirm that while the protein is well preorganised, some changes take place along the reaction that favour the hydride transfer and the product release. Coordinates defining these movements are, in fact, part of the real reaction coordinate.

Promoting Vibrations and the Function of Enzymes. Emerging Theoretical and Experimental Convergence

A complete understanding of enzyme catalysis requires knowledge of both transition state features and the detailed motions of atoms that cause reactant molecules to form and traverse the transition state. The seeming intractability of the problem arises from the femtosecond lifetime of chemical transition states, preventing most experimental access. Computational chemistry is admirably suited to short time scale analysis but can be misled by inappropriate starting points or by biased assumptions. Kinetic isotope effects provide an experimental approach to transition state structure and a method for obtaining transition state analogues but, alone, do not inform how that transition state is reached. Enzyme structures with transition state analogues provide computational starting points near the transition state geometry. These well-conditioned starting points, combined with the unbiased computational method of transition path sampling, provide realistic atomistic motions involved in transition state formation and passage. In many, but not all, enzymatic systems, femtosecond local protein motions near the catalytic site are linked to transition state formation. These motions are not inherently revealed by most approaches of transition state theory, because transition state theory replaces dynamics with the statistics of the transition state. Experimental and theoretical convergence of the link between local catalytic site vibrational modes and catalysis comes from heavy atom ("Born-Oppenheimer") enzymes. Fully labeled and catalytic site local heavy atom labels perturb the probability of finding enzymatic transition states in ways that can be analyzed and predicted by transition path sampling. Recent applications of these experimental and computational approaches reveal how subpicosecond local catalytic site protein modes play important roles in creating the transition state.

Trehalose Mediated Inhibition of Lactate Dehydrogenase from Rabbit Muscle. The Application of Kramers' Theory in Enzyme Catalysis

Lactate dehydrogenase (LDH) catalyzes the reduction of pyruvate to lactate by using NADH. LDH kinetics has been proposed to be dependent on the dynamics of a loop over the active site. Kramers' theory has been useful in the study of enzyme catalysis dependent on large structural dynamics. In this work, LDH kinetics was studied in the presence of trehalose and at different temperatures. In the absence of trehalose, temperature increase raised exponentially the LDH V_max and revealed a sigmoid transition of K_m toward a low-affinity state similar to protein unfolding. Notably, LDH V_max diminished when in the presence of trehalose, while pyruvate affinity increased and the temperature-mediated binding site transition was hindered. The effect of trehalose on k_cat was viscosity dependent as described by Kramers' theory since V_max correlated inversely with the viscosity of the medium. As a result, activation energy ( E_a) for pyruvate reduction was dramatically increased by trehalose presence. This work provides experimental evidence that the dynamics of a structural component in LDH is essential for catalysis, i.e., the closing of the loop on the active site. While the trehalose mediated-increased of pyruvate affinity is proposed to be due to the compaction and/or increase of structural order at the binding site.

Reaction Mechanism of Organocatalytic Michael Addition of Nitromethane to Cinnamaldehyde: A Case Study on Catalyst Regeneration and Solvent Effects

The Michael addition of nitromethane to cinnamaldehyde has been computationally studied in the absence of a catalyst and the presence of a biotinylated secondary amine by a combined computational and experimental approach. The calculations were performed at the density functional theory (DFT) level with the M06-2X hybrid functional, and a polarizable continuum model has been employed to mimic the effect of two different solvents: dichloromethane (DCM) and water. Contrary to common assumption, the product-derived iminium intermediate was absent in both of the solvents tested. Instead, hydrating the C1–C2 double bond in the enamine intermediate directly yields the tetrahedral intermediate, which is key for forming the product and regenerating the catalyst. Enamine hydration is concerted and found to be rate-limiting in DCM but segregated into two non-rate-limiting steps when the solvent is replaced with water. However, further analysis revealed that the use of water as solvent also raises the energy barriers for other chemical steps, particularly the critical step of C–C bond formation between the iminium intermediate and nucleophile; this consequently lowers both the reaction yield and enantioselectivity of this LUMO-lowering reaction, as experimentally detected. These findings provide a logical explanation to why water often enhances organocatalysis when used as an additive but hampers the reaction progress when employed as a solvent.

Single Chain Polymeric Nanoparticles to Promote Selective Hydroxylation Reactions of Phenol Catalyzed by Copper

Metal-containing single chain polymeric nanoparticles (SCPNs) can be used as synthetic mimics of metalloenzymes. Currently, the role of the folded polymer backbones on the activity and selectivity of metal sites is not clear. Herein, we report our findings on how polymeric frameworks modulate the coordination of Cu sites and the catalytic activity/selectivity of Cu-containing SCPNs mimicking monophenol hydroxylation reactions. Imidazole-functionalized copolymers of poly(methyl methacrylate-co-3-imidazolyl-2-hydroxy propyl methacrylate) were used for intramolecular Cu-imidazole binding that triggered the self-folding of polymers. Polymer chains imposed steric hindrance which yielded unsaturated Cu sites with an average coordination number of 3.3. Cu-containing SCPNs showed a high selectivity for the hydroxylation reaction of phenol to catechol, >80%, with a turnover frequency of >870 h^-1 at 60 °C. The selectivity was largely influenced by the flexibility of the folded polymer backbone where a more flexible polymer backbone allows the cooperative catalysis of two Cu sites. The second coordination sphere provided by the folded polymer that has been less studied is therefore critical in the design of active mimics of metalloenzymes.

Benchmarking Quantum Mechanics/Molecular Mechanics (QM/MM) Methods on the Thymidylate Synthase-Catalyzed Hydride Transfer

Given the ubiquity of hydride-transfer reactions in enzyme-catalyzed processes, identifying the appropriate computational method for evaluating such biological reactions is crucial to perform theoretical studies of these processes. In this paper, the hydride-transfer step catalyzed by thymidylate synthase (TSase) is studied by examining hybrid quantum mechanics/molecular mechanics (QM/MM) potentials via multiple semiempirical methods and the M06-2X hybrid density functional. Calculations of protium and tritium transfer in these reactions across a range of temperatures allowed calculation of the temperature dependence of kinetic isotope effects (KIE). Dynamics and quantum-tunneling effects are revealed to have little effect on the reaction rate, but are significant in determining the KIEs and their temperature dependence. A good agreement with experiments is found, especially when computed for RM1/MM simulations. The small temperature dependence of quantum tunneling corrections and the quasiclassical contribution term cancel each other, while the recrossing transmission coefficient seems to be temperature-independent over the interval of 5-40 °C.

CADEE: Computer-Aided Directed Evolution of Enzymes

The tremendous interest in enzymes as biocatalysts has led to extensive work in enzyme engineering, as well as associated methodology development. Here, a new framework for computer-aided directed evolution of enzymes (CADEE) is presented which allows a drastic reduction in the time necessary to prepare and analyze in silico semi-automated directed evolution of enzymes. A pedagogical example of the application of CADEE to a real biological system is also presented in order to illustrate the CADEE workflow.

Dynamic and Electrostatic Effects on the Reaction Catalyzed by HIV-1 Protease

HIV-1 Protease (HIV-1 PR) is one of the three enzymes essential for the replication process of HIV-1 virus, which explains why it has been the main target for design of drugs against acquired immunodeficiency syndrome (AIDS). This work is focused on exploring the proteolysis reaction catalyzed by HIV-1 PR, with special attention to the dynamic and electrostatic effects governing its catalytic power. Free energy surfaces for all possible mechanisms have been computed in terms of potentials of mean force (PMFs) within hybrid QM/MM potentials, with the QM subset of atoms described at semiempirical (AM1) and DFT (M06-2X) level. The results suggest that the most favorable reaction mechanism involves formation of a gem-diol intermediate, whose decomposition into the product complex would correspond to the rate-limiting step. The agreement between the activation free energy of this step with experimental data, as well as kinetic isotope effects (KIEs), supports this prediction. The role of the protein dynamic was studied by protein isotope labeling in the framework of the Variational Transition State Theory. The predicted enzyme KIEs, also very close to the values measured experimentally, reveal a measurable but small dynamic effect. Our calculations show how the contribution of dynamic effects to the effective activation free energy appears to be below 1 kcal·mol^-1. On the contrary, the electric field created by the protein in the active site of the enzyme emerges as being critical for the electronic reorganization required during the reaction. These electrostatic properties of the active site could be used as a mold for future drug design.

Triple Isotope Effects Support Concerted Hydride and Proton Transfer and Promoting Vibrations in Human Heart Lactate Dehydrogenase

Transition path sampling simulations have proposed that human heart lactate dehydrogenase (LDH) employs protein promoting vibrations (PPVs) on the femtosecond (fs) to picosecond (ps) time scale to promote crossing of the chemical barrier. This chemical barrier involves both hydride and proton transfers to pyruvate to form l-lactate, using reduced nicotinamide adenine dinucleotide (NADH) as the cofactor. Here we report experimental evidence from three types of isotope effect experiments that support coupling of the promoting vibrations to barrier crossing and the coincidence of hydride and proton transfer. We prepared the native (light) LDH and a heavy LDH labeled with ¹³C, ¹⁵N, and nonexchangeable ²H (D) to perturb the predicted PPVs. Heavy LDH has slowed chemistry in single turnover experiments, supporting a contribution of PPVs to transition state formation. Both the [4-²H]NADH (NADD) kinetic isotope effect and the D₂O solvent isotope effect were increased in dual-label experiments combining both NADD and D₂O, a pattern maintained with both light and heavy LDHs. These isotope effects support concerted hydride and proton transfer for both light and heavy LDHs. Although the transition state barrier-crossing probability is reduced in heavy LDH, the concerted mechanism of the hydride-proton transfer reaction is not altered. This study takes advantage of triple isotope effects to resolve the chemical mechanism of LDH and establish the coupling of fs-ps protein dynamics to barrier crossing.

The influence of active site conformations on the hydride transfer step of the thymidylate synthase reaction mechanism

The hydride transfer from C6 of tetrahydrofolate to the reaction's exocyclic methylene-dUMP intermediate is the rate limiting step in thymidylate synthase (TSase) catalysis. This step has been studied by means of QM/MM molecular dynamics simulations to generate the corresponding free energy surfaces. The use of two different initial X-ray structures has allowed exploring different conformational spaces and the existence of chemical paths with not only different reactivities but also different reaction mechanisms. The results confirm that this chemical conversion takes place preferentially via a concerted mechanism where the hydride transfer is conjugated to thiol-elimination from the product. The findings also confirm the labile character of the substrate-enzyme covalent bond established between the C6 of the nucleotide substrate and a conserved cysteine residue. The calculations also reproduce and rationalize a normal H/T 2° kinetic isotope effect measured for that step. From a computational point of view, the results demonstrate that the use of an incomplete number of coordinates to describe the real reaction coordinate can render biased results.

Hydrostatic Pressure Studies Distinguish Global from Local Protein Motions in C-H Activation by Soybean Lipoxygenase-1

The proposed contributions of distinct classes of local versus global protein motions during enzymatic bond making/breaking processes has been difficult to verify. We employed soybean lipoxygenase-1 as a model system to investigate the impact of high pressure at variable temperatures on the hydrogen-tunneling properties of the wild-type protein and three single-site mutants. For all variants, pressure dramatically elevates the enthalpies of activation for the C-H activation. In contrast, the primary kinetic isotope effects (KIEs) for C-H activation and their corresponding temperature dependencies remain unchanged up to ca. 700 bar. The differential impact of elevated hydrostatic pressure on the temperature dependencies of rate constants versus substrate KIEs provides direct evidence for two distinct classes of protein motions: local, isotope-dependent donor-acceptor distance-sampling modes, and a more global, isotope-independent search for productive protein conformational sub-states.

QM/MM study of l-lactate oxidation by flavocytochrome b2

Free energy surfaces calculated from a state-of-the-art computational methodology highlight the role of active site residues in l-lactate oxidation by flavocytochrome b₂.

QM/MM study of the mechanism of reduction of 3-hydroxy-3-methylglutaryl coenzyme A catalyzed by human HMG-CoA reductase

Detailing with atomistic resolution the reaction mechanism of human HMG-CoA reductase (HMG-CoA-R) might provide valuable insights for the development of new cholesterol-lowering drugs.

Transition state theory for enzyme kinetics

This article is an essay that discusses the concepts underlying the application of modern transition state theory to reactions in enzymes. Issues covered include the potential of mean force, the quantization of vibrations, the free energy of activation, and transmission coefficients to account for nonequilibrium effect, recrossing, and tunneling.