Evolution of the Chemical Step in Enzyme Catalysis

Engineering Glucosamine-6-Phosphate Synthase to Achieve Efficient One-Step Biosynthesis of Glucosamine

As an important functional monosaccharide, glucosamine (GlcN) is widely used in fields such as medicine, food nutrition, and health care. Here, we report a distinct GlcN biosynthesis method that utilizes engineered Bacillus subtilis glucosamine-6-phosphate synthase (BsGlmS) to convert D-fructose to directly generate GlcN. The best variant obtained by using a combinatorial active-site saturation test/iterative saturation mutagenesis (CAST/ISM) strategy was a quadruple mutant S596D/V597G/S347H/G299Q (BsGlmS-BK19), which has a catalytic activity 1736-fold that of the wild type toward D-fructose. Upon using mutant BK19 as a whole-cell catalyst, D-fructose was converted into GlcN with 65.32% conversion in 6 h, whereas the wild type only attained a conversion rate of 0.31% under the same conditions. Molecular docking and molecular dynamics simulations were implemented to provide insights into the mechanism underlying the enhanced activity of BK19. Importantly, the BsGlmS-BK19 variant specifically catalyzes D-fructose without the need for phosphorylated substrates, representing a significant advancement in GlcN biosynthesis.

Experimental Evidence for the Role of Dynamics in pH-Dependent Enzymatic Activity

Enzymatic activity is heavily influenced by pH, but the rationale for the dynamical mechanism of pH-dependent enzymatic activity has not been fully understood. In this work, combined neutron scattering techniques, including quasielastic neutron scattering (QENS) and small angle neutron scattering (SANS), are used to study the structural and dynamic changes of a model enzyme, xylanase, under different pH and temperature environments. The QENS results reveal that xylanase at optimal pH exhibits faster relaxational dynamics and a lower energy barrier between conformational substates. The SANS results demonstrate that pH affects both xylanase's stability and monodispersity. Our findings indicate that enzymes have optimized stability and function under their optimal pH conditions, with both structure and dynamics being affected. The current study offers valuable insights into enzymatic functionality mechanisms, allowing for broad industrial applications.

Impact of N-Glycosylation on Protein Structure and Dynamics Linked to Enzymatic C-H Activation in the M. oryzae Lipoxygenase

Lipoxygenases (LOXs) from pathogenic fungi are potential therapeutic targets for defense against plant and select human diseases. In contrast to the canonical LOXs in plants and animals, fungal LOXs are unique in having appended N-linked glycans. Such important post-translational modifications (PTMs) endow proteins with altered structure, stability, and/or function. In this study, we present the structural and functional outcomes of removing or altering these surface carbohydrates on the LOX from the devastating rice blast fungus, M. oryzae, MoLOX. Alteration of the PTMs did notinfluence the active site enzyme-substrate ground state structures as visualized by electron-nuclear double resonance (ENDOR) spectroscopy. However, removal of the eight N-linked glycans by asparagine-to-glutamine mutagenesis nonetheless led to a change in substrate selectivity and an elevated activation energy for the reaction with substrate linoleic acid, as determined by kinetic measurements. Comparative hydrogen-deuterium exchange mass spectrometry (HDX-MS) analysis of wild-type and Asn-to-Gln MoLOX variants revealed a regionally defined impact on the dynamics of the arched helix that covers the active site. Guided by these HDX results, a single glycan sequon knockout was generated at position 72, and its comparative substrate selectivity from kinetics nearly matched that of the Asn-to-Gln variant. The cumulative data from model glyco-enzyme MoLOX showcase how the presence, alteration, or removal of even a single N-linked glycan can influence the structural integrity and dynamics of the protein that are linked to an enzyme's catalytic proficiency, while indicating that extensive glycosylation protects the enzyme during pathogenesis by protecting it from protease degradation.

Study of the Effects of Remote Heavy Group Vibrations on the Temperature Dependence of Hydride Kinetic Isotope Effects of the NADH/NAD+ Model Reactions

It has recently been observed that the temperature(T)-dependence of KIEs in H-tunneling reactions, characterized by isotopic activation energy difference (ΔEa = EaD - EaH), is correlated to the rigidity of the tunneling ready states (TRSs) in enzymes. A more rigid system with narrowly distributed H-donor-acceptor distances (DADs) at the TRSs gives rise to a weaker T-dependence of KIEs (i.e., a smaller ΔEa). Theoreticians have attempted to develop new H-tunneling models to explain this, but none has been universally accepted. In order to further understand the observations in enzymes and provide useful data to build new theoretical models, we have studied the electronic and solvent effects on ΔEa's for the hydride-tunneling reactions of NADH/NAD+ analogues. We found that a tighter charge-transfer (CT) complex system gives rises to a smaller ΔEa, consistent with the enzyme observations. In this paper, we use the remote heavy group (R) vibrational effects to mediate the system rigidity to study the rigidity-ΔEa relationship. The specific hypothesis is that slower vibrations of a heavier remote group would broaden the DAD distributions and increase the ΔEa value. Four NADH/NAD+ systems were studied in acetonitrile but most of such heavy group vibrations do not appear to significantly increase the ΔEa. The remote heavy group vibrations in these systems may have not affected the CT complexation rigidity to a degree that can significantly increase the DADs, and further, the ΔEa values.

Identification of the Thermal Activation Network in Human 15-Lipoxygenase-2: Divergence from Plant Orthologs and Its Relationship to Hydrogen Tunneling Activation Barriers

The oxidation of polyunsaturated fatty acids by lipoxygenases (LOXs) is initiated by a C-H cleavage step in which the hydrogen atom is transferred quantum mechanically (i.e., via tunneling). In these reactions, protein thermal motions facilitate the conversion of ground-state enzyme-substrate complexes to tunneling-ready configurations and are thus important for transferring energy from the solvent to the active site for the activation of catalysis. In this report, we employed temperature-dependent hydrogen-deuterium exchange mass spectrometry (TDHDX-MS) to identify catalytically linked, thermally activated peptides in a representative animal LOX, human epithelial 15-LOX-2. TDHDX-MS of wild-type 15-LOX-2 was compared to two active site mutations that retain structural stability but have increased activation energies (Ea) of catalysis. The Ea value of one variant, V427L, is implicated to arise from suboptimal substrate positioning by increased active-site side chain rotamer dynamics, as determined by X-ray crystallography and ensemble refinement. The resolved thermal network from the comparative Eas of TDHDX-MS between wild-type and V426A is localized along the front face of the 15-LOX-2 catalytic domain. The network contains a clustering of isoleucine, leucine, and valine side chains within the helical peptides. This thermal network of 15-LOX-2 is different in location, area, and backbone structure compared to a model plant lipoxygenase from soybean that exhibits a low Ea value of catalysis compared to the human ortholog. The presented data provide insights into the divergence of thermally activated protein motions in plant and animal LOXs and their relationships to the enthalpic barriers for facilitating hydrogen tunneling.

Structural Effects on the Temperature Dependence of Hydride Kinetic Isotope Effects of the NADH/NAD+ Model Reactions in Acetonitrile: Charge-Transfer Complex Tightness Is a Key

It has recently frequently been found that the kinetic isotope effect (KIE) is independent of temperature (T) in H-tunneling reactions in enzymes but becomes dependent on T in their mutants. Many enzymologists found that the trend is related to different donor-acceptor distances (DADs) at tunneling-ready states (TRSs), which could be sampled by protein dynamics. That is, a more rigid system of densely populated short DADs gives rise to a weaker T dependence of KIEs. Theoreticians have attempted to develop H-tunneling theories to explain the observations, but none have been universally accepted. It is reasonable to assume that the DAD sampling concept, if it exists, applies to the H-transfer reactions in solution, as well. In this work, we designed NADH/NAD+ model reactions to investigate their structural effects on the T dependence of hydride KIEs in acetonitrile. Hammett correlations together with N-CH3/CD3 secondary KIEs were used to provide the electronic structure of the TRSs and thus the rigidity of their charge-transfer complexation vibrations. In all three pairs of reactions, a weaker T dependence of KIEs always corresponds to a steeper Hammett slope on the substituted hydride acceptors. It was found that a tighter/rigid charge-transfer complexation system corresponds with a weaker T dependence of KIEs, consistent with the observations in enzymes.

Ratcheting synthesis

Synthetic chemistry has traditionally relied on reactions between reactants of high chemical potential and transformations that proceed energetically downhill to either a global or local minimum (thermodynamic or kinetic control). Catalysts can be used to manipulate kinetic control, lowering activation energies to influence reaction outcomes. However, such chemistry is still constrained by the shape of one-dimensional reaction coordinates. Coupling synthesis to an orthogonal energy input can allow ratcheting of chemical reaction outcomes, reminiscent of the ways that molecular machines ratchet random thermal motion to bias conformational dynamics. This fundamentally distinct approach to synthesis allows multi-dimensional potential energy surfaces to be navigated, enabling reaction outcomes that cannot be achieved under conventional kinetic or thermodynamic control. In this Review, we discuss how ratcheted synthesis is ubiquitous throughout biology and consider how chemists might harness ratchet mechanisms to accelerate catalysis, drive chemical reactions uphill and programme complex reaction sequences.

Effect of solvent viscosity on the activation barrier of hydrogen tunneling in the lipoxygenase reaction

Hydrogen tunneling in enzyme reactions has played an important role in linking protein thermal motions to the chemical steps of catalysis. Lipoxygenases (LOXs) have served as model systems for such reactions, showcasing deep hydrogen tunneling mechanisms associated with enzymatic C-H bond cleavage from polyunsaturated fatty acids. Here, we examined the effect of solvent viscosity on the protein thermal motions associated with LOX catalysis using trehalose and glucose as viscogens. Kinetic analysis of the reaction of the paradigm plant orthologue, soybean lipoxygenase (SLO), with linoleic acid revealed no effect on the first-order rate constants, kcat, or activation energy, Ea. Further studies of SLO active site mutants displaying varying Eas, which have been used to probe catalytically relevant motions, likewise provided no evidence for viscogen-dependent motions. Kinetic analyses were extended to a representative fungal LOX from M. oryzae, MoLOX, and a human LOX, 15-LOX-2. While MoLOX behaved similarly to SLO, we show that viscogens inhibit 15-LOX-2 activity. The latter implicates viscogen sensitive, conformational motions in animal LOX reactions. The data provide insight into the role of water hydration layers in facilitating hydrogen (quantum) tunneling in LOX.

Correlation of temperature dependence of hydride kinetic isotope effects with donor-acceptor distances in two solvents of different polarities

Recently observed nearly temperature (T)-independent kinetic isotope effects (KIEs) in wild-type enzymes and T-dependent KIEs in variants were used to suggest that H-tunneling in enzymes is assisted by the fast protein vibrations that help sample short donor-acceptor distances (DADs). This supports the recently proposed role of protein vibrations in DAD sampling catalysis. However, use of T-dependence of KIEs to suggest DAD sampling associated with protein vibrations is debated. We have formulated a hypothesis regarding the correlation and designed experiments in solution to investigate it. The hypothesis is, a more rigid system with shorter DADTRS's at the tunneling ready states (TRSs) gives rise to a weaker T-dependence of KIEs, i.e., a smaller ΔEa (= EaD - EaH). In a former work, the solvent effects of acetonitrile versus chloroform on the ΔEa of NADH/NAD+ model reactions were determined, and the DADPRC's of the productive reactant complexes (PRCs) were computed to substitute the DADTRS for the DADTRS-ΔEa correlation study. A smaller ΔEa was found in the more polar acetonitrile where the positively charged PRC is better solvated and has a shorter DADPRC, indirectly supporting the hypothesis. In this work, the TRS structures of different DADTRS's for the hydride tunneling reaction from 1,3-dimethyl-2-phenylimidazoline to 10-methylacridinium were computed. The N-CH3/CD3 secondary KIEs on both reactants were calculated and fitted to the observed values to find the DADTRS order in both solutions. It was found that the equilibrium DADTRS is shorter in acetonitrile than in chloroform. Results directly support the DADTRS-ΔEa correlation hypothesis as well as the explanation that links T-dependence of KIEs to DAD sampling catalysis in enzymes.

Second-Shell Residues Contribute to Catalysis by Predominately Preorganizing the Apo State in PafA

Residues beyond the first coordination shell are often observed to make considerable cumulative contributions in enzymes. Due to typically indirect perturbations of multiple physicochemical properties of the active site, however, their individual and specific roles in enzyme catalysis and disease-causing mutations remain difficult to predict and understand at the molecular level. Here we analyze the contributions of several second-shell residues in phosphate-irrepressible alkaline phosphatase of flavobacterium (PafA), a representative system as one of the most efficient enzymes. By adopting a multifaceted approach that integrates quantum-mechanical/molecular-mechanical free energy computations, molecular-mechanical molecular dynamics simulations, and density functional theory cluster model calculations, we probe the rate-limiting phosphoryl transfer step and structural properties of all relevant enzyme states. In combination with available experimental data, our computational results show that mutations of the studied second-shell residues impact catalytic efficiency mainly by perturbation of the apo state and therefore substrate binding, while they do not affect the ground state or alter the nature of phosphoryl transfer transition state significantly. Several second-shell mutations also modulate the active site hydration level, which in turn influences the energetics of phosphoryl transfer. These mechanistic insights also help inform strategies that may improve the efficiency of enzyme design and engineering by going beyond the current focus on the first coordination shell.

Solvent Effects on the Temperature Dependence of Hydride Kinetic Isotope Effects: Correlation to the Donor-Acceptor Distances

Protein structural effects on the temperature (T) dependence of kinetic isotope effects (KIEs) in H-tunneling reactions have recently been used to discuss about the role of enzyme thermal motions in catalysis. Frequently observed nearly T-independent KIEs in the wild-type enzymes and T-dependent KIEs in variants suggest that H-tunneling in the former is assisted by the naturally evolved protein constructive vibrations that help sample short donor-acceptor distances (DADs) needed. This explanation that correlates the T-dependence of KIEs with DAD sampling has been highly debated as simulations following other H-tunneling models sometimes gave alternative explanations. In this paper, solvent effects on the T-dependence of KIEs of two hydride tunneling reactions of NADH/NAD⁺ analogues (represented by ΔE_a = E_aD - E_aH) were determined in attempts to replicate the observations in enzymes and test the protein vibration-assisted DAD sampling concept. Effects of selected aprotic solvents on the DAD_PRC's of the productive reactant complexes (PRCs) and the DAD_TRS's of the activated tunneling ready states (TRSs) were obtained through computations and analyses of the kinetic data, including 2° KIEs, respectively. A weaker T-dependence of KIEs (i.e., smaller ΔE_a) was found in a more polar aprotic solvent in which the system has a shorter average DAD_PRC and DAD_TRS. Further results show that a charge-transfer (CT) complexation made of a stronger donor/acceptor gives rise to a smaller ΔE_a. Overall, the shorter and less broadly distributed DADs resulting from the stronger CT complexation vibrations give rise to a smaller ΔE_a. Our results appear to support the explanation that links the T-dependence of KIEs to the donor-acceptor rigidity in enzymes.

Functional roles of enzyme dynamics in accelerating active site chemistry: Emerging techniques and changing concepts

With the growing acceptance of the contribution of protein conformational ensembles to enzyme catalysis and regulation, research in the field of protein dynamics has shifted toward an understanding of the atomistic properties of protein dynamical networks and the mechanisms and time scales that control such behavior. A full description of an enzymatic reaction coordinate is expected to extend beyond the active site and include site-specific networks that communicate with the protein/water interface. Advances in experimental tools for the spatial resolution of thermal activation pathways are being complemented by biophysical methods for visualizing dynamics in real time. An emerging multidimensional model integrates the impacts of bound substrate/effector on the distribution of protein substates that are in rapid equilibration near room temperature with reaction-specific protein embedded heat transfer conduits.

Enabling Role of Ligand-Driven Conformational Changes in Enzyme Evolution

Many enzymes that show a large specificity in binding the enzymatic transition state with a higher affinity than the substrate utilize substrate binding energy to drive protein conformational changes to form caged substrate complexes. These protein cages provide strong stabilization of enzymatic transition states. Using part of the substrate binding energy to drive the protein conformational change avoids a similar strong stabilization of the Michaelis complex and irreversible ligand binding. A seminal step in the development of modern enzyme catalysts was the evolution of enzymes that couple substrate binding to a conformational change. These include enzymes that function in glycolysis (triosephosphate isomerase), the biosynthesis of lipids (glycerol phosphate dehydrogenase), the hexose monophosphate shunt (6-phosphogluconate dehydrogenase), and the mevalonate pathway (isopentenyl diphosphate isomerase), catalyze the final step in the biosynthesis of pyrimidine nucleotides (orotidine monophosphate decarboxylase), and regulate the cellular levels of adenine nucleotides (adenylate kinase). The evolution of enzymes that undergo ligand-driven conformational changes to form active protein-substrate cages is proposed to proceed by selection of variants, in which the selected side chain substitutions destabilize a second protein conformer that shows compensating enhanced binding interactions with the substrate. The advantages inherent to enzymes that incorporate a conformational change into the catalytic cycle provide a strong driving force for the evolution of flexible protein folds such as the TIM barrel. The appearance of these folds represented a watershed event in enzyme evolution that enabled the rapid propagation of enzyme activities within enzyme superfamilies.

Hydride Transfer Mechanism of Enzymatic Sugar Nucleotide C2 Epimerization Probed with a Loose-Fit CDP-Glucose Substrate

Transient oxidation–reduction through hydride transfer with tightly bound NAD coenzyme is used by a large class of sugar nucleotide epimerases to promote configurational inversion of carbon stereocenters in carbohydrate substrates. A requirement for the epimerases to coordinate hydride abstraction and re-addition with substrate rotation in the binding pocket poses a challenge for dynamical protein conformational selection linked to enzyme catalysis. Here, we studied the thermophilic C2 epimerase from Thermodesulfatator atlanticus (TaCPa2E) in combination with a slow CDP-glucose substrate (k_cat ≈ 1.0 min^–1; 60 °C) to explore the sensitivity of the enzymatic hydride transfer toward environmental fluctuations affected by temperature (20–80 °C). We determined noncompetitive primary kinetic isotope effects (KIE) due to ²H at the glucose C2 and showed that a normal KIE on the k_cat (^Dk_cat) reflects isotope sensitivity of the hydrogen abstraction to enzyme-NAD⁺ in a rate-limiting transient oxidation. The ^Dk_cat peaked at 40 °C was 6.1 and decreased to 2.1 at low (20 °C) and 3.3 at high temperature (80 °C). The temperature profiles for k_cat with the ¹H and ²H substrate showed a decrease in the rate below a dynamically important breakpoint (∼40 °C), suggesting an equilibrium shift to an impaired conformational landscape relevant for catalysis in the low-temperature region. Full Marcus-like model fits of the rate and KIE profiles provided evidence for a high-temperature reaction via low-frequency conformational sampling associated with a broad distribution of hydride donor–acceptor distances (long-distance population centered at 3.31 ± 0.02 Å), only poorly suitable for quantum mechanical tunneling. Collectively, dynamical characteristics of TaCPa2E-catalyzed hydride transfer during transient oxidation of CDP-glucose reveal important analogies to mechanistically simpler enzymes such as alcohol dehydrogenase and dihydrofolate reductase. A loose-fit substrate (in TaCPa2E) resembles structural variants of these enzymes by extensive dynamical sampling to balance conformational flexibility and catalytic efficiency.

Engineered Tryptophan Synthase Balances Equilibrium Effects and Fast Dynamic Effects

Creating efficient and stable enzymes for catalysis in pharmaceutical and industrial laboratories is an important research goal. Arnold et al. used directed evolution to engineer a natural tryptophan synthase to create a mutant that is operable under laboratory conditions without the need for a natural allosteric effector. The use of directed evolution allows researchers to improve enzymes without understanding the structure-activity relationship. Here, we present a transition path sampling study of a key chemical transformation in the tryptophan synthase catalytic cycle. We observed that while directed evolution does mimic the natural allosteric effect from a stability perspective, fast protein dynamics associated with chemistry has been dramatically altered. This work provides further evidence of the role of protein dynamics in catalysis and clearly demonstrates the multifaceted complexity of mutations associated with protein engineering. This study also demonstrates a fascinating contrast between allosteric and stand-alone functions at the femtosecond time scale.

Evolution of Optimized Hydride Transfer Reaction and Overall Enzyme Turnover in Human Dihydrofolate Reductase

Evolution of dihydrofolate reductase (DHFR) has been studied using the enzyme from Escherichia coli DHFR (ecDHFR) as a model, but less studies have used the enzyme from Homo sapiens DHFR (hsDHFR). Each enzyme maintains a short and narrow distribution of hydride donor-acceptor distances (DAD) at the tunneling ready state (TRS). Evolution of the enzyme was previously studied in ecDHFR where three key sites were identified as important to the catalyzed reaction. The corresponding sites in hsDHFR are F28, 62-PEKN, and 26-PPLR. Each of these sites was studied here through the creation of mutant variants of the enzyme and measurements of the temperature dependence of the intrinsic kinetic isotope effects (KIEs) on the reaction. F28 is mutated first to M (F28M) and then to the L of the bacterial enzyme (F28L). The KIEs of the F28M variant are larger and more temperature-dependent than wild-type (WT), suggesting a broader and longer average DAD at the TRS. To more fully mimic ecDHFR, we also study a triple mutant of the human enzyme (F32L-PP26N-PEKN62G). Remarkably, the intrinsic KIEs, while larger in magnitude, are temperature-independent like the WT enzymes. We also construct deletion mutations of hsDHFR removing both the 62-PEKN and 26-PPLR sequences. The results mirror those described previously for insertion mutants of ecDHFR. Taken together, these results suggest a balancing act during DHFR evolution between achieving an optimal TRS for hydride transfer and preventing product inhibition arising from the different intercellular pools of NADPH and NADP⁺ in prokaryotic and eukaryotic cells.