「超级细菌」魔咒或将打破,MIT 利用深度学习发现新型抗生素

2024-01-11 05:52

本文主要是介绍「超级细菌」魔咒或将打破,MIT 利用深度学习发现新型抗生素,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

作者:加零

编辑:李宝珠、三羊

MIT 利用图神经网络 Chemprop 识别潜在抗生素,特异性杀死鲍曼不动杆菌。

自然界中充满了各种各样的微生物,例如结核杆菌(导致肺结核)、霍乱弧菌(导致霍乱)等严重危害了人们的健康。而在人类历史上,感染这些致病细菌,除了依赖人体自身的免疫系统,几乎无解。直到 1928 年,青霉素的发现,让人类第一次拥有了战胜致病细菌的有力武器。

然而,抗生素的广泛应用也带来了一个巨大的危机 —— 抗生素耐药性 (AMR),据世界卫生组织 (WHO) 统计, 2019 年全球约有 120 万人死于抗生素耐药性 (AMR) 所加剧的细菌感染,这已经高于艾滋病导致的死亡人数。滥用抗生素催生某些「超级细菌」的出现,成为了 21 世纪临床重要致病因。要解决这一难题,迫切需要开发全新的抗生素。

针对特异细菌,深度学习模型可以提高对化合物药效和安全性的预测准确性,有效减少实验室实验和临床试验的时间和资源消耗,对有效且安全的抗生素发现至关重要。

为此,来自 MIT 的研究者们开发了一种深度学习方法来发现抗生素,利用图神经网络 Chemprop 从大型化学库中识别潜在的抗生素,并发现了一类新型抗生素。 他们能够特异性杀死重要的耐药菌 —— 鲍曼不动杆菌 (Acinetobacter baumannii),相关论文已发表于 Nature。

在这里插入图片描述

论文已发表于 Nature

论文地址:
https://www.nature.com/articles/s41586-023-06887-8
关注公众号,回复「发现抗生素」下载完整论文

实验方法:利用图神经网络的深度学习

数据集:多种化合物

初步筛选: 该研究初步筛选了 39,312 种化合物,对其抗生素活性和人类细胞毒性特征进行分析。
扩大预测: 为进一步扩大预测范围,这些模型对 12,076,365 种化合物进行测试,其中 11,277,225 种来自 Mcule 数据库,799,140 种来自 Broad Institute 数据库。

在这里插入图片描述

39,312 个化合物的分子量分布

算法训练:利用图神经网络进行训练

使用图神经网络 Chemprop 对筛选出的 39,312 种化合物进行训练,预测其抗生素活性和人类细胞毒性,训练过程如下。

在这里插入图片描述

分子表示: 使用 RDKit 从每个化合物的 SMILES(Simplified Molecular Input Line Entry Specification,简化分子输入线性入口系统)字符串生成基于图的分子表示。

特征向量生成: 为每个原子和键生成特征向量,包括原子特征(如原子序数、键数、正电荷等)和键特征(如键类型、共轭、环成员资格等)。

信息传递: 实现基于键的信息传递卷积神经网络,更新键的消息,通过神经网络层并应用非线性激活函数。

模型输出: 模型在固定的信息传递步骤后,对整个分子的消息进行汇总,通过前馈神经网络预测化合物的活性,如抗生素活性、细胞毒性或改变质子动力势的活性。

优化措施: 包括添加额外的分子级特征、使用超参数优化选择最佳性能参数,以及通过集成学习增强模型的鲁棒性。

模型架构:保证化合物的有效性和安全性

抗生素活性模型

研究人员在化合物浓度 50 μM 的培养基中预测其对金黄色葡萄球菌 (S. aureus) 的抗生素活性,以 80% 标准化生长抑制截止值 (normalized growth inhibition cut-off) 区分为活性或无活性。共 10 个图神经网络模型在同一 80% - 20% 的训练数据集上进行了训练、验证和测试。

结果显示,具有 RDKit 特征的 Chemprop 模型表现出更优越的预测能力,在 39,312 种化合物中识别出 512 种活性化合物。

在这里插入图片描述

预测抗生素活性的深度学习模型比较

人类细胞毒性模型

研究人员利用 39,312 种化合物对人类肝癌细胞 (HepG2)、人类初级骨骼肌细胞 (HSkMCs) 和人类肺成纤维细胞 (IMR-90) 进行毒性筛选。在用每种化合物 10 μM 浓度处理 2-3 天后,评估细胞活力,以 90% 细胞活力截止值 (cell viability cut-off) 将化合物活性分类。

同样对 10 个 Chemprop 模型集合进行了训练、验证和测试,对比结果如下图:

在这里插入图片描述

预测人类细胞毒性的深度学习模型比较

结果显示,分别有 3,341 (8.5%)、1,490 (3.8%) 和 3,447 (8.8%) 的化合物对 HepG2 细胞、HSkMCs 和 IMR-90 细胞具有毒性。在上步筛选出的 512 种活性抗菌化合物中,306 种对这三种细胞类型均无毒性。

总结来看,尽管与抗生素活性模型相比存在一定局限性,这一模型平衡了药物的有效性和对人体的无害性,展示了在药物发现中使用先进计算方法的潜力。

实验结果:抗生素的筛选识别

精炼和大范围应用模型:全化学空间的筛选和可视化

在这一研究阶段中,重点在于精炼和应用模型,以便在庞大的化学空间中识别潜在的抗生素化合物,并评估它们的细胞毒性。研究人员重新训练了 20 个 Chemprop 模型用于预测抗生素活性和 HepG2、HSkMC、IMR-90 细胞的细胞毒性,改进后的模型被应用于 12,076,365 种化合物预测上。

化合物筛选

抗生素活性筛选:从 Mcule 数据库中筛选出抗生素预测得分超过 0.4 的 3,004 种化合物;从 Broad Institute 数据库中筛选出得分超过 0.2 的 7,306 种化合物。

人类细胞毒性筛选:

保留细胞毒性预测得分低于 0.2 的化合物,最终筛选出 3,646 种化合物(1,210 种来自Mcule 数据库,2,436 种来自 Broad Institute 数据库)——占所有评估化合物的 0.03%。

在这里插入图片描述

化合物筛选

a:计算机模拟过滤程序b-e:抗生素活性和 HepG2、HSkMC、IMR-90 细胞的细胞毒性预测

化学空间的可视化

以摩根指纹 (Morgan fingerprints) 为分子表示,采用 t 分布随机相邻嵌入 (t-SNE, t-distributed stochastic neighbour embedding) 方法对化学空间进行可视化。

如下图所示,t-SNE 的可视化揭示了 hits(通过筛选的化合物)和 non-hits(抗生素预测得分低的化合物)之间的明显差异。

在这里插入图片描述

对化合物的 t-SNE 分析

进一步筛选:识别两个有效化合物

在 3,646 个化合物中筛选出 2 个对金黄色葡萄球菌表现出高活性且对人类细胞具有良好选择性的化合物(1 号和 2 号),这两种化合物在各种测试条件下的表现、特别是在含血清的培养基中的生长抑制能力很优越,值得进一步研究。

这些化合物的研究表明,深度学习模型预测的结构类别可以有效地指导实验筛选,从而发现新的抗生素候选物。

筛选过程

化合物筛选:从最初的 3,646 个 hits 化合物中去除了含有可能反应活性、致突变或药代动力学不利的 PAINS 和 Brenk 警报的化合物,筛选至 2,209 个。

结构筛选:进一步筛选与训练集中化合物结构不同的化合物,使用最大 Tanimoto 相似性得分 ≤0.5 作为初步切断点,并排除含 β-内酰胺环或喹诺酮双环核心的化合物,最终得到 1,261 个化合物。

有效化合物的识别

生长抑制测试: 在 9 个与理性组 G1-G5 相关的 hits 化合物中,发现 4 个化合物 (44%) 对金黄色葡萄球菌 (S. aureus) 有活性,最小抑制浓度 (MIC) ≤32 μg/ml。

结构类别和有效性:这些有效化合物与理性组 G1、G2 和 G5 相关,其中 G2 组的两个化合物(1号和2号)被证实为活性。这两个化合物同时符合 Lipinski 的规则和 Ghose 标准,表明它们具有良好的口服生物利用度和药物样特性,值得进一步研究。

在这里插入图片描述

筛选出的两个化合物

进一步研究:两个化合物的特性

通过深入的机制研究和体内外实验,化合物 1 和 2 显示出作为新型抗生素候选物的潜力。它们不仅针对多重耐药菌株有效,而且具有较低的耐药性发展趋势和良好的安全性。

这些发现表明,这两种化合物可作为抗生素药物发展的有希望的化学系列。

作用机制和耐药性

共同结构: 化合物 1 和 2 共享 N-[2-(2-chlorophenoxy)ethyl]aniline core 这一结构,预测表明这可能是抗生素活性的基本结构。

生长抑制实验: 在对金黄色葡萄球菌和枯草杆菌的时间杀菌实验中,这两种化合物显示出了抗生素活性,与范科霉素类似但杀菌力较低。

耐药性研究: 在对抗生素耐药的金黄色葡萄球菌株进行的实验中,这两种化合物的最小抑制浓度 (MIC) 只略有增加,表明它们可能具有与常见抗生素不同的作用机制。

耐药性的发展: 经过 30 天的连续培养,这两种化合物的 MIC 几乎没有变化,显示出较低的耐药性发展趋势。

抗多重耐药菌效果

广泛抗性: 这两种化合物对包括万古霉素耐药性肠球菌在内的 40 种不同细菌种类表现出活性,其中 MIC 的中位数分别为 4 和 3 μg/ml。

对静止期细菌的有效性: 这两种化合物对枯草杆菌的静止期细胞也显示出活性。

毒理学、化学特性及功效

安全性研究: 这两种化合物在体外实验中表现出良好的安全性,包括非溶血性、不结合金属离子、无基因毒性、化学稳定性以及在小鼠体内的安全性。

体内功效实验: 在小鼠金黄色葡萄球菌皮肤感染和大腿感染模型中,化合物 1 表现出显著的抗菌活性。

在这里插入图片描述

化合物的体内疗效

深度学习:应对抗生素耐药性的利器

研究人员多年来一直在探索应对抗生素耐药性的有效、易推广的方法,在这一过程中,深度学习的出现为研究人员提供了新的解题思路,深度学习在应对抗生素耐药性方面的应用价值在于:

超越传统抗生素发现方法: 传统的抗生素发现方法往往依赖于已知的活性结构,这限制了新药发现的范围。深度学习方法可以识别与传统抗生素结构不同的新化合物,可能对当前耐药的菌株有效。

个性化和精准医疗: 深度学习可以用于分析特定病原体的遗传和表型特征,从而有助于开发针对特定病原体或感染类型的个性化抗生素。

道阻且长,行之将至。深度学习在药物研发中的应用仍处于相对早期阶段,可能会面临数据质量和可解释性等挑战,但作为人类对抗细菌的重要防线,相关研究意义重大,相信势必会在技术迭代的加持下持续向前。

参考文献:
https://www.nature.com/articles/s41586-023-06887-8

这篇关于「超级细菌」魔咒或将打破,MIT 利用深度学习发现新型抗生素的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/593360

相关文章

HarmonyOS学习(七)——UI(五)常用布局总结

自适应布局 1.1、线性布局(LinearLayout) 通过线性容器Row和Column实现线性布局。Column容器内的子组件按照垂直方向排列,Row组件中的子组件按照水平方向排列。 属性说明space通过space参数设置主轴上子组件的间距,达到各子组件在排列上的等间距效果alignItems设置子组件在交叉轴上的对齐方式,且在各类尺寸屏幕上表现一致,其中交叉轴为垂直时,取值为Vert

Ilya-AI分享的他在OpenAI学习到的15个提示工程技巧

Ilya(不是本人,claude AI)在社交媒体上分享了他在OpenAI学习到的15个Prompt撰写技巧。 以下是详细的内容: 提示精确化:在编写提示时,力求表达清晰准确。清楚地阐述任务需求和概念定义至关重要。例:不用"分析文本",而用"判断这段话的情感倾向:积极、消极还是中性"。 快速迭代:善于快速连续调整提示。熟练的提示工程师能够灵活地进行多轮优化。例:从"总结文章"到"用

【前端学习】AntV G6-08 深入图形与图形分组、自定义节点、节点动画(下)

【课程链接】 AntV G6:深入图形与图形分组、自定义节点、节点动画(下)_哔哩哔哩_bilibili 本章十吾老师讲解了一个复杂的自定义节点中,应该怎样去计算和绘制图形,如何给一个图形制作不间断的动画,以及在鼠标事件之后产生动画。(有点难,需要好好理解) <!DOCTYPE html><html><head><meta charset="UTF-8"><title>06

学习hash总结

2014/1/29/   最近刚开始学hash,名字很陌生,但是hash的思想却很熟悉,以前早就做过此类的题,但是不知道这就是hash思想而已,说白了hash就是一个映射,往往灵活利用数组的下标来实现算法,hash的作用:1、判重;2、统计次数;

Andrej Karpathy最新采访:认知核心模型10亿参数就够了,AI会打破教育不公的僵局

夕小瑶科技说 原创  作者 | 海野 AI圈子的红人,AI大神Andrej Karpathy,曾是OpenAI联合创始人之一,特斯拉AI总监。上一次的动态是官宣创办一家名为 Eureka Labs 的人工智能+教育公司 ,宣布将长期致力于AI原生教育。 近日,Andrej Karpathy接受了No Priors(投资博客)的采访,与硅谷知名投资人 Sara Guo 和 Elad G

零基础学习Redis(10) -- zset类型命令使用

zset是有序集合,内部除了存储元素外,还会存储一个score,存储在zset中的元素会按照score的大小升序排列,不同元素的score可以重复,score相同的元素会按照元素的字典序排列。 1. zset常用命令 1.1 zadd  zadd key [NX | XX] [GT | LT]   [CH] [INCR] score member [score member ...]

【机器学习】高斯过程的基本概念和应用领域以及在python中的实例

引言 高斯过程(Gaussian Process,简称GP)是一种概率模型,用于描述一组随机变量的联合概率分布,其中任何一个有限维度的子集都具有高斯分布 文章目录 引言一、高斯过程1.1 基本定义1.1.1 随机过程1.1.2 高斯分布 1.2 高斯过程的特性1.2.1 联合高斯性1.2.2 均值函数1.2.3 协方差函数(或核函数) 1.3 核函数1.4 高斯过程回归(Gauss

【学习笔记】 陈强-机器学习-Python-Ch15 人工神经网络(1)sklearn

系列文章目录 监督学习:参数方法 【学习笔记】 陈强-机器学习-Python-Ch4 线性回归 【学习笔记】 陈强-机器学习-Python-Ch5 逻辑回归 【课后题练习】 陈强-机器学习-Python-Ch5 逻辑回归(SAheart.csv) 【学习笔记】 陈强-机器学习-Python-Ch6 多项逻辑回归 【学习笔记 及 课后题练习】 陈强-机器学习-Python-Ch7 判别分析 【学

系统架构师考试学习笔记第三篇——架构设计高级知识(20)通信系统架构设计理论与实践

本章知识考点:         第20课时主要学习通信系统架构设计的理论和工作中的实践。根据新版考试大纲,本课时知识点会涉及案例分析题(25分),而在历年考试中,案例题对该部分内容的考查并不多,虽在综合知识选择题目中经常考查,但分值也不高。本课时内容侧重于对知识点的记忆和理解,按照以往的出题规律,通信系统架构设计基础知识点多来源于教材内的基础网络设备、网络架构和教材外最新时事热点技术。本课时知识

线性代数|机器学习-P36在图中找聚类

文章目录 1. 常见图结构2. 谱聚类 感觉后面几节课的内容跨越太大,需要补充太多的知识点,教授讲得内容跨越较大,一般一节课的内容是书本上的一章节内容,所以看视频比较吃力,需要先预习课本内容后才能够很好的理解教授讲解的知识点。 1. 常见图结构 假设我们有如下图结构: Adjacency Matrix:行和列表示的是节点的位置,A[i,j]表示的第 i 个节点和第 j 个