Hydrogen-Deuterium Exchange within Adenosine Deaminase, a TIM Barrel Hydrolase, Identifies Networks for Thermal Activation of Catalysis

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.

Recent advances in structure-based enzyme engineering for functional reconstruction

Structural information can help engineer enzymes. Usually, specific amino acids in particular regions are targeted for functional reconstruction to enhance the catalytic performance, including activity, stereoselectivity, and thermostability. Appropriate selection of target sites is the key to structure-based design, which requires elucidation of the structure-function relationships. Here, we summarize the mutations of residues in different specific regions, including active center, access tunnels, and flexible loops, on fine-tuning the catalytic performance of enzymes, and discuss the effects of altering the local structural environment on the functions. In addition, we keep up with the recent progress of structure-based approaches for enzyme engineering, aiming to provide some guidance on how to take advantage of the structural information.

Mutexa: A Computational Ecosystem for Intelligent Protein Engineering

Protein engineering holds immense promise in shaping the future of biomedicine and biotechnology. This Review focuses on our ongoing development of Mutexa, a computational ecosystem designed to enable "intelligent protein engineering". In this vision, researchers will seamlessly acquire sequences of protein variants with desired functions as biocatalysts, therapeutic peptides, and diagnostic proteins through a finely-tuned computational machine, akin to Amazon Alexa's role as a versatile virtual assistant. The technical foundation of Mutexa has been established through the development of a database that combines and relates enzyme structures and their respective functions (e.g., IntEnzyDB), workflow software packages that enable high-throughput protein modeling (e.g., EnzyHTP and LassoHTP), and scoring functions that map the sequence-structure-function relationship of proteins (e.g., EnzyKR and DeepLasso). We will showcase the applications of these tools in benchmarking the convergence conditions of enzyme functional descriptors across mutants, investigating protein electrostatics and cavity distributions in SAM-dependent methyltransferases, and understanding the role of nonelectrostatic dynamic effects in enzyme catalysis. Finally, we will conclude by addressing the future steps and fundamental challenges in our endeavor to develop new Mutexa applications that assist the identification of beneficial mutants in protein engineering.

Engineering the Active Site Lid Dynamics to Improve the Catalytic Efficiency of Yeast Cytosine Deaminase

Conformational dynamics is important for enzyme catalysis. However, engineering dynamics to achieve a higher catalytic efficiency is still challenging. In this work, we develop a new strategy to improve the activity of yeast cytosine deaminase (yCD) by engineering its conformational dynamics. Specifically, we increase the dynamics of the yCD C-terminal helix, an active site lid that controls the product release. The C-terminal is extended by a dynamical single α-helix (SAH), which improves the product release rate by up to ~8-fold, and the overall catalytic rate kcat by up to ~2-fold. It is also shown that the kcat increase is due to the favorable activation entropy change. The NMR H/D exchange data indicate that the conformational dynamics of the transition state analog complex increases as the helix is extended, elucidating the origin of the enhanced catalytic entropy. This study highlights a novel dynamics engineering strategy that can accelerate the overall catalysis through the entropy-driven mechanism.

Non-active Site Residue in Loop L4 Alters Substrate Capture and Product Release in d-Arginine Dehydrogenase

Numerous studies demonstrate that enzymes undergo multiple conformational changes during catalysis. The malleability of enzymes forms the basis for allosteric regulation: residues located far from the active site can exert long-range dynamical effects on the active site residues to modulate catalysis. The structure of Pseudomonas aeruginosa d-arginine dehydrogenase (PaDADH) shows four loops (L1, L2, L3, and L4) that span the substrate and the FAD-binding domains. Loop L4 comprises residues 329-336, spanning over the flavin cofactor. The I335 residue on loop L4 is ∼10 Å away from the active site and ∼3.8 Å from N(1)-C(2)═O atoms of the flavin. In this study, we used molecular dynamics and biochemical techniques to investigate the effect of the mutation of I335 to histidine on the catalytic function of PaDADH. Molecular dynamics showed that the conformational dynamics of PaDADH are shifted to a more closed conformation in the I335H variant. In agreement with an enzyme that samples more in a closed conformation, the kinetic data of the I335H variant showed a 40-fold decrease in the rate constant of substrate association (k₁), a 340-fold reduction in the rate constant of substrate dissociation from the enzyme-substrate complex (k₂), and a 24-fold decrease in the rate constant of product release (k₅), compared to that of the wild-type. Surprisingly, the kinetic data are consistent with the mutation having a negligible effect on the reactivity of the flavin. Altogether, the data indicate that the residue at position 335 has a long-range dynamical effect on the catalytic function in PaDADH.

Temporal and spatial resolution of distal protein motions that activate hydrogen tunneling in soybean lipoxygenase

The enzyme soybean lipoxygenase (SLO) provides a prototype for deep tunneling mechanisms in hydrogen transfer catalysis. This work combines room temperature X-ray studies with extended hydrogen-deuterium exchange experiments to define a catalytically-linked, radiating cone of aliphatic side chains that connects an active site iron center of SLO to the protein-solvent interface. Employing eight variants of SLO that have been appended with a fluorescent probe at the identified surface loop, nanosecond fluorescence Stokes shifts have been measured. We report a remarkable identity of the energies of activation (Ea) for the Stokes shifts decay rates and the millisecond C-H bond cleavage step that is restricted to side chain mutants within an identified thermal network. These findings implicate a direct coupling of distal protein motions surrounding the exposed fluorescent probe to active site motions controlling catalysis. While the role of dynamics in enzyme function has been predominantly attributed to a distributed protein conformational landscape, the presented data implicate a thermally initiated, cooperative protein reorganization that occurs on a timescale faster than nanosecond and represents the enthalpic barrier to the reaction of SLO.

Data-driven enzyme engineering to identify function-enhancing enzymes

Identifying function-enhancing enzyme variants is a 'holy grail' challenge in protein science because it will allow researchers to expand the biocatalytic toolbox for late-stage functionalization of drug-like molecules, environmental degradation of plastics and other pollutants, and medical treatment of food allergies. Data-driven strategies, including statistical modeling, machine learning, and deep learning, have largely advanced the understanding of the sequence-structure-function relationships for enzymes. They have also enhanced the capability of predicting and designing new enzymes and enzyme variants for catalyzing the transformation of new-to-nature reactions. Here, we reviewed the recent progresses of data-driven models that were applied in identifying efficiency-enhancing mutants for catalytic reactions. We also discussed existing challenges and obstacles faced by the community. Although the review is by no means comprehensive, we hope that the discussion can inform the readers about the state-of-the-art in data-driven enzyme engineering, inspiring more joint experimental-computational efforts to develop and apply data-driven modeling to innovate biocatalysts for synthetic and pharmaceutical applications.

Biochemical and hydrogen-deuterium exchange studies of the single nucleotide polymorphism Y649C in human platelet 12-lipoxygenase linked to a bleeding disorder

Human platelet 12-lipoxygenase (h12-LOX) is responsible for the formation of oxylipin products that play an important role in platelet aggregation. Single nucleotide polymorphisms (SNPs) of h12-LOX have been implicated in several diseases. In this study, we investigate the structural, dynamical, and functional impact of a h12-LOX SNP that generates a tyrosine-to-cysteine mutation at a buried site (Y649C h12-LOX) and was previously ascribed with reduced levels of 12(S)-hydroxyeicosatetraenoic acid (12S-HETE) production in isolated platelets. Herein, in vitro Michaelis-Menten kinetics show reduced catalytic rates for Y649C compared to WT h12-LOX at physiological or lower temperatures. Both proteins exhibited similar melting temperatures, metal content, and oligomerization state. Liposome binding for both proteins was also dependent upon the presence of calcium, temperature, and liposome composition; however, the Y649C variant was found to have lowered binding capacity to liposomes compared to WT at physiological temperatures. Further, hydrogen-deuterium exchange mass spectrometry (HDX-MS) experiments revealed a regional defined enhancement in the peptide mobility caused by the mutation. This increased instability for the mutation stemmed from a change in an interaction with an arched helix that lines the substrate binding site, located ≥15 Å from the mutation site. Finally, differential scanning calorimetry demonstrated a reduced protein (un)folding enthalpy, consistent with the HDX results. Taken together, these results demonstrate remarkable similarity between the mutant and WT h12-LOX, and yet, subtle changes in activity, membrane affinity and protein stability may be responsible for the significant physiological changes that the Y649C SNP manifests in platelet biology.

Dynamical activation of function in metalloenzymes

Formulations of hydrogen tunneling in enzyme-catalysed C-H activation reactions indicate enthalpic barriers to reaction that are independent of chemical steps and dependent on the protein scaffold. A tool to identify catalytically relevant site-specific protein thermal networks has emerged from temperature-dependent hydrogen deuterium exchange (TDHDX). Focusing on mutant enzyme forms with altered activation energies for catalysis, TDHDX provides a comparative analysis of the impact of mutation on Ea for local protein unfolding. Identified thermal networks appear unrelated to protein scaffold conservation and track to the dictates of the catalysed reaction, including sites for metal binding. The positions of thermal networks provide a framework for further understanding of time-dependent, functionally relevant protein motions. Measurement of nanosecond Stokes shifts at the surface of the thermal network in soybean lipoxygenase yields activation energies that are identical to Ea values measured for kcat . This finding identifies a rapid (> nanosecond), long-range and cooperative structural reorganization as the thermal barrier to catalysis. A model for protein dynamics is put forward that integrates broadly distributed protein conformational sampling with protein embedded thermal networks.

DHFR Mutants Modulate Their Synchronized Dynamics with the Substrate by Shifting Hydrogen Bond Occupancies

Antibiotic resistance is a global health problem in which mutations occurring in functional proteins render drugs ineffective. The working mechanisms of the arising mutants are seldom apparent; a methodology to decipher these mechanisms systematically would render devising therapies to control the arising mutational pathways possible. Here we utilize C_α-C_β bond vector relaxations obtained from moderate length MD trajectories to determine conduits for functionality of the resistance conferring mutants of Escherichia coli dihydrofolate reductase. We find that the whole enzyme is synchronized to the motions of the substrate, irrespective of the mutation introducing gain-of-function or loss-of function. The total coordination of the motions suggests changes in the hydrogen bond dynamics with respect to the wild type as a possible route to determine and classify the mode-of-action of individual mutants. As a result, nine trimethoprim-resistant point mutations arising frequently in evolution experiments are categorized. One group of mutants that display the largest occurrence (L28R, W30G) work directly by modifying the dihydrofolate binding region. Conversely, W30R works indirectly by the formation of the E139-R30 salt bridge which releases energy resulting from tight binding by distorting the binding cavity. A third group (D27E, F153S, I94L) arising as single, resistance invoking mutants in evolution experiment trajectories allosterically and dynamically affects a hydrogen bonding motif formed at residues 59-69-71 which in turn modifies the binding site dynamics. The final group (I5F, A26T, R98P) consists of those mutants that have properties most similar to the wild type; these only appear after one of the other mutants is fixed on the protein structure and therefore display clear epistasis. Thus, we show that the binding event is governed by the entire enzyme dynamics while the binding site residues play gating roles. The adjustments made in the total enzyme in response to point mutations are what make quantifying and pinpointing their effect a hard problem. Here, we show that hydrogen bond dynamics recorded on sub-μs time scales provide the necessary fingerprints to decipher the various mechanisms at play.

Convergence in determining enzyme functional descriptors across Kemp eliminase variants

Molecular simulations have been extensively employed to accelerate biocatalytic discoveries. Enzyme functional descriptors derived from molecular simulations have been leveraged to guide the search for beneficial enzyme mutants. However, the ideal active-site region size for computing the descriptors over multiple enzyme variants remains untested. Here, we conducted convergence tests for dynamics-derived and electrostatic descriptors on 18 Kemp eliminase variants across six active-site regions with various boundary distances to the substrate. The tested descriptors include the root-mean-square deviation of the active-site region, the solvent accessible surface area ratio between the substrate and active site, and the projection of the electric field (EF) on the breaking C-H bond. All descriptors were evaluated using molecular mechanics methods. To understand the effects of electronic structure, the EF was also evaluated using quantum mechanics/molecular mechanics methods. The descriptor values were computed for 18 Kemp eliminase variants. Spearman correlation matrices were used to determine the region size condition under which further expansion of the region boundary does not substantially change the ranking of descriptor values. We observed that protein dynamics-derived descriptors, including RMSD_{active_site} and SASA_ratio, converge at a distance cutoff of 5 Å from the substrate. The electrostatic descriptor, EF_C-H, converges at 6 Å using molecular mechanics methods with truncated enzyme models and 4 Å using quantum mechanics/molecular mechanics methods with whole enzyme model. This study serves as a future reference to determine descriptors for predictive modeling of enzyme engineering.

Temperature-dependent hydrogen deuterium exchange shows impact of analog binding on adenosine deaminase flexibility but not embedded thermal networks

The analysis of hydrogen deuterium exchange by mass spectrometry as a function of temperature and mutation (TDHDX-MS) has emerged as a generic and efficient tool for the spatial resolution of protein networks that are proposed to function in the thermal activation of catalysis. In this work, we extend TDHDX from apo-enzyme structures to protein-ligand complexes. Using adenosine deaminase as a prototype, we compared the impacts of a substrate analog (1-deaza-adenosine or DAA) and a very tight-binding inhibitor/transition state analog (pentostatin) at single and multiple temperatures. At a single temperature, we observed different HDX-MS properties for the two ligands, as expected from their 10⁶-fold differences in strength of binding. By contrast, analogous patterns for TDHDX-MS emerge in the presence of both DAA and pentostatin, indicating similar impacts of either ligand on the enthalpic barriers for local protein unfolding. We extended TDHDX to a function-altering mutant of adenosine deaminase in the presence of pentostatin and revealed a protein thermal network that is highly similar to that previously reported for the apo-enzyme (Gao et al., 2020, JACS 142, 19936-19949). Finally, we discuss the differential impacts of pentostatin binding on overall protein flexibility vs. site-specific thermal transfer pathways in the context of models for substrate-induced changes to a distributed protein conformational landscape that act in synergy with embedded protein thermal networks to achieve efficient catalysis.

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.

Molecular Dynamics-Derived Descriptor Informs the Impact of Mutation on the Catalytic Turnover Number in Lactonase Across Substrates

Molecular dynamics simulations have been extensively employed to reveal the roles of protein dynamics in mediating enzyme catalysis. However, simulation-derived predictive descriptors that inform the impacts of mutations on catalytic turnover numbers remain largely unexplored. In this work, we report the identification of molecular modeling-derived descriptors to predict mutation effect on the turnover number of lactonase SsoPox with both native and non-native substrates. The study consists of 10 enzyme-substrate complexes resulting from a combination of five enzyme variants with two substrates. For each complex, we derived 15 descriptors from molecular dynamics simulations and applied principal component analysis to rank the predictive capability of the descriptors. A top-ranked descriptor was identified, which is the solvent-accessible surface area (SASA) ratio of the substrate to the active site pocket. A uniform volcano-shaped plot was observed in the distribution of experimental activation free energy against the SASA ratio. To achieve efficient lactonase hydrolysis, a non-native substrate-bound enzyme variant needs to involve a similar range of the SASA ratio to the native substrate-bound wild-type enzyme. The descriptor reflects how well the enzyme active site pocket accommodates a substrate for reaction, which has the potential of guiding optimization of enzyme reaction turnover for non-native chemical transformations.

Advances in Hydrogen/Deuterium Exchange Mass Spectrometry and the Pursuit of Challenging Biological Systems

Solution-phase hydrogen/deuterium exchange (HDX) coupled to mass spectrometry (MS) is a widespread tool for structural analysis across academia and the biopharmaceutical industry. By monitoring the exchangeability of backbone amide protons, HDX-MS can reveal information about higher-order structure and dynamics throughout a protein, can track protein folding pathways, map interaction sites, and assess conformational states of protein samples. The combination of the versatility of the hydrogen/deuterium exchange reaction with the sensitivity of mass spectrometry has enabled the study of extremely challenging protein systems, some of which cannot be suitably studied using other techniques. Improvements over the past three decades have continually increased throughput, robustness, and expanded the limits of what is feasible for HDX-MS investigations. To provide an overview for researchers seeking to utilize and derive the most from HDX-MS for protein structural analysis, we summarize the fundamental principles, basic methodology, strengths and weaknesses, and the established applications of HDX-MS while highlighting new developments and applications.

Identification of Thermal Conduits That Link the Protein-Water Interface to the Active Site Loop and Catalytic Base in Enolase

We report here on the salient role of protein mobility in accessing conformational landscapes that enable efficient enzyme catalysis. We are focused on yeast enolase, a highly conserved lyase with a TIM barrel domain and catalytic loop, as part of a larger study of the relationship of site selective protein motions to chemical reactivity within superfamilies. Enthalpically hindered variants were developed by replacement of a conserved hydrophobic side chain (Leu 343) with smaller side chains. Leu343 is proximal to the active site base in enolase, and comparative pH rate profiles for the valine and alanine variants indicate a role for side chain hydrophobicity in tuning the pK_a of the catalytic base. However, the magnitude of a substrate deuterium isotope effect is almost identical for wild-type (WT) and Leu343Ala, supporting an unchanged rate-determining proton abstraction step. The introduced hydrophobic side chains at position 343 lead to a discontinuous break in both activity and activation energy as a function of side chain volume. Hydrogen-deuterium exchange mass spectrometry (HDX-MS) experiments were performed as a function of time and temperature for WT and Leu343Ala, and provide a spatially resolved map of changes in protein flexibility following mutation. Impacts on protein flexibility are localized to specific networks that arise at the protein-solvent interface and terminate in a loop that has been shown by X-ray crystallography to close over the active site. These interrelated effects are discussed in the context of long-range, solvent-accessible and thermally activated networks that play key roles in tuning the precise distances and interactions among reactants.