WO2021103491A1 - 一种测试和拟合力场二面角参数的方法 - Google Patents
一种测试和拟合力场二面角参数的方法 Download PDFInfo
- Publication number
- WO2021103491A1 WO2021103491A1 PCT/CN2020/096122 CN2020096122W WO2021103491A1 WO 2021103491 A1 WO2021103491 A1 WO 2021103491A1 CN 2020096122 W CN2020096122 W CN 2020096122W WO 2021103491 A1 WO2021103491 A1 WO 2021103491A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- dihedral angle
- force field
- dihedral
- parameters
- fitting
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000012360 testing method Methods 0.000 title claims abstract description 9
- 239000012634 fragment Substances 0.000 claims abstract description 15
- 229920002521 macromolecule Polymers 0.000 claims abstract description 13
- 238000012795 verification Methods 0.000 claims abstract description 3
- 238000001311 chemical methods and process Methods 0.000 claims description 4
- 238000012617 force field calculation Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 15
- 238000001514 detection method Methods 0.000 abstract description 8
- 150000002605 large molecules Chemical class 0.000 abstract description 5
- 150000003384 small molecules Chemical class 0.000 description 5
- 239000000126 substance Substances 0.000 description 4
- 238000000324 molecular mechanic Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000004057 DFT-B3LYP calculation Methods 0.000 description 2
- 238000005452 bending Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000003077 quantum chemistry computational method Methods 0.000 description 2
- 238000003775 Density Functional Theory Methods 0.000 description 1
- 101100425949 Mus musculus Tnfrsf13c gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 238000009510 drug design Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000005624 perturbation theories Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N10/00—Quantum computing, i.e. information processing based on quantum-mechanical phenomena
- G06N10/20—Models of quantum computing, e.g. quantum circuits or universal quantum computers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N10/00—Quantum computing, i.e. information processing based on quantum-mechanical phenomena
- G06N10/60—Quantum algorithms, e.g. based on quantum optimisation, quantum Fourier or Hadamard transforms
Definitions
- the present invention belongs to the field of molecular mechanics, and specifically relates to a method for testing and fitting dihedral angle parameters of a force field, which is suitable for evaluating dihedral angle parameters in a molecular force field and simulating parameters that perform poorly. He revised. Background technique
- Molecular mechanics is widely used in many fields such as drug design due to its speed advantage and reliable accuracy. Molecular mechanics is based on formulas that describe molecular properties (such as energy) and corresponding parameters, and the parameter part is called the molecular force field.
- E E bond + E angle + E dlhedral + E i ni-CT + E ele + E vdw
- E bcmd the energy determined by the bond length of the two connected atoms
- E angle I the energy of the angle determined by the three connected atoms
- E di ( te ( tal S is the energy of the dihedral angle determined by the four connected atoms
- E inipropra ⁇
- E ⁇ is the energy between the charges of two atoms
- E vdw is the van der Waals energy of the two atoms.
- calculate the bond, angle, dihedral angle, out-of-plane bending, charge and Van der Waals parameters form the force field.
- the dihedral angle parameter exceeds other parameters in quantity and flexibility, so the quality of the dihedral angle parameter is very important to the overall quality of the force field.
- the development of a molecular force field is generally based on smaller molecular fragments, and then a quantum chemistry calculation is performed and the calculation result (usually using energy) as a target is fitted to obtain a series of force field parameters.
- the evaluation of the performance of the force field parameters usually takes the quantum chemistry calculation data of some other small molecules that are not within the fitting range as the standard. Small molecules can perform a large number of accurate quantum chemical calculations.
- the selected quantum chemistry method is generally density functional theory or a high-precision method based on perturbation theory. These small molecules generally have 1-2 flexible dihedral angles. Usually dihedral scanning is to rotate these dihedral angles and scan at regular angles within the range of -180 degrees -180 degrees. During the scanning process, the specific dihedral angle is fixed at a specific angle (generally 15 degrees) for structural optimization. Finally, compare the force field with the energy obtained by quantum chemistry, and evaluate the force field.
- the present invention provides a detection and fitting method suitable for the dihedral angle parameters of the macromolecular force field, which requires less calculation than the traditional method.
- these structures represent different angles of the flexible dihedral angle in the molecule; compare these structures with force field and quantum chemical methods, if they meet the standards, it is deemed that the force field parameters are satisfactory , End the process;
- the macromolecule is further cut into molecular fragments containing only one flexible dihedral angle, and the dihedral angle scan is performed.
- the quantum chemistry results of each flexible dihedral angle with the force field results, find Find out the flexible dihedral angle that does not meet the standard, and fit its parameters; after obtaining the new dihedral angle parameter, return it to the original series of structures of the entire molecule for verification. If it meets the standard, end the entire process , Complete the detection and fitting of poorly performing dihedral parameters; if it does not meet the standard, perform a dihedral scan of the entire molecule for the poorly performing flexible dihedral angle.
- the present invention is preferably R greater than 0.7 and d E á 2.0 kcal/mol (the two The threshold standard can be lowered or raised according to the specific requirements of the user, and the standard in the following steps is the same setting), terminate the process, and the force field parameters perform better in this molecule. Otherwise, go to step (4);
- each fragment contains a flexible dihedral angle (applicable to rdkit's rdkit. Chem. rdmolops. FragmentOnBonds() function). Perform the usual dihedral scan of these fragments, compare the quantum chemistry data and the force field data, for the dihedral that does not perform well (the correlation coefficient R is less than 0.7 or the energy deviation dE is greater than 2.0 kcal/mol) Angle, use the obtained quantum chemistry data to fit the dihedral angle parameter to obtain a new dihedral angle parameter.
- step (5) Using the new dihedral angle parameter obtained in step (4), repeat the force field calculation in step (2) to obtain a new energy E MM .
- the quantum chemistry previously calculated in step (2) Data EQM does not need to be calculated repeatedly, and then E _ and £ are linearly fitted to obtain the correlation coefficient R'and energy deviation dE' of the two sets of data
- a method for testing and fitting dihedral angle parameters of a force field provided by the present invention has the following technical advantages: [0018]
- the present invention proposes a new method for detecting and fitting dihedral angle parameters for macromolecules This method saves more calculations than the usual dihedral scanning of the entire molecule. By combining with the usual process, the detection and fitting of the dihedral angle parameters of the force field with macromolecules can be completed with lower computational resource consumption.
- FIG. 1 is a general flow of detection and fitting of macromolecules.
- FIG. 2 is a flowchart of the present invention.
- the lower left part is a general method based on dihedral scanning of the entire molecule. After being combined with the new method of the present invention, these steps only need to be applied to the dihedral angle that does not meet the set standard in the previous step.
- Inventive embodiment of the present invention is a general method based on dihedral scanning of the entire molecule. After being combined with the new method of the present invention, these steps only need to be applied to the dihedral angle that does not meet the set standard in the previous step.
- the force field method here the published GAFF2 force field is used
- the time can be ignored . So the CPU time required for this step is 60 hours.
- the calculation required by the new detection and fitting method of the present invention is far less than the usual dihedral scan of the entire molecule.
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Artificial Intelligence (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种测试和拟合力场二面角参数的方法，产生一些具有代表性的构象，对结构进行力场和量子化学方法的比较，如符合标准，就认定力场参数表现满意，结束流程；如不符合标准，将大分子切成只含一个柔性二面角的分子碎片，进行二面角扫描，通过对每个柔性二面角的量子化学结果和力场结果对比，找出不符合标准的柔性二面角，对其参数进行拟合；得到新的二面角参数后，返回到最初产生的整个分子的一系列结构进行验证，如符合标准，结束整个流程，完成检测和拟合表现不好的二面角参数；如不符合标准，将表现不好的柔性二面角进行整个分子的二面角扫描。上述方法以更低的计算资源消耗完成用大分子对力场二面角参数进行检测和拟合。
Description
一种 测试 和拟 合力 场二 面 角参 数的 方法 技术领域
[0001] 本发明属于分子力学领域， 具体涉及一种测试和拟合力场二面角参数的方法， 适用于对于对分子力场中二面角参数进行评价和对于表现不好的参数进行拟合 修正。 背景技术
[0002] 分子力学因其速度优势和可靠的精度被广泛应用于很多领域比如药物设计。 分 子力学基于描述分子性质 （比如能量） 的公式和相对应的参数， 其中的参数部 分被成为分子力场。 常用的力场能量的定义为： E=E bond+E angle+E dlhedral+ E ini-CT+E ele+E vdw，其中 E bcmd是相连的两个原子的键长决定的能量， E angle是相 连的三个原子决定的角度的能量， E di（te（talS相连的四个原子决定的二面角的能 量， E inipropra^维持一个平面的四个原子在同一平面的面外弯曲能量， E ^是两 个原子的电荷之间的能量， E vdw是两个原子的范德华能量。 与之相对应的计算 各项能量所需要的键， 角， 二面角， 面外弯曲， 电荷和范德华参数组成力场。 二面角参数在数量上和灵活性上超过其他参数， 因此二面角参数的质量对于力 场的整体质量至关重要。
[0003] 发展分子力场一般基于较小的分子碎片， 然后进行量子化学计算并将计算结果 （一般用能量） 作为目标进行拟合， 从而得到一系列力场参数。 在此之后， 力 场参数表现的评价通常将一些不在拟合范围内的其他小分子的量子化学的计算 数据作为标准。 小分子可以进行大量准确的量子化学计算。 选取的量子化学方 法一般为密度泛函或者基于微扰理论的高精度的方法。 这些小的分子一般有 1-2 个柔性二面角。 通常二面角扫描是旋转这些二面角， 并且在 -180度 -180度范围 内每隔一定角度进行扫描。 在扫描过程中， 固定特定二面角在特定角度（一般 15 度）， 进行结构优化。 最后比较力场和量子化学得到的能量， 对力场进行评价。
[0004] 如上所述， 发展力场一般基于小的分子， 但在实际应用比如药物分子设计中， 分子往往比较大 （柔性角在 3个以上） ， 这就需要力场参数特别是二面角参数有
好的从小分子到大分子的迀移性。 通常的检测二面角参数和拟合新的二面角参 数的方法是上述对整个分子二面角扫描的方法 (流程见图 1 ) ， 由于大的分子含 有柔性二面角比较多， 而且二面角之间可能存在耦合， 需要对耦合的二面角进 行组合扫描， 所以对于大分子来说， 传统的方法需要进行大量的高精度量子化 学计算。 发明概 述 技术问题 问题的解 决方案 技术解决方案
[0005] 针对上述技术问题， 本发明提供一种适用于大分子力场二面角参数的检测和拟 合方法， 这种方法比传统方法需要更少的计算量。
[0006] 具体技术方案为：
[0007] 一种测试和拟合力场二面角参数的方法
[0008] 首先产生一些具有代表性的一部分构象， 这些结构代表分子中柔性二面角的不 同角度； 对这些结构进行力场和量子化学方法的比较， 如果符合标准， 就认定 力场参数表现满意， 结束流程；
[0009] 如果不符合标准， 进一步将大分子切成只含一个柔性二面角的分子碎片， 进行 二面角扫描， 通过对每个柔性二面角的量子化学结果和力场结果对比， 找出不 符合标准的柔性二面角， 对其参数进行拟合； 得到新的二面角参数后， 将其返 回到最初的产生的整个分子的一系列结构进行验证， 如果符合标准， 结束整个 流程， 完成检测和拟合表现不好的二面角参数； 如果不符合标准， 将表现不好 的柔性二面角进行整个分子的二面角扫描。
[0010] 具体步骤为：
[0011] ( 1 ) 对于一个大分子， 首先使用 rdkit对于每一个分子产生 500个构象 (具体 命令为 rdkit. EmbedMultipleConfs (mol, 500) )， 将这些分子用 rdkit自带的 UFF 力场进行结构优化 (rdkit. UFFGetMoleculeForceField (mol) . Minimize 0 )， 计 算每个结构的每个柔性二面角的角度。 根据柔性二面角的角度分布， 选取 30个 结构， 涵盖从 -180度到 180度不同的区域， 因为实际应用中， 能量较低的结构占
有更大的比重， 所以柔性二面角角度相同时， 优先选择能量较低的结构；
[0012] （2） 使用量子化学计算软件 （比如 PSI4） ，将步骤 （1） 中的结构使用高精度 量子化学方法 （比如 B3LYP/6-31G（d））进行进一步结构优化， 得到相应的能量为 E _， 同时用需要检测的力场将这些结构进行优化， 得到对应的能量为 E
[0013] （3） 将步骤 （2） 中得到的两组能量以分子为单位进行线性拟合 （python的 sc ipy模块， R, A E=scipy. stats, linregress（E QM, E丽 ）， 得到两组数据的 Pearson相关系数 R和能量偏差 dE， 如果 R大于第一阈值而且 d E小于第二阈值， 本发明优选为 R大于 0. 7而且 d E á2. 0千卡 /摩尔 （这两个阈 值标准可以根据用户的具体要求降低或者升高， 下面的步骤里的标准同此设定 ） ， 终止流程， 力场参数在此分子表现较好。 否则， 进入步骤 （4） ;
[0014] （4） 将步骤 （3） 中进入此步骤的分子切成较小的碎片， 每个碎片包含一个柔 性二面角（适用 rdkit的 rdkit. Chem. rdmolops. FragmentOnBonds（）函数） 。 对这 些碎片进行通常的二面角扫描， 对比量子化学数据和力场数据， 对于表现不好 （两者相关系数 R小于 0. 7或者能量偏差 dE大于 2. 0千卡 /摩尔） 的二面角， 利用 得到的量子化学数据对二面角参数进行拟合， 得到新的二面角参数。
[0015] （5） 用步骤 （4） 得到的新的二面角参数， 重复步骤 （2） 中的力场计算得到 新的能量 E MM，. 此时步骤 （2） 中以前计算的量子化学数据 EQM无需重复计算， 然后将 E _和£ 进行线性拟合， 得到两组数据的相关系数 R’ 和能量偏差 dE’
。 如果 R ’ 大于 0. 7而且 d E ’ á2. 0千卡 /摩尔， 终止流程， 新拟合的二面角参数 在此分子表现较好。 否则， 进入步骤 （6） ;
[0016] （6） 将经过步骤 （4） - 步骤 （5） 中仍然表现不好的参数进行通常的整个分 子的二面角扫描进行拟合过程。 可以先扫描步骤 （4） 中表现不好的二面角， 然 后进行二面角参数拟合， 将新的拟合参数用最初产生的结构进行力场的计算， 和步骤 （2） 中已有的量子化学结果 E QM进行比较， 如果 R大于 0. 7而且 d E á2. 0 千卡 /摩尔， 终止流程。 否者， 接着进行其他柔性二面角的扫描， 拟合相关二面 角参数。 发明的 有益效果 有益效果
[0017] 本发明提供的一种测试和拟合力场二面角参数的方法， 具有以下技术优势： [0018] 本发明提出的运用于大分子的新的检测和拟合二面角参数的方法， 比通常的对 于整个分子的二面角扫描更加节省计算量。 通过与通常的流程结合使用， 以更 低的计算资源消耗完成用大分子对力场二面角参数进行检测和拟合。 对附图 的简要说 明 附图说明
[0019] 图 1为现有通用的大分子检测和拟合的流程。
[0020] 图 2为本发明流程图。 左下部分为基于整个分子二面角扫描的通用方法。 与本 发明的新方法结合后， 这些步骤只需要应用于其上步中表现不符合设定标准的 二面角。 发明实施 例 本发明的实施方 式
[0021] 结合实施例说明本发明的具体技术方案。
[0022] 如图 2所示的流程：
[0023] 以大分子（SMILES: C0C1=C（0CCCC2=NC3=CC=CC=C3 [N] 2C） C=C4N=CN=CC4=C1）为例 ， 该分子有 6个柔性二面角， 如果用通常的整个分子二面角扫描， 需要优化 600 个构象， 每个结构优化 （量子化学方法使用 B3LYP/6-31G（d）） 需要的 CPU时间 （ 单个 CPU核心需要计算的时间） 为 2小时， 相对于量子化学计算的时间， 力场计 算和力场拟合的时间可以忽略不计 （在下面的方法里， 也将只对比量子化学计 算的时间） ， 所以共需要 CPU时间 1200小时。
[0024] 根据本发明设计的方法， 先用 rdkit产生 500构象， 根据结构筛选， 确定 30个最 终构象。 这个过程所需的计算量可以忽略不计。
[0025] 然后将这 30个构象用量子化学的方法结构优化， 所需时间为 2小时 *30=60小时. 力场方法 （此处采用公开的 GAFF2力场） 进行结构优化， 其时间可以忽略。 所以 此步所需 CPU时间为 60小时。
[0026] 对比量子化学的能量 EQM和力场的能量 EMM， 得到 Pearson相关系数 R为 0. 76和方 差 dE为 1. 6 kcal/mol o 这个符合本发明预设的标准。 所以检测力场的表现满足 预定标准， 无需进一步拟合。 整个方法计算量为 60个 CPU时间， 比传统的 1200
小时大幅度降低。
[0027] 为了对比下面步骤的计算量， 本实施例继续下一步。 将这个分子切成 6个分子 碎片。 每个进行从 -180度到 180度每个 15度进行二面角扫描， 这样对于每个分子 碎片需要计算 24个构象， 总共需要计算的构象为 24*6=144个。 每个分子碎片计 算所需的平均 CHJ时间为 0. 16小时， 所以此步所需的总时间为 144*0. 16=23小时 。 本实施例利用这些分子碎片量子化学数据进行拟合得到新的六个柔性二面角 的参数。 将新的参数应用第一步中产生的 30个构象， 与第一步中的量子化学数 据对比， 重新计算能量相关系数和偏差分别为 0. 83 和 1. 3kCal/mol. 可见基于 分子碎片的拟合提升了二面角参数的质量。 到此步， 总共需要的计算时间为 60+ 23=83小时， 远远小于通常所需的 1200小时。
[0028] 为了展示整个流程， 本实施例选取 6个二面角中相对较差的 2个二面角进行整个 分子的二面角扫描， 在不考虑耦合的情况下， 需要的 CPU时间为 2 （二面角个数 ） *24 （每个二面角扫描所需的构象数） *2小时 =96小时， 然后进行这两个二面 角参数的拟合。 以此时所需的总时间为 96+60+23=179小时。
[0029] 根据以上实例的对比， 本发明新的检测和拟合方法所需的计算量远远小于通常 的整个分子的二面角扫描。
Claims
( 1 ) 对于一个大分子， 首先使用 rdkit对于每一个分子产生 500个构 象， 将这些分子用 rdkit自带的 UFF力场进行结构优化， 计算每个结构 的每个柔性二面角的角度； 根据柔性二面角的角度分布， 选取 30个结 构， 涵盖从 -180度到 180度不同的区域， 优先选择能量较低的结构；
(2) 使用量子化学计算软件，将步骤 ( 1 ) 中的结构使用高精度量子 化学方法进行进一步结构优化， 得到相应的能量为 E QM， 同时用需要 检测的力场将这些结构进行优化， 得到对应的能量为 EMM;
(3) 将步骤 (2) 中得到的两组能量以分子为单位进行线性拟合， 得 到两组数据的 Pearson相关系数 R和能量偏差 dE， 如果 R大于第一阈值 而且 d E小于第二阈值， 终止流程， 力场参数在此分子表现较好； 否 则， 进入步骤 (4) ;
(4) 将步骤 (3) 中进入此步骤的分子切成较小的碎片， 每个碎片包 含一个柔性二面角； 对碎片进行通常的二面角扫描， 对比量子化学数 据和力场数据， 对于表现不好的二面角， 利用得到的量子化学数据对
二面角参数进行拟合， 得到新的二面角参数； 所述的对于表现不好， 为两者相关系数 R小于第一阈值或者能量偏差 d E大于第二阈值；
（5） 用步骤 （4） 得到的新的二面角参数， 重复步骤 （2） 中的力场 计算得到新的能量 EMM’； 此时步骤 （2） 中以前计算的量子化学数据 EQM无需重复计算， 然后将 E QM和 EMM’进行线性拟合， 得到两组数据 的相关系数 R’ 和能量偏差 dE’ ； 如果 R’ 大于第一阈值而且 dE’ 小于 第二阈值， 终止流程， 新拟合的二面角参数在此分子表现较好； 否则 ， 进入步骤 （6） ;
（6） 将经过步骤 （4） - 步骤 （5） 中仍然表现不好的参数进行通常 的整个分子的二面角扫描进行拟合过程； 先扫描步骤 （4） 中表现不 好的二面角， 然后进行二面角参数拟合， 将新的拟合参数用最初产生 的结构进行力场的计算， 和步骤 （2） 中已有的量子化学结果 EQM进行 比较， 如果 R大于第一阈值而且 dE小于第二阈值， 终止流程； 否者
， 接着进行其他柔性二面角的扫描， 拟合相关二面角参数。
[权利要求 3] 根据权利要求 2所述的一种测试和拟合力场二面角参数的方法， 其特 征在于， 所述的第一阈值为 0. 7, 第二阈值为 2. 0千卡 /摩尔。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2020/096122 WO2021103491A1 (zh) | 2020-06-15 | 2020-06-15 | 一种测试和拟合力场二面角参数的方法 |
| US17/293,923 US20220310210A1 (en) | 2020-06-15 | 2020-06-15 | Method of testing and fitting the dihedral angle parameters in force field |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2020/096122 WO2021103491A1 (zh) | 2020-06-15 | 2020-06-15 | 一种测试和拟合力场二面角参数的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2021103491A1 true WO2021103491A1 (zh) | 2021-06-03 |
Family
ID=76129143
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/CN2020/096122 WO2021103491A1 (zh) | 2020-06-15 | 2020-06-15 | 一种测试和拟合力场二面角参数的方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US20220310210A1 (zh) |
| WO (1) | WO2021103491A1 (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102779239A (zh) * | 2011-05-09 | 2012-11-14 | 中国科学院研究生院 | 一种用于建立蛋白质体系分子模拟力场的方法 |
| CN106355025A (zh) * | 2016-09-06 | 2017-01-25 | 北京理工大学 | 生命体系中等位基因竞争性反应qm/mm方法 |
| WO2017196963A1 (en) * | 2016-05-10 | 2017-11-16 | Accutar Biotechnology Inc. | Computational method for classifying and predicting protein side chain conformations |
| CN108763852A (zh) * | 2018-05-09 | 2018-11-06 | 深圳晶泰科技有限公司 | 类药有机分子的自动化构象分析方法 |
-
2020
- 2020-06-15 WO PCT/CN2020/096122 patent/WO2021103491A1/zh active Application Filing
- 2020-06-15 US US17/293,923 patent/US20220310210A1/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102779239A (zh) * | 2011-05-09 | 2012-11-14 | 中国科学院研究生院 | 一种用于建立蛋白质体系分子模拟力场的方法 |
| WO2017196963A1 (en) * | 2016-05-10 | 2017-11-16 | Accutar Biotechnology Inc. | Computational method for classifying and predicting protein side chain conformations |
| CN106355025A (zh) * | 2016-09-06 | 2017-01-25 | 北京理工大学 | 生命体系中等位基因竞争性反应qm/mm方法 |
| CN108763852A (zh) * | 2018-05-09 | 2018-11-06 | 深圳晶泰科技有限公司 | 类药有机分子的自动化构象分析方法 |
Non-Patent Citations (2)
| Title |
|---|
| ALEXANDROS ALTIS, PHUONG H. NGUYEN, RAINER HEGGER, GERHARD STOCK: "Dihedral angel principal component analysis of molecular dynamics simulations", THE JOURNAL OF CHEMICAL PHYSICS, vol. 126, no. 24, 28 June 2007 (2007-06-28), pages 244111.1 - 244111.10, XP009528282, ISSN: 0021-9606, DOI: 10.1063/1.2746330 * |
| WANG JUNMEI , ZHAO ZHULIU , YE XUEQI: "Parameterization Procedures in Molecular Mechanisc Calculation", ACTA PHYSICO-CHIMICA SINICA, vol. 11, no. 5, 31 May 1995 (1995-05-31), pages 424 - 428, XP055816578 * |
Also Published As
| Publication number | Publication date |
|---|---|
| US20220310210A1 (en) | 2022-09-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gupta et al. | Quantifying similarity between motifs | |
| Terwilliger | SOLVE and RESOLVE: automated structure solution, density modification and model building | |
| JP4597853B2 (ja) | 堅牢な統計的方法を使用する質量スペクトルの自己較正 | |
| Richterich | Estimation of errors in “raw” DNA sequences: a validation study | |
| US8158436B2 (en) | Methods for characterizing heparin-like glycosaminoglycan mixtures | |
| Knudsen et al. | Sequence alignments and pair hidden Markov models using evolutionary history | |
| WO2023168812A1 (zh) | 一种搜索系统的优化方法、装置、存储介质及计算机设备 | |
| Zhang et al. | A continuous threshold expectile model | |
| WO2021103491A1 (zh) | 一种测试和拟合力场二面角参数的方法 | |
| Empereur-Mot et al. | Automatic multi-objective optimization of coarse-grained lipid force fields using SwarmCG | |
| CN111863140B (zh) | 一种测试和拟合力场二面角参数的方法 | |
| CN112966435B (zh) | 一种桥梁变形实时预测方法 | |
| CN108491690B (zh) | 一种蛋白质组学中肽段的肽段定量效率预测方法 | |
| CN109346128B (zh) | 一种基于残基信息动态选择策略的蛋白质结构预测方法 | |
| Cavattoni et al. | Classgraph: improving metagenomic read classification with overlap graphs | |
| US20140379270A1 (en) | System and method for aligning genome sequence considering mismatch | |
| CN112861308B (zh) | 一种基于时序的弹道曲线平移判断方法 | |
| Chakrabarti et al. | State of the art: refinement of multiple sequence alignments | |
| CN112381456A (zh) | 基于冲突数据融合的装备研制阶段测试性评估及系统 | |
| US20140214338A1 (en) | Polymer side chain analysis method and polymer side chain analysis device | |
| TWI302246B (zh) | ||
| CN117238368B (zh) | 分子遗传标记分型方法和装置、生物个体识别方法和装置 | |
| Cao et al. | Using Variable-Length Aligned Fragment Pairs and an Improved Transition Function for Flexible Protein Structure Alignment | |
| CN116486895A (zh) | 基于残基坐标系的距离信息引导折叠蛋白质结构预测方法 | |
| JP2002350310A (ja) | 応力拡大係数の誤差評価方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 20894560 Country of ref document: EP Kind code of ref document: A1 |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| 32PN | Ep: public notification in the ep bulletin as address of the adressee cannot be established |
Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 25.05.2023) |
|
| 32PN | Ep: public notification in the ep bulletin as address of the adressee cannot be established |
Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 20894560 Country of ref document: EP Kind code of ref document: A1 |