[发明专利]靶点选择方法及装置、计算机设备

说明书

本公开提供了一种靶点选择方法及装置、计算机设备。方法包括：获取mRNA靶点序列和NGS数据；基于所述mRNA靶点序列生成n1×N的序列特征矩阵；将所述序列矩阵和所述序列特征矩阵进行矩阵叠加，得到(n1+1)×N的一级结构矩阵；将所述一级结构矩阵降维处理成1×N的矩阵，并对所述1×N的矩阵进行对角化，得到N×N的二级结构矩阵；基于最优路径搜索算法从所述二级结构矩阵中选取一个或多个最优路径，并基于所述最优路径确定所有候选靶标序列。该方案在选择靶点时不仅考虑了NGS实验的数据，还考虑了各个位点的特征，使得选出的位点结合率高。并且，该方法无需对每个位点进行生物实验，耗时短，成本低，效果好。

技术领域

本公开涉及生物技术领域，特别涉及一种靶点选择方法及装置、计算机设备。

背景技术

基于核糖核酸(Ribonucleic Acid，RNA)的疗法，需要确定药物RNA分子和靶向RNA小分子。各种形式的RNA可用于选择性地作用于常规小分子或蛋白质无法靶向的蛋白质、转录物和基因。靶向RNA的疗法有小干扰RNA(siRNA)和反义寡核苷酸(AntisenseOligonucleotide，ASO)，虽然这两类药物特异性较高，但在屈指可数的靶器官范围外，仍有着不小的递送难题。

因此，选择RNA靶点成为当下RNA疗法的一个难题。相关技术中，通常采用生物实验的方法确定RNA序列上的各个位点能否作为靶点。但是，这种通过实验的方法耗时较长，并且成本高昂。

发明内容

本公开实施例提供了一种靶点选择方法及装置、计算机设备，能够通过计算机对RNA数据处理的方法确定靶点，减少耗时和成本。所述技术方案如下：

本公开至少一实施例提供了一种靶点选择方法，所述方法包括：

获取信使核糖核酸mRNA靶点序列和新一代测序NGS数据，所述mRNA靶点序列的长度为N，所述NGS数据为1×N的序列矩阵，所述序列矩阵的每个元素对应mRNA靶点序列的一个位点，所述序列矩阵的每个元素表示对应的位点在NGS实验中被N₂₀片段杂合的次数；

基于所述mRNA靶点序列生成n1×N的序列特征矩阵，所述序列特征矩阵的每一列元素对应mRNA靶点序列的一个位点，所述序列特征矩阵的每一列元素表示对应的位点的n1个特征；

将所述序列矩阵和所述序列特征矩阵进行矩阵叠加，得到(n1+1)×N的一级结构矩阵；

将所述一级结构矩阵降维处理成1×N的矩阵，并对所述1×N的矩阵进行对角化，得到N×N的二级结构矩阵，所述二级结构矩阵的每个元素表示对应二级结构的打开概率；

基于最优路径搜索算法从所述二级结构矩阵中选取一个或多个最优路径，并基于所述最优路径确定所有候选靶标序列。

可选地，所述特征包括如下至少一项：

回文结构、序列多态性、重复序列、翻译区起始点、编码区。

可选地，所述将所述一级结构矩阵降维处理成1×N的矩阵，并对所述1×N的矩阵进行对角化，得到N×N的二级结构矩阵，包括：

将所述一级结构矩阵的每一列元素相加，得到所述1×N的矩阵；

按照公式Pij＝Pi×Pj的方式计算所述二级结构矩阵中的各个元素，将计算结果填入所述二级结构矩阵中；其中，Pij为所述二级结构矩阵中的第i行第j个元素，Pi为所述1×N的矩阵中的第i个元素，Pj为所述1×N的矩阵中的第j个元素，所述i、j均为[1，N]中的正整数。

可选地，所述基于最优路径搜索算法从所述二级结构矩阵中选取一个或多个最优路径，并基于所述最优路径确定所有候选靶标序列，包括：