sp3 -Rich Heterocycle Synthesis on DNA: Application to DNA-Encoded Library Production.

DNA encoded library (DEL) synthesis represents a convenient means to produce, annotate and store large collections of compounds in a small volume. While DELs are well suited for drug discovery campaigns, the chemistry used in their production must be compatible with the DNA tag, which can limit compound class accessibility. As a result, most DELs are heavily populated with peptidomimetic and sp2 -rich molecules. Herein, we show that sp3 -rich mono- and bicyclic heterocycles can be made on DNA from ketochlorohydrin aldol products through a reductive amination and cyclization process. The resulting hydroxypyrrolidines possess structural features that are desirable for DELs and target a distinct region of pharmaceutically relevant chemical space.

Porin-independent accumulation in Pseudomonas enables antibiotic discovery.

Gram-negative antibiotic development has been hindered by a poor understanding of the types of compounds that can accumulate within these bacteria1,2. The presence of efflux pumps and substrate-specific outer-membrane porins in Pseudomonas aeruginosa renders this pathogen particularly challenging3. As a result, there are few antibiotic options for P. aeruginosa infections4 and its many porins have made the prospect of discovering general accumulation guidelines seem unlikely5. Here we assess the whole-cell accumulation of 345 diverse compounds in P. aeruginosa and Escherichia coli. Although certain positively charged compounds permeate both bacterial species, P. aeruginosa is more restrictive compared to E. coli. Computational analysis identified distinct physicochemical properties of small molecules that specifically correlate with P. aeruginosa accumulation, such as formal charge, positive polar surface area and hydrogen bond donor surface area. Mode of uptake studies revealed that most small molecules permeate P. aeruginosa using a porin-independent pathway, thus enabling discovery of general P. aeruginosa accumulation trends with important implications for future antibiotic development. Retrospective antibiotic examples confirmed these trends and these discoveries were then applied to expand the spectrum of activity of a gram-positive-only antibiotic, fusidic acid, into a version that demonstrates a dramatic improvement in antibacterial activity against P. aeruginosa. We anticipate that these discoveries will facilitate the design and development of high-permeating antipseudomonals.

LipMetE (Lipophilic Metabolism Efficiency) as a Simple Guide for Half-Life and Dosing Regimen Prediction of Oral Drugs.

The in vivo half-life is a key property of every drug molecule, as it determines dosing regimens, peak-to-trough ratios and often dose. However, half-life optimization can be challenging due to its multifactorial nature, with in vitro metabolic turnover, plasma protein binding and volume of distribution all impacting half-life. We here propose that the medicinal chemistry design parameter Lipophilic Metabolism Efficiency (LipMetE) can greatly simplify half-life optimization of neutral and basic compounds. Using mathematical transformations, examples from preclinical GABAA projects and clinical compounds with human pharmacokinetic data, we show that LipMetE is directly proportional to the logarithm of half-life. As the design parameter LipMetE can be swiftly calculated using the readily available parameters LogD, intrinsic clearance and fraction unbound in human liver microsomes or hepatocytes, this approach enables rational half-life optimization from the early stages of drug discovery projects.

The Playbooks of Medicinal Chemistry Design Moves.

Large databases of biologically relevant molecules, such as ChEMBL, SureChEMBL, or compound collections of pharmaceutical or agrochemical companies, are invaluable sources of medicinal chemistry information, albeit implicit. We developed a modified matched molecular pair approach to systematically and exhaustively extract the transformations in these databases and distill them into snippets of explicit design knowledge that are easily interpretable and directly applicable. The resulting "playbooks of medicinal chemistry design moves" capture the collective pharmaceutical and agrochemical research expertise across multiple chemists, companies, targets, and projects. They can be queried in an automated fashion for systematic prospective design and compound generation. The ChEMBL playbook and an application to exploit it are available at https://github.com/mahendra-awale/medchem_moves.

Tautomer Database: A Comprehensive Resource for Tautomerism Analyses.

We report a database of tautomeric structures that contains 2819 tautomeric tuples extracted from 171 publications. Each tautomeric entry has been annotated with experimental conditions reported in the respective publication, plus bibliographic details, structural identifiers (e.g., NCI/CADD identifiers FICTS, FICuS, uuuuu, and Standard InChI), and chemical information (e.g., SMILES, molecular weight). The majority of tautomeric tuples found were pairs; the remaining 10% were triples, quadruples, or quintuples, amounting to a total number of structures of 5977. The types of tautomerism were mainly prototropic tautomerism (79%), followed by ring-chain (13%) and valence tautomerism (8%). The experimental conditions reported in the publications included about 50 pure solvents and 9 solvent mixtures with 26 unique spectroscopic or nonspectroscopic methods. 1H and 13C NMR were the most frequently used methods. A total of 77 different tautomeric transform rules (SMIRKS) are covered by at least one example tuple in the database. This database is freely available as a spreadsheet at https://cactus.nci.nih.gov/download/tautomer/.

Epigenetic re-wiring of breast cancer by pharmacological targeting of C-terminal binding protein.

The C-terminal binding protein (CtBP) is an NADH-dependent dimeric family of nuclear proteins that scaffold interactions between transcriptional regulators and chromatin-modifying complexes. Its association with poor survival in several cancers implicates CtBP as a promising target for pharmacological intervention. We employed computer-assisted drug design to search for CtBP inhibitors, using quantitative structure-activity relationship (QSAR) modeling and docking. Functional screening of these drugs identified 4 compounds with low toxicity and high water solubility. Micro molar concentrations of these CtBP inhibitors produces significant de-repression of epigenetically silenced pro-epithelial genes, preferentially in the triple-negative breast cancer cell line MDA-MB-231. This epigenetic reprogramming occurs through eviction of CtBP from gene promoters; disrupted recruitment of chromatin-modifying protein complexes containing LSD1, and HDAC1; and re-wiring of activating histone marks at targeted genes. In functional assays, CtBP inhibition disrupts CtBP dimerization, decreases cell migration, abolishes cellular invasion, and improves DNA repair. Combinatorial use of CtBP inhibitors with the LSD1 inhibitor pargyline has synergistic influence. Finally, integrated correlation of gene expression in breast cancer patients with nuclear levels of CtBP1 and LSD1, reveals new potential therapeutic vulnerabilities. These findings implicate a broad role for this class of compounds in strategies for epigenetically targeted therapeutic intervention.

Experimental and Chemoinformatics Study of Tautomerism in a Database of Commercially Available Screening Samples.

We investigated how many cases of the same chemical sold as different products (at possibly different prices) occurred in a prototypical large aggregated database and simultaneously tested the tautomerism definitions in the chemoinformatics toolkit CACTVS. We applied the standard CACTVS tautomeric transforms plus a set of recently developed ring-chain transforms to the Aldrich Market Select (AMS) database of 6 million screening samples and building blocks. In 30 000 cases, two or more AMS products were found to be just different tautomeric forms of the same compound. We purchased and analyzed 166 such tautomer pairs and triplets by 1H and 13C NMR to determine whether the CACTVS transforms accurately predicted what is the same "stuff in the bottle". Essentially all prototropic transforms with examples in the AMS were confirmed. Some of the ring-chain transforms were found to be too "aggressive", i.e. to equate structures with one another that were different compounds.

Characterizing the Covalent Targets of a Small Molecule Inhibitor of the Lysine Acetyltransferase P300.

C646 inhibits the lysine acetyltransferases (KATs) p300 and CBP and represents the most potent and selective small molecule KAT inhibitor identified to date. To gain insights into the cellular activity of this epigenetic probe, we applied chemoproteomics to identify covalent targets of the C646 chemotype. Modeling and synthetic derivatization was used to develop a clickable analogue (C646-yne) that inhibits p300 similarly to the parent compound and enables enrichment of bound proteins. LC-MS/MS identified the major covalent targets of C646-yne as highly abundant cysteine-containing proteins, and follow-up studies found that C646 can inhibit tubulin polymerization in vitro. Finally, we provide evidence that thiol reactivity of C646 may limit its ability to antagonize acetylation in cells. These findings should enable a more precise interpretation of studies utilizing C646 as a chemical probe of KAT activity and suggest that an underappreciated liability of electrophile-containing inhibitors is a reduction in their cellular potency due to consumption by abundant protein and metabolite thiol sinks.

Novel HIV-1 Integrase Inhibitor Development by Virtual Screening Based on QSAR Models.

HIV-1 integrase (IN) plays an important role in the life cycle of HIV and is responsible for integration of the virus into the human genome. We present computational approaches used to design novel HIV-1 IN inhibitors. We created an IN inhibitor database by collecting experimental data from the literature. We developed quantitative structure-activity relationship (QSAR) models of HIV-1 IN strand transfer (ST) inhibitors using this database. The prediction accuracy of these models was estimated by external 5-fold cross-validation as well as with an additional validation set of 308 structurally distinct compounds from the publicly accessible BindingDB database. The validated models were used to screen a small combinatorial library of potential synthetic candidates to identify hits, with a subsequent docking approach applied to further filter out compounds to arrive at a small set of potential HIV-1 IN inhibitors. As result, 236 compounds with good druglikeness properties and with correct docking poses were identified as potential candidates for synthesis. One of the six compounds finally chosen for synthesis was experimentally confirmed to inhibit the ST reaction with an IC50(ST) of 37 µM. The IN inhibitor database is available for download from http://cactus.nci.nih.gov/download/iidb/.

Tautomerism of Warfarin: Combined Chemoinformatics, Quantum Chemical, and NMR Investigation.

Warfarin, an important anticoagulant drug, can exist in solution in 40 distinct tautomeric forms through both prototropic tautomerism and ring-chain tautomerism. We have investigated all warfarin tautomers with computational and NMR approaches. Relative energies calculated at the B3LYP/6-311G++(d,p) level of theory indicate that the 4-hydroxycoumarin cyclic hemiketal tautomer is the most stable tautomer in aqueous solution, followed by the 4-hydroxycoumarin open-chain tautomer. This is in agreement with our NMR experiments where the spectral assignments indicate that warfarin exists mainly as a mixture of cyclic hemiketal diastereomers, with an open-chain tautomer as a minor component. We present a diagram of the interconversion of warfarin created taking into account the calculated equilibrium constants (pK(T)) for all tautomeric reactions. These findings help with gaining further understanding of proton transfer and ring closure tautomerization processes. We also discuss the results in the context of chemoinformatics rules for handling tautomerism.

Metabolic Regulation of Histone Acetyltransferases by Endogenous Acyl-CoA Cofactors.

The finding that chromatin modifications are sensitive to changes in cellular cofactor levels potentially links altered tumor cell metabolism and gene expression. However, the specific enzymes and metabolites that connect these two processes remain obscure. Characterizing these metabolic-epigenetic axes is critical to understanding how metabolism supports signaling in cancer, and developing therapeutic strategies to disrupt this process. Here, we describe a chemical approach to define the metabolic regulation of lysine acetyltransferase (KAT) enzymes. Using a novel chemoproteomic probe, we identify a previously unreported interaction between palmitoyl coenzyme A (palmitoyl-CoA) and KAT enzymes. Further analysis reveals that palmitoyl-CoA is a potent inhibitor of KAT activity and that fatty acyl-CoA precursors reduce cellular histone acetylation levels. These studies implicate fatty acyl-CoAs as endogenous regulators of histone acetylation, and suggest novel strategies for the investigation and metabolic modulation of epigenetic signaling.

Peroxisome Proliferator-Activated Receptor γ (PPARγ) and Ligand Choreography: Newcomers Take the Stage.

Thiazolidinediones (TZDs), such as rosiglitazone and pioglitazone, are peroxisome proliferator-activated receptor γ (PPARγ) full agonists that have been widely used in the treatment of type 2 diabetes mellitus. Despite the demonstrated beneficial effect of reducing glucose levels in the plasma, TZDs also induce several adverse effects. Consequently, the search for new compounds with potent antidiabetic effects but fewer undesired effects is an active field of research. Interestingly, the novel proposed mechanisms for the antidiabetic activity of PPARγ agonists, consisting of PPARγ Ser273 phosphorylation inhibition, ligand and receptor mutual dynamics, and the presence of an alternate binding site, have recently changed the view regarding the optimal characteristics for the screening of novel PPARγ ligands. Furthermore, transcriptional genomics could bring essential information about the genome-wide effects of PPARγ ligands. Consequently, facing the new mechanistic scenario proposed for these compounds is essential for resolving the paradoxes among their agonistic function, antidiabetic activities, and side effects and should allow the rational development of better and safer PPARγ-mediated antidiabetic drugs.

Enumeration of ring-chain tautomers based on SMIRKS rules.

A compound exhibits (prototropic) tautomerism if it can be represented by two or more structures that are related by a formal intramolecular movement of a hydrogen atom from one heavy atom position to another. When the movement of the proton is accompanied by the opening or closing of a ring it is called ring-chain tautomerism. This type of tautomerism is well observed in carbohydrates, but it also occurs in other molecules such as warfarin. In this work, we present an approach that allows for the generation of all ring-chain tautomers of a given chemical structure. Based on Baldwin's Rules estimating the likelihood of ring closure reactions to occur, we have defined a set of transform rules covering the majority of ring-chain tautomerism cases. The rules automatically detect substructures in a given compound that can undergo a ring-chain tautomeric transformation. Each transformation is encoded in SMIRKS line notation. All work was implemented in the chemoinformatics toolkit CACTVS. We report on the application of our ring-chain tautomerism rules to a large database of commercially available screening samples in order to identify ring-chain tautomers.

Identification of PPARgamma partial agonists of natural origin (II): in silico prediction in natural extracts with known antidiabetic activity.

BACKGROUND: Natural extracts have played an important role in the prevention and treatment of diseases and are important sources for drug discovery. However, to be effectively used in these processes, natural extracts must be characterized through the identification of their active compounds and their modes of action. METHODOLOGY/PRINCIPAL FINDINGS: From an initial set of 29,779 natural products that are annotated with their natural source and using a previously developed virtual screening procedure (carefully validated experimentally), we have predicted as potential peroxisome proliferators-activated receptor gamma (PPARγ) partial agonists 12 molecules from 11 extracts known to have antidiabetic activity. Six of these molecules are similar to molecules with described antidiabetic activity but whose mechanism of action is unknown. Therefore, it is plausible that these 12 molecules could be the bioactive molecules responsible, at least in part, for the antidiabetic activity of the extracts containing them. In addition, we have also identified as potential PPARγ partial agonists 10 molecules from 16 plants with undescribed antidiabetic activity but that are related (i.e., they are from the same genus) to plants with known antidiabetic properties. None of the 22 molecules that we predict as PPARγ partial agonists show chemical similarity with a group of 211 known PPARγ partial agonists obtained from the literature. CONCLUSIONS/SIGNIFICANCE: Our results provide a new hypothesis about the active molecules of natural extracts with antidiabetic properties and their mode of action. We also suggest plants with undescribed antidiabetic activity that may contain PPARγ partial agonists. These plants represent a new source of potential antidiabetic extracts. Consequently, our work opens the door to the discovery of new antidiabetic extracts and molecules that can be of use, for instance, in the design of new antidiabetic drugs or functional foods focused towards the prevention/treatment of type 2 Diabetes Mellitus.

Identification of PPARgamma partial agonists of natural origin (I): development of a virtual screening procedure and in vitro validation.

BACKGROUND: Although there are successful examples of the discovery of new PPARγ agonists, it has recently been of great interest to identify new PPARγ partial agonists that do not present the adverse side effects caused by PPARγ full agonists. Consequently, the goal of this work was to design, apply and validate a virtual screening workflow to identify novel PPARγ partial agonists among natural products. METHODOLOGY/PRINCIPAL FINDINGS: We have developed a virtual screening procedure based on structure-based pharmacophore construction, protein-ligand docking and electrostatic/shape similarity to discover novel scaffolds of PPARγ partial agonists. From an initial set of 89,165 natural products and natural product derivatives, 135 compounds were identified as potential PPARγ partial agonists with good ADME properties. Ten compounds that represent ten new chemical scaffolds for PPARγ partial agonists were selected for in vitro biological testing, but two of them were not assayed due to solubility problems. Five out of the remaining eight compounds were confirmed as PPARγ partial agonists: they bind to PPARγ, do not or only moderately stimulate the transactivation activity of PPARγ, do not induce adipogenesis of preadipocyte cells and stimulate the insulin-induced glucose uptake of adipocytes. CONCLUSIONS/SIGNIFICANCE: We have demonstrated that our virtual screening protocol was successful in identifying novel scaffolds for PPARγ partial agonists.

Identification of novel human dipeptidyl peptidase-IV inhibitors of natural origin (Part II): in silico prediction in antidiabetic extracts.

BACKGROUND: Natural extracts play an important role in traditional medicines for the treatment of diabetes mellitus and are also an essential resource for new drug discovery. Dipeptidyl peptidase IV (DPP-IV) inhibitors are potential candidates for the treatment of type 2 diabetes mellitus, and the effectiveness of certain antidiabetic extracts of natural origin could be, at least partially, explained by the inhibition of DPP-IV. METHODOLOGY/PRINCIPAL FINDINGS: Using an initial set of 29,779 natural products that are annotated with their natural source and an experimentally validated virtual screening procedure previously developed in our lab (Guasch et al.; 2012) [1], we have predicted 12 potential DPP-IV inhibitors from 12 different plant extracts that are known to have antidiabetic activity. Seven of these molecules are identical or similar to molecules with described antidiabetic activity (although their role as DPP-IV inhibitors has not been suggested as an explanation for their bioactivity). Therefore, it is plausible that these 12 molecules could be responsible, at least in part, for the antidiabetic activity of these extracts through their inhibitory effect on DPP-IV. In addition, we also identified as potential DPP-IV inhibitors 6 molecules from 6 different plants with no described antidiabetic activity but that share the same genus as plants with known antidiabetic properties. Moreover, none of the 18 molecules that we predicted as DPP-IV inhibitors exhibits chemical similarity with a group of 2,342 known DPP-IV inhibitors. CONCLUSIONS/SIGNIFICANCE: Our study identified 18 potential DPP-IV inhibitors in 18 different plant extracts (12 of these plants have known antidiabetic properties, whereas, for the remaining 6, antidiabetic activity has been reported for other plant species from the same genus). Moreover, none of the 18 molecules exhibits chemical similarity with a large group of known DPP-IV inhibitors.

Identification of novel human dipeptidyl peptidase-IV inhibitors of natural origin (part I): virtual screening and activity assays.

BACKGROUND: There has been great interest in determining whether natural products show biological activity toward protein targets of pharmacological relevance. One target of particular interest is DPP-IV whose most important substrates are incretins that, among other beneficial effects, stimulates insulin biosynthesis and secretion. Incretins have very short half-lives because of their rapid degradation by DPP-IV and, therefore, inhibiting this enzyme improves glucose homeostasis. As a result, DPP-IV inhibitors are of considerable interest to the pharmaceutical industry. The main goals of this study were (a) to develop a virtual screening process to identify potential DPP-IV inhibitors of natural origin; (b) to evaluate the reliability of our virtual-screening protocol by experimentally testing the in vitro activity of selected natural-product hits; and (c) to use the most active hit for predicting derivatives with higher binding affinities for the DPP-IV binding site. METHODOLOGY/PRINCIPAL FINDINGS: We predicted that 446 out of the 89,165 molecules present in the natural products subset of the ZINC database would inhibit DPP-IV with good ADMET properties. Notably, when these 446 molecules were merged with 2,342 known DPP-IV inhibitors and the resulting set was classified into 50 clusters according to chemical similarity, there were 12 clusters that contained only natural products for which no DPP-IV inhibitory activity has been previously reported. Nine molecules from 7 of these 12 clusters were then selected for in vitro activity testing and 7 out of the 9 molecules were shown to inhibit DPP-IV (where the remaining two molecules could not be solubilized, preventing the evaluation of their DPP-IV inhibitory activity). Then, the hit with the highest activity was used as a lead compound in the prediction of more potent derivatives. CONCLUSIONS/SIGNIFICANCE: We have demonstrated that our virtual-screening protocol was successful in identifying novel lead compounds for developing more potent DPP-IV inhibitors.

Development of docking-based 3D-QSAR models for PPARgamma full agonists.

Peroxisome proliferator-activated receptor gamma (PPARγ) has become an attractive molecular target for drugs that aim to treat diabetes mellitus type II, and its therapeutic potency against skin cancer and other skin diseases is also currently being explored. To study the relationship between the structure of several PPARγ full agonists and the trans-activation activity of PPARγ, we have performed a three-dimensional quantitative structure-activity relationship (3D-QSAR) study of tyrosine-based derivatives, based on the 3D alignment of conformations obtained by docking. Highly predictive 3D-QSAR models, with Pearson-R values of 0.86 and 0.90, were obtained for the transactivation activity and binding affinity of PPARγ, respectively. These models are in good agreement with the structural characteristics of the binding pocket of PPARγ and provide some structural insights for the improvement of PPARγ full agonist bioactivities.

DecoyFinder: an easy-to-use python GUI application for building target-specific decoy sets.

UNLABELLED: Decoys are molecules that are presumed to be inactive against a target (i.e. will not likely bind to the target) and are used to validate the performance of molecular docking or a virtual screening workflow. The Directory of Useful Decoys database (http://dud.docking.org/) provides a free directory of decoys for use in virtual screening, though it only contains a limited set of decoys for 40 targets.To overcome this limitation, we have developed an application called DecoyFinder that selects, for a given collection of active ligands of a target, a set of decoys from a database of compounds. Decoys are selected if they are similar to active ligands according to five physical descriptors (molecular weight, number of rotational bonds, total hydrogen bond donors, total hydrogen bond acceptors and the octanol-water partition coefficient) without being chemically similar to any of the active ligands used as an input (according to the Tanimoto coefficient between MACCS fingerprints). To the best of our knowledge, DecoyFinder is the first application designed to build target-specific decoy sets. AVAILABILITY: A complete description of the software is included on the application home page. A validation of DecoyFinder on 10 DUD targets is provided as Supplementary Table S1. DecoyFinder is freely available at http://URVnutrigenomica-CTNS.github.com/DecoyFinder.