align_to_unip: Align protein segnent around the hotspot to the...

Description Usage Arguments Value Author(s) Examples

View source: R/align_to_unip.R

Description

This function alignes the Ensembl protein region with a hotspot to the UniProt sequence. The Ensembl region encompasses 15 amino acids where the hotspot is in the middle. If the hotspot was at the protein start or end the region is still 15 amino acids long, but the hotspot position is shifted.

Usage

1
  align_to_unip(ens.seq, uni.seq, ensembl_mut_position)

Arguments

ens.seq

AAString object with the Ensembl protein sequence corresponding to the representative transcript.

uni.seq

AAString with the UniProt sequence for the identifier matching the Ensembl gene name.

ensembl_mut_position

Numeric vector defining the hotspot position in the Ensembl sequence, i.e. in the ens.seq

Value

Returns a list where the first element is a character vector defining whether there was a significant alignment and the second element provides the hotspot position in the UniProt sequence.

Author(s)

Marija Buljan <buljan@imsb.biol.ethz.ch> Peter Blattmann <blattmann@imsb.biol.ethz.ch>

Examples

1
2
3
4
5
6
7
library(Biostrings)

ens.seq <- AAString("MDLSALREEVQNVINAMQKILECPICLELIKEPVSTKCDHIFCKFCMLK")
uni.seq <- AAString("MDLSALRVEEVQNVINAMQKILECPICLELIKEPVSTKCDHIFCKFCMLA")
ensembl_mut_position <- 25

align_to_unip(ens.seq, uni.seq, ensembl_mut_position)

DominoEffect documentation built on Nov. 8, 2020, 5:46 p.m.