类风湿性关节炎（Rheumatoid Arthritis, RA）是一种让患者关节逐渐畸形的自身免疫疾病，它主要困扰女性和老年人，导致持续性疼痛、功能障碍甚至过早死亡。尽管现有药物如来氟米特（Leflunomide）、甲氨蝶呤（Methotrexate）能够缓解症状，但往往伴随严重的神经系统、皮肤和胃肠道毒副作用，更棘手的是，它们只能控制病情进展，无法实现根治。面对这一临床困境，科学家们将目光投向了一个关键靶点——人源二氢乳清酸脱氢酶（Human Dihydroorotate Dehydrogenase, DHODH）。这种位于线粒体内膜的黄素依赖酶，是嘧啶核苷酸从头合成（de novo pyrimidine biosynthesis）途径的限速酶，而快速增殖的免疫细胞正依赖此途径获取遗传物质合成原料。抑制DHODH就能切断异常活化淋巴细胞的"补给线"，从而有效控制RA疾病进程。然而，现有DHODH抑制剂如brequinar因药代动力学缺陷止步临床，这促使摩洛哥研究团队在《Scientific African》发表了一项创新研究，通过计算机模拟技术设计新一代高效低毒的DHODH抑制剂。
研究团队运用多层次计算生物学方法展开攻关：首先基于35个丙烯酰胺衍生物的活性数据建立3D-QSAR模型，通过分子场分析揭示影响抑制活性的关键化学特征；接着利用分子对接技术将设计的22个新分子与DHODH蛋白（PDB: 1D3G）结合，筛选出结合模式最优的候选分子；进一步采用MM-GBSA方法精确计算结合自由能，并通过100纳秒分子动力学模拟验证复合物稳定性；最后利用瑞士ADME和pkCSM在线平台全面评估候选物的药代动力学特性和毒性风险。
分子对接与结合自由能分析
通过将候选分子与DHODH活性位点对接，发现P4、P11和P18的结合得分（-9.801至-10.216 kcal/mol）显著优于先导化合物M35（-9.744 kcal/mol）。MM-GBSA计算显示这些分子与蛋白的结合自由能（ΔGbind 介于-48.19至-51.41 kcal/mol）远低于M35（-44.57 kcal/mol），其中范德华相互作用和疏水效应对结合稳定性贡献最大。
ADMET性质预测
药代动力学评估表明，所有候选分子均符合利宾斯基五规则（Lipinski's Rule of Five），具有良好口服生物利用度。特别是P4、P11和P18表现出适宜的水溶性（logS约-2.9）、肠道吸收率（>30%）和组织分布能力（VDss>0.15），且无血液脑屏障穿透性，预示中枢神经副作用风险较低。AMES毒性测试显示所有分子无致突变性，安全性特征与参考抑制剂brequinar相当。
分子动力学模拟验证
100纳秒动力学轨迹分析显示，P4-DHODH等复合物的均方根偏差（RMSD）稳定在1.5-2.1Å之间，关键氨基酸Gln147和Arg136通过氢键与配体保持持续相互作用。配体自身构象参数（分子表面积、回转半径等）波动幅度微小，证明复合物在生理环境下能维持稳定结合构象。
合成可行性分析
研究人员还设计了以邻苯二胺为起始物的合成路线，通过环化、氧化、Wittig烯化等六步反应可高效制备目标分子。合成可及性评分（2.56-2.75/10）证实实验可行性，核磁共振（NMR）和红外光谱（IR）特征峰为后续结构验证提供明确指引。
本研究通过计算生物学手段成功设计出P4、P11和P18三个新型丙烯酰胺类DHODH抑制剂。这些分子不仅表现出优于先导化合物的结合亲和力和特异性，更具备良好的类药性与安全性特征。分子动力学模拟从原子层面揭示了配体-受体复合物的稳定机制，特别是Gln147和Arg136残基在结合过程中的关键作用。该研究不仅为RA治疗提供了极具潜力的候选药物，更展示了计算机辅助药物设计在优化药物效能与安全性方面的强大能力，为自身免疫疾病药物开发开辟了新途径。
打赏
