Abstract: BACKGROUND: LuxI is a component of the quorum sensing signaling pathway in Vibrio fischeriresponsible for the inducer synthesis that is essential for bioluminescence. METHODS: Homology modeling of LuxI was carried out using Phyre2 and refined with the GalaxyWEB server. Five models were generated and evaluated by ERRAT, ANOLEA, QMEAN6, and Procheck. RESULTS: Five refined models were generated by the GalaxyWEB server, with Model 4 having the greatest quality based on the QMEAN6 score of 0.732. ERRAT analysis revealed an overall quality of 98.9%, while the overall quality of the initial model was 54%. The mean force potential energy, as analyzed by ANOLEA, were better compared to the initial model. Sterochemical quality estimation by Procheck showed that the refined Model 4 had a reliable structure, and was therefore submitted to the protein model database. Drug Discovery Workbench V.2 was used to screen 2700 experimental compounds from the DrugBank database to identify inhibitors that can bind to the active site between amino acids 24 and 110. Ten compounds with high negative scores were selected as the best in binding. CONCLUSION: The model produced, and the predicted acteyltransferase binding site, could be useful in modeling homologous sequences from other microorganisms and the design of new antimicrobials.

Abstract: Chagas disease is a parasitic infection caused by the protozoa Trypanosoma cruzi that affects about 6 million people in Latin America. Despite its sanitary importance, there are currently only two drugs available for treatment: benznidazole and nifurtimox, both exhibiting serious adverse effects and limited efficacy in the chronic stage of the disease. Polyamines are ubiquitous to all living organisms where they participate in multiple basic functions such as biosynthesis of nucleic acids and proteins, proliferation and cell differentiation. T. cruzi is auxotroph for polyamines, which are taken up from the extracellular medium by efficient transporters and, to a large extent, incorporated into trypanothione (bis-glutathionylspermidine), the major redox cosubstrate of trypanosomatids. From a 268-compound database containing polyamine analogs with and without inhibitory effect on T. cruzi we have inferred classificatory models that were later applied in a virtual screening campaign to identify anti-trypanosomal compounds among drugs already used for other therapeutic indications (i.e. computer-guided drug repositioning) compiled in the DrugBank and Sweetlead databases. Five of the candidates identified with this strategy were evaluated in cellular models from different pathogenic trypanosomatids (T. cruzi wt, T. cruzi PAT12, T. brucei and Leishmania infantum), and in vitro models of aminoacid/polyamine transport assays and trypanothione synthetase inhibition assay. Triclabendazole, sertaconazole and paroxetine displayed inhibitory effects on the proliferation of T. cruzi (epimastigotes) and the uptake of putrescine by the parasite. They also interfered with the uptake of others aminoacids and the proliferation of infective T. brucei and L. infantum (promastigotes). Trypanothione synthetase was ruled out as molecular target for the anti-parasitic activity of these compounds.

Abstract: High-throughput, 'omic' methods provide sensitive measures of biological responses to perturbations. However, inherent biases in high-throughput assays make it difficult to interpret experiments in which more than one type of data is collected. In this work, we introduce Omics Integrator, a software package that takes a variety of 'omic' data as input and identifies putative underlying molecular pathways. The approach applies advanced network optimization algorithms to a network of thousands of molecular interactions to find high-confidence, interpretable subnetworks that best explain the data. These subnetworks connect changes observed in gene expression, protein abundance or other global assays to proteins that may not have been measured in the screens due to inherent bias or noise in measurement. This approach reveals unannotated molecular pathways that would not be detectable by searching pathway databases. Omics Integrator also provides an elegant framework to incorporate not only positive data, but also negative evidence. Incorporating negative evidence allows Omics Integrator to avoid unexpressed genes and avoid being biased toward highly-studied hub proteins, except when they are strongly implicated by the data. The software is comprised of two individual tools, Garnet and Forest, that can be run together or independently to allow a user to perform advanced integration of multiple types of high-throughput data as well as create condition-specific subnetworks of protein interactions that best connect the observed changes in various datasets. It is available at http://fraenkel.mit.edu/omicsintegrator and on GitHub at https://github.com/fraenkel-lab/OmicsIntegrator.

Abstract: BACKGROUND: Therapeutic treatments for schizophrenia do not alleviate symptoms for all patients and efficacy is limited by common, often severe, side-effects. Genetic studies of disease can identify novel drug targets, and drugs for which the mechanism has direct genetic support have increased likelihood of clinical success. Large-scale genetic studies of schizophrenia have increased the number of genes and gene sets associated with risk. We aimed to examine the overlap between schizophrenia risk loci and gene targets of a comprehensive set of medications to potentially inform and improve treatment of schizophrenia. METHODS: We defined schizophrenia risk loci as genomic regions reaching genome-wide significance in the latest Psychiatric Genomics Consortium schizophrenia genome-wide association study (GWAS) of 36 989 cases and 113 075 controls and loss of function variants observed only once among 5079 individuals in an exome-sequencing study of 2536 schizophrenia cases and 2543 controls (Swedish Schizophrenia Study). Using two large and orthogonally created databases, we collated drug targets into 167 gene sets targeted by pharmacologically similar drugs and examined enrichment of schizophrenia risk loci in these sets. We further linked the exome-sequenced data with a national drug registry (the Swedish Prescribed Drug Register) to assess the contribution of rare variants to treatment response, using clozapine prescription as a proxy for treatment resistance. FINDINGS: We combined results from testing rare and common variation and, after correction for multiple testing, two gene sets were associated with schizophrenia risk: agents against amoebiasis and other protozoal diseases (106 genes, p=0.00046, pcorrected =0.024) and antipsychotics (347 genes, p=0.00078, pcorrected=0.046). Further analysis pointed to antipsychotics as having independent enrichment after removing genes that overlapped these two target sets. We noted significant enrichment both in known targets of antipsychotics (70 genes, p=0.0078) and novel predicted targets (277 genes, p=0.019). Patients with treatment-resistant schizophrenia had an excess of rare disruptive variants in gene targets of antipsychotics (347 genes, p=0.0067) and in genes with evidence for a role in antipsychotic efficacy (91 genes, p=0.0029). INTERPRETATION: Our results support genetic overlap between schizophrenia pathogenesis and antipsychotic mechanism of action. This finding is consistent with treatment efficacy being polygenic and suggests that single-target therapeutics might be insufficient. We provide evidence of a role for rare functional variants in antipsychotic treatment response, pointing to a subset of patients where their genetic information could inform treatment. Finally, we present a novel framework for identifying treatments from genetic data and improving our understanding of therapeutic mechanism. FUNDING: US National Institutes of Health.

Abstract: A novel functional-genomics based concept of pharmacology that uses artificial intelligence techniques for mining and knowledge discovery in "big data" providing comprehensive information about the drugs' targets and their functional genomics is proposed. In "process pharmacology", drugs are associated with biological processes. This puts the disease, regarded as alterations in the activity in one or several cellular processes, in the focus of drug therapy. In this setting, the molecular drug targets are merely intermediates. The identification of drugs for therapeutic or repurposing is based on similarities in the high-dimensional space of the biological processes that a drug influences. Applying this principle to data associated with lymphoblastic leukemia identified a short list of candidate drugs, including one that was recently proposed as novel rescue medication for lymphocytic leukemia. The pharmacological data science approach provides successful selections of drug candidates within development and repurposing tasks.

Abstract: Drug discovery has produced many successful therapeutic agents; however, most of these drugs were developed without a deep understanding of the system-wide mechanisms of action responsible for their indications. Gene-disease associations produced by molecular and genetic studies of complex diseases provide great opportunities for a system-level understanding of drug activity. In this study, we focused on acute myocardial infarction (MI) and conducted an integrative network analysis to illuminate drug actions. We integrated MI drugs, MI drug interactors, drug targets, and MI disease genes into the human interactome and showed that MI drug targets are significantly proximate to MI disease proteins. We then constructed a bipartite network of MI-related drug targets and MI disease proteins and derived 12 drug-target-disease (DTD) modules. We assessed the biological relevance of these modules and demonstrated the benefits of incorporating disease genes. The results indicate that DTD modules provide insights into the mechanisms of action of MI drugs and the cardiovascular (side) effects of non-MI drugs.

Abstract: BACKGROUND: Drug repositioning refers to the process of developing new indications for existing drugs. As a phenotypic indicator of drug response in humans, clinical side effects may provide straightforward signals and unique opportunities for drug repositioning. OBJECTIVE: We aimed to identify drugs frequently associated with hypotension adverse reactions (ie, the opposite condition of hypertension), which could be potential candidates as antihypertensive agents. METHODS: We systematically searched the electronic records of the US Food and Drug Administration (FDA) Adverse Event Reporting System (FAERS) through the openFDA platform to assess the association between hypotension incidence and antihypertensive therapeutic effect regarding a list of 683 drugs. RESULTS: Statistical analysis of FAERS data demonstrated that those drugs frequently co-occurring with hypotension events were more likely to have antihypertensive activity. Ranked by the statistical significance of frequent hypotension reporting, the well-known antihypertensive drugs were effectively distinguished from others (with an area under the receiver operating characteristic curve > 0.80 and a normalized discounted cumulative gain of 0.77). In addition, we found a series of antihypertensive agents (particularly drugs originally developed for treating nervous system diseases) among the drugs with top significant reporting, suggesting the good potential of Web-based and data-driven drug repositioning. CONCLUSIONS: We found several candidate agents among the hypotension-related drugs on our list that may be redirected for lowering blood pressure. More important, we showed that a pharmacovigilance system could alternatively be used to identify antihypertensive agents and sustainably create opportunities for drug repositioning.

Abstract: The advances of genomics, sequencing, and high throughput technologies have led to the creation of large volumes of diverse datasets for drug discovery. Analyzing these datasets to better understand disease and discover new drugs is becoming more common. Recent open data initiatives in basic and clinical research have dramatically increased the types of data available to the public. The past few years have witnessed successful use of big data in many sectors across the whole drug discovery pipeline. In this review, we will highlight the state of the art in leveraging big data to identify new targets, drug indications, and drug response biomarkers in this era of precision medicine.

Abstract: Large-scale computing technologies have enabled high-throughput virtual screening involving thousands to millions of drug candidates. It is not trivial, however, for biochemical scientists to evaluate the technical alternatives and their implications for running such large experiments. Besides experience with the molecular docking tool itself, the scientist needs to learn how to run it on high-performance computing (HPC) infrastructures, and understand the impact of the choices made. Here, we review such considerations for a specific tool, AutoDock Vina, and use experimental data to illustrate the following points: (1) an additional level of parallelization increases virtual screening throughput on a multi-core machine; (2) capturing of the random seed is not enough (though necessary) for reproducibility on heterogeneous distributed computing systems; (3) the overall time spent on the screening of a ligand library can be improved by analysis of factors affecting execution time per ligand, including number of active torsions, heavy atoms and exhaustiveness. We also illustrate differences among four common HPC infrastructures: grid, Hadoop, small cluster and multi-core (virtual machine on the cloud). Our analysis shows that these platforms are suitable for screening experiments of different sizes. These considerations can guide scientists when choosing the best computing platform and set-up for their future large virtual screening experiments.

Abstract: The adverse outcome pathway (AOP) concept links molecular perturbations with organism and population-level outcomes to support high-throughput toxicity (HTT) testing. International efforts are underway to define AOPs and store the information supporting these AOPs in a central knowledge base; however, this process is currently labor-intensive and time-consuming. Publicly available data sources provide a wealth of information that could be used to define computationally predicted AOPs (cpAOPs), which could serve as a basis for creating expert-derived AOPs in a much more efficient way. Computational tools for mining large datasets provide the means for extracting and organizing the information captured in these public data sources. Using cpAOPs as a starting point for expert-derived AOPs should accelerate AOP development. Coupling this with tools to coordinate and facilitate the expert development efforts will increase the number and quality of AOPs produced, which should play a key role in advancing the adoption of HTT testing, thereby reducing the use of animals in toxicity testing and greatly increasing the number of chemicals that can be tested.