title: "longevityDrugs: Database of Longevity-Associated Drugs (LAD)" author: "Authors: Thomas Girke, Tyler Backman, Dan Evans" date: "Last update: 25 May, 2016" package: "longevityDrugs 1.0.0" output: BiocStyle::html_document: toc: true toc_depth: 3 fig_caption: yes
fontsize: 14pt bibliography: bibtex.bib
document.addEventListener("DOMContentLoaded", function() { document.querySelector("h1").className = "title"; }); document.addEventListener("DOMContentLoaded", function() { var links = document.links; for (var i = 0, linksLength = links.length; i < linksLength; i++) if (links[i].hostname != window.location.hostname) links[i].target = '_blank'; });This vignette is part of the NIA funded Longevity Genomics project. For more information on this project please visit its website here. The GitHub repository of the corresponding R package is available here and the most recent version of this vignette can be found here.
The data package longevityDrugs
contains molecular structures and annotation information
of longevity-associated drugs (LADs). An iterative process will be used to assemble, update and
curate this data set long-term [@Backman2011-uw; @Cao2008-zo]. Most small molecules will be identified by integrating several
large-scale community data sources such as data from drug databases (e.g.
Drugbank), bioassay databases (e.g. PubChem Bioassay, ChEMBL), literature,
drug-related gene expression fingerprints (e.g. LINCS), genetic (e.g. GWAS) and
phenotype resources. This will also include nearest neighbors of LADs sharing
with them structural, physicochemical or biological properties.
The R software for running longevityDrugs
and longevityTools
can be downloaded from CRAN. The longevityTools
package can be installed from the R console using the following biocLite
install command.
source("http://bioconductor.org/biocLite.R") # Sources the biocLite.R installation script
biocLite("tgirke/longevityTools", build_vignettes=FALSE, dependencies=FALSE) # Installs package from GitHub
biocLite("tgirke/longevityDrugs", build_vignettes=FALSE, dependencies=FALSE)
library(RSQLite); library(ChemmineR); library(longevityDrugs)
## Loading required package: DBI
mypath <- system.file("extdata", "cmap.db", package="longevityDrugs")
conn <- initDb(mypath)
results <- getAllCompoundIds(conn)
sdfset <- getCompounds(conn, results, keepOrder=TRUE)
sdfset
## An instance of "SDFset" with 1309 molecules
plot(sdfset[1:4], print=FALSE)
as.data.frame(datablock2ma(datablock(sdfset)))[1:4,]
## instance_id batch_id cmap_name INN1 concentration..M. duration..h. cell2 array3
## 201 1 1 metformin INN 1.00e-05 6 MCF7 HG-U133A
## 202 21 2 phenformin INN 1.00e-05 6 MCF7 HG-U133A
## 203 22 2 phenyl biguanide 1.00e-05 6 MCF7 HG-U133A
## 204 23 2 valproic acid INN 1.00e-03 6 MCF7 HG-U133A
## perturbation_scan_id vehicle_scan_id4 scanner vehicle vendor catalog_number
## 201 EC2003090503AA EC2003090502AA HP GeneArray Scanner medium Sigma-Aldrich D5035
## 202 EC2003091104AA EC2003091102AA HP GeneArray Scanner medium Sigma-Aldrich P7045
## 203 EC2003091105AA EC2003091102AA HP GeneArray Scanner medium Sigma-Aldrich P19906
## 204 EC2003091106AA EC2003091102AA HP GeneArray Scanner medium Sigma-Aldrich P4543
## catalog_name SOURCE_DRUG UNIPROT P_SCORE DIRECTIONALITY
## 201 1,1-dimethylbiguanide hydrochloride Metformin ZN396_HUMAN 1.53601286 Stimulatory
## 202 phenformin hydrochloride Phenformin KI2L3_HUMAN -6.33081681 Inhibitory
## 203 1-phenylbiguanide hydrochloride PHENYL BIGUANIDE RHOC_HUMAN -4.21084571 Inhibitory
## 204 2-propylpentanoic acid Valproic Acid GSK3B_HUMAN -58.47057966 Inhibitory
## PUBCHEM_ID DRUGBANK_ID DRUGBANK_GROUP ATCCODES chembank_id chembank_name
## 201 CID4091 DB00331 approved A10BD11|A10BA02 1714 metformin
## 202 CID8249 DB00914 approved|withdrawn A10BA01 1018627 phenformin
## 203 CID4780 32656 phenylbiguanide
## 204 CID3121 DB00313 approved|investigational N03AG01 471 valproic acid
## match_distance smiles
## 201 0 CN(C)C(=N)NC(=N)N
## 202 0 NC(=N)NC(=N)NCCc1ccccc1
## 203 0 NC(=N)NC(=N)Nc1ccccc1
## 204 0 CCCC(CCC)C(=O)O
myfeat <- listFeatures(conn)
feat <- getCompoundFeatures(conn, results, myfeat)
feat[1:4,]
## compound_id aromatic cansmi cansmins formula hba1 hba2 hbd
## 1 201 0 CN(C(=N)NC(=N)N)C CN(C(=N)NC(=N)N)C C4H11N5 5 5 4
## 2 202 1 N=C(NC(=N)N)NCCc1ccccc1 N=C(NC(=N)N)NCCc1ccccc1 C10H15N5 5 5 5
## 3 203 1 N=C(Nc1ccccc1)NC(=N)N N=C(Nc1ccccc1)NC(=N)N C8H11N5 5 5 5
## 4 204 0 CCCC(C(=O)O)CCC CCCC(C(=O)O)CCC C8H16O2 2 2 1
## inchi logp
## 1 InChI=1S/C4H11N5/c1-9(2)4(7)8-3(5)6/h1-2H3,(H5,5,6,7,8) 0.2565
## 2 InChI=1S/C10H15N5/c11-9(12)15-10(13)14-7-6-8-4-2-1-3-5-8/h1-5H,6-7H2,(H6,11,12,13,14,15) 1.9181
## 3 InChI=1S/C8H11N5/c9-7(10)13-8(11)12-6-4-2-1-3-5-6/h1-5H,(H6,9,10,11,12,13) 1.8800
## 4 InChI=1S/C8H16O2/c1-3-5-7(6-4-2)8(9)10/h7H,3-6H2,1-2H3,(H,9,10) 2.2874
## mr mw ncharges nf r2nh r3n rcch rcho rcn rcooh rcoor rcor rings rnh2 roh ropo3 ror
## 1 36.9285 129.1636 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0
## 2 61.3212 205.2596 0 0 2 0 0 0 0 0 0 0 1 0 0 0 0
## 3 53.2452 177.2064 0 0 2 0 0 0 0 0 0 0 1 0 0 0 0
## 4 42.3418 144.2114 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0
## title tpsa
## 1 metformin 88.99
## 2 phenformin 97.78
## 3 phenyl biguanide 97.78
## 4 valproic acid 37.30
This project is funded by NIH grant U24AG051129 awarded by the National Intitute on Aging (NIA).
sessionInfo()
## R version 3.2.2 (2015-08-14)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: CentOS Linux 7 (Core)
##
## locale:
## [1] C
##
## attached base packages:
## [1] stats graphics utils datasets grDevices methods base
##
## other attached packages:
## [1] longevityDrugs_1.0.0 ChemmineR_2.22.3 RSQLite_1.0.0 DBI_0.3.1
## [5] ggplot2_2.0.0 longevityTools_1.0.6 BiocStyle_1.8.0
##
## loaded via a namespace (and not attached):
## [1] Rcpp_0.12.3 knitr_1.12.3 magrittr_1.5 munsell_0.4.3 colorspace_1.2-6
## [6] rjson_0.2.15 stringr_1.0.0 plyr_1.8.3 tools_3.2.2 grid_3.2.2
## [11] gtable_0.1.2 htmltools_0.3 yaml_2.1.13 digest_0.6.9 formatR_1.2.1
## [16] bitops_1.0-6 codetools_0.2-14 RCurl_1.95-4.7 evaluate_0.8 rmarkdown_0.9.2
## [21] stringi_1.0-1 scales_0.3.0
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.