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.

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.

13C Electron Nuclear Double Resonance Spectroscopy-Guided Molecular Dynamics Computations Reveal the Structure of the Enzyme-Substrate Complex of an Active, N-Linked Glycosylated Lipoxygenase.

Lipoxygenase (LOX) enzymes produce important cell-signaling mediators, yet attempts to capture and characterize LOX-substrate complexes by X-ray co-crystallography are commonly unsuccessful, requiring development of alternative structural methods. We previously reported the structure of the complex of soybean lipoxygenase, SLO, with substrate linoleic acid (LA), as visualized through the integration of 13C/1H electron nuclear double resonance (ENDOR) spectroscopy and molecular dynamics (MD) computations. However, this required substitution of the catalytic mononuclear, nonheme iron by the structurally faithful, yet inactive Mn2+ ion as a spin probe. Unlike canonical Fe-LOXs from plants and animals, LOXs from pathogenic fungi contain active mononuclear Mn2+ metallocenters. Here, we report the ground-state active-site structure of the native, fully glycosylated fungal LOX from rice blast pathogen Magnaporthe oryzae, MoLOX complexed with LA, as obtained through the 13C/1H ENDOR-guided MD approach. The catalytically important distance between the hydrogen donor, carbon-11 (C11), and the acceptor, Mn-bound oxygen, (donor-acceptor distance, DAD) for the MoLOX-LA complex derived in this fashion is 3.4 ± 0.1 Å. The difference of the MoLOX-LA DAD from that of the SLO-LA complex, 3.1 ± 0.1 Å, is functionally important, although is only 0.3 Å, despite the MoLOX complex having a Mn-C11 distance of 5.4 Å and a "carboxylate-out" substrate-binding orientation, whereas the SLO complex has a 4.9 Å Mn-C11 distance and a "carboxylate-in" substrate orientation. The results provide structural insights into reactivity differences across the LOX family, give a foundation for guiding development of MoLOX inhibitors, and highlight the robustness of the ENDOR-guided MD approach to describe LOX-substrate structures.

Investigating membrane-binding properties of lipoxygenases using surface plasmon resonance.

Lipoxygenases (LOXs) catalyze the oxidation of polyunsaturated fatty acids and synthesize oxylipin products that drive important cellular signaling processes in plants and animals. While there has been indirect evidence presented for the interaction of mammalian LOXs with membranes, a quantitative study of the molecular details of LOX-membrane interactions is lacking. Here, we mimicked biological membranes using surface plasmon resonance (SPR) sensor chips derivatized with 2-D planar lipophilic anchors (2D LP) to capture liposomes of varying phospholipid compositions that self-assemble into lipid bilayers on the SPR chip. The sensor chip surfaces were then used to investigate the membrane-binding properties of model LOX enzymes. SPR binding assays displayed reproducible and stable liposome capture to the sensor chip surface that allowed for the detailed characterization of LOX-membrane interactions. Our studies demonstrate a calcium-dependence for the membrane binding activities of coral 8R-LOX and human 15-LOX-2. Furthermore, our data confirm the importance of key membrane insertion loop residues in each of these LOX enzymes for membrane binding activity. Experiments utilizing model plant and human LOXs reveal differences in membrane-binding specificities. Our study establishes and validates a robust SPR-based platform using 2D LP sensor chips that allows for the detailed study of LOX-membrane interactions under different experimental conditions, including altered membrane compositions. Collectively, this investigation improves our overall understanding of LOX-membrane interaction properties, and our SPR-based approach holds potential for future use in the development of LOX-based therapeutics.

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.

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.

Membrane-Mimicking Reverse Micelles for High-Resolution Interfacial Study of Proteins and Membranes.

Despite substantial advances, the study of proteins interacting with membranes remains a significant challenge. While integral membrane proteins have been a major focus of recent efforts, peripheral membrane proteins (PMPs) and their interactions with membranes and lipids have far less high-resolution information available. Their small size and the dynamic nature of their interactions have stalled detailed interfacial study using structural methods like cryo-EM and X-ray crystallography. A major roadblock for the structural analysis of PMP interactions is limitations in membrane models to study the membrane recruited state. Commonly used membrane mimics such as liposomes, bicelles, nanodiscs, and micelles are either very large or composed of non-biological detergents, limiting their utility for the NMR study of PMPs. While there have been previous successes with integral and peripheral membrane proteins, currently employed reverse micelle (RM) compositions are optimized for their inertness with proteins rather than their ability to mimic membranes. Applying more native, membrane-like lipids and surfactants promises to be a valuable advancement for the study of interfacial interactions between proteins and membranes. Here, we describe the development of phosphocholine-based RM systems that mimic biological membranes and are compatible with high-resolution protein NMR. We demonstrate new formulations that are able to encapsulate the model soluble protein, ubiquitin, with minimal perturbations of the protein structure. Furthermore, one formula, DLPC:DPC, allowed the encapsulation of the PMPs glutathione peroxidase 4 (GPx4) and phosphatidylethanolamine-binding protein 1 (PEBP1) and enabled the embedment of these proteins, matching the expected interactions with biological membranes. Dynamic light scattering and small-angle X-ray scattering characterization of the RMs reveals small, approximately spherical, and non-aggregated particles, a prerequisite for protein NMR and other avenues of study. The formulations presented here represent a new tool for the study of elusive PMP interactions and other membrane interfacial investigations.

Development of Transient Recombinant Expression and Affinity Chromatography Systems for Human Fibrinogen.

Fibrin forms the structural scaffold of blood clots and has great potential for biomaterial applications. Creating recombinant expression systems of fibrinogen, fibrin's soluble precursor, would advance the ability to construct mutational libraries that would enable structure-function studies of fibrinogen and expand the utility of fibrin as a biomaterial. Despite these needs, recombinant fibrinogen expression systems, thus far, have relied on the time-consuming creation of stable cell lines. Here we present tests of a transient fibrinogen expression system that can rapidly generate yields of 8-12 mg/L using suspension HEK Expi293TM cells. We report results from two different plasmid systems encoding the fibrinogen cDNAs and two different transfection reagents. In addition, we describe a novel, affinity-based approach to purifying fibrinogen from complex media such as human plasma. We show that using a high-affinity peptide which mimics fibrin's knob 'A' sequence enables the purification of 50-75% of fibrinogen present in plasma. Having robust expression and purification systems of fibrinogen will enable future studies of basic fibrin(ogen) biology, while paving the way for the ubiquitous use of fibrin as a biomaterial.

Mutagenesis, Hydrogen-Deuterium Exchange, and Molecular Docking Investigations Establish the Dimeric Interface of Human Platelet-Type 12-Lipoxygenase.

It was previously shown that human platelet 12S-lipoxygenase (h12-LOX) exists as a dimer; however, the specific structure is unknown. In this study, we create a model of the dimer through a combination of computational methods, experimental mutagenesis, and hydrogen-deuterium exchange (HDX) investigations. Initially, Leu183 and Leu187 were replaced by negatively charged glutamate residues and neighboring aromatic residues were replaced with alanine residues (F174A/W176A/L183E/L187E/Y191A). This quintuple mutant disrupted both the hydrophobic and π-π interactions, generating an h12-LOX monomer. To refine the determinants for dimer formation further, the L183E/L187E mutant was generated and the equilibrium shifted mostly toward the monomer. We then submitted the predicted monomeric structure to protein-protein docking to create a model of the dimeric complex. A total of nine of the top 10 most energetically favorable docking conformations predict a TOP-to-TOP dimeric arrangement of h12-LOX, with the α-helices containing a Leu-rich region (L172, L183, L187, and L194), corroborating our experimental results showing the importance of these hydrophobic interactions for dimerization. This model was supported by HDX investigations that demonstrated the stabilization of four, non-overlapping peptides within helix α2 of the TOP subdomain for wt-h12-LOX, consistent with the dimer interface. Most importantly, our data reveal that the dimer and monomer of h12-LOX have distinct biochemical properties, suggesting that the structural changes due to dimerization have allosteric effects on active site catalysis and inhibitor binding.

Kinetic and structural investigations of novel inhibitors of human epithelial 15-lipoxygenase-2.

Human epithelial 15-lipoxygenase-2 (h15-LOX-2, ALOX15B) is expressed in many tissues and has been implicated in atherosclerosis, cystic fibrosis and ferroptosis. However, there are few reported potent/selective inhibitors that are active ex vivo. In the current work, we report newly discovered molecules that are more potent and structurally distinct from our previous inhibitors, MLS000545091 and MLS000536924 (Jameson et al, PLoS One, 2014, 9, e104094), in that they contain a central imidazole ring, which is substituted at the 1-position with a phenyl moiety and with a benzylthio moiety at the 2-position. The initial three molecules were mixed-type, non-reductive inhibitors, with IC50 values of 0.34 ±â€¯0.05 µM for MLS000327069, 0.53 ±â€¯0.04 µM for MLS000327186 and 0.87 ±â€¯0.06 µM for MLS000327206 and greater than 50-fold selectivity versus h5-LOX, h12-LOX, h15-LOX-1, COX-1 and COX-2. A small set of focused analogs was synthesized to demonstrate the validity of the hits. In addition, a binding model was developed for the three imidazole inhibitors based on computational docking and a co-structure of h15-LOX-2 with MLS000536924. Hydrogen/deuterium exchange (HDX) results indicate a similar binding mode between MLS000536924 and MLS000327069, however, the latter restricts protein motion of helix-α2 more, consistent with its greater potency. Given these results, we designed, docked, and synthesized novel inhibitors of the imidazole scaffold and confirmed our binding mode hypothesis. Importantly, four of the five inhibitors mentioned above are active in an h15-LOX-2/HEK293 cell assay and thus they could be important tool compounds in gaining a better understanding of h15-LOX-2's role in human biology. As such, a suite of similar pharmacophores that target h15-LOX-2 both in vitro and ex vivo are presented in the hope of developing them as therapeutic agents.

The Soybean Lipoxygenase-Substrate Complex: Correlation between the Properties of Tunneling-Ready States and ENDOR-Detected Structures of Ground States.

Hydrogen tunneling in enzymatic C-H activation requires a dynamical sampling among ground-state enzyme-substrate (E-S) conformations, which transiently generates a tunneling-ready state (TRS). The TRS is characterized by a hydrogen donor-acceptor distance (DAD) of 2.7 Å, ∼0.5 Å shorter than the dominant DAD of optimized ground states. Recently, a high-resolution, 13C electron-nuclear double-resonance (ENDOR) approach was developed to characterize the ground-state structure of the complex of the linoleic acid (LA) substrate with soybean lipoxygenase (SLO). The resulting enzyme-substrate model revealed two ground-state conformers with different distances between the target C11 of LA and the catalytically active cofactor [Fe(III)-OH]: the active conformer "a", with a van der Waals DAD of 3.1 Å between C11 and metal-bound hydroxide, and an inactive conformer "b", with a distance that is almost 1 Å longer. Herein, the structure of the E-S complex is examined for a series of six variants in which subtle structural modifications of SLO have been introduced either at a hydrophobic side chain near the bound substrate or at a remote residue within a protein network whose flexibility influences hydrogen transfer. A remarkable correlation is found between the ENDOR-derived population of the active ground-state conformer a and the kinetically derived differential enthalpic barrier for D versus H transfer, ΔEa, with the latter increasing as the fraction of conformer a decreases. As proposed, ΔEa provides a "ruler" for the DAD within the TRS. ENDOR measurements further corroborate the previous identification of a dynamical network coupling the buried active site of SLO to the surface. This study shows that subtle imperfections within the initial ground-state structures of E-S complexes are accompanied by compromised geometries at the TRS.

Comparative kinetic isotope effects on first- and second-order rate constants of soybean lipoxygenase variants uncover a substrate-binding network.

Lipoxygenases are widespread enzymes found in virtually all eukaryotes, including fungi, and, more recently, in prokaryotes. These enzymes act on long-chain polyunsaturated fatty acid substrates (C18 to C20), raising questions regarding how the substrate threads its way from solvent to the active site. Herein, we report a comparison of the temperature dependence of isotope effects on first- and second-order rate constants among single-site variants of the prototypic plant enzyme soybean lipoxygenase-1 substituted at amino acid residues inferred to impact substrate binding. We created 10 protein variants including four amino acid positions, Val-750, Ile-552, Ile-839, and Trp-500, located within a previously proposed substrate portal. The conversion of these bulky hydrophobic side chains to smaller side chains is concluded to increase the mobility of flanking helices, giving rise to increased off rates for substrate dissociation from the enzyme. In this manner, we identified a specific "binding network" that can regulate movement of the substrate from the solvent to the active site. Taken together with our previous findings on C-H and O2 activation of soybean lipoxygenase-1, these results support the emergence of multiple complementary networks within a single protein scaffold that modulate different steps along the enzymatic reaction coordinate.

Fatty Acid Allosteric Regulation of C-H Activation in Plant and Animal Lipoxygenases.

Lipoxygenases (LOXs) catalyze the (per) oxidation of fatty acids that serve as important mediators for cell signaling and inflammation. These reactions are initiated by a C-H activation step that is allosterically regulated in plant and animal enzymes. LOXs from higher eukaryotes are equipped with an N-terminal PLAT (Polycystin-1, Lipoxygenase, Alpha-Toxin) domain that has been implicated to bind to small molecule allosteric effectors, which in turn modulate substrate specificity and the rate-limiting steps of catalysis. Herein, the kinetic and structural evidence that describes the allosteric regulation of plant and animal lipoxygenase chemistry by fatty acids and their derivatives are summarized.

Kinetic Characterization of the C-H Activation Step for the Lipoxygenase from the Pathogenic Fungus Magnaporthe oryzae: Impact of N-Linked Glycosylation.

Lipoxygenases from pathogenic fungi belong to the lipoxygenase family of enzymes, which catalyze C-H activation of polyunsaturated fatty acids to form a diverse set of cell-signaling hydroperoxides. While the lipoxygenase catalytic domains are structurally and functionally similar, these fungal enzymes are decorated with N-linked glycans. The impact of N-linked glycans on the structure and function of these enzymes remains largely unknown. One exemplary system is MoLOX, a lipoxygenase from the fungus Magnaporthe oryzae, that is emerging as an important target for the devastating rice blast disease. Herein, we demonstrate that hydrogen transfer, associated with C-H cleavage of the substrate linoleic acid by MoLOX, is rate-determining and occurs by a hydrogen tunneling mechanism. Using the differential enthalpic barrier for hydrogen and deuterium transfer, ΔEa, as a kinetic reporter of tunneling efficiency, a disproportionate increase in the activation energy for deuterium transfer is observed upon treatment of MoLOX with a peptide:N-glycosidase that cleaves N-linked carbohydrates from the protein. This increased ΔEa is consistent with an impairment of substrate positioning in the enzyme-substrate complex for both the tunneling ready state and the ground state. These results provide new insight into the functional consequences of N-linked glycosylation on lipoxygenase C-H activation and have important implications for MoLOX inhibitor design.

Hydrogen-deuterium exchange reveals long-range dynamical allostery in soybean lipoxygenase.

In lipoxygenases, the topologically conserved C-terminal domain catalyzes the oxidation of polyunsaturated fatty acids, generating an assortment of biologically relevant signaling mediators. Plant and animal lipoxygenases also contain a 100-150-amino acid N-terminal C2-like domain that has been implicated in interactions with isolated fatty acids and at the phospholipid bilayer. These interactions may lead to increased substrate availability and contribute to the regulation of active-site catalysis. Because of a lack of structural information, a molecular understanding of this lipid-protein interaction remains unresolved. Herein, we employed hydrogen-deuterium exchange MS (HDXMS) to spatially resolve changes in protein conformation upon interaction of soybean lipoxygenase with a fatty acid surrogate, oleyl sulfate (OS), previously shown to act at a site separate from the substrate-binding site. Specific, OS-induced conformational changes are detected both at the N-terminal domain and within the substrate portal nearly 30 Å away. Combining previously measured kinetic properties in the presence of OS with its impact on the Kd for linoleic acid substrate binding, we conclude that OS binding brings about an increase in rate constants for both the ingress and egress of substrate. We discuss the role of OS-induced changes in protein flexibility in the context of changes in the mechanism of substrate acquisition.

Biophysical Characterization of a Disabled Double Mutant of Soybean Lipoxygenase: The "Undoing" of Precise Substrate Positioning Relative to Metal Cofactor and an Identified Dynamical Network.

Soybean lipoxygenase (SLO) has served as a prototype for understanding the molecular origin of enzymatic rate accelerations. The double mutant (DM) L546A/L754A is considered a dramatic outlier, due to the unprecedented size and near temperature-independence of its primary kinetic isotope effect, low catalytic efficiency, and elevated enthalpy of activation. To uncover the physical basis of these features, we herein apply three structural probes: hydrogen-deuterium exchange mass spectrometry, room-temperature X-ray crystallography and EPR spectroscopy on four SLO variants (wild-type (WT) enzyme, DM, and the two parental single mutants, L546A and L754A). DM is found to incorporate features of each parent, with the perturbation at position 546 predominantly influencing thermally activated motions that connect the active site to a protein-solvent interface, while mutation at position 754 disrupts the ligand field and solvation near the cofactor iron. However, the expanded active site in DM leads to more active site water molecules and their associated hydrogen bond network, and the individual features from L546A and L754A alone cannot explain the aggregate kinetic properties for DM. Using recently published QM/MM-derived ground-state SLO-substrate complexes for WT and DM, together with the thorough structural analyses presented herein, we propose that the impairment of DM is the combined result of a repositioning of the reactive carbon of linoleic acid substrate with regard to both the iron cofactor and a catalytically linked dynamic region of protein.

Understanding Biological Hydrogen Transfer Through the Lens of Temperature Dependent Kinetic Isotope Effects.

Hydrogen atom transfer (HAT) is a salient feature of many enzymatic C-H cleavage mechanisms. In systems where kinetic isolation of HAT is achieved, selective labeling of substrate with hydrogen isotopes, such as deuterium, enables the determination of intrinsic kinetic isotope effects (KIEs). While the magnitude of the KIE is itself informative, ultimately the size of the temperature dependence of the KIE, Δ Ea = Ea(D) - Ea(H), serves as a critical, and often misinterpreted (or even ignored) descriptor of the reaction coordinate. As will be highlighted in this Account, Δ Ea is one of the most robust parameters to emerge from studies of enzyme catalyzed hydrogen transfer. Kinetic parameters for C-H reactions via HAT can appear consistent with either classical "over-the-barrier" or "Bell-like tunneling correction" models. However, neither of these models is able to explain the observation of near-zero Δ Ea values with many native enzymes that increase upon extrinsic or intrinsic perturbations to function. Instead, a full tunneling model has been developed that can account for the aggregate trends in the temperature dependence of the KIE. This model is reminiscent of Marcus-like theory for electron tunneling, with the additional incorporation of an H atom donor-acceptor distance (DAD) sampling term for effective wave function overlap; the role of the latter term is manifested in the experimentally determined Δ Ea. Three enzyme systems from this laboratory that illustrate different aspects of HAT are presented: taurine dioxygenase, the dual copper ß-monooxygenases, and soybean lipoxygenase (SLO). The latter provides a particularly compelling system for understanding the properties of hydrogen tunneling, showing systematic increases in Δ Ea upon reduction in the size of hydrophobic residues both proximal and distal from the active site iron cofactor. Of note, recent ENDOR-based studies of enzyme-substrate complexes with SLO indicate an increase in DAD for mutants with increased Δ Ea, observations that are inconsistent with "Bell-like correction" models. Overall, the surmounting kinetic and biophysical evidence corroborates a multidimensional approach for understanding HAT, offering a robust mechanistic explanation for the magnitude and trends of the KIE and Δ Ea. Recent DFT and QM/MM computations on SLO are compared to the developed nonadiabatic analytical constructs, providing considerable insight into ground state structures and reactivity. However, QM/MM is unable to readily reproduce the small Δ Ea values characteristic of native enzymes. Future theoretical developments to capture these experimental observations may necessitate a parsing of protein motions for local, substrate deuteration-sensitive modes from isotope-insensitive modes within the larger conformational landscape, in the process providing deeper understanding of how native enzymes have evolved to transiently optimize their active site configurations.

Origins of Enzyme Catalysis: Experimental Findings for C-H Activation, New Models, and Their Relevance to Prevailing Theoretical Constructs.

The physical basis for enzymatic rate accelerations is a subject of great fundamental interest and of direct relevance to areas that include the de novo design of green catalysts and the pursuit of new drug regimens. Extensive investigations of C-H activating systems have provided considerable insight into the relationship between an enzyme's overall structure and the catalytic chemistry at its active site. This Perspective highlights recent experimental data for two members of distinct, yet iconic C-H activation enzyme classes, lipoxygenases and prokaryotic alcohol dehydrogenases. The data necessitate a reformulation of the dominant textbook definition of biological catalysis. A multidimensional model emerges that incorporates a range of protein motions that can be parsed into a combination of global stochastic conformational thermal fluctuations and local donor-acceptor distance sampling. These motions are needed to achieve a high degree of precision with regard to internuclear distances, geometries, and charges within the active site. The available model also suggests a physical framework for understanding the empirical enthalpic barrier in enzyme-catalyzed processes. We conclude by addressing the often conflicting interface between computational and experimental chemists, emphasizing the need for computation to predict experimental results in advance of their measurement.

13C ENDOR Spectroscopy of Lipoxygenase-Substrate Complexes Reveals the Structural Basis for C-H Activation by Tunneling.

In enzymatic C-H activation by hydrogen tunneling, reduced barrier width is important for efficient hydrogen wave function overlap during catalysis. For native enzymes displaying nonadiabatic tunneling, the dominant reactive hydrogen donor-acceptor distance (DAD) is typically ca. 2.7 Å, considerably shorter than normal van der Waals distances. Without a ground state substrate-bound structure for the prototypical nonadiabatic tunneling system, soybean lipoxygenase (SLO), it has remained unclear whether the requisite close tunneling distance occurs through an unusual ground state active site arrangement or by thermally sampling conformational substates. Herein, we introduce Mn2+ as a spin-probe surrogate for the SLO Fe ion; X-ray diffraction shows Mn-SLO is structurally faithful to the native enzyme. 13C ENDOR then reveals the locations of 13C10 and reactive 13C11 of linoleic acid relative to the metal; 1H ENDOR and molecular dynamics simulations of the fully solvated SLO model using ENDOR-derived restraints give additional metrical information. The resulting three-dimensional representation of the SLO active site ground state contains a reactive (a) conformer with hydrogen DAD of ∼3.1 Å, approximately van der Waals contact, plus an inactive (b) conformer with even longer DAD, establishing that stochastic conformational sampling is required to achieve reactive tunneling geometries. Tunneling-impaired SLO variants show increased DADs and variations in substrate positioning and rigidity, confirming previous kinetic and theoretical predictions of such behavior. Overall, this investigation highlights the (i) predictive power of nonadiabatic quantum treatments of proton-coupled electron transfer in SLO and (ii) sensitivity of ENDOR probes to test, detect, and corroborate kinetically predicted trends in active site reactivity and to reveal unexpected features of active site architecture.

Synthesis of Site-Specifically 13C Labeled Linoleic Acids.

Soybean lipoxygenase-1 (SLO-1) catalyzes the C-H abstraction from the reactive carbon (C-11) in linoleic acid as the first and rate-determining step in the formation of alkylhydroperoxides. While previous labeling strategies have focused on deuterium labeling to ascertain the primary and secondary kinetic isotope effects for this reaction, there is an emerging interest and need for selectively enriched 13C isotopologues. In this report, we present synthetic strategies for site-specific 13C labeled linoleic acid substrates. We take advantage of a Corey-Fuchs formyl to terminal 13C-labeled alkyne conversion, using 13CBr4 as the labeling source, to reduce the number of steps from a previous fatty acid 13C synthetic labeling approach. The labeled linoleic acid substrates are useful as nuclear tunneling markers and for extracting active site geometries of the enzyme-substrate complex in lipoxygenase.