一种分子对接药物筛选方法和装置制造方法及图纸

本发明专利技术涉及一种分子对接药物筛选方法和装置，属于量子计算与药物设计技术领域，解决现有方法对指数大的分子及化合物空间，复杂度高且在训练过程中不稳定的问题。方法包括：将药物分子按其结合基团的特征抽象成结合部位节点和结合原件节点的集合；根据结合基团之间的亲和度打分函数值构建亲和度矩阵并列出约束条件和优化目标；对需要优化设计的相关变量进行简化；根据亲和度矩阵基于横向独热模型构建量子比特分配与算子映射；根据Mixer和Phaser算子分层构建QAOA参数化线路；设置迭代次数优化步长，通过损失函数的梯度下降训练QAOA参数化线路；和多次测量训练后的QAOA参数化线路并获取目标函数最大的比特串。引入辅助优化变量，提高优化效率，减少量子比特。减少量子比特。减少量子比特。

【技术实现步骤摘要】
一种分子对接药物筛选方法和装置


[0001]本专利技术涉及量子计算与药物设计
，尤其涉及一种分子对接药物筛选方法和装置。

技术介绍

[0002]新型药物的研发过程主要由目标发现，分子设计，临床前研究和临床试验这几个阶段组成。食品药品监督管理局批准的大部分新药都是小分子药物，其结构和功能的多样性使其与生物结合位点的匹配成为可能。长期以来，药物研发主要通过人工反复实验探索、对比为主，不仅需要时间、劳力和财力，还存在较大的误差与不确定性。为了解决上述问题，物理和化学、计算机和信息科学、材料科学等多种学科越来越多地参与到创新药物的研发中，具有重大潜力的新技术、新方法不断涌现。
[0003]2014年，Edward Farhi等人提出了量子近似优化算法(QAOA：Quantum Approximate Optimization Algorithm)，并将其运用在了有界约束问题中。量子近似优化算法是一种广泛研究的方法，用于解决NISQ设备上的组合优化问题。量子近似优化算法通过时间演化在哈密顿量和量子电路之间建立联系，并构建层级的电路来实现强大的算法功能。这是一种启发式方法，可以被认为是绝热量子计算的时间离散化。与量子计算的许多近期方法一样，QAOA是一种混合经典量子算法，它将量子电路和经典优化问题相结合。量子QAOA的应用广泛而深远，其算法的性能也是量子计算研究界非常感兴趣的。
[0004]综合上述背景，在医药领域，分子对接作为计算基辅助药物设计的重要应用，是通过研究小分子配体与受体生物大分子相互作用，预测其亲和度和结合模式的过程，实现了基于受体结构的药物设计。然而基于传统优化算法及计算机的分子对接算法，面对指数大的分子及化合物空间，其复杂度较高，且在训练过程中会出现不稳定的现象。

技术实现思路

[0005]鉴于上述的分析，本专利技术实施例旨在提供一种分子对接药物筛选方法和装置，用以解决现有优化算法及计算机的分子对接算法，面对指数大的分子及化合物空间，其复杂度较高，且在训练过程中会出现不稳定的现象的问题。
[0006]一方面，本专利技术实施例提供了一种分子对接药物筛选方法，包括：将药物分子按其结合基团的特征抽象成结合部位节点和结合原件节点的集合，其中，受体分子包含结合部位节点，配体分子包含结合原件节点；根据结合基团之间的亲和度打分函数值构建结合部位和结合原件的亲和度矩阵，并列出约束条件和优化目标；根据亲和度矩阵和约束条件，对需要优化设计的相关变量进行简化；根据亲和度矩阵，基于横向独热模型构建量子比特分配与算子映射；根据映射得到的Mixer算子和Phaser算子，分层构建QAOA参数化线路；设置迭代次数，优化步长，通过损失函数的梯度下降对QAOA参数化线路进行训练；以及多次测量训练后的QAOA参数化线路，并获取目标函数最大的比特串作为结果。
[0007]上述技术方案的有益效果如下：利用量子计算在生物化学分子体系中进行准确、
快速计算的量子优势。分子对接药物筛选方法，将传统的0
‑
1规划问题转换为独热问题的组合，引入辅助优化变量，使得优化效率大大提高，减少量子比特，避免了引入过多额外的量子比特求解提升复杂度。针对分子对接算法，在指数大的分子及化合物空间中计算时复杂度较高、不稳定的现象，该方法不仅减少了训练参数，提升了表达性，并且充分运用量子电路并行处理数据的能力，是量子计算与药物发现相结合的重要应用。
[0008]基于上述方法的进一步改进，亲和度打分函数是配体分子与受体分子结合位点间的相互作用能的线性表示，其中，列出约束条件和优化目标包括：约束条件为配位相关约束条件；以及优化目标为使总亲和度最大，其中，当亲和度打分函数值为
‑
∞时，结合部位节点与结合原件节点无法结合。
[0009]基于上述方法的进一步改进，根据亲和度矩阵和约束条件，对优化设计的相关变量进行简化包括：
[0010]根据亲和度矩阵、目标函数以及约束条件，抽象出以下数学模型：给定矩阵W∈R
M
×
N
、集合B：＝{X1，X2，...}，对二值矩阵变量X∈{0，1}
M
×
N
进行求解，使得：
[0011][0011][0012][0013][0014][0015]其中，B表示额外限制条件的集合，X表示不允许同时配对的结合方式集合，w
ij
表示配体分子的第i结合原件与受体分子的第j结合部位基团间亲和度的打分函数结果，x
ij
表示配体分子的第i结合原件与受体分子的第j结合部位配对；以及当亲和度打分函数值为
‑
∞时，去除其对应逻辑量子比特并短路其横向连接。
[0016]基于上述方法的进一步改进，根据亲和度矩阵，基于横向独热模型构建量子比特分配与算子映射包括：通过增加全0列，将原始M
×
N维亲和度矩阵W重写为M
×
(N+1)维矩阵W
′
，根据重写矩阵W
′
的维度设置M
×
(N+1)位逻辑量子比特；将所有判定x
m，N+1
＝0的对应辅助量子比特q
m，N+1
删除；Mixer算子和Phaser算子分别为和H
P
＝H
S
+H
PV
+H
PB
，其中，为表示横向连接的集合，D
k，l；m，n
＝X
k，l
X
m，n
+Y
k，l
Y
m，n
，X
k，l
、X
m，n
、Y
k，l
、Y
m，n
是下标编号为k，l和m，n的哈密顿量Pauli
‑
X和Pauli
‑
Y算符，H
S
为编码原问题的目标函数，H
PV
为编码纵向惩罚项，H
PB
为编码额外限制条件；以及获得单层参数化量子线路U(β，γ)＝U
M
(β)U
P
(γ)，其中，U
M
(β)＝exp(
‑
iβH
M
)，U
P
(γ)＝exp(
‑
iγH
P
)。
[0017]基于上述方法的进一步改进，当亲和度矩阵W∈R
M
×
N
的维度满足列数大于等于行数时，使得额外增加的优化变量尽可能小；以及当亲和度矩阵W∈R
M
×
N
的维度不满足列数大于等于行数时，对亲和度矩阵进行转置，使得转置矩阵W
T
的维度满足列数大于等于行数。
[0018]基于上述方法的进一步改进，将Mixer算子改写为2组M个独立量子系统的哈密顿量的和：
[0019][0020]当实现exp(...

【技术保护点】

【技术特征摘要】
1.一种分子对接药物筛选方法，其特征在于，包括：将药物分子按其结合基团的特征抽象成结合部位节点和结合原件节点的集合，其中，受体分子包含所述结合部位节点，配体分子包含所述结合原件节点；根据所述结合基团之间的亲和度打分函数值构建结合部位和结合原件的亲和度矩阵，并列出约束条件和优化目标；根据所述亲和度矩阵和所述约束条件，对需要优化设计的相关变量进行简化；根据所述亲和度矩阵，基于横向独热模型构建量子比特分配与算子映射；根据映射得到的Mixer算子和Phaser算子，分层构建QAOA参数化线路；设置迭代次数，优化步长，通过损失函数的梯度下降对所述QAOA参数化线路进行训练；以及多次测量训练后的QAOA参数化线路，并获取目标函数最大的比特串作为结果。2.根据权利要求1所述的分子对接药物筛选方法，其特征在于，所述亲和度打分函数是所述配体分子与所述受体分子结合位点间的相互作用能的线性表示，其中，列出约束条件和优化目标包括：所述约束条件为配位相关约束条件；以及所述优化目标为使总亲和度最大，其中，当所述亲和度打分函数值为
‑
∞时，所述结合部位节点与所述结合原件节点无法结合。3.根据权利要求1所述的分子对接药物筛选方法，其特征在于，根据所述亲和度矩阵和所述约束条件，对优化设计的相关变量进行简化包括：根据所述亲和度矩阵、所述目标函数以及所述约束条件，抽象出以下数学模型：给定矩阵W∈R
M
×
N
、集合B：＝{X1，X2，
…
}，对二值矩阵变量X∈{0，1}
M
×
N
进行求解，使得：进行求解，使得：进行求解，使得：进行求解，使得：进行求解，使得：其中，B表示额外限制条件的集合，X表示不允许同时配对的结合方式集合，w
ij
表示所述配体分子的第i结合原件与所述受体分子的第j结合部位基团间亲和度的打分函数结果，x
ij
表示所述配体分子的第i结合原件与所述受体分子的第j结合部位配对；以及当所述亲和度打分函数值为
‑
∞时，去除其对应逻辑量子比特并短路其横向连接。4.根据权利要求1所述的分子对接药物筛选方法，其特征在于，根据所述亲和度矩阵，基于横向独热模型构建量子比特分配与算子映射包括：通过增加全0列，将原始M
×
N维亲和度矩阵W重写为M
×
(N+1)维矩阵W
′
，根据重写矩阵W
′
的维度设置M
×
(N+1)位逻辑量子比特；将所有判定x
m，N+1
＝0的对应辅助量子比特q
m，N+1
删除；
Mixer算子和Phaser算子分别为和H
P
＝H
S
+H
PV
+H
PB
，其中，为表示横向连接的集合，D
k，l；m，n
＝X
k，l
X
m，n
+Y
k，l
Y
m，n
，X
k，l
、X
m，n
、Y
k，l
、Y
m，n
是下标编号为k，l和m，n的哈密顿量Pauli
‑
X和Pauli
‑
Y算符，H
S
为编码原问题的目标函数，H
PV
为编码纵向惩罚项，H
PB
为编码额外限制条件；以及获得单层参数化量子线路U(β，γ)＝U
M
(β)U
P
(γ)，其中，U
M
(β)＝exp(
‑
iβH
M
)，U
P
(γ)＝exp(
‑
iγH
P
)。5.根据权利要求4所述的分子对接药物筛选方法，其特征在于，当所述亲和度矩阵W∈R
M
×
N
的维度满足列数大于等于行数时，使得额外增加的优化变量尽可能小；以及当所述亲和度矩阵W∈R
M
×
N
的维度不满足列数大于等于行数时，对所述亲和度矩阵进行转置，使得转置矩阵W
T
的维度满足列数大于等于行数。6.根据权利要求4所述的分子对接药物筛选方法，其特征在于，将所述Mixer算子改写为2组个独立量子系统的哈密顿量的和：当实现exp(
‑
iβH
M
)时，每个独立量子系统的哈密顿模拟并行实现，线路深度仅为O(1)；所述Phaser算子为H
P
＝H
S
+H
PV
+H
PB
，其中，H
S
编码原问题的目标函数为：所述编码纵向惩罚项H
PV
为：所述编码额外限制条件H
P...

【专利技术属性】
技术研发人员：张先超，张在琛，余旭涛，拱翟锐，
申请(专利权)人：嘉兴学院，
类型：发明
国别省市：