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.

Halting the FGF/FGFR axis leads to antitumor activity in Waldenström macroglobulinemia by silencing MYD88.

The human fibroblast growth factor/fibroblast growth factor receptor (FGF/FGFR) axis deregulation is largely involved in supporting the pathogenesis of hematologic malignancies, including Waldenström macroglobulinemia (WM). WM is still an incurable disease, and patients succumb because of disease progression. Therefore, novel therapeutics designed to specifically target deregulated signaling pathways in WM are required. We aimed to investigate the role of FGF/FGFR system blockade in WM by using a pan-FGF trap molecule (NSC12). Wide-transcriptome profiling confirmed inhibition of FGFR signaling in NSC12-treated WM cells; unveiling a significant inhibition of MYD88 was also confirmed at the protein level. Importantly, the NSC12-dependent silencing of MYD88 was functionally active, as it led to inhibition of MYD88-driven pathways, such as BTK and SYK, as well as the MYD88-downstream target HCK. Of note, both canonical and noncanonical NF-κB cascades were downregulated in WM cells upon NSC12 treatment. Functional sequelae exerted by NSC12 in WM cells were studied, demonstrating significant inhibition of WM cell growth, induction of WM cell apoptosis, halting MAPK, JAK/STAT3, and PI3K-Akt pathways. Importantly, NSC12 exerted an anti-WM effect even in the presence of bone marrow microenvironment, both in vitro and in vivo. Our studies provide the evidence for using NSC12 as a specific FGF/FGFR system inhibitor, thus representing a novel therapeutic strategy in WM.

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.

Systematic Modification of the Substitution Pattern of the 7-Hydroxy-5-oxopyrazolo[4,3-b]pyridine-6-carboxamide Scaffold Enabled the Discovery of New Ligands with High Affinity and Selectivity for the Cannabinoid Type 2 Receptor.

Selective ligands of the CB2 receptor are receiving considerable attention due to their potential as therapeutic agents for a variety of diseases. Recently, 7-hydroxy-5-oxopyrazolo[4,3-b]pyridine-6-carboxamide derivatives were shown to act at the CB2 receptor either as agonists or as inverse agonists/antagonists in vitro and to have anti-osteoarthritic activity in vivo. In this article, we report the synthesis, pharmacological profile, and molecular modeling of a series of twenty-three new 7-hydroxy-5-oxopyrazolo[4,3-b]pyridine-6-carboxamides with the aim of further developing this new class of selective CB2 ligands. In addition to these compounds, seven other analogs that had been previously synthesized were included in this study to better define the structure-activity relationship (SAR). Ten of the new compounds studied were found to be potent and selective ligands of the CB2 receptor, with Ki values ranging from 48.46 to 0.45 nM and CB1/CB2 selectivity indices (SI) ranging from >206 to >4739. In particular, compounds 54 and 55 were found to be high-affinity CB2 inverse agonists that were not active at all at the CB1 receptor, whereas 57 acted as an agonist. The functional activity profile of the compounds within this structural class depends mainly on the substitution pattern of the pyrazole ring.

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.

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.

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.

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.

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.

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.

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.