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

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

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.

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.

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.

Conformational diversity induces nanosecond-timescale chemical disorder in the HIV-1 protease reaction pathway

The role of conformational diversity in enzyme catalysis has been a matter of analysis in recent studies. Pre-organization of the active site has been pointed out as the major source for enzymes' catalytic power. Following this line of thought, it is becoming clear that specific, instantaneous, non-rare enzyme conformations that make the active site perfectly pre-organized for the reaction lead to the lowest activation barriers that mostly contribute to the macroscopically observed reaction rate. The present work is focused on exploring the relationship between structure and catalysis in HIV-1 protease (PR) with an adiabatic mapping method, starting from different initial structures, collected from a classical MD simulation. The first, rate-limiting step of the HIV-1 PR catalytic mechanism was studied with the ONIOM QM/MM methodology (B3LYP/6-31G(d):ff99SB), with activation and reaction energies calculated at the M06-2X/6-311++G(2d,2p):ff99SB level of theory, in 19 different enzyme:substrate conformations. The results showed that the instantaneous enzyme conformations have two independent consequences on the enzyme's chemistry: they influence the barrier height, something also observed in the past in other enzymes, and they also influence the specific reaction pathway, which is something unusual and unexpected, challenging the "one enzyme-one substrate-one reaction mechanism" paradigm. Two different reaction mechanisms, with similar reactant probabilities and barrier heights, lead to the same gem-diol intermediate. Subtle nanosecond-timescale rearrangements in the active site hydrogen bonding network were shown to determine which reaction the enzyme follows. We named this phenomenon chemical disorder. The results make us realize the unexpected mechanistic consequences of conformational diversity in enzymatic reactivity.