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.

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, Δ E_a = E_a(D) - E_a(H), serves as a critical, and often misinterpreted (or even ignored) descriptor of the reaction coordinate. As will be highlighted in this Account, Δ E_a 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 Δ E_a 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 Δ E_a. 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 Δ E_a 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 Δ E_a, 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 Δ E_a. 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 Δ E_a 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.

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.

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 K_d 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.

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.

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 ¹³C isotopologues. In this report, we present synthetic strategies for site-specific ¹³C labeled linoleic acid substrates. We take advantage of a Corey-Fuchs formyl to terminal ¹³C-labeled alkyne conversion, using ¹³CBr₄ as the labeling source, to reduce the number of steps from a previous fatty acid ¹³C 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.

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.

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 O₂ 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.

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, ¹³C 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, ΔE_a, with the latter increasing as the fraction of conformer a decreases. As proposed, ΔE_a 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.

The rate-limiting chemical reaction catalyzed by soybean lipoxygenase (SLO) involves quantum mechanical tunneling of a hydrogen atom from substrate to its active site ferric-hydroxide cofactor. SLO has emerged as a prototypical system for linking the thermal activation of a protein scaffold to the efficiency of active site chemistry. Significantly, hydrogen-deuterium exchange-mass spectrometry (HDX-MS) experiments on wild type and mutant forms of SLO have uncovered trends in the enthalpic barriers for HDX within a solvent-exposed loop (positions 317-334) that correlate well with trends in the corresponding enthalpic barriers for kcat. A model for this behavior posits that collisions between water and loop 317-334 initiate thermal activation at the protein surface that is then propagated 15-34 Å inward toward the reactive carbon of substrate in proximity to the iron catalyst. In this study, we have prepared protein samples containing cysteine residues either at the tip of the loop 317-334 (Q322C) or on a control loop, 586-603 (S596C). Chemical modification of cysteines with the fluorophore 6-bromoacetyl-2-dimethylaminonaphthalene (Badan, BD) provides site-specific probes for the measurement of fluorescence relaxation lifetimes and Stokes shift decays as a function of temperature. Computational studies indicate that surface water structure is likely to be largely preserved in each sample. While both loops exhibit temperature-independent fluorescence relaxation lifetimes as do the Stokes shifts for S596C-BD, the activation enthalpy for the nanosecond solvent reorganization at Q322C-BD (Ea(ksolv) = 2.8(0.9) kcal/mol)) approximates the enthalpy of activation for catalytic C-H activation (Ea(kcat) = 2.3(0.4) kcal/mol). This study establishes and validates the methodology for measuring rates of rapid local motions at the protein/solvent interface of SLO. These new findings, when combined with previously published correlations between protein motions and the rate-limiting hydride transfer in a thermophilic alcohol dehydrogenase, provide experimental evidence for thermally induced protein quakes as the origin of enthalpic barriers in catalysis.