Computational Methods for the Discovery and Optimization of TAAR1 and TAAR5 Ligands.

G-protein-coupled receptors (GPCRs) represent a family of druggable targets when treating several diseases and continue to be a leading part of the drug discovery process. Trace amine-associated receptors (TAARs) are GPCRs involved in many physiological functions with TAAR1 having important roles within the central nervous system (CNS). By using homology modeling methods, the responsiveness of TAAR1 to endogenous and synthetic ligands has been explored. In addition, the discovery of different chemo-types as selective murine and/or human TAAR1 ligands has helped in the understanding of the species-specificity preferences. The availability of TAAR1-ligand complexes sheds light on how different ligands bind TAAR1. TAAR5 is considered an olfactory receptor but has specific involvement in some brain functions. In this case, the drug discovery effort has been limited. Here, we review the successful computational efforts developed in the search for novel TAAR1 and TAAR5 ligands. A specific focus on applying structure-based and/or ligand-based methods has been done. We also give a perspective of the experimental data available to guide the future drug design of new ligands, probing species-specificity preferences towards more selective ligands. Hints for applying repositioning approaches are also discussed.

Discovery of Novel Thiazole-Based SIRT2 Inhibitors as Anticancer Agents: Molecular Modeling, Chemical Synthesis and Biological Assays.

The search and development of effective sirtuin small molecule inhibitors (SIRTIs) continues to draw great attention due to their wide range of pharmacological applications. Based on SIRTs' involvement in different biological pathways, their ligands were investigated for many diseases, such as cancer, neurodegenerative disorders, diabetes, cardiovascular diseases and autoimmune diseases. The elucidation of a substantial number of SIRT2-ligand complexes is steering the identification of novel and more selective modulators. Among them, SIRT2 in the presence of the SirReal2 analog series was the most studied. On this basis, we recently reported structure-based analyses leading to the discovery of thiazole-based compounds acting as SIRT2 inhibitors (T1, SIRT2 IC50 = 17.3 µM). Herein, ligand-based approaches followed by molecular docking simulations allowed us to evaluate in silico a novel small series of thiazoles (3a-3d and 5a, 5d) as putative SIRT2 inhibitors. Results from the computational studies revealed comparable molecular interaction fields (MIFs) and docking positionings of most of these compounds with respect to reference SIRT2Is. Biochemical and biological assays validated this study and pointed to compound 5a (SIRT2 IC50 = 9.0 µM) as the most interesting SIRT2I that was worthy of further development as an anticancer agent.

Discovery of a Novel Chemo-Type for TAAR1 Agonism via Molecular Modeling.

The search for novel effective TAAR1 ligands continues to draw great attention due to the wide range of pharmacological applications related to TAAR1 targeting. Herein, molecular docking studies of known TAAR1 ligands, characterized by an oxazoline core, have been performed in order to identify novel promising chemo-types for the discovery of more active TAAR1 agonists. In particular, the oxazoline-based compound S18616 has been taken as a reference compound for the computational study, leading to the development of quite flat and conformationally locked ligands. The choice of a "Y-shape" conformation was suggested for the design of TAAR1 ligands, interacting with the protein cavity delimited by ASP103 and aromatic residues such as PHE186, PHE195, PHE268, and PHE267. The obtained results allowed us to preliminary in silico screen an in-house series of pyrimidinone-benzimidazoles (1a-10a) as a novel scaffold to target TAAR1. Combined ligand-based (LBCM) and structure based (SBCM) computational methods suggested the biological evaluation of compounds 1a-10a, leading to the identification of derivatives 1a-3a (hTAAR1 EC50 = 526.3-657.4 nM) as promising novel TAAR1 agonists.

Novel thienocycloalkylpyridazinones as useful scaffolds for acetylcholinesterase inhibition and serotonin 5-HT6 receptor interaction.

A library of eighteen thienocycloalkylpyridazinones was synthesized for human acetylcholinesterase (hAChE) and butyrylcholinesterase (hBChE) inhibition and serotonin 5-HT6 receptor subtype interaction by following a multitarget-directed ligand approach (MTDL), as a suitable strategy for treatment of Alzheimer's disease (AD). The novel compounds featured a tricyclic scaffold, namely thieno[3,2-h]cinnolinone, thienocyclopentapyridazinone and thienocycloheptapyridazinone, connected through alkyl chains of variable length to proper amine moieties, most often represented by N-benzylpiperazine or 1-(phenylsulfonyl)-4-(piperazin-1-ylmethyl)-1H-indole as structural elements addressing AChE and 5-HT6 interaction, respectively. Our study highlighted the versatility of thienocycloalkylpyridazinones as useful architectures for AChE interaction, with several N-benzylpiperazine-based analogues emerging as potent and selective hAChE inhibitors with IC50 in the 0.17-1.23 µM range, exhibiting low to poor activity for hBChE (IC50 = 4.13-9.70 µM). The introduction of 5-HT6 structural moiety phenylsulfonylindole in place of N-benzylpiperazine, in tandem with a pentamethylene linker, gave potent 5-HT6 thieno[3,2-h]cinnolinone and thienocyclopentapyridazinone-based ligands both displaying hAChE inhibition in the low micromolar range and unappreciable activity towards hBChE. While docking studies provided a rational structural explanation for AChE/BChE enzyme and 5-HT6 receptor interaction, in silico prediction of ADME properties of tested compounds suggested further optimization for development of such compounds in the field of MTDL for AD.

Virtual Screening Combined with Enzymatic Assays to Guide the Discovery of Novel SIRT2 Inhibitors.

Sirtuin isoform 2 (SIRT2) is one of the seven sirtuin isoforms present in humans, being classified as class III histone deacetylases (HDACs). Based on the high sequence similarity among SIRTs, the identification of isoform selective modulators represents a challenging task, especially for the high conservation observed in the catalytic site. Efforts in rationalizing selectivity based on key residues belonging to the SIRT2 enzyme were accompanied in 2015 by the publication of the first X-ray crystallographic structure of the potent and selective SIRT2 inhibitor SirReal2. The subsequent studies led to different experimental data regarding this protein in complex with further different chemo-types as SIRT2 inhibitors. Herein, we reported preliminary Structure-Based Virtual Screening (SBVS) studies using a commercially available library of compounds to identify novel scaffolds for the design of new SIRT2 inhibitors. Biochemical assays involving five selected compounds allowed us to highlight the most effective chemical features supporting the observed SIRT2 inhibitory ability. This information guided the following in silico evaluation and in vitro testing of further compounds from in-house libraries of pyrazolo-pyrimidine derivatives towards novel SIRT2 inhibitors (1-5). The final results indicated the effectiveness of this scaffold for the design of promising and selective SIRT2 inhibitors, featuring the highest inhibition among the tested compounds, and validating the applied strategy.

Discovery of Novel Trace Amine-Associated Receptor 5 (TAAR5) Antagonists Using a Deep Convolutional Neural Network.

Trace amine-associated receptor 5 (TAAR5) is a G protein-coupled receptor that belongs to the TAARs family (TAAR1-TAAR9). TAAR5 is expressed in the olfactory epithelium and is responsible for sensing 3-methylamine (TMA). However, recent studies showed that TAAR5 is also expressed in the limbic brain regions and is involved in the regulation of emotional behaviour and adult neurogenesis, suggesting that TAAR5 antagonism may represent a novel therapeutic strategy for anxiety and depression. We used the AtomNet® model, the first deep learning neural network for structure-based drug discovery, to identify putative TAAR5 ligands and tested them in an in vitro BRET assay. We found two mTAAR5 antagonists with low to submicromolar activity that are able to inhibit the cAMP production induced by TMA. Moreover, these two compounds also inhibited the mTAAR5 downstream signalling, such as the phosphorylation of CREB and ERK. These two hits exhibit drug-like properties and could be used to further develop more potent TAAR5 ligands with putative anxiolytic and antidepressant activity.

Virtual Screening in the Identification of Sirtuins' Activity Modulators.

Sirtuins are NAD+-dependent deac(et)ylases with different subcellular localization. The sirtuins' family is composed of seven members, named SIRT-1 to SIRT-7. Their substrates include histones and also an increasing number of different proteins. Sirtuins regulate a wide range of different processes, ranging from transcription to metabolism to genome stability. Thus, their dysregulation has been related to the pathogenesis of different diseases. In this review, we discussed the pharmacological approaches based on sirtuins' modulators (both inhibitors and activators) that have been attempted in in vitro and/or in in vivo experimental settings, to highlight the therapeutic potential of targeting one/more specific sirtuin isoform(s) in cancer, neurodegenerative disorders and type 2 diabetes. Extensive research has already been performed to identify SIRT-1 and -2 modulators, while compounds targeting the other sirtuins have been less studied so far. Beside sections dedicated to each sirtuin, in the present review we also included sections dedicated to pan-sirtuins' and to parasitic sirtuins' modulators. A special focus is dedicated to the sirtuins' modulators identified by the use of virtual screening.

Molecular Docking and QSAR Studies as Computational Tools Exploring the Rescue Ability of F508del CFTR Correctors.

Cystic fibrosis (CF) is the autosomal recessive disorder most recurrent in Caucasian populations. Different mutations involving the cystic fibrosis transmembrane regulator protein (CFTR) gene, which encodes the CFTR channel, are involved in CF. A number of life-prolonging therapies have been conceived and deeply investigated to combat this disease. Among them, the administration of the so-called CFTR modulators, such as correctors and potentiators, have led to quite beneficial effects. Recently, based on QSAR (quantitative structure activity relationship) studies, we reported the rational design and synthesis of compound 2, an aminoarylthiazole-VX-809 hybrid derivative exhibiting promising F508del-CFTR corrector ability. Herein, we explored the docking mode of the prototype VX-809 as well as of the aforementioned correctors in order to derive useful guidelines for the rational design of further analogues. In addition, we refined our previous QSAR analysis taking into account our first series of in-house hybrids. This allowed us to optimize the QSAR model based on the chemical structure and the potency profile of hybrids as F508del-CFTR correctors, identifying novel molecular descriptors explaining the SAR of the dataset. This study is expected to speed up the discovery process of novel potent CFTR modulators.

Recent Strategic Advances in CFTR Drug Discovery: An Overview.

Cystic fibrosis transmembrane conductance regulator (CFTR)-rescuing drugs have already transformed cystic fibrosis (CF) from a fatal disease to a treatable chronic condition. However, new-generation drugs able to bind CFTR with higher specificity/affinity and to exert stronger therapeutic benefits and fewer side effects are still awaited. Computational methods and biosensors have become indispensable tools in the process of drug discovery for many important human pathologies. Instead, they have been used only piecemeal in CF so far, calling for their appropriate integration with well-tried CF biochemical and cell-based models to speed up the discovery of new CFTR-rescuing drugs. This review will give an overview of the available structures and computational models of CFTR and of the biosensors, biochemical and cell-based assays already used in CF-oriented studies. It will also give the reader some insights about how to integrate these tools as to improve the efficiency of the drug discovery process targeted to CFTR.

Synthesis and Biological Evaluation of Novel (thio)semicarbazone-Based Benzimidazoles as Antiviral Agents against Human Respiratory Viruses.

Respiratory RNA viruses are responsible for recurrent acute respiratory illnesses that still represent a major medical need. Previously we developed a large variety of benzimidazole derivatives able to inhibit these viruses. Herein, two series of (thio)semicarbazone- and hydrazone-based benzimidazoles have been explored, by derivatizing 5-acetyl benzimidazoles previously reported by us, thereby evaluating the influence of the modification on the antiviral activity. Compounds 6, 8, 16 and 17, bearing the 5-(thio)semicarbazone and 5-hydrazone functionalities in combination with the 2-benzyl ring on the benzimidazole core structure, acted as dual inhibitors of influenza A virus and human coronavirus. For respiratory syncytial virus (RSV), activity is limited to the 5-thiosemicarbazone (25) and 5-hydrazone (22) compounds carrying the 2-[(benzotriazol-1/2-yl)methyl]benzimidazole scaffold. These molecules proved to be the most effective antiviral agents, able to reach the potency profile of the licensed drug ribavirin. The molecular docking analysis explained the SAR of these compounds around their binding mode to the target RSV F protein, revealing the key contacts for further assessment. The herein-investigated benzimidazole-based derivatives may represent valuable hit compounds, deserving subsequent structural improvements towards more efficient antiviral agents for the treatment of pathologies caused by these human respiratory viruses.

Novel sulfenamides and sulfonamides based on pyridazinone and pyridazine scaffolds as CB1 receptor ligand antagonists.

A series of sulfenamide and sulfonamide derivatives was synthesized and evaluated for the affinity at CB1 and CB2 receptors. The N-bornyl-S-(5,6-di-p-tolylpyridazin-3-yl)-sulfenamide, compound 11, displayed good affinity and high selectivity for CB1 receptors (Ki values of 44.6 nM for CB1 receptors and >40 µM for CB2 receptors, respectively). The N-isopinocampheyl-sulfenamide 12 and its sulfonamide analogue 22 showed similar selectivity for CB1 receptors with Ki values of 75.5 and 73.2 nM, respectively. These novel compounds behave as antagonists/inverse agonists at CB1 receptor in the [35S]-GTPγS binding assays, and none showed adequate predictive blood-brain barrier permeation, exhibiting low estimated LD50. However, testing compound 12 in a supraspinal analgesic test (hot-plate) revealed that it was as effective as the classic CB1 receptor antagonist rimonabant, in reversing the analgesic effect of a cannabinoid agonist.

Speeding Up the Identification of Cystic Fibrosis Transmembrane Conductance Regulator-Targeted Drugs: An Approach Based on Bioinformatics Strategies and Surface Plasmon Resonance.

Cystic fibrosis (CF) is mainly caused by the deletion of Phe 508 (ΔF508) in the cystic fibrosis transmembrane conductance regulator (CFTR) protein that is thus withheld in the endoplasmic reticulum and rapidly degraded by the ubiquitin/proteasome system. New drugs able to rescue ΔF508-CFTR trafficking are eagerly awaited. An integrated bioinformatics and surface plasmon resonance (SPR) approach was here applied to investigate the rescue mechanism(s) of a series of CFTR-ligands including VX809, VX770 and some aminoarylthiazole derivatives (AAT). Computational studies tentatively identified a large binding pocket in the ΔF508-CFTR nucleotide binding domain-1 (NBD1) and predicted all the tested compounds to bind to three sub-regions of this main pocket. Noticeably, the known CFTR chaperone keratin-8 (K8) seems to interact with some residues located in one of these sub-pockets, potentially interfering with the binding of some ligands. SPR results corroborated all these computational findings. Moreover, for all the considered ligands, a statistically significant correlation was determined between their binding capability to ΔF508-NBD1 measured by SPR and the pockets availability measured by computational studies. Taken together, these results demonstrate a strong agreement between the in silico prediction and the SPR-generated binding data, suggesting a path to speed up the identification of new drugs for the treatment of cystic fibrosis.

Exhaustive CoMFA and CoMSIA analyses around different chemical entities: a ligand-based study exploring the affinity and selectivity profiles of 5-HT1A ligands.

The 5-hydroxytryptamine (5-HT1A) receptors represent an attractive target in drug discovery. In particular, 5-HT1A agonists and partial agonists are deeply investigated for their potential role in the treatment of anxiety, depression, ischaemic brain disorder and more recently, of pain. On the other hand, 5-HT1A antagonists have been revealed promising compounds in cognition disorders and, lately, in cancer. Thus, the discovery of 5HT1A ligands is nowadays an appealing research activity in medicinal chemistry. In this work, Comparative Molecular Fields Analysis (CoMFA) and Comparative Molecular Similarity Index Analysis (CoMSIA) were applied on an in-house library of 5-HT1A ligands bearing different chemical scaffolds in order to elucidate their affinity and selectivity for the target. Following this procedure, a number of structural modifications have been drawn for the development of much more effective 5-HT1AR ligands. [Formula: see text].

Benzimidazole-based derivatives as privileged scaffold developed for the treatment of the RSV infection: a computational study exploring the potency and cytotoxicity profiles.

Respiratory syncytial virus (RSV) has been identified as a main cause of hospitalisation in infants and children. To date, the current therapeutic arsenal is limited to ribavirin and palivizumab with variable efficacy. In this work, starting from a number of in-house series of previously described anti-RSV agents based on the benzimidazole scaffold, with the aim at gaining a better understanding of the related chemical features involved in potency and safety profiles, we applied a computational study including two focussed comparative molecular fields analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The results allowed us to derive useful suggestions for the design of derivatives and also to set up statistical models predicting the potency and selectivity index (SI = CC50/EC50) of any new analogue prior to synthesis. Accordingly, here, we discuss preliminary results obtained through the applied exhaustive QSAR analyses, leading to design and synthesise more effective anti-RSV agents.

New arylsparteine derivatives as positive inotropic drugs.

Positive inotropic agents are fundamental in the treatment of heart failure; however, their arrhythmogenic liability and the increased myocardial oxygen demand strongly limit their therapeutic utility. Pursuing our study on cardiovascular activities of lupin alkaloid derivatives, several 2-(4-substituted-phenyl)-2-dehydrosparteines and 2-(4-substituted-phenyl)sparteines were prepared and tested for inotropic and chronotropic activities on isolated guinea pig atria. Four compounds (6b, 6e, 7b, and 7f) exhibited significant inotropism that, at the higher concentrations, was followed by negative inotropism or toxicity. Compound 7e (2-(4-tolyl)sparteine) exhibited a steep dose-depending inotropic activity up to the highest concentration tested (300 µM) with an Emax of 116.5 ± 3.4% of basal force, proving less potent but much more active in comparison to the highest concentrations tested of digoxin and milrinone having Emax of 87.5 ± 3.1% and 52.2 ± 1.1%, respectively. Finally, docking studies suggested that the relevant sparteine derivatives could target the sigma-1 receptor, whose involvement in cardiac activity is well documented.

Molecular Chaperones in the Pathogenesis of Amyotrophic Lateral Sclerosis: The Role of HSPB1.

Genetic discoveries in amyotrophic lateral sclerosis (ALS) have a significant impact on deciphering molecular mechanisms of motor neuron degeneration but, despite recent advances, the etiology of most sporadic cases remains elusive. Several cellular mechanisms contribute to the motor neuron degeneration in ALS, including RNA metabolism, cellular interactions between neurons and nonneuronal cells, and seeding of misfolded protein with prion-like propagation. In this scenario, the importance of protein turnover and degradation in motor neuron homeostasis gained increased recognition. In this study, we evaluated the role of the candidate gene HSPB1, a molecular chaperone involved in several proteome-maintenance functions. In a cohort of 247 unrelated Italian ALS patients, we identified two variants (c.570G>C, p.Gln190His and c.610dupG, p.Ala204Glyfs* 6). Functional characterization of the p.Ala204Glyfs* 6 demonstrated that the mutant protein alters HSPB1 dynamic equilibrium, sequestering the wild-type protein in a stable dimer and resulting in a loss of chaperone-like activity. Our results underline the relevance of identifying rare but pathogenic variations in sporadic neurodegenerative diseases, suggesting a possible correlation between specific pathomechanisms linked to HSPB1 mutations and the associated neurological phenotype. Our study provides additional lines of evidence to support the involvement of HSPB1 in the pathogenesis of sporadic ALS.

Synthesis, molecular modeling and SAR study of novel pyrazolo[5,1-f][1,6]naphthyridines as CB2 receptor antagonists/inverse agonists.

Pyrazolo[5,1-f][1,6]naphthyridine-carboxamide derivatives were synthesized and evaluated for the affinity at CB1 and CB2 receptors. Based on the AgOTf and proline-cocatalyzed multicomponent methodology, the ethyl 5-(p-tolyl)pyrazolo[5,1-f][1,6]naphthyridine-2-carboxylate (12) and ethyl 5-(2,4-dichlorophenyl)pyrazolo[5,1-f][1,6]naphthyridine-2-carboxylate (13) intermediates were synthesized from the appropriate o-alkynylaldehydes, p-toluenesulfonyl hydrazide and ethyl pyruvate. Most of the novel compounds feature a p-tolyl (8a-i) or a 2,4-dichlorophenyl (8j) motif at the C5-position of the tricyclic pyrazolo[5,1-f][1,6]naphthyridine scaffold. Structural variation on the carboxamide moiety at the C2-position includes basic monocyclic, terpenoid and adamantine-based amines. Among these derivatives, compound 8h (N-adamant-1-yl-5-(p-tolyl)pyrazolo[5,1-f][1,6]naphthyridine-2-carboxamide) exhibited the highest CB2 receptor affinity (Ki=33nM) and a high degree of selectivity (KiCB1/KiCB2=173:1), whereas a similar trend in the near nM range was seen for the bornyl analogue (compound 8f, Ki=53nM) and the myrtanyl derivative 8j (Ki=67nM). Effects of 8h, 8f and 8j on forskolin-stimulated cAMP levels were determined, showing antagonist/inverse agonist properties for such compounds. Docking studies conducted for these derivatives and the reference antagonist/inverse agonist compound 4 (SR144528) disclosed the specific pattern of interactions probably related to the pyrazolo[5,1-f][1,6]naphthyridine scaffold as CB2 inverse agonists.

New insights into PDE4B inhibitor selectivity: CoMFA analyses and molecular docking studies.

PDE4 inhibitors have been largely studied because of their promising therapeutic effects concerning inflammation and neurodegenerative dysfunctions, such as depression, schizophrenia and Alzheimer's diseases. In this context, the PDE4B isoform proved to be particularly involved in the activation of inflammatory responses, while the PDE4D subfamily is more associated with neuropathologies. The clinical use of PDE4 inhibitors was restricted by the presence of prominent side effects probably due to their non-specific action across the different isoforms. Therefore, this work deals with the development of 3D-QSAR models, supported by molecular docking studies, to identify the key requirements underlying selective PDE4B or PDE4D inhibition. The results highlighted the ligand-based approach as a promising tool to guide the rational design of novel PDE4 inhibitors endowed with high affinity and selectivity profiles. The alignment of compound 1-85 and the model A statistical results are depicted.

Impaired Chaperone Activity of Human Heat Shock Protein Hsp27 Site-Specifically Modified with Argpyrimidine.

Non-enzymatic posttranslational modifications (nPTMs) affect at least ∼30 % of human proteins, but our understanding of their impact on protein structure and function is limited. Studies of nPTMs are difficult because many modifications are not included in common chemical libraries or protein expression systems and should be introduced site-specifically. Herein, we probed the effect of the nPTM argpyrimidine on the structure and function of human protein Hsp27, which acquires argpyrimidine at residue 188 in vivo. We developed a synthetic approach to an argpyrimidine building block, which we then incorporated at position 188 of Hsp27 through protein semisynthesis. This modification did not affect the protein secondary structure, but perturbed the oligomeric assembly and impaired chaperone activity. Our work demonstrates that protein function can be altered by a single nPTM and opens up a new area of investigation only accessible by methods that allow site-selective protein modification.

In silico evaluation of human small heat shock protein HSP27: homology modeling, mutation analyses and docking studies.

Small heat-shock proteins, possessing chaperone-like activity, represented crucial proteins actively involved in maintain protein homeostasis, which act to prevent improper polypeptide aggregation and deposition of misfolded proteins. In this context, a number of mutations concerning the HspB1 protein proved to be associated with the development of several neuropathologies. Unfortunately, molecular mechanisms underlying the onset of these diseases and in particular the changes induced by the mutations in HspB1 structure, remain poorly characterized. On the other hand, more recent studies demonstrated that HspB1 overexpression leads to an overactive chaperone activity, which in turn contributes to the anticancer agent resistance. On these basis, Hsp27 could represent a good innovative target for development of novel cancer therapy. Therefore, in this work a computational study, based on the homology model of the complete Hsp27 protein and of several pathological mutant forms, was developed. Finally, the derived model was employed to perform, for the first time, docking simulations on a recently identified Hsp27 inhibitor, disclosing a new useful panorama to be exploited for the further development of new compounds.