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.

Molecular Insights into the Regulation of 3-Phosphoinositide-Dependent Protein Kinase 1: Modeling the Interaction between the Kinase and the Pleckstrin Homology Domains.

The 3-phosphoinositide-dependent protein kinase 1 (PDK1) K465E mutant kinase can still activate protein kinase B (PKB) at the membrane in a phosphatidylinositol-3,4,5-trisphosphate (PIP3, PtdIns(3,4,5)P3) independent manner. To understand this new PDK1 regulatory mechanism, docking and molecular dynamics calculations were performed for the first time to simulate the wild-type kinase domain-pleckstrin homology (PH) domain complex with PH-in and PH-out conformations. These simulations were then compared to the PH-in model of the KD-PH(mutant K465E) PDK1 complex. Additionally, three KD-PH complexes were simulated, including a substrate analogue bound to a hydrophobic pocket (denominated the PIF-pocket) substrate-docking site. We find that only the PH-out conformation, with the PH domain well-oriented to interact with the cellular membrane, is active for wild-type PDK1. In contrast, the active conformation of the PDK1 K465E mutant is PH-in, being ATP-stable at the active site while the PIF-pocket is more accessible to the peptide substrate. We corroborate that both the docking-site binding and the catalytic activity are in fact enhanced in knock-in mouse samples expressing the PDK1 K465E protein, enabling the phosphorylation of PKB in the absence of PIP3 binding.

Theoretical Characterization of the Step-by-Step Mechanism of Conversion of Leukotriene A4 to Leukotriene B4 Catalysed by the Enzyme Leukotriene A4 Hydrolase.

LTA4H is a bifunctional zinc metalloenzyme that converts leukotriene A4 (LTA4) into leukotriene B4 (LTB4), one of the most potent chemotactic agents involved in acute and chronic inflammatory diseases. In this reaction, LTA4H acts as an epoxide hydrolase with a unique and fascinating mechanism, which includes the stereoselective attachment of one water molecule to the carbon backbone of LTA4 several methylene units away from the epoxide moiety. By combining Molecular Dynamics simulations and Quantum Mechanics/Molecular Mechanics calculations, we obtained a very detailed molecular picture of the different consecutive steps of that mechanism. By means of a rather unusual 1,7-nucleophilic substitution through a clear SN1 mechanism, the epoxide opens and the triene moiety of the substrate twists in such a way that the bond C6-C7 adopts its cis (Z) configuration, thus exposing the R face of C12 to the addition of a water molecule hydrogen-bonded to ASP375. Thus, the two stereochemical features that are required for the bioactivity of LTB4 appear to be closely related. The noncovalent π-π stacking interactions between the triene moiety and two tyrosines (TYR267 and, especially, TYR378) that wrap the triene system along the whole reaction explain the preference for the cis configuration inside LTA4H.

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 H16-abstraction from EPA, required for forming 18-HpEPE, are higher than for the H13-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 H16proS-abstraction energy barriers are somewhat lower than the H13proS energy barriers and much smaller than the H16-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 O2 addition to C18 is preferred versus the addition to C14 which also agrees with 18R-HEPE and 18S-HEPE being the main products from EPA in aspirin-acetylated COX-2.

N-Substituted 5-(1H-Indol-2-yl)-2-methoxyanilines Are Allosteric Inhibitors of the Linoleate Oxygenase Activity of Selected Mammalian ALOX15 Orthologs: Mechanism of Action.

Here, we describe the first systematic study on the mechanism of substrate-selective inhibition of mammalian ALOX15 orthologs. For this purpose, we prepared a series of N-substituted 5-(1H-indol-2-yl)anilines and found that (N-(5-(1H-indol-2-yl)-2-methoxyphenyl)sulfamoyl)carbamates and their monofluorinated analogues are potent and selective inhibitors of the linoleate oxygenase activity of rabbit and human ALOX15. Introduction of a 2-methoxyaniline moiety into the core pharmacophore plays a crucial role in substrate-selective inhibition of ALOX15-catalyzed oxygenation of linoleic acid at submicromolar concentrations without affecting arachidonic acid oxygenation. Steady-state kinetics, mutagenesis studies, and molecular dynamics (MD) simulations suggested an allosteric mechanism of action. Using a dimer model of ALOX15, our MD simulations suggest that the binding of the inhibitor at the active site of one monomer induces conformational alterations in the other monomer so that the formation of a productive enzyme-linoleic acid complex is energetically compromised.

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 C14 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 lipoxinB4 after the final action of a reductase. In contrast, our calculations are in agreement with experiments that lipoxinA4 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.

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

The significance of tunneling contributions to the rate constants of enzymatic reactions has been described at length using experimental evidence as well as theoretical computations. Within the frame of variational transition state theory (VTST), tunneling corrections are included using the so-called ground-state tunneling transmission coefficient. For the calculation of those enzymatic rate constants using the ensemble-averaged extension of VTST on QM/MM potential energy surfaces, the transmission coefficient at a given temperature is averaged over a set of coefficient values, each one obtained from an individual minimum energy path (MEP). However, the calculation of accurate QM/MM MEPs for tunneling calculations, also using a reliable QM method like DFT, is highly costly in enzyme models. For this reason, more affordable methodologies have been used. In this paper, we validate a feasible computational strategy to compute multidimensional tunneling corrections that describes better than cheaper alternatives the physics of the hydrogen abstraction from linoleic acid catalyzed by the enzyme 15-rLOX-1. Our recommendations to obtain better values of kinetic isotope effects and, especially, of rate constants are based on multidimensional small-curvature tunneling (SCT) coefficients derived from electrostatic embedding QM(DFT)/MM MEPs. The MEPs used must be calculated with a small enough step-size. Also, the number of gradients and Hessians along the reaction path must be checked to cover the whole tunneling region and to obtain converged adiabatic potential energy profiles. Distinguished reaction coordinates (DCPs) that are commonly used to describe enzyme reaction mechanisms are not adequate for tunneling calculations in such biological systems.

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.

Unraveling how the Gly526Ser mutation arrests prostaglandin formation from arachidonic acid catalyzed by cyclooxygenase-2: a combined molecular dynamics and QM/MM study.

Cyclooxygenases (COXs) are the enzymes responsible for the biosynthesis of prostaglandins, eicosanoids that play a major role in many physiological processes. Particularly, prostaglandins are known to trigger inflammation, and COX-2, the enzyme isoform associated with this inflammatory response, catalyzes the cyclooxidation of arachidonic acid, leading to prostaglandin G2. For this reason, COX-2 has been a very important pharmacological target for several decades now. The catalytic mechanism of COX-2, a so-called all-radical mechanism, consists of six chemical steps. One of the most intriguing aspects of this mechanism is how COX-2 manages to control the regio- and stereospecificity of the products formed at each step. Mutagenesis experiments have previously been performed in an attempt to find those hot-spot residues that make such control possible. In this context, it is worth mentioning that in experiments with the Gly526Ser COX-2 mutant, prostaglandins were not detected. In this paper, we have combined molecular dynamics simulations and quantum mechanics/molecular mechanics calculations to analyze how the COX-2 catalytic mechanism is modified in the Gly526Ser mutant. Therefore, this study provides new insights into the COX-2 catalytic function.

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.

Synthetic Photoswitchable Neurotransmitters Based on Bridged Azobenzenes.

Photoswitchable neurotransmitters of ionotropic kainate receptors were synthesized by tethering a glutamate moiety to disubstituted C2-bridged azobenzenes, which were prepared through a novel methodology that allows access to diazocines with higher yields and versatility. Because of the singular properties of these photochromes, photoisomerizable compounds were obtained with larger thermal stability for their inert cis isomer than for their biologically activity trans state. This enabled selective neuronal firing upon irradiation without background activity in the dark.

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.

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.

Inhibition of Mammalian 15-Lipoxygenase by Three Ebselen-like Drugs. A QM/MM and MM/PBSA Comparative Study.

Ebselen is a potent competitive inhibitor of the active form of rabbit 15-lipoxygenase, an enzyme involved in many inflammatory diseases. Light-induced Z-to-E isomerization of the ebselen-like 2-(3-benzylidene)-3-oxo-2,3-dihydrobenzo[b]thiophene-7-carboxylic acid methyl ester (BODTCM) molecule was used to convert the weak (Z)-BOTDCM inhibitor into the (E)-isomer with much higher inhibitory capacity. In this study, the binding modes of ebselen, (E)-BOTDCM and (Z)-BOTDCM, have been analyzed to provide molecular insights on the inhibitory potency of ebselen and on the geometric-isomer specificity of (E)- and (Z)-BOTDCM inhibitors. The inhibitor-enzyme structures obtained from docking and molecular dynamics simulations as well as from QM/MM calculations show that the inhibitor molecules are not coordinated to the nonheme iron in the active site. Thermal motion allows ebselen and (E)-BOTDCM to visit a wide range of the configurational space competing with the polyunsaturated fatty acid for binding at the active site. Both molecules present similar MM/PBSA binding free energies. The energy penalty for the bigger geometric deformation undergone by (E)-BODTCM would explain its lower inhibitor potency. The (Z)-isomer is the weakest inhibitor because thermal motion moves it to a region very far from the first coordination sphere of Fe, where it could not compete with the fatty acid substrate.

Understanding how cAMP-dependent protein kinase can catalyze phosphoryl transfer in the presence of Ca2+ and Sr2+: a QM/MM study.

Recent experimental results have challenged conventional views on the role metals play in the chemistry of protein kinases because it has been shown that (cAMP)-dependent protein kinase (PKA) is active in the presence of other divalent alkaline earth metal cations besides physiological Mg2+ ions. This has raised the important possibility that Ca2+ may also be a physiological cofactor of protein kinases. In this work, QM/MM calculations, at the DFT and MP2 levels for the QM part, on complete solvated models of PKAc-M2ATP-substrate ternary complexes, with PKAc as the catalytic subunit of PKA, M denoting Ca2+ or Sr2+ and substrate denoting SP20 or Kemptide, have been carried out for the overall phosphoryl transfer reaction. In accordance with the experimental data, our theoretical results show for the first time at the molecular level how the overall PKAc-catalyzed phosphorylation of SP20, via a dissociative mechanism, is plausible with Ca2+ and Sr2+. The viability of the catalytic reaction with Kemptide and Ca2+ is also verified here. The energy barrier of the rate-limiting phosphoryl-transfer step does not depend on different coordination environments of the alkaline earth metal cations whereas the proton-transfer step region is metal dependent making the global chemical process more exoergic on going from Mg2+ to Sr2+. This trend is in agreement with the less effective release of the phosphorylated product observed experimentally in the presence of Ca2+versus Mg2+, and would explain also the lower activity of PKAc with Ca2+, since phospho-substrate and ADP releases are rate limiting for catalytic turnover. For the same reason, we predict an even lower activity of PKAc with Sr2+. Moreover, the active sites of the in silico reactant and product complexes and the available X-ray crystallographic structures show good agreement.

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.