生物医药 × 计算科学 × AI：蛋白质折叠与药物发现

1
# 从一个DNA序列开始——这是生命的"源代码"
2
dna = "ATGGCCATTGTAATGGGCCGCTGAAAGGGTGCCCGATAG"
3

4
# 计算碱基组成
5
print("🧬 DNA序列分析：")
6
print(f"   序列: {dna}")
7
print(f"   长度: {len(dna)} 个碱基")
8
print(f"   腺嘌呤 A: {dna.count('A')} ({dna.count('A')/len(dna)*100:.1f}%)")
9
print(f"   胸腺嘧啶 T: {dna.count('T')} ({dna.count('T')/len(dna)*100:.1f}%)")
10
print(f"   鸟嘌呤 G: {dna.count('G')} ({dna.count('G')/len(dna)*100:.1f}%)")
11
print(f"   胞嘧啶 C: {dna.count('C')} ({dna.count('C')/len(dna)*100:.1f}%)")
12
print(f"   GC含量: {(dna.count('G') + dna.count('C'))/len(dna)*100:.1f}%")
13

14
# 转录：DNA → mRNA
15
transcription_table = str.maketrans('ATGC', 'UACG')
16
mrna = dna.translate(transcription_table)
17
print(f"\n📝 转录结果 (mRNA): {mrna}")
18

19
# 翻译：mRNA → 蛋白质
20
# （这里用一个简化版的密码子表）
21
genetic_code = {
22
    'AUG': 'M(起始)', 'GCC': 'A', 'AUU': 'I', 'GUA': 'V',
23
    'GGG': 'G', 'CGC': 'R', 'UGA': '终止', 'UAG': '终止'
24
}
25
protein = []
26
for i in range(0, len(mrna)-2, 3):
27
    codon = mrna[i:i+3]
28
    aa = genetic_code.get(codon, '?')
29
    protein.append(aa)
30
    if aa == '终止':
31
        break
32

33
print(f"\n🧪 翻译结果 (氨基酸链): {'-'.join(protein)}")

1
# 场景：你想快速了解"mRNA疫苗的最新进展"
2
# 以前：去PubMed搜，看几十篇摘要
3
# 现在：让AI帮你读、帮你总结
4

5
question = """
6
请帮我查阅关于mRNA疫苗在肿瘤治疗领域的最新进展（2024-2026），
7
总结：
8
1. 目前处于临床试验阶段的mRNA癌症疫苗有哪些？
9
2. 与传统疗法相比，mRNA疫苗的主要优势是什么？
10
3. 目前面临的最大挑战是什么？
11
4. 未来2-3年最值得关注的突破方向？
12

13
请用通俗语言解释，让非生物专业的人也能看懂。
14
"""
15

16
# response = openai_client.chat.completions.create(
17
#     model="gpt-4o",
18
#     messages=[{"role": "user", "content": question}]
19
# )