做药物靶点研究、小分子药物筛选的小伙伴，大概率都遇到过同一个难题：辛辛苦苦做完分子对接，拿到看似完美的蛋白 - 配体结合构象，一深入研究就发现静态对接结果根本经不起推敲。真实生理环境中，蛋白会动态波动、配体可能脱落、复合物稳定性随时变化，仅凭静态结构根本无法判断药物结合的有效性。而能解决这一问题的分子动力学模拟，却因配置繁琐、命令行复杂、参数调试难度大、分析流程冗杂，让无数科研人员望而却步，白白耗费大量时间精力。
今天我们就以经典免疫药物靶点体系——FKBP12–他克莫司复合物（PDB：1FKJ）为实战案例，带大家零门槛上手完整的分子动力学模拟流程，轻松搞定蛋白 - 药物复合物稳定性分析、动态相互作用探究，让药物机理研究更高效、更省心。
AsterFire Go 是一个用聊天就能做科学计算的平台。分子对接、蛋白设计、动力学模拟、量子化学、DFT 计算……打字告诉它你想算什么，30 分钟内给你结果。平台致力于让科研伙伴可以轻松调用各类计算能力，不用琢磨配置集群、挑选方法、调试参数这些事，把精力留给科研本身。现阶段封装了 SPONGE（对标 GROMACS）、ABACUS（对标 VASP）、ORCA（学术用途）、ProteinMPNN、Vina、Bindcraft、RFdiffusion、文献检索总结、pdf 保真翻译等百种基础工具及其环境，均可通过自然语言调用或构建工作流。
AsterFire Go 用户交流群
产品内测期间加入官方用户群可得 1 个月专业版会员，含 4 个流动星币 + 50 GB 存储空间 + 百余种高级模型 + 全算力消费 9 折优惠。
分子对接可以预测候选小分子在蛋白结合区域中的可能构象，但对接结果本质上是一张静态结构，不能直接说明该构象在水溶液、温度和压力条件下是否能够保持，也不能反映蛋白和配体随时间发生的构象变化。
对于一个已经获得的蛋白–小分子对接复合物，我们通常还需要回答两个问题：
分子动力学模拟（molecular dynamics, MD）是进一步研究这些问题的常用方法。它通过连续计算体系中原子的运动，获得随时间变化的轨迹，并结合 RMSD、RMSF、回转半径、SASA、氢键和相互作用能等指标，对复合物的动态行为进行分析。MD 的作用是对候选结合构象进行动态检查和比较，而不是仅凭一条短轨迹直接证明高亲和力或实验活性。
本教程选用经典免疫抑制药案例 FKBP12–他克莫司。PDB 1FKJ 是人源 FKBP12 与 FK506（Tacrolimus，他克莫司；结构中配体编号为 FK5）的 1.70 Å 晶体结构。原始 PDB 非对称单元中包含一条 FKBP12 蛋白链 A（1–107）和一个 FK5 配体（A 108）。本教程保留链 A 和 FK5 108，仅删除结晶水。结构来源：https://www.rcsb.org/structure/1FKJ
经典分子动力学使用力场描述原子之间的键合和非键相互作用，根据势能函数计算原子受力，再以飞秒量级时间步长积分运动方程，从而得到体系随时间变化的构象轨迹。显式水和离子用于近似溶液环境，周期性边界条件用于减小有限模拟盒带来的表面效应。
一个常规蛋白–小分子 MD 流程通常包括结构准备、蛋白和配体参数化、建盒、加水和离子、能量最小化、平衡模拟、生产期模拟以及轨迹分析。力场和配体电荷决定相互作用的近似方式，模拟时长决定有限时间内能够采样到的状态。
本教程使用的平台工具名称为“GROMACS 2023.2_蛋白和小分子复合物体系”。该 Kit 面向由标准氨基酸组成的单蛋白或多链蛋白与单一小分子配体复合物，可自动完成结构处理、配体参数化、体系构建、GROMACS 模拟和标准化分析。
蛋白部分固定采用 amber14SB_ROC_lipid17 力场和 TIP3P 水模型；小分子通过 ACPYPE/AmberTools 生成 GAFF2 参数，并使用 AM1-BCC 电荷。支持同时上传含正确键级的原始配体 mol2 或 PDB 文件，Kit 会将其作为参考修复对接文件中的配体连接信息。
图 1：GROMACS 2023.2_蛋白和小分子复合物体系使用教程
任务完成后，用户可以直接查看 report.md、综合分析图、末态结构和处理后轨迹，同时下载 RMSD、RMSF、SASA、回转半径、温度、压力、总能量、氢键、短程 LJ/Coulomb 相互作用、RMSD–Rg 自由能形貌、MM/GBSA 结果及量化统计表。
图 2：GROMACS 传统模拟流程
在传统操作流程下，GROMACS 分子动力学模拟需要配置相应的软件环境，并准备力场、水模型及相关参数文件。完成蛋白–配体复合物的结构预处理后，需要进行拓扑生成、溶剂盒构建、离子添加和能量最小化，并通过命令行设置温度、压力、模拟步长及模拟时间等参数。计算完成后，还需要对轨迹、能量和结构变化等数据分别进行整理，并借助可视化软件开展进一步分析。对于不熟悉分子动力学理论、参数设置和命令行操作的用户，整个流程具有一定的使用门槛。
以上介绍的 GROMACS 2023.2_蛋白和小分子复合物体系，AsterFire Go 平台已将上述流程整合到工具百宝箱（Kit）模式中，用户只需上传蛋白和配体文件，填写关键模拟参数（如模拟温度、模拟时间等），即可提交计算并获取可视化结果。相比传统方式，Kit 模式省去了软件安装、环境配置和命令行操作，但仍需用户自行准备结构文件并确认计算参数。
在此基础上，平台进一步提供智算智能体（Agent）模式，将文件准备和参数配置环节也纳入自动化流程。该模式采用对话式交互，用户只需用自然语言描述模拟任务（如“对 FKBP12–他克莫司复合物进行模拟”），平台即可根据需求自动准备结构文件、配置计算参数、调用相应工具执行模拟，并返回结果。整个过程中用户无需上传任何文件，也无需填写任何参数。
图 3：智算 Agent 实现分子动力学模拟
蛋白质–小分子分子动力学模拟是研究复合物构象稳定性、动态相互作用和结构变化的常用计算方法。以 FKBP12–他克莫司复合物（PDB ID：1FKJ）为例，GROMACS 可在显式溶剂环境中模拟蛋白质与药物分子随时间的运动，并通过 RMSD、RMSF、回转半径、SASA、氢键和最小距离等指标，评估复合物的整体稳定性及关键相互作用，为结合模式分析和后续研究提供参考。
传统分子动力学模拟通常涉及软件安装、运行环境配置、复合物结构清理、配体参数化、溶剂盒构建、离子添加、能量最小化、平衡模拟、生产期模拟以及轨迹分析等多个环节，对用户的计算环境和专业经验具有较高要求。AsterFire Go 将 GROMACS 2023.2 蛋白质–小分子模拟流程封装为统一的计算 Kit：在 Kit 模式下，用户可自主上传 FKBP12–他克莫司复合物和配体结构，并控制模拟时间、温度、压力、配体电荷等关键参数；在智算 Agent 模式下，平台自研 Agent 可围绕 1FKJ 进行结构调研和输入核对，辅助完成文件准备、参数配置、Kit 调用及结果整理，从而降低分子动力学模拟的使用门槛。
通过本教程的 FKBP12–他克莫司案例，可以看到 Kit 与智算 Agent 能够协同完成从结构确认、体系构建和模拟执行到轨迹分析与结果解读的完整流程。需要注意的是，本教程采用的短时间模拟主要用于验证流程是否能够正确运行，并初步观察复合物的结构稳定性，不能直接证明他克莫司的结合亲和力、药理活性或作用机制，也不能代表体系已经得到充分采样。更可靠的结论仍需结合更长时间的重复模拟、实验结构数据及生化实验进一步验证。
