Unraveling How Enzymes Can Use Bulky Residues To Drive Site-Selective C–H Activation: The Case of Mammalian Lipoxygenases Catalyzing Arachidonic Acid Oxidation

Bony Fish Arachidonic Acid 15-Lipoxygenases Exhibit Different Catalytic Properties than Their Mammalian Orthologs, Suggesting Functional Enzyme Evolution during Vertebrate Development

The human genome involves six functional arachidonic acid lipoxygenase (ALOX) genes and the corresponding enzymes (ALOX15, ALOX15B, ALOX12, ALOX12B, ALOXE3, ALOX5) have been implicated in cell differentiation and in the pathogenesis of inflammatory, hyperproliferative, metabolic, and neurological disorders. In other vertebrates, ALOX-isoforms have also been identified, but they occur less frequently. Since bony fish represent the most abundant subclass of vertebrates, we recently expressed and characterized putative ALOX15 orthologs of three different bony fish species (Nothobranchius furzeri, Pundamilia nyererei, Scleropages formosus). To explore whether these enzymes represent functional equivalents of mammalian ALOX15 orthologs, we here compared a number of structural and functional characteristics of these ALOX-isoforms with those of mammalian enzymes. We found that in contrast to mammalian ALOX15 orthologs, which exhibit a broad substrate specificity, a membrane oxygenase activity, and a special type of dual reaction specificity, the putative bony fish ALOX15 orthologs strongly prefer C20 fatty acids, lack any membrane oxygenase activity and exhibit a different type of dual reaction specificity with arachidonic acid. Moreover, mutagenesis studies indicated that the Triad Concept, which explains the reaction specificity of all mammalian ALOX15 orthologs, is not applicable for the putative bony fish enzymes. The observed functional differences between putative bony fish ALOX15 orthologs and corresponding mammalian enzymes suggest a targeted optimization of the catalytic properties of ALOX15 orthologs during vertebrate development.

Different Structures-Similar Effect: Do Substituted 5-(4-Methoxyphenyl)-1H-indoles and 5-(4-Methoxyphenyl)-1H-imidazoles Represent a Common Pharmacophore for Substrate Selective Inhibition of Linoleate Oxygenase Activity of ALOX15?

Mammalian 15-lipoxygenases (ALOX15) are lipid peroxidizing enzymes that exhibit variable functionality in different cancer and inflammation models. The pathophysiological role of linoleic acid- and arachidonic acid-derived ALOX15 metabolites rendered this enzyme a target for pharmacological research. Several indole and imidazole derivatives inhibit the catalytic activity of rabbit ALOX15 in a substrate-specific manner, but the molecular basis for this allosteric inhibition remains unclear. Here, we attempt to define a common pharmacophore, which is critical for this allosteric inhibition. We found that substituted imidazoles induce weaker inhibitory effects when compared with the indole derivatives. In silico docking studies and molecular dynamics simulations using a dimeric allosteric enzyme model, in which the inhibitor occupies the substrate-binding pocket of one monomer, whereas the substrate fatty acid is bound at the catalytic center of another monomer within the ALOX15 dimer, indicated that chemical modification of the core pharmacophore alters the enzyme-inhibitor interactions, inducing a reduced inhibitory potency. In our dimeric ALOX15 model, the structural differences induced by inhibitor binding are translated to the hydrophobic dimerization cluster and affect the structures of enzyme-substrate complexes. These data are of particular importance since substrate-specific inhibition may contribute to elucidation of the putative roles of ALOX15 metabolites derived from different polyunsaturated fatty acids in mammalian pathophysiology.

Functional Characterization of Mouse and Human Arachidonic Acid Lipoxygenase 15B (ALOX15B) Orthologs and of Their Mutants Exhibiting Humanized and Murinized Reaction Specificities

The arachidonic acid lipoxygenase 15B (ALOX15B) orthologs of men and mice form different reaction products when arachidonic acid is used as the substrate. Tyr603Asp+His604Val double mutation in mouse arachidonic acid lipoxygenase 15b humanized the product pattern and an inverse mutagenesis strategy murinized the specificity of the human enzyme. As the mechanistic basis for these functional differences, an inverse substrate binding at the active site of the enzymes has been suggested, but experimental proof for this hypothesis is still pending. Here we expressed wildtype mouse and human arachidonic acid lipoxygenase 15B orthologs as well as their humanized and murinized double mutants as recombinant proteins and analyzed the product patterns of these enzymes with different polyenoic fatty acids. In addition, in silico substrate docking studies and molecular dynamics simulation were performed to explore the mechanistic basis for the distinct reaction specificities of the different enzyme variants. Wildtype human arachidonic acid lipoxygenase 15B converted arachidonic acid and eicosapentaenoic acid to their 15-hydroperoxy derivatives but the Asp602Tyr+Val603His exchange murinized the product pattern. The inverse mutagenesis strategy in mouse arachidonic acid lipoxygenase 15b (Tyr603Asp+His604Val exchange) humanized the product pattern with these substrates, but the situation was different with docosahexaenoic acid. Here, Tyr603Asp+His604Val substitution in mouse arachidonic acid lipoxygenase 15b also humanized the specificity but the inverse mutagenesis (Asp602Tyr+Val603His) did not murinize the human enzyme. With linoleic acid Tyr603Asp+His604Val substitution in mouse arachidonic acid lipoxygenase 15b humanized the product pattern but the inverse mutagenesis in human arachidonic acid lipoxygenase 15B induced racemic product formation. Amino acid exchanges at critical positions of human and mouse arachidonic acid lipoxygenase 15B orthologs humanized/murinized the product pattern with C20 fatty acids, but this was not the case with fatty acid substrates of different chain lengths. Asp602Tyr+Val603His exchange murinized the product pattern of human arachidonic acid lipoxygenase 15B with arachidonic acid, eicosapentaenoic acid, and docosahexaenoic acid. An inverse mutagenesis strategy on mouse arachidonic acid lipoxygenase 15b (Tyr603Asp+His604Val exchange) did humanize the reaction products with arachidonic acid and eicosapentaenoic acid, but not with docosahexaenoic acid.

Hydroperoxidation of Docosahexaenoic Acid by Human ALOX12 and pigALOX15-mini-LOX

Human lipoxygenase 12 (hALOX12) catalyzes the conversion of docosahexaenoic acid (DHA) into mainly 14S-hydroperoxy-4Z,7Z,10Z,12E,16Z,19Z-docosahexaenoic acid (14S-H(p)DHA). This hydroperoxidation reaction is followed by an epoxidation and hydrolysis process that finally leads to maresin 1 (MaR1), a potent bioactive specialized pro-resolving mediator (SPM) in chronic inflammation resolution. By combining docking, molecular dynamics simulations, and quantum mechanics/molecular mechanics calculations, we have computed the potential energy profile of DHA hydroperoxidation in the active site of hALOX12. Our results describe the structural evolution of the molecular system at each step of this catalytic reaction pathway. Noteworthy, the required stereospecificity of the reaction leading to MaR1 is explained by the configurations adopted by DHA bound to hALOX12, along with the stereochemistry of the pentadienyl radical formed after the first step of the mechanism. In pig lipoxygenase 15 (pigALOX15-mini-LOX), our calculations suggest that 14S-H(p)DHA can be formed, but with a stereochemistry that is inadequate for MaR1 biosynthesis.

C-N Bond Forming Radical Rebound Is the Enantioselectivity-Determining Step in P411-Catalyzed Enantioselective C(sp3)-H Amination: A Combined Computational and Experimental Investigation

Engineered metalloenzymes represent promising catalysts for stereoselective C-H functionalization reactions. Recently, P450 enzymes have been evolved to allow for new-to-nature intramolecular C(sp3)-H amination reactions via a nitrene transfer mechanism, giving rise to diamine derivatives with excellent enantiocontrol. To shed light on the origin of enantioselectivity, a combined computational and experimental study was carried out. Hybrid quantum mechanics/molecular mechanics calculations were performed to investigate the activation energies and enantioselectivities of both the hydrogen atom transfer (HAT) and the subsequent C-N bond forming radical rebound steps. Contrary to previously hypothesized enantioinduction mechanisms, our calculations show that the radical rebound step is enantioselectivity-determining, whereas the preceding HAT step is only moderately stereoselective. Furthermore, the selectivity in the initial HAT is ablated by rapid conformational change of the radical intermediate prior to C-N bond formation. This finding is corroborated by our experimental study using a set of enantiomerically pure, monodeuterated substrates. Furthermore, classical and ab initio molecular dynamics simulations were carried out to investigate the conformational flexibility of the carbon-centered radical intermediate. This key radical species undergoes a facile conformational change in the enzyme active site from the pro-(R) to the pro-(S) configuration, whereas the radical rebound is slower due to the spin-state change and ring strain of the cyclization process, thereby allowing stereoablative C-N bond formation. Together, these studies revealed an underappreciated enantioinduction mechanism in biocatalytic C(sp3)-H functionalizations involving radical intermediates, opening up new avenues for the development of other challenging asymmetric C(sp3)-H functionalizations.

The role of acetylated cyclooxygenase-2 in the biosynthesis of resolvin precursors derived from eicosapentaenoic acid

Specialized pro-resolving lipid mediators (SPMs) are natural bioactive agents actively involved in inflammation resolution. SPMs act when uncontrolled inflammatory processes are developed, for instance, in patients of COVID-19 or other diseases. The so-called resolution pharmacology aims at developing new treatments based on the use of SPMs as agonists, which promote inflammation resolution without unwanted side effects. It has been shown that the biosynthesis of SPMs called eicosapentaenoic acid (EPA)-derived E-series resolvins is initiated by aspirin-acetylated COX-2 from EPA, leading to 18-hydroperoxy-eicosapentaenoic acid (18-HpEPE). However, there are many open questions concerning the intriguing role of aspirin in the molecular mechanism of resolvin formation. Our MD simulations, combined with QM/MM calculations, show that the potential energy barriers for the H₁₆-abstraction from EPA, required for forming 18-HpEPE, are higher than for the H₁₃-abstraction, thus explaining why 18-HpEPE is a marginal product of COX-2 catalysis. By contrast, in the aspirin-acetylated COX-2/EPA complex, the H_16proS-abstraction energy barriers are somewhat lower than the H_13proS energy barriers and much smaller than the H₁₆-transfer barriers in the wild type COX-2/EPA system. Those results agree with the experimental observation that aspirin favours the synthesis of several SPMs known as aspirin-triggered resolvins. In the following step of the catalytic mechanism, the calculated O₂ addition to C₁₈ is preferred versus the addition to C₁₄ which also agrees with 18R-HEPE and 18S-HEPE being the main products from EPA in aspirin-acetylated COX-2.

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.

Conformational Heterogeneity and Cooperative Effects of Mammalian ALOX15

Arachidonic acid lipoxygenases (ALOXs) have been suggested to function as monomeric enzymes, but more recent data on rabbit ALOX15 indicated that there is a dynamic monomer-dimer equilibrium in aqueous solution. In the presence of an active site ligand (the ALOX15 inhibitor RS7) rabbit ALOX15 was crystalized as heterodimer and the X-ray coordinates of the two monomers within the dimer exhibit subtle structural differences. Using native polyacrylamide electrophoresis, we here observed that highly purified and predominantly monomeric rabbit ALOX15 and human ALOX15B are present in two conformers with distinct electrophoretic mobilities. In silico docking studies, molecular dynamics simulations, site directed mutagenesis experiments and kinetic measurements suggested that in aqueous solutions the two enzymes exhibit motional flexibility, which may impact the enzymatic properties.

Deciphering the Molecular Details of the Lipoxin Formation Mechanism in the 5(S),15(S)-DiHpETE Biosynthetic Pathway Catalyzed by Reticulocyte 15-Lipoxygenase-1

Chronic inflammation is now widely recognized to play important roles in many commonly occurring diseases, including COVID-19. The resolution response to this chronic inflammation is an active process governed by specialized pro-resolving mediators (SPMs) like the lipid mediators known as lipoxins. The biosynthesis of lipoxins is catalyzed by several lipoxygenases (LOXs) from arachidonic acid. However, the molecular details of the mechanisms involved are not well known yet. In this paper, we have combined molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations to analyze how reticulocyte 15-LOX-1 catalyzes the production of lipoxins from 5(S),15(S)-diHpETE. Our results indicate that the dehydration mechanism from 5(S),15(S)-diHpETE, via the formation of an epoxide, presents huge energy barriers even though it was one of the two a priori synthetic proposals. This result is compatible with the fact that no epoxide has been directly detected as an intermediate in the catalytic formation of lipoxins from 5(S),15(S)-diHpETE. Conversely, the oxygenation of 5(S),15(S)-diHpETE at C₁₄ is feasible because there is an open channel connecting the protein surface with this carbon atom, and the energy barrier for oxygen addition through this channel is small. The analysis of the following steps of this mechanism, leading to the corresponding hydroperoxide at the 15-LOX-1 active site, indicates that the oxygenation mechanism will lead to the formation of lipoxinB₄ after the final action of a reductase. In contrast, our calculations are in agreement with experiments that lipoxinA₄ cannot derive from 5(S),15(S)-diHpETE by either of the two proposed mechanisms and that 5(S),15(S)-diHETE is not an intermediate of lipoxin biosynthesis catalyzed by 15-LOX-1.

Quantum Chemical and QM/MM Models in Biochemistry

Quantum chemical (QC) calculations provide a basis for deriving a microscopic understanding of enzymes and photobiological systems. Here we describe how QC models can be used to explore the electronic structure, dynamics, and energetics of biomolecules. We introduce the hybrid quantum mechanics/classical mechanics (QM/MM) approach, where a quantum mechanically described system of interest is embedded in a classically described force field representation of the biochemical surroundings. We also discuss the QM cluster model approach, as well as embedding theories, that provide complementary methodologies to model quantum mechanical effects in biomolecules. The chapter also provides some practical guides for building quantum biochemical models using the quinone reduction catalysis in respiratory complex I and a model reaction in solution as examples.

A protocol to obtain multidimensional quantum tunneling corrections derived from QM(DFT)/MM calculations for an enzyme reaction

The multidimensional small-curvature tunneling (SCT) method with Electrostatic Embedding calculations is a compromise between an accessible computational cost and the attainment of an accurate enough estimation of tunneling for an enzyme reaction.

Mutations of Triad Determinants Changes the Substrate Alignment at the Catalytic Center of Human ALOX5

For the specificity of ALOX15 orthologs of different mammals, the geometry of the amino acids Phe353, Ile418, Met419, and Ile593 ("triad determinants") is important, and mutagenesis of these residues altered the reaction specificity of these enzymes. Here we expressed wild-type human ALOX5 and its F359W/A424I/N425M/A603I mutant in Sf9 insect cells and characterized the catalytic differences of the two enzyme variants. We found that wild-type ALOX5 converted arachidonic acid mainly to 5(S)-hydroperoxyeicosatetraenoic acid (HpETE). In contrast, 15(S)- and 8(S)-H(p)ETE were formed by the mutant enzyme. In addition to arachidonic acid, wild-type ALOX5 accepted eicosapentaenoic acid (EPA) as substrate, but C18 fatty acids were not oxygenated. The quadruple mutant also accepted linoleic acid and α- and γ-linolenic acid as substrate. Structural analysis of the oxygenation products and kinetic studies with stereospecifically labeled 11(S)- and 11(R)-deutero-linoleic acid suggested alternative ways of substrate orientation at the active site. In silico docking studies, molecular dynamics simulations, and quantum mechanics/molecular mechanics (QM/MM) calculations confirmed this hypothesis. These data indicate that "triad determinant" mutagenesis alters the catalytic properties of ALOX5 abolishing its leukotriene synthase activity but improving its biosynthetic capacity for pro-resolving lipoxins.

Insights into the Mechanism and Enantioselectivity in the Biosynthesis of Ergot Alkaloid Cycloclavine Catalyzed by Aj_EasH from Aspergillus japonicus

Cycloclavine is a complex ergot alkaloid containing an unusual cyclopropyl moiety, which has a wide range of biological activities and pharmaceutical applications. The biosynthesis of cycloclavine requires a series of enzymes, one of which is a nonheme Fe^II/α-ketoglutarate-dependent (aKG) oxidase (Aj_EasH). According to the previous proposal, the cyclopropyl ring formation catalyzed by Aj_EasH follows an unprecedented oxidative mechanism; however, the reaction details are unknown. In this article, on the basis of the recently obtained crystal structure of Aj_EasH (EasH from Aspergillus japonicas), the reactant models were built, and the reaction details were investigated by performing QM-only and combined QM and MM calculations. Our calculation results reveal that the biosynthesis of cyclopropyl moiety involves a radical intermediate rather than a carbocationic or carbanionic intermediate as in the biosynthesis of terpenoid family. The iron(IV)-oxo first abstracts a hydrogen atom from the substrate to trigger the reaction, and then the generated radical intermediate undergoes ring rearrangement to form the fused 5-3 ring system of cycloclavine. On the basis of our calculations, the absolute configuration of the cycloclavine catalyzed by Aj_EasH from Aspergillus japonicus should be (5R,8R,10R), which is different from the product isolated from Ipomoea hildebrandtii (5R,8S,10S). Residues at the active site play an important role in substrate binding, ring rearrangement, and enantioselectivity.

Computational Understanding of the Selectivities in Metalloenzymes

Metalloenzymes catalyze many different types of biological reactions with high efficiency and remarkable selectivity. The quantum chemical cluster approach and the combined quantum mechanics/molecular mechanics methods have proven very successful in the elucidation of the reaction mechanism and rationalization of selectivities in enzymes. In this review, recent progress in the computational understanding of various selectivities including chemoselectivity, regioselectivity, and stereoselectivity, in metalloenzymes, is discussed.

Computational evidence for the degradation mechanism of haloalkane dehalogenase LinB and mutants of Leu248 to 1-chlorobutane

The catalytic degradation ability of the haloalkane dehalogenase LinB toward 1-chlorobutane (1-CB) was studied using a combined quantum mechanics/molecular mechanics (QM/MM) approach. Two major processes are involved in the LinB-catalyzed removal of halogens: dechlorination and hydrolyzation. The present study confirmed the experimentally proposed reaction path at the molecular level. Moreover, based on nucleophilic substitution mechanism (SN2 reaction), dechlorination was found to be the rate-determining step of the entire reaction process. In this study, the Boltzmann-weighted average barrier for dechlorination was determined to be 17.0 kcal mol-1, which is fairly close to the experimental value (17.4 kcal mol-1). The state of His107 and the influence of Leu248 on the dechlorination process were also explored. In addition, an intriguing phenomenon was discovered: the potential energy barrier decreased by 7.5 kcal mol-1 when the Leu248 residue was mutated into Phe248. This discovery might be of great help to design new mutant enzymes or novel biocatalysts.

Enzymatic browning in avocado (Persea americana) revisited: History, advances, and future perspectives

Considering nearly 80 years of research regarding one of the enzymes responsible for catalyzing the formation of pigments in higher animals, plants, fungi and bacteria, this review will focus on collecting and categorizing the existing information about polyphenol oxidase (PPO) in fruits, with particular emphasis on the information in relation to avocado, which is one of the hardiest species in terms of inactivation, has documented dual activity (EC 1.14.18.1/EC 1.10.3.1), and represents one of the oldest challenges for food science research and fruit processors. It is expected that this review will contribute to the further development of the field by highlighting the questions that have arisen during the characterization of PPO, the progress that has been made and the questions that remain today, in addition to new methodologies that are being applied to study this system. Holistic methodologies offer unexplored potential for advancing our understanding of the complex phenomena that govern PPO activity in fruits, because these methodologies will enable the characterization of this family of enzymes in all of its complexity. Subsequently, it will be possible to develop better techniques for controlling enzymatic browning in this valuable fruit.

Mutagenesis of Sequence Determinants of Truncated Porcine ALOX15 Induces Changes in the Reaction Specificity by Altering the Catalytic Mechanism of Initial Hydrogen Abstraction

The reaction specificity of lipoxygenases is of physiological relevance since the various oxygenation products exhibit different biological activities. Among mammalian ALOX15 orthologs there are arachidonic acid 12- and 15-lipoxygenating enzymes and recent studies suggested an evolutionary switch in that reaction specificity during late primate development. Previous reports showed that 12-lipoxygenating ALOX15 orthologs can be converted to 15-lipoxygenating enzymes by site-directed mutagenesis of some sequence determinants. Unfortunately, the molecular basis for those alterations are not well understood. Here, the arachidonic acid 12-lipoxygenating N-terminal truncation variant of pig ALOX15, for which a crystal structure is available, was used to explore the catalytic mechanism of the specificity switch induced by mutagenesis of Val418 and Val419 sequence determinants. We found that Val418Ile+Val419Met double mutant is dominantly 15-lipoxygenating. Docking and MD simulations, and quantum mechanics/molecular mechanics calculations indicated that the wildtype energy barrier for arachidonic acid 15-lipoxygenation is 3.4 kcal mol^-1 higher than for 12-lipoxygenation. In contrast, for the Val418Ile+Val419Met double mutant the energy barrier for 12-lipoxygenation is 6.0 kcal mol^-1 higher than for 15-lipoxygenation. Our data suggest that enzyme-substrate complex geometries determine the value of these energy barriers and, as a consequence, the reaction specificity of ALOX15 orthologs.

Insights into the unprecedented epoxidation mechanism of fumitremorgin B endoperoxidase (FtmOx1) from Aspergillus fumigatus by QM/MM calculations

QM/MM calculations indicate that the quintet of the Fe^IVO complex firstly abstracts the hydrogen from Tyr228 to initiate the reaction, then the generated Tyr228 radical extracts the hydrogen from C₂₁ to form the C₂₁ radical, which binds the second dioxygen to complete the epoxidation.

How Many Conformations of Enzymes Should Be Sampled for DFT/MM Calculations? A Case Study of Fluoroacetate Dehalogenase

The quantum mechanics/molecular mechanics (QM/MM) method (e.g., density functional theory (DFT)/MM) is important in elucidating enzymatic mechanisms. It is indispensable to study "multiple" conformations of enzymes to get unbiased energetic and structural results. One challenging problem, however, is to determine the minimum number of conformations for DFT/MM calculations. Here, we propose two convergence criteria, namely the Boltzmann-weighted average barrier and the disproportionate effect, to tentatively address this issue. The criteria were tested by defluorination reaction catalyzed by fluoroacetate dehalogenase. The results suggest that at least 20 conformations of enzymatic residues are required for convergence using DFT/MM calculations. We also tested the correlation of energy barriers between small QM regions and big QM regions. A roughly positive correlation was found. This kind of correlation has not been reported in the literature. The correlation inspires us to propose a protocol for more efficient sampling. This saves 50% of the computational cost in our current case.

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.

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.

Catalytic mechanism of C–F bond cleavage: insights from QM/MM analysis of fluoroacetate dehalogenase

The catalytic mechanisms of fluoroacetate dehalogenase (FAcD) toward substrates fluoroacetate and chloroacetate were studied by a combined quantum mechanics/molecular mechanics (QM/MM) method.

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

Using a multi-scale approach to search for the arachidonic acid binding modes that determine the catalytic specificity of human 5-LOX.

Insight into the catalytic mechanism of meta-cleavage product hydrolase BphD: a quantum mechanics/molecular mechanics study

The catalytic mechanism of BphD (the fourth enzyme of the biphenyl catabolic pathway) toward its natural substrate 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoic acid (HOPDA) was investigated in atomistic detail by QM/MM approach.

Insights into the catalytic mechanism of chlorophenol 4-monooxygenase: a quantum mechanics/molecular mechanics study

The degradation mechanism of chlorophenol 4-monooxygenase toward pollutants 2,4,5-trichlorophenol, 2,4,6-trichlorophenol, and 2,5-dichloro-p-hydroquinone was studied by QM/MM investigations.