《Chemical Biology & Drug Design》：Cruzain Inhibitors for Chagas Disease: Anticorrelated Optimisation Landscapes and the Multiparametric Path to Clinical Candidates
三十余年来，克氏锥虫（Trypanosoma cruzi）主要半胱氨酸蛋白酶cruzain（CZ）作为恰加斯病（Chagas disease, CD）的验证药物靶点，已积累了215个IC50 ≤ 1 μM的抑制剂、超过30个晶体结构，却未能产生任何临床候选药物。该综述认为，这一转化悖论并非源于靶点无效，而是源于cruzain抑制剂优化方式与细胞内疗效所需条件之间的持续性战略错配。在肽基乙烯基砜（peptidyl vinyl sulfones）、非肽骨架（non-peptidic scaffolds）及硫代半卡巴腙（thiosemicarbazones）中，酶促效力被系统性优先于细胞内暴露的耦合决定因素：渗透性（permeability）、代谢稳定性（metabolic stability）和组织蛋白酶（cathepsin, CTS）选择性。每类抑制剂因不同的近因失败，而这些失败汇聚于同一种优化偏差。计算方法（分子对接（docking）、分子动力学（molecular dynamics, MD）、自由能微扰（free-energy perturbation, FEP）和机器学习定量构效关系（machine-learning quantitative structure–activity relationship, ML-QSAR））的应用日趋精密，但均针对以亲和力为中心的目标，未能充分代表暴露约束。作为其核心贡献，本研究提出了一种多参数优化（multiparametric optimisation, MPO）框架、临时可行性基准以及可行性包络（feasibility-envelope）概念，将cruzain抑制剂的发现重新定义为约束满足问题而非效力最大化问题。该框架概述了下一代cruzain先导化合物可能需要满足的临时定量标准，并解释了结构可及性良好的细胞内半胱氨酸蛋白酶靶点为何仍未能实现临床转化。
1 引言：Cruzain靶向药物发现中的三十年悖论
恰加斯病（CD）药物发现中，极少有靶点能像cruzain（CZ）——克氏锥虫（T. cruzi）的主要半胱氨酸蛋白酶——如此清晰地展示生物学合理性与转化成功之间的差距。三十年间，CZ积累了高可信度临床前靶点的标志性特征：无直接人类直系同源物、遗传必需性跨所有寄生虫生活史阶段、超过30个高分辨率晶体结构，以及K777的药理学概念验证——该亚纳摩尔抑制剂（K
i=0.2 nM）在急性感染小鼠模型中实现寄生虫学治愈。研究人员系统整理了ChEMBL靶点CHEMBL3563数据：在pChEMBL ≥ 6（IC
50 ≤ 1 μM）阈值下，215个独特抑制剂通过Molecule ChEMBL ID去重后获得资格；将阈值提高至pChEMBL ≥ 7（IC
50 ≤ 100 nM）则产生132个化合物。然而，无一进入临床开发。该综述认为，障碍不在于靶点合法性，而在于反复出现的优化错配：三十年以酶促效力为优先的单参数优化，系统性忽视了渗透性、代谢稳定性和CTS选择性这些耦合约束条件——这些参数最终决定细胞内药物暴露和临床可行性。
CZ抑制剂景观的回顾性分析揭示了三种化学不同骨架类别中的汇聚性失败模式。肽基乙烯基砜最大化靶点结合，但因暴露、安全性和选择性缺陷而未能超越临床前模型。非肽骨架部分满足Lipinski类药性标准，但酶促IC
50与细胞EC
50之间存在显著脱节，跨骨架类别约为一到两个数量级。苯并咪唑衍生物的酶促IC
50为5–50 μM，而细胞活性仍较弱（>100 μM；酶促-细胞脱节约6–20倍）。类似地，硫代半卡巴腙的酶促IC
50低至380 nM，但仅在低微摩尔范围保留杀锥虫活性，表明约一个数量级的酶-细胞脱节——Du等将其归因于快速代谢淬灭而非渗透性失败。这些差异反映了寄生虫纳虫空泡（parasitophorous vacuole）的穿透不足，以及主动外排介导的清除。值得注意的是，优化的环状酰亚胺实现了有利的酶促-细胞转化，表明该渗透性障碍可通过理性设计克服。源自天然产物的化合物（IC
50 = 15–100 μM）贡献结构多样性，但缺乏先导优化所需内在效力。所有三类的统一解释是顺序单参数优化（sequential single-parameter optimisation, SPO）：每个骨架被驱动朝向单一可测量性质，而细胞内药物暴露的相互关联需求被系统性降格。结果是属性景观中，沿一个参数的改善常以牺牲另一个为代价——迭代单参数药物化学已被证明难以驾驭的约束。
2 Cruzain作为药物靶点：卓越的可及性与固有限制
CZ在克氏锥虫主要生活史阶段（上鞭毛体、锥鞭毛体和细胞内无鞭毛体）均功能活跃，但治疗决定性阶段为细胞内无鞭毛体，其驱动慢性感染中的疾病持续。这一生物学现实对CZ抑制剂设计施加核心药理学约束：化合物必须达到并维持抑制浓度，不仅在血浆中，而且在寄居于宿主细胞内的寄生虫的酸性溶酶体区室中——这一递送挑战远比细胞外靶点复杂。
木瓜蛋白酶样（papain-like, C1）折叠组织了一个双叶活性位点，包含α-螺旋L结构域和β-折叠富集的R结构域，围绕保守的Cys25-His162-Asn182催化三联体。这一亲核几何构型在弱酸性条件下被激活。超过30个高分辨率CZ晶体结构共同在亚2.0 ?分辨率定义了这一结合结构。S2亚位点异常深、疏水且耐受 bulky P2取代基，是结构导向设计的主要药效团锚点。CZ的遗传必需性得到趋同性证据支持：cruzain基因阵列的遗传破坏和功能减弱导致严重上鞭毛体生长缺陷、分化受损和宿主细胞侵袭减少。化学验证进一步佐证这一必需性：肽模拟物乙烯基砜K777不可逆灭活CZ，具有亚纳摩尔酶促效力，并在急性T. cruzi感染小鼠模型中经口服给药实现稳健寄生虫学治愈。
CZ定位随生活史阶段变化，具有直接药理学后果。该酶在约pH 5.5（溶酶体腔）至7.4（宿主细胞细胞质界面）的梯度中运作，要求抑制剂在超过2个pH单位范围内保持结合能力。此外，化合物必须依次穿越宿主细胞膜和寄生虫纳虫空泡膜——累积渗透性需求远非标准Caco-2或PAMPA（平行人工膜渗透性实验）测定所能捕捉。这一pH分配约束引入设计悖论：在酸性溶酶体中积累的理想弱碱性化合物，可能因催化三联体几何构型的质子化依赖性变化而同时遭受结合亲和力降低。
除药代动力学约束外，进化保守性带来选择性挑战。CZ与人类CTSs仅共享部分总体序列同一性，但其活性位点区域保守性显著更高，与人CTSs的活性位点残基报告有>70%序列同一性。催化机制在结构和功能层面相应高度保守。Cys25-His162-Asn182三联体和Gln19介导的氧阴离子孔在CA族半胱氨酸蛋白酶中结构可叠合。虽此保守性确认CZ为机制上可成药靶点，但也施加重约束：S2和S3亚位点中的效力驱动优化常导致脱靶CTS活性。通过非保守残基的相互作用工程化选择性必须利用S2和S3亚位点——恰为CZ同工型家族间变异性较大的区域。
3 Cruzain抑制剂
3.1 为何单独效力失败：K777范式
K777是最高级的CZ抑制剂，其未能超越临床前开发揭示了更广泛的转化病理。K777通过Michael加成不可逆灭活CZ，形成抗水解裂解的β-硫醚砜加合物，酶促效力明确：K
i = 0.2 nM。该效力转化为药理学概念验证：口服K777在急性T. cruzi感染小鼠模型中产生稳健寄生虫学治愈。
尽管疗效令人信服，K777未能推进为临床候选药物。失败源于多重骨架固有限质的汇聚。肽模拟骨架严重限制且非线性的口服暴露，部分由于转运蛋白介导的外排；同时，乙烯基砜弹头作为CYP3A4的机制性抑制剂，引入高风险药物-药物相互作用（drug-drug interaction）可能性。重复剂量毒理学研究进一步揭示高级别物种中的剂量限制性肝脏和胃肠道不良反应，有效缩小治疗窗。
这些固有限质并非独立或顺序障碍，而是同一分子结构的结构纠缠特征。对任何单一局限的缓解尝试均保留其余局限，例证了顺序单参数优化（SPO）的根本约束。从K777到WRR-483、WRR-662再到WRR-669的轨迹——每次修饰针对单一可测量参数，但均未逃脱肽骨架施加的根本暴露约束——是整个乙烯基砜类别的缩影。
此外，K777范式的另一维度由酶-有机体转化间隙揭示。K777本身保留对细胞内无鞭毛体的低微摩尔活性，但更广泛系列的系统分析记录从重组CZ测定到全寄生虫感染测定的效力侵蚀达一至三个数量级。这一脱节反映寄生虫纳虫空泡穿透、pH依赖性细胞内分配和主动外排对寄生虫复制部位有效药物浓度的综合影响。因此，分离酶测得的效力是生物体水平疗效的不良预测因子——这一事实在K777系列中已 evident，并在所有后续肽基乙烯基砜类似物中被放大。
3.2 渗透性墙：当类药性未能转化为细胞内暴露
K777僵局之后，领域转向非肽、类药性骨架（腙、硫代半卡巴腙、三嗪和苯并咪唑），旨在满足肽架构固有违反的分子量、极性表面积和化学稳定性标准。
然而，非肽时代暴露了第二个同等限制性的转化障碍。跨化学不同非肽系列，对CZ的酶促效力常未能转化为细胞或感染相关测定中的持续抗寄生虫活性。这一酶-细胞脱节已多次记录，其一致性指向系统性暴露问题而非随机测定伪影。
有效抑制剂因此必须依次穿越宿主细胞膜、寄生虫纳虫空泡膜和寄生虫细胞膜，在感染宿主细胞中以足以靶点结合的浓度积累。这一多层障碍已被认识为抗寄生虫效力的主要决定因素，且常规渗透性测定或Lipinski启发式方法难以捕捉。
非肽CZ抑制剂频繁说明这一局限。为生理pH下膜渗透性优化的弱碱性或中性杂环，可能在酸性或区室化细胞内微环境中未能保持充分质子化，限制空泡保留和有效靶点结合。反之，为利用pH依赖性积累而设计的化合物在遭遇近中性胞质条件时可能重新分布或外排。cLogP增加>1单位虽改善被动渗透性，但同时加剧CTS交叉反应性和代谢清除——这一权衡在CZ-CTS选择性景观中呈现为可量化约束。
3.3 代谢悬崖：硫代半卡巴腙作为警示故事
渗透性墙促使进一步的药物化学转向。若非肽骨架主要因细胞内暴露不足而失败，则结合充分渗透性与对CZ高内在反应性的化学型可能恢复抗寄生虫效价。硫代半卡巴腙作为对此逻辑的回应出现。这些化合物放弃肽识别元素，保留相对紧凑的分子结构，并纳入可与催化Cys25残基反应的亲电功能基团。
确实，硫代半卡巴腙是报道的对CZ最效的非肽不可逆类别之一。筛选鉴定了多个具有亚微摩尔酶促抑制和可测量的人CTS选择性的成员。抑制机制明确：Cys25硫醇负离子对硫羰基中心的亲核攻击产生共价硫脲加合物，在生化测定条件下产生持续靶点结合。
然而，这一表面成功掩盖了第三个决定性转化障碍。硫代半卡巴腙常表现不利的代谢稳定性。使Cys25结合得以实现的硫羰基亲电性，也使这些化合物易受谷胱甘肽（glutathione）缀合和非特异性蛋白加合的影响。Du等注意到尽管酶促效力强，杀锥虫活性在细胞系统中快速丧失——归因于代谢淬灭而非渗透性失败。通过增量化学修饰缓解这一局限的尝试收效有限：屏蔽亲电中心减少对脱靶亲核体的反应性，但同时降低对靶效力；增强亲电性则损害酶促抑制，同时可能增加代谢不稳定性。这些权衡反映了第3.1和3.2节中遇到的反相关优化景观。
4 计算方法：精确而无转化
该领域的计算投入可观：分子对接、分子动力学（MD）、自由能微扰（FEP）和基于机器学习的QSAR（ML-QSAR）均已应用于CZ且技术日趋精密。然而，日益增长的计算精密性未伴随临床转化。局限可能较少在于计算精度，而在于目标设定：所有四种方法均 predominantly 部署用于优化结合自由能（ΔG
bind或其替代指标），而细胞内暴露的多参数决定因素仍未被建模。
尽管方法多样，所有四种方法共享两个关键盲点：（i）无一模拟膜穿越、细胞内运输、溶酶体积累或代谢周转；（ii）均依赖隐式物理化学假设（最关键的是酶促质子化状态的选择），这些假设在感染相关条件下系统性失效。
4.1 共同失败模式：结合亲和力作为转化的误设替代指标
对于细胞内寄生虫靶点，结合亲和力是必要但次要的条件。治疗结果的决定因素是化合物能否到达寄生虫、在相关细胞内区室持续存在，并在不招致禁止性毒性的情况下维持足够游离浓度。优化ΔG
bind而忽视这些变量的计算方法将结构性优先注定失败的化合物，无论其在界定目标空间内实现的精度如何。
结构对接虽成功合理化了配体结合模式并识别S2亚位点关键相互作用，但其固有条件为理想化：刚性或半柔性蛋白、隐式溶剂环境，且无膜穿越、细胞内运输或pH依赖性分配的表达。分子动力学（MD）扩展了对接，纳入蛋白柔性、溶剂动力学和时间分辨配体-蛋白相互作用，但根本仍属局部范围：在纳秒至微秒时间尺度的均一水环境中进行模拟，不建模膜渗透、溶酶体积累、主动外排或代谢周转。自由能微扰（FEP）在预测同源系列相对结合自由能方面定量最严谨，在蛋白酶系统中与实验亲和力数据高度一致，但FEP在狭窄界定的化学空间内运作，优化单一目标ΔΔG
bind——在CZ背景下，FEP引导优化获得的增量亲和力增益常被渗透性、选择性或代谢稳定性的损失所掩盖。基于机器学习的QSAR（ML-QSAR）在酶促抑制数据集上应用时内部验证性能合理，但底层数据集通常较小（最大CZ特异性训练集<350个化合物）、化学异质且由生化而非细胞终点主导。
4.2 pH-质子化偏差：计算工作流程中的系统性错误
除共同的目标函数误设外，应用于CZ的计算工作流程嵌入了第二个常被忽视的系统性误差来源：对催化三联体中性pH质子化状态的假设。在pH ~7.0建模CZ（对接、MD和FEP工作流程的标准默认值）未能捕捉细胞内无鞭毛体经历的酸性微环境，其中酶在pH 5.0–5.5运作。
机制上，中性pH下His162仅部分质子化，Cys25保留大量亲核特征；而在溶酶体条件下，His162主要质子化，Cys25有效pKa预期升高，降低反应性硫醇负离子形式的种群。这些pH依赖性催化三联体质子化转移改变支配抑制剂结合和共价反应性的电子环境，削弱标准中性pH条件下测量或建模的亲和力可向感染相关酸性区室转移的假设。最清晰的实验例证为WRR-669系列，其在pH 5.5经历从不可逆共价抑制到pH 8.0完全可逆非共价结合的转变。
4.3 重新定位计算工具：从优化器至多参数过滤器
上述诊断并非计算方法与CZ药物发现无关，而是其角色必须根本重构：从单目标优化器至多参数过滤器和整合器。
在MPO框架内，每种方法在正确目的下保留价值。对接和MD模拟最适合重新定位为排除性过滤器：它们可有效消除无法合理结合CZ或依赖与选择性约束不兼容的相互作用的化合物。FEP最富成效地部署于严格界定的同源系列内，优先选择在不降低已满足渗透性或稳定性指标的情况下改善亲和力的取代。ML-QSAR模型在训练于多参数终点（整合生化活性与细胞效力、代谢稳定性和转运蛋白易感性）而非仅酶促效力时最为有用。
5 多参数优化（MPO）框架用于Cruzain抑制剂
诊断三元组要求根本不同的优化逻辑。图4形式化从顺序单参数优化（SPO）向MPO框架的过渡，其中效力、渗透性、代谢稳定性和CTS选择性被视为耦合而非顺序设计约束。
5.1 从顺序优化到约束满足
对于CZ，改善任何单一参数可预测地收紧其他参数约束。这一反相关景观要求根本不同的优化逻辑。MPO框架将优化重新构建为约束满足问题：候选化合物必须同时满足跨越正交维度的定义最低标准集合，而非最大化单一目标。任何单一约束失败的化合物均降优先，无论其他轴上的卓越表现如何。
5.2 定义Cruzain抑制剂的生存力阈值
酶促效力保持必要，但必须有界。建议亚100 nM生化IC
50范围作为效力筛选的实际起点，确保合理细胞内浓度可转化为靶点结合。对人CTSs L、B和S的选择性必不可少以避免机制性毒性；鉴于催化三联体和氧阴离子孔在CA族半胱氨酸蛋白酶中结构可叠合，至少两个数量级的选择性 margin 似乎是合理最小值。细胞内暴露为优势转化变量；化合物必须证明在感染宿主细胞中积累和持续存在的能力。代谢稳定性定义暴露的时间维度；微粒体或细胞t
1/2超过约2小时提供实际暴露持续时间基准。
这些标准非独立目标而系耦合约束。三维优秀但第四维失败的化合物不太可能转化。
5.3 未来Cruzain项目的战略意涵
采用MPO框架强制项目战略中的纪律性转变。针对效力最大化的大规模合成活动变得适得其反。反之，暴露和稳定性筛选的早期整合是强制性的，即使以舍弃高活性分子为代价。进展不以创纪录IC
50值衡量，而以跨所有约束保持活力的化合物的出现衡量。
5.4 恰加斯病药物设计背景下的Cruzain特异性挑战
CZ抑制剂的转化失败是否源于CZ的特异性，抑或在细胞内寄生虫靶点的药物设计中具有普遍性？答案微妙——问题部分通用但 critically 被CZ及其生物学背景的特异性特征所加剧。
通用组分已被充分认识：任何设计用于抑制寄居于宿主细胞内寄生虫中的酶的化合物，必须穿越多重生物膜、抵抗肝脏首过代谢、并在作用部位维持足够游离药物浓度。CYP51抑制剂用于CD的临床失败（泊沙康唑和fosravuconazole均显示初始寄生虫清除后频繁寄生虫学复发）确认CD中的转化障碍并非CZ独有。
但三个结构和生物学特征使反相关问题对CZ尤为严峻。第一，亚细胞定位：CZ主要位于细胞内无鞭毛体的溶酶体中，口服抑制剂因此必须穿越六个连续生物屏障——远大于疟原虫falcipain抑制剂面临的障碍，后者靶点寄居于去核红细胞的食物泡中，缺乏主动胞吞机制且存在寄生虫衍生渗透途径。第二，与人CTSs的结构同源性：CZ与人CTSs仅部分总体序列同一，但活性位点区域保守性显著更高。因此，有意义的选择性必须从S2和S3亚位点的细微差异中提取，而实践中改善选择性的取代基常增加极性和分子大小，直接对抗细胞内靶点可达所需的物理化学特征——形成选择性-渗透性权衡循环。第三，靶点冗余性：T. cruzi表达至少四种cruzipain（CZP）亚型，分为两个基因组家族，仅针对重组CZ优化的抑制剂因此对 in situ 遇到的同工型群体可能显示减弱疗效。
5.5 驾驭反相关：Cruzain抑制剂优化的可行性包络
若效力和类药性在CZ抑制剂空间中反相关，MPO能否真正解决这一张力？研究人员认为，记载的 anticorrelation 主要为经验性，反映单参数范式下实际采样的化学空间，而非绝对热力学约束——尽管不能完全排除部分热力学成分。
三条证据支持经验性解释。第一，CZ抑制的化学空间相对于其他类似年龄蛋白酶靶点仍较狭窄：215个独特抑制剂中，多数聚集于肽基骨架或少数非肽化学型，而类片段空间（MW ≤ 300 Da，cLogP ≤ 3）基本未探索。第二，配体效率（ligand efficiency, LE）分析显示，效力追踪分子量远强于配体效率（Pearson r(MW, pIC
50) = +0.610 vs. r(LE, pIC
50) = ?0.213），高效力通过分子增长而非更有效的活性位点结合实现。第三，相邻领域成功先例证明表观anticorrelation可在MPO ab initio 应用时被克服：HIV-1蛋白酶抑制剂从首代肽模拟物到达芦那韦（darunavir, K
i = 16 pM，口服生物利用度>80%，含利托那韦增强）的演化说明蛋白酶抑制剂空间在从头多参数设计下可航行。
关键概念转变从寻求各维度最优化合物，转向识别满足所有关键参数同时最低阈值——即可行性包络——的化合物。对于CZ抑制剂，现实的可行性包络可能由表3中的并发标准定义。关键洞察力在于：这些个体阈值无一要求极端优化；挑战完全在于其同时满足。硫代砜TSO-3和喹唑啉1s等实现中等效力与可接受选择性及细胞转化的化合物的存在，表明可行性包络非空。
6 未决问题与转化议程
将诊断框架转化为可操作进展的未决问题包括：pH依赖性细胞内分配是否优于被动渗透性预测细胞效力；早期多参数筛选能否在不减损疗效的情况下降低损耗；K777毒性是否由脱靶CTS抑制或骨架固有限质驱动；何种程度的CZ同工型结合对寄生虫学治愈是必需的；可调谐亲电体能否平衡效力与代谢稳定性；多参数计算管道能否在预期命中筛选中超越以亲和力为中心的管道。
7 结论与展望
三十年215个强效抑制剂未产生CZ临床候选药物；这一结果难以主要归因于靶点生物学、结构知识或计算精密性的缺陷，而反映优化策略与细胞内寄生虫感染疗效生物学决定因素之间的系统性错配。核心贡献在于将此问题重新构建于MPO框架内，将效力、选择性、细胞内积累的物理化学性质和代谢稳定性视为同时设计约束而非顺序里程碑。该框架不要求新实验技术；它需要优先级的重新排序，以及对早期权衡决定后期生存力的明确认识。对ChEMBL3563中已有215个效力合格抑制剂的数据集应用MPO标准进行系统重新筛选——需要新化学但仅需重新界定的目标函数——可识别约10–50个符合可行性包络的子集。这一预测可通过回顾性重新筛选直接验证。
