CN109461475B - 一种基于人工神经网络的分子属性预测方法 - Google Patents
一种基于人工神经网络的分子属性预测方法 Download PDFInfo
- Publication number
- CN109461475B CN109461475B CN201811258268.XA CN201811258268A CN109461475B CN 109461475 B CN109461475 B CN 109461475B CN 201811258268 A CN201811258268 A CN 201811258268A CN 109461475 B CN109461475 B CN 109461475B
- Authority
- CN
- China
- Prior art keywords
- characterization
- atom
- molecule
- molecular
- representation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种基于人工神经网络的分子属性预测方法,包括:S1)对分子数据进行预处理:通过图的数据结构表征的方法,得到原子空间表征与原子构成表征;S2)建立模型:将原子空间表征与原子构成表征通过多层卷积神经网络,得到分子各级的表征,并将分子各级的表征进行组合,得到模型;S3)根据模型预测分子属性。与现有技术相比,本发明利用多层级卷积神经网络,能够利用已有数据的信息以及分子的多层级结构,从中学出分子属性和空间构成的关系,并用来预测未知分子的相关属性,因此具有较好的速度与精度。
Description
技术领域
本发明属于材料学技术领域,尤其涉及一种基于人工神经网络的分子属性预测方法。
背景技术
从药物研发,到材料开发,都离不开分子发现。为了寻找到具有特定属性的分子来满足应用上的需求,通用的方法是遍历一个未知的可能分子的集合(称之为化学空间),在遍历过程中,研究人员通过各种方法来预测分子的属性,如果发现某分子符合要求,则记录下来做进一步研究。举例而言,对分子能量属性的预测能够帮助研究人员找到稳定的分子。
但是,化学空间往往非常庞大,一个被广泛应用的化学空间有一千六百多亿个分子。因此,一种快速的分子属性测定方法能够极大加速寻求特定分子的进程。但是传统物理上常用的密度泛函分析方法(Density Function Theory,下面简称为DFT),在时间效率上并不理想,无法对大量数据进行处理。
围绕该问题,研究者们提出了很多方法,但其中大部分仍然是基于DFT开发的。与此同时,该研究领域已经积累了许多相关的数据,但是大部分方法都无法利用这些已有的数据。
发明内容
有鉴于此,本发明要解决的技术问题在于提供一种基于人工神经网络的分子属性预测方法,该分子属性预测方法可充分利用已有的数据,且具有较好的速度与精度。
本发明提供了一种基于人工神经网络的分子属性预测方法,包括:
S1)对分子数据进行预处理:通过图的数据结构表征的方法,得到原子空间表征与原子构成表征;
S2)建立模型:将原子空间表征与原子构成表征通过多层卷积神经网络,得到分子各级的表征,并将分子各级的表征进行组合,得到模型;
S3)根据模型预测分子属性。
优选的,所述步骤S1)具体为:
分子数据包括分子的原子构成与原子的三维空间坐标;
将原子的三维空间坐标转换为原子之间的距离矩阵,然后通过径向基函数扩展为距离张量,得到原子空间表征;
将分子的原子构成进行嵌入表示,得到原子构成表征。
优选的,距离矩阵通过径向基函数扩展为距离张量,具体按照以下模型进行:
通过径向基函数得到张量D∈RN×N×K,N为距离矩阵中每对距离的数量。
优选的,将分子的原子构成进行嵌入表示,具体按照以下步骤进行:
用向量表示分子中的节点和边,将原子看做分子中的节点,其表示为向量a0∈RD,进而得到整个分子的节点表示矩阵A0∈RN×D,经嵌入表示得到边表征矩阵为E∈RN×D×D。
优选的,将原子空间表征与原子构成表征通过多层卷积神经网络,得到分子各级的表征具体为:
优选的,其特征在于,
优选的,所述步骤S2)中将分子各级的表征进行组合,具体为:
优选的,所述模型为:
优选的,得到模型后,还包括:用已有的分子属性数据训练模型中的参数。
优选的,所述步骤S2)中多层卷积神经网络中卷积层的个数为4~5。
本发明提供了一种基于人工神经网络的分子属性预测方法,包括:S1)对分子数据进行预处理:通过图的数据结构表征的方法,得到原子空间表征与原子构成表征;S2)建立模型:将原子空间表征与原子构成表征通过多层卷积神经网络,得到分子各级的表征,并将分子各级的表征进行组合,得到模型;S3)根据模型预测分子属性。与现有技术相比,本发明利用多层级卷积神经网络,能够利用已有数据的信息以及分子的多层级结构,从中学出分子属性和空间构成的关系,并用来预测未知分子的相关属性,因此具有较好的速度与精度。
在预测的速度上有了很大的提高:对于一个中等大小的分子,在同样的计算资源下,本发明仅需2.4×10-2秒,而传统基于密度泛函分析方法的模型需要3.6×103秒。相比之下本发明的效率高了大约1.5×105倍。
在预测的精度方面,本发明在多项属性的预测上都达到了目前最高的精度。
附图说明
图1为本发明提供的分子属性预测方法的流程示意图。
具体实施方式
下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供了一种基于人工神经网络的分子属性预测方法,包括:
S1)对分子数据进行预处理:通过图的数据结构表征的方法,得到原子空间表征与原子构成表征;
S2)建立模型:将原子空间表征与原子构成表征通过多层卷积神经网络,得到分子各级的表征,并将分子各级的表征进行组合,得到模型;
S3)根据模型预测分子属性。
参见图1,图1为本发明提供的分子属性预测方法的流程示意图。
分子数据包括分子的原子构成和每个原子的三维空间坐标;所述分子数据为本领域技术人员熟知的分子数据即可,并无特殊的限制,其可优选为由密度泛函分析方法计算得到的分子数据;对分子数据进行预处理,优选包括分子空间信息的处理和分子结构中节点和边的嵌入表征,前者在原始分子数据中表现为分子中每个原子的三维空间坐标,后者在原始分子数据中表现为分子的原子构成,本发明通过图的数据结构表征的方法,得到原子空间表征与原子构成表征。
其中,所述原子空间表征优选按照以下步骤进行:将原子的三维空间坐标转换为距离矩阵,然后通过径向基函数扩展为距离张量,得到原子空间表征。将原子的三维空间坐标转换为原子之间的距离矩阵,可消除坐标系选择对分子观测的影响,因此使本发明提供的模型具有旋转和平移不变性,不会被不同视角下的同个分子迷惑;然后用径向基函数将该距离矩阵拓展为一个距离张量。径向基函数是一种被广泛使用的核方法,给定一个集合中K个中心点{μ1,...μk},对于单个数据点x(即距离矩阵中的每对距离),按照以下模型进行处理:
其中,x为距离矩阵中的每对距离;∩表示对元素的拼接;μi为中心点;||x-μi||表示x到第i个中心点的欧几里得距离;h为径向基,本发明中优选为高斯核,即使用高斯核的好处在于可以避免训练初期较长的平台期;K为最短距离到最长距离选取中心点的数量,本发明中优选为在最短距离到最长距离中平均选K个点,这样,所有的距离都能被覆盖在其中。
通过径向基函数得到张量D∈RN×N×K,dij表示第i个原子与第j个原子之间的距离张量,N为距离矩阵中每对距离的数量。
预处理中另一部分的操作为嵌入操作,将分子的原子构成进行嵌入表示,得到原子构成表征,可分为原子表征与边表征。在嵌入表示时,用向量表示分子中的节点和边,将原子看做分子中的节点,其表示为向量a0∈RD,进而得到整个分子的节点表示矩阵即原子表征A0∈RN×D。在嵌入表示时,分子中同类原子优选共享一个嵌入值,如一个分子中所有氢原子的向量表示相同,而且该分子中的氢原子向量表示和其他分子中的氢原子表示也相同,因此,向量表示即原子表征仅仅与原子的种类相关,而和原子所在的分子构成以及空间结构没有关系。基于上述同样的技术,对分子中的边进行嵌入,优选当且仅当一个边的两个节点种类相同时,它们的向量表示相同,如所有连接着氢原子与氧原子的边都用同一个向量进行表示,经嵌入表示得到边表征矩阵为E∈RN×D×D,其中,eij∈R为第i个原子与第j个原子之间的边表征。
将原子空间表征与原子构成表征通过多层卷积神经网络,得到分子各级的表征。在此过程中,首选需要构建若干个卷积层,卷积层的个数是一个可以调节的超参数,通常在预测中小型分子时,优选4~5层。在本发明中,每层卷积的输入优选为边嵌入、空间信息以及上一层的相互作用表征,如在第一层的输入时,相互作用表征就是原子表征,而经过第一层卷积后,得到二级相互作用表征来作为第二个卷积层的输入,在每一层嵌入之后,可以得到更高一层的相互作用的表征。考虑到分子中的量子作用力可以被转换为多个级别,如仅基于原子的,两两原子之间的,以及每三个原子之间等等,因此通过多层卷积神经网络,得到分子各级的表征优选具体为:
其中,he是用来更新边表示的,hv是用来生成更高层的原子相互作用表示的。根据此多层级的建模,能够高效的保留分子的结构信息并描述其量子相互作用力。具体而言,在第一个相互作用层即第一卷积层中,输入的是预处理得到的原子表征a0,此表征代表了不同化学元素的原子的固有属性,在第一相互作用层中,模型给改表征加入了一阶邻居节点、边以及空间信息,则输出了a1,a1表示二级的相互作用,即每两个原子之间的相互作用表征。以类似的方式,可计算出a2表示三阶相互作用,即每三个原子之间的表征,以此类推,还可以得到a3、a4直到aT,T代表相互作用层的个数,即卷积层的层数。
在本发明中,更新函数he优选具体为:
所述更新函数hv优选具体为:
其中,σ为tanh激活函数,Wuv为权重矩阵,Mfa、Mfd与Mfe为全连接层,即M(x)=Wx+b。通过该函数,模型实行了卷积操作以产生更高阶的原子表示,这里距离张量用以控制相互作用影响的强度,而边嵌入提供了无法被映射到节点上的额外信息。
本发明采用多层级的结构能够给模型带来可迁移的优点。由于传入第一级相互作用层的原子表征和边嵌入仅仅和原子及边的种类有关而和特定的分子结构以及空间信息无关,因此,本发明得到的模型在原子表征和边嵌入中学到的化学领域的知识在分子系统中是通用的。接着,在本发明实例多层级的模型中,这些表征被用来生成更深层次的相互作用表征。虽然大分子和小分子在原子和边的分布上有不同,但是由于分子结构的局部性和可分解性,它们相互作用的方式是类似的。通过这些通用的表征和相似的相互作用机制,本发明提供的模型能够以较高的精度推断出高阶的相互作用表征。因此,本发明提供的模型能够将从小分子数据上学习到的知识迁移到大分子的属性预测上去。此外,另一种迁移知识的方法是,利用小分子数据来学习原子和边的表征,并用学习到的值作为另一个模型的初始化来取代上述的随机初始化方法。
另一方面,由于本发明提供的模型中嵌入层的操作赋予了表征以平移不变性和旋转不变性。加上模型在相互作用层的上述操作是基于元素的,因此,模型的预测结果和原子的顺序无关。这种次序不变性增强了模型的泛化性。
在得到分子各级的表征之后,将其进行组合,优选按照下式进行:
通过分子各级的表征组合,得到模型,优选具体为:
由于分子的属性具有可加性和局部性,模型能够通过对每个原子各自进行预测最终加和来获得预测值。在本发明中,模型的前一项表示映射到每个原子的量子相互作用,而额外的和边相关的信息由最后一项提供。由于和原子相关的相互作用占了分子内相互作用的大部分,而与边相关的部分很少,因此,当分子较小时,可以忽略模型公式中的后一项。
按照本发明本发明,优选使用小批次随机梯度下降算法来更新模型参数;所述小批次随机梯度下降算法为本领域技术人员熟知的算法即可,并无特殊的限制,本发明中优选算法中使用的优化器为自适应力矩优化器(Adam optimizer)。在本发明提供的一些实施例中,小批量的数值优选为被设定64,而初始学习率为10-5;在本发明提供的一些实施例中,处理时可以根据数值设定一个阈值,超过这个阈值的距离都用掩码为0,忽略这对距离对应原子的相互作用,由此来减少模型的复杂度并增强鲁棒性。
在本发明提供的一些实施例中,在训练时,使用80%的数据作为训练集,使用10%的数据作为验证集,使用剩余10%的数据作为测试集。验证集用以进行调节超参数,确定好超参数后,如果损失函数连续十次全量数据迭代后都没有变小,则停止训练,并选择在测试集上表现最好的模型。需要注意的是,对于分子的不同属性,需要单独训练不同的模型可以取得比联合训练更好的效果。
本发明依照上述方案,可以充分利用已有数据,学习到分子系统内多层级相互作用的表征,并用以高效地预测未知分子的属性,从而能加速材料学,医学等领域上特性属性分子的搜寻。
通过本发明提供的预测方法能够预测的属性包括但不限于:U0(温度为0K时分子内能量)、U(温度为298.15K时分子内能量)、H(温度为298.15K时分子热焓)、G(温度为298.15K时分子自由能)、Cv(温度为298.15K时分子热容量)、zpve(基本振动能)、R2(电子空间广度)、EHOMO(最高占据分子轨道能量)、ELUMO(最低未占分子轨道)、μ(偶极矩)与α(等向极化率)
为了进一步说明本发明,以下结合实施例对本发明提供的一种基于人工神经网络的分子属性预测方法进行详细描述。
以下实施例中所用的试剂均为市售。
实施例1
以下面三个分子的U0预测为例,它们都来自国际上常用的QM9数据集,单位为eV。我们以QM9数据集作为训练集,训练上方式如上所述,然后利用训练后得到的模型对以下分子作预测。衡量误差的标准采取绝对误差,即预测值和真实值差的绝对值。
(1)CH4的真实值为-17.1717476eV,预测值为-17.1681695eV,误差为0.0035781eV。(2)NH3的真实值为-12.0055513eV,预测值为-12.0187658eV,误差为0.0132145eV。
(3)HOH的真实值为-9.2401279eV,预测值为-9.2371538eV,误差为0.0029741eV。
并且在下表中给出了本方法在整个QM9数据集上各个属性预测的平均误差。
Claims (8)
1.一种基于人工神经网络的分子属性预测方法,其特征在于,包括以下步骤:
S1)对分子数据进行预处理:通过图的数据结构表征的方法,得到原子空间表征与原子构成表征;
S2)建立模型:将原子空间表征与原子构成表征通过多层卷积神经网络,得到分子各级的表征,并将分子各级的表征进行组合,得到模型;
S3)根据模型预测分子属性;
所述步骤S1)具体为:
分子数据包括分子的原子构成与原子的三维空间坐标;
将原子的三维空间坐标转换为原子之间的距离矩阵,然后通过径向基函数扩展为距离张量,得到原子空间表征;
将分子的原子构成进行嵌入表示,得到原子构成表征;
距离矩阵通过径向基函数扩展为距离张量,具体按照以下模型进行:
通过径向基函数得到张量D∈RN×N×K,N为距离矩阵中每对距离的数量。
2.根据权利要求1所述的分子属性预测方法,其特征在于,将分子的原子构成进行嵌入表示,具体按照以下步骤进行:
用向量表示分子中的节点和边,将原子看做分子中的节点,其表示为向量a0∈RD,进而得到整个分子的节点表示矩阵A0∈RN×D,经嵌入表示得到边表征矩阵为E∈RN×D×D。
7.根据权利要求6所述的分子属性预测方法,其特征在于,得到模型后,还包括:用已有的分子属性数据训练模型中的参数。
8.根据权利要求1所述的分子属性预测方法,其特征在于,所述步骤S2)中多层卷积神经网络中卷积层的个数为4~5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811258268.XA CN109461475B (zh) | 2018-10-26 | 2018-10-26 | 一种基于人工神经网络的分子属性预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811258268.XA CN109461475B (zh) | 2018-10-26 | 2018-10-26 | 一种基于人工神经网络的分子属性预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109461475A CN109461475A (zh) | 2019-03-12 |
CN109461475B true CN109461475B (zh) | 2022-05-13 |
Family
ID=65608537
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811258268.XA Active CN109461475B (zh) | 2018-10-26 | 2018-10-26 | 一种基于人工神经网络的分子属性预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109461475B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111798935A (zh) * | 2019-04-09 | 2020-10-20 | 南京药石科技股份有限公司 | 基于神经网络的普适性化合物结构-性质相关性预测方法 |
US11798655B2 (en) | 2019-05-15 | 2023-10-24 | International Business Machines Corporation | Feature vector feasibility estimation |
US11817184B2 (en) | 2019-05-16 | 2023-11-14 | Robert Bosch Gmbh | Graph neural network force field computational algorithms for molecular dynamics computer simulations |
CN110767271B (zh) * | 2019-10-15 | 2021-01-08 | 腾讯科技(深圳)有限公司 | 化合物性质预测方法、装置、计算机设备及可读存储介质 |
CN110957012B (zh) * | 2019-11-28 | 2021-04-09 | 腾讯科技(深圳)有限公司 | 化合物的性质分析方法、装置、设备及存储介质 |
CN112162182A (zh) * | 2020-09-28 | 2021-01-01 | 哈尔滨理工大学 | 一种基于神经网络的气体介电强度预测方法 |
CN112147473A (zh) * | 2020-09-28 | 2020-12-29 | 哈尔滨理工大学 | 一种高绝缘强度气体的筛选方法 |
CN113284565B (zh) * | 2021-05-18 | 2023-09-22 | 百度时代网络技术(北京)有限公司 | 信息处理的方法和装置 |
CN113409893B (zh) * | 2021-06-25 | 2022-05-31 | 成都职业技术学院 | 一种基于图像卷积的分子特征提取及性能预测方法 |
CN115171807B (zh) * | 2022-09-07 | 2022-12-06 | 合肥机数量子科技有限公司 | 一种分子编码模型训练方法、分子编码方法和系统 |
CN115713986B (zh) * | 2022-11-11 | 2023-07-11 | 中南大学 | 基于注意力机制的材料晶体属性预测方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104978498A (zh) * | 2015-04-16 | 2015-10-14 | 上海大学 | 生物分子网络拓扑结构比对的自适应方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0718027D0 (en) * | 2007-09-14 | 2007-10-24 | Univ Manchester | Method for determining three-dimensional structures of dynamic molecules |
CN101587510A (zh) * | 2008-05-23 | 2009-11-25 | 中国科学院上海药物研究所 | 基于复杂抽样和改进决策森林算法的化合物致癌毒性预测方法 |
EP2596352B1 (en) * | 2010-07-22 | 2016-08-31 | Ramot at Tel Aviv University, Ltd. | Systems and methods for detection of cellular stress |
WO2015017798A2 (en) * | 2013-08-02 | 2015-02-05 | CRIXlabs, Inc. | Method and system for predicting spatial and temporal distributions of therapeutic substance carriers |
-
2018
- 2018-10-26 CN CN201811258268.XA patent/CN109461475B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104978498A (zh) * | 2015-04-16 | 2015-10-14 | 上海大学 | 生物分子网络拓扑结构比对的自适应方法 |
Non-Patent Citations (3)
Title |
---|
A deep learning architecture for molecules and materials;K. T. Schütt 等;《The Journal of Chemical Physics》;20180329;第148卷(第24期);第1-11页 * |
Edge Attention-based Multi-Relational Graph Convolutional Networks;Chao Shang 等;《arXiv:1802.04944v1》;20180214;第1-10页 * |
多尺度复合材料力学研究进展;陈玉丽 等;《固体力学学报》;20180228;第39卷(第1期);第1-68页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109461475A (zh) | 2019-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109461475B (zh) | 一种基于人工神经网络的分子属性预测方法 | |
Zhang et al. | An end-to-end deep learning architecture for graph classification | |
Zhang et al. | Feature selection algorithm based on bare bones particle swarm optimization | |
CN107862173B (zh) | 一种先导化合物虚拟筛选方法和装置 | |
CN111724867B (zh) | 分子属性测定方法、装置、电子设备及存储介质 | |
Young et al. | Optimizing deep learning hyper-parameters through an evolutionary algorithm | |
Ayvaz | Simultaneous determination of aquifer parameters and zone structures with fuzzy c-means clustering and meta-heuristic harmony search algorithm | |
Mansinghka et al. | Structured priors for structure learning | |
Peng et al. | Accelerating minibatch stochastic gradient descent using typicality sampling | |
CN103116324A (zh) | 基于指标预测和在线学习的微电子生产线调度方法 | |
Yan | Weighted K-nearest neighbor classification algorithm based on Genetic Algorithm | |
Bej et al. | A multi-schematic classifier-independent oversampling approach for imbalanced datasets | |
Ma | An Efficient Optimization Method for Extreme Learning Machine Using Artificial Bee Colony. | |
CN109063418A (zh) | 疾病预测分类器的确定方法、装置、设备及可读存储介质 | |
Naik et al. | A global-best harmony search based gradient descent learning FLANN (GbHS-GDL-FLANN) for data classification | |
Tuli et al. | FlexiBERT: Are current transformer architectures too homogeneous and rigid? | |
Chen et al. | A weighted bagging LightGBM model for potential lncRNA-disease association identification | |
Yan et al. | A review about RNA–protein-binding sites prediction based on deep learning | |
Rahman et al. | IDMIL: an alignment-free Interpretable Deep Multiple Instance Learning (MIL) for predicting disease from whole-metagenomic data | |
Bharathi et al. | Optimal feature subset selection using differential evolution and extreme learning machine | |
Bichat et al. | Hierarchical correction of p-values via an ultrametric tree running Ornstein-Uhlenbeck process | |
Liu et al. | Improving genetic classifiers with a boosting algorithm | |
CN110162704B (zh) | 基于多因子遗传算法的多规模关键用户提取方法 | |
JP6993250B2 (ja) | コンテンツ特徴量抽出装置、方法、及びプログラム | |
CN113223622A (zh) | 基于元路径的miRNA-疾病关联预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |