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.

Targeting glucosylceramide synthase induces antiproliferative and proapoptotic effects in osimertinib-resistant NSCLC cell models.

The EGFR tyrosine kinase inhibitor osimertinib has been approved for the first-line treatment of EGFR-mutated Non-Small Cell Lung Cancer (NSCLC) patients. Despite its efficacy, patients develop resistance. Mechanisms of resistance are heterogeneous and not fully understood, and their characterization is essential to find new strategies to overcome resistance. Ceramides are well-known regulators of apoptosis and are converted into glucosylceramides (GlcCer) by glucosylceramide synthase (GCS). A higher content of GlcCers was observed in lung pleural effusions from NSCLC patients and their role in osimertinib-resistance has not been documented. The aim of this study was to determine the therapeutic potential of inhibiting GCS in NSCLC EGFR-mutant models resistant to osimertinib in vitro and in vivo. Lipidomic analysis showed a significant increase in the intracellular levels of glycosylceramides, including GlcCers in osimertinib resistant clones compared to sensitive cells. In resistant cells, the GCS inhibitor PDMP caused cell cycle arrest, inhibition of 2D and 3D cell proliferation, colony formation and migration capability, and apoptosis induction. The intratumoral injection of PDMP completely suppressed the growth of OR xenograft models. This study demonstrated that dysregulation of ceramide metabolism is involved in osimertinib-resistance and targeting GCS may be a promising therapeutic strategy for patients progressed to osimertinib.

Peptidyl nitroalkene inhibitors of main protease rationalized by computational and crystallographic investigations as antivirals against SARS-CoV-2.

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 exhibiting enzyme inhibitory activity against Mpro (Ki: 1-10 µM) 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. Structural analysis shows the binding mode of FGA146 and FGA147 to the active site of the protein. Furthermore, our results illustrate that peptidyl nitroalkenes are effective covalent reversible inhibitors of the Mpro and cathepsin L, and that inhibitors FGA145, FGA146 and FGA147 prevent infection against SARS-CoV-2.

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?

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.

A Pharmacological Investigation of Eph-Ephrin Antagonism in Prostate Cancer: UniPR1331 Efficacy Evidence.

The Eph kinases are the largest receptor tyrosine kinases (RTKs) family in humans. PC3 human prostate adenocarcinoma cells are a well-established model for studying Eph-ephrin pharmacology as they naturally express a high level of EphA2, a promising target for new cancer therapies. A pharmacological approach with agonists did not show significant efficacy on tumor growth in prostate orthotopic murine models, but reduced distal metastasis formation. In order to improve the comprehension of the pharmacological targeting of Eph receptors in prostate cancer, in the present work, we investigated the efficacy of Eph antagonism both in vitro and in vivo, using UniPR1331, a small orally bioavailable Eph-ephrin interaction inhibitor. UniPR1331 was able to inhibit PC3 cells' growth in vitro in a dose-dependent manner, affecting the cell cycle and inducing apoptosis. Moreover, UniPR1331 promoted the PC3 epithelial phenotype, downregulating epithelial mesenchymal transition (EMT) markers. As a consequence, UniPR1331 reduced in vitro PC3 migration, invasion, and vasculomimicry capabilities. The antitumor activity of UniPR1331 was confirmed in vivo when administered alone or in combination with cytotoxic drugs in PC3-xenograft mice. Our results demonstrated that Eph antagonism is a promising strategy for inhibiting prostate cancer growth, especially in combination with cytotoxic drugs.

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.

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.

Pharmacological characterization of second generation FXR agonists as effective EphA2 antagonists: A successful application of target hopping approach.

It is well demonstrated the key role of Eph-ephrin system, specifically of EphA2 receptor, in supporting tumor growth, invasion, metastasis and neovascularization. We previously identified FXR agonists as eligible antagonists of Eph-ephrin system. Herein we characterize new commercially available FXR (Farnesoid X Receptor) agonists as potential Eph ligands including Cilofexor, Nidufexor, Tropifexor, Turofexorate isopropyl and Vonafexor. Our exploration based on molecular modelling investigations and binding assays shows that Cilofexor binds specifically and reversibly to EphA2 receptor with a Ki value in the low micromolar range. Furthermore, Cilofexor interferes with the phosphorylation of EphA2 and the cell retraction and rounding in PC3 prostate cancer cells, both events depending on EphA2 activation. In conclusion, we can confirm that target hopping can be a successful approach to discover new moiety of protein-protein inhibitors.

Development of potent and selective FAAH inhibitors with improved drug-like properties as potential tools to treat neuroinflammatory conditions.

The neuroprotective performance against neuroinflammation of the endocannabinoid system (ECS) can be remarkably improved by indirect stimulation mediated by the pharmacological inhibition of the key ECS catabolic enzyme fatty acid amide hydrolase (FAAH). Based on our previous works and aiming to discover new selective FAAH inhibitors , we herein reported a new series of carbamate-based FAAH inhibitors (4a-t) which showed improved drug disposition properties compared to the previously reported analogues 2a-b. The introduction of ionizable functions allowed us to obtain new FAAH inhibitors of nanomolar potency characterized by good water solubility and chemical stability at physiological pH. Interesting structure-activity relationships (SARs), deeply analyzed by molecular docking and molecular dynamic (MD) simulations, were obtained. All the newly developed inhibitors showed an excellent selectivity profile evaluated against monoacylglycerol lipase and cannabinoid receptors. The reversible mechanism of action was determined by a rapid dilution assay. Absence of toxicity was confirmed in mouse fibroblasts NIH3T3 (for compounds 4e, 4g, 4n-o, and 4s) and in human astrocytes cell line 1321N1 (for compounds 4e, 4n, and 4s). The absence of undesired cardiac effects was also confirmed for compound 4n. Selected analogues (compounds 4e, 4g, 4n, and 4s) were able to reduce oxidative stress in 1321N1 astrocytes and exhibited notable neuroprotective effects when tested in an ex vivo model of neuroinflammation.

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.

In silico drug discovery of melatonin receptor ligands with therapeutic potential.

INTRODUCTION: The neurohormone melatonin (N-acetyl-5-methoxytryptamine) regulates circadian rhythms exerting a variety of effects in the central nervous system and in periphery. These activities are mainly mediated by activation of MT1 and MT2 GPCRs. MT1/MT2 agonist compounds are used clinically for insomnia, depression, and circadian rhythm disturbances. AREA COVERED: The following review describes the design strategies that have led to the identification of melatonin receptor ligands, guided by in silico approaches and molecular modeling. Initial ligand-based design, mainly relying on pharmacophore modeling and 3D-QSAR, has been flanked by structure-based virtual screening, given the recent availability of MT1 and MT2 crystal structures. Receptor ligands with different activity profiles, agonist/antagonist and subtype-selective compounds, are available. EXPERT OPINION: An insight on the pharmacological characterization and therapeutic perspectives for relevant ligands is provided. In silico drug discovery has been instrumental in the design of novel ligands targeting melatonin receptors. Ligand-based approaches has led to the construction of a solid framework defining structure-activity relationships to obtain compounds with a tailored pharmacological profile. Structure-based techniques could integrate previous knowledge and provide compounds with novel chemotypes and pharmacological activity as drug candidates for disease conditions in which melatonin receptor ligands are currently being investigated, including cancer and pain.

Protein-Protein Interaction Inhibitors Targeting the Eph-Ephrin System with a Focus on Amino Acid Conjugates of Bile Acids.

The role of the Eph-ephrin system in the etiology of pathological conditions has been consolidated throughout the years. In this context, approaches directed against this signaling system, intended to modulate its activity, can be strategic therapeutic opportunities. Currently, the most promising class of compounds able to interfere with the Eph receptor-ephrin protein interaction is composed of synthetic derivatives of bile acids. In the present review, we summarize the progresses achieved, in terms of chemical expansions and structure-activity relationships, both in the steroidal core and the terminal carboxylic acid group, along with the pharmacological characterization for the most promising Eph-ephrin antagonists in in vivo settings.

Impact of Warhead Modulations on the Covalent Inhibition of SARS-CoV-2 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.

Cholenic acid derivative UniPR1331 impairs tumor angiogenesis via blockade of VEGF/VEGFR2 in addition to Eph/ephrin.

Angiogenesis, the formation of new blood vessels from preexisting ones, is crucial for tumor growth and metastatization, and is considered a promising therapeutic target. Unfortunately, drugs directed against a specific proangiogenic growth factor or receptor turned out to be of limited benefit for oncology patients, likely due to the high biochemical redundancy of the neovascularization process. In this scenario, multitarget compounds that are able to simultaneously tackle different proangiogenic pathways are eagerly awaited. UniPR1331 is a 3ß-hydroxy-Δ5-cholenic acid derivative, which is already known to inhibit Eph-ephrin interaction. Here, we employed an analysis pipeline consisting of molecular modeling and simulation, surface plasmon resonance spectrometry, biochemical assays, and endothelial cell models to demonstrate that UniPR1331 directly interacts with the vascular endothelial growth factor receptor 2 (VEGFR2) too. The binding of UniPR1331 to VEGFR2 prevents its interaction with the natural ligand vascular endothelial growth factor and subsequent autophosphorylation, signal transduction, and in vitro proangiogenic activation of endothelial cells. In vivo, UniPR1331 inhibits tumor cell-driven angiogenesis in zebrafish. Taken together, these data shed light on the pleiotropic pharmacological effect of UniPR1331, and point to Δ5-cholenic acid as a promising molecular scaffold for the development of multitarget antiangiogenic compounds.

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.

Discovery of SARS-CoV-2 Mpro peptide inhibitors from modelling substrate and ligand binding.

The main protease (Mpro) of SARS-CoV-2 is central to viral maturation and is a promising drug target, but little is known about structural aspects of how it binds to its 11 natural cleavage sites. We used biophysical and crystallographic data and an array of biomolecular simulation techniques, including automated docking, molecular dynamics (MD) and interactive MD in virtual reality, QM/MM, and linear-scaling DFT, to investigate the molecular features underlying recognition of the natural Mpro substrates. We extensively analysed the subsite interactions of modelled 11-residue cleavage site peptides, crystallographic ligands, and docked COVID Moonshot-designed covalent inhibitors. Our modelling studies reveal remarkable consistency in the hydrogen bonding patterns of the natural Mpro substrates, particularly on the N-terminal side of the scissile bond. They highlight the critical role of interactions beyond the immediate active site in recognition and catalysis, in particular plasticity at the S2 site. Building on our initial Mpro-substrate models, we used predictive saturation variation scanning (PreSaVS) to design peptides with improved affinity. Non-denaturing mass spectrometry and other biophysical analyses confirm these new and effective 'peptibitors' inhibit Mpro competitively. Our combined results provide new insights and highlight opportunities for the development of Mpro inhibitors as anti-COVID-19 drugs.

A sulfonyl fluoride derivative inhibits EGFRL858R/T790M/C797S by covalent modification of the catalytic lysine.

The emergence of the C797S mutation in EGFR is a frequent mechanism of resistance to osimertinib in the treatment of non-small cell lung cancer (NSCLC). In the present work, we report the design, synthesis and biochemical characterization of UPR1444 (compound 11), a new sulfonyl fluoride derivative which potently and irreversibly inhibits EGFRL858R/T790M/C797S through the formation of a sulfonamide bond with the catalytic residue Lys745. Enzymatic assays show that compound 11 displayed an inhibitory activity on EGFRWT comparable to that of osimertinib, and it resulted more selective than the sulfonyl fluoride probe XO44, recently reported to inhibit a significant part of the kinome. Neither compound 11 nor XO44 inhibited EGFRdel19/T790M/C797S triple mutant. When tested in Ba/F3 cells expressing EGFRL858R/T790M/C797S, compound 11 resulted significantly more potent than osimertinib at inhibiting both EGFR autophosphorylation and proliferation, even if the inhibition of EGFR autophosphorylation by compound 11 in Ba/F3 cells was not long lasting.

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.