Mechanistic insights into glycoside 3-oxidases involved in C-glycoside metabolism in soil microorganisms.

C-glycosides are natural products with important biological activities but are recalcitrant to degradation. Glycoside 3-oxidases (G3Oxs) are recently identified bacterial flavo-oxidases from the glucose-methanol-coline (GMC) superfamily that catalyze the oxidation of C-glycosides with the concomitant reduction of O2 to H2O2. This oxidation is followed by C-C acid/base-assisted bond cleavage in two-step C-deglycosylation pathways. Soil and gut microorganisms have different oxidative enzymes, but the details of their catalytic mechanisms are largely unknown. Here, we report that PsG3Ox oxidizes at 50,000-fold higher specificity (kcat/Km) the glucose moiety of mangiferin to 3-keto-mangiferin than free D-glucose to 2-keto-glucose. Analysis of PsG3Ox X-ray crystal structures and PsG3Ox in complex with glucose and mangiferin, combined with mutagenesis and molecular dynamics simulations, reveal distinctive features in the topology surrounding the active site that favor catalytically competent conformational states suitable for recognition, stabilization, and oxidation of the glucose moiety of mangiferin. Furthermore, their distinction to pyranose 2-oxidases (P2Oxs) involved in wood decay and recycling is discussed from an evolutionary, structural, and functional viewpoint.

Corrigendum to "Unveiling molecular details behind improved activity at neutral to alkaline pH of an engineered DyP-type peroxidase" Comput Struct Biotechnol J, 20 (2022), 3899-3910.

[This corrects the article DOI: 10.1016/j.csbj.2022.07.032.].

Distal Mutations Shape Substrate-Binding Sites during Evolution of a Metallo-Oxidase into a Laccase.

Laccases are in increasing demand as innovative solutions in the biorefinery fields. Here, we combine mutagenesis with structural, kinetic, and in silico analyses to characterize the molecular features that cause the evolution of a hyperthermostable metallo-oxidase from the multicopper oxidase family into a laccase (k cat 273 s-1 for a bulky aromatic substrate). We show that six mutations scattered across the enzyme collectively modulate dynamics to improve the binding and catalysis of a bulky aromatic substrate. The replacement of residues during the early stages of evolution is a stepping stone for altering the shape and size of substrate-binding sites. Binding sites are then fine-tuned through high-order epistasis interactions by inserting distal mutations during later stages of evolution. Allosterically coupled, long-range dynamic networks favor catalytically competent conformational states that are more suitable for recognizing and stabilizing the aromatic substrate. This work provides mechanistic insight into enzymatic and evolutionary molecular mechanisms and spots the importance of iterative experimental and computational analyses to understand local-to-global changes.

Getting Deeper into the Molecular Events of Heme Binding Mechanisms: A Comparative Multi-level Computational Study of HasAsm and HasAyp Hemophores.

Many biological systems obtain their activity by the inclusion of metalloporphyrins into one or several binding pockets. However, decoding the molecular mechanism under which these compounds bind to their receptors is something that has not been widely explored and is a field with open questions. In the present work, we apply computational techniques to unravel and compare the mechanisms of two heme-binding systems, concretely the HasA hemophores from Gram negative bacteria Serratiamarcescens (HasAsm) and Yersinia pestis (HasAyp). Despite the high sequence identity between both systems, the comparison between the X-ray structures of their apo and holo forms suggests different heme-binding mechanisms. HasAyp has extremely similar structures for heme-free and heme-bound forms, while HasAsm presents a very large displacement of a loop that ultimately leads to an additional coordination to the metal with respect to HasAyp. We combined Gaussian accelerated molecular dynamics simulations (GaMDs) in explicit solvent and protein-ligand docking optimized for metalloligands. GaMDs were first carried out on heme-free forms of both hemophores. Then, protein-ligand dockings of the heme were performed on cluster representatives of these simulations and the best poses were then subjected to a new series of GaMDs. A series of analyses reveal the following: (1) HasAyp has a conformational landscape extremely similar between heme-bound and unbound states with no to limited impact on the binding of the cofactor, (2) HasAsm presents as a slightly broader conformational landscape in its apo state but can only visit conformations similar to the X-ray of the holo form when the heme has been bound. Such behavior results from a complex cascade of changes in interactions that spread from the heme-binding pocket to the flexible loop previously mentioned. This study sheds light on the diversity of molecular mechanisms of heme-binding and discusses the weight between the pre-organization of the receptor as well as the induced motions resulting in association.

The evolutionary advantage of an aromatic clamp in plant family 3 glycoside exo-hydrolases.

In the barley ß-D-glucan glucohydrolase, a glycoside hydrolase family 3 (GH3) enzyme, the Trp286/Trp434 clamp ensures ß-D-glucosides binding, which is fundamental for substrate hydrolysis during plant growth and development. We employ mutagenesis, high-resolution X-ray crystallography, and multi-scale molecular modelling methods to examine the binding and conformational behaviour of isomeric ß-D-glucosides during substrate-product assisted processive catalysis that operates in GH3 hydrolases. Enzyme kinetics reveals that the W434H mutant retains broad specificity, while W434A behaves as a strict (1,3)-ß-D-glucosidase. Investigations of reactant movements on the nanoscale reveal that processivity is sensitive to mutation-specific alterations of the tryptophan clamp. While wild-type and W434H utilise a lateral cavity for glucose displacement and sliding of (1,3)-linked hydrolytic products through the catalytic site without dissociation, consistent with their high hydrolytic rates, W434A does not adopt processive catalysis. Phylogenomic analyses of GH3 hydrolases disclose the evolutionary advantage of the tryptophan clamp that confers broad specificity, high catalytic efficiency, and processivity.

Unveiling molecular details behind improved activity at neutral to alkaline pH of an engineered DyP-type peroxidase.

DyP-type peroxidases (DyPs) are microbial enzymes that catalyze the oxidation of a wide range of substrates, including synthetic dyes, lignin-derived compounds, and metals, such as Mn2+ and Fe2+, and have enormous biotechnological potential in biorefineries. However, many questions on the molecular basis of enzyme function and stability remain unanswered. In this work, high-resolution structures of PpDyP wild-type and two engineered variants (6E10 and 29E4) generated by directed evolution were obtained. The X-ray crystal structures revealed the typical ferredoxin-like folds, with three heme access pathways, two tunnels, and one cavity, limited by three long loops including catalytic residues. Variant 6E10 displays significantly increased loops' flexibility that favors function over stability: despite the considerably higher catalytic efficiency, this variant shows poorer protein stability compared to wild-type and 29E4 variants. Constant-pH MD simulations revealed a more positively charged microenvironment near the heme pocket of variant 6E10, particularly in the neutral to alkaline pH range. This microenvironment affects enzyme activity by modulating the pK a of essential residues in the heme vicinity and should account for variant 6E10 improved activity at pH 7-8 compared to the wild-type and 29E4 that show optimal enzymatic activity close to pH 4. Our findings shed light on the structure-function relationships of DyPs at the molecular level, including their pH-dependent conformational plasticity. These are essential for understanding and engineering the catalytic properties of DyPs for future biotechnological applications.

Accounting for the instantaneous disorder in the enzyme-substrate Michaelis complex to calculate the Gibbs free energy barrier of an enzyme reaction.

Many enzyme reactions present instantaneous disorder. These dynamic fluctuations in the enzyme-substrate Michaelis complexes generate a wide range of energy barriers that cannot be experimentally observed, but that determine the measured kinetics of the reaction. These individual energy barriers can be calculated using QM/MM methods, but then the problem is how to deal with this dispersion of energy barriers to provide kinetic information. So far, the most usual procedure has implied the so-called exponential average of the energy barriers. In this paper, we discuss the foundations of this method, and we use the free energy perturbation theory to derive an alternative equation to get the Gibbs free energy barrier of the enzyme reaction. In addition, we propose a practical way to implement it. We have chosen four enzyme reactions as examples. In particular, we have studied the hydrolysis of a glycosidic bond catalyzed by the enzyme Thermus thermophilus ß-glycosidase, and the mutant Y284P Ttb-gly, and the hydrogen abstraction reactions from C13 and C7 of arachidonic acid catalyzed by the enzyme rabbit 15-lipoxygenase-1.

Computational modeling of carbohydrate processing enzymes reactions.

Carbohydrate processing enzymes are of biocatalytic interest. Glycoside hydrolases and the recently discovered lytic polysaccharide monooxygenase for their use in biomass degradation to obtain biofuels or valued chemical entities. Glycosyltransferases or engineered glycosidases and phosphorylases for the synthesis of carbohydrates and glycosylated products. Quantum mechanics-molecular mechanics (QM/MM) methods are highly contributing to establish their different chemical reaction mechanisms. Other computational methods are also used to study enzyme conformational changes, ligand pathways, and processivity, e.g. for processive glycosidases like cellobiohydrolases. There is still a long road to travel to fully understand the role of conformational dynamics in enzyme activity and also to disclose the variety of reaction mechanisms these enzymes employ. Additionally, computational tools for enzyme engineering are beginning to be applied to evaluate substrate specificity or aid in the design of new biocatalysts with increased thermostability or tailored activity, a growing field where molecular modeling is finding its way.

Synthesis of polyfluoroalkyl sp2-iminosugar glycolipids and evaluation of their immunomodulatory properties towards anti-tumor, anti-leishmanial and anti-inflammatory therapies.

Immunomodulatory glycolipids, among which α-galactosylceramide (KRN7000) is an iconic example, have shown strong therapeutic potential in a variety of conditions ranging from cancer and infection to autoimmune or neurodegenerative diseases. A main difficulty for those channels is that they often provoke a cytokine storm comprising both pro- and anti-inflammatory mediators that antagonize each other and negatively affect the immune response. The synthesis of analogues with narrower cytokine secretion-inducing capabilities is hampered by the intrinsic difficulty at controlling the stereochemical outcome in glycosidation reactions, particularly if targeting the α-anomer, which seriously hampers drug optimization strategies. Here we show that replacing the monosaccharide glycone by a sp2-iminosugar glycomimetic moiety allows accessing N-linked sp2-iminosugar glycolipids (sp2-IGLs) with total α-stereocontrol in a single step with no need of protecting groups or glycosidation promotors. The lipid tail has been then readily tailored by incorporating polyfluoroalkyl segments of varied lengths in view of favouring binding to the lipid binding site of the master p38 mitogen activated protein kinase (p38 MAPK), thereby polarizing the immune response in a cell-context dependent manner. The compounds have been evaluated for their antiproliferative, anti-leishmanial and anti-inflammatory activities in different cell assays. The size of the fluorous segment was found to be critical for the biological activity, probably by regulating the aggregation and membrane-crossing properties, whereas the hydroxylation profile (gluco or galacto-like) was less relevant. Biochemical and computational data further support a mechanism of action implying binding to the allosteric lipid binding site of p38 MAPK and subsequent activation of the noncanonical autophosphorylation route. The ensemble of results provide a proof of concept of the potential of sp2-IGLs as immunoregulators.

Defective AMH signaling disrupts GnRH neuron development and function and contributes to hypogonadotropic hypogonadism.

Congenital hypogonadotropic hypogonadism (CHH) is a condition characterized by absent puberty and infertility due to gonadotropin releasing hormone (GnRH) deficiency, which is often associated with anosmia (Kallmann syndrome, KS). We identified loss-of-function heterozygous mutations in anti-Müllerian hormone (AMH) and its receptor, AMHR2, in 3% of CHH probands using whole-exome sequencing. We showed that during embryonic development, AMH is expressed in migratory GnRH neurons in both mouse and human fetuses and unconvered a novel function of AMH as a pro-motility factor for GnRH neurons. Pathohistological analysis of Amhr2-deficient mice showed abnormal development of the peripheral olfactory system and defective embryonic migration of the neuroendocrine GnRH cells to the basal forebrain, which results in reduced fertility in adults. Our findings highlight a novel role for AMH in the development and function of GnRH neurons and indicate that AMH signaling insufficiency contributes to the pathogenesis of CHH in humans.

Discovery of processive catalysis by an exo-hydrolase with a pocket-shaped active site.

Substrates associate and products dissociate from enzyme catalytic sites rapidly, which hampers investigations of their trajectories. The high-resolution structure of the native Hordeum exo-hydrolase HvExoI isolated from seedlings reveals that non-covalently trapped glucose forms a stable enzyme-product complex. Here, we report that the alkyl ß-D-glucoside and methyl 6-thio-ß-gentiobioside substrate analogues perfused in crystalline HvExoI bind across the catalytic site after they displace glucose, while methyl 2-thio-ß-sophoroside attaches nearby. Structural analyses and multi-scale molecular modelling of nanoscale reactant movements in HvExoI reveal that upon productive binding of incoming substrates, the glucose product modifies its binding patterns and evokes the formation of a transient lateral cavity, which serves as a conduit for glucose departure to allow for the next catalytic round. This path enables substrate-product assisted processive catalysis through multiple hydrolytic events without HvExoI losing contact with oligo- or polymeric substrates. We anticipate that such enzyme plasticity could be prevalent among exo-hydrolases.

Comparing Hydrolysis and Transglycosylation Reactions Catalyzed by Thermus thermophilus ß-Glycosidase. A Combined MD and QM/MM Study.

The synthesis of oligosaccharides and other carbohydrate derivatives is of relevance for the advancement of glycosciences both at the fundamental and applied level. For many years, glycosyl hydrolases (GHs) have been explored to catalyze the synthesis of glycosidic bonds. In particular, retaining GHs can catalyze a transglycosylation (T) reaction that competes with hydrolysis (H). This has been done either employing controlled conditions in wild type GHs or by engineering new mutants. The goal, which is to increase the T/H ratio, has been achieved with moderate success in several cases despite the fact that the molecular basis for T/H modulation are unclear. Here we have used QM(DFT)/MM calculations to compare the glycosylation, hydrolysis and transglycosylation steps catalyzed by wild type Thermus thermophilus ß-glycosidase (family GH1), a retaining glycosyl hydrolase for which a transglycosylation yield of 36% has been determined experimentally. The three transition states have a strong oxocarbenium character and ring conformations between 4H3 and 4E. The atomic charges at the transition states for hydrolysis and transglycosylation are very similar, except for the more negative charge of the oxygen atom of water when compared to that of the acceptor Glc. The glycosylation transition state has a stronger S N 2 character than the deglycosylation ones and the proton transfer is less advanced. At the QM(PBE0/TZVP)/MM level, the TS for transglycosylation has shorter O4GLC-C1FUC (forming bond) distance and longer OE2GLU338-C1FUC (breaking) distance than the hydrolysis one, although the HACC proton is closer to the Glu164 base in the hydrolysis TS. The QM(SCC-DFTB)/MM free energy maxima show the inverted situation, although the hydrolysis TS presents significant structural fluctuations. The 3-OHGLC group of the acceptor Glc (transglycosylation) and WAT432 (neighbor water in hydrolysis) are identified to stabilize the oxocarbenium transition states through interaction with O5FUC and O4FUC. The analysis of interaction suggests that perturbing the Glu392-Fuc interaction could increase the T/H ratio, either by direct mutation of this residue or indirectly as reported experimentally in the Asn390I and Phe401S cases. The molecular understanding of similarities and differences between hydrolysis and transglycosylation steps may be of help in the design of new biocatalysts for glycan synthesis.

The sp2-iminosugar glycolipid 1-dodecylsulfonyl-5N,6O-oxomethylidenenojirimycin (DSO2-ONJ) as selective anti-inflammatory agent by modulation of hemeoxygenase-1 in Bv.2 microglial cells and retinal explants.

Neuroinflammation is an early event during diabetic retinopathy (DR) that impacts the dynamics of microglia polarization. Gliosis is a hallmark of DR and we have reported the beneficial effects of 1R-DSO-ONJ, a member of the sp2-iminosugar glycolipid (sp2-IGL) family, in targeting microglia and reducing gliosis in diabetic db/db mice. Herein, we analyzed the effect of DSO2-ONJ, another family compound incorporating a sulfone group that better mimics the phosphate group of phosphatidylinositol ether lipid analogues (PIAs), in Bv.2 microglial cells treated with bacterial lipopolysaccaride (LPS) and in retinal explants from db/db mice. In addition to decreasing iNOS and inflammasome activation, the anti-inflammatory effect of DSO2-ONJ was mediated by direct p38α MAPK activation. Computational docking experiments demonstrated that DSO2-ONJ binds to p38α MAPK at the same site where PIAs and the alkyl phospholipid perifosine activators do, suggesting similar mechanism of action. Moreover, treatment of microglial cells with DSO2-ONJ increased both heme-oxygenase (HO)-1 and Il10 expression regardless the presence of LPS. In retinal explants from db/db mice, DSO2-ONJ also induced HO-1 and reduced gliosis. Since IL-10-mediated induction of HO-1 expression is mediated by p38α MAPK activation, our results suggest that this molecular mechanism is involved in the anti-inflammatory effects of DSO2-ONJ in microglia.

Ab Initio QM/MM Modeling of the Rate-Limiting Proton Transfer Step in the Deamination of Tryptamine by Aromatic Amine Dehydrogenase.

Aromatic amine dehydrogenase (AADH) and related enzymes are at the heart of debates on the roles of quantum tunneling and protein dynamics in catalysis. The reaction of tryptamine in AADH involves significant quantum tunneling in the rate-limiting proton transfer step, shown by large H/D primary kinetic isotope effects (KIEs), with unusual temperature dependence. We apply correlated ab initio combined quantum mechanics/molecular mechanics (QM/MM) methods, at levels up to local coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), to calculate accurate potential energy surfaces for this reaction, which are necessary for quantitative analysis of tunneling contributions and reaction dynamics. Different levels of QM/MM treatment are tested. Multiple pathways are calculated with fully flexible transition state optimization by the climbing-image nudged elastic band method at the density functional QM/MM level. The average LCCSD(T) potential energy barriers to proton transfer are 16.7 and 14.0 kcal/mol for proton transfer to the two carboxylate atoms of the catalytic base, Asp128ß. The results show that two similar, but distinct pathways are energetically accessible. These two pathways have different barriers, exothermicity and curvature, and should be considered in analyses of the temperature dependence of reaction and KIEs in AADH and other enzymes. These results provide a benchmark for this prototypical enzyme reaction and will be useful for developing empirical models, and analyzing experimental data, to distinguish between different conceptual models of enzyme catalysis.

Computational insights into active site shaping for substrate specificity and reaction regioselectivity in the EXTL2 retaining glycosyltransferase.

Glycosyltransferases are enzymes that catalyze a monosaccharide transfer reaction from a donor to an acceptor substrate with the synthesis of a new glycosidic bond. They are highly substrate specific and regioselective, even though the acceptor substrate often presents multiple reactive groups. Currently, many efforts are dedicated to the development of biocatalysts for glycan synthesis and, therefore, a better understanding of how natural enzymes achieve this goal can be of valuable help. To gain a deeper insight into the catalytic strategies used by retaining glycosyltransferases, the wild type EXTL2 (CAZy family GT64) and four mutant forms (at positions 293 and 246) were studied using QM(DFT)/MM calculations and molecular dynamics simulations. Existing hypotheses on the roles of Arg293, an enigmatic residue in the CAZy family GT64 that seemed to contradict a mechanism through an oxocarbenium intermediate, and of Asp246 have been tested. We also provide a molecular interpretation for the results of site-directed mutagenesis experiments. Moreover, we have investigated why an Asp, and not a Glu like in the family GT6, is found on the ß-face of the transferred GlcNAc. It is predicted that an Asp246Glu mutant of EXTL2 would be unable to catalyze the α-1,4 transfer. The results herein presented clarify the roles that Arg293, Asp246 and Leu213 have at different stages of the catalytic process (for binding but also for efficient chemical reaction). Altogether, we provide a molecular view that connects the identity and conformation of these residues to the substrate specificity and regioselectivity of the enzyme, illustrating a delicate interplay between all these aspects.

Computational insight into the catalytic implication of head/tail-first orientation of arachidonic acid in human 5-lipoxygenase: consequences for the positional specificity of oxygenation.

In the present work we have combined homology modeling, protein-ligand dockings, quantum mechanics/molecular mechanics calculations and molecular dynamics simulations to generate human 5-lipoxygenase (5-LOX):arachidonic acid (AA) complexes consistent with the 5-lipoxygenating activity (which implies hydrogen abstraction at the C7 position). Our results suggest that both the holo and the apo forms of human Stable 5-LOX could accommodate AA in a productive form for 5-lipoxygenation. The former, in a tail-first orientation, with the AA carboxylate end interacting with Lys409, gives the desired structures with C7 close to the Fe-OH(-) cofactor and suitable barrier heights for H7 abstraction. Only when using the apo form structure, a head-first orientation with the AA carboxylate close to His600 (a residue recently proposed as essential for AA positioning) is obtained in the docking calculations. However, the calculated barrier heights for this head-first orientation are in principle consistent with 5-LOX specificity, but also with 12/8 regioselectivity. Finally, long MD simulations give support to the recent hypothesis that the Phe177 + Tyr181 pair needs to close the active site access during the chemical reaction, and suggest that in the case of a head-first orientation Phe177 may be the residue interacting with the AA carboxylate.

Evolutionary alteration of ALOX15 specificity optimizes the biosynthesis of antiinflammatory and proresolving lipoxins.

ALOX15 (12/15-lipoxygenase) orthologs have been implicated in maturational degradation of intracellular organelles and in the biosynthesis of antiinflammatory and proresolving eicosanoids. Here we hypothesized that lower mammals (mice, rats, pigs) express 12-lipoxygenating ALOX15 orthologs. In contrast, 15-lipoxygenating isoforms are found in higher primates (orangutans, men), and these results suggest an evolution of ALOX15 specificity. To test this hypothesis we first cloned and characterized ALOX15 orthologs of selected Catarrhini representing different stages of late primate evolution and found that higher primates (men, chimpanzees) express 15-lipoxygenating orthologs. In contrast, lower primates (baboons, rhesus monkeys) express 12-lipoxygenating enzymes. Gibbons, which are flanked in evolution by rhesus monkeys (12-lipoxygenating ALOX15) and orangutans (15-lipoxygenating ALOX15), express an ALOX15 ortholog with pronounced dual specificity. To explore the driving force for this evolutionary alterations, we quantified the lipoxin synthase activity of 12-lipoxygenating (rhesus monkey, mouse, rat, pig, humIle418Ala) and 15-lipoxygenating (man, chimpanzee, orangutan, rabbit, ratLeu353Phe) ALOX15 variants and found that, when normalized to their arachidonic acid oxygenase activities, the lipoxin synthase activities of 15-lipoxygenating ALOX15 variants were more than fivefold higher (P < 0.01) [corrected]. Comparative molecular dynamics simulations and quantum mechanics/molecular mechanics calculations indicated that, for the 15-lipoxygenating rabbit ALOX15, the energy barrier for C13-hydrogen abstraction (15-lipoxygenation) was 17 kJ/mol lower than for arachidonic acid 12-lipoxygenation. In contrast, for the 12-lipoxygenating Ile418Ala mutant, the energy barrier for 15-lipoxygenation was 10 kJ/mol higher than for 12-lipoxygenation. Taken together, our data suggest an evolution of ALOX15 specificity, which is aimed at optimizing the biosynthetic capacity for antiinflammatory and proresolving lipoxins.

Is Regioselectivity in the Enzyme-Catalyzed Hydroperoxidation of Arachidonic Acid Necessarily Determined by Hydrogen Abstraction? The Case of Rabbit Leu597Ala/Ile663Ala ALOX15 Mutant.

Molecular dynamics simulations and quantum mechanics/molecular mechanics calculations were performed on the in silico Leu597Ala/Ile663Ala double mutant of rabbit ALOX15 (12/15 lipoxygenase). The computational results suggested that subtle steric hindrance by the conserved Leu597 and C-terminal Ile663 residues disturbed H10 abstractions in wildtype ALOX15 (which abstracts H13), but if these two bulky residues were mutated to smaller ones, H10 abstraction was no longer impeded and the regioselectivity of the initial H-abstraction step was changed. However, site-directed mutagenesis with HPLC analysis of the products of the whole oxidation process showed that the regioselectivity of the hydroperoxidation was not altered. This disagreement may be explained by the conformational reorganization of the system needed to rotate the -OO. group from an antarafacial to a suprafacial arrangement prior to back-hydrogen transfer. After H10 abstraction and O2 insertion, the evolution of the peroxy radical at C12 was sterically impeded, whereas peroxyl group rotation at C15 (after H13 abstraction) could easily evolve to a suprafacial arrangement, which thus led to the final product. For this reason, the global regiospecificity was not affected in the mutant. These findings exemplify that the regioselectivity of initial hydrogen abstraction and the regioselectivity of the final product do not necessarily coincide (in fact, they can be opposite) for the hydroperoxidation of arachidonic acid catalyzed by a lipoxygenase.

Understanding the Mechanism of the Hydrogen Abstraction from Arachidonic Acid Catalyzed by the Human Enzyme 15-Lipoxygenase-2. A Quantum Mechanics/Molecular Mechanics Free Energy Simulation.

Lipoxygenases (LOXs) are a family of enzymes involved in the biosynthesis of several lipid mediators. In the case of human 15-LOX, the 15-LOX-1 and 15-LOX-2 isoforms show slightly different reaction regiospecificity and substrate specificity, indicating that substrate binding and recognition may be different, a fact that could be related to their different biological role. Here, we have used long molecular dynamics simulations, QM(DFT)/MM potential energy and free energy calculations (using the newly developed DHAM method), to investigate the binding mode of the arachidonic acid (AA) substrate into 15-LOX-2 and the rate-limiting hydrogen-abstraction reaction 15-LOX-2 catalyzes. Our results strongly indicate that hydrogen abstraction from C13 in 15-LOX-2 is only consistent with the "tail-first" orientation of AA, with its carboxylate group interacting with Arg429, and that only the pro-S H13 hydrogen will be abstracted (being the pro-R H13 and H10 too far from the acceptor oxygen atom). At the B3LYP/6-31G(d) level the potential and free energy barriers for the pro-S H13 abstraction of AA by 15-LOX-2 are 18.0 and 18.6 kcal/mol, respectively. To analyze the kinetics of the hydrogen abstraction process, we determined a Markov model corresponding to the unbiased simulations along the state-discretized reaction coordinate. The calculated rates based on the second largest eigenvalue of the Markov matrices agree well with experimental measurements, and also provide the means to directly determine the pre-exponential factor for the reaction by comparing with the free energy barrier height. Our calculated pre-exponential factor is close to the value of kBT/h. On the other hand, our results suggest that the spin inversion of the complete system (including the O2 molecule) that is required to happen at some point along the full process to lead to the final hydroperoxide product, is likely to take place during the hydrogen transfer, which is a proton coupled electron transfer. Overall, a different binding mode from the one accepted for 15-LOX-1 is proposed, which provides a molecular basis for 15-LOX-2 exclusive 15-HPETE production in front of the double (although highly 15-) 12/15 regiospecificity of 15-LOX-1. Understanding how these different isoenzymes achieve their regiospecificity is expected to help in specific inhibitor design.

How Can Linoleic Acid Be the Preferential Substrate of the Enzyme 15-Lipoxygenase-1? A QM/MM Approach.

The most common substrate of mammalian lipoxygenases (LOXs) is arachidonic acid (AA). However, 15-LOXs can present dual substrate specificity. These LOXs catalyze the peroxidation of AA, initiated by a H-abstraction step (mainly H13-abstraction) by the Fe(III)-OH(-) cofactor, and the peroxidation of linoleic acid (LA) after H11-abstraction. In this paper, QM(B3LYP)/MM(CHARMM) calculations of the rate-limiting H11-abstraction process of LA catalyzed by rabbit 15-LOX-1 (15-rLOX-1) have been carried out using a complete model of the solvated 15-rLOX-1:LA complex. A total of 26 QM/MM potential energy profiles as a function of the H-transfer reaction coordinate have been computed along with one QM/MM free energy profile obtained using the Free Energy Perturbation method. The molecular origin of substrate specificity of 15-rLOX-1 for LA in comparison with AA has been analyzed. In many of the QM/MM reactive 15-rLOX-1:LA energy minima, LA adopts more elongated conformations than AA, although having a shorter carbon chain, because LA has one double bond between C1 and C11 whereas AA has three double bonds between C1 and C13. Consequently, C11 of LA can be located in the same region of the active site as C13 of AA, a zone where H11-abstraction from LA as well as H13-abstraction from AA is not hindered by bulky residue side chains. This explains at a molecular level how 15-LOXs might accommodate and recognize for catalysis two substrates that are different in length by two carbons. Our results also explain why (9Z,11E)-13-hydro(pero)xyoctadeca-9,11-dienoic acid is the major product of the peroxidation and why LA is the preferential substrate of 15-rLOX-1.