-> 启动多尺度增强采样算法…并行遍历构象空间…
-> 计算中…预计剩余时间:…
进度条飞速跳动。
12.8秒!
仅仅 12.8秒,计算完成!
-> 模拟完成。总有效采样时间等效:1.52 μs。
-> 生成结合模式分析报告、关键相互作用动态图、结合自由能估算值。
林枫立刻查看结果。报告详尽得令人发指:
稳定结合模式: 清晰地显示A-73分子并非如其设计者预想的那样嵌入某个深口袋。
而是以一种巧妙的“表面锚定”方式,与蛋白两个柔性loop区(环区)的关键残基(K311, I354)形成稳定的疏水簇和氢键网络。
关键相互作用: 指出了一个之前被忽略的、与残基D283形成的水分子介导的氢键,此氢键寿命长达 ~850 ns,对结合稳定性贡献显着。
结合自由能 (ΔG): 估算值为 -9.8 kcal/mol(误差范围 ±0.3),表明结合能力尚可,但有优化空间。
不稳定性根源: 明确指出了分子结构中一个甲氧基的旋转自由度太大,在模拟中频繁撞击蛋白表面,造成局部构象扰动,是影响结合稳定性和特异性的主要负面因素。
优化建议: 系统甚至基于此分析,推演了三个具体的结构修饰方案,并给出了每个方案的预期结合自由能提升预测(分别约为 +1.2, +0.7, +1.5 kcal/mol)。
林枫将这份核心发现和优化建议,整理成一份简洁明了、图文并茂的技术摘要,隐去了所有涉及计算方法和速度的信息,通过“火鸡”的渠道,匿名回复给了提问者。
本小章还未完,请点击下一页继续阅读后面精彩内容!
---
位于华东某市高新区的“安华生物科技”研发中心内,药物化学部门负责人张建明博士正对着电脑屏幕上的分子动力学模拟进度条发愁。
已经运行了两周多,才完成了不到0.4微秒,数据还噪声重重。
他几天前在某个小圈子技术论坛上,偶然看到了关于“幽灵愿望池”的传说,抱着死马当活马医的心态,匿名提交了关于A-73分子的问题。
他并没抱太大希望,只当是个心理安慰。
然而,就在提交问题后的第二天,他收到了加密回复!
当他点开那份技术摘要时,整个人如同被电流击中,僵在了座位上。
“这……这不可能!”他喃喃自语。
报告里指出的稳定结合模式,与他团队基于晶体结构猜测的完全不同!
那个水分子介导的氢键,他们从未在分析中考虑到!
而那个被指责为“祸根”的甲氧基,正是他一位手下爱将的“得意设计”!
更让他震惊的是报告数据的详实程度和内在逻辑的自洽性。