北京哪里有专门做网站的地方,前端做的好的网站,大数据工程技术,12380网站建设的意见建议文 | 李梅、黄楠源 | AI科技评论2018 年#xff0c;Deepmind 首次发布基于深度神经网络的蛋白质结构预测数据库 AlphaFold#xff0c;在蛋白质预测中实现了最先进的性能#xff1b;去年#xff0c;AlphaFold 2 获得了 98.5% 的蛋白质预测率#xff1b;前段时间#xff0c…文 | 李梅、黄楠源 | AI科技评论2018 年Deepmind 首次发布基于深度神经网络的蛋白质结构预测数据库 AlphaFold在蛋白质预测中实现了最先进的性能去年AlphaFold 2 获得了 98.5% 的蛋白质预测率前段时间Deepmind 又重磅发布了数据集更新称目前的 AlphaFold 已经预测了几乎所有已知的蛋白质。如何有效识别药物作用机制在今天仍然是一个巨大挑战计算对接的方法已被广泛用于预测药物结合靶点。有了大规模蛋白质结构预测技术药物发现将变得更容易。所以自 AlphaFold 问世以来称其将引发一场结构生物学的革命、彻底改变药物发现的声音就不绝于耳。本质上AlphaFold 是一个工具我们目前真的能利用好这个工具吗近日来自 MIT 的研究团队给出了否定的回答。他们对使用 AlphaFold2 的分子对接模拟的模型性能进行了评估发现模型在识别真正的蛋白质-配体相互作用方面的预测能力较弱并证明需要使用基于机器学习的方法进行建模来提高模型性能以更好地利用AlphaFold2 进行药物发现。该论文“Benchmarking AlphaFold-enabled molecular docking predictions for antibiotic discovery”发表在了Molecular Systems Biology 期刊上。论文地址https://www.embopress.org/doi/epdf/10.15252/msb.202211081使用 AlphaFold 2 预测的分子对接所谓化合物的对接计算是将候选化合物列表中的每一个对接到目标蛋白质中生成最有可能结合的化合物的粗略排序。这个过程可以在化合物集合上完成获得庞大的虚拟库这种虚拟筛选已经成为计算化学领域的长期目标。筛选 218 种大肠杆菌活性化合物研究团队首先筛选了一组化合物包含大约 39128 种其中包括已知药物已知抗生素、活性天然产物和一系列其他不同结构并在针对大肠杆菌的筛选中发现了 218 种化合物培养物。仅仅有 218 个阳性这个结果是令人惊讶的但考虑到抗菌药物发现工作的难度这个数字也算比较难得了。在 218 种活性化合物中有大约 80% 是已知抗生素类别中的成员剩下的部分则是已知细胞毒性化合物和一些新的通配类型的混合。这为实验的进行提供了一个很好的背景因为在大多数情况下我们可以预测从反向对接筛选中获得什么结果。将活性化合物与预测的蛋白质结构对接接着团队研究了这些活性化合物的潜在结合靶标。多年来在大肠杆菌中进行的大量基因组敲除扫描的共识评分已经得出了 296 项基本蛋白质所以可以合理推断所有真正抑制生长的靶标蛋白质都可能在这些列表当中。作者将 218 种活性化合物中与 AlphaFold 2 预测出的 296 种基本大肠杆菌蛋白质结构进行对接并通过几种不同的计算方法对 218 种化合物与 296 种蛋白质的组合进行计算预测了 64000 多个蛋白质-配体对的结合位姿与结合亲和力预测。▲图注在 AlphaFold 结构上进行分子对接从计算量来看这是一个可靠度相当高的测试尤其是考虑到内部控制的数量具有已知靶标的化合物和在这些靶标内具有已知结合构象的化合物这项计算是很有价值的。作为对照研究团队还从一组化合物中随机选择了 100 种对细菌生长完全没有抑制作用的化合物进行相同的计算从而获得对 29600 个蛋白质-配体对的结合位姿与亲和力预测。基于 AlphaFold 2 预测结构的模型性能很弱虽然这项工作预测了包括活性和非活性化合物的化合物与蛋白质混杂性但问题是这些预测中有多少是假阳性将模型预测与已知的抗生素结合目标进行比较为了评估所用模型方法的性能作者将模型预测与常用抗生素类别的已知相互作用进行比较。作者搜集了先前文献中的抗生素-蛋白质靶对组成一个包含 142 种抗生素-蛋白质相互作用的数据集。结果发现他们的模型仅仅正确预测了 3 种具有强结合性即结合亲和力阈值为 -7 kcal/mol 的相互作用以及 43 种具有一般结合性即结合亲和力阈值为 -5 kcal/mol 的相互作用。所以模型预测的真阳性率分别为 2.1% 和 30.3%。这种比较表明基于 AlphaFold 2 预测结构的建模平台性能很弱。测量 12 种基本蛋白质的酶抑制作者接着选取了 12 种基本蛋白质它们可以用于酶促测定通过测量 218 种活性化合物对这些蛋白质的酶抑制作者对模型预测的子集进行进一步的评估。▲图注所有 218 种活性化合物的平均相对活性12 种蛋白质都经过抑制实验测试。结合相互作用命中是蛋白质-配体相互作用红点所有其他交互都被指定为非命中灰点。结果表明所有测试中的基本蛋白质都被至少四种不同的化合物所抑制涵盖从强到弱的一系列结合亲和力阈值基于 AlphaFold 2 的模型所预测的对接显示出了广泛的混杂性。模型的基准测试最后作者对建模平台的性能进行了统计基准测试。基于酶抑制测量数据作者将实验观察到的结合相互作用命中与他们预测的相互作用进行比较结合亲和力阈值设定为小于 -5 kcal/mol 和 -7 kcal/mol。结果发现更强结合性的亲和力阈值会导致更少的结合相互作用预测此时真阳性率更低准确性更高。模型性能根据所设定的结合亲和力阈值从弱到中等变化。作者又在独立于结合亲和力阈值的条件下使用接受者操作特征曲线 (ROC) 和精确召回曲线 (PR) 再次进行评估也都表明模型性能较弱。总之模型的预测结果中不仅存在大量假阳性即非活性化合物被预测为与关键细菌蛋白质的活性位点结合而且还有大量的假阴性即已知存在相互作用但没有被发现。只有达到最严格的结合亲和力阈值时模型才会比随机预测表现得略好一些。AlphaFold 本身没错 用好机器学习方法是关键接下来的问题是模型的弱性能是由 AlphaFold2 所提供的蛋白质结构质量导致的吗问题出自对接方法而非蛋白质结构质量为了验证这个问题作者将 218 种活性化合物与八种实验确定的蛋白质结构中的每一种对接进行了重复的对接模拟并同样对模型性能进行了基准测试。结果是 auROC 值在数量上与先前相似范围从 0.25 ( glmU ) 到 0.69 ( gyrAB )平均值为 0.46。auPRC 值也发现了类似的结果范围从 0.03 ( ligA ) 到 0.56 ( gyrAB )平均值为 0.22。这些发现表明使用 AlphaFold2 预测结构的分子对接与使用实验确定的结构是类似的。这也与之前对 AlphaFold 对实验确定的蛋白质结构的保真度评估一致由此可以得出模型的性能弱是因为对接方法的原因而不是蛋白质结构的质量差。使用机器学习方法可改进模型性能基于分子对接的弱性能问题研究团队探索了可以提高性能的方法。研究中使用了四种不同的基于机器学习的评分函数分别是 RF-Score 、RF-Score-VS、PLEC score 和 NNScore以对模型性能进行基准测试和改进。相比于 RF-Score 和 RF-Score-VS - RF-Score 的虚拟筛选适应性--利用随机森林或决策树的组合来预测蛋白质与配体的结合亲和力PLEC score 采用了蛋白质-配体对之间的扩展连接指纹NNScore 是基于神经网络的集合。作者在研究中采用了评分函数使用 PDBbind v2016 或有用的诱饵目录对增强DUD-E数据库进行训练以重新评估 AutoDock Vina 预测的对接姿势。此外研究使用 DOCK6.9 和应用于 AutoDock Vina 姿势的每个基于机器学习的评分函数还预测了每种抗菌化合物与 12 种经验测试必需蛋白中每一种之间的结合亲和力并对每种方法的性能进行基准测试。测试结果发现平均 auROC 值在 0.46 和 0.63 之间下图 A。其中与 DOCK6.9 对接并使用 PLEC score 对 AutoDock Vina 姿势进行重新评分平均导致 auROC 值低于单独使用 AutoDock Vina 的结果DOCK6.9 的 auROC 值为为 0.46范围为 0.25 至 0.61和 0.47范围 PLEC score 为 0.28 至 0.63下图 A相比之下使用 RF-Score、RF-Score-VS 或 NNScore 对 AutoDock Vina 姿势进行重新评分可提高模型性能平均 auROC 值分别为 0.62范围为 0.53 至 0.69、0.63范围为 0.46 至 0.75和 0.58范围为 0.41 到 0.69。研究结果也与 auPRC 相似当使用 RF-Score 重新评分时其平均值高达 0.24。这些模型性能评估表明某些基于机器学习的评分函数提高了预测准确性。▲图注使用机器学习对模型性能进行基准测试和改进A. 在不同的分子对接程序和不同的基于机器学习的姿势评分函数。白点表示平均值灰色条表第25-75个百分位值的范围灰色箱线图须线表示不被视为异常值的值范围0.5 处的水平线表示随机预测生成的基准。B. 通过在 AutoDock Vina 应用基于机器学习的重新评分函数建模的蛋白质-配体对的排序结合亲和力 。曲线根据 (A) 中使用的重新评分函数着色阴影区域表示 7 的结合亲和力阈值。C-E. 预测准确性、预测阳性数蛋白质-配体相互作用和真阳性率/假阳性率对所用模型数量的依赖性群体智慧方法可提高预测准确性由于某些基于机器学习的评分函数会增加 auROC 和 auPRC研究还探讨了在严格限制结合亲和力阈值的情况下结合“群体智慧”方法使用重新评分模型是否可以提高预测准确性和真阳性率。作者将预测的蛋白质-配体相互作用定义为满足所有模型的结合亲和力阈值并将 AutoDock Vina 预测与上述四种基于机器学习的评分函数的预测相结合研究通过使用这种共识方法发现预测准确性可随着使用的模型数量而提高上图 C这同预测的蛋白质-配体相互作用数量的相应减少预期一致上图 D。与此同时真阳性率与假阳性率的比率则是随使用模型数量的增加而增加在预期之外上图 E。可以看到该结果同使用某些基于机器学习的评分函数提高预测能力的发现一致这也进一步表明了将分子对接与基于机器学习的模型结合起来可以让人们更好地利用 AlphaFold2 预测的蛋白质结构进行药物筛选。所以一些机器学习方法确实可以提高了预测的准确性。不过这只是部分的成功当前研究所用的数据集中有很多已经确定的蛋白质和化合物的实验事实如果涉及那些较少被关注的领域这些方法是否仍然奏效就不可知了。虽然 AlphaFold 为我们提供了大量且合理的蛋白质结构但我们实现它的价值的能力还非常有限。所以至少在目前看来“AlphaFold 将彻底改变药物发现”的说法还尚待证实成功还在未来。后台回复关键词【入群】加入卖萌屋NLP、CV、搜广推与求职讨论群[1] https://www.science.org/content/blog-post/not-alphafold-s-fault
