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.

Unbinding Kinetics of Muscarinic M3 Receptor Antagonists Explained by Metadynamics Simulations.

The residence time (RT), the time for which a drug remains bound to its biological target, is a critical parameter for drug design. The prediction of this key kinetic property has been proven to be challenging and computationally demanding in the framework of atomistic simulations. In the present work, we setup and applied two distinct metadynamics protocols to estimate the RTs of muscarinic M3 receptor antagonists. In the first method, derived from the conformational flooding approach, the kinetics of unbinding is retrieved from a physics-based parameter known as the acceleration factor α (i.e., the running average over time of the potential deposited in the bound state). Such an approach is expected to recover the absolute RT value for a compound of interest. In the second method, known as the tMETA-D approach, a qualitative estimation of the RT is given by the time of simulation required to drive the ligand from the binding site to the solvent bulk. This approach has been developed to reproduce the change of experimental RTs for compounds targeting the same target. Our analysis shows that both computational protocols are able to rank compounds in agreement with their experimental RTs. Quantitative structure-kinetics relationship (SKR) models can be identified and employed to predict the impact of a chemical modification on the experimental RT once a calibration study has been performed.

Molecular Determinants of EphA2 and EphB2 Antagonism Enable the Design of Ligands with Improved Selectivity.

With the aim of identifying novel antagonists selective for the EphA receptor family, a combined experimental and computational approach was taken to investigate the molecular basis of the recognition between a prototypical Eph-ephrin antagonist (UniPR1447) and two representative receptors of the EphA and EphB subfamilies, namely, EphA2 and EphB2 receptors. The conformational free-energy surface (FES) of the binding state of UniPR1447 within the ligand binding domain of EphA2 and EphB2, reconstructed from molecular dynamics (MD) simulations performed on the microsecond time scale, was exploited to drive the design and synthesis of a novel antagonist selective for EphA2 over the EphB2 receptor. The availability of compounds with this pharmacological profile will help discriminate the importance of these two receptors in the insurgence and progression of cancer.

Free-Energy Simulations Support a Lipophilic Binding Route for Melatonin Receptors.

The effects of the neurohormone melatonin are mediated by the activation of the GPCRs MT1 and MT2 in a variety of tissues. Crystal structures suggest ligand access to the orthosteric binding site of MT1 and MT2 receptors through a lateral channel between transmembrane (TM) helices IV and V. We investigated the feasibility of this lipophilic entry route for 2-iodomelatonin, a nonselective agonist with a slower dissociation rate from the MT2 receptor, applying enhanced sampling simulations and free-energy calculations. 2-Iodomelatonin unbinding was investigated with steered molecular dynamics simulations which revealed different trajectories passing through the gap between TM helices IV and V for both receptors. For one of these unbinding trajectories from the MT1 receptor, an umbrella-sampling protocol with path-collective variables provided a calculated energy barrier consistent with the experimental dissociation rate. The side-chain flexibility of Tyr5.38 was significantly different in the two receptor subtypes, as assessed by metadynamics simulations, and during ligand unbinding it frequently assumes an open conformation in the MT1 but not in the MT2 receptor, favoring 2-iodomelatonin egress. Taken together, our simulations are consistent with the possibility that the gap between TM IV and V is a way of connecting the orthosteric binding site and the membrane core for lipophilic melatonin receptor ligands. Our simulations also suggest that the open state of Tyr5.38 generates a small pocket on the surface of MT1 receptor, which could participate in the recognition of MT1-selective ligands and may be exploited in the design of new selective compounds.

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 (IC50 = 0.2 nM) and two recently described inhibitors featuring a pyrazole (IC50 = 1800 nM) or a 4-cyanopyrazole (IC50 = 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.

Different roles for the acyl chain and the amine leaving group in the substrate selectivity of N-Acylethanolamine acid amidase.

N-acylethanolamine acid amidase (NAAA) is an N-terminal nucleophile (Ntn) hydrolase that catalyses the intracellular deactivation of the endogenous analgesic and anti-inflammatory agent palmitoylethanolamide (PEA). NAAA inhibitors counteract this process and exert marked therapeutic effects in animal models of pain, inflammation and neurodegeneration. While it is known that NAAA preferentially hydrolyses saturated fatty acid ethanolamides (FAEs), a detailed profile of the relationship between catalytic efficiency and fatty acid-chain length is still lacking. In this report, we combined enzymatic and molecular modelling approaches to determine the effects of acyl chain and polar head modifications on substrate recognition and hydrolysis by NAAA. The results show that, in both saturated and monounsaturated FAEs, the catalytic efficiency is strictly dependent upon fatty acyl chain length, whereas there is a wider tolerance for modifications of the polar heads. This relationship reflects the relative stability of enzyme-substrate complexes in molecular dynamics simulations.

Analysis of ADAM12-Mediated Ephrin-A1 Cleavage and Its Biological Functions.

Accumulating evidence indicates that an elevated ephrin-A1 expression is positively correlated with a worse prognosis in some cancers such as colon and liver cancer. The detailed mechanism of an elevated ephrin-A1 expression in a worse prognosis still remains to be fully elucidated. We previously reported that ADAM12-cleaved ephrin-A1 enhanced lung vascular permeability and thereby induced lung metastasis. However, it is still unclear whether or not cleaved forms of ephrin-A1 are derived from primary tumors and have biological activities. We identified the ADAM12-mediated cleavage site of ephrin-A1 by a Matrix-assisted laser desorption ionization mass spectrometry and checked levels of ephrin-A1 in the serum and the urine derived from the primary tumors by using a mouse model. We found elevated levels of tumor-derived ephrin-A1 in the serum and the urine in the tumor-bearing mice. Moreover, inhibition of ADAM-mediated cleavage of ephrin-A1 or antagonization of the EphA receptors resulted in a significant reduction of lung metastasis. The results suggest that tumor-derived ephrin-A1 is not only a potential biomarker to predict lung metastasis from the primary tumor highly expressing ephrin-A1 but also a therapeutic target of lung metastasis.

New classes of potent heparanase inhibitors from ligand-based virtual screening.

Heparanase is a validated target in cancer therapy and a potential target for several inflammatory pathologies. A ligand-based virtual screening of commercial libraries was performed to expand the chemical space of small-molecule inhibitors. The screening was based on similarity with known inhibitors and was performed in several runs, starting from literature compounds and progressing through newly discovered inhibitors. Among the fifty-five tested compounds, nineteen had IC50 values lower than 5 µM and some showed remarkable potencies. Importantly, tere- and isophthalamides derivatives belong to new structural classes of heparanase inhibitors and some of them showed enzyme affinities (61 and 63, IC50 = 0.32 and 0.12 µM, respectively) similar to those of the most potent small-molecule inhibitors reported so far. Docking studies provided a comprehensive binding hypothesis shared by compounds with significant structural diversity. The most potent inhibitors reduced cell invasiveness and inhibited the expression of proangiogenic factors in tumour cell lines.

N-tert-butyloxycarbonyl-Phe-Leu-Phe-Leu-Phe (BOC2) inhibits the angiogenic activity of heparin-binding growth factors.

The peptides N-tert-butyloxycarbonyl-Phe-Leu-Phe-Leu-Phe (BOC2) and BOC-Met-Leu-Phe (BOC1) are widely used antagonists of formyl peptide receptors (FPRs), BOC2 acting as an FPR1/FPR2 antagonist whereas BOC1 inhibits FPR1 only. Extensive investigations have been performed by using these FPR antagonists as a tool to assess the role of FPRs in physiological and pathological conditions. Based on previous observations from our laboratory, we assessed the possibility that BOC2 may exert also a direct inhibitory effect on the angiogenic activity of vascular endothelial growth factor-A (VEGF-A). Our data demonstrate that BOC2, but not BOC1, inhibits the angiogenic activity of heparin-binding VEGF-A165 with no effect on the activity of the non-heparin-binding VEGF-A121 isoform. Endothelial cell-based bioassays, surface plasmon resonance analysis, and computer modeling indicate that BOC2 may interact with the heparin-binding domain of VEGF-A165, thus competing for heparin interaction and preventing the binding of VEGF-A165 to tyrosine kinase receptor VEGFR2, its phosphorylation and downstream signaling. In addition, BOC2 inhibits the interaction of a variety of heparin-binding angiogenic growth factors with heparin, including fibroblast growth factor 2 (FGF2) whose angiogenic activity is blocked by the compound. Accordingly, BOC2 suppresses the angiogenic potential of human tumor cell lines that co-express VEGF-A and FGF2. Thus, BOC2 appears to act as a novel multi-heparin-binding growth factor antagonist. These findings caution about the interpretation of FPR-focusing experimental data obtained with this compound and set the basis for the design of novel BOC2-derived, FPR independent multi-target angiogenesis inhibitors.

Theoretical Model of EphA2-Ephrin A1 Inhibition.

This work aims at the theoretical description of EphA2-ephrin A1 inhibition by small molecules. Recently proposed ab initio-based scoring models, comprising long-range components of interaction energy, is tested on lithocholic acid class inhibitors of this protein⁻protein interaction (PPI) against common empirical descriptors. We show that, although limited to compounds with similar solvation energy, the ab initio model is able to rank the set of selected inhibitors more effectively than empirical scoring functions, aiding the design of novel compounds.

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.

Exploiting Free-Energy Minima to Design Novel EphA2 Protein-Protein Antagonists: From Simulation to Experiment and Return.

The free-energy surface (FES) of protein-ligand binding contains information useful for drug design. Here we show how to exploit a free-energy minimum of a protein-ligand complex identified by metadynamics simulations to design a new EphA2 antagonist with improved inhibitory potency.

Pushing the boundaries of vinylogous reactivity: catalytic enantioselective mukaiyama aldol reactions of highly unsaturated 2-silyloxyindoles.

The first example of catalytic, enantioselective hypervinylogous Mukaiyama aldol reaction (HVMAR) involving multiply unsaturated 2-silyloxyindoles is reported. The reaction utilizes a chiral Lewis base-catalyzed Lewis acid-mediated technology to deliver homoallylic 3-polyenylidene 2-oxindoles with extraordinary levels of regio-, enantio-, and geometrical selectivity. This work highlights a subtle yet decisive influence of the indole N-substituents on the propagation of the vinylogous reactivity space of the donor substrates up to ten bonds away from the origin of the vinylogy effect. Analysis of the (13) C NMR chemical shifts of the C-ω remote site within homologous silyloxyindole donors enabled rationalization of the results and easy qualitative prediction of the HVMAR reactivity/inertia toward a given aldehyde acceptor.

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.

Comparative Analysis of Virtual Screening Approaches in the Search for Novel EphA2 Receptor Antagonists.

The EphA2 receptor and its ephrin-A1 ligand form a key cell communication system, which has been found overexpressed in many cancer types and involved in tumor growth. Recent medicinal chemistry efforts have identified bile acid derivatives as low micromolar binders of the EphA2 receptor. However, these compounds suffer from poor physicochemical properties, hampering their use in vivo. The identification of compounds able to disrupt the EphA2-ephrin-A1 complex lacking the bile acid scaffold may lead to new pharmacological tools suitable for in vivo studies. To identify the most promising virtual screening (VS) protocol aimed at finding novel EphA2 antagonists, we investigated the ability of both ligand-based and structure-based approaches to retrieve known EphA2 antagonists from libraries of decoys with similar molecular properties. While ligand-based VSs were conducted using UniPR129 and ephrin-A1 ligand as reference structures, structure-based VSs were performed with Glide, using the X-ray structure of the EphA2 receptor/ephrin-A1 complex. A comparison of enrichment factors showed that ligand-based approaches outperformed the structure-based ones, suggesting ligand-based methods using the G-H loop of ephrin-A1 ligand as template as the most promising protocols to search for novel EphA2 antagonists.

Combining ligand- and structure-based approaches for the discovery of new inhibitors of the EPHA2-ephrin-A1 interaction.

The EPH receptor A2 (EPHA2) represents an attractive anticancer target. With the aim to identify novel EPHA2 receptor antagonists, a virtual screening campaign, combining shape-similarity and docking calculations, was conducted on a set of commercially available compounds. A combined score, taking into account both ligand- and structure-based results, was then used to identify the most promising candidates. Two compounds, selected among the best-ranked ones, were identified as EPHA2 receptor antagonists with micromolar affinity.

Synthesis and characterization of new bivalent agents as melatonin- and histamine H3-ligands.

Melatonin is an endogenous molecule involved in many pathophysiological processes. In addition to the control of circadian rhythms, its antioxidant and neuroprotective properties have been widely described. Thus far, different bivalent compounds composed by a melatonin molecule linked to another neuroprotective agent were synthesized and tested for their ability to block neurodegenerative processes in vitro and in vivo. To identify a novel class of potential neuroprotective compounds, we prepared a series of bivalent ligands, in which a prototypic melatonergic ligand is connected to an imidazole-based H3 receptor antagonist through a flexible linker. Four imidazolyl-alkyloxy-anilinoethylamide derivatives, characterized by linkers of different length, were synthesized and their binding affinity for human MT1, MT2 and H3 receptor subtypes was evaluated. Among the tested compounds, 14c and 14d, bearing a pentyl and a hexyl linker, respectively, were able to bind to all receptor subtypes at micromolar concentrations and represent the first bivalent melatonergic/histaminergic ligands reported so far. These preliminary results, based on binding affinity evaluation, pave the way for the future development of new dual-acting compounds targeting both melatonin and histamine receptors, which could represent promising therapeutic agents for the treatment of neurodegenerative pathologies.

Drug discovery: In silico dry data can bypass biological wet data?

The recent and extraordinary increase in computer power, along with the availability of efficient algorithms based on artificial intelligence, has prompted a large number of inexperienced scientists to challenge the complex and yet competitive world of drug discovery, by pretending to identify new hits through the sole use of computer aided drug design (CADD). Does the golden era of dry data run the risk of overshadowing the importance of wet data and, in doing so, forget that in silico and biological data need each other in successful preclinical drug discovery programmes?

Targeting metabolic adaptive responses induced by glucose starvation inhibits cell proliferation and enhances cell death in osimertinib-resistant non-small cell lung cancer (NSCLC) cell lines.

Osimertinib, a tyrosine kinase inhibitor targeting mutant EGFR, has received approval for initial treatment in patients with Non-Small Cell Lung Cancer (NSCLC). While effective in both first- and second-line treatments, patients eventually develop acquired resistance. Metabolic reprogramming represents a strategy through which cancer cells may resist and adapt to the selective pressure exerted by the drug. In the current study, we investigated the metabolic adaptations associated with osimertinib-resistance in NSCLC cells under low glucose culture conditions. We demonstrated that, unlike osimertinib-sensitive cells, osimertinib-resistant cells were able to survive under low glucose conditions by increasing the rate of glucose and glutamine uptake and by shifting towards mitochondrial metabolism. Inhibiting glucose/pyruvate contribution to mitochondrial respiration, glutamine deamination to glutamate, and oxidative phosphorylation decreased the proliferation and survival abilities of osimertinib-resistant cells to glucose starvation. Our findings underscore the remarkable adaptability of osimertinib-resistant NSCLC cells in a low glucose environment and highlight the pivotal role of mitochondrial metabolism in mediating this adaptation. Targeting the metabolic adaptive responses triggered by glucose shortage emerges as a promising strategy, effectively inhibiting cell proliferation and promoting cell death in osimertinib-resistant cells.

Discovery of a new 1-(phenylsulfonyl)-1H-indole derivative targeting the EphA2 receptor with antiproliferative activity on U251 glioblastoma cell line.

In our continuing effort devoted at developing agents targeting the EphA2 receptor by means of protein-protein interaction (PPI) inhibitors, we report here the design and synthesis of a new class of l-ß-homotryptophan conjugates of 3-ß-hydroxy-Δ5-cholenic acid bearing a set of arylsulfonyl substituents at the indole nitrogen atom. An extensive structure-activity relationship (SAR) analysis indicates that the presence of a bulky lipophilic moiety at the indole nitrogen is fundamental for improving potency on the EphA2 receptor, while abrogating activity on the EphB1-EphB3 receptor subtypes. A rational exploration, guided by the combined application of an experimental design on σp and π physicochemical descriptors and docking simulations, led to the discovery of UniPR1454, a 1-(4-(trifluoromethyl)phenyl)sulfonyl derivative acting as potent and competitive EphA2 antagonist able to inhibit ephrin-A1 dependent signals and to reduce proliferation of glioblastoma (U251) cell line at micromolar concentration.