NAME

Bio::DB::NextProt - Object interface to NextProt REST API.

SYNOPSIS

my $np = Bio::DB::NextProt->new();

my @result_1 = $np->search_cv(-query => "kinase#");

my @result_2 = $np->get_isoform_info(-id => "NX_O00142-2");

my @result_3 = $np->get_protein_cv_info(-id => "PTM-0205", -format => "html");

DESCRIPTION

The module allows the dynamic retrieval of information from the NextProt Database through its API service. All the information below was extracted from the API webpage. For the moment the results obtained from the API are in pure HTML, XML or JSON, so you will have to parse them yourself.

Search functionalities

Search Protein

Search proteins matching the query or search proteins for which the filter is true. Available filter values are: structure, disease, expression, mutagenesis or proteomics. Note: only one filter parameter at a time is possible for the moment.

@result = $np->search_protein(-query => "kinase");
@result = $np->search_protein(-query => "kinase", -filter => "disease");

Control Vocabulary Terms

Search control vocabulary terms matching the query or search control vocabulary terms in the category specified by the filter. Available filter values are: enzyme, go, mesh, cell, domain, family, tissue, metal, pathway, disease, keyword, ptm, subcell. Note: only one category at a time is possible.

@result = $np->search_cv(-query => "colon");
@result = $np->search_cv(-query => "colon", -filter => "keyword");

Format:

Output format maybe in JSON (default), HTML or XML.

@result = $np->search_protein(-query => "kinase", -filter => "disease", -format => "html");

Find information by protein entry

Protein ID

ID is neXtProt identifier. Retrieve gene name as well as main identifier and isoform sequences

@result = $np->get_protein_info(-query => "NX_P13051");

Post-translational modifications

For each isoform of the specified entry, retrieve all post-translational modifications.

@result = $np->get_protein_info(-query => "NX_P13051", -retrieve => "ptm");

Variant

For each isoform of the specified entry, retrieve all variants.

@result = $np->get_protein_info(-query => "NX_P13051", -retrieve => "variant");

Localisation

For each isoform of the specified entry, retrieve all subcellular location.

@result = $np->get_protein_info(-query => "NX_P13051", -retrieve => "localisation");

Expression

Retrieve all expression information by tissue for the specified entry.

@result = $np->get_protein_info(-query => "NX_P13051", -retrieve => "expression");

Format:

Output format maybe in JSON (default), HTML or XML.

@result = $np->get_protein_info(-query => "NX_P13051", -retrieve => "expression", -format => "html");

Find information by isoform

Protein ID

ID is neXtProt identifier. Retrieve gene name as well as main identifier and isoform sequences

@result = $np->get_isoform_info(-query => "NX_O00142-2");

Post-translational modifications

For each isoform of the specified entry, retrieve all post-translational modifications.

@result = $np->get_isoform_info(-query => "NX_P01116-2", -retrieve => "ptm");

Variant

For each isoform of the specified entry, retrieve all variants.

@result = $np->get_isoform_info(-query => "NX_P01116-2", -retrieve => "variant");

Localisation

For each isoform of the specified entry, retrieve all subcellular location.

@result = $np->get_isoform_info(-query => "NX_P01116-2", -retrieve => "localisation");

Find information by controlled vocabulary term

Protein ID

ID is neXtProt identifier. Retrieve the accession, the name and the category of the CV term.

@result = $np->get_protein_cv_info(-query => "PTM-0205");

Protein List

List all the proteins associated with the term in neXtProt.

@result = $np->get_protein_cv_info(-query => "PTM-0205", -retrieve => "proteins");

Format:

Output format maybe in JSON (default), HTML or XML.

@result = $np->get_protein_cv_info(-query => "PTM-0205", -retrieve => "proteins", -format => "html");

Retrieving Accession Lists

Allows the retrieval of all accession codes from individual chromossomes or from the entire NextProt database.

@result = $np->get_accession_list(-chromosome => "10");

@result = $np->get_accession_list(-chromosome => "all");

entries with a protein existence "at protein level" (PE 1)

@result = $np->get_accession_list(-evidence => "protein_level");

entries with a protein existence "at transcript level" (PE 2)

@result = $np->get_accession_list(-evidence => "transcript_level");

entries with a protein existence "by homology" (PE 3)

@result = $np->get_accession_list(-evidence => "homology");

entries with a protein existence "predicted" (PE 4)

@result = $np->get_accession_list(-evidence => "predicted");

entries with a protein existence "uncertain" (PE 5)

@result = $np->get_accession_list(-evidence => "uncertain");

Customized Report Files for the HUPO Human Proteome Project (HPP)

Individual report files for each chromosomes (1 to 22, X, Y and MT)

my @list = $np->get_hpp_report(-chromosome => 10);

Annotated phosphorylated residues per chromosome

my @list = $np->get_hpp_report(-phospho => "true");

Annotated N-Acetyl residues per chromosome

my @list = $np->get_hpp_report(-nacetyl => "true");

NextProt Mapping

Mapping of neXtProt accession numbers to external resources.

Ensembl gene identifiers

@list = $np->get_mapping(-map => "ensembl_gene");

Ensembl protein identifiers

@list = $np->get_mapping(-map => "ensembl_protein");

protein ids that cannot be mapped to any isoform in neXtProt

@list = $np->get_mapping(-map => "ensembl_unmapped");

Ensembl transcript identifiers

@list = $np->get_mapping(-map => "ensembl_transcript");

transcript is considered as coding by Ensembl that cannot be mapped to any isoform in neXtProt

@list = $np->get_mapping(-map => "ensembl_transcript_unmapped");

NCBI GeneID gene accession numbers

@list = $np->get_mapping(-map => "geneid");

HGNC gene accession numbers

@list = $np->get_mapping(-map => "hgnc");

MGI mouse gene accession numbers

@list = $np->get_mapping(-map => "mgi");

NCBI RefSeq gene accession numbers

@list = $np->get_mapping(-map => "refseq");

Chromosome Report

The module also allows the programatic access to chromosome information by accessing and formatting the chr_report tables from the nextprot ftp server. The retrieved structure is a hash of hashes, being the first key the NextProt Accession Number. The internal hashes have the following values:

* Gene Name * Chromosomal position * Start position * Stop position * Protein existence * Proteomics * Antibody * 3D * Disease * Isoforms * Variants * PTMs * Description

This is how the data is representes in the hashes:

NX_A7E2V4                   {
    antibody         "yes",
    description      "Zinc finger SWIM domain-containing protein 8",
    disease          "no",
    existence        "protein level",
    has_3d           "no",
    isoforms         5,
    gene_name        "ZSWIM8",
    position         "10q22.2",
    proteomics       "yes",
    ptms             6,
    start_position   75545340,
    stop_position    75561551,
    variants         67
}

Loading the Chromosome table.

Loas all the information from tha table.

my %data = $np->get_chromosome(-chromosome => 10);

Accessing Protein information:

say $data{ZSWIM8}->{isoforms};
say $data{ZSWIM8}->{proteomics};
say $data{ZSWIM8}->{description};

Counting the number of Proteins in the Chromosome

$sum = (keys %data);
say $sum;		 

Retrieve all Gene Names from a giving Chromosome

for my $prot (keys %data) {
	say $prot;
}

FEEDBACK

Mailing Lists

User feedback is an integral part of the evolution of this and other Bioperl modules. Send your comments and suggestions preferably to one of the Bioperl mailing lists. Your participation is much appreciated.

bioperl-l@bioperl.org                  - General discussion
http://bioperl.org/wiki/Mailing_lists  - About the mailing lists

Support

Please direct usage questions or support issues to the mailing list:

bioperl-l@bioperl.org

rather than to the module maintainer directly. Many experienced and reponsive experts will be able look at the problem and quickly address it. Please include a thorough description of the problem with code and data examples if at all possible.

Reporting Bugs

Report bugs to the Bioperl bug tracking system to help us keep track the bugs and their resolution. Bug reports can be submitted via the web:

https://github.com/bioperl/bioperl-live

AUTHOR - Felipe da Veiga Leprevost

Email leprevost@cpan.org