Role of Arg403 for thermostability and catalytic activity of rabbit 12/15-lipoxygenase

Development of novel ALOX15 inhibitors combining dual machine learning filtering and fragment substitution optimisation approaches, molecular docking and dynamic simulation methods

The oxidation of unsaturated lipids, facilitated by the enzyme Arachidonic acid 15-lipoxygenase (ALOX15), is an essential element in the development of ferroptosis. This study combined a dual-score exclusion strategy with high-throughput virtual screening, naive Bayesian and recursive partitioning machine learning models, the already established ALOX15 inhibitor i472, and a docking-based fragment substitution optimisation approach to identify potential ALOX15 inhibitors, ultimately leading to the discovery of three FDA-approved drugs that demonstrate optimal inhibitory potential against ALOX15. Through fragment substitution-based optimisation, seven new inhibitor structures have been developed. To evaluate their practicality, ADMET predictions and molecular dynamics simulations were performed. In conclusion, the compounds found in this study provide a novel approach to combat conditions related to ferroptosis-related injury by inhibiting ALOX15.

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.

Computational Insights and In Silico Characterization of a Novel Mini-Lipoxygenase from Nostoc Sphaeroides and Its Application in the Quality Improvement of Steamed Bread

Lipoxygenase (EC1.13.11.12, LOX) has been potentially used in the food industry for food quality improvement. However, the low activity, poor thermal stability, narrow range of pH stability, as well as undesirable isoenzymes and off-flavors, have hampered the application of current commercial LOX. In this study, a putative mini-lipoxygenase gene from cyanobacteria, Nostoc sphaeroides (NsLOX), was cloned and expressed in E. coli BL21. NsLOX displayed only 26.62% structural identity with the reported LOX from Cyanothece sp., indicating it as a novel LOX. The purified NsLOX showed the maximum activity at pH 8.0 and 15 °C, with superior stability at a pH range from 6.0 to 13.0, retaining about 40% activity at 40 °C for 90 min. Notably, NsLOX exhibited the highest specific activity of 78,080 U/mg towards linoleic acid (LA), and the kinetic parameters-K_m, k_cat, and k_cat/K_m-attain values of 19.46 μM, 9199.75 s^-1, and 473.85 μM^-1 s^-1, respectively. Moreover, the activity of NsLOX was obviously activated by Ca²⁺, but it was completely inhibited by Zn²⁺ and Cu²⁺. Finally, NsLOX was supplied in steamed bread and contributed even better improved bread quality than the commercial LOX. These results suggest NsLOX as a promising substitute of current commercial LOX for application in the food industry.

One for All, All for One: The Peculiar Dynamics of TNF-Receptor-Associated Factor (TRAF2) Subunits

TNF Receptor-Associated Factor 2 (TRAF2) is a homo-trimer belonging to the TNF-receptor-associated factor family (TRAFs). The TRAF2 oligomeric state is crucial for receptor binding, the interaction with other proteins (involved in the TNFR signaling), and the interaction with biological membranes. In this study, we present a computational analysis of the Molecular Dynamics of TRAF2-C (a truncated and soluble TRAF2 form) to identify patterns in the interactions between the three chains. We have performed a canonical analysis of the motion applied to molecular dynamics starting from the available crystal structure to identify correlated motions in TRAF2 dynamics. We have computed the displacement matrix, providing a frame-by-frame displacement for each residue in the dynamic. We provide the results in terms of the correlation matrix, which represents a detailed map of the correlated motions of residues. Eventually, we computed the so-called dynamical clusters, based on the Principal Component Analysis (PCA) of the motion (displacement) and the k means application on the first two principal components space. The results clearly indicate that, most of the time, two chains move in a strongly correlated motion, while the third chain follows a freer motion. A detailed analysis of the correlation matrix also shows that a few specific interface residues characterize the interaction of the more independent subunit with the other two. These findings suggest that the equilibrium between the trimer and the dissociated species (dimers and monomers) might be finely tuned by controlling a few critical residues in the protein quaternary structure, probably facilitating the regulation of oligomerization and dissociation in vivo.

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.

A role of Gln596 in fine-tuning mammalian ALOX15 specificity, protein stability and allosteric properties

His596 of human ALOX12 has been suggested to interact with the COO^--group of arachidonic acid during ALOX catalysis. In mammalian ALOX15 orthologs Gln596 occupies this position and this amino acid exchange might contribute to the functional differences between the two ALOX-isoforms. To explore the role of Gln596 for ALOX15 functionality we mutated this amino acid to different residues in rabbit and human ALOX15 and investigated the impact of these mutations on structural, catalytic and allosteric enzyme properties. To shed light on the molecular basis of the observed functional alterations we performed in silico substrate docking studies and molecular dynamics simulations and also explored the impact of Gln596 exchange on the protein structure. The combined theoretical and experimental data suggest that Gln596 may not directly interact with the COO^--group of arachidonic acid. In contrast, mutations at Gln596 destabilize the secondary and tertiary structure of ALOX15 orthologs, which may be related to a disturbance of the electrostatic interaction network with other amino acids in the immediate surrounding. Moreover, our MD-simulations suggest that the geometry of the dimer interface depends on the structure of substrate bound inside the substrate-binding pocket and that Gln596Ala exchange impairs the allosteric properties of the enzyme. Taken together, these data indicate the structural and functional importance of Gln596 for ALOX15 catalysis.

Gly188Arg substitution eliminates substrate inhibition in arachidonate 11R-lipoxygenase

Lipoxygenases (LOXs) are dioxygenases that catalyze the oxygenation of polyunsaturated fatty acids to hydroperoxyl derivates. These products are precursors for different lipid mediators which are associated with pathogenesis of various diseases such as asthma, atherosclerosis and cancer. Several LOXs suffer from substrate inhibition, a potential regulatory mechanism, yet it is unclear what is the cause of this phenomenon. One such enzyme is the coral 11R-LOX which displays a significant decrease in turnover rate at arachidonic acid concentrations above 30 μM. In this report, site-directed mutagenesis and inhibition assays were employed to shed light on the mechanism of substrate inhibition in 11R-LOX. We found that introduction of a positive charge to the active site entrance with Gly188Arg substitution completely eliminates the slow-down at higher substrate concentrations. Inhibition of 11R-LOX by its catalysis product, 11(R)-hydroperoxyeicosatetraenoic acid, suggests an uncompetitive mechanism. We reason that substrate inhibition in 11R-LOX is due to additional fatty acid binding by the enzyme:substrate complex at an allosteric site situated in the very vicinity of the active site entrance.

Probing the Electrostatic and Steric Requirements for Substrate Binding in Human Platelet-Type 12-Lipoxygenase

Human platelet ALOX12 (hALOX12 or h12-LOX) has been implicated in a variety of human diseases. The present study investigates the active site of hALOX12 to more thoroughly understand how it positions the substrate and achieves nearly perfect regio- and stereospecificities (i.e., 100 ± 5% of the 12(S)-hydroperoxide product), utilizing site-directed mutagenesis. Specifically, we have determined that Arg402 is not as important in substrate binding as previously seen for hALOX15 but that His596 may play a role in anchoring the carboxy terminal of the arachidonic acid during catalysis. In addition, Phe414 creates a π-stacking interaction with a double bond of arachidonic acid (Δ¹¹), and Ala417/Val418 define the bottom of the cavity. However, the influence of Ala417/Val418 on the profile is markedly less for hALOX12 than that seen in hALOX15. Mutating these two residues to larger amino acids (Ala417Ile/Val418Met) only increased the generation of 15-HpETE by 24 ± 2%, but conversely, smaller residues at these positions converted hALOX15 to almost 100% hALOX12 reactivity [Gan et al. (1996) J. Biol. Chem. 271, 25412-25418]. However, we were able to increase 15-HpETE to 46 ± 3% by restricting the width of the active site with the Ala417Ile/Val418Met/Ser594Thr mutation, indicating both depth and width of the active site are important. Finally, residue Leu407 is shown to play a critical role in positioning the substrate correctly, as seen by the increase of 15-HpETE to 21 ± 1% for the single Leu407Gly mutant. These results outline critical differences between the active site requirements of hALOX12 relative to hALOX15 and explain both their product specificity and inhibitory differences.

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.

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.

Structural and functional biology of arachidonic acid 15-lipoxygenase-1 (ALOX15)

Lipoxygenases (LOX) form a family of lipid peroxidizing enzymes, which have been implicated in a number of physiological processes and in the pathogenesis of inflammatory, hyperproliferative and neurodegenerative diseases. They occur in two of the three domains of terrestrial life (bacteria, eucarya) and the human genome involves six functional LOX genes, which encode for six different LOX isoforms. One of these isoforms is ALOX15, which has first been described in rabbits in 1974 as enzyme capable of oxidizing membrane phospholipids during the maturational breakdown of mitochondria in immature red blood cells. During the following decades ALOX15 has extensively been characterized and its biological functions have been studied in a number of cellular in vitro systems as well as in various whole animal disease models. This review is aimed at summarizing the current knowledge on the protein-chemical, molecular biological and enzymatic properties of ALOX15 in various species (human, mouse, rabbit, rat) as well as its implication in cellular physiology and in the pathogenesis of various diseases.

Molecular basis for the catalytic inactivity of a naturally occurring near-null variant of human ALOX15

Mammalian lipoxygenases belong to a family of lipid-peroxidizing enzymes, which have been implicated in cardiovascular, hyperproliferative and neurodegenerative diseases. Here we report that a naturally occurring mutation in the hALOX15 gene leads to expression of a catalytically near-null enzyme variant (hGly422Glu). The inactivity may be related to severe misfolding of the enzyme protein, which was concluded from CD-spectra as well as from thermal and chemical stability assays. In silico mutagenesis experiments suggest that most mutations at hGly422 have the potential to induce sterical clash, which might be considered a reason for protein misfolding. hGly422 is conserved among ALOX5, ALOX12 and ALOX15 isoforms and corresponding hALOX12 and hALOX5 mutants also exhibited a reduced catalytic activity. Interestingly, in the hALOX5 Gly429Glu mutants the reaction specificity of arachidonic acid oxygenation was shifted from 5S- to 8S- and 12R-H(p)ETE formation. Taken together, our data indicate that the conserved glycine is of functional importance for these enzyme variants and most mutants at this position lose catalytic activity.