如何用Python从DNA链分离蛋白密码子？代码返回‘无蛋白’问题排查

你的代码存在几个关键逻辑错误，导致始终返回"No protein"：

• 第一个循环完全无效：第一个for i in range(0,len(seqDNA),3)循环里，每次都重新执行pi=seqDNA.find(codaoi)，最后pi只是ATG第一次出现的索引，循环结束后i会变成序列末尾的索引值，直接打乱了后续判断逻辑。
• 条件判断逻辑错误：第二个循环里的if条件要求同一个3碱基片段里同时存在起始密码子和终止密码子，这根本不可能（起始是ATG，终止是TAG/TAA/TGA，没有重叠可能），所以第一次循环就会进入else分支返回"No protein"。
• 终止密码子查找逻辑错误：就算条件成立，seqDNA.find(codaof1)会返回整个序列中第一个终止密码子的位置，不管它是否在起始密码子之后，这会导致切割出错误片段。

下面是实现需求的正确代码：

def isolarprot(seqDNA):
    start_codon = "ATG"
    stop_codons = {"TAG", "TAA", "TGA"}
    
    # 找到第一个起始密码子的位置
    start_pos = seqDNA.find(start_codon)
    if start_pos == -1:
        return "No protein"
    
    # 从起始密码子之后，按3碱基步长查找终止密码子
    current_pos = start_pos + 3  # 跳过起始密码子本身
    while current_pos <= len(seqDNA) - 3:
        current_codon = seqDNA[current_pos:current_pos+3]
        if current_codon in stop_codons:
            # 返回起始密码子之后到终止密码子之前的片段
            return seqDNA[start_pos+3 : current_pos]
        current_pos += 3
    
    # 没找到有效的终止密码子
    return "No protein"

代码说明 #

1. 先定位第一个起始密码子ATG的位置，如果找不到直接返回"No protein"。
2. 从起始密码子的下一个3碱基位置开始，按密码子的步长（每次+3）遍历序列，检查每个密码子是否是终止密码子。
3. 找到第一个在起始密码子之后的终止密码子后，返回起始密码子结束位置到终止密码子开始位置之间的片段（也就是示例中的YYYYYY）。
4. 如果遍历完序列都没找到有效的终止密码子，返回"No protein"。

测试示例：

print(isolarprot("XXXATGYYYYYYTAGXXX"))  # 输出: YYYYYY

内容的提问来源于stack exchange，提问作者Alícia Magalhães