Application of a SCC-DFTB QM/MM approach to the investigation of the catalytic mechanism of fatty acid amide hydrolase

Mechanistic Modeling of Monoglyceride Lipase Covalent Modification Elucidates the Role of Leaving Group Expulsion and Discriminates Inhibitors with High and Low Potency

Inhibition of monoglyceride lipase (MGL), also known as monoacylglycerol lipase (MAGL), has emerged as a promising approach for treating neurological diseases. To gain useful insights in the design of agents with balanced potency and reactivity, we investigated the mechanism of MGL carbamoylation by the reference triazole urea SAR629 (IC₅₀ = 0.2 nM) and two recently described inhibitors featuring a pyrazole (IC₅₀ = 1800 nM) or a 4-cyanopyrazole (IC₅₀ = 8 nM) leaving group (LG), using a hybrid quantum mechanics/molecular mechanics (QM/MM) approach. Opposite to what was found for substrate 2-arachidonoyl-sn-glycerol (2-AG), covalent modification of MGL by azole ureas is controlled by LG expulsion. Simulations indicated that changes in the electronic structure of the LG greatly affect reaction energetics with triazole and 4-cyanopyrazole inhibitors following a more accessible carbamoylation path compared to the unsubstituted pyrazole derivative. The computational protocol provided reaction barriers able to discriminate between MGL inhibitors with different potencies. These results highlight how QM/MM simulations can contribute to elucidating structure–activity relationships and provide insights for the design of covalent inhibitors.

A Review of Recent Developments in Molecular Dynamics Simulations of the Photoelectrochemical Water Splitting Process

In this review, we provide a short overview of the Molecular Dynamics (MD) method and how it can be used to model the water splitting process in photoelectrochemical hydrogen production. We cover classical non-reactive and reactive MD techniques as well as multiscale extensions combining classical MD with quantum chemical and continuum methods. Selected examples of MD investigations of various aqueous semiconductor interfaces with a special focus on TiO2 are discussed. Finally, we identify gaps in the current state-of-the-art where further developments will be needed for better utilization of MD techniques in the field of water splitting.

Synthesis, Molecular Modeling and Biological Evaluation of Metabolically Stable Analogues of the Endogenous Fatty Acid Amide Palmitoylethanolamide

Palmitoylethanolamide (PEA) belongs to the class of N-acylethanolamine and is an endogenous lipid potentially useful in a wide range of therapeutic areas; products containing PEA are licensed for use in humans as a nutraceutical, a food supplement, or food for medical purposes for its analgesic and anti-inflammatory properties demonstrating efficacy and tolerability. However, the exogenously administered PEA is rapidly inactivated; in this process, fatty acid amide hydrolase (FAAH) plays a key role both in hepatic metabolism and in intracellular degradation. So, the aim of the present study was the design and synthesis of PEA analogues that are more resistant to FAAH-mediated hydrolysis. A small library of PEA analogues was designed and tested by molecular docking and density functional theory calculations to find the more stable analogue. The computational investigation identified RePEA as the best candidate in terms of both synthetic accessibility and metabolic stability to FAAH-mediated hydrolysis. The selected compound was synthesized and assayed ex vivo to monitor FAAH-mediated hydrolysis and to confirm its anti-inflammatory properties. ¹H-NMR spectroscopy performed on membrane samples containing FAAH in integral membrane protein demonstrated that RePEA is not processed by FAAH, in contrast with PEA. Moreover, RePEA retains PEA’s ability to inhibit LPS-induced cytokine release in both murine N9 microglial cells and human PMA-THP-1 cells.

Benchmarking quantum mechanical methods for calculating reaction energies of reactions catalyzed by enzymes

To assess the accuracy of different quantum mechanical methods for biochemical modeling, the reaction energies of 20 small model reactions (chosen to represent chemical steps catalyzed by commonly studied enzymes) were calculated. The methods tested included several popular Density Functional Theory (DFT) functionals, second-order Møller Plesset perturbation theory (MP2) and its spin-component scaled variant (SCS-MP2), and coupled cluster singles and doubles and perturbative triples (CCSD(T)). Different basis sets were tested. CCSD(T)/aug-cc-pVTZ results for all 20 reactions were used to benchmark the other methods. It was found that MP2 and SCS-MP2 reaction energy calculation results are similar in quality to CCSD(T) (mean absolute error (MAE) of 1.2 and 1.3 kcal mol−1, respectively). MP2 calculations gave a large error in one case, and are more subject to basis set effects, so in general SCS-MP2 calculations are a good choice when CCSD(T) calculations are not feasible. Results with different DFT functionals were of reasonably good quality (MAEs of 2.5–5.1 kcal mol−1), whereas popular semi-empirical methods (AM1, PM3, SCC-DFTB) gave much larger errors (MAEs of 11.6–14.6 kcal mol−1). These results should be useful in guiding methodological choices and assessing the accuracy of QM/MM calculations on enzyme-catalyzed reactions.

Design and SAR Analysis of Covalent Inhibitors Driven by Hybrid QM/MM Simulations

Quantum mechanics/molecular mechanics (QM/MM) hybrid technique is emerging as a reliable computational method to investigate and characterize chemical reactions occurring in enzymes. From a drug discovery perspective, a thorough understanding of enzyme catalysis appears pivotal to assist the design of inhibitors able to covalently bind one of the residues belonging to the enzyme catalytic machinery. Thanks to the current advances in computer power, and the availability of more efficient algorithms for QM-based simulations, the use of QM/MM methodology is becoming a viable option in the field of covalent inhibitor design. In the present review, we summarized our experience in the field of QM/MM simulations applied to drug design problems which involved the optimization of agents working on two well-known drug targets, namely fatty acid amide hydrolase (FAAH) and epidermal growth factor receptor (EGFR). In this context, QM/MM simulations gave valuable information in terms of geometry (i.e., of transition states and metastable intermediates) and reaction energetics that allowed to correctly predict inhibitor binding orientation and substituent effect on enzyme inhibition. What is more, enzyme reaction modelling with QM/MM provided insights that were translated into the synthesis of new covalent inhibitor featured by a unique combination of intrinsic reactivity, on-target activity, and selectivity.

Quantum mechanics/molecular mechanics modeling of covalent addition between EGFR-cysteine 797 and N-(4-anilinoquinazolin-6-yl) acrylamide

Irreversible epidermal growth factor receptor (EGFR) inhibitors can circumvent resistance to first-generation ATP-competitive inhibitors in the treatment of nonsmall-cell lung cancer. They covalently bind a noncatalytic cysteine (Cys797) at the surface of EGFR active site by an acrylamide warhead. Herein, we used a hybrid quantum mechanics/molecular mechanics (QM/MM) potential in combination with umbrella sampling in the path-collective variable space to investigate the mechanism of alkylation of Cys797 by the prototypical covalent inhibitor N-(4-anilinoquinazolin-6-yl) acrylamide. Calculations show that Cys797 reacts with the acrylamide group of the inhibitor through a direct addition mechanism, with Asp800 acting as a general base/general acid in distinct steps of the reaction. The obtained reaction free energy is negative (ΔA = -12 kcal/mol) consistent with the spontaneous and irreversible alkylation of Cys797 by N-(4-anilinoquinazolin-6-yl) acrylamide. Our calculations identify desolvation of Cys797 thiolate anion as a key step of the alkylation process, indicating that changes in the intrinsic reactivity of the acrylamide would have only a minor impact on the inhibitor potency.

Computational insights into function and inhibition of fatty acid amide hydrolase

The Fatty Acid Amide Hydrolase (FAAH) enzyme is a membrane-bound serine hydrolase responsible for the deactivating hydrolysis of a family of naturally occurring fatty acid amides. FAAH is a critical enzyme of the endocannabinoid system, being mainly responsible for regulating the level of its main cannabinoid substrate anandamide. For this reason, pharmacological inhibition of FAAH, which increases the level of endogenous anandamide, is a promising strategy to cure a variety of diseases including pain, inflammation, and cancer. Much structural, mutagenesis, and kinetic data on FAAH has been generated over the last couple of decades. This has prompted several informative computational investigations to elucidate, at the atomic-level, mechanistic details on catalysis and inhibition of this pharmaceutically relevant enzyme. Here, we review how these computational studies - based on classical molecular dynamics, full quantum mechanics, and hybrid QM/MM methods - have clarified the binding and reactivity of some relevant substrates and inhibitors of FAAH. We also discuss the experimental implications of these computational insights, which have provided a thoughtful elucidation of the complex physical and chemical steps of the enzymatic mechanism of FAAH. Finally, we discuss how computations have been helpful for building structure-activity relationships of potent FAAH inhibitors.

Unraveling the differences of the hydrolytic activity of Trypanosoma cruzi trans-sialidase and Trypanosoma rangeli sialidase: a quantum mechanics-molecular mechanics modeling study

Chagas' disease, also known as American trypanosomiasis, is a lethal, chronic disease that currently affects more than 10 million people in Central and South America. The trans-sialidase from Trypanosoma cruzi (T. cruzi, TcTS) is a crucial enzyme for the survival of this parasite: sialic acids from the host are transferred to the cell surface glycoproteins of the trypanosome, thereby evading the host's immune system. On the other hand, the sialidase of T. rangeli (TrSA), which shares 70% sequence identity with TcTS, is a strict hydrolase and shows no trans-sialidase activity. Therefore, TcTS and TrSA represent an excellent framework to understand how different catalytic activities can be achieved with extremely similar structures. By means of combined quantum mechanics-molecular mechanics (QM/MM, SCC-DFTB/Amberff99SB) calculations and umbrella sampling simulations, we investigated the hydrolysis mechanisms of TcTS and TrSA and computed the free energy profiles of these reactions. The results, together with our previous computational investigations, are able to explain the catalytic mechanism of sialidases and describe how subtle differences in the active site make TrSA a strict hydrolase and TcTS a more efficient trans-sialidase.

Comparison of ab Initio, DFT, and Semiempirical QM/MM Approaches for Description of Catalytic Mechanism of Hairpin Ribozyme

We have analyzed the capability of state-of-the-art multiscale computational approaches to provide atomic-resolution electronic structure insights into possible catalytic scenarios of the hairpin ribozyme by evaluating potential and free energy surfaces of the reactions by various hybrid QM/MM methods. The hairpin ribozyme is a unique catalytic RNA that achieves rate acceleration similar to other small self-cleaving ribozymes but without direct metal ion participation. Guanine 8 (G8) and adenine 38 (A38) have been identified as the catalytically essential nucleobases. However, their exact catalytic roles are still being investigated. In line with the available experimental data, we considered two reaction scenarios involving protonated A38H(+) as a general acid which is further assisted by either canonical G8 or deprotonated G8(-) forms. We used the spin-component scaled Møller-Plesset (SCS-MP2) method at the complete basis set limit as the reference method. The semiempirical AM1/d-PhoT and SCC-DFTBPR methods provided acceptable activation barriers with respect to the SCS-MP2 data but predicted significantly different reaction pathways. DFT functionals (BLYP and MPW1K) yielded the same reaction pathway as the SCS-MP2 method. The activation barriers were slightly underestimated by the GGA BLYP functional, although with accuracy comparable to the semiempirical methods. The SCS-MP2 method and hybrid MPW1K functional gave activation barriers that were closest to those derived from experimentally measured rate constants.

Quantum mechanics/molecular mechanics modeling of fatty acid amide hydrolase reactivation distinguishes substrate from irreversible covalent inhibitors

Carbamate and urea derivatives are important classes of fatty acid amide hydrolase (FAAH) inhibitors that carbamoylate the active-site nucleophile Ser241. In the present work, the reactivation mechanism of carbamoylated FAAH is investigated by means of a quantum mechanics/molecular mechanics (QM/MM) approach. The potential energy surfaces for decarbamoylation of FAAH covalent adducts, derived from the O-aryl carbamate URB597 and from the N-piperazinylurea JNJ1661610, were calculated and compared to that for deacylation of FAAH acylated by the substrate oleamide. Calculations show that a carbamic group bound to Ser241 prevents efficient stabilization of transition states of hydrolysis, leading to large increments in the activation barrier. Moreover, the energy barrier for the piperazine carboxylate was significantly lower than that for the cyclohexyl carbamate derived from URB597. This is consistent with experimental data showing slowly reversible FAAH inhibition for the N-piperazinylurea inhibitor and irreversible inhibition for URB597.

Wagging the Tail: Essential Role of Substrate Flexibility in FAAH Catalysis

The serine hydrolase, fatty acid amide hydrolase (FAAH), is responsible for the intracellular degradation of anandamide and other bioactive fatty acid ethanolamides involved in the regulation of pain, inflammation, and other pathophysiological processes. The catalytic site of FAAH is composed of multiple cavities with mixed hydrophobic and hydrophilic properties, the role of which remains incompletely understood. Anandamide is thought to enter the active site through a "membrane-access" (MA) channel and position its flexible fatty acyl chain in a highly hydrophobic "acyl chain-binding" (AB) cavity to allow for hydrolysis to occur. Using microsecond molecular dynamics (MD) simulations of FAAH embedded in a realistic membrane/water environment, we show now that anandamide may not lock itself into the AB cavity but may rather assume catalytically significant conformations required for hydrolysis by moving its flexible arachidonoyl tail between the MA and AB cavities. This process is regulated by a phenylalanine residue (Phe432) located at the boundary between the two cavities, which may act as a "dynamic paddle." The results identify structural flexibility as a key determinant by which FAAH recognizes its primary lipid substrate.

A catalytic mechanism for cysteine N-terminal nucleophile hydrolases, as revealed by free energy simulations

The N-terminal nucleophile (Ntn) hydrolases are a superfamily of enzymes specialized in the hydrolytic cleavage of amide bonds. Even though several members of this family are emerging as innovative drug targets for cancer, inflammation, and pain, the processes through which they catalyze amide hydrolysis remains poorly understood. In particular, the catalytic reactions of cysteine Ntn-hydrolases have never been investigated from a mechanistic point of view. In the present study, we used free energy simulations in the quantum mechanics/molecular mechanics framework to determine the reaction mechanism of amide hydrolysis catalyzed by the prototypical cysteine Ntn-hydrolase, conjugated bile acid hydrolase (CBAH). The computational analyses, which were confirmed in water and using different CBAH mutants, revealed the existence of a chair-like transition state, which might be one of the specific features of the catalytic cycle of Ntn-hydrolases. Our results offer new insights on Ntn-mediated hydrolysis and suggest possible strategies for the creation of therapeutically useful inhibitors.

The increasing role of QM/MM in drug discovery

Since its first appearance in 1976, the quantum mechanics/molecular mechanics (QM/MM) approach has mostly been used to study the chemical reactions of enzymes, which are frequently the target of drug discovery programs. In principle, a detailed understanding of the enzymatic mechanism should help researchers to design a potent enzyme inhibitor or new drug. However, QM/MM has not yet had a widespread impact on structure-based drug design. This is mostly due to its high computational cost. We expect this to change with the recent and extraordinary increases in computational power, and with the availability of more efficient algorithms for QM/MM calculations. Here, we report on some representative examples of QM/MM studies, including our own research, of pharmaceutically relevant enzymes, such as ribonuclease H and fatty acid amide hydrolase (FAAH). We aim to show how QM/MM has traditionally been used to study enzymatic catalysis. In this regard, we discuss its potential to become a routinely used drug design tool. To support this, we also discuss selected computational studies where QM/MM insights have been helpful in improving the potency of covalent inhibitors of FAAH.

Mechanism of C-terminal intein cleavage in protein splicing from QM/MM molecular dynamics simulations

Protein splicing is a post-translational process in which a biologically inactive protein is activated by the release of a segment denoted as an intein. The process involves four steps. In the third, the scission of the intein takes place after the cyclization of the last amino acid of the segment, an asparagine. Little is known about the chemical reaction necessary for this cyclization. Experiments demonstrate that two histidines (the penultimate amino acid of the intein, and a histidine located 10 amino acids upstream) are relevant in the cyclization of the asparagine. We have investigated the mechanism and determinants of reaction in the GyrA intein focusing on the requirements for asparagine activation for its cyclization. First, the influence that the protonation states of these two histidines have on the orientation of the asparagine side chain is investigated by means of molecular dynamics simulation. Molecular dynamics simulations using the CHARMM27 force field were carried out on the three possible protonation states for each of these two histidines. The results indicate that the only protonation state in which the conformation of the system is suitable for cyclization is when the penultimate histidine is fully protonated (positively charged), and the upstream histidine is in the His(ε) neutral tautomeric form. The free energy profile for the reaction in which the asparagine is activated by a proton transfer to the upstream histidine is presented, computed by hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics at the SCCDFTB/CHARMM27 level of theory. The calculated free energy barrier for the reaction is 19.0 kcal mol(-1). B3LYP/6-31+G(d) QM/MM single-point calculations give a qualitatively a similar energy profile, although with somewhat higher energy barriers, in good agreement with the value derived from experiment of 25 kcal mol(-1) at 60 °C. QM/MM molecular dynamics simulations of the reactant, activated reactant and intermediate states highlight the importance of the Arg181-Val182-Asp183 segment in catalysing the reaction. Overall, the results indicate that nucleophilic activation of the asparagine for its cyclization by the upstream histidine acting as the base is a plausible mechanism for the C-terminal cleavage in protein splicing.