Discovery of novel FGF trap small molecules endowed with anti-myeloma activity.

Fibroblast growth factors (FGFs) act as proangiogenic and mitogenic cytokines in several cancers, including multiple myeloma (MM). Indeed, corrupted FGF autocrine and paracrine secretion induces an aberrant activation of the FGF receptor (FGFR) signaling sustaining cancer cell spreading and resistance to pharmacological treatments. Thus, FGF traps may represent a promising anti-cancer strategy to hamper the ligand-dependent activation of the FGF/FGFR system. We previously identified NSC12 as the first orally available small molecule FGF trap able to inhibit the growth and progression of several FGF-dependent tumor models. NSC12 is a pregnenolone derivative carrying a 1,1-bis-trifluoromethyl-1,3-propanediol chain in position 17 of the steroid nucleus. Investigation of structure-activity relationships (SARs) provided more potent and specific NSC12 steroid derivatives and highlighted that the C17-side chain is pivotal for the FGF trap activity. Here, a scaffold hopping approach allowed to obtain two FGF trap compounds (22 and 57) devoid of the steroid nucleus and able to efficiently bind FGF2 and to inhibit FGFR activation in MM cells. Accordingly, these compounds exert a potent anti-tumor activity on MM cell lines both in vitro and in vivo and on MM patient-derived primary cells, strongly affecting the survival of both proteasome-inhibitor sensitive and resistant MM cells. These results propose a new therapeutic option for relapsed/refractory MM patients and set the bases for the development of novel FGF traps prone to chemical diversification to be used in the clinic for the treatment of those tumors in which the FGF/FGFR system plays a pivotal role, including MM.

Melatonin receptor structure and signaling.

Melatonin (5-methoxy-N-acetyltryptamine) binds with high affinity and specificity to membrane receptors. Several receptor subtypes exist in different species, of which the mammalian MT1 and MT2 receptors are the best-characterized. They are members of the G protein-coupled receptor superfamily, preferentially coupling to Gi/o proteins but also to other G proteins in a cell-context-depending manner. In this review, experts on melatonin receptors will summarize the current state of the field. We briefly report on the discovery and classification of melatonin receptors, then focus on the molecular structure of human MT1 and MT2 receptors and highlight the importance of molecular simulations to identify new ligands and to understand the structural dynamics of these receptors. We then describe the state-of-the-art of the intracellular signaling pathways activated by melatonin receptors and their complexes. Brief statements on the molecular toolbox available for melatonin receptor studies and future perspectives will round-up this review.

Binding and unbinding of potent melatonin receptor ligands: Mechanistic simulations and experimental evidence.

The labeled ligand commonly employed in competition binding studies for melatonin receptor ligands, 2-[125I]iodomelatonin, showed slow dissociation with different half-lives at the two receptor subtypes. This may affect the operational measures of affinity constants, which at short incubation times could not be obtained in equilibrium conditions, and structure-activity relationships, as the Ki values of tested ligands could depend on either interaction at the binding site or the dissociation path. To address these issues, the kinetic and saturation binding parameters of 2-[125I]iodomelatonin as well as the competition constants for a series of representative ligands were measured at a short (2 h) and a long (20 h) incubation time. Concurrently, we simulated by molecular modeling the dissociation path of 2-iodomelatonin from MT1 and MT2 receptors and investigated the role of interactions at the binding site on the stereoselectivity observed for the enantiomers of the subtype-selective ligand UCM1014. We found that equilibrium conditions for 2-[125I]iodomelatonin binding can be reached only with long incubation times, particularly for the MT2 receptor subtype, for which a time of 20 h approximates this condition. On the other hand, measured Ki values for a set of ligands including agonists, antagonists, nonselective, and subtype-selective compounds were not significantly affected by the length of incubation, suggesting that structure-activity relationships based on data collected at shorter time reflect different interactions at the binding site. Molecular modeling simulations evidenced that the slower dissociation of 2-iodomelatonin from the MT2 receptor can be related to the restricted mobility of a gatekeeper tyrosine along a lipophilic path from the binding site to the membrane bilayer. The enantiomers of the potent, MT2-selective agonist UCM1014 were separately synthesized and tested. Molecular dynamics simulations of the receptor-ligand complexes provided an explanation for their stereoselectivity as due to the preference shown by the eutomer at the binding site for the most abundant axial conformation adopted by the ligand in solution. These results suggest that, despite the slow-binding kinetics occurring for the labeled ligand, affinity measures at shorter incubation times give robust results consistent with known structure-activity relationships and with interactions taken at the receptor binding site.

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.

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.

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.

New Coumarin Derivatives as Cholinergic and Cannabinoid System Modulators.

In the last years, the connection between the endocannabinoid system (eCS) and neuroprotection has been discovered, and evidence indicates that eCS signaling is involved in the regulation of cognitive processes and in the pathophysiology of Alzheimer's disease (AD). Accordingly, pharmacotherapy targeting eCS could represent a valuable contribution in fighting a multifaceted disease such as AD, opening a new perspective for the development of active agents with multitarget potential. In this paper, a series of coumarin-based carbamic and amide derivatives were designed and synthesized as multipotent compounds acting on cholinergic system and eCS-related targets. Indeed, they were tested with appropriate enzymatic assays on acetyl and butyryl-cholinesterases and on fatty acid amide hydrolase (FAAH), and also evaluated as cannabinoid receptor (CB1 and CB2) ligands. Moreover, their ability to reduce the self-aggregation of beta amyloid protein (Aß42) was assessed. Compounds 2 and 3, bearing a carbamate function, emerged as promising inhibitors of hAChE, hBuChE, FAAH and Aß42 self-aggregation, albeit with moderate potencies, while the amide 6 also appears a promising CB1/CB2 receptors ligand. These data prove for the new compounds an encouraging multitarget profile, deserving further evaluation.

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.

The Control of Heparanase Through the Use of Small Molecules.

Despite the enormous progress made in recent years with antibodies, vaccines, antisense oligonucleotides, etc., the so-called "biological" approaches for tackling the control of various diseases, medicinal chemistry remains a bulwark to refer to for the development of new drugs. Also in the case of heparanase, medicinal chemistry has always been in the forefront to identify new inhibitors, through modification of natural macromolecules, e.g., sulfated polysaccharides like heparin, or of natural compounds isolated from bacteria or plants, or through rational design. In this chapter, the reader will find a detailed description of the most relevant small-molecule heparanase inhibitors reported so far in the scientific literature and in patent applications, with mention to the design strategy and to structure-activity relationships. Starting from heparanase inhibitors of natural origin and the attempts to improve their potency and selectivity, the reader will be guided through the major chemical classes of synthetic inhibitors, with representation of the structure of the most relevant compounds. The last paragraph is dedicated to a brief description of inhibitors that have reached clinical trials, highlighting their structure, mechanism, and improved derivatives.

Chiral Recognition of Flexible Melatonin Receptor Ligands Induced by Conformational Equilibria.

N-anilinoethylamides are a class of melatoninergic agents with the aniline portion mimicking the indole ring of the natural ligand and the ethylamide chain reproducing that of melatonin. The simplest compound in this class, N-{2-[(3-methoxyphenyl)methylamino]ethyl}acetamide (UCM793), has nanomolar binding affinity for MT1 and MT2 membrane receptors. To explore the effect of chain conformation on receptor binding, a methyl group was inserted on the methylene alpha or beta to the amide nitrogen and conformational equilibria were investigated by NMR spectroscopy and molecular dynamics simulations. Receptor affinity was conserved only for the beta-methyl derivative, which also showed significant stereoselectivity, with the (S) enantiomer being the eutomer. Molecular dynamics simulations, validated by NMR spectroscopy, showed that the beta-methyl group affects the conformational preferences of the ethylamide chain. Docking into the receptor crystal structure provides a rationale for the observed chiral recognition, suggesting that the (S)-beta-methyl group favors the conformation that better fits the receptor binding site.

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.

Kinetic analysis and molecular modeling of the inhibition mechanism of roneparstat (SST0001) on human heparanase.

Heparanase is a ß-d-glucuronidase which cleaves heparan sulfate chains in the extracellular matrix and on cellular membranes. A dysregulated heparanase activity is intimately associated with cell invasion, tumor metastasis and angiogenesis, making heparanase an attractive target for the development of anticancer therapies. SST0001 (roneparstat; Sigma-Tau Research Switzerland S.A.) is a non-anticoagulant 100% N-acetylated and glycol-split heparin acting as a potent heparanase inhibitor, currently in phase I in advanced multiple myeloma. Herein, the kinetics of heparanase inhibition by roneparstat is reported. The analysis of dose-inhibition curves confirmed the high potency of roneparstat (IC50 ≈ 3 nM) and showed, at higher concentrations, a Hill coefficient consistent with the engagement of two molecules of inhibitor. A homology model of human heparanase GS3 construct was built and used for docking experiments with inhibitor fragments. The model has high structural similarity with the recently reported crystal structure of human heparanase. Different interaction schemes are proposed, which support the hypothesis of a complex binding mechanism involving the recruitment of one or multiple roneparstat chains, depending on its concentration. In particular, docking solutions were obtained in which (i) a single roneparstat molecule interacts with both heparin-binding domains (HBDs) of heparanase or (ii) two fragments of roneparstat interact with either HBD-1 or HBD-2, consistent with the possibility of different inhibitor:enzyme binding stoichiometries. This study provides unique insights into the mode of action of roneparstat as well as clues of its interaction with heparanase at a molecular level, which could be exploited to design novel potential inhibitor molecules.

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.

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.

Synthesis, structure and inhibitory activity of a stereoisomer of oseltamivir carboxylate.

A stereodivergent plan is presented leading to all eight stereoisomers of oseltamivir carboxylate (OC). Key chemical manoeuvers are (1) a three-component vinylogous Mukaiyama-Mannich reaction, which sets the whole carbon skeleton and heteroatom substituents, and (2) an intramolecular, silylative Mukaiyama aldol reaction, which creates the targeted carbocycle. The viability of the plan was demonstrated by the first total synthesis of 4-epi-oseltamivir carboxylate (6), accessed in 15 steps from glyceraldehyde, o-anisidine and pyrrole siloxydiene precursors. Compound 6 inhibits influenza A virus strains H1N1 and H3N2 at the µM level, about 150 000-fold less than the OC reference, testifying that the stereodisposition of the C4 acetamido function is key for enzyme recognition. Guided by in-depth structural evaluation including NMR solution studies, molecular mechanics simulations, docking analyses and X-ray crystallography, rationalization of the biological verdict was established.

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.

Social stress-induced depressive-like symptoms and changes in gut microbial and lipidomic profiles are prevented by pharmacological inhibition of FAAH activity in male rats.

Pharmacological inhibition of fatty acid amide hydrolase (FAAH) activity has antidepressant-like effects in preclinical models of stress. In this study, we investigated whether the antidepressant-like effects of FAAH inhibition are associated with corresponding changes in gut microbial and lipidomic profiles, which are emerging as critical components in the pathophysiology of depression. Adult male Wistar rats experienced five weeks of repeated social defeat or control procedure and were treated with the FAAH inhibitor URB694 (0.3 mg/kg/day, i.p.) or vehicle starting from the third week. Repeated social defeat induced the emergence of depressive-like behavioral (sucrose preference reduction and passive coping behaviors in the forced swim test) and neuroendocrine (increased corticosterone levels) changes, which were prevented by URB694 treatment. Repeated social defeat also provoked a significant variation in gut microbiota (changes in the relative abundance of 14 bacterial taxa) and lipidic (e.g., glycerophospholipids) composition. These stress-induced changes were prevented by URB694 treatment. These findings indicate that inhibition of FAAH activity with URB694 blocks the co-occurrence of depressive-like behavioral and neuroendocrine changes and alterations in gut microbial and lipid composition in rats exposed to repeated social defeat. In conclusion, these results suggest that the gut microbiota-lipid crosstalk may represent a novel biological target for FAAH inhibitors to enhance stress resilience.

Long-lasting inhibition of EGFR autophosphorylation in A549 tumor cells by intracellular accumulation of non-covalent inhibitors.

In the present study, a small set of reversible or irreversible 4-anilinoquinazoline EGFR inhibitors was tested in A549 cells at early (1h) and late (8h) time points after inhibitor removal from culture medium. A combination of assays was employed to explain the observed long-lasting inhibition of EGFR autophosphorylation. We found that EGFR inhibition at 8h can be due, besides to the covalent interaction of the inhibitor with Cys797, as for PD168393 (2) and its prodrug 4, to the intracellular accumulation of non-covalent inhibitors by means of an active cell uptake, as for 5 and 6. Compounds 5-6 showed similar potency and duration of inhibition of EGFR autophosphorylation as the covalent inhibitor 2, while being devoid of reactive groups forming covalent bonds with protein thiols.