Similarity Mapplet: Interactive Visualization of the Directory of Useful Decoys and ChEMBL in High Dimensional Chemical Spaces

An Internet portal accessible at www.gdb.unibe.ch has been set up to automatically generate color-coded similarity maps of the ChEMBL database in relation to up to two sets of active compounds taken from the enhanced Directory of Useful Decoys (eDUD), a random set of molecules, or up to two sets of user-defined reference molecules. These maps visualize the relationships between the selected compounds and ChEMBL in six different high dimensional chemical spaces, namely MQN (42-D molecular quantum numbers), SMIfp (34-D SMILES fingerprint), APfp (20-D shape fingerprint), Xfp (55-D pharmacophore fingerprint), Sfp (1024-bit substructure fingerprint), and ECfp4 (1024-bit extended connectivity fingerprint). The maps are supplied in form of Java based desktop applications called “similarity mapplets” allowing interactive content browsing and linked to a “Multifingerprint Browser for ChEMBL” (also accessible directly at www.gdb.unibe.ch) to perform nearest neighbor searches. One can obtain six similarity mapplets of ChEMBL relative to random reference compounds, 606 similarity mapplets relative to single eDUD active sets, 30 300 similarity mapplets relative to pairs of eDUD active sets, and any number of similarity mapplets relative to user-defined reference sets to help visualize the structural diversity of compound series in drug optimization projects and their relationship to other known bioactive compounds.

Discovery and characterization of a novel non-competitive inhibitor of the divalent metal transporter DMT1/SLC11A2

Divalent metal transporter-1 (SLC11A2/DMT1) uses the H+ electrochemical gradient as the driving force to transport divalent metal ions such as Fe2+, Mn2+ and others metals into mammalian cells. DMT1 is ubiquitously expressed, most notably in proximal duodenum, immature erythroid cells, brain and kidney. This transporter mediates H+-coupled transport of ferrous iron across the apical membrane of enterocytes. In addition, in cells such as to erythroid precursors, following transferrin receptor (TfR) mediated endocytosis; it mediates H+-coupled exit of ferrous iron from endocytic vesicles into the cytosol. Dysfunction of human DMT1 is associated with several pathologies such as iron deficiency anemia hemochromatosis, Parkinson's disease and Alzheimer's disease, as well as colorectal cancer and esophageal adenocarcinoma, making DMT1 an attractive target for drug discovery. In the present study, we performed a ligand-based virtual screening of the Princeton database (700,000 commercially available compounds) to search for pharmacophore shape analogs of recently reported DMT1 inhibitors. We discovered a new compound, named pyrimidinone 8, which mediates a reversible linear non-competitive inhibition of human DMT1 (hDMT1) transport activity with a Ki of ∼20 μM. This compound does not affect hDMT1 cell surface expression and shows no dependence on extracellular pH. To our knowledge, this is the first experimental evidence that hDMT1 can be allosterically modulated by pharmacological agents. Pyrimidinone 8 represents a novel versatile tool compound and it may serve as a lead structure for the development of therapeutic compounds for pre-clinical assessment.

Stereoselective virtual screening of the ZINC database using atom pair 3D-fingerprints

Background Tools to explore large compound databases in search for analogs of query molecules provide a strategically important support in drug discovery to help identify available analogs of any given reference or hit compound by ligand based virtual screening (LBVS). We recently showed that large databases can be formatted for very fast searching with various 2D-fingerprints using the city-block distance as similarity measure, in particular a 2D-atom pair fingerprint (APfp) and the related category extended atom pair fingerprint (Xfp) which efficiently encode molecular shape and pharmacophores, but do not perceive stereochemistry. Here we investigated related 3D-atom pair fingerprints to enable rapid stereoselective searches in the ZINC database (23.2 million 3D structures). Results Molecular fingerprints counting atom pairs at increasing through-space distance intervals were designed using either all atoms (16-bit 3DAPfp) or different atom categories (80-bit 3DXfp). These 3D-fingerprints retrieved molecular shape and pharmacophore analogs (defined by OpenEye ROCS scoring functions) of 110,000 compounds from the Cambridge Structural Database with equal or better accuracy than the 2D-fingerprints APfp and Xfp, and showed comparable performance in recovering actives from decoys in the DUD database. LBVS by 3DXfp or 3DAPfp similarity was stereoselective and gave very different analogs when starting from different diastereomers of the same chiral drug. Results were also different from LBVS with the parent 2D-fingerprints Xfp or APfp. 3D- and 2D-fingerprints also gave very different results in LBVS of folded molecules where through-space distances between atom pairs are much shorter than topological distances. Conclusions 3DAPfp and 3DXfp are suitable for stereoselective searches for shape and pharmacophore analogs of query molecules in large databases. Web-browsers for searching ZINC by 3DAPfp and 3DXfp similarity are accessible at www.gdb.unibe.ch webcite and should provide useful assistance to drug discovery projects.

Molecular determinants of tissue selectivity in estrogen receptor modulators

Interaction of the estrogen receptor/ligand complex with a DNA estrogen response element is known to regulate gene transcription. In turn, specific conformations of the receptor-ligand complex have been postulated to influence unique subsets of estrogen-responsive genes resulting in differential modulation and, ultimately, tissue-selective outcomes. The estrogen receptor ligands raloxifene and tamoxifen have demonstrated such tissue-specific estrogen agonist/antagonist effects. Both agents antagonize the effects of estrogen on mammary tissue while mimicking the actions of estrogen on bone. However, tamoxifen induces significant stimulation of uterine tissue whereas raloxifene does not. We postulate that structural differences between raloxifene and tamoxifen may influence the conformations of their respective receptor/ligand complexes, thereby affecting which estrogen-responsive genes are modulated in various tissues. These structural differences are 4-fold: (A) the presence of phenolic hydroxyls, (B) different substituents on the basic amine, (C) incorporation of the stilbene moiety into a cyclic benzothiophene framework, and (D) the imposition of a carbonyl “hinge” between the basic amine-containing side chain and the olefin. A series of raloxifene analogs that separately exemplify each of these differences have been prepared and evaluated in a series of in vitro and in vivo assays. This strategy has resulted in the development of a pharmacophore model that attributes the differences in effects on the uterus between raloxifene and tamoxifen to a low-energy conformational preference imparting an orthogonal orientation of the basic side chain with respect to the stilbene plane. This three-dimensional array is dictated by a single carbon atom in the hinge region of raloxifene. These data indicate that differences in tissue selective actions among benzothiophene and triarylethylene estrogen receptor modulators can be ascribed to discrete ligand conformations.

Role of glutamic acid 216 in cytochrome P450 2D6 substrate binding and catalysis

Human cytochrome P450 (P450) 2D6 is an important enzyme involved in the metabolism of drugs, many of which are amines or contain other basic nitrogen atoms. Asp301 has generally been considered to be involved in electrostatic docking with the basic substrates, on the basis of previous modeling studies and site-directed mutagenesis. Substitution of Glu216 with a residue other than Asp strongly attenuated the binding of quinidine, bufuralol, and several other P450 2D6 ligands. Catalytic activity with the substrates bufuralol and 4-methoxyphenethylamine was strongly inhibited by neutral or basic mutations at Glu216 (>95%), to the same extent as the substitution of Asn at Asp301. Unlike the Asp301 mutants, the Gln216 mutant (E216Q) retained 40% enzyme efficiency with the substrate spirosulfonamide, devoid of basic nitrogen, suggesting that the substitutions at Glu216 affect binding of amine substrates more than other catalytic steps. Attempts to induce catalytic specificity toward new substrates by substitutions at Asp301 and Glu216 were unsuccessful. Collectively, the results provide evidence for electrostatic interaction of amine substrates with Glu216, and we propose that both of these acidic residues plus at least another residue(s) is (are) involved in binding the repertoire of P450 2D6 ligands.

Nematocidal thiocyanatins from a southern Australian marine sponge Oceanapia sp.

Investigations of a southern Australian marine sponge, Oceanapia sp., have yielded two new methyl branched bisthiocyanates, thiocyanatins D-1 (3a) and D-2 (3b), along with two new thiocarbamate thiocyanates, thiocyanatins E-l (4a) and E-2 (4b). The new thiocyanatins belong to a rare class of bioactive marine metabolite previously only represented by thiocyanatins A-C (1, 2a/b). Structures were assigned on the basis of detailed spectroscopic analysis, with comparisons to the known bisthiocyanate thiocyanatin A (1) and synthetic model compounds (5-7). The thiocyanatins exhibit potent nematocidal activity, and preliminary structure-activity relationship investigations have confirmed key characteristics of the thiocyanatin pharmacophore.

Potent cyclic antagonists of the complement C5a receptor on human polymorphonulcear leukocytes. Relationships between structures and activity

Human C5a is a plasma protein with potent chemoattractant and pro-inflammatory properties, and its overexpression correlates with severity of inflammatory diseases. C5a binds to its G protein-coupled receptor (C5aR) on polymorphonuclear leukocytes (PMNLs) through a high-affinity helical bundle and a low-affinity C terminus, the latter being solely responsible for receptor activation. Potent and selective C5a antagonists are predicted to be effective anti-inflammatory drugs, but no pharmacophore for small molecule antagonists has yet been developed, and it would significantly aid drug design. We have hypothesized that a turn conformation is important for activity of the C terminus of C5a and herein report small cyclic peptides that are stable turn mimics with potent antagonism at C5aR on human PMNLs. A comparison of solution structures for the C terminus of C5a, small acyclic peptide ligands, and cyclic antagonists supports the importance of a turn for receptor binding. Competition between a cyclic antagonist and either C5a or an acyclic agonist for C5aR on PMNLs supports a common or overlapping binding site on the C5aR. Structure-activity relationships for 60 cyclic analogs were evaluated by competitive radioligand binding with C5a (affinity) and myeloperoxidase release (antagonist potency) from human PMNLs, with 20 compounds having high antagonist potencies (IC50, 20 nM(-1) muM). Computer modeling comparisons reveal that potent antagonists share a common cyclic backbone shape, with affinity-determining side chains of defined volume projecting from the cyclic scaffold. These results define a new pharmacophore for C5a antagonist development and advance our understanding of ligand recognition and receptor activation of this G protein-coupled receptor.

Carbohydrate-based scaffolds in drug discovery

Carbohydrates have been proven as valuable scaffolds to display pharmocophores and the resulting molecules have demonstrated useful biological activity towards various targets including the somatostatin receptors (SSTR), integrins, HIV-1 protease, matrix metalloproteinases (MMP), multidrug resistance-associated protein (MRP), and as RNA binders. Carbohydrate-based compounds have also shown antibacterial and herbicidal activity.

Isokibdelones: Novel heterocyclic polyketides from a Kibdelosporangium sp.

The isokibdelones are an unprecedented family of polyketides produced by an Australian isolate of a rare actinomycete, Kibdelosporangium sp. The structures of the isokibdelones were assigned by spectroscopic analysis and chemical interconversion. A proposed biosynthesis requires a novel molecular twist that generates an unprecedented heterocyclic system and differentiates the isokibdelones from their kibdelone co-metabolites. SAR analysis on the isokibdelones further defines the anticancer pharmacophore of these novel polyketides.

Multi-scale approaches to Food Sciences: computational QM and MM explorations at the microscopic level

L'obiettivo principale della politica di sicurezza alimentare è quello di garantire la salute dei consumatori attraverso regole e protocolli di sicurezza specifici. Al fine di rispondere ai requisiti di sicurezza alimentare e standardizzazione della qualità, nel 2002 il Parlamento Europeo e il Consiglio dell'UE (Regolamento (CE) 178/2002 (CE, 2002)), hanno cercato di uniformare concetti, principi e procedure in modo da fornire una base comune in materia di disciplina degli alimenti e mangimi provenienti da Stati membri a livello comunitario. La formalizzazione di regole e protocolli di standardizzazione dovrebbe però passare attraverso una più dettagliata e accurata comprensione ed armonizzazione delle proprietà globali (macroscopiche), pseudo-locali (mesoscopiche), ed eventualmente, locali (microscopiche) dei prodotti alimentari. L'obiettivo principale di questa tesi di dottorato è di illustrare come le tecniche computazionali possano rappresentare un valido supporto per l'analisi e ciò tramite (i) l’applicazione di protocolli e (ii) miglioramento delle tecniche ampiamente applicate. Una dimostrazione diretta delle potenzialità già offerte dagli approcci computazionali viene offerta nel primo lavoro in cui un virtual screening basato su docking è stato applicato al fine di valutare la preliminare xeno-androgenicità di alcuni contaminanti alimentari. Il secondo e terzo lavoro riguardano lo sviluppo e la convalida di nuovi descrittori chimico-fisici in un contesto 3D-QSAR. Denominata HyPhar (Hydrophobic Pharmacophore), la nuova metodologia così messa a punto è stata usata per esplorare il tema della selettività tra bersagli molecolari strutturalmente correlati e ha così dimostrato di possedere i necessari requisiti di applicabilità e adattabilità in un contesto alimentare. Nel complesso, i risultati ci permettono di essere fiduciosi nel potenziale impatto che le tecniche in silico potranno avere nella identificazione e chiarificazione di eventi molecolari implicati negli aspetti tossicologici e nutrizionali degli alimenti.

Molecular modelling assisted design and synthesis of cyclothialidine derivatives as DNA gyrase inhibitors

Since cyclothialidine was discovered as the most active DNA gyrase inhibitor in 1994, enormous efforts have been devoted to make it into a commercial medicine by a number of pharmaceutical companies and research groups worldwide. However, no serious breakthrough has been made up to now. An essential problem involved with cyclothialidine is that though it demonstrated the potent inhibition of DNA gyrase, it showed little activity against bacteria. This probably is attributable to its inability to penetrate bacterial cell walls and membranes. We applied the TSAR programme to generate a QSAR equation to the gram-negative organisms. In that equation, LogP is profoundly indicated as the key factor influencing the cyclothialidine activity against bacteria. However, the synthesized new analogues have failed to prove that. In the structure based drug design stage, we designed a group of open chain cyclothialidine derivatives by applying the SPROUT programme and completed the syntheses. Improved activity is found in a few analogues and a 3D pharmacophore of the DNA gyrase B is proposed to lead to synthesis of the new derivatives for development of potent antibiotics.

Design, synthesis and evaluation of cyclothialidine analogues as DNA gyrase inhibitors

Cyclothialidine, a natural product isolated from Streptomyces .filipinensis NR0484, has been proven to be a potent and selective inhibitor of the bacterial enzyme DNA gyrase. Gyrase inhibition results in cell death, the enzyme being the target of several currently used antibiotics. Cyclothialidine showed poor activity against whole bacterial cells, highlighting scope for improvement regarding cell membrane pemeability in order for the full potential of this new class of antibiotics to be realised, Structurally, cyclothialidine contains a 12-membered lactone ring which is partly integrated into a pentapeptide chain, with a substituted aromatic moiety bordering the lactone, Retrosynthetically it can be traced back to cis-3-hydroxyproline, 3,5-dihydroxy-2,6-dimethylbenzoic acid and four commercially available amino acids; two serine, one cysteine and one alanine. In this work, a model of cyclothialidine was synthesised in order to establish the methodology for more complex compounds. Analogues with hydroxy, dihydroxy and dihydroxymethyl substituted aromatic moieties were then prepared to ensure successful protection methods could be performed and the pharmacophore synthesised. The key aromatic moiety, 2,6-dimethyl-3,5-dihydroxybenzoic acid was produced via two successive Mannich reaction/reduction steps. Acid protection using 4-nitrobenzyl bromide and TBDMS hydroxyl protection followed by bromination of one methyl afforded the desired intermediate. Reaction with a serine/cysteine dipeptide, followed by deprotection and cyclisation under Mitsunobu conditions lead to the 12-membered lactone. An amine substituted aromatic analogue and also replacement of the cysteine sulphur by oxygen were attempted but without success. In an effort to improve cell permeability, a conjugate was synthesised between the pharmacophore and a cholesterol moiety. It was hoped the steroid fragment would serve to increase potency by escorting the molecule through the lipid environment of the cell membrane. The pharmacophore and conjugate were tested against a variety of bacterial strains but the conjugate failed to improve activity.

Investigation of the antiproliferative properties of tumour promoting phorbol esters and related compounds

Tumour promoting phorbol esters such as 12-0-tetradecanoylphorbol-13-acetate (TPA) exert a multitude of biological effects on many cellular systems, many of which are believed to be mediated via the activation of the enzyme protein kinase C (PKC). TPA and other biologically active phorbol esters inhibited the proliferation of the A549 human lung carcinoma cell line. However, after 5-6 days culture in the continued presence of the phorbol ester cells began to proliferate at a rate similar to that of untreated cells. Resistance to TPA was lost following subculturing, although subculture in the presence of 10 nM TPA for more than 9 weeks resulted in a more resistant phenotype. The selection of a TPA-resistant subpopulation was not responsible for the observed resistance. The antiproliferative properties of other PKC activators were investigated. Mezerein induced the same antiproliferative effects as TPA but synthetic diacylglycerols (DAGs), the presumed physiological ligands of PKC, exerted only a non-specific cytotoxic influence on growth. Bryostatins 1 and 2 were able to induce transient growth arrest of A549 cells in a manner similar to phorbol esters at nanomolar concentrations, but at higher concentrations blocked both their own antiproliferative action and also that of phorbol esters and mezerein. Fourteen compounds synthesized to mimic features of the phorbol ester pharmacophore and/or DAGs did not mimic the antiproliferative properties of TPA in A549 cells and exerted only a DAG-like non-specific cytotoxicity at high concentrations. The subcellular distribution and activity of PKC was determined following partial purification by non-denaturing polyacrylamide gel electrophoresis. Treatment with TPA, mezerein or bryostatins resulted in a concentration-dependent shift of PKC activity from the cytosol to cellular membranes within 30 min. Significant translocation was not observed on treatment with DAGs. Chronic exposure of cells to TPA caused a time- and concentration dependent down-regulation of functional PKC activity. A complete loss of PKC activity was also observed on treatment with growth-inhibitory concentrations of bryostatins. No PKC activity was detected in cells resistant to the growth-inhibitory influence of TPA. Measurement of intracellular Ca2+ concentrations using A549 cells cultured on Cytodex 1 microcarrier beads revealed that TPA, mezerein and the bryostatins induced a similar rapid rise in intracellular Ca2+ levels.

Free fatty acid receptor 1 (GPR40) agonists containing spirocyclic periphery inspired by LY2881835

The free fatty acid receptor 1 (FFA1), a G protein-coupled receptor (GPCR) naturally activated by long-chain fatty acids is a novel target for the treatment of metabolic diseases. The basic amine spirocyclic periphery of Eli Lilly's drug candidate LY2881835 for treatment of type 2 diabetes mellitus (which reached phase I clinical trials) inspired a series of novel FFA1 agonists. These were designed to incorporate the 3-[4-(benzyloxy)phenyl]propanoic acid pharmacophore core decorated with a range of spirocyclic motifs. The latter were prepared via the Prins cyclization and subsequent modification of the 4-hydroxytetrahydropyran moiety in the Prins product. Here, we synthesize 19 compounds and test for FFA1 activity. Within this pilot set, a nanomolar potency (EC50=55nM) was reached. Four lead compounds (EC50 range 55-410nM) were characterized for aqueous solubility, metabolic stability, plasma protein binding and Caco-2 permeability. While some instability in the presence of mouse liver microsomes was noted, mouse pharmacokinetic profile of the compound having the best overall ADME properties was evaluated to reveal acceptable bioavailability (F=10.3%) and plasma levels achieved on oral administration.

Design and synthesis of novel quinolones directed to the Plasmodium falciparum bc1 protein complex

Tese de Doutoramento, Química, Especialização em Química Orgânica, Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2016