[发明专利]采用k近邻算法筛选人甲状腺素运载蛋白干扰物的方法

权利要求书

1.采用k近邻算法筛选人甲状腺素运载蛋白hTTR干扰物的方法，其特征在于，具体步骤如下：

（1）收集有机化学品干扰效应数据

收集有机化学品的干扰效应数据，所述的干扰效应数据为有机化学品与¹²⁵I-T₄竞争hTTR结合位点的能力即半数竞争效应浓度IC₅₀；

（2）计算描述符

采用基于形态修正的量化描述符来表征可电离基团解离的影响：采用Gaussian 16软件优化有机化学品分子态和离子态的结构，然后基于Gaussian 16的输出文件，直接提取或计算有机化学品分子态和离子态的量化描述符，再根据式（1）计算基于形态修正的量化描述符X_修正

(1)

， (2)

其中，X_M和X_I分别是有机化学品分子态和离子态的描述符值，δ_M和δ_I分别是分子态和离子态的比例分数；并采用Dragon 6.0软件计算官能团和分子碎片描述符以表征有机化学品各种基团对干扰效应的影响；

（3）二元分类模型的构建及表征

采用收集的有机化学品的有无活性的定性数据，根据基于欧几里德距离的kNN算法构建二元分类模型，模型采用经济合作与发展组织关于模型构建与验证的导则进行表征，确定最优模型，所述的最优模型含三个描述符即形态修正的平均分子静电势V_aver-adj、连接到sp3杂化的碳原子的氟原子F-083和连接到sp3杂化或sp2杂化的碳原子的氢原子H-047，邻近数k为3，所述的二元分类模型的应用域为欧几里德距离小于0.928；

（4）定量预测模型的构建及表征

选取采用相同测试方法和测试条件得到的定量数据，根据基于欧几里德距离的kNN算法构建定量预测模型，建模时用取对数的相对效应势RP表征有机化学品与¹²⁵I-T₄竞争hTTR结合位点的能力，RP定义为：

(3)

其中，IC₅₀(T₄)和IC₅₀(有机化学品)分别代表甲状腺素和有机化学品的IC₅₀；确定最优模型，所述的最优模型含四个描述符：Sp2杂化的取代苯碳原子数nCb-、酚羟基数nArOH、分子内氢键数nHBonds和形态修正的平均分散度V_adj，邻近数k为3；所述的定量预测模型的应用域为欧几里德距离小于1.11；

(5) 人甲状腺素运载蛋白干扰物的筛选

①计算分类模型所需描述符，即形态修正的平均分子静电势V_aver-adj、连接到sp3杂化的碳原子的氟原子F-083、连接到sp3杂化或sp2杂化的碳原子的氢原子H-047，评估目标有机化学品是否在二元分类模型应用域内；

若目标有机化学品在二元分类模型应用域范围内，则根据二元分类模型计算目标有机化学品是否具有hTTR干扰活性；若目标有机化学品无活性，则不需要进一步评估；若目标有机化学品有活性，则依据定量预测模型预测其干扰效应数值大小；若目标有机化学品不在二元分类模型的应用域范围内，则不能用二元分类模型进行预测；

②对于有活性的目标有机化学品，根据定量预测模型的要求计算所需描述符，即Sp2杂化的取代苯碳原子数nCb-、酚羟基数nArOH、分子内氢键数nHBonds和形态修正的平均分散度V_adj，评估其是否在定量预测模型的应用域范围内；

若目标有机化学品在定量预测模型的应用域范围内，则依据选择的定量预测模型计算目标有机化学品对hTTR的logRP值；若目标有机化学品不在定量预测模型的应用域范围内，则不能用定量预测模型进行预测；

③根据定量预测模型预测的logRP值判断目标有机化学品是否具有干扰hTTR转运甲状腺素的能力：

若有机化学品logRP 0，则表明目标有机化学品与hTTR结合能力强于甲状腺素；

若有机化学品logRP= 0，则表明目标有机化学品与hTTR结合能力与甲状腺素相近；

若有机化学品logRP0，则表明目标有机化学品与hTTR结合能力弱于甲状腺素。