▊药物研发中结合自由能计算应用现状药物分子通过对靶蛋白的识别与结合作用能够调控靶蛋白功能进而实现治疗疾病的效果。蛋白质的许多关键生理和药理活动是通过与小分子相互作用来实现比如酶的催化特性是由其与底物的相互作用所体现的。在药物研发中众多药物和其他生物活性分子的功效实质体现在它们与生物大分子受体的相互协作中这意味着药物与受体间的结合强度直接影响到药物的活性表现。结合自由能是评价药物与受体相互作用强度的关键参数之一它直接关系到药物分子能否有效结合至靶标受体并进一步影响其生物活性和药理效应。高精度的结合自由能计算一直都是现代基于结构的药物发现工作中关键的组成部分这类计算通常以统计力学为基础譬如基于分子动力学模拟的自由能计算。自由能微扰(FEP)方法是其中最常用的一种计算方法它能以与实验精度(1kcal/mol)相当的计算精度预测配体和受体之间的相对结合自由能(RBFE)和绝对结合自由能(ABFE)该方法的核心思想就是从一个已知的体系出发通过一系列微小变化的累积变成另一个体系其中每一步的变化都是一次分子动力学模拟最后再利用自由能估算算法从多次的分子动力学模拟中得出自由能的变化。目前有效的应用FEP展开药物研发并不简单。即便不考虑初始结构的合理性仅使用FEP方法进行模拟就需要考虑诸多要素以RBFE计算为例需要考虑的要素就包括拓扑方案、原子间映射、微扰区域、微扰窗口、软核参数等技术要素以及电荷微扰和骨架跃迁等场景要素。这些要素的组合决定了预测的速度、精度、通量和成本。基于FEP方法的开源软件包诸如PyAutoFEP往往需要进行精细复杂的配置以及稳固的计算环境构建然而对于许多化合物pair而言建立一个既能保持稳定性又能提供可靠计算结果的体系并非易事这使得它们难以满足工业化应用级别的标准要求。此外FEP方法存在诸多经验性操作无疑给那些意图通过定制FEP模拟协议以精确计算各类不同体系结合自由能的研究人员带来了显著的挑战。▊ 高效精确灵活的解决方案MaXFlow的Neo-BFE模块引入了一种名为ATM(Alchemical Transfer Method,炼金术转移法)的创新方法来执行高精度的结合自由能计算。ATM方法提供了一种通过使用坐标变换交换配体位置来计算绝对/相对结合自由能的协议。与常用的FEP方法不同ATM方法在单个溶剂化盒子中就能完成结合自由能的计算可避免配体净电荷变化的影响不仅减少了计算量还极大简化了计算流程。此外由于该方法无需像FEP那样通过开、关闭配体与受体之间的静电和非静电相互作用项来实现配体与受体之间的解耦进而模拟计算出结合自由能因此该方法能兼容任何势能函数(如可极化量子力学以及机器学习势)同时还能避免FEP执行过程中的那些经验性较强的技术和场景因素带来的问题。▊打造结合自由能计算的“灰盒”工具Neo-BFE模块的设计目标是将结合自由能计算作为一个“灰盒” 即允许用户根据计算任务的实际需求对过程进行精细化的控制同时也可以直接采取它提供的“最佳实践”默认值以减少个人干预和实现自动化。因此无论是专家或非专家用户均可使用Neo-BFE模块来进行稳健、高效、准确的结合自由能计算。▊计算性能精度方面Neo-BFE支持机器学习力场和哈密顿副本交换增强采样的算法可在兼顾计算精度的同时做到充分采样速度方面经过优化后的分子动力学模拟积分器- Langevin Middle Integrator 会使得模拟更稳定在添加氢质量再分配情况下可将默认时间步长从2fs增加到4fs而不损失稳定性使得模拟所需要的GPU时间减少近两倍在推荐的参数下一次结合自由能计算能够在5-10小时内完成基于MXF的云计算平台可以调用海量资源一天完成成百上千次的自由能计算。易用性方面因为不需要像FEP那样定义复杂的原子映射关系微扰区域软核或软键势等经验性较强的参数所以和常规动力学类似简单的设置便能得到高精度的计算结果同时MaXFlow平台实现了计算流程的极大自动化可轻易的进行高通量计算灵活分配计算资源。目前Neo-BFE中的ATM方法已达到与实验精度相当的结合自由能预测值(1 kcal/mol)。在已发表的论文中研究人员通过选择Schrödinger论文中发表的数据集进行测试(包括8 个不同蛋白质系统中的330个配体对)并将ATM方法与其他常见方法(例如 FEP、Amber 和 pmx)进行比较验证它们在相对结合自由能预测方面的性能[1]。结果表明ATM 在Pearson相关性方面与更复杂的自由能扰动 (FEP) 方法的性能相匹配甚至超越。尽管其平均误差略高但考虑到ATM方法在实现上的简便性、灵活性以及在势函数通用性方面的优势这一方法在药物发现中有着巨大的应用前景和价值。表1:自由能方法的性能比较8个被测蛋白的Pearson相关性(r)和平均绝对误差(MAE) 单位kcal/mol▊MaXFlow新一代分子模拟与人工智能平台在MaXFlow找到该模块通过拖拽组件的形式搭建如下工作流选择要计算的结构设置相应的参数即可运行该工作流。本次案例使用的结构是PDB ID4HBV。表中展示了计算得到的自由能变化ΔG的值该值为-7.7 Kcal/mol表示有利于驱动化学反应的进行。标准差为0.4 Kcal/mol相较于ΔG的绝对值来说较小所以可以认为计算出的ΔG-7.7 Kcal/mol是比较准确和可靠的。使用的约束类型是FloatBottom。采用哈密顿副本交换增强采样每个副本采样时间为1ns采样间隔为10ps因此每个副本采样数为100个。分析采用未分箱加权直方图分析方法UWHAM。表2: ATM结合自由能计算结果下图展示的是两条炼金术路径的扰动能量概率密度分布图。从图中可以看出在扰动能量比较小的时候采样很充分重叠也很好。但是当扰动逐渐增大时出现了不够重叠的概率密度分布这预示着1ns的采样时间可能不太够需要更加充分的采样以获得更良好的自由能估计。▊总结在药物研发中结合自由能计算对于评估药物分子与受体蛋白的相互作用至关重要。传统自由能微扰(FEP)方法虽能精确预测结合自由能但实施复杂且受限于多种技术因素。MaXFlow的Neo-BFE模块创新地引入了ATM方法通过坐标变换简化并加速结合自由能计算克服了FEP方法的部分局限性提供更高兼容性、稳定性和计算效率。Neo-BFE设计为“灰盒”工具既方便专家精细化调控也适合非专家用户快速准确地进行大规模计算。实验证明ATM方法在结合自由能预测上达到与FEP相当的精度展现出强大的应用潜力尤其在药物发现领域。[1] Sabanés Zariquiey F,Pérez A,Majewski M,et al.Validation of the Alchemical Transfer Method for the Estimation of Relative Binding Affinities of Molecular Series.[J].J Chem Inf Model,2023,63,8:2438-2444.