Code Monkey home page Code Monkey logo

maxabpp's Introduction

maxabpp

R package for augmented visualization of MaxQuant identification and label-free quantitation data in peptide-centric competitive activity-based protein profiling or other untargeted proteomics platforms.

Visualizing ABPP Data

Volcano Plot

Conventional Venn Diagram of Three Protein ID lists

Conventional Venn Diagram of Four Protein ID lists

Target Diagram of Three Protein ID lists

Target Diagram of Four Protein ID lists

Getting Started

How to install this package?

  1. You need to install the devtools package. You can do this from CRAN. Invoke R and then type
install.packages("devtools")
  1. Load the devtools package.
library(devtools)
  1. Install this package directly from github.
install_github("devradiumking/maxabpp")

Usage

  1. Install dependent R packages.
install.packages("tidyverse")
install.packages("stringdist")
install.packages("stringr")
install.packages("ggplot2")
install.packages("ggrepel")
install.packages("RColorBrewer")
install.packages("grid")
  1. Load dependent R packages.
library(tidyverse)
library(stringdist)
library(stringr)
library(ggplot2)
library(ggforce)
library(ggrepel)
library(grid)
library(RColorBrewer)
library(plyr)

  1. Call function pairwise_LFQ() on raw MaxQuant output ("modificationSpecificPeptides.txt" and a customized metadata file "metadata.txt" must be put in the folder set as the working directory) to obtain output1. Example metadata.txt: Raw data file name Replicate group Intensity A0288_klc_20140708m_HEK293_Redoxome_H2O2_BR1_C1 H2O2_20 Intensity A0288_klc_20140708m_HEK293_Redoxome_H2O2_BR1_C5 H2O2_2000 Intensity A0288_klc_20140708m_HEK293_Redoxome_H2O2_BR2_C1 H2O2_20 Intensity A0288_klc_20140708m_HEK293_Redoxome_H2O2_BR2_C5 H2O2_2000 Intensity A0288_klc_20140708m_HEK293_Redoxome_H2O2_BR3_C1 H2O2_20 Intensity A0288_klc_20140708m_HEK293_Redoxome_H2O2_BR3_C5 H2O2_2000
output1 <- pairwise_LFQ(
raw = read.delim("modificationSpecificPeptides.txt", header = TRUE, sep = "\t"), 
metadata = read.delim("metadata.txt", header=TRUE, sep = "\t"), 
name_probe_mod = c("Mod"), 
max_each_mod = 1, 
max_total_mods = 1, 
quantitation_level = "peptide", 
background_check = FALSE,
normalize_to = "sum_all")

Note: Multiple modification forms of a single chemical probe can be used as name_probe_mod = c("Mod1", "Mod2"). For instance, original (+ probe mass) and hydrolyzed (+ probe mass + 18 Da). These should be previously set on MaxQuant.

  1. Call function append_ec_sites() on output1 to obtain output2, for example:
output2 <- append_ec_sites(output1, quantitation_level = "peptide")
  1. Call function plot_volcano(), on output2 to obtain a volcano plot, for example:
plot_volcano(output2, "InhibitorHigh _vs_ InhibitorLow _log2fold_change", "InhibitorHigh _vs_ InhibitorLow _-log10p-value", xlim = c(-6, 2), ylim = c(0, 5), "Gene Names", 1, 0, "InhibitorName/ProbeName")
  1. New feature of v1.1, you can plot all volcano plots by calling multi_volcano_plots() functions
multi_volcano_plots(raw = raw, meta = meta, name_probe_mod = c("Mod"),
                    max_each_mod = 1, max_total_mods = 1, quantitation_level = "peptide" , background_check = FALSE,
                    xlim = c(-10, 3), ylim = c(0, 5), label_col_name = "Gene Names", pCutoff = 0.05, FCcutoff = 0)
  1. New features of v2.3: visualization of identified proteins groups from MaxQuant proteinGroups.txt with Venn Diagram and Target Diagram of tiered intersection. User-renamed proteinGroups.txt files must be put in the designated folder (default folder name is "proteinGroups"). Create one if needed. Call the functions below to make the plots.
setList <- make_proteinGroups_setList(folderName = "proteinGroups")
plot_Max_Venn(Max_Venn(setList, IndividualAnalysis = FALSE))
plot_target(make_tiers(setList))

Citation

maxabpp was developed at the Yao Lab at Chemistry Department, University of Connecticut

If you use this package please cite as:

Lei Wang and Xudong Yao (2020). maxabpp: R package for augmented visualization of peptide-centric competitive activity-based protein profiling data from MaxQuant protein identification and label-free quantitation output. package version 2.5. https://github.com/devradiumking/maxabpp

maxabpp's People

Contributors

devradiumking avatar

Watchers

 avatar

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    ๐Ÿ–– Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. ๐Ÿ“Š๐Ÿ“ˆ๐ŸŽ‰

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google โค๏ธ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.