Home > Blockchain >  Codon generation in Python
Codon generation in Python

Time:09-12

I have this code for converting a DNA string into a list of codons and after that to convert this list into a string with their respective aminoacids. However, when I ran the code and the DNA string ends in a pair of nucleotides(like CT for instance) and not in a triplet, then the code doesn't generate the aminoacid sequence. As you can see in the output.

from collections import defaultdict
from collections import Counter


dna_sequence = "GAGCGTCTGCTCCGTGTATAAGCCACGTCGGAGCT"
codons = [dna_sequence[i:i 3]
      for i in range (0, len(dna_sequence), 3)]
      print(codons)

genetic_code = {
   "GCG" :"A","GCA" :"A","GCT" :"A","GCC" :"A",
   "AGG" :"R","AGA" :"R","CGG" :"R", "CGA" :"R","CGT" :"R","CGC" :"R",
   "AAT" :"N","AAC" :"N",
   "GAT" :"D", "GAC":"D", "TGT" :"C","TGC" :"C",
   "TGA" :"*","TAG" :"*","TAA" :"*", # * Stop codon
   "CAG" :"Q","CAA" :"Q",
   "GAG" :"E","GAA" :"E",
   "GGG" :"G","GGA" :"G","GGT" :"G","GGC" :"G",
   "CAT" :"H","CAC" :"H",
   "ATA" :"I","ATT" :"I","ATC" :"I",
   "TTG" :"L","TTA" :"L","CTG" : "L","CTA" :"L","CTT" :"L","CTC" :"L",
   "AAG" :"K","AAA" :"K",
   "ATG" :"M" , # Start codon
   "TTT" :"F" ,"TTC" :"F" ,
   "CCG" :"P" ,"CCA" :"P" ,"CCT" :"P" ,"CCC" :"P" ,
   "AGT" :"S" ,"AGC" :"S" ,"TCG" :"S" ,"TCA" :"S" ,"TCT" :"S" ,"TCC" :"S" ,
   "ACG" :"T" ,"ACA" :"T" ,"ACT" :"T" ,"ACC" :"T" ,
   "TGG" :"W" ,
   "TAT" :"Y" ,"TAC" :"Y" ,
   "GTG" :"V" ,"GTA" :"V" ,"GTT" :"V" ,"GTC" :"V" 


}
def codon_seq(seq):
tmpList = []
for i in range(0, len(seq) - 2, 3):
    if genetic_code [seq[i:i   3]]:
        tmpList.append(seq[i:i   3])
        print(tmpList)

def amino_seq(seq):
 protein =""
 if len(seq)%3 == 0:
  for i in range(0, len(seq), 3): 
   codon = seq[i:i   3]
   protein = genetic_code[codon]

  return protein

print("Aminoacids: ")
amino_seq(dna_sequence)

OUTPUT:

Codons: ['GAG', 'CGT', 'CTG', 'CTC', 'CGT', 'GTA', 'TAA', 'GCC', 'ACG', 'TCG', 'GAG', 'CT']

Aminoacids: ''

I'd like to find a solution in order that I can use the two last nucleotides of the string to predict the next aminoacid, selecting the codons (genetic_code) that start with these two nucleotides and choosing one ramdonly. How could I do that? Any advice please.

CodePudding user response:

You can define a function random_protein to get all the codons that start with the given prefix and random.choice one of those and use that if the len of the current codon is != 3.

import random
def random_protein(prefix):
    codons = [c for c in genetic_code if c.startswith(prefix)]
    return genetic_code[random.choice(codons)]

dna_sequence = "GAGCGTCTGCTCCGTGTATAAGCCACGTCGGAGCT"
codons = [dna_sequence[i:i 3] for i in range (0, len(dna_sequence), 3)]
proteins = [genetic_code[c] if len(c) == 3 else random_protein(c) for c in codons]
  • Related