CN116092598A - 基于流形正则化非负矩阵分解的抗病毒药物筛选方法 - Google Patents

Abstract

本发明公开了一种基于流形正则化非负矩阵分解的抗病毒药物筛选方法，属于人工智能与生物学交叉领域。具体包括如下步骤：构建病毒‑药物关联的邻接矩阵，分别计算病毒高斯距离相似性和药物高斯距离相似性；输入病毒基因组序列计算病毒基因序列相似性，输入药物的化学结构计算药物化学结构相似性；使用快速核学习法计算整合病毒相似性、整合药物相似性；构造流形正则非负矩阵分解损失函数，使用迭代算法求得基矩阵与系数矩阵；计算关联对预测得分矩阵，筛选、排序给出最终结果。本发明能有效利用病毒‑药物关联数据中的阴性样本信息，高效、精确、快速地筛选出病毒有效治疗药物，解决生物实验方法的耗时长且成本高的问题。

Description

基于流形正则化非负矩阵分解的抗病毒药物筛选方法

技术领域

本发明涉及人工智能与生物学交叉领域，更具体地说，本发明涉及一种基于流形正则化非负矩阵分解的抗病毒药物筛选方法。

背景技术

寻找有前景的COVID-19抗病毒药物是当务之急，然而，传统的药物筛选包括目标发现、目标验证、先导化合物识别、先导优化、临床前药物开发、推进临床试验和临床试验等诸多阶段，昂贵且耗时的缺陷十分显著，因此亟需有效的机器学习建模计算方法实现“老药新用”的目的。已报道的方法存在许多不足之处，如使用随机采样方法生成负样本，引入了人为的不确切干扰信息，亦无法充分利用阴性样本信息，最终导致模型预测性能有限。

发明内容

为了克服现有技术的上述缺陷，本发明提供一种基于流形正则化的非负矩阵分解的病毒药物筛选方法（Manifold Regularized Non-negative Matrix Factorization forVirus Drug Association prediction，简称MRNMF-VDA），其步骤大体为，首先通过构建病毒-药物关联的邻接矩阵分别计算病毒高斯距离相似性和药物高斯距离相似性，然后根据病毒基因组序列计算病毒序列相似性，根据药物化学结构计算药物化学结构相似性，再使用快速核学习法计算整合病毒相似性、整合药物相似性；接下来使用病毒-药物关联邻接矩阵、整合病毒相似性和整合药物相似性构建目标函数并求解基矩阵和系数矩阵，得到病毒-药物关联对预测得分矩阵，最后排序给出最终预测结果。本方法能根据病毒-药物关联数据，精确高效地预测抗病毒相关药物。

为实现上述目的，本发明提供如下技术方案，一种基于流形正则化非负矩阵分解的病毒药物筛选方法，具体包括如下步骤：

步骤一，输入已知的病毒-药物关联对，构建邻接矩阵X，此矩阵的行数为病毒数量nv，列数为药物数量nd；

步骤二，分别计算药物间高斯距离相似性和病毒间高斯距离相似性：若药物d(i)与某个病毒之间存在关联，则对应位置记为1否则记为0，形成一个1×nv大小的0或1构成的向量，记之为药物d(i)的向量谱IP(d(i))，然后计算药物d(i)和d(j)之间的高斯距离相似性

上式中，参数γ_d用于控制核带宽，通过归一化新带宽参数γ’_d获得：

以类似的方式定义病毒v(i)和v(j)之间的高斯距离相似性，1×nd大小的0或1构成的向量，记之为病毒v(i)的向量谱IP(v(i))；

以上γ’_d和γ’_v都是常数；

步骤三，输入病毒基因组序列，使用既有较为成熟的多序列比方法计算病毒序列相似性，输入药物的化学结构得到药物MACCS指纹，采用谷本系数（Tanimoto Coefficient，也即Jaccard相似度，关联的交集除于关联的并集）计算药物化学结构相似性；

步骤四，使用快速核学习方法整合病毒基因组序列相似矩阵和病毒高斯距离相似矩阵，具体是通过求解下面的半正定规划式：

式中，第一项为范数项表示相似矩阵的整合误差大小，第二项为为正则化项，作用是避免过拟合，其中X为病毒-药物关联邻接矩阵，S_j ^v（j=1,2）表示病毒基因组序列相似矩阵、病毒高斯距离相似矩阵，μ^v为正则化参数，λ^v∈R^1×2为待求解的系数，通过λ^v得到病毒的整合相似矩阵，类似方法可获得药物结构相似性与药物高斯距离相似性集成参数λ^d∈R^1×2，然后计算药物整合相似矩阵，其中S_j ^d（j=1,2）表示药物结构相似矩阵、药物高斯距离相似矩阵；

步骤五，基于非负矩阵分解方法，引入拉普拉斯正则化和L_2，1范数约束项构建目标损失函数，具体为：

上式中X表示病毒药物关联，U和V分别是待求的基矩阵和系数矩阵，是非负矩阵分解的重构误差项，也称为最小二乘损失函数；和是刻画流形信息的拉普拉斯正则化项，用来保留局部结构信息；判别约束项，其中I_k是大小为k×k的单位矩阵，其目的是使V近似正交以获得判别信息；是矩阵V的L_2,1范数项，其中v_j.表示矩阵V的第j行。行稀疏正则化的目的是保留重要特征，并通过将V中的某些行向量缩小至0来删除不相关的特征；L₂正则化项确保矩阵光滑度防止过拟合，tr()表示求矩阵的迹，‖·‖_F表示矩阵的F范数。

步骤六，使用迭代方法求解目标函数得到基矩阵U和系数矩阵V，具体更新公式为：

，其中，

，其中R是一个对角矩阵，其第i个元素为：；类似地，

步骤七，计算评分矩阵F_m：

根据算出的病毒-药物关联对关系分数排序，给出最终预测结果。

本发明的技术效果和优点：

1、本发明在非负矩阵分解框架中引入流形正则化（拉普拉斯正则化）项，能有效捕获数据内在结构信息，只需要阳性样本和未标记样本而不依赖于负样本，大大降低了模型构建的难度，提升了预测性能。

2、L_2,1范数约束保证了系数矩阵能保留重要特征同时删除不相关的特征，故能减弱噪声数据的不利影响，进而获得较佳预测结果。

3、本发明合理融合了流形学习理论、非负矩阵分解和机器学习方法，该模型是半监督模型，能够高效、精确、快速地给出病毒-药物预测结果，且具有较好的可扩展性与健壮性。

附图说明

图1为本发明的总体流程图。

图2为本发明与几种已报道方法在同一数据集上五倍交叉验证的结果图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行完整地描述。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明的实施例中所使用的已知人类药物-病毒关联数据是从相关文献中收集的，使用文本挖掘技术对文献中经过实验验证的药物-病毒相互作用对进行整理后，获得455个已证实的人类病毒-药物相互作用，涉及34种病毒与219种药物（文献DOI: 10.1016/j.asoc.2021.107135）；药物化学结构从DrugBank数据库中下载，病毒基因组核苷酸序列从美国国家生物技术信息中心NCBI数据库获得，然后执行如图1所示的一种基于流形正则化非负矩阵分解的病毒药物筛选方法，具体包括如下步骤：

步骤一，输入已知的药物-病毒关联对，构建邻接矩阵X：

得到的邻接矩阵X元素为0或1，大小为34行×219列，i与j的取值范围满足1≤i≤34，1≤j≤219；

步骤二，分别计算病毒高斯距离相似性和药物高斯距离相似性：

若某一个病毒v(i)与某药物之间存在关联，则对应位置记为1否则记为0，形成一个1×219大小的0或1构成的行向量，记之为病毒v(i)的向量谱IP(v(i))，然后计算病毒v(i)和v(j)之间的高斯距离相似性：

上式中，参数γ_v用于控制核带宽，通过归一化新带宽参数γ’_v获得：

以类似的方式定义药物d(i)和d(j)之间的高斯距离相似性：

取；

其中nv表示病毒的数量，此例中为34，nd表示药物的数量，此例中为219，此步计算后得到大小为34×34的对称矩阵S₁ ^v（病毒高斯距离相似性）和大小为219×219的对称矩阵S₁ ^d（药物高斯距离相似性），且这两个矩阵全部元素大小都在0到1之间；

步骤三，输入病毒基因组序列，基于多序列比方法使用多序列比对工具MAFFT计算病毒的基因序列信息得到病毒序列相似矩阵S₁ ^v；输入SMILES编码表示的药物化学结构，然后用化学信息学软件RDKit或Open Babel获得药物的分子访问系统指纹（MACCS），再使用R包RxnSim计算Tanimoto相似度矩阵来测量两种药物之间的化学结构相似性S₁ ^d，具体计算方法是，对d(i)和d(j)两种药物，将此两种药物的MACCS片段二进制表示的字符串集分别记为D(i)和D(j)，d(i)和d(j)间的相似度S^d _ij值可以用下面公式计算：

步骤四，使用快速核学习方法整合病毒基因组序列相似矩阵和病毒高斯距离相似矩阵，具体通过求解下面的半正定规划式：

式中，第一项为范数项表示相似矩阵的整合误差大小，第二项为为正则化项，作用是避免过拟合，其中X为病毒-药物关联邻接矩阵，S_j ^v是病毒基因组序列相似矩阵、病毒高斯距离相似矩阵，μ^v为正则化参数，λ^v∈R^1×2为待求解的系数，可使用Matlab软件中的CVX工具箱求解此优化问题，得到病毒的整合相似矩阵，类似地，使用快速核学习方法可获得药物结构相似性与药物高斯距离相似性集成参数λ^d∈R^1×2，然后计算药物整合相似矩阵定义为；

步骤五，基于非负矩阵分解方法，引入拉普拉斯正则化和L₂范数约束项构建目标损失函数，具体为：

上式中X表示病毒药物关联，U和V分别是待求的基矩阵和系数矩阵，是非负矩阵分解的重构误差项，也称为最小二乘损失函数；和是刻画流形信息的拉普拉斯正则化项，用来保留局部结构信息；判别约束项，其中I_k是大小为k×k的单位矩阵，其目的是使V近似正交以获得判别信息；是矩阵V的L_2,1范数项，其中v_j.表示矩阵V的第j行。行稀疏正则化的目的是保留重要特征，并通过将V中的某些行向量缩小至0来删除不相关的特征；L₂正则化项确保矩阵光滑度防止过拟合，tr()表示求矩阵的迹，‖·‖_F表示矩阵的F范数。

步骤六，根据拉格朗日乘数法求偏导再使用Karush-Kuhn-Tucker（KKT）条件将矩阵U和V反解出来，便可得到其更新公式：

，其中，

其中R是一个对角矩阵，其第i个元素为：；类似地，有；L_v计算方法为病毒整合相似矩阵S^v按列求和后再对角化，L_d计算方法为药物整合相似矩阵S^d按列求和后再对角化；

在上述的实施例中，经初步优化后，选取参数值λ₁=λ₂=1，μ=0.4，γ₁=γ₂=1；

使用Matlab编程实现上述算法时，将矩阵U初始化为34行×30列的随机矩阵，V初始化为30行×219列的随机矩阵（对角矩阵R大小同矩阵V的行数，此处为30行×30列），U和V的全部元素都在(0,1)区间范围，设定迭代循环次数为50次或者满足退出迭代循环，运行结束后得到矩阵U、V。

步骤七，计算评分矩阵F_m：

根据算出的药物-病毒关联对关系分数排序给出最终预测结果。

本发明的有效性验证：

如图1所示的一种基于流形正则化非负矩阵分解的病毒药物筛选方法，采用五重交叉验证进行预测性能评估，具体实施方式为：先将所有已知的药物-病毒关联随机平均分成5组，再将5组中的每一组分别设为测试样本，其他组作为训练样本。使用训练样本作为本方法的输入得到预测结果，最后将该组中每个测试样本的预测分数与候选的分数进行比较。为了减少在获得测试样品的过程中随机划分对结果造成的影响，进行了100次五折交叉验证。

使用Matlab计算后获得了如下数据，如图2所示为本方法MNMFVDA与现已报道的的几种病毒-药物筛选模型之间的AUROC（ROC曲线下面积）值比较。本方法在5折交叉验证中取得了0.7686±0.0044的AUROC值，表现出了比几种经典模型更加出色的预测性能。

另外一方面，对具体某种病毒，如新型冠状病毒（SARS-CoV-2）使用本方法做预测，筛选评分矩阵F_m中SARS-CoV-2对应的行即取得相关药物的预测得分，将其降序排列后中前20个药物有17个能够得到已报道文献的支持。

下表展示了预测结果前20个药物名称和支持的文献PMID号。

最后：以上所述仅为本发明的优选实施例之一，并不用于限制本发明，凡在本发明的精神和原则之内，所作的细微修改、等同替换、简单改进等，均应包含在本发明的保护范围之内。

Claims (2)

1.基于流形正则化非负矩阵分解的病毒药物筛选方法，其特征在于，具体包括如下步骤：

步骤一，输入已知的病毒-药物关联对，构建邻接矩阵；

步骤二，根据邻接矩阵分别计算药物间高斯距离相似性和病毒间高斯距离相似性；

步骤三，输入病毒基因组序列，基于既有的多序列比对法计算病毒序列相似性；输入药物的化学结构的SMILES编码，使用谷本系数计算药物化学结构相似性矩阵；

步骤四，使用快速核学习方法整合病毒基因组序列相似矩阵和病毒高斯距离相似矩阵、药物结构相似矩阵和药物高斯距离相似矩阵；

步骤五，基于非负矩阵分解整合拉普拉斯正则化项和L_2,1范数约束项构建目标损失函数；

步骤六，求解目标函数，得到基矩阵U和系数矩阵V；

步骤七，计算评分矩阵，根据算出的病毒-药物关联对关系分数排序，给出最终预测结果。

2.根据权利要求1所述的基于流形正则化非负矩阵分解的病毒药物筛选方法，其特征在于，上述步骤五中的目标函数以非负矩阵分解为主体且加入了拉普拉斯正则化和L_2,1范数约束项，具体为：

式中X表示病毒-药物关联邻接矩阵，U和V分别是待求的基矩阵和系数矩阵，是非负矩阵分解的重构误差项，也称为最小二乘损失函数；和是刻画流形信息的拉普拉斯正则化项，用来保留局部结构信息；判别约束项，其中I_k是大小为k×k的单位矩阵，其目的是使V近似正交以获得判别信息；是矩阵V的L_2,1范数项，其中v_j.表示矩阵V的第j行；行稀疏正则化的目的是保留重要特征，并通过将V中的某些行向量缩小至0来删除不相关的特征；L₂正则化项确保矩阵光滑度防止过拟合，tr( )表示求矩阵的迹，‖·‖_F表示矩阵的F范数。

Images