计算机辅助药物设计通过在分子水平预测配体与靶蛋白的结合模式来加速先导化合物发现。分子对接可以快速搜索配体在蛋白活性口袋内的结合构象并给出打分，但其打分函数通常基于简化能量模型，忽略了蛋白与配体的构象柔性以及溶剂效应，预测的结合模式需要进一步验证。本项目的科学问题是：小分子抑制剂以何种构象结合在靶蛋白的活性位点？对接预测的最佳构象能否在动态环境中保持稳定？通过将分子对接与分子动力学模拟相结合，可以在纳秒至微秒时间尺度上评估配体-蛋白复合物的稳定性与结合模式的可靠性。
分子对接算法由构象搜索与打分函数两部分构成。AutoDock Vina采用基于经验的打分函数，综合范德华、静电、氢键、去溶剂化与扭转熵等贡献，以kcal/mol为单位的打分越低表示预测结合越强。对接时受体通常被视为刚性，配体的可旋转键在搜索空间内被采样，Vina输出多个排序后的结合构象供选择。分子动力学模拟则在显式溶剂中对蛋白-配体复合物进行全原子模拟，以力场描述非键相互作用。模拟中通过配体与蛋白的RMSD监控复合物稳定性：若配体从初始结合位点显著漂移，说明对接预测的构象在热涨落下不稳定；若配体RMSD在平台期保持较小值，则验证了结合模式合理。
以靶蛋白的晶体结构为受体，去除原共晶配体与水分子后加氢并分配电荷。使用AutoDock Tools设置对接盒子，盒子中心对准实验已知的活性口袋并覆盖整个结合位点，盒子边长约20埃以保证配体有足够采样空间。用AutoDock Vina对配体进行对接，输出打分最高的9个构象并按打分排序。取最优构象构建蛋白-配体复合物，用AMBER力场（配体参数经GAFF力场与RESP电荷拟合）将其置于TIP3P水盒子中，加入抗衡离子后在300 K下进行50纳秒分子动力学模拟，分析蛋白骨架与配体重原子的RMSD随时间的变化。
图10的左侧面板（a）给出Vina打分排名前9的对接构象。最佳构象打分约-9.6 kcal/mol，随后各构象打分依次上升至约-6.9 kcal/mol，相邻构象的分数差异逐渐缩小。最佳构象与第二名相差约0.7 kcal/mol，说明存在一个明显的优势结合模式，配体在活性口袋内有明确的优先取向。较低的负打分表明配体与口袋之间形成了较强的范德华与氢键互补，对接构象的聚类分析也显示前几个构象在空间上彼此接近，均占据相同的结合子口袋，进一步增强了结合模式预测的可信度。
图10的右侧面板（b）给出最优对接构象在50纳秒分子动力学模拟中的稳定性。蓝色实线为蛋白骨架RMSD，曲线在约5纳秒内上升并稳定于约1.9埃的平台，说明蛋白整体构象在结合配体后保持稳定，没有发生显著形变或去折叠。红色虚线为配体重原子RMSD，快速上升至约0.9埃后长期保持在该水平，说明配体在热涨落中仅发生小幅的取向调整而未脱离结合位点。配体RMSD明显小于蛋白RMSD表明配体被活性口袋稳定锚定，其结合模式在动力学时间尺度上具有良好的一致性，验证了对接预测结果可靠。RMSD曲线没有出现配体跳跃式位移的迹象，排除了配体从口袋解离或结合到非特异位点的可能。
两面板结合，左面板给出了结合模式的对接打分排序，右面板提供了复合物动力学的稳定性验证，共同完成了分子对接与动力学模拟的联合分析，为后续结合自由能计算与先导化合物优化提供了可靠起点。
本项目以AutoDock Vina对接预测了配体与靶蛋白的结合构象，最佳打分-9.6 kcal/mol，并通过50纳秒分子动力学模拟验证了复合物的稳定性，蛋白骨架RMSD约1.9埃、配体RMSD约0.9埃。案例图以对接打分柱状图与复合物RMSD双面板呈现。后续可拓展：(1)采用MM-PBSA或自由能微扰方法定量计算结合自由能；(2)对系列衍生物开展虚拟筛选与构效关系分析；(3)延长模拟时间并结合增强采样考察配体的结合-解离路径。
