自由能微扰能够以化学精度计算化合物与靶标的结合自由能，主要用于先导化合物优化阶段。本期Uni-FEP靶点“探月”计划的FXR激动剂优化，选取了文献中具有代表性的化合物作为配体，使用Uni-FEP计算其与FXR的结合自由能，探索并验证了Uni-FEP在FXR体系中的适应性和准确性。
FXR靶点
法尼醇X受体（FXR）是核受体超家族和胆汁酸受体成员之一。FXR作为单体与类视黄醇X受体结合形成异源二聚体，并结合FXR应答元件，调控众多下游基因的转录。FXR主要在肝、肠、肾和肾上腺中大量表达[1,2]，通过复杂的信号通路网络，FXR与机体内相关激素或各代谢途径限速酶相互作用，参与调控机体胆汁酸代谢、脂质代谢、葡萄糖代谢等多种代谢途径，与糖尿病、非酒精性脂肪肝和肥胖症等多种代谢性疾病的发生发展有关[3,4]。
FXR作为治疗非酒精性脂肪性肝病的热门靶点，激动剂主要有甾体类和非甾体类两大类。甾体类FXR激动剂最初是从天然配体中研究发现，如鹅去氧胆酸。非甾体类FXR激动剂大多是人工合成的，结构类型多样，主要包括异恶唑类、苯并哌喃类、氮杂并吲哚类及苯并咪唑类等。本文以D3RCG2的FreeEnergySet2数据集的6个IC50跨度将近倍的苯并咪唑类FXR激动剂为例，探究了Hermite?Uni-FEP预测结合亲和力及指导先导化合物优化的能力。
Uni-FEP准确计算已发表
化合物与FXR的结合自由能
·计算体系
本文使用含有苯并咪唑骨架共晶分子的FXR蛋白晶体结构（PDBID：3OKI）快速验证文献报道[5]已有活性的6个苯并咪唑类化合物（图1）与FXR的结合能力。图1用于Uni-FEP计算的6个化合物
·计算过程
Hermite?Uni-FEP提供完整的、自动化的FEP计算流程和结果分析指导，整个流程包括：蛋白准备、配体准备、化合物对齐（支持刚性对齐和柔性对齐）、化合物对的自动Mapping/Pertubation、构建微扰图、结果分析。图2化合物对齐后与FXR靶点的结合
图3Uni-FEP构建的微扰图示意图。Exp.为实验ΔG值，Calc.为计算ΔG值，ΔG单位为kcal/mol，箭头上数字为化合物间的相似度
·计算结果
表1为Hermite?Uni-FEP计算ΔG与实验ΔG的比较表格。6个化合物的实验值和预测值具备良好的对应关系，除化合物FXR_13之外，其余化合物预测的结合自由能和实验值排序一致。单独分析每个化合物的计算结果和实验结果，所有化合物的预测结果与实验结果偏差均小于1kcal/mol，计算精度很高（备注：对于FEP的计算而言，预测值与实验值间的差值小于1.4kcal/mol，化合物活性相差10倍左右，属于精度可以接受的范围）。
表1：Hermite?Uni-FEP计算结果与实验结果比较
进一步分析预测ΔΔG与实验ΔΔG的相关性（图4），两者之间的相关性R2=0.89，RMSE=1.05kcal/mol，表明Uni-FEP在FXR的预测体系中，已经具备完美反映化合物实验亲和力的能力。
图4FEP计算结果图。横坐标为实验值，纵坐标为计算值。(a)Cycleclosure所得绝对结合自由能ΔG预测值和实验值的相关性；(b)相对结合自由能ΔΔG(pair)预测值和实验值的相关性
结论
本文以FXR靶点为案例，计算了FXR激动剂的结合自由能，结果表明预测的结合亲和力与实验的结合亲和力值具有较高的相关性，证实了Hermite?Uni-FEP能够通过计算一系列化合物与靶标的相对结合自由能，从而实现指导先导化合物优化的工作。
Hermite?Uni-FEP靶点“探月”计划更多靶点案例，请参考：
?
“骨架跃迁”不难算！Uni-FEP破解化合物优化计算难题
?
Uni-FEP靶点“探月”计划
聚焦热休克蛋白HSP90
?
Uni-FEP如何“破解”难成药靶标KRASG12D？
靶点“探月”计划
?
药物设计“手术刀”Uni-FEP解析Eg5靶点变构口袋分子
靶点“探月”计划
?
从虚拟到现实，Uni-FEP助力HNE靶点先导化合物优化
靶点“探月”计划
致谢
感谢中国药科大学的梁礼同学，在FXR体系中对Uni-FEP计算做出的重要贡献。
实习生招聘
Uni-FEP靶点“探月”计划招新啦~感兴趣的同学可以进入下方链接，投递简历：
