Traduction de l'ADN en séquence d'acides aminés

<sxh python; title : dictionary_adn_protein.py> #!/usr/bin/env python # -*- coding: iso-8859-1 -*-

# python code to translante DNA sequences into proteins # traduction de l'ADN en séquence d'acides aminés (protéine)

def translate_dna(sequence): # to translate a DNA sequence in equivalent protein 

  # see http://fr.wikipedia.org/wiki/Code_g%C3%A9n%C3%A9tique
  # http://en.wikipedia.org/wiki/Genetic_code
  gencode = {                   # définition d'un dictionnaire dont chaque clé (un codon de 3 nucléotides) donne l'acide aminé correpondant
  'ATA':'I', 'ATC':'I', 'ATT':'I', 'ATG':'M',  
  'ACA':'T', 'ACC':'T', 'ACG':'T', 'ACT':'T',  
  'AAC':'N', 'AAT':'N', 'AAA':'K', 'AAG':'K',  
  'AGC':'S', 'AGT':'S', 'AGA':'R', 'AGG':'R',  
  'CTA':'L', 'CTC':'L', 'CTG':'L', 'CTT':'L',  
  'CCA':'P', 'CCC':'P', 'CCG':'P', 'CCT':'P',  
  'CAC':'H', 'CAT':'H', 'CAA':'Q', 'CAG':'Q',  
  'CGA':'R', 'CGC':'R', 'CGG':'R', 'CGT':'R',  
  'GTA':'V', 'GTC':'V', 'GTG':'V', 'GTT':'V',  
  'GCA':'A', 'GCC':'A', 'GCG':'A', 'GCT':'A',  
  'GAC':'D', 'GAT':'D', 'GAA':'E', 'GAG':'E',  
  'GGA':'G', 'GGC':'G', 'GGG':'G', 'GGT':'G',  
  'TCA':'S', 'TCC':'S', 'TCG':'S', 'TCT':'S',  
  'TTC':'F', 'TTT':'F', 'TTA':'L', 'TTG':'L',  
  'TAC':'Y', 'TAT':'Y', 'TAA':'_', 'TAG':'_',  
  'TGC':'C', 'TGT':'C', 'TGA':'_', 'TGG':'W',
  }
  print sequence  
  proteinseq = ''  
  for n in range(0,len(sequence),3):  
      if gencode.has_key(sequence[n:n+3]) == True:  # check that key exists in dictionary
          proteinseq += gencode[sequence[n:n+3]]  
  return proteinseq

# DNA sequences : GenBank : http://www.ncbi.nlm.nih.gov/

se='GTGACTTTGTTCAACGGCCGCGGTATCCTAACCGTGCGAAGGTAGCGTAATCACTTGTTCTTTAAATAAGGACTAGTATGAATGGCATCACGAGGGCTTTACTGTCTCCTTTTTCTAATCAGTGAA' #se='GTGACTTTG' print len(se) print translate_dna(se)

# extensions : lire la séquence à partir d'un fichier # générer des mutations aléatoires de l'ADN et examiner les modifications éventuelles de la protéine # faire correspondre des noms complets ou abbréviations sur 3 lettres des acides aminés # générer l'ARN correspondant </sxh>

références