Using Biojava 5, I am trying to convert a sequence of Amino Acids to a Nucleotide sequence. Any idea how to do this?
Asked
Active
Viewed 81 times
0
-
Welcome to SE.bioinformatics. What have you tried so far? – user438383 Sep 14 '20 at 15:24
-
Just looked in BioJava documentation - I don't think they have anything to do this. There is http://www.bioinformatics.org/sms2/rev_trans.html, but I need to write code for this - so can't use that tool – user2689782 Sep 14 '20 at 16:46
-
"convert a sequence of Amino Acids to a Nucleotide sequence". Like you've mentioned in a comment, this is called "reverse translation". Avoid the word "convert" as much as you can - it is as generic as a word can get in our context. It doesn't look like BioJava has a reverse translate provision - you might need to write some code yourself. – Ram RS Sep 14 '20 at 18:08
-
3Does this answer your question? How to translate amino acid sequences to Nucleotide sequences – Chris_Rands Sep 14 '20 at 20:22
-
Thanks Chris for that link. Interesting. I ended up using files from NCBI, that worked. – user2689782 Sep 20 '20 at 03:31
1 Answers
1
Note that unlike the reverse (nucleotide sequence to protein sequence), which is uniquely defined given a reading frame, the mapping from protein sequence to nucleotide sequence is not unique (it is not a function in the mathematical sense). The method you use to select which of the many (exponential in input sequence length) possible sequences to output depends critically on what the question is.
If the goal is to design a synthetic gene that will produce a desired protein, then you need to optimize it for expression in the organism into which the gene will be inserted, subject to the constraint of the encoded amino acid sequence.
If the goal is to search for nucleotide sequences in one or more genomes that could give rise to a certain polypeptide, there is a much better way than to do reverse translation at all. In particular, most putative protein-coding open reading frames in most sequenced genomes have already been run through forward translation code and the expected amino acid sequence encoded in them can be found in databases. So, this problem is reduced to a substring search in these databases.
biohacker
- 131
- 3
-
Very interesting! Thanks. That is what I ended up doing. Matching up FASTA files to get sequences. – user2689782 Sep 20 '20 at 03:26