---
title: "longevityDrugs: Database of Longevity-Associated Drugs (LAD)"
author: "Authors: Thomas Girke, Tyler Backman, Dan Evans"
date: "Last update: `r format(Sys.time(), '%d %B, %Y')`"
package: "`r pkg_ver('longevityDrugs')`"
output:
BiocStyle::html_document:
toc: true
toc_depth: 3
fig_caption: yes
fontsize: 14pt
bibliography: bibtex.bib
---
```{r style, echo = FALSE, results = 'asis'}
BiocStyle::markdown()
options(width=100, max.print=1000)
knitr::opts_chunk$set(
eval=as.logical(Sys.getenv("KNITR_EVAL", "TRUE")),
cache=as.logical(Sys.getenv("KNITR_CACHE", "TRUE")))
```
```{r setup, echo=FALSE, messages=FALSE, warnings=FALSE}
suppressPackageStartupMessages({
library(longevityTools)
library(ggplot2) })
```
# Introduction
This vignette is part of the NIA funded Longevity Genomics project. For more information on this project please visit its
website [here](http://www.longevitygenomics.org/projects/). 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.
[Back to Table of Contents]()
# Getting Started
## Installation
The R software for running [_`longevityDrugs`_](https://github.com/tgirke/longevityDrugs) and [_`longevityTools`_](https://github.com/tgirke/longevityTools) can be downloaded from [_CRAN_](http://cran.at.r-project.org/). The _`longevityTools`_ package can be installed from the R console using the following _`biocLite`_ install command.
```{r install, eval=FALSE}
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)
```
[Back to Table of Contents]()
## Loading of required packages
```{r documentation, eval=TRUE}
library(RSQLite); library(ChemmineR); library(longevityDrugs)
```
[Back to Table of Contents]()
# Load database
```{r load_database, eval=TRUE}
mypath <- system.file("extdata", "cmap.db", package="longevityDrugs")
conn <- initDb(mypath)
```
# Query database
## Retrieve compound structures
```{r query_structures, eval=TRUE}
results <- getAllCompoundIds(conn)
sdfset <- getCompounds(conn, results, keepOrder=TRUE)
sdfset
plot(sdfset[1:4], print=FALSE)
as.data.frame(datablock2ma(datablock(sdfset)))[1:4,]
```
## Retrieve compound properties
```{r query_properties, eval=TRUE}
myfeat <- listFeatures(conn)
feat <- getCompoundFeatures(conn, results, myfeat)
feat[1:4,]
```
[Back to Table of Contents]()
# Funding
This project is funded by NIH grant U24AG051129 awarded by the National Intitute on Aging (NIA).
[Back to Table of Contents]()
# Version information
```{r sessionInfo}
sessionInfo()
```
[Back to Table of Contents]()
# References