[Bioperl-l] Degenerate primer calculation
cjfields at illinois.edu
Mon Dec 8 11:41:04 EST 2008
On Dec 8, 2008, at 9:59 AM, Samantha Thompson wrote:
> I also have another similar sequence analysis/primer problem.
> What I'd like to do is produce degenerate primers from amino acid
> What I did initially was take the codon usage table and rewrite it
> in a
> hash in perl in the form of degenerate codon usage e.g Lysine/K
> would be
> AAR, its reverse complement would be YTT. So my form then takes an
> acid sequence (derived as a consensus from multiple the alignment of
> homologous proteins) and converts them into degenerate codons and then
> that degenerate primer (actually several primers synthesised with
> different bases pooled together), in order to search for homologues to
> the protein in unsequenced organisms.
> I would like to improve this by being able to take a consensus
> more in the form of a Prosite motif (I think thats the right one) such
> as [TS]YW[RKSD] and then develop a degenerate nucleotide sequence
> corresponding to this.
> So I'm wondering if bioperl contains anything like this (both prosite
> motif format parsing and degenerate code from multiple alignments or
> such a motif), or if I need to write this myself (which I want to if
> doesn't exist already).
> Thanks again,
Bio::Tools::CodonTable reverse translates, but I don't think it
accepts patterns. Maybe a pipeline including Bio::Tools::SeqPattern?
Might be an interesting programming challenge if it isn't already set
up for that.
More information about the Bioperl-l