R/uniprot_get.R

In egenn/rtemis: Machine Learning and Visualization

# uniprot_get.R
# ::rtemis::
# 2022 E.D. Gennatas lambdamd.org

#' Get protein sequence from UniProt
#'
#' @param accession Character: UniProt Accession number - e.g. "Q9UMX9"
#' @param baseURL Character: UniProt rest API base URL.
#' Default = "https://rest.uniprot.org/uniprotkb"
#' @param verbose Logical: If TRUE, print messages to console
#'
#' @return List with two elements: Annotation & Sequence
#' @author E.D. Gennatas
#' @export
#' @examples
#' \dontrun{
#' matp <- uniprot_get("Q9UMX9")
#' }
uniprot_get <- function(accession = "Q9UMX9",
                        baseURL = "https://rest.uniprot.org/uniprotkb",
                        verbose = TRUE) {
  path <- paste0(baseURL, "/", accession, ".fasta")
  dat <- seqinr::read.fasta(path, seqtype = "AA")
  Annotation <- attr(dat[[1]], "Annot")
  Identifier <- gsub(" .*", "", Annotation)
  if (verbose) {
    msg2("Got:", hilite(Annotation))
  }

  list(
    Identifier = Identifier,
    Annotation = Annotation,
    Sequence = as.character(dat[[1]])
  )
} # rtemis::uniprot_get