《Journal of Molecular Structure》：Indole-Tethered Thiazoles: Synthesis, Computational Evaluation as BACE-1 Inhibitor for Alzheimer’s Disease, and Photophysical Studies
编辑推荐：
阿尔茨海默病（Alzheimer's disease, AD）是一种以神经元进行性丧失为特征的神经退行性疾病，可导致学习能力缺失、记忆丧失及其他认知功能损害。目前AD已成为全球常见的痴呆形式，约4400万人患有痴呆，预计到2050年患者人数将达1.35亿，其中
阿尔茨海默病（Alzheimer's disease, AD）是一种以神经元进行性丧失为特征的神经退行性疾病，可导致学习能力缺失、记忆丧失及其他认知功能损害。目前AD已成为全球常见的痴呆形式，约4400万人患有痴呆，预计到2050年患者人数将达1.35亿，其中半数以上为AD病例。AD是一种多因素疾病，由乙酰胆碱耗竭、Aβ纤维聚集、tau蛋白过度磷酸化、N-甲基-D-天冬氨酸（NMDA）受体过度刺激、生物金属稳态改变、单胺氧化酶-B活性升高及氧化应激等因素引起。在AD发病机制涉及的酶中，β-分泌酶（即β-位点淀粉样前体蛋白裂解酶1，BACE-1）发挥关键作用。淀粉样前体蛋白（APP）经α-、β-和γ-分泌酶依次裂解，产生不溶性Aβ蛋白（Aβ40和Aβ42），其中Aβ42是AD病理生理学的重要因子。BACE-1催化的限速步骤产生不溶性Aβ，聚集导致斑块沉积并驱动神经退行性变，使其成为治疗AD的可溶性靶点。因此，开发有效的BACE-1抑制剂对AD治疗至关重要。
杂环化合物在药物化学中发挥重要作用，超过85%的生物活性化学实体含有杂环，超过60%的小分子药物含有氮原子。含氮和硫的杂环在AD治疗中具有重要意义。1,3-噻唑是由1位硫原子和3位氮原子组成的唑类化合物，是多种生物活性化合物的重要杂环骨架，如噻唑呋林（tiazofurin）、达沙替尼（dasatinib）、利托那韦（ritonavir）、雷夫康唑（ravuconazole）和硝唑尼特（nitazoxanide）等。虽然噻唑骨架作为胆碱酯酶抑制剂已被广泛研究，但其作为AD多靶点治疗的潜力研究有限。同样，多种吲哚骨架已被研究，但先前研究集中在单一药效团单元或替代靶点，尚无研究评估吲哚-噻唑杂合体作为BACE-1抑制剂的能力。本研究采用分子杂交策略，将吲哚和噻唑药效团单元结合到单一分子中，设计合成结构新颖的吲哚-四氢噻唑类化合物，靶向AD发病机制中涉及的β-分泌酶。研究通过合成、生物学筛选、分子对接、ADMET分析、DFT计算和光物理研究对新骨架进行了全面评估。
**论文解读：新型吲哚-四氢噻唑类BACE-1抑制剂的合理设计、合成与多维评估**
阿尔茨海默病（Alzheimer's disease, AD）是一种全球常见的神经退行性疾病，以进行性神经元丧失、学习能力缺失和记忆功能障碍为主要特征。目前全球约有4400万痴呆患者，预计2050年将达1.35亿。AD的病理机制复杂，涉及乙酰胆碱耗竭、Aβ纤维聚集、tau蛋白过度磷酸化、NMDA受体过度刺激、氧化应激等多因素。在AD发病机制中，β-分泌酶（β-site amyloid precursor protein-cleaving enzyme 1, BACE-1）催化淀粉样前体蛋白（amyloid precursor protein, APP）裂解的限速步骤，产生不溶性Aβ蛋白，聚集形成老年斑并驱动神经退行性变，因此BACE-1是治疗AD的关键药物靶点。现有治疗药物仅局限于胆碱酯酶抑制剂和NMDA受体调节剂，无法从根本上延缓疾病进展。杂环化合物在药物发现中占据核心地位，超过85%的生物活性分子含有杂环，其中噻唑骨架已在多种上市药物中展现出重要价值，但噻唑-吲哚杂合体作为BACE-1抑制剂的研究尚属空白。基于此，研究人员采用分子杂交策略，将吲哚与噻唑两种药效团单元整合于同一分子骨架，设计并合成了一系列结构新颖的吲哚-四氢噻唑衍生物，旨在开发新型BACE-1抑制剂以治疗AD。该研究发表于《Journal of Molecular Structure》。
研究人员通过Hantzsch噻唑合成法，以吲哚-2-甲醛与氨基硫脲经酸催化Schiff碱缩合生成硫代缩氨基脲中间体，再与系列溴代苯乙酮衍生物环化缩合，合成了8个目标化合物（5a-h）。所有化合物经FT-IR、1H NMR、13C NMR和质谱进行结构确证。计算方法方面，采用密度泛函理论（Density Functional Theory, DFT）在B3LYP/TZVP水平下进行基态几何优化和前沿分子轨道分析，激发态计算采用含时密度泛函理论（Time-Dependent DFT, TD-DFT）。分子对接使用Schr?dinger套件的Glide模块（Xtra Precision模式），靶蛋白选用BACE-1晶体结构（PDB ID: 2B8L，分辨率1.7 ?），对接协议通过共晶配体再对接验证（RMSD = 0.1864 ?）。结合自由能采用MM-GBSA方法计算。分子动力学（Molecular Dynamics, MD）模拟采用Desmond模块，OPLS4力场，TIP3P显式水模型，0.15 M NaCl生理离子强度，NPT系综，300 K条件下运行200 ns。ADME性质通过QikProp模块预测，毒性使用pkCSM网络工具评估。光物理性质采用紫外-可见（UV-Vis）吸收光谱（DMF中测定）和荧光发射光谱（DMSO中0.01 mM溶液，激发波长380 nm）表征。
**3.2.2 分子对接研究**
分子对接结果表明，所有噻唑衍生物（5a-h）均能进入BACE-1活性位点，展现出与共晶配体相似的结合模式。化合物5a通过离子相互作用直接锚定催化天冬氨酸双联体（Asp32和Asp228），并与Gly230形成氢键，其对接评分为-6.206 kcal/mol，结合自由能为-39.83 kcal/mol。其他衍生物（5b-h）的吲哚NH基团与Gly230形成氢键，苯环与Tyr71形成π-π堆积相互作用。其中，化合物5d（3,4-二氟取代）展现出最优的对接评分（-6.786 kcal/mol）和结合自由能（-53.61 kcal/mol），其次为5h（-6.542 kcal/mol）。所有化合物均能识别Gly230和Tyr71这两个关键残基，与共晶配体的结合模式具有良好相关性。DFT计算揭示的电子描述符与对接结果相互印证：卤素取代衍生物（5d、5f、5h）具有适中的能隙和硬度，有利于电荷转移和与Tyr71/Gly230的相互作用；硝基衍生物（5a）虽亲电性最高，但硝基的脱溶剂惩罚导致其对接评分较低；供电子取代基（5e、5g）降低亲电性，主要通过π-π堆积和疏水作用结合。
**3.2.3 MD模拟**
200 ns的分子动力学模拟显示，化合物5d的蛋白Cα骨架RMSD稳定在约2.0-2.8 ?范围内，而参比抑制剂verubecestat的结合体系RMSD维持在较高水平（约2.8-3.2 ?）。配体RMSD方面，5d在整个模拟过程中保持约3-4.2 ?的稳定值，而verubecestat在模拟后期（165-200 ns）出现高达约6.0 ?的位移，表明5d在结合口袋中的稳定性优于verubecestat。RMSF分析显示两种复合物的整体构象波动较小（约0.5-1.0 ?），5d在活性位点附近的残基78处表现出较高柔性（约2.6 ?），而verubecestat在230-240口袋区域具有更低柔性。蛋白-配体接触分析表明，5d主要通过疏水相互作用（与Tyr71、Phe108、Trp115）以及氢键和水桥（与Gln12、Gly230）维持稳定的接触网络；verubecestat则与催化残基Asp32和Gln73形成密集的极性接触。接触时间线显示，5d从75 ns起与Asp32、Phe108、Ser113和Gln12形成高度一致的稳定接触，而verubecestat的接触网络在周边残基区域表现出动态位移特征。
**3.2.4 ADME预测**
QikProp预测结果表明，所有化合物均满足Lipinski五规则，分子量低于500 Da，氢键供体不超过5个，氢键受体不超过10个。极性表面积（PSA）在46-94 ?2之间，显著低于verubecestat（118.9 ?2），预示更好的膜通透性。所有化合物的口服生物利用度预测值高达93-100%，优于verubecestat（82.3%）。QPPCaco值显示优异的肠通透性，QPlogBB值处于可接受范围（-3.0至1.2），血浆蛋白结合适中，表明具有良好血脑屏障穿透能力。pkCSM毒性预测显示所有化合物无hERG抑制风险，多数为AMES阴性（非致突变），LOAEL和MTD均在安全范围内。所有化合物均为P-糖蛋白（P-gp）底物，可能通过外排机制影响脑组织渗透，但其低PSA和高通透性提示被动扩散可能克服外排效应。
**SAR分析**
构效关系（Structure-Activity Relationship, SAR）分析表明：硝基取代类似物（5a）通过强吸电子效应降低HOMO能量并提高亲电性（ω=5.77 eV），与Asp32/Asp228形成离子相互作用，但脱溶剂惩罚降低其结合能。单卤素衍生物（5b、5c、5f）硬度增加、柔软度降低，通过π-π堆积和氢键实现中等结合能力。二卤素类似物（5d、5h）兼具吸电子和位阻效应，展示出优异的对接评分。供电子基团类似物（5e、5g）稳定HOMO、降低亲电性，增强分子亲核性，促进π-π和疏水相互作用但降低结合能。位阻效应影响取代基在活性位点的取向。
**光物理研究**
UV-Vis吸收光谱显示，化合物5a-h在DMF中于350-400 nm区间呈现最低能量吸收带，250-300 nm区间出现次级吸收峰，归属为π-π*分子内电子跃迁和分子内电荷转移（Intramolecular Charge Transfer, ICT）。吸电子取代基使吸收峰发生红移，表明这些取代基更好地稳定LUMO能级；供电子取代基（5e、5g）则导致轻微蓝移。荧光光谱显示所有衍生物的发射峰位于440-560 nm范围，呈现蓝绿色发光特性。其中5b-f和5h的最大发射强度位于约515 nm，5a和5g位于约450 nm。色度图证实发射坐标落在蓝绿区域。TD-DFT计算的吸收光谱与实验吻合良好，除5a外，约400 nm吸收带源于HOMO→LUMO激发（π→π*特性，高振子强度）；5a的该吸收带源自第三激发态（HOMO-1→LUMO）。发射计算显示，除5a外，所有分子的绝热发射能在约495-506 nm范围，与实验的约500 nm发射一致；5a的实验发射（约450 nm）可能源自更高激发态而非第一激发态。
**结论**
本研究成功通过Hantzsch噻唑合成法制备了一系列新型吲哚-四氢噻唑衍生物（5a-h），并经1H NMR、13C NMR、IR和质谱充分表征。分子对接揭示化合物对BACE-1酶活性位点具有强结合亲和力，MD模拟验证5d与蛋白复合物在200 ns内保持稳定。DFT计算阐明了取代基电子效应对反应性和稳定性的调控规律。体外ADME预测和pkCSM毒性评估证实了化合物具有良好的类药性和药代动力学特征，展现出作为β-分泌酶抑制剂的潜力。光物理研究显示所有噻唑衍生物的发射峰位于440-560 nm，呈蓝绿色发光，具备光电子应用前景。综合计算与结构研究结果，该噻唑-吲哚分子骨架为未来AD治疗研究提供了有前景的先导化合物方向。
