Python 中的密码子生成-解网

问：

我有这个代码，用于将 DNA 字符串转换为密码子列表，然后将此列表转换为具有各自氨基酸的字符串。但是，当我运行代码并且 DNA 字符串以一对核苷酸（例如 CT）而不是三联体结尾时，代码不会生成氨基酸序列。正如您在输出中看到的那样。

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)

输出：

密码子： ['GAG'， 'CGT'， 'CTG'， 'CTC'， 'CGT'， 'GTA'， 'TAA'， 'GCC'， 'ACG'， 'TCG'， 'GAG'， 'CT']

氨基酸： ''

我想找到一个解决方案，以便我可以使用字符串的最后两个核苷酸来预测下一个氨基酸，选择以这两个核苷酸开头的密码子（genetic_code），然后选择一个。我该怎么做？请提供任何建议。

python-3.x 字符串列表 DNA 序列

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]

上一个：如何打印包含部分输入的所有行？

下一个：生成具有不同长度的序列

Python 中的密码子生成

Codon generation in Python

评论