《In Silico Pharmacology》：In silico identification of?novel PTP1B inhibitors by the investigation of molecular mechanism, QSAR, and DFT studies
编辑推荐：
蛋白酪氨酸磷酸酶1B（PTP1B）是一个公认且具有前景的药物靶点，参与胰岛素和瘦素信号的负调控，治疗2型糖尿病的新型抗糖尿病分子直接依赖于它。研究人员利用计算机模拟技术，包括虚拟筛选、分子对接、ADMET分析、毒性预测、VEGA-QSAR分析、分子动力学（MD
蛋白酪氨酸磷酸酶1B（PTP1B）是一个公认且具有前景的药物靶点，参与胰岛素和瘦素信号的负调控，治疗2型糖尿病的新型抗糖尿病分子直接依赖于它。研究人员利用计算机模拟技术，包括虚拟筛选、分子对接、ADMET分析、毒性预测、VEGA-QSAR分析、分子动力学（MD）模拟和密度泛函理论（DFT）计算，发现并表征了一系列PTP1B抑制剂。从1000个抗糖尿病化合物的初始库中筛选，鉴定出PubChem CID 44560696为最成功的命中化合物。随后通过基于化学类似性的优化，鉴定出最有效的抑制剂为CID 44560744，其结合自由能为?9.13 kcal/mol，并能与关键催化残基形成多个稳定氢键和疏水相互作用。DFT优化几何结构和对接揭示，该先导化合物显示出高稳定性和强PTP1B抑制活性，与熊果酸（结合能为?6.34 kcal mol?1）相比，具有更优的结合亲和力。MD模拟显示，在200 ns模拟过程中，PTP1B-配体复合物结构稳定且紧凑，在RMSD、RMSF、Rg、SASA和氢键方面均表现出理想特征。DFT和前线分子轨道（FMO）分析表明，该化合物因其低HOMO-LUMO能隙、高柔软度和中等亲电性，具有极高的化学反应活性和电子转移能力。ADMET和毒性评估（SwissADME、ProTox-II和VEGA-QSAR）的结果表明，其具有满意的溶解度、无致癌性、无致突变性以及高药理学安全性。总体而言，计算方法的协同作用使CID 44560744成为开发和优化用于T2DM治疗的选择性PTP1B抑制剂的先导骨架。
**研究背景与问题**
2型糖尿病（Type 2 Diabetes Mellitus, T2DM）是最常见的异质性代谢紊乱疾病，以慢性高血糖和碳水化合物、脂质、蛋白质代谢紊乱为特征。蛋白酪氨酸磷酸酶1B（Protein Tyrosine Phosphatase 1B, PTP1B）是胰岛素和瘦素信号通路的关键负调控因子，被认为是治疗T2DM和肥胖的潜在药物靶点。然而，PTP1B抑制剂的开发面临重大挑战：其催化位点与其他蛋白酪氨酸磷酸酶（如T细胞蛋白酪氨酸磷酸酶TCPTP，同源性约80%）高度保守，导致选择性难以实现；此外，还需平衡细胞通透性和口服生物利用度。尽管过去二十年进行了大量研究，但仅有少数化合物（如trodusquemine、ertiprotafib、JTT-551）进入临床试验，且多数因副作用或疗效不足而终止。因此，迫切需要利用计算机辅助药物设计（Computer-Aided Drug Design, CADD）方法快速、低成本地发现新型PTP1B抑制剂。本研究旨在通过整合虚拟筛选、分子对接、ADMET分析、毒性预测、VEGA-QSAR、分子动力学（Molecular Dynamics, MD）模拟和密度泛函理论（Density Functional Theory, DFT）等多种计算机模拟技术，系统筛选并优化PTP1B抑制剂，为T2DM治疗提供有前景的先导化合物。该论文发表在《In Silico Pharmacology》。
**主要技术方法**
研究人员采用以下关键计算技术：从PubChem和ChEMBL数据库获取1000个抗糖尿病化合物作为初始库；使用AutoDock Vina（通过PyRx界面）进行基于结构的虚拟筛选和分子对接；利用SwissADME、ProTox-II和VEGA-QSAR平台进行ADMET（吸收、分布、代谢、排泄、毒性）及毒性预测；采用GROMACS软件，结合AMBER力场，进行200 ns的分子动力学模拟，分析均方根偏差（RMSD）、均方根波动（RMSF）、回转半径（Rg）、溶剂可及表面积（SASA）和氢键；使用Materials Studio 2020的DMol3模块，在M06-2X/DNP水平上进行DFT计算，包括前线分子轨道（FMO）分析和分子静电势（MEP）映射。
**研究结果**
- **虚拟筛选**：从1000个化合物中筛选，初始命中化合物为PubChem CID 44560696（结合能-8.34 kcal/mol）。通过基于化学类似性的优化，鉴定出类似物CID 44560744，其结合能显著提升至-9.13 kcal/mol，RMSD为2.02 ?，表明形成更稳定的复合物。
- **DFT：前线分子轨道（FMOs）和分子静电势（MEP）**：DFT计算显示，CID 44560744的HOMO-LUMO能隙为7.168 eV，表明高稳定性；化学硬度（η）为3.584 eV，柔软度（S）为0.139 eV
?1，亲电性指数（ω）为2.13 eV，揭示了该化合物具有受控的化学反应活性和电子转移能力。MEP图显示富电子（亲核）和缺电子（亲电）区域，为反应位点预测提供依据。
- **态密度（DOS）和部分态密度（PDOS）**：DOS图显示在费米能级附近有尖锐峰，表明活跃分子轨道参与电子跃迁；PDOS分析表明p轨道对费米能级贡献最大，主导分子化学反应活性。
- **对接研究**：对接显示CID 44560744与PTP1B催化口袋中的关键残基（LYS116、ASP181、GLY183、ARG221、GLN262、GLN266等）形成常规氢键，并与LYS120、MET258等产生疏水及π-硫相互作用。该化合物同时与A、B、C、D四个亚口袋残基结合，表现出多口袋结合模式，有利于选择性抑制。
- **熊果酸（标准化合物）的对接研究**：标准抑制剂熊果酸与PTP1B的结合能为-6.34 kcal/mol，低于先导化合物，证实CID 44560744具有更优的结合潜力。
- **天然配体对接验证**：将共结晶天然配体重新对接，其预测构象与晶体结构高度重叠，验证了对接方法的可靠性。
- **MD模拟**：200 ns的MD模拟中，PTP1B-CID 44560744复合物的RMSD、RMSF、Rg、SASA和氢键分析均显示结构稳定且紧凑。RMSD在初始波动后趋于稳定，RMSF表明配体结合对残基波动影响小，Rg和SASA值显示复合物比游离蛋白更紧凑。氢键分析显示配体与溶剂持续形成1-2个氢键，支持良好水溶性。配体RMSD稳定在4.5-5 ?，表明构象稳定。主成分分析（PCA）显示配体结合后PTP1B构象自由度降低。MM/GBSA计算得到总结合自由能为-13.78 ± 0.18 kcal/mol，范德华力（-35.25 kcal/mol）为主要驱动力。氢键占有率分析显示与ASP240的氢键短暂（0.79%），芳香相互作用图显示多个芳香接触。
- **药物毒性**：ProTox-II预测CID 44560744无肝毒性、神经毒性、肾毒性、呼吸毒性、致癌性、免疫毒性、致突变性和细胞毒性，仅临床毒性预测为“活性”（概率0.51），但属边缘预测。
- **VEGA QSAR和ADMET分析**：VEGA QSAR预测无皮肤刺激、无致突变性、无P-糖蛋白抑制，血浆蛋白结合低（logK=0.0164），消除半衰期1.158 h。SwissADME显示该化合物水溶性好，但存在多项Lipinski规则违反（如分子量大、极性高），胃肠道吸收低，提示需后续结构优化。
**讨论与结论**
讨论部分强调，通过虚拟筛选、DFT优化、分子对接、MD模拟等计算方法的协同作用，成功鉴定出CID 44560744为高效PTP1B抑制剂。该化合物具有高结合亲和力、动态稳定性、良好电子结构特征和低毒性，但存在Lipinski规则违反和低胃肠道吸收的问题，需通过结构修饰或剂型改进解决。结论部分翻译如下：本研究采用系统、高效的计算机模拟研究流程，全面筛选、优化并阐明蛋白酪氨酸磷酸酶1B（PTP1B）抑制剂的作用机制，最终鉴定出有前景的候选化合物CID 44560744。从大规模虚拟筛选开始，对来自PubChem和ChEMBL的1000个抗糖尿病相关分子进行结构基础筛选，鉴定出初始高亲和力分子CID 44560696。随后以CID 44560744为模板进行结构类似物搜索，并进一步进行分子对接筛选。CID 44560744展现出优异的结合能（-9.13 kcal/mol）和稳定的结合构象。使用DFT优化结构对先导化合物进行分子重新对接。它与关键催化残基（如Asp181、Arg221和Gln262等）形成多个氢键和疏水相互作用，表现出强活性位点亲和力和结构适应性。与参考抑制剂熊果酸的对比对接显示，其结合亲和力（-6.34 kcal mol
?1）低于先导化合物，证实了筛选化合物对PTP1B的更优结合潜力。对接回测实验显示预测构象与晶体结构高度重叠，证明了当前对接过程的可靠性。为进一步确认其结合稳定性，随后的200 ns分子动力学模拟显示，PTP1B-CID 44560744复合物在RMSD、RMSF、Rg和SASA指标上均表现出稳定趋势。复合物结构更紧凑，证明其在动态环境中的结合耐受性，有利于其抑制效果。氢键时序分析也表明，该分子能够长期维持关键相互作用，充分支持其稳定结合能力。在电子结构层面，基于DFT和FMO的计算进一步揭示了分子电子行为。结果表明，该化合物具有极高的HOMO-LUMO能隙（7.168 eV），表明高电子迁移率和稳定性，有利于与酶催化口袋形成稳定的静电和极性相互作用。同时，其低柔软度和中等电荷转移能力反映了其良好的电子调控特性和潜在稳定抑制效果。在药代动力学和毒理学方面，使用SwissADME、ProTox-II和VEGA-QSAR的联合评估表明，CID 44560744具有良好的水溶性、无致癌性、无显著肝毒性或遗传毒性，总体安全性高。尽管其与Lipinski规则有所偏离且胃肠道吸收低，但这种高极性、高分子量的结构在PTP1B抑制剂中较为常见，可通过后续结构优化、修饰或药物递送策略进一步改善。总之，本研究成功鉴定出CID 44560744作为潜在高效PTP1B抑制剂，在结合亲和力、动力学稳定性、电子结构特征和毒理学安全性方面均表现出优势。这些结果为其作为治疗2型糖尿病先导化合物的潜力奠定了坚实基础。本研究也证明了计算机模拟多技术融合策略在发现新型PTP1B抑制剂中的重要性。未来应结合体外酶学实验、细胞模型和体内药效评估，进一步验证其抑制活性、特异性和药代动力学性质，并开展必要的结构修饰，以推动其向临床候选药物迈进。
