Conformational effects on the pro-S hydrogen abstraction reaction in cyclooxygenase-1: an integrated QM/MM and MD study

A key step in the cyclooxygenase reaction cycle of cyclooxygenase 1 (COX-1) is abstraction of the pro-S hydrogen atom of the arachidonic acid by a radical that is formed at the protein residue Tyr-385. Here we investigate this reaction step by a quantum-mechanics/molecular-mechanics approach in combination with molecular-dynamics simulations. The simulations identify the hydrogen abstraction angle as a crucial geometric determinant of the reaction, thus revealing the importance of the cyclooxygenase active site for calculating the potential energy surface of the reaction

Impact of Warhead Modulations on the Covalent Inhibition of SARS-CoV2 Mpro Explored by QM\/MM Simulations

The COVID-19 pandemic, caused by the severe acute respiratory syndrome coronavirus-2, SARS-CoV-2, shows the need for effective antiviral treatments. Here, we present a simulation study of the inhibition of the SARS-CoV-2 main protease (Mpro), a cysteine hydrolase essential for the life cycle of the virus. The free energy landscape for the mechanism of the inhibition process is explored by QM/MM umbrella sampling and free energy perturbation simulations at the M06-2X/MM level of theory for two proposed peptidyl covalent inhibitors that share the same recognition motif but feature distinct cysteine-targeting warheads. Regardless of the intrinsic reactivity of the modeled inhibitors, namely a Michael acceptor and a hydroxymethyl ketone activated carbonyl, our results confirm that the inhibitory process takes place by means of a two-step mechanism, in which the formation of an ion pair C145/H41 dyad precedes the protein–inhibitor covalent bond formation. The nature of this second step is strongly dependent on the functional groups in the warhead: while the nucleophilic attack of the C145 sulfur atom on the Cα of the double bond of the Michael acceptor takes place concertedly with the proton transfer from H41 to Cβ, in the compound with an activated carbonyl, the sulfur attacks the carbonyl carbon concomitant with a proton transfer from H41 to the carbonyl oxygen via the hydroxyl group. An analysis of the free energy profiles, structures along the reaction path, and interactions between the inhibitors and the different pockets of the active site on the protein shows a measurable effect of the warhead on the kinetics and thermodynamics of the process. These results and QM/MM methods can be used as a guide to select warheads to design efficient irreversible and reversible inhibitors of SARS-CoV-2 Mpro

Metadynamics Simulations Distinguish Short- and Long-Residence-Time Inhibitors of Cyclin-Dependent Kinase 8.

The duration of drug efficacy in vivo is a key aspect primarily addressed during the lead optimization phase of drug discovery. Hence, the availability of robust computational approaches that can predict the residence time of a compound at its target would accelerate candidate selection. Nowadays the theoretical prediction of this parameter is still very challenging. Starting from methods reported in the literature, we set up and validated a new metadynamics (META-D)-based protocol that was used to rank the experimental residence times of 10 arylpyrazole cyclin-dependent kinase 8 (CDK8) inhibitors for which target-bound X-ray structures are available. The application of reported methods based on the detection of the escape from the first free energy well gave a poor correlation with the experimental values. Our protocol evaluates the energetics of the whole unbinding process, accounting for multiple intermediates and transition states. Using seven collective variables (CVs) encoding both roto-translational and conformational motions of the ligand, a history-dependent biasing potential is deposited as a sum of constant-height Gaussian functions until the ligand reaches an unbound state. The time required to achieve this state is proportional to the integral of the deposited potential over the CV hyperspace. Average values of this time, for replicated META-D simulations, provided an accurate classification of CDK8 inhibitors spanning short, medium, and long residence times

Mechanistic Modeling of Lys745 Sulfonylation in EGFR C797S Reveals Chemical Determinants for Inhibitor Activity and Discriminates Reversible from Irreversible Agents

Targeted covalent inhibitors hold promise for drug discovery, particularly for kinases. Targeting the catalytic lysine of epidermal growth factor receptor (EGFR) has attracted attention as a new strategy to overcome resistance due to the emergence of C797S mutation. Sulfonyl fluoride derivatives able to inhibit EGFRL858R/T790M/C797S by sulfonylation of Lys745 have been reported. However, atomistic details of this process are still poorly understood. Here, we describe the mechanism of inhibition of an innovative class of compounds that covalently engage the catalytic lysine of EGFR, through a sulfur(VI) fluoride exchange (SuFEx) process, with the help of hybrid quantum mechanics/molecular mechanics (QM/MM) and path collective variables (PCVs) approaches. Our simulations identify the chemical determinants accounting for the irreversible activity of agents targeting Lys745 and provide hints for the further optimization of sulfonyl fluoride agents

Molecular modelling studies on PEA hydrolysis by the enzyme <em>N</em>-acylethanolamine acid amidase

n/

Mechanistic Insights into the Reaction of Chlorination of Tryptophan Catalyzed by Tryptophan 7-Halogenase

Tryptophan 7-halogenase catalyzes chlorination of free tryptophan to 7-chlorotryptophan, which is the first step in the antibiotic pyrrolnitrin biosynthesis. Many biologically and pharmaceutically active natural products contain chlorine and thus, an understanding of the mechanism of its introduction into organic molecules is important. Whilst enzyme-catalyzed chlorination is accomplished with ease, it remains a difficult task for the chemists. Therefore, utilizing enzymes in the synthesis of chlorinated organic compounds is important, and providing atomistic mechanistic insights about the reaction mechanism of tryptophan 7-halogenase is vital and timely. In this work, we examined a mechanism for the reaction of tryptophan chlorination, performed by tryptophan 7-halogenase, by calculating potential energy and free energy surfaces using two different Combined Quantum Mechanical/Molecular Mechanical (QM/MM) methods both employing Density Functional Theory (DFT) for the QM region. Both computational strategies agree on the nature of the rate-limiting step and provided close results for the reaction barriers of the two reaction steps. The calculations for both the potential energy and the free energy profiles showed very similar geometric features and hydrogen bonding interactions for the characterized stationary points.Peer ReviewedPostprint (published version

Phenotype Screening of an Azole-bisindole Chemical Library Identifies URB1483 as a New Antileishmanial Agent Devoid of Toxicity on Human Cells

none11sìWe report the evaluation of a small library of azole-bisindoles for their antileishmanial potential, in terms of efficacy on Leishmania infantum promastigotes and intracellular amastigotes. Nine compounds showed good activity on L. infantum MHOM/TN/80/IPT1 promastigotes with IC50 values ranging from 4 to 10 μM. These active compounds were also tested on human (THP-1, HEPG2, HaCaT, and human primary fibroblasts) and canine (DH82) cell lines. URB1483 was selected as the best compound, with no quantifiable cytotoxicity in mammalian cells, to test the efficacy on intracellular amastigotes. URB1483 significantly reduced the infection index of both human and canine macrophages with an effect comparable to the clinically used drug pentamidine. URB1483 emerges as a new anti-infective agent with remarkable antileishmanial activity and no cytotoxic effects on human and canine cells.openDiotallevi, Aurora; Scalvini, Laura; Buffi, Gloria; Pérez-Pertejo, Yolanda; De Santi, Mauro; Verboni, Michele; Favi, Gianfranco; Magnani, Mauro; Lodola, Alessio; Lucarini, Simone; Galluzzi, LucaDiotallevi, Aurora; Scalvini, Laura; Buffi, Gloria; Pérez-Pertejo, Yolanda; De Santi, Mauro; Verboni, Michele; Favi, Gianfranco; Magnani, Mauro; Lodola, Alessio; Lucarini, Simone; Galluzzi, Luc

Metadynamics for perspective drug design: Computationally driven synthesis of new protein-protein interaction inhibitors targeting the EphA2 receptor

Metadynamics (META-D) is emerging as a powerful method for the computation of the multidimensional freeenergy surface (FES) describing the protein-ligand binding process. Herein, the FES of unbinding of the antagonist N-(3α-hydroxy-5β-cholan-24-oyl)-L-β-homotryptophan (UniPR129) from its EphA2 receptor was reconstructed by META-D simulations. The characterization of the free-energy minima identified on this FES proposes a binding mode fully consistent with previously reported and new structure-activity relationship data. To validate this binding mode, new N-(3α-hydroxy-5β-cholan-24-oyl)-L-β-homotryptophan derivatives were designed, synthesized, and tested for their ability to displace ephrin-A1 from the EphA2 receptor. Among them, two antagonists, namely compounds 21 and 22, displayed high affinity versus the EphA2 receptor and resulted endowed with better physicochemical and pharmacokinetic properties than the parent compound. These findings highlight the importance of free-energy calculations in drug design, confirming that META-D simulations can be used to successfully design novel bioactive compounds

Benzisothiazolinone Derivatives as Potent Allosteric Monoacylglycerol Lipase Inhibitors That Functionally Mimic Sulfenylation of Regulatory Cysteines

We describe a set of benzisothiazolinone (BTZ) derivatives that are potent inhibitors of monoacylglycerol lipase (MGL), the primary degrading enzyme for the endocannabinoid 2-arachidonoyl-sn-glycerol (2-AG). Structure-activity relationship studies evaluated various substitutions on the nitrogen atom and the benzene ring of the BTZ nucleus. Optimized derivatives with nanomolar potency allowed us to investigate the mechanism of MGL inhibition. Site-directed mutagenesis and mass spectrometry experiments showed that BTZs interact in a covalent reversible manner with regulatory cysteines, Cys201 and Cys208, causing a reversible sulfenylation known to modulate MGL activity. Metadynamics simulations revealed that BTZ adducts favor a closed conformation of MGL that occludes substrate recruitment. The BTZ derivative 13 protected neuronal cells from oxidative stimuli and increased 2-AG levels in the mouse brain. The results identify Cys201 and Cys208 as key regulators of MGL function and point to the BTZ scaffold as a useful starting point for the discovery of allosteric MGL inhibitors

Peptidyl Nitroalkene Inhibitors of Main Protease (Mpro) rationalized by Computational/1 Crystallographic 2 Investigations as Antivirals against SARS-CoV-2

34 p.-10 fig.-2 tab.-2 schem.The coronavirus disease 2019 (COVID-19) pandemic continues to represent a global public health issue. The viral main protease (Mpro) represents one of the most attractive targets for the development of antiviral drugs. Herein we report peptidyl nitroalkenes exhibited enzyme inhibitory activity against Mpro (Ki: 1-10 μM) and three of them good anti-SARS-CoV-2 infection activity in the low micromolar range (EC50: 1-12 μM) without significant toxicity. Additional kinetic studies of compounds FGA145, FGA146 and FGA147 show that all three compounds inhibit Cathepsin L, denoting a possible multitarget effect of these compounds in the antiviral activity. QM/MM computer simulations assisted in the design and in elucidating the way of action. Finally, structural analysis shows, in agreement with the computer predictions, the binding mode of FGA146 and FGA147 to the active site of the protein. Our results illustrate that peptidyl nitroalkenes are potent covalent reversible inhibitors of the Mpro and cathepsin L, and that inhibitors FGA145, FGA146 and FGA147 prevent infection becoming promising drugs against SARS-CoV-2.This research was funded by the Consejo Superior de Investigaciones Científicas, grant number PIE- 202020E224, the Spanish Ministerio de Ciencia e Innovación (ref. PID2021-123332OB-C21 and PID2019- 107098RJ-I00), the Generalitat Valenciana (PROMETEO with ref. CIPROM/2021/079, and SEJI/2020/007), Universitat Jaume I (UJI-B2020-03, UJI-B2021-71 and SomUJIcontracovid crowdfunding campaign).K.Ś.thanks to Ministerio de Ciencia e Innovación and Fondo Social Europeo for a Ramon y Cajal contract (Ref. RYC2020-030596-I). The authors wish to thank the staff of beamlines ID30B (ESRF Synchrotron) and BL13- XALOC (ALBA Synchrotron) for their generous and much appreciated support, and the Serveis Centrals d’Instrumentació Científica of Universitat Jaume I for technical support. Finally, the authors acknowledge the computer resources at Mare Nostrum of the Barcelona Supercomputing Center (QH-2022-2-0004 and QH- 2022-3-0008), as well as the local computational resources founded by Generalitat Valenciana - European Regional Development Fund (REF: IDIFEDER/2021/02).Peer reviewe