卤代维拉帕米衍生物作为预测性自噬调节剂的理性设计：网络药理学与分子动力学揭示其对AKT/mTOR和MAPK级联的靶向潜力

《Letters in Drug Design & Discovery》：Rational design of halogenated predicted autophagy modulators: network pharmacology and molecular dynamics reveal verapamil derivatives targeting AKT/mTOR and MAPK cascades

　　

维拉帕米类似物的生物电子等排设计产生了新型自噬调节剂。VA7类似物通过卤化显示出对MTOR（mechanistic target of rapamycin，机制性雷帕霉素靶蛋白）和MAPK8（mitogen-activated protein kinase 8，丝裂原活化蛋白激酶8）的更优结合。100 ns分子动力学（molecular dynamics，MD）模拟将MAPK3（mitogen-activated protein kinase 3，丝裂原活化蛋白激酶3）过滤为静态对接假阳性。VA7通过预测的双命中相互作用谱发挥作用：抑制MTOR并稳定MAPK8。ADMET（absorption, distribution, metabolism, excretion, and toxicity，吸收、分布、代谢、排泄和毒性）分析证实VA7的药物相似性且缺乏遗传毒性风险。

自噬（autophagy）是高度有序、进化保守的分解代谢过程，对真核细胞稳态维持至关重要。其分子机器受mTORC₁（mechanistic target of rapamycin complex 1，机制性雷帕霉素靶蛋白复合物1）和AMPK（AMP-activated protein kinase，AMP活化蛋白激酶）等复杂信号网络调控，涉及ULK1、PI3K、膜延伸、自噬溶酶体降解和营养回收。该通路失调与心血管疾病、神经退行性疾病和癌症相关，故药理学调节自噬是药物发现优先策略。二甲双胍、地奥司明和维拉帕米等临床化合物可干预自噬级联不同节点，但其临床应用受亚优药代动力学和特定胞内靶标选择性不足限制，因此需要理性结构优化。计算机药理学可模拟分子相互作用并预测药代性质，减少研究时间和成本，符合3Rs原则。该研究旨在提出并评估二甲双胍、地奥司明和维拉帕米的新型衍生物，通过生物电子等排和分子动力学模拟获得对自噬主调控因子选择性和结合稳定性增强的分子。

关键技术方法：研究人员通过文献和数据库筛选候选药物，从PubChem获取SMILES；用MolOpt进行生物电子等排生成类似物；用SwissADME评估类药性；用ADMETlab 3.0评估毒性、合成可及性和ADMET；从PubChem与CTD获取靶基因并以InteractiVenn取交集，用STRING构建PPI网络；从RCSB PDB和AlphaFold获取蛋白结构，经UniProt确认定位；用Avogadro、UCSF Chimera、AutoDock Tools和AutoDock Vina对接；用GROMACS 2025.2、AMBER99SB-ILDN、ACPYPE和TIP3P进行100 ns分子动力学与MM-PBSA。无人类或动物样本队列，数据均来自公开数据库。

研究结果：3.1 化学空间扩展与理性先导优化。研究人员用MolOpt对三个FDA批准支架进行数据挖掘转化，生成80个类似物（排除母药）。3.1.1 二甲双胍衍生类似物：生成12个，引入末端苯基、吡啶和酰胺基团，以提高亲脂性和分子复杂性。3.1.2 地奥司明衍生类似物：生成18个，策略包括B环羟基化、多甲氧基化、氟原子和侧链羧酸基团。3.1.3 维拉帕米衍生类似物：生成50个，包括腈到酰胺生物电子等排、三氟甲氧基（-OCF₃）、侧链截短和羟基引入。3.2 理化筛选与类药性评价：SwissADME五规则筛选中，二甲双胍衍生物因分子量低而常失败Muegge过滤；地奥司明B环羟基化类似物零违反，羧酸衍生物因TPSA>131.6 ?²失败Egan和Muegge；维拉帕米衍生物最平衡，三氟甲氧基和酰胺取代通过Lipinski、Ghose和Egan，部分因>10可旋转键触发Veber，截短侧链和羟基化衍生物生物利用度评分0.55。3.3 毒理分析和合成可及性：ADMETlab 3.0显示多数二甲双胍和地奥司明类似物在致癌性、致突变性或遗传毒性上为“Suspect”；维拉帕米衍生物VA7、VA9、VA12、VA14、VA16、VA17和VA44通过全部毒理过滤且合成可及性容易。3.4 药代动力学与安全性评估：候选物Caco-2渗透性为-4.805至-5.459；VA7和VA12血浆蛋白结合高（85.9%和88.6%），VA14和VA17低（17.7%）；VA7和VA12血脑屏障渗透高，短链类似物低或无；均为CYP1A2、CYP2C19、CYP2D6和CYP3A4底物，T_1/2为1.713–2.542 h；AMES和遗传毒性阴性，致癌风险安全，但VA7的hERG阻断概率为0.89。3.4.1 吸收与分布：候选物具中等至高肠道吸收潜力，PPB和BBB差异显著。3.4.2 代谢与排泄：候选物经主要CYP亚型代谢，半衰期适合标准给药。3.4.3 毒理学评价（分子安全性）：优先化合物无致突变和遗传毒性，VA16致癌风险最低0.082，人肝毒性中等，VA7 hERG风险需全细胞膜片钳验证。3.5 靶基因识别与数据整合：PubChem和CTD整合得到88个共同基因（41.5%），PubChem独有65个，CTD独有59个。3.6 PPI网络与枢纽识别：STRING网络识别MTOR、AKT1、MAPK1、MAPK3、MAPK8和BCL2为高连接枢纽，选作对接靶点。3.7 分子对接：AutoDock Vina显示VA7对MTOR为-7.0 kcal/mol（维拉帕米-6.2），对MAPK8为-7.5 kcal/mol（维拉帕米-6.8），对MAPK3达-8.0 kcal/mol；VA12对MAPK1为-6.9 kcal/mol；VA7与MTOR形成ARG2224氢键、ALA912和VAL915卤键及ALA2210 Pi-Sigma接触，与MAPK8形成ARG69 Pi-Cation、GLU73氢键和ASN114卤键。AKT1-维拉帕米存在ASN53不利受体-受体相互作用，新类似物避免该冲突。3.8 分子动力学模拟：100 ns轨迹显示MTOR-VA7和MAPK8-VA7稳定。3.8.1 构象稳定性（RMSD）与结构紧凑性（R_g）：MTOR-VA7 RMSD低于1 nm，MAPK8-VA7近零波动；MAPK3体系偏差5–20 nm，R_g峰值达8.0 nm，MAPK8约2.2 nm。3.8.2 残基柔性与结合位点组织（RMSF）：MTOR和MAPK8波动低，MAPK3-维拉帕米波动超过4.5 nm，表明VA7可抑制不利构象。3.8.3 相互作用持续性：MTOR多数轨迹维持至少一个氢键，MAPK8在一到两个氢键间交替，MAPK3氢键稀疏短暂。3.8.4 每残基结合自由能分解：MM-PBSA显示MTOR-VA7的1583TYR贡献-16.5 kJ/mol，1582SER -7.3，1580GLY -5.9；维拉帕米为829GLN -15.4和926SER -10.1。MAPK8-VA7的290LYS贡献-9.8，126GLU -6.9；维拉帕米96LYS仅-4.3。MAPK3贡献极低。总结合自由能中MTOR-VA7为-189.453 kJ/mol，MTOR-维拉帕米-82.331，MAPK8-VA7 -83.299，MAPK8-维拉帕米-69.887。

讨论部分总结：原文指出，研究核心目标是评估维拉帕米理性结构修饰能否增强其对自噬机器的多靶点效力。计算流程优先VA7，其预测热力学和动力学优于母药。MD轨迹提示VA7与MTOR和MAPK8形成稳定复合物，达到维拉帕米无法复制的能量锚定。MTOR是自噬主负调控因子，抑制MTOR模拟营养剥夺并诱导自噬；VA7绕过AKT1空间冲突，可能形成MTOR抑制和MAPK8稳定的双命中作用谱。MAPK8为应激响应自噬正调控因子，VA7与MAPK8近刚性稳定，原文假设VA7可能调节自噬通量，但需生化验证。与传统认为维拉帕米为mTOR非依赖性自噬诱导剂不同，卤化后获得直接靶向mTOR/MAPK催化区的新相互作用。ADMET支持VA7转化可行性，但hERG风险需实验处理。局限包括纯计算机研究、无合成和生化及细胞自噬通量测定；单条100 ns轨迹；无阳性对照；AlphaFold模型无法redocking；卤素以点电荷建模；ADMET为理论预测；MM-PBSA未解析溶剂惩罚，总结合自由能反映气相分子力学相互作用。

研究结论部分翻译：5.1 背景与依据：药理学调节自噬是细胞和蛋白毒性失调疾病的关键治疗策略。但传统调节剂常有选择性不足和不利药代动力学。本研究应用多阶段计算机药物化学和网络药理学框架，设计并计算评估二甲双胍、地奥司明和维拉帕米支架的新型生物电子等排衍生物。5.2 关键计算结果：基于结构的流程优先卤化类似物VA7，其预测热力学优于母药维拉帕米。100 ns分子动力学轨迹显示，VA7与核心自噬调节因子形成稳定、低波动复合物，并由MTOR（1583TYR；-16.5 kJ/mol）和MAPK8（290LYS；-9.8 kJ/mol）的强能量锚定驱动。核心假设是VA7可能作为自噬通量的推定性调节剂。通过竞争性结合MTOR催化区并稳定应激响应激酶MAPK8构象，该化合物呈现与囊泡成核相容的理论谱。预测ADMET模型提示其具有利于胞内靶向的分布体积，无遗传毒性，但显示特定hERG阻断概率。5.3 未来展望与实验验证：尽管结果有希望，该研究完全基于预测性生物物理模型。所用计算机流程成功过滤不稳定候选，支持3Rs原则，减少动物依赖，但这些结构假设严格需要实验验证。因此，VA7类似物未来任何转化推进均依赖体外生化试验，尤其是LC3-I向LC3-II的免疫印迹以确认功能性自噬通量诱导，以及全细胞膜片钳试验以严格排除预测的心脏毒性。