Mendelevium
Diary
Drug Design
Field Knowledge
Academia
Yang
Biology
Physics
Free Energy
Machine Learning & AI
Active Learning
Basics
Boltz-2
Data
Generation
Interpretability
QSAR application
Representations
Mol2Image
Workflow & Agent
Molecular Dynamics
FF & Algorithm
Small Molecule
martini
water
Interaction
Modeling & Tools
QM
Sampling & Analysis
Allostery
Fundamental
Other
Specific Sytems
Enzyme Engineering
Fiber & LLPS
Membrane
orientation_penetration
Metal
Nano Polymers
Skin Permeation
Techniques
Linux
Python
Research
Web
about
Home
Contact
Copyright © 2025 Xufan Gao | Academic Research Blog
Home
>
Field Knowledge
> Biology
A Bunch of Biophysics is Loading ...
Biology
把重复长度写成方程:HttEx1 与 polyQ 疾病的长度依赖模型汇总(临床—分子—凝聚态)
把重复长度写成方程:HttEx1 与 polyQ 疾病的长度依赖模型汇总(临床—分子—凝聚态) 多种AI调研的综合内容,请自行甄别信息的正确性 摘要 多聚谷氨酰胺(polyQ)疾病是一组由CAG重复序列扩增引起的神经退行性疾病,包括亨廷顿病(HD)、脊髓小脑共济失调(SCA)、齿状核红核苍白球路易体萎缩症(DRPLA)和脊髓延髓肌萎缩症(SBMA)等。这些疾病的共同特征是致病性与polyQ长度高度相关——重复序列越长,发病越早、症状越重。然而,现有文献中关于长度依赖性的描述往往是定性的,缺乏可直接用于定量建模的数学公式和参数范围。本文档系统整理了polyQ疾病的长度依赖模型,从临床预测、分子机制到(亚)细胞三个尺度,展示重复长度这一核心参数如何决定疾病表型。在临床尺度,我们总结了Weibull分布、对数线性模型等发病年龄预测公式,并量化了体细胞扩增、HTT1a异常剪接等修饰因子的影响;在分子尺度,我们整合了Httex1的成核-延伸模型、四聚体平衡方程和纤维延伸速率等动力学参数;在细胞尺度,我们分析了液-液相分离(LLPS)的Flory-Huggins相互作用参数$\chi(Q)$、临界浓度$C_{\mathrm{sat}}(Q)$以及膜相互作用和轴突运输缺陷的定量关系。总体上,以 $Q$ 为自变量的项可以归纳为多种函数族,常见函数族包括线性或分段线性、指数、幂律标度与S形函数,以及用于阈值刻画的变号判据。我们发现35-40个谷氨酰胺的临界阈值是polyQ疾病的普适特征,临床风险随Q长度呈连续陡峭的S形曲线,而分子层面在Q23-26和Q36-40存在类相变特征。本文总结的可复用定量框架为建立更准确的polyQ疾病预测模型和早期诊断干预策略提供了理论基础。 1. 研究背景 1.1 PolyQ疾病的共同分子特征 polyQ 疾病是一组由编码区 $\text{CAG}$ 重复扩增引起的遗传性神经退行性疾病。经典意义上的polyQ疾病通常指至少九种:Huntington 病(HD)、脊髓小脑共济失调 SCA1/2/3/6/7/17、DRPLA 与 SBMA。123 长度阈值与毒性正相关的普适规律 尽管致病蛋白在 polyQ 区段以外彼此并不相同,但近年来的系统综述揭示了一个高度普适的长度-毒性关系:当 polyQ 长度超过35-40个谷氨酰胺的临界阈值时,致病性才显著显现,且毒性与polyQ长度呈正相关——更长重复导致更严重的疾病表型和更早的发病年龄。2 不同 polyQ 疾病的致病阈值存在差异(32Q至54Q不等),这反映了不同蛋白背景对 polyQ 毒性的调制作用(表1)。 表1 不同polyQ疾病的致病阈值对比 疾病类型 致病阈值 备注 HD ≥36Q 典型致病阈值为36-40Q SCA1 ≥39Q 39Q以上外显率显著增加 SCA2 ≥32Q 致病阈值相对较低 SCA3 ≥52Q 致病阈值较高 SCA6 ≥19-21Q 较短,特殊类型 SCA7 ≥35Q 与HD类似 DRPLA ≥49Q 致病阈值较高 SBMA ≥38Q X连锁遗传 这种阈值效应和长度依赖性为理解 polyQ 疾病的共同机制提供了重要线索:polyQ长度可以作为预测疾病风险和进展的核心参数。 它们在分子层面有三条高度一致的共性主线,也是本文选择把长度作为统一自变量来整理模型的原因: 1.1.1 动态突变与代际提前:CAG不稳定把“长度”变成一个随机过程 polyQ 疾病的 $\text{CAG}$ 重复属于不稳定重复序列,可在代际传递中发生进一步扩增,并表现为代际提前(anticipation):子代更早发病且更重。HD 的人群研究显示,父系传递时重复更易净扩增,且父代到子代的重复变化与 AO 的变化相关,这为“代际提前”提供了直接的遗传学证据。43 因此,从建模角度,$\text{CAG}$ 不应只被视作一个固定参数,而应被视作随个体、组织与年龄演化的随机变量;本文第 2 章与第 2.1.2 节把体细胞扩增写成可进入 AO 方程的噪声项,作为跨尺度连接的第一步。56 更重要的是,“聚集”不仅是形态学现象,还会与蛋白稳态网络发生耦合:蛋白酶体、分子伴侣与自噬等系统既是清除通路,也可能在某些条件下促进包涵体形成或被聚集体重塑。以 HD 与 ataxin-3 为例,有研究显示 19S 相关的蛋白酶体伴侣亚基可与包涵体共定位,并能促进 mutant $\mathrm{HTT}$ 与 ataxin-3 的聚集,提示“PQC 只会抑制聚集”的直觉并不总成立。7 Enroll-HD 等队列显示,36–39 重复的“降低外显性”人群发病呈缓慢上升的累积风险(70 岁时,38 重复携带者的累积发病概率约 32%,39 重复约 68%),未见绝对突变式跳变,更像陡峭的 S 形曲线。8 40–42 重复的“晚发”患者与常规 40–42 重复群体在症状和进展上无显著差异,提示在阈值以上风险主要是连续变化,环境和修饰基因决定方差。9 1.2.2 分子层面:存在阶跃式相变特征 成核阶数跃迁:体外简单 polyQ 肽的关键核团簇大小在 Q23–Q26 出现阶跃:Q 从 26 降到 23 时,$N^{*}$ 由单体核变为二聚体核再到四聚体核,使短 Q 的聚集显著变慢;反过来 Q 增至 ≥26 时可由单体直接成核并明显加速,带来“相变感”。1011 聚集速率陡增:β-折叠概率和聚集速率在 Q36–40 附近提升 1–2 个数量级,与临床阈值共振;同时正常长度 polyQ 还能作为“助核剂”加速长链聚集,放大该阶跃效应。1213 浓度驱动的相分离:mHttex1 形成 M/S/F 三相(单体/球形寡聚体/纤维)存在尖锐浓度边界,并被 profilin 等配体左右,提示相变式分区是浓度与长度共同决定的。14 临床风险随 CAG 主要呈连续陡峭曲线,分子聚集动力学在 Q23–26、Q36–40 存在阶跃或类相变特征,两层信息应同时纳入模型。 本文主题:系统整理polyQ疾病的长度依赖模型 多聚谷氨酰胺(polyQ)疾病是一组由CAG重复序列扩增引起的神经退行性疾病,其核心特征是致病性与polyQ长度高度相关。从临床预测到分子机制,长度作为唯一明确的遗传参数,贯穿了疾病风险、发病年龄、蛋白构象、聚集动力学和细胞毒性等多个尺度。 然而,现有文献中关于长度依赖性的描述往往是定性的(如重复越长则发病越早),缺乏可直接用于定量建模的数学公式、参数范围和函数形式。本文档的目标是填补这一gap:系统整理polyQ疾病的长度依赖模型,将分散在临床研究、生物物理和细胞生物学中的定量关系,统一写成可复用的函数、速率方程和参数建议。 本文特别关注以下几个核心问题: 长度如何作为自变量进入预测模型? 从临床的Weibull分布到分子动力学的指数增长 哪些参数具有明确的Q长度依赖? 如成核速率、相分离临界浓度、降解速率等 如何实现跨尺度耦合? 将分子层面的长度依赖性质映射到临床时间轴 当前的研究缺口在哪里? 哪些函数关系尚未建立,哪些参数尚未测定 通过整理这些模型,我们希望为研究者提供一套可直接用于仿真、参数拟合和实验设计的工具集,并推动polyQ疾病研究从定性描述向定量预测转变。 2. 临床与群体级模型 2.1 Huntington病(HD) 2.1.1 发病年龄(AO)预测模型 Weibull 生存模型(Langbehn 2004) Langbehn 模型假设 AO(age of onset,发病年龄)服从 Weibull 分布(一种用于生存分析的参数分布,可描述随时间变化的风险),log-scale 参数 $\mu = 0.053 \times (\text{CAG} - 35.5)$($\mu$ 为位置参数,控制分布中心位置),形状参数 $k = 20.2$($k$ 控制分布形状和危害率随时间的变化),允许根据 CAG 重复数计算任一年龄的累积分布函数(cumulative distribution function, CDF)并驱动危害函数(hazard function,描述在某一时刻发病的瞬时风险)。15 用更标准的Weibull生存分析写法,可把“给定CAG长度的发病年龄分布”写成生存函数 $S(t\mid Q)$、累积发病概率 $F(t\mid Q)$ 与瞬时危害率 $h(t\mid Q)$。如果用 $\mu(Q)$ 表示 log-scale 参数,则可取尺度参数 $\lambda(Q)=\exp[\mu(Q)]$:15 \[\mu(Q)=0.053\left(\text{CAG}-35.5\right), \quad \lambda(Q)=\exp\!\bigl[\mu(Q)\bigr], \quad k=20.2.\] 在Weibull分布下: \[S(t\mid Q)=\exp\!\left[-\left(\frac{t}{\lambda(Q)}\right)^k\right], \quad F(t\mid Q)=1-S(t\mid Q).\] 对应的危害函数为: \[h(t\mid Q)=\frac{k}{\lambda(Q)}\left(\frac{t}{\lambda(Q)}\right)^{k-1}.\] 公式的通俗解释 这组公式的用途是把遗传长度 $Q$ 直接映射成任意年龄 $t$ 的发病风险:$F(t\mid Q)$ 给出到 $t$ 岁之前已经发病的累计概率,$S(t\mid Q)$ 给出到 $t$ 岁仍未发病的概率,$h(t\mid Q)$ 则描述在 $t$ 附近的瞬时发病风险随年龄如何变化。参数 $k$ 主要控制风险曲线是否随年龄变陡,$\lambda(Q)$ 则把整条风险曲线沿年龄轴平移并随 $Q$ 改变。 概率密度函数: \(\begin{aligned} &f(x ; \lambda, k)= \begin{cases}\frac{k}{\lambda}\left(\frac{x}{\lambda}\right)^{k-1} e^{-\left(\frac{x}{\lambda}\right)^k} & x \geq 0 \\ 0 & x<0\end{cases}\\ &\text { 其中,} x \text { 是随机变量,} \lambda>0 \text { 是比例参数(scale parameter),} k>0 \text { 是形状参数 } \end{aligned}\) https://zh.wikipedia.org/zh-cn/%E9%9F%A6%E4%BC%AF%E5%88%86%E5%B8%83 分段 YTO 公式 同一研究还提供了基于TRACK-HD/Enroll-HD队列拟合的距发病年份预测模型:15 \[\mathrm{YTO}(Q)= \begin{cases} -20.854 - 0.886\, (Q - 44), & 40 \le Q \le 44, \\ -9.653 - 1.494\, (Q - 44), & Q \ge 45 . \end{cases}\] 公式的通俗解释 方程以CAG重复数 $Q$ 为自变量,输出预估的距发病年份 $\mathrm{YTO}$。 分段特征 在 $Q\le 44$ 时,斜率约为 $-0.89$;当 $Q\ge 45$ 时,斜率陡增到约 $-1.49$,意味着每多1个重复,平均发病年龄额外提前约1.5年。系数来源于 TRACK-HD/Enroll-HD 队列拟合,可直接用于人群风险模拟和个体化随访。 对数线性 AO 回归(Lee 2012) 在纳入正常等位基因后可写成:16 \[\ln(AO) = \beta_0 - 0.049\,\text{CAG}_{\mathrm{exp}} + 0.013\,\text{CAG}_{\mathrm{norm}}\] 公式的通俗解释 这个公式预测HD患者的发病年龄。其中 $\beta_0$ 为截距常数,$\text{CAG}{\mathrm{exp}}$ 为扩增的病理性等位基因的CAG重复数,$\text{CAG}{\mathrm{norm}}$ 为正常等位基因的CAG重复数。 核心规律 系数 $-0.049$ 表示病理性CAG每增加1个重复,$\ln(\mathrm{AO})$ 降低 $0.049$(即发病提前);系数 $+0.013$ 表示正常等位基因具有缓冲作用,$\text{CAG}_{\mathrm{norm}}$ 越长,AO 略微推迟。 体细胞扩增噪声项(Swami 2009) 纹状体(striatum)与皮质(cortex)样本显示,独立测得的体细胞扩增量(somatic expansion, SE,指CAG重复在有丝分裂后神经元中的进一步扩增)与 AO 呈线性关系(每 SD 扩增约对应 AO 提前 3.3 年),可作为 length-dependent 噪声项:5 \[AO \approx f(\text{CAG}_{\mathrm{germline}}) - 3.3 \times \mathrm{z}(SE)\] 公式的通俗解释 这个公式解释为什么具有相同遗传CAG长度的患者发病年龄会有差异。其中 $\text{CAG}_{\mathrm{germline}}$ 为生殖细胞遗传的CAG长度,$\mathrm{z}(SE)$ 为体细胞扩增量的标准化值(z-score,以标准差SD为单位)。 核心规律 脑组织中的CAG会持续扩增,每多1个标准差的体细胞扩增,发病年龄会额外提前约3.3年(文献报道为约$2$–$4\,\mathrm{年/SD}$,$3.3$为中值近似)。这解释了同一家族内、相同CAG长度患者AO的个体差异。 指数衰减备选模型(Poirier 2002,待验证) 除上述主流模型外,文献中还提出了指数衰减模型:17 \[AO_{\mathrm{alt}}(Q) = \exp\!\left(6.657 - 0.0662\,Q\right)\] 公式的通俗解释 该式把发病年龄写成随CAG长度指数衰减的形式:每多1个重复,$AO_{\mathrm{alt}}$ 按因子 $\mathrm{e}^{-0.0662}$ 缩短。 注意 因缺乏同行评审,仅作模型对照,不替代上述主流分段/对数线性模型。 另一版本为: \[AO_{\mathrm{alt}}(Q) = \exp\!\left(5.34 - 0.0363\,Q\right)\] 斜率更缓(-0.0363);若使用,可与分段/对数线性模型交叉验证,评估是否过度拟合或低估高 Q 区域风险。 2.1.2 遗传修饰因子 HTT1a 转录本(Hoschek 2024) HTT1a 是一种通过不完全剪接(aberrant splicing)产生的仅含 exon 1 的短转录本,由于exon 1到exon 2的剪接受阻,随后在intron 1中切割并加尾,最终翻译成高度聚集和毒性的exon 1蛋白片段。HTT1a/HTT(全长huntingtin转录本)比例随 CAG 每增加 10 个重复上升约 0.15,可作为毒性剂量代理,纳入多层贝叶斯模型:18 \[\mathrm{HTT1a}/\mathrm{HTT} = 0.015 \times \text{CAG} + \epsilon\] 公式的通俗解释 这个公式预测异常剪接产物HTT1a相对于全长HTT的比例。HTT1a是一种仅含exon 1的短转录本,翻译后形成高度聚集和毒性的片段,是HD病理的关键驱动因素。 核心规律 HTT1a/HTT比例随CAG长度线性增加,系数$0.015$表示CAG每增加10个重复,HTT1a比例上升约$15\%$。其中$\epsilon$为个体间变异的误差项,反映不同患者剪接效率的差异。 m6A RNA修饰调控 METTL3(m6A甲基转移酶复合物的核心催化亚基)和YTHDF1/3(m6A阅读蛋白,识别并调控m6A修饰的mRNA)上调会抑制 HTT1a 的翻译或稳定性,等效于在 HTT1a 方程添加负反馈项 $-k_{m6A} \times \text{METTL3}{\mathrm{act}}$,其中 $k{m6A}$ 为m6A介导的HTT1a抑制速率常数,$\text{METTL3}_{\mathrm{act}}$ 为METTL3的活性水平。该机制提示m6A修饰是调控HTT1a毒性的潜在治疗靶点。19 2.1.3 疾病进展与神经影像 神经元层读数:皮层组织中体细胞扩增与萎缩速度的共变,提示在 state-space 模型中需把 SER(somatic expansion ratio)与 MRI 体积变化耦合。6 2.2 其他polyQ疾病 SCA(spinocerebellar ataxia,遗传性脊髓小脑性共济失调)是一组常染色体显性遗传的小脑变性病,编号 1、2、3、6、7、17 等对应不同致病基因,其中多种属于 polyQ 扩增疾病。 2.2.1 SCA1 / SCA2 / SCA3 / SCA6 对数 AO 回归:四种 SCA 的最佳模型均为:20 \[\ln(AO) = \beta_0 + \beta_1 \text{CAG}_{\mathrm{exp}} + \beta_2 \text{CAG}_{\mathrm{norm}}\] 公式的通俗解释 这个公式形式与HD的对数线性模型一致,但不同SCA疾病的系数差异显著,反映了不同polyQ蛋白对CAG长度的敏感性不同。 具体系数 SCA1($\beta_1=-0.049$,$\beta_2=+0.013$):系数数值与HD模型巧合一致,但来自 Tezenas du Montcel 2014 的独立回归;正常等位基因同样表现为缓冲作用。2016 SCA2($\beta_1=-0.105$):对CAG长度最敏感,病理性CAG每多1个重复的影响约为亨廷顿病的2.1倍 SCA3($\beta_1=-0.056$):对CAG长度的敏感性略高于亨廷顿病 SCA6($\beta_1=-0.090$,$\beta_2=-0.029$):特殊的是正常等位基因系数为负,表明正常等位基因的CAG重复越长,发病越早 这些系数可直接用来比较不同polyQ蛋白对AO的敏感度。 功能退变速率与生存危害(EUROSCA 队列):2122 疾病 SARA 年进展率(分/年) 每多1个 CAG 增加的SARA 斜率 每多1个 CAG 的死亡危害比 备注 SCA1 2.11 +0.06 1.06 进展最快 SCA2 1.49 +0.04 1.16 CAG 对 hazard 最敏感 SCA3 1.56 +0.03 1.08 中等 SCA6 0.80 +0.02 1.05 进展最慢,hazard 影响最小 说明:SARA(Scale for the Assessment and Rating of Ataxia)是评估共济失调严重度的 8 项量表,总分 0–40 分,分数越高表示症状越重,常用于 SCA 等小脑变性疾病的纵向随访。 解释:$\mathrm{Pr}(\text{juvenile})$ 表示少年型表型概率(青少年起病,常见肌阵挛性癫痫及快速退化),$\sigma$ 为 logistic 函数;当 $Q=63$ 时概率为 0.5,$Q$ 每增加 1,概率按 sigmoid 曲线陡增。 2.2.2 SCA7 AO vs CAG:在 25 家族 131 名患者中,AO 与 CAG 呈强负相关($r=-0.84,\,p<0.001$),拟合线约 $AO = 102 - 1.7 Q$;强调 retina/脑干受累的急剧斜率。23 疾病阶段模型:将 AO、一年内病程与呼吸功能分阶段建模显示,CAG 长度决定3个阶段间的转移概率,尤以 >60Q 者迅速进入 Stage 2–3,可构建 Markov chain。24 眼科指标:角膜内皮细胞密度(ECD)与 CAG 呈线性下降($\mathrm{ECD} \approx 3171 - 48\times Q$),可把视网膜病程量化为长度函数。25 2.2.3 SCA17 线性穿透力:多中心队列显示 AO 与 CAG 呈负线性,约 47 重复是成人与少年表型分界,可写成 \(AO = 119 - 1.4 Q\) 近似式。26 脑结构关联:体素基形态学(VBM)分析发现,小脑灰质体积与 CAG 呈线性递减($R^2=0.33$),为三维结构模型提供约束。27 2.2.4 Dentatorubral-pallidoluysian atrophy (DRPLA) 表型分类:重复数 62–79 对应少年肌阵挛性癫痫表型,49–71 对应痴呆/共济失调,提供了piecewise 逻辑映射28 \(\mathrm{Pr}(\text{juvenile}) = \sigma(0.37(Q-63))\) AO 与 CAG:AO 与 CAG 负相关($r=-0.696,\,p<0.001$),回归式约 \(AO = 132 - 1.7 Q\)。29 疾病里程碑:CAG 长度越高,步行/轮椅/死亡的转换时间越短;每增加一个重复,步行→轮椅时间缩短 0.26 年,可直接用于多状态模型。30 解释:$\mathrm{Pr}(\text{juvenile})$ 表示少年型表型概率,$\sigma$ 为 logistic 函数;当 $Q=63$ 时 $\mathrm{Pr}=0.5$,$Q$ 每增加 1,概率按 sigmoid 曲线陡增。 2.2.5 Spinal and bulbar muscular atrophy (SBMA) Meta 模型:系统回顾 1,317 名患者显示 AO 与 CAG 呈线性(slope ≈ -1.3 年/Q),$R^{2} \approx 0.34$。31 人群数据:韩国157例回归式$AO = 92.7 - 1.21 Q$($r=-0.407$),肌力(MRC)和功能(ALSFRS-R)与CAG亦呈显著负相关。32 内分泌调制:多元回归揭示血清睾酮、SHBG 与 CAG 共同解释握力/行走时间差异,可建模为:33 \[\text{Strength} = \alpha - 0.35 Q + \beta T + \gamma \text{SHBG}\] 公式的通俗解释 这个公式揭示了SBMA疾病特有的激素调控机制。SBMA由雄激素受体(AR)基因中的polyQ扩增引起,因此激素水平会显著影响疾病表型。 其中 $\alpha$ 为截距常数,$Q$ 为CAG重复数,$T$ 为血清睾酮水平,$\text{SHBG}$ 为性激素结合球蛋白水平。 核心规律 CAG每多1个重复,肌力下降0.35个单位,但这个效应可以被激素水平调控。系数 $\beta$ 和 $\gamma$ 分别表征睾酮和SHBG对肌力的影响,提示内分泌治疗可能是SBMA的潜在干预靶点。 2.2.6 跨疾病比较 | 疾病 | 模型类型 | 长度效应示例 | | — | — | — | | SCA1 | $\ln(AO)$ 线性、多态 hazard | -0.049 per repeat (AO), HR 1.06 | | SCA2 | $\ln(AO)$ 线性 | -0.105 per repeat, HR 1.16 | | SCA3 | logistic AO + progression | -0.056 per repeat, SARA +0.03/yr/repeat | | SCA6 | 超线性(正常等位基因也显著) | $\beta_{\exp}=-0.090$, $\beta_{\mathrm{norm}}=-0.029$ | | SCA7 | AO & 眼科线性 | $r=-0.84$, ECD -48 cells/Q | | SCA17 | AO 线性 + VBM | AO slope ≈ -1.4 年/Q | | DRPLA | piecewise phenotype, milestone hazard | juvenile σ(0.37(Q-63)) | | SBMA | AO 线性 + 激素交互 | AO slope -1.21 年/Q,strength 公式见上 | 2.3 参数速查表 模块 关系 用途 参考 AO 分布 Weibull: $\mu = 0.053 (Q-35.5), k=20.2$ 生存/危害仿真 15 对数 AO $\ln(AO) = \beta_0 -0.049 Q_{\mathrm{exp}} + 0.013 Q_{\mathrm{norm}}$ 个体化预测 16 Somatic 噪声 $AO = f(Q) - 3.3 \times \mathrm{z}(SE)$ 解释残差 5 HTT1a 剂量 $0.015 \times Q$ 分子剂量代理 18 神经影像 SER与MRI体积共变 State-space模型 6 3. 分子尺度模型 3.1 Httex1结构与构象动力学 3.1.1 结构域组成与功能 Huntingtin exon 1(Httex1)由三个主要结构域组成: N17结构域:前17个氨基酸,形成两亲性α-螺旋结构。该结构域在脊椎动物中高度保守,可插入脂双层膜,促进mHttex1聚集。N17的两亲性螺旋性质使其既能与膜相互作用,也能介导蛋白-蛋白相互作用。343536 PolyQ核心区:可变长度的谷氨酰胺重复序列,是HD及其他polyQ疾病的致病核心。PolyQ区在病理性长度(Q>35)时采用长α-螺旋构象,由谷氨酰胺侧链到骨架氢键传播和稳定。该区域的长度直接决定聚集倾向和疾病严重程度。37 PRD(proline-rich domain):富含脯氨酸的C端结构域,包含P10和P11两段polyproline序列,由短的富脯氨酸序列连接。PRD通常形成不完美的polyproline II(PPII)螺旋构象,在纤维表面占据显著比例,参与蛋白-蛋白相互作用和疾病病理。34 这三个结构域的协同作用决定了Httex1的整体构象、聚集动力学和细胞毒性。 3.1.2 单体构象特征 α-螺旋稳定性 局部构象能量:固体 NMR + 分子动力学表明,N17(huntingtin exon 1的前17个氨基酸,形成两亲性α-螺旋结构域)α-螺旋稳定性与 polyQ 长度呈正相关,PRD(proline-rich domain,富含脯氨酸的C端结构域,包含P10和P11两段polyproline序列)的构象熵与毒性表型强相关,提示能量函数可写成:34 \[E = E_0 + \alpha_{N17} \times Q + \beta_{PRD} \times \Delta S_{PRD}\] 公式的通俗解释 这个能量函数描述了Httex1蛋白单体的构象稳定性如何随polyQ长度和flanking区域的构象变化。 其中 $E_0$ 为基础能量,$\alpha_{N17}$ 为N17螺旋稳定性对Q长度的敏感系数(在本符号约定下,若Q变长使N17稳定化并降低能量,则 $\alpha_{N17}<0$),$\beta_{PRD}$ 为PRD构象熵的能量权重,$\Delta S_{PRD}$ 为PRD相对于参考态的构象熵变化。 核心规律 随着polyQ长度增加,N17结构域的α-螺旋更稳定,同时PRD的构象灵活性变化也会贡献能量,这两个flanking区域的构象变化共同影响蛋白的整体稳定性和聚集倾向。 参数符号提示:在上式约定下,若N17稳定化使能量降低,则 $\alpha_{N17}<0$;若Q变长增强聚集驱动力并降低 $\Delta G_{\mathrm{agg}}$,则 $\beta_Q<0$。 跨域互作 $k_{\mathrm{off}}$:多域 NMR(DEST/CPMG)测得 N17–PRD 相互作用的解离速率常数从 Q25 的约 $20\,\mathrm{s^{-1}}$ 降到 Q46 的约 $8\,\mathrm{s^{-1}}$,说明 polyQ 拉长强化跨域瞬时互作,使单体更紧凑,可用于多状态马尔可夫模型。已公开实验仅有 Q25 与 Q46 两点;Q30–Q80 尚无实测 $k_{\mathrm{off}}$。长时程 MD(Mohanty 2025)与 EPR 模拟均预测更长 Q 会进一步降低 $k_{\mathrm{off}}$,但缺少实验数值,属待补数据。38 α-螺旋稳定性作为主导因子(来自2023年研究) Elena-Real等(2023)通过位点特异性同位素标记NMR揭示,病理性httex1(Q46和Q66)的polyQ区采用长α-螺旋构象,由谷氨酰胺侧链到骨架氢键传播和稳定。整合数据分析表明,α-螺旋稳定性是比Q数量更强的聚集动力学和纤维结构特征。该发现提示在能量函数中应优先考虑螺旋稳定性项: \[\Delta G_{\mathrm{agg}} = \Delta G_{\mathrm{helix}} + \beta_Q \times Q + \Delta G_{\mathrm{flank}}\] 聚集驱动力随 Q 增强,故在 $\Delta G_{\mathrm{agg}} = \Delta G_{\mathrm{helix}} + \beta_Q \times Q + \Delta G_{\mathrm{flank}}$ 中,$\beta_Q<0$,可用 −0.02 至 −0.08 $\,\mathrm{kcal\cdot mol^{-1}\cdot Q^{-1}}$ 作为初始量级。此符号方向与固体 NMR/SMFRET 及 Mohanty 2025 长时程 MD 预测一致。 公式的通俗解释 该式表明聚集驱动力由三部分组成:其中 $\Delta G_{\mathrm{helix}}$ 为螺旋稳定化能(负值表示稳定),对聚集速率的贡献大于线性Q项。37 核心规律 聚集驱动力 = 螺旋稳定能(越稳定越易聚集)+ Q长度线性项 + flanking序列贡献。螺旋稳定能主导,意味着即使Q数相同,螺旋更稳定的变体也更易聚集。 超紧密单体态 smFRET(单分子荧光共振能量转移,可测量蛋白末端间距离)+ 分子动力学显示,Httex1 的回转半径(radius of gyration, $R_g$,衡量蛋白整体尺寸)与 Q 长度遵循标度律 $R_g \propto Q^{\nu}$,其中 $\nu \approx 0.22$,明显低于典型球状蛋白($\nu \approx 0.33$)或无序链($\nu \approx 0.5$,理想链),源于谷氨酰胺侧链的高密度内聚氢键网络形成超紧密构象。该异常低的标度指数表明polyQ蛋白并非典型的无序蛋白。参数 $\alpha \approx 2.62$Å 的来源待补充验证。39 \[\langle R_g\rangle = \alpha\,N^{\nu}\] 单体 $\beta$ 倾向与能垒 Vitalis 等的分子模拟指出,单体 polyQ 形成 $\beta$-rich 状态的自由能代价可达 $10$–$20\,\mathrm{kcal\cdot mol^{-1}}$,且长度依赖并不支持“单体 $\beta$ 核”随长度显著变得容易的简单图景,更合理的解释是 $\beta$-sheet 更像“肽富集相”中的属性。40 Jakubek 等用 UVCD 等手段显示,polyQ 的溶液态构象分布会随重复长度偏移,长度增加可提高 $\beta$-strand 倾向的占比,从统计意义上降低“找到成核构象”的等待时间,可作为“滞后期随 Q 指数缩短”的分子起点(推断)。41 3.1.3 浓度依赖的结构转变 非病理性HttEx1-17Q的浓度依赖结构转变(Yoo 2025) Yoo等(2025)通过NMR、CD、TEM和AFM系统研究了非病理性HttEx1-17Q在不同浓度下的结构转变:单体在低浓度下以无序为主,随浓度升高经历无序→螺旋→β结构的多重转变,这种重排可加速短淀粉样纤维成核。该发现为理解聚集早期事件提供了浓度阈值参数,可写成分段函数: \[\text{Structure}(c) = \begin{cases} \text{Disordered}, & c < c_1, \\ \text{Helical}, & c_1 \le c < c_2, \\ \beta\text{-sheet}, & c \ge c_2 . \end{cases}\] 其中 $c_1$ 和 $c_2$ 为实验测定的临界浓度(HttEx1-17Q约在 μM 量级)。42该式描述蛋白构象随浓度的分段变化,$c_1$ 和 $c_2$ 是两个临界浓度阈值,可从CD光谱和ThT荧光实验确定。 需要强调:这种分段写法是对给定实验条件下的经验近似,并不必然对应严格意义上的数学不连续。多种polyQ体系中确实常见在某段长度或浓度区间变得更陡的现象,但更常见的解释是连续但高度非线性的加速叠加机制切换(例如成核阶数变化),从而呈现类似阈值的观感。12 3.2 聚集动力学 PolyQ蛋白的聚集是一个多步骤、多尺度的复杂过程,从单体的构象转变开始,经历成核、寡聚化、纤维延伸,最终形成成熟的淀粉样纤维。这一过程的动力学特征高度依赖于polyQ长度,并在Q23-26和Q35-40两个关键阈值处表现出显著的非线性转变。 本节系统梳理聚集动力学的核心机制,整合实验测量、计算模拟和临床数据,构建从分子到疾病的定量连接。 3.2.1 成核机制:从单体到临界核 成核是聚集过程的限速步骤,决定了聚集的滞后期和整体速率。PolyQ成核的独特之处在于其长度依赖的阶数跃迁和自毒性机制,这些特征共同塑造了疾病的阈值效应。 成核动力学常数与能垒 对Q47肽段的热力学/动力学解析测得关键参数:43 二级聚集速率常数:$k_{+} \approx 1.14 \times 10^{4}\,\mathrm{M}^{-1}\,\mathrm{s}^{-1}$ 核形成平衡常数:$K_{n^*} = 2.6 \times 10^{-9}$ 核形成自由能:$\Delta G_{n^*} \approx +12.2\,\mathrm{kcal\cdot mol^{-1}}$ 正值的自由能表明成核是热力学不利过程,需要跨越显著能垒。该体系支持单体核($n^*\approx 1$)的解释框架。 成核时间的指数依赖关系 HD患者数据分析和polyQ肽实验表明,成核滞后期与Q长度呈指数关系。在半对数坐标图上,成核时间对Q长度呈线性下降趋势: \[\log(t_{\mathrm{lag}}) \propto -Q\] 这一现象表明成核速率随polyQ长度指数增长。Perutz等人最早提出假说:单体形成临界β-结构核的概率随链长指数增长,因为每个额外的谷氨酰胺都增加了形成稳定β-sheet核心的机会。 实验观察 病理性长度polyQ(如HTTex1 Q46):聚集时间尺度为小时至天 正常长度polyQ(Q<35):可保持溶解状态数周至数月 发病年龄在Q35-40附近急剧下降,反映了聚集倾向的陡增 这种指数关系解释了为何Q长度微小的增加(如从36到40)会导致发病年龄的巨大差异。 经典成核-延伸方程 这些参数可直接参数化到经典成核-延伸模型: \[\begin{aligned} \frac{\mathrm{d}[M]}{\mathrm{d}t} &= -k_n[M]^{n^{*}} - k_{+}[M][F], \\ \frac{\mathrm{d}[F]}{\mathrm{d}t} &= 2k_n[M]^{n^{*}} + k_{+}[M][F]. \end{aligned}\] 其中: $k_n$:成核速率常数,与$K_{n^}$相关:$k_n \propto k_{+}^2 K_{n^}$ $n^*$:临界核尺寸(Q47约为1) $[M]$:单体浓度 $[F]$:纤维末端浓度 第一式描述单体通过成核和延伸消耗,第二式描述纤维末端通过成核产生(因子2表示每个核有两端)和延伸增加。 成核阶数的长度依赖跃迁 Wetzel团队对双赖氨酸封端的polyQ肽(如$K_{2}Q_{n}K_{2}$)的系统研究揭示了成核阶数在Q23-26区间的离散跃迁:1011 Q≤23时需要四聚体核:短链polyQ(Q≤23)需要形成四聚体核($N^{*}=4$)才能启动淀粉样生长,多个单链必须协同聚集才能克服成核能垒,导致聚集速率显著降低 Q25附近过渡为二聚体核:在Q25附近,成核阶数从四聚体平滑过渡到二聚体($N^{*}\approx 2$),所需协同的单链数量减少,成核效率相应提高 Q≥26以单体核为主:当polyQ长度达到或超过26个谷氨酰胺时,单体即可作为成核核心($N^{*}\approx 1$),无需多链协同,成核速率大幅提升并呈现明显的加速效应 即使在$N^{}\approx 1$的长链区间,聚集速率仍随Q继续上升,因为成核效率项$k_{+}^{2}K_{N^{}}$对Q呈非线性增强,产生类似相变的加速效应。 此外,临床数据驱动的动力学映射也把“重复长度依赖的成核”与“成核依赖的延伸”区分开来,支持长度效应主要由成核步骤主导的建模框架。44 重要区分 Q23-26的阶跃主要描述体外简单polyQ肽的物理转折区,并不等同于HD的临床致病阈值(Q36-40)。更高的临床阈值很可能来自flanking区域(N17、PRD)、细胞环境与体内稳态网络对能垒的再塑形。 多长度 Httex1 的实测动力学数据(Vieweg 2016) Vieweg 等使用 intein 纯化获得 tag-free Httex1,并对多种 Q 长度在统一条件下做 ThT 动力学与电镜表征,给出了一组可直接用于拟合的长度依赖数据矩阵。45 以 $15\,\mathrm{\mu M}$、$37\,^{\circ}\mathrm{C}$ 的典型条件为例: 构建 完全聚集时间(小时) 平均纤维长度(nm) 备注 Httex1–23Q 1176 522 作为慢聚集基准 Httex1–29Q 124.5 284 完全聚集快约 9 倍 Httex1–37Q 242 268 完全聚集快约 4.9 倍 Httex1–43Q 5.98 183 完全聚集快约 200 倍 所以这是突变?我没看原文 Q 增长不仅缩短聚集时间,还缩短终态纤维平均长度,提示长度同时改变成核与碎裂/封端动力学(推断)。45 临界浓度与能量分解(Sahoo 2014,Httex1–Q23/Q42) Sahoo 等用化学合成的全长 Httex1(含 PRD)直接比较 Q23 与 Q42 的聚集过程,并把平衡态单体浓度写成临界浓度 $C_r$。该工作给出了多个构建的 $C_r$,使得 flanking 区域对稳定性的贡献可以转成 $\Delta\Delta G$:46 Httex1–Q23 的 $C_r = 0.44\pm0.13$ μM。 HTT$^{NT}$Q23P10K2 的 $C_r = 0.28\pm0.11$ μM。 K2Q23K2 的 $C_r = 3.0$ μM(文中引用既往结果)。 HTT$^{NT}$Q23K2 的 $C_r \le 0.1$ μM(受检测限约束)。 基于这些 $C_r$,作者估算:HTT$^{NT}$ 对 Q23 纤维稳定性贡献至少 $2.2\,\mathrm{kcal\cdot mol^{-1}}$,而 PRD 对稳定性贡献至少 $0.9\,\mathrm{kcal\cdot mol^{-1}}$ 的去稳定化。46 同时,该工作还提供了两个对“长度→动力学”建模很关键的事实: 早期出现球形寡聚体(约 100–600 个分子)与短 protofibril(约 500–2600 个分子)。46 成熟纤维的解聚与重新达到平衡是一个月量级的过程。46 β-螺旋几何约束与Q35-40阈值 计算模拟与结构研究提示,Q35-40临床阈值可能源于β-螺旋形成的几何要求: β-螺旋每转包含约18.5±2个残基,因此要形成稳定β-螺旋,通常需要约33-40个谷氨酰胺残基 模拟对比显示,Q25在所测试温度范围内难以形成稳定β-螺旋,而Q45在更宽温度范围内可稳定形成β-螺旋 随机线团→平行β-螺旋的构象转变更倾向发生在超过约37个谷氨酰胺残基的肽段 这一几何约束与临床表型吻合:$Q<35$通常不致病,$35\le Q\le 39$可能致病,$40\le Q\le 60$多为成人发病,$Q>60$与少年型HD相关。来源 建模意义 β-螺旋形成可作为构象转变的order parameter,临界长度$Q_c \approx 37$对应β-螺旋稳定性的相变点。成核速率可表示为: \[k_{\mathrm{nuc}}(Q) = k_0 \exp\bigl[\gamma (Q - Q_c)\bigr]\] 其中$k_0$为参考速率,$\gamma$控制长度敏感度,$Q_c$为经验阈值(约35-40)。指数型形式表示一旦Q超过阈值,成核速率呈倍数级飙升,对应临床上更长CAG重复往往更早发病的现象。47 自毒性机制与steric zipper结构 Kandola等(2023)使用细胞内直接报告系统揭示了polyQ成核的自毒性机制:来源 核心发现 病理性polyQ成核涉及steric zipper结构:polyQ核形成涉及每隔一个位置的三个谷氨酰胺残基段,编码四链steric zipper结构,侧链交叉指状排列形成稳定的疏水核心,这种特殊结构是病理性聚集的关键分子基础 临界长度对应结构完成:polyQ需要折叠成由四个互锁链组成的zipper形状才能成核,形成该完整steric zipper结构所需的长度正好对应引起神经退行性疾病的临界长度,解释了为何存在明确的致病阈值 自毒效应抑制纤维生长:polyQ倾向于以阻碍生长的方式结合到已形成的核上,这种自毒作用在短链时占主导地位,抑制聚集体的进一步生长和延伸 这种自毒性可表示为核生长的抑制项: \[\dfrac{\mathrm{d}[N]}{\mathrm{d}t} = k_{\mathrm{nuc}}[M] - k_{\mathrm{poison}}[M][N]\] 其中$[N]$为核浓度,$[M]$为单体浓度,$k_{\mathrm{poison}}$为自毒速率常数。在$Q < Q_c$时,自毒项占主导,抑制聚集;在$Q \ge Q_c$时,成核项占主导,聚集加速。47 能量景观的长度依赖转变 能量景观理论的计算预测揭示了Q40附近的热力学相变:来源 短链片段(Q20或Q30)聚集热力学不利:Q20或Q30的Httex1片段grand canonical自由能曲线向上(uphill),表示聚集过程需要克服显著的自由能壁垒,因此在热力学上不利于自发聚集 临界长度片段(Q40)聚集变为自发过程:Q40片段的聚集景观转变为向下(downhill),意味着聚集体比单体更稳定,聚集过程在热力学上成为自发过程,这一临界长度与HD发病的临床阈值高度一致 这为阈值效应提供了热力学基础:Q≥40时,聚集从动力学控制转为热力学驱动。自由能可表示为: \[\Delta G_{\mathrm{agg}}(Q) = \Delta G_0 + \alpha (Q - Q_c)\] 其中$Q_c \approx 37$,当$\Delta G_{\mathrm{agg}} < 0$时聚集downhill。良性与病理长度的核形成自由能差$<1\,\mathrm{kcal\cdot mol^{-1}}$,解释了为何阈值附近的微小长度变化会导致显著的表型差异。48 毒性饱和效应 患者数据的数学建模提示,随polyQ长度增加,polyQ蛋白水平的增加对毒性的贡献递减,表明毒性饱和:来源 \[\text{Toxicity} = T_{\max} \times \dfrac{[\mathrm{polyQ}]^n}{EC_{50}^n + [\mathrm{polyQ}]^n}\] 其中$EC_{50}$随Q长度降低,Hill系数$n$描述陡峭度。该模型解释了为何更长polyQ在较低表达水平即引起毒性。 热力学机制:熵焓平衡与溶剂效应 polyQ聚集的热力学反映了熵-焓平衡的精细博弈: 熵的代价 从无规卷曲形成有序β-sheet导致巨大的构象熵损失($-T\Delta S > 0$) 柔性链变成刚性结构化链,熵损失大致与链长成正比 对于短polyQ,熵损失主导自由能,使聚集不利 焓的增益 在polyQ淀粉样中,每个谷氨酰胺可形成多个氢键(骨架-骨架和侧链-侧链) 高度互插的polar zipper排列释放能量,降低聚集体相对于溶剂化单体的焓 谷氨酰胺侧链的酰胺基团从与水氢键转变为彼此氢键:以水换谷氨酰胺 溶剂熵的贡献 释放有序水分子产生有利的溶剂熵($+\Delta S_{\mathrm{solvent}}$) 进一步补偿构象熵损失 净自由能变化 \[\Delta G_{\mathrm{agg}} = \Delta H_{\mathrm{HB}} - T\Delta S_{\mathrm{conf}} + T\Delta S_{\mathrm{solvent}}\] 其中: $\Delta H_{\mathrm{HB}} < 0$:氢键形成焓增益(有利) $-T\Delta S_{\mathrm{conf}} \gg 0$:构象熵代价(不利) $T\Delta S_{\mathrm{solvent}} > 0$:溶剂熵增益(有利) 关键物理图像 需要达到临界尺寸才能在能量上获益——在此之前,小组装或单链无法补偿熵代价;在此之后,每个添加的单体实际上进一步降低自由能(使更大聚集体更稳定)。这一概念体现在Q40的downhill自由能曲线与Q20的uphill曲线对比中。 3.2.2 从寡聚体到纤维:生长动力学 成核之后,聚集过程进入纤维生长阶段,包括四聚体转化、纤维延伸和二级成核三个关键步骤。这些过程的速率常数已通过多种实验技术精确测定,为定量建模提供了坚实基础。 统一动力学模型与关键速率常数 针对Q35 Httex1的系统研究建立了统一动力学模型,给出三个关键速率常数:49 四聚体转化速率:$k_c = 0.07 \pm 0.01\,\mathrm{h}^{-1}$(描述预成核四聚体向聚集活性态转化) 二级成核速率:$k_s = 0.30 \pm 0.04\,\mathrm{M}^{-1}\mathrm{h}^{-1}$(纤维表面催化新核形成) 纤维延伸速率:$k_{+} = 6.4 \pm 0.6 \times 10^5\,\mathrm{M}^{-1}\mathrm{h}^{-1}$(单体添加到纤维末端) 该模型还包含自毒性剪枝项$k_{\mathrm{poison}}$,描述polyQ单体以生长抑制方式结合到纤维核心。 完整动力学方程组 \(\begin{aligned} \frac{\mathrm{d}[M]}{\mathrm{d}t} &= -k_c[T_4] - 2k_{+}[M][F] - k_s[M][F] - k_{\mathrm{poison}}[M][N], \\ \frac{\mathrm{d}[T_4]}{\mathrm{d}t} &= K_{\mathrm{tetra}}[M]^4 - k_c[T_4], \\ \frac{\mathrm{d}[N]}{\mathrm{d}t} &= k_c[T_4] + k_s[M][F] - k_{\mathrm{poison}}[M][N], \\ \frac{\mathrm{d}[F]}{\mathrm{d}t} &= k_{+}[M][F]. \end{aligned}\) 其中: $[M]$:单体浓度 $[T_4]$:四聚体浓度 $[N]$:成核种子浓度 $[F]$:纤维末端浓度(可延伸末端) $K_{\mathrm{tetra}}$:四聚化平衡常数 方程解释 单体消耗(第一式):单体通过四条路径减少 四聚体转化为活性核($k_c[T_4]$) 纤维两端延伸($2k_{+}[M][F]$,因子2表示两端) 纤维表面二级成核($k_s[M][F]$) 自毒性结合到核上($k_{\mathrm{poison}}[M][N]$) 四聚体动力学(第二式):四聚体由单体快速平衡形成($K_{\mathrm{tetra}}[M]^4$),缓慢转化为活性核($k_c[T_4]$) 成核种子生成(第三式):核由四聚体转化和二级成核产生,被自毒性抑制 纤维生长(第四式):纤维末端数量随延伸增加 Native httex1Q35通过四级一级成核聚集(与预成核四聚化一致),耦合一级二级成核,形成复杂的聚集网络。这些参数可直接用于常微分方程(ODE)或动力学蒙特卡罗(KMC)模拟。来源 N17自组装界面与高阶多聚体 Mishra等(2024)通过丙氨酸扫描和蛋白对接(ClusPro)鉴定了httN17四聚体上的自组装界面。50 该对称界面可能介导四聚体进一步组装为八聚体乃至更高阶多聚体,从而提供可被动力学模型显式表示的寡聚化路径。 该界面把新生polyQ链带到更接近的位置,提高局部有效浓度与空间约束,有利于跨链氢键网络形成并加速聚集。 多个Ala替换会增强螺旋度和/或聚集,提示这些残基可能构成自组装界面,可作为突变验证与配体设计的候选靶点。 已知exon-1聚集抑制剂的对接预测显示,配体可能接触该界面残基,从而通过干扰寡聚体装配来抑制聚集(推断)。 该发现为靶向早期寡聚体提供了结构基础,N17结构域通过促进大型结构稳定寡聚体加速聚集动力学,同时降低纤维异质性。5152 纤维延伸的实时观察 高速原子力显微镜(HS-AFM)首次实现了Htt淀粉样形成中二级成核的单颗粒实时观察:来源 纤维elongation和secondary nucleation路径的直接可视化 纤维延伸显示快速生长与停滞期交替,表明复杂的生长动力学 纤维表面的二级成核事件可被实时追踪 这些观察揭示了纤维生长并非简单的单体逐个添加,而是涉及构象重排、表面催化和多路径竞争的复杂过程。 延伸速率与Q长度的关系 seeding实验的一个重要发现是:纤维延伸速率几乎与polyQ长度无关。 这一观察表明: 长度依赖性主要影响成核步骤,而非已有纤维的生长速率 一旦临界核形成,无论polyQ长度如何,单体添加到纤维末端的速率相似 不同长度的polyQ种子可以以相似速率催化聚集 这与成核时间的指数依赖形成鲜明对比,进一步强调了成核步骤是长度依赖的主要来源。 两步成核机制:从液-液相分离到有序淀粉样 模拟和实验证据支持polyQ聚集的两步机制: 步骤1:液-液相分离(LLPS),详见第四章 polyQ链形成液态样球形寡聚体(浓密簇),这些寡聚体通过疏水相互作用和链间熵驱动自发形成,类似于油水分离中的相分离过程 这些寡聚体是亚稳态的,无定形且无β-结构,内部动力学性质类似液体,分子可以快速交换和重组 在过饱和溶液中类似于相分离成核,形成浓密相作为后续有序化的前体状态 步骤2:有序化转变 在浓密簇内部,缓慢的分子内重排可将某些链转变为β-sheet构象,这一构象转变是速率限制步骤,涉及氢键网络的重组 这一种子可启动有序纤维核的形成,一旦形成稳定的β-sheet核,纤维延伸迅速进行并最终形成成熟的淀粉样纤维 概念上类似于无定形寡聚体先于纤维的实验观察,解释了为何在病理条件下能观察到多种中间态结构 物理图像 \[\text{单体} \xrightarrow{\text{LLPS}} \text{液态寡聚体} \xrightarrow{\text{构象重排}} \text{β-核} \xrightarrow{\text{延伸}} \text{纤维}\] 这种两步机制解释了为何短polyQ寡聚体(如果存在)可能off-pathway或无害——它们无法转化为β-核,最终解离;而超过阈值的寡聚体可转化为毒性β-丰富核。 长度依赖的生长速率 尽管延伸速率本身与Q长度无关,但粗粒化MD模拟显示,Q48比Q23生长速度显著更快,且生长主要沿β-sheet延伸方向,也观察到通过steric zippering生长。53 这与实验观察一致: 更长的polyQ不仅成核更快,纤维延伸也更快 生长速率的长度依赖性可能源于β-sheet稳定性和侧链相互作用强度的增强 Seeded聚集实验证实,Q48种子催化的聚集速率远高于Q23种子 3.2.3 Q长度依赖的动力学模型(补充:多为综述性转述,需逐条回原文复核) Wetzel 的成核–延伸模型(定性趋势) Wetzel 等通过多谷氨酰肽的体外纤维化实验提出成核–延伸框架,强调限速的成核步骤可视作单链进入某种有序构象的平衡,并指出该成核平衡常数会随 $Q$ 增加而上升。54 定量测量:Chen 等的半衰期数据(待核对原文) Chen 等测定不同长度 polyQ 肽的聚集半衰期1:从 $Q20$ 延长到 $Q24$ 时,滞后期由约 $5$ 天缩短至约 $3$ 天,提示少量残基增加即可显著加速。54 Dokholyan 的 $Q\approx 37$ 临界长度模拟(待核对原文) 离散分子动力学等模拟工作常提出3:当 $Q$ 超过某个临界长度(常被写作约 $Q37$)时,单链更可能自发形成稳定的 $\beta$ 结构核(例如平行 $\beta$-螺旋),从而降低成核所需的分子数并提高单分子成核几率。54 短于阈值时:单链难以形成完整有序核,需多个链协同聚集才能成核。 长于阈值时:更长链降低成核所需分子数,并提高单分子成核概率与速率。 Crick:FCS 的构象标度指数(待核对原文) 有实验观测支持“高度塌缩”的单链图像:Crick 等的 FCS 研究报道55 polyQ 肽的水动力半径 $R_h$ 随长度增长的指数约为 $0.32$(小于理想线团的 $0.5$),对应更塌缩的构象。54 更塌缩的单链意味着更易自我关联并进入成核路径,这可解释为何超过临界长度后,某些动力学读数会呈现超线性加速(推断)。54 $k(Q)$ 的经验拟合框架(占位) 在更统一的表述下,可把某个动力学参数写成 $k(Q)$ 的增函数,并用最小可辨识的函数族做拟合,例如指数或幂指数: \[k(Q)=k_0\,\exp(\lambda Q), \quad \text{或} \quad k(Q)=k_0\,Q^{\alpha}.\] 哪一种形式更合理,需要用同一体系下的多长度扫描数据来判别;仅有两三个长度点时,通常无法区分指数与幂律。54 3.2.4 计算模拟方法与多尺度建模 计算模拟在揭示聚集机制、预测长度依赖性和连接分子与临床尺度方面发挥了关键作用。 Multi-eGO多尺度模拟框架 Kulshrestha等(2025)开发的Multi-eGO混合多态结构模型实现了聚集动力学与纤维多态性的统一描述:5152 聚集模拟显示polyQ纤维通过β-turn、β-arc、β-strand组合形成高度异质性形态 N17结构域通过促进大型结构稳定寡聚体加速聚集动力学,同时降低纤维异质性 早期聚集涉及两种机制:骨架相互作用驱动β-sheet形成、侧链交叉指状(interdigitation) 该模型可用于预测不同flanking序列对聚集路径的影响 粗粒化MD的disorder-to-order相变 Dekker等(2025)构建的校准粗粒化MD模型系统探索了从核化生长到液-固相转变的聚集路径:53 通过调节侧链相互作用强度和氢键强度,可覆盖多种聚集机制 Seeded聚集模拟显示,Q48比Q23生长速度显著更快 生长主要沿β-sheet延伸方向,也观察到通过steric zippering生长 模型参数可调,为探索序列变异和更广泛聚集机制提供了通用框架 聚集动力学与临床数据的定量连接 Takahashi等分析了HD患者的polyQ长度与发病年龄(AO)数据,测试多种数学模型:来源 最佳拟合模型:平方和关系 \(t_A^2 = t_N^2 + \Delta t^2\) 其中$t_A$为聚集时间(对应发病年龄),$t_N$为成核时间,$\Delta t$为延伸时间。 物理意义 该模型反映了成核生长聚合分为长度依赖的成核和成核依赖的延伸两个阶段。从成核-延伸动力学可推导: 成核时间呈指数依赖浓度:$t_N \propto [M]^{-n^}$,成核时间强烈依赖于polyQ长度,因为成核阶数$n^$随长度变化(长链$n^* \approx 1$,短链$n^* \approx 4$),这种指数依赖关系使得微小的Q长度变化导致成核时间的数量级差异 延伸时间呈线性依赖浓度:$\Delta t \propto (k_{+}[M])^{-1}$,延伸时间对浓度的依赖相对较弱,呈简单的反比关系,对polyQ长度变化的敏感度远低于成核时间 平方和形式表明成核和延伸是独立的随机过程,总时间为两者的均方根。对于长polyQ($n^* \approx 1$),成核时间主导;对于短polyQ($n^* \approx 4$),成核时间显著延长。 与分子参数的连接 结合前述动力学参数,可建立从分子到临床的定量关系: \(\text{AO}(Q) \approx \sqrt{\left(\dfrac{C_1}{k_n(Q)}\right)^2 + \left(\dfrac{C_2}{k_{+}(Q)}\right)^2}\) 其中$C_1$、$C_2$为与细胞环境相关的常数,$k_n(Q)$和$k_{+}(Q)$分别为长度依赖的成核和延伸速率常数。 其他测试模型(线性、倒数、指数)拟合较差,支持聚集动力学的核化-延伸框架是连接分子机制与临床表型的正确范式。 3.2.5 小结:从分子聚集到疾病发生 PolyQ聚集动力学的研究揭示了多层次的长度依赖阈值效应: Q23-26阈值:成核阶数从四聚体→二聚体→单体的跃迁,体外polyQ肽的物理转折点 Q35-40阈值:β-螺旋几何约束、能量景观uphill→downhill转变、临床致病阈值 自毒性机制:在阈值以下抑制聚集,在阈值以上加速聚集,放大阈值效应 Flanking区域调控:N17自组装界面和PRD构象熵共同调节聚集速率,解释体外与体内阈值的差异 这些机制共同塑造了HD的陡峭S形剂量-反应曲线:在Q35-40附近,微小的长度变化导致发病年龄的巨大差异。 多尺度建模策略 为建立从分子到临床的定量预测,需要整合多尺度数据和模型。 统一建模框架:用多元回归或层级贝叶斯把 $\ln(AO)$ 与 CAG 长度关联,使正常等位基因、somatic expansion、修饰基因(如FAN1、RAI1)作为协变量或随机效应进入同一模型;并把 CAG 长度写入功能读数的状态转移率(如SARA评分、眼科指标、运动里程碑、ALSFRS-R),用半马尔可夫模型捕捉长度到进展速度的差异。 多尺度耦合:将HD分子参数(如 $R_g \propto Q^{\nu}$ 中的 $\nu \approx 0.22$、N17-PRD解离速率常数 $k_{\mathrm{off}}$、四聚体转化速率常数 $k_c$、LLPS临界浓度、膜结合自由能 $\Delta G_{\mathrm{mem}}$)或SBMA的激素依赖项嵌入细胞或组织尺度的ODE或PDE,再用AO或功能读数做数据同化,把体外与体内数据连接到同一参数集。 跨蛋白比较:把不同疾病的回归斜率或hazard归一化成长度灵敏度(如年/重复、hazard per repeat),用于跨蛋白对照、机制对照与靶点优先级排序。 Hazard-聚集耦合模型:把危害函数 $h(t)$ 与聚集概率 $P_{\mathrm{agg}}(t)$ 耦合,例如 $P_{\mathrm{agg}}(t)=1-\exp[-k_{\mathrm{nuc}}(Q)t]$,并用 $k_h$ 作为人群尺度调节参数,以便把体外速率常数映射到人群生存模型。 参数化的落点:聚集动力学的定量参数($k_c = 0.07\,\mathrm{h}^{-1}$、$k_s = 0.30\,\mathrm{M}^{-1}\,\mathrm{h}^{-1}$、$k_{+} = 6.4 \times 10^5\,\mathrm{M}^{-1}\,\mathrm{h}^{-1}$、$k_{\mathrm{poison}}$)可直接作为ODE或KMC的初值或先验,使从分子聚集到临床表型的定量预测更可操作。 3.3 膜相互作用与定位 N17区域因其两亲性螺旋结构可插入脂双层膜,线粒体/ER膜富集度与polyQ长度相关,并影响聚集动力学。35 3.3.1 N17膜结合能 N17的两亲性α-螺旋结构使其能够插入脂双层膜,这种膜相互作用对polyQ聚集有多重影响。 膜结合的Q长度依赖 模型可在能量项中加入膜结合贡献: \[-\Delta G_{\mathrm{mem}} = \sigma (Q - Q_0)\] 其中$\sigma$为膜结合能的Q长度敏感系数(估计值:$0.05-0.1\,\mathrm{kcal\cdot mol^{-1}\cdot Q^{-1}}$),负号表示膜结合降低自由能(热力学有利),$Q_0$为参考长度(通常取23)。35 Nt17稳定化能 N17使polyQ内部有利相互作用减少约$5\,\mathrm{kcal\cdot mol^{-1}}$,降低$\beta$-sheet概率。36 3.3.2 膜富集度的Q长度依赖 Atwal等(2014)通过荧光显微镜定量分析发现,mHtt在线粒体和ER膜上的富集度随polyQ长度线性增加。虽然原文未给出精确的定量函数,但基于实验观察可建立以下关系: \[F_{\mathrm{mem}}(Q) = F_0 + \beta_{\mathrm{mem}} \times (Q - Q_0)\] 其中: $F_{\mathrm{mem}}(Q)$:膜富集度(可用荧光强度比或共定位系数表征) $F_0$:正常长度polyQ(如Q23)的基线膜结合水平 $\beta_{\mathrm{mem}}$:膜富集度对Q长度的敏感系数(单位:富集度/Q) $Q_0$:参考Q长度(通常取23) 实验观察 Q23 Httex1显示低水平的膜定位 Q73 Httex1显示显著增强的线粒体和ER膜富集 膜富集与聚集倾向正相关,提示膜表面可能作为聚集的成核位点 3.3.3 膜定位与聚集的耦合 膜表面可能作为聚集成核的催化界面,形成正反馈循环: \[\text{polyQ} \xrightarrow{\text{N17}} \text{膜定位} \xrightarrow{\text{局部浓缩}} \text{成核加速} \xrightarrow{\text{聚集体生长}} \text{更强膜结合}\] 这种正反馈解释了为何膜富集度与疾病严重程度强相关。定量模型需要耦合膜结合动力学与聚集动力学: \[\dfrac{\mathrm{d}[M_{\mathrm{mem}}]}{\mathrm{d}t} = k_{\mathrm{on}}^{\mathrm{mem}}(Q)[M]_{\mathrm{cyto}} - k_{\mathrm{off}}^{\mathrm{mem}}[M_{\mathrm{mem}}] - k_{\mathrm{nuc}}^{\mathrm{mem}}(Q)[M_{\mathrm{mem}}]^{n^*}\] 其中$[M_{\mathrm{mem}}]$为膜结合的单体浓度,$k_{\mathrm{nuc}}^{\mathrm{mem}}(Q)$为膜表面的成核速率常数(可能比溶液中高数倍)。 3.3.4 膜表面成核动力学 膜表面可能作为聚集成核的催化界面。相对于溶液中的成核,膜表面成核具有以下特点: 成核速率增强 \[k_{\mathrm{nuc}}^{\mathrm{mem}}(Q) = \eta_{\mathrm{mem}} \times k_{\mathrm{nuc}}(Q)\] 其中$\eta_{\mathrm{mem}} > 1$为膜表面成核增强因子(估计值:2-10倍)。 物理机制 二维浓缩效应提高局部有效浓度:膜表面把相关组分限制在二维空间,提高局部有效浓度与相遇频率,从而提升成核概率。 界面诱导的构象偏置:膜界面可能富集或稳定更易成核的构象子集,从而降低有效能垒并提高有效成核速率。 N17锚定带来定位与取向效应:N17插入膜后将polyQ核心定位在膜表面,促进分子间相互作用与取向匹配,加速成核与后续生长。 建模意义 在完整的多尺度模型中,需要同时考虑溶液和膜表面的成核路径: \[\dfrac{\mathrm{d}[N]}{\mathrm{d}t} = k_{\mathrm{nuc}}^{\mathrm{sol}}(Q)[M]_{\mathrm{cyto}}^{n^*} + k_{\mathrm{nuc}}^{\mathrm{mem}}(Q)[M]_{\mathrm{mem}}^{n^*} - k_{\mathrm{poison}}[M][N]\] 3.4 蛋白质降解与清除 细胞通过两条主要途径清除polyQ蛋白:泛素-蛋白酶体系统(UPS)和自噬,包括选择性自噬聚集体的aggrephagy。 3.4.1 双路径降解模型 polyQ蛋白根据其聚集状态被分配到不同降解路径: 泛素-蛋白酶体系统(UPS) 处理可溶性、较小的聚集体 降解速率:$k_{\mathrm{UPS}} \approx 0.1-1\,\mathrm{h}^{-1}$(取决于polyQ长度和聚集体大小) 主要清除单体和小寡聚体 巨自噬/Aggrephagy 清除较大的寡聚体和聚集体(无法进入蛋白酶体) 降解速率:$k_{\mathrm{agg}} \approx 0.01-0.1\,\mathrm{h}^{-1}$(比UPS慢10-100倍) 针对已形成的包涵体和纤维 关键发现:扩展polyQ(Q>35-40)导致蛋白酶体功能受损,形成双稳态开关——当降解能力被overwhelmed时聚集体累积。 3.4.2 Aggrephagy定量模型 2022年Autophagy期刊的研究建立了aggrephagy速率常数的数学模型: \[k_{\mathrm{agg}}(Q, R_{\mathrm{agg}}, \Phi_{\mathrm{auto}}) = k_0 \times f_Q(Q) \times g_R(R_{\mathrm{agg}}) \times h_{\Phi}(\Phi_{\mathrm{auto}})\] 其中: $k_0$:基础aggrephagy速率常数 $f_Q(Q)$:polyQ长度依赖的增强因子(扩展polyQ识别更快) $g_R(R_{\mathrm{agg}})$:聚集体尺寸依赖因子(更大聚集体清除更快) $\Phi_{\mathrm{auto}}$:自噬通量(受mTOR、AMPK等调控) 实验观察 Aggrephagy清除聚集体半衰期:数小时至数天 自噬通量增强10倍可显著减少polyQ累积 扩展polyQ(Q>40)可被选择性识别,清除速率提高2-3倍 3.4.3 分子伴侣协同(Hsp70系统) Hsp70与协同因子(Hsp40、CHIP)在triage决策中起关键作用: \[\text{Substrate} \xrightarrow{\text{Hsp70}} \begin{cases} \text{Refolding} & \text{if } t_{\mathrm{bind}} < t_{\mathrm{crit}} \\ \text{Degradation} & \text{if } t_{\mathrm{bind}} \ge t_{\mathrm{crit}} \end{cases}\] 其中$t_{\mathrm{bind}}$为底物与Hsp70结合时间,$t_{\mathrm{crit}}$为临界时间阈值。 定量参数(来自Hsp70与其他底物的研究,polyQ数据缺失) ATP态亲和力:$K_d^{\mathrm{ATP}} \approx 50-100\,\mu\mathrm{M}$ ADP态亲和力:$K_d^{\mathrm{ADP}} \approx 5-10\,\mu\mathrm{M}$(约10倍增强) ATP水解速率:$k_{\mathrm{hyd}} \approx 0.1-1\,\mathrm{s}^{-1}$ polyQ-Hsp70直接结合的$K_d$值和速率常数仍缺失,阻碍了准确的聚集与清除耦合动力学建模。 3.4.4 双稳态切换与阈值行为 当聚集速率超过清除速率时,系统出现双稳态切换: \[\dfrac{\mathrm{d}[A]}{\mathrm{d}t} = k_{\mathrm{agg}}[M] - (k_{\mathrm{UPS}} + k_{\mathrm{agg}})[A]\] 稳态分析显示: 当 $k_{\mathrm{agg}} < k_{\mathrm{UPS}} + k_{\mathrm{agg}}$ 时,聚集体浓度保持低水平(健康态) 当 $k_{\mathrm{agg}} > k_{\mathrm{UPS}} + k_{\mathrm{agg}}$ 时,聚集体指数累积(疾病态) 扩展polyQ降低$k_{\mathrm{UPS}}$(蛋白酶体捕获),促进切换 这种正反馈循环解释了为何polyQ疾病进展中的突然恶化和不可逆性。 3.5 分子参数速查表 3.5.1 结构与构象参数 参数 数值/关系 来源 N17-PRD解离速率 $k_{\mathrm{off}}$: $20\to 8\,\mathrm{s^{-1}}$(Q25→Q46) 38 回转半径标度 $R_g \propto Q^{\nu}$, $\nu \approx 0.22$ 39 浓度阈值 $c_1, c_2$(HttEx1-17Q: μM级) 42 3.5.2 聚集动力学参数 参数 数值/关系 来源 四聚体转化 $k_c = 0.07 \pm 0.01\,\mathrm{h}^{-1}$ 49 二级成核 $k_s = 0.30 \pm 0.04\,\mathrm{M}^{-1}\mathrm{h}^{-1}$ 49 纤维延伸 $k_{+} = 6.4 \pm 0.6 \times 10^5\,\mathrm{M}^{-1}\mathrm{h}^{-1}$ 49 成核平衡常数 $K_{n^*} = 2.6 \times 10^{-9}$(Q47) 11 成核自由能 $\Delta G_{n^*} \approx +12.2\,\mathrm{kcal\cdot mol^{-1}}$(Q47) 11 3.5.3 膜相互作用参数 参数类别 参数 数值/关系 参考 膜结合能 $\Delta G_{\mathrm{mem}}(Q)$ 建议模型:$\Delta G_{\mathrm{mem}}(Q)=\Delta G_0-\sigma(Q-Q_0)$,$\sigma$ 待拟合 35 膜富集度 $F_{\mathrm{mem}}(Q)$ 建议模型:$F_{\mathrm{mem}}(Q)=F_0+\beta_{\mathrm{mem}}(Q-Q_0)$,$\beta_{\mathrm{mem}}$ 待拟合 35 Nt17稳定化能 $\Delta G_{\mathrm{N17}}$ 约$-5\,\mathrm{kcal\cdot mol^{-1}}$ (待验证) 36 膜表面成核增强 $\eta_{\mathrm{mem}}$ 2-10倍(估计) 推断 3.5.4 其他分子参数(建议移出速查表) 与“长度依赖的分子主线”相比,蛋白质降解与清除、翻译后修饰与抑制剂更适合放在第 5 章集中整理;这些条目往往缺乏跨长度的系统扫描,或只在少数长度点给出趋势,放入速查表容易造成“可跨长度复用”的误解。 蛋白质降解与清除的动力学建模与数据缺口见第 3.4 节与第 5.1 节。 翻译后修饰与蛋白质量控制的具体实例见第 5.1 节。 小分子与抑制剂(SeNPs、姜黄素、共溶质)相关的模型与参数见第 5.2 节。 4. 相分离与凝聚态(LLPS) 4.1 定量理论框架:Flory-Huggins与Cahn-Hilliard 液-固相转变 Httex1(Q43)在生理盐缓冲液下的LLPS(liquid-liquid phase separation,液-液相分离,形成无膜液滴凝聚物)到固化(gelation/solidification)可由Cahn-Hilliard方程(描述相分离动力学的偏微分方程)描述,实验给出液滴黏度从初始的2.9 Pa·s 增至 17 Pa·s(约6小时),表明液滴逐渐老化、硬化。关键临界浓度(critical concentration,相分离的最低蛋白浓度阈值)随 Q 升高而下降,意味着更长的polyQ更容易发生相分离。56 Flory-Huggins理论的Q长度依赖预言 经典的高分子热力学模型——Flory–Huggins溶液理论——提供了定量框架,将聚谷氨酰胺链长Q作为参数来预测相分离行为。在该模型中,高分子链长度$N$(在此相当于PolyQ序列的Q数)会显著影响相分离的临界条件。54 Starikov等(1999):Flory–Huggins与长度阈值 Starikov 等在纯 polyQ–水的简化体系中,使用 Flory–Huggins 平均场晶格模型讨论随机线团→$\beta$-hairpin($\beta$-sheet-hairpin)转变。其混合的单位格点 Gibbs 自由能(每个残基)写为:57 \[\tilde{G}(\phi,N) = RT\left[ \frac{\phi}{N}\ln\phi +(1-\phi)\ln(1-\phi) -X\,\phi(1-\phi) \right].\] 其中,$\phi$ 为聚合物体积分数,$N$ 为链长(残基数),$X$ 为 Flory 参数(其符号约定与常见 $\chi$ 写法不同,因此这里保持作者记号)。该工作把 hairpin 的形成近似为链长从 $N$ 有效变为 $N/2$,并定义转变的自由能差:57 \[\Delta\tilde{G}_{c\to h}(N)=\tilde{G}(\phi_c,N/2)-\tilde{G}(\phi_c,N).\] 作者在差溶剂(文中采用 $X\approx 1.133$ 的 glutamine–water 估计)下,引用聚合物溶液理论指出存在临界体积分数 $\phi_c$,其随链长呈 $\phi_c\propto N^{-1/2}$ 的标度(原文在该处写作 $\phi_c\sim \cdots/\sqrt{N}$)。57 在达到该 $\phi_c$ 后,作者评估 $\Delta\tilde{G}{c\to h}(N)$ 随 $N$ 的变化,并得到一个与链长 $Q$ 直接相关的定量判据:当 $N$(可视为 $Q$)增大到约 $40$ 左右时,$\Delta\tilde{G}{c\to h}(N)$ 从正变负,提示随机线团对 hairpin($\beta$-sheet)构象出现热力学不稳定。57 \[Q_c \approx N_0 \approx 40, \quad \Delta\tilde{G}_{c\to h}(Q_c)=0.\] 因此,在该简化模型下,$Q<Q_c$ 对应 $\Delta\tilde{G}{c\to h}(Q)>0$(不自发),而 $Q>Q_c$ 对应 $\Delta\tilde{G}{c\to h}(Q)<0$(倾向自发形成 hairpin),这就是其长度阈值预测的数学形式。57 公式的通俗解释 该模型的核心是把构象转变等效成链长变化,然后比较在同一溶液条件($\phi=\phi_c$)下的自由能: 若 $\Delta\tilde{G}_{c\to h}(N)>0$,说明从随机线团变成 hairpin 在热力学上不划算,不会自发发生; 若 $\Delta\tilde{G}_{c\to h}(N)<0$,说明 hairpin 更有利,随机线团会倾向于转向 $\beta$-相关构象。 因此,$N\approx 40$ 附近的变号给出一种早期的、基于聚合物热力学的解释:为什么 polyQ 体系会在 $40$ 左右出现显著加速或类似阈值的现象。但作者也明确指出,纯 polyQ–水体系的临界浓度对溶剂组成高度敏感,因此该模型更适合用来理解趋势与尺度,而不是逐点定量预测。57 分子力学的补充:polyQ40 的两类折叠模式(力场依赖) 同一篇工作还对 polyQ40 做了分子力学构象搜索,发现两类拓扑上不同的折叠模式:57 在 CHARMM 力场下,polyQ40 倾向形成较经典的 $\beta$-hairpin,并观察到约 $57$ 个氢键;除主链间氢键外,侧链间氢键也大量出现并可能稳定 $\beta$-sheet 的褶皱与弯曲。 在 AMBER/OPLSAA/TINKER 等力场下,更倾向得到高度紧密的随机线团,并且氢键数量更高。 该结果提示:polyQ40 的折叠基序可能具有明显的力场敏感性与构象多形性。作者据此提出两个临界浓度的图景($C_1<C_2$):先发生随机线团→hairpin,再在更高浓度下发生进一步病理性折叠或致密化,随后进入成核与聚集(文中称 in-quarto 折叠)。这一部分更像机制性假说,需与现代 Httex1 的结构学与动力学数据对照后再用于定量建模。57 Pappu课题组Crick在其博士论文(2011)中首次系统测定了polyQ溶液的相分离参数。他们通过测量聚谷氨酰胺肽(30Q vs 40Q)在水溶液中的饱和浓度随温度变化,构建了polyQ溶液的参考相图。Flory-Huggins拟合清晰显示:40Q肽的饱和浓度曲线显著低于30Q肽,意味着链长增加显著降低了临界浓度、提高了聚集/相分离倾向。正如作者指出,固定其他条件时,polyQ链越长,驱动自组装和相分离的内作用越强。58 Zeng等(2020)提出了单链塌缩–相分离相当性框架,将单链线团-球体转变参数映射到相分离相图。他们利用模拟提取了polyQ的Θ温度$T_θ$和二阶、三阶维里系数,代入改进的Flory-Huggins随机相近似(RPA)计算相图。测试发现,随着Q从低到高延长,临界温度$T_c$单调升高,而临界体积分数$\phi_c$显著降低(例如70Q的$T_c$比短序列更高,$\phi_c$更低)。这定量验证了Flory-Huggins理论预言:polyQ长度与相分离参数呈系统相关。58 Q长度依赖的定性趋势:虽然Peskett等(2018)明确指出临界浓度会随Q升高而下降,但目前文献中缺乏系统的定量测量。Suarez等(2022)研究了mHttex1的浓度依赖相转变,发现M/S/F三相(单体/球形寡聚体/纤维)存在尖锐浓度边界,并被profilin等配体调控,但同样未提供不同Q长度的临界浓度数值。5614 为了把是否相分离与相分离后如何演化写成可计算模型,通常需要一个平衡态自由能与一个动力学方程: Flory–Huggins 自由能(聚合物溶液的最简模型) 令 $\phi$ 为聚合物体积分数,$N$ 为聚合度(链长),$\chi$ 为 Flory–Huggins 相互作用参数,则混合自由能密度可写为: \[\frac{f(\phi)}{k_B T} = \frac{\phi}{N}\ln\phi + (1-\phi)\ln(1-\phi) + \chi\,\phi(1-\phi)\] 公式的通俗解释 该式把相分离驱动力拆成三部分:链长相关的聚合物熵项、溶剂熵项与相互作用项。对于固定 $\chi$,链越长($N$ 越大),熵惩罚越小,因此更容易进入两相区并发生 LLPS。 补充:$\chi$ 是 Flory–Huggins 相互作用参数(无量纲),用于把大量微观相互作用(氢键、疏水、静电等)压缩成一个平均场强度。直觉上,$\chi$ 越大表示聚合物–溶剂越“不相容”,链更倾向彼此吸引、更容易塌缩或分相;$\chi$ 越小表示更相容、更倾向均匀混合。在简单近似与长链极限下,$\theta$ 条件常对应 $\chi\approx 1/2$。 临界点条件(从“混合”到“自发分相”) 在上述模型中,临界点满足: \[\frac{\partial^2 f}{\partial \phi^2}=0,\qquad \frac{\partial^3 f}{\partial \phi^3}=0\] 由此可得到经典结果: \[\phi_c = \frac{1}{1+\sqrt{N}},\qquad \chi_c = \frac{1}{2}\left(1+\frac{1}{\sqrt{N}}\right)^2\] 公式的通俗解释 这里的 $\chi$(Flory–Huggins 相互作用参数)是一个无量纲的平均场参数,用来概括“聚合物链段–溶剂”的有效相互作用强度。直觉上,$\chi$ 越大表示聚合物与溶剂越“不相容”,链更倾向彼此吸引,因此更易塌缩或发生相分离;$\chi$ 越小表示更相容,更倾向均匀混合。在简单近似与长链极限下,$\theta$ 条件常对应 $\chi\approx 1/2$。 这些公式揭示了polyQ长度对相分离的核心影响: 随着PolyQ长度$N$(即Q数)增加,发生液液相分离所需的相互作用强度 $\chi_{cr}$ 逐渐接近 $\theta$ 条件($\chi=0.5$) 相分离所需的聚合物临界浓度会显著降低,近似按 $Q^{-1/2}$ 的规律降低 这意味着更长的PolyQ链在较低浓度下就更容易发生相分离,从而具备更高的LLPS倾向 理论意义 Starikov等人基于Flory-Huggins建立的简单模型就明确以PolyQ长度为自变量,预测当Q序列超过约40时无规线团更易失稳并偏向β折叠聚集核。这一理论阈值(约Q≈40)与体外聚集实验和临床发病阈值非常吻合。54 Cahn–Hilliard 动力学(描述液滴生成、长大与熟化) 若用 $\phi(\mathbf{r},t)$ 表示空间与时间依赖的组分场,常用 Cahn–Hilliard 方程描述保守型相分离动力学: \[\frac{\partial \phi}{\partial t} = \nabla\cdot\left(M\nabla\mu\right),\qquad \mu=\frac{\partial f}{\partial \phi}-\kappa\nabla^2\phi\] 其中 $M$ 为迁移率,$\mu$ 为化学势,$\kappa$ 为梯度能系数。 公式的通俗解释 这组方程回答两个问题:什么时候开始分相、以及分相以后结构如何演化。 第一式是“质量守恒扩散方程”,右侧写成 $\nabla\cdot(M\nabla\mu)$ 表示物质沿化学势梯度扩散;$M$ 越大,相分离和粗化越快。 第二式给出化学势 $\mu$ 的来源:$\partial f/\partial\phi$ 决定局部混合是否稳定;$-\kappa\nabla^2\phi$ 惩罚过陡的界面,等价于赋予界面有限厚度与表面张力。 更“化学直觉”的理解是:$\mu$ 就像“把 1 个分子从当前位置搬走所付出的自由能代价”。如果某个地方 $\mu$ 高,体系就倾向于把该处的分子搬走,流向 $\mu$ 低的地方,让整体自由能下降。 之所以不是“沿浓度梯度扩散”,而是“沿化学势梯度扩散”,是因为在相分离体系里,驱动力不只是“哪边更浓”,还包括“哪边更不稳定、哪边界面代价更大”。$\mu$ 把这些因素合并成一个驱动力量。 $\kappa$ 可以理解为“界面惩罚系数”:$\kappa$ 越大,体系越不喜欢出现很尖锐的界面,因此液滴边界会更平滑、更厚,并在动力学上更倾向于减少界面总面积(这对应化学上常说的“表面张力驱动液滴并并/粗化”)。 在扩散控制的粗化阶段,典型液滴尺度 $R(t)$ 常满足近似的幂律增长 $R(t)\propto t^{1/3}$,对应 Ostwald ripening 或扩散控制的相域粗化,能用来把显微尺度的液滴尺寸随时间变化映射到参数 $M$ 与自由能曲面。 对Httex1这类体系,液滴老化→固化往往意味着 $M$ 随时间下降,或体系出现粘弹性与不可逆键合,从而偏离简单的 $t^{1/3}$ 粗化并转入凝胶化/玻璃化动力学,与实验观察到的黏度随小时尺度上升相一致。56 Q长度依赖的理论参数:当前缺口 尽管Flory-Huggins和Cahn-Hilliard理论为LLPS提供了成熟的理论框架,但将这些理论参数与Httex1的Q长度定量关联的系统研究仍然缺失。 已知的定性趋势 临界浓度$C_{\mathrm{sat}}$随Q长度增加而下降(Peskett 2018)56 更长polyQ更容易相分离,与Flory-Huggins预测的长链降低$\phi_c$一致 液滴老化速率可能与Q长度相关(Q43的黏度6小时内增加6倍) 缺失的定量参数 Flory-Huggins相互作用参数:$\chi(Q)$与Q长度的函数关系未知。理论预测$\chi_c \to 1/2$(长链极限),但实际polyQ的$\chi$值及其Q依赖性未测量。 临界浓度的定量关系:虽知$C_{\mathrm{sat}}(Q)$递减,但缺乏具体数值。例如: Q25的$C_{\mathrm{sat}}$是多少? Q35、Q46、Q60的$C_{\mathrm{sat}}$如何变化? 是否遵循指数、幂律或其他函数形式? Cahn-Hilliard动力学参数: 迁移率$M(Q)$:液滴粗化速率与Q长度的关系 梯度能系数$\kappa(Q)$:界面张力与Q长度的关系 这些参数可从液滴尺寸随时间演化$R(t)$拟合,但缺乏系统测量 相图数据:binodal(共存曲线)和spinodal(失稳曲线)的Q长度依赖未绘制。完整相图需要测量不同Q长度在多个温度和浓度下的相边界。 为何缺乏这些数据? 实验挑战包括:polyQ蛋白聚集快速,难以达到热力学平衡;液滴快速老化固化,偏离经典LLPS假设;Httex1含N17和PRD flanking区域,非简单均聚物;聚集与相分离过程耦合,难以分离。理论上,可能需要超越Flory-Huggins的框架,考虑多组分、非平衡和构象转变耦合效应。 建模策略 在缺乏精确参数的情况下,可采用以下近似: 使用Flory-Huggins框架的定性预测:$\phi_c \propto N^{-1/2}$,$C_{\mathrm{sat}}(Q) \propto Q^{-\alpha}$($\alpha \approx 0.5$) 从已知的Q43数据(黏度、时间尺度)外推其他Q长度 结合聚集动力学参数($k_c$、$k_s$、$k_{+}$)与LLPS模型,构建耦合方程 使用粗粒化模拟(如Baidya 2025的方法)预测$\chi(Q)$和$c_\theta(Q)$ 4.2 实验验证的定量关系 Posey等(2018)的三相模型 Posey等利用纯化的Htt N端片段进行体外逐步饱和实验59,观测到Htt-NTF在体外可存在三种相态,每一相对应一个特征饱和浓度范围。58 M相(Monomer/Oligomer相):可溶性单体/低聚相 S相(Solution相):液滴状浓缩相,动态可逆 F相(Fiber相):不可溶纤维相,固态聚集体 对于致病扩展长度的polyQ片段,研究者测定了出现液滴相的起始浓度$c_S$以及出现纤维沉淀的浓度$c_F$。这些相边界提供了类似binodal线的信息。更重要的是,他们发现N17片段能够显著降低纤维相的饱和浓度$c_F$,即加入N17提高了聚集倾向,使纤维相更早出现;反之,去除N17则需要更高浓度才出现固相。58 Peskett等(2018)的液-固转变研究 Peskett等在哺乳细胞和酵母中对比了Httex1-25Q、43Q和97Q的行为,发现: 所有长度(包括25Q)均能形成液滴状的动态聚集体 但仅有含扩展polyQ的43Q和97Q进一步经历液-固转变,生成稳定的固样包涵体 25Q的液滴保持液态且可逆,而扩增的polyQ使体系越过某临界点发生固化 Krobitsch和Lindquist的早期研究也支持这一观察:野生型Htt(Q~25)在酵母中不形成可见包涵体,而超长polyQ(如Q>70)才会形成聚集斑。近期的工作将这一差异解释为相分离相图的改变:短Q的Httex1即使浓度升高也多以可溶/液态存在,而长Q的Httex1更容易跨过相界限形成浓缩相甚至固相。58 Liu等(2025)Ataxin-2研究的定量观察 Liu等(2025)以Ataxin-2蛋白的N端片段为模型,系统研究了polyQ扩展如何调控液滴性质。他们构建了Ataxin-2的两个片段(含N端polyQ且缺失C端结构域),分别引入不同长度的Q扩展(9Q、23Q、33Q、41Q),在细胞和体外体系中观察相分离行为。60 液滴形成能力 荧光显微和颗粒计数显示,随着polyQ重复数增加,形成的液滴数量显著增加且平均液滴尺寸也变大 特别是在体外5 µM蛋白浓度+PEG诱导条件下,9Q和23Q样品仅产生少量小液滴,而33Q和41Q样品产生大量且更大的液滴 研究人员据此绘制了液滴数量和直径随Q长度的关系曲线,呈现出在Q约23–33之间液滴形成能力出现跃变 这一转折点与Ataxin-2相关疾病的中间扩增区间(ALS相关的中等扩增>23Q,SCA2致病阈值>33Q)相符 液滴内部流动性(FRAP实验) 短Q(如Atx2-9Q、23Q)液滴表现出较高的动态性,在漂白后数分钟内可恢复约60%的荧光 长Q(33Q、41Q)液滴几乎不恢复,只有不到10%信号回升,表明后者内部分子几乎被固化 这一差异明确了polyQ扩增使凝聚相由液态转向固态的趋势: 短polyQ液滴多保持液态流动性,相分离可逆,分子扩散与交换相对更快 长polyQ液滴更易凝胶化或固化,内部分子扩散受限,提示凝聚态内部可能已形成高阶网络 两步成核机制的热力学依据 Crick等利用Flory-Huggins拟合得到的相图计算了不稳定边界(spinodal)和饱和曲线(binodal)之间的间隔,发现在生理温度附近,两者间浓度跨度可达约两位数量级。这意味着对于一定浓度范围,体系处于热力学亚稳区:单链已倾向塌缩但整体尚未自发相分离。在此区域内,体系可能通过形核形成亚稳寡聚液滴,再经构象重排走向不可逆聚集。55 这一分析为两步成核机制提供了热力学依据: polyQ链先经历液液相分离形成浓缩液滴(由Flory-Huggins自由能面预测的亚稳相) 随后在液滴内部或表面发生转变生成固态聚集体 通过将Q长度纳入Flory-Huggins模型,研究者不仅定量比较了不同Q长度下的相图差异,还推演了聚合驱动力与链长的定量关系:链长增加相当于有效相互作用参数$\chi$增大或临界$\chi_c$降低,使体系更易进入相分离甚至聚集区域。55 从动力学角度看,长polyQ导致的液滴内部固化现象(FRAP不可逆)实际反映了凝聚态结构形成。Liu等通过Ataxin-2的FRAP比较指出,polyQ扩增使液滴分子扩散受阻,暗示链在液滴中可能已经部分堆积为不溶原纤维或高阶网络。60 Ataxin-2 polyQ扩展错误隔离TDP-43 Wijegunawardana等(2024)报告,Ataxin-2 polyQ扩展错误地隔离TDP-43于RNP凝聚物内,破坏其沿轴突的运动性和液体样性质。Ataxin-2控制神经元RNP凝聚物的运动性和翻译,polyQ扩展从根本上扰乱mRNA空间定位并抑制局部翻译。该研究支持一个模型:Ataxin-2 polyQ扩展破坏关键轴突和细胞骨架mRNA的稳定性、定位和/或翻译,对运动神经元完整性特别重要。为polyQ毒性的RNA翻译调控机制提供了定量框架。61 FOXP2:PolyQ长度依赖的LLPS与功能调控 FOXP2转录因子天然含有两个长PolyQ序列,其不同物种的长度变异影响FOXP2在细胞核内形成凝聚体的性质,为理解PolyQ长度对LLPS的影响提供了另一个重要案例。54 研究发现 FOXP2的PolyQ区可形成coiled-coil结构驱动LLPS 不同长度组合改变凝聚体液态与纤维态平衡以及转录活性 PolyQ长度作为分子开关调控FOXP2的相分离行为 生物学意义 这一发现表明PolyQ扩展不仅与疾病相关,在正常生理条件下也可作为调控蛋白相分离和功能的分子开关。FOXP2通过PolyQ长度变化调控其LLPS行为和转录活性,为理解PolyQ序列的生理功能提供了新视角。54 理论框架 这些案例共同支持从序列长度预测相图的框架。在已知环境条件下,通过输入序列中的Q数,模型可估算其相图参数(如$\chi(Q)$、$\phi_c(Q)$、$T_c(Q)$)以及聚集路径(如成核方式、速率)。54 4.3 内在刚度与溶剂效应 Baidya等(2025):θ溶剂区的长度依赖研究 Baidya等(2025)通过计算机模拟系统研究了具有极性侧链(polyQ)和疏水侧链(polyL)的IDPs在不同共溶质浓度下的θ溶剂区。由于内在刚度,这些IDPs在短长度尺度上总是扩展的,与溶剂质量无关。因此,对于短IDP序列(<25残基),其LLPS倾向无法从单链性质推断。进一步,对于有限尺寸IDPs,从结构因子(模拟SAXS)和配对距离(模拟smFRET)提取的θ溶剂共溶质浓度 $c_\theta$ 不同,仅在大 $N$ 时收敛。研究表明,θ溶剂区的回转半径满足标度关系 $R_g(N) \propto N^{\nu_\theta}$,其中 $\nu_\theta \approx 0.5$,可用于准确提取 $c_\theta$。该研究强调有限尺寸校正在分析IDP性质以识别θ溶剂区时的重要性。62 公式解释 θ溶剂区是刚好不好不坏的溶剂条件,IDP既不过度收缩也不过度扩展。$\nu_\theta \approx 0.5$ 是理想链标度指数,但短链因刚性而偏离。SAXS和smFRET测到的 $c_\theta$ 只在长链时一致,提示短polyQ分析需谨慎。 θ 溶剂区是什么 θ 溶剂区可以理解为排斥体积效应被抵消的条件,单链统计接近理想链,等价表述包括 $B_2 \approx 0$(其中 $B_2$ 是第二维里系数;$B_2>0$ 对应好溶剂使链更伸展,$B_2<0$ 对应差溶剂使链更收缩)。在 Flory–Huggins 框架下,长链极限的 θ 条件常对应 $\chi_\theta \approx 1/2$。 因此,Baidya 等把共溶质浓度 $c_\theta$ 视为把体系从好溶剂调到 θ 条件的控制参量。对 polyQ 这类具有内在氢键内聚的极性同聚物来说,$c_\theta$ 与链长、刚度都会共同影响 LLPS 的出现与临界浓度(推断)。63 Murphy等(2009)的Q长度依赖溶剂质量实验 Murphy等(J. Mol. Biol. 2009)系统考察了Q8~Q24多聚谷氨酰胺肽的构象和聚集性质,为理解Q长度依赖的溶剂质量转变提供了直接实验证据。54 短链表现为良好溶剂条件(Q8、Q12):链较为扩展,水分子对polyQ链的溶剂化作用强,链内氢键被水分子竞争抑制,链呈现扩展的随机线团构象,此时排斥体积效应占主导 临界长度接近Θ溶剂(Q16):Q16肽在37℃下尚能保持单体态,被认为是良溶剂和差溶剂的分界长度,链内与链外相互作用达到微妙平衡,为链塌缩的临界点 长链体现差溶剂条件(Q20、Q24):链明显塌缩,链内氢键逐渐占据主导地位,链间疏水相互作用增强,驱动链采取更紧凑的构象,Q20和Q24肽一稀释即自发形成可溶性寡聚体(液滴态前体),而更短的Q8、Q12则不出现聚集 研究通过链构象半径的测定进一步量化了这一趋势:随Q增长,多聚Q肽的末端距离显著减小,有效刚度(持续长度)从约11Å降至7Å,表示链逐渐塌缩。54 理想链公式 在理想$\Theta$溶剂下,PolyQ末端距满足公式: \(R_{\Theta} = 5.7\sqrt{n+1}\AA\) 其中$n$为Gln数。54 与理论对比 结果显示: Q8、Q12的实际尺寸比理想链更扩展($R>R_{\Theta}$,说明溶剂偏好) Q20、Q24则更紧凑($R<R_{\Theta}$,溶剂不良) 这证明随着PolyQ长度增加,溶剂对链的有效良溶性下降,链间相互吸引占优,定量上相当于Flory-Huggins的相互作用参数$\chi$随Q升高而增大、接近并超过相分离阈值0.5。54 理论意义 这些结果与Flory-Huggins理论预言的$\phi_{cr}(Q)$下降趋势一致。简而言之,PolyQ链越长,体系越接近相分离条件:临界浓度$\phi_c$随Q增加而降低,临界温度$T_c$则升高趋近$\Theta$温度。54 4.4 其他细胞因子对LLPS的影响 Rad23B异型缓冲延迟Ataxin-3相变 Prasad等(2025)发现,Rad23B(蛋白酶体凝聚物主要成分)与Ataxin-3的异型相互作用抑制Ataxin-3液滴成熟,但不抑制稀释条件下的淀粉样形成。表明Ataxin-3通过错误折叠路径的聚集不同于凝聚路径。Ataxin-3在arsenite应激下被整合到液体样应激颗粒中。该发现为polyQ蛋白聚集动力学与相分离的解耦提供了证据,提示在模型中应分别处理两种路径: \[\dfrac{\mathrm{d}[\mathrm{Agg}]}{\mathrm{d}t} = k_{\mathrm{misfold}}[M] - k_{\mathrm{buffer}}[M][\mathrm{Rad23B}], \quad \dfrac{\mathrm{d}[\mathrm{LLPS}]}{\mathrm{d}t} = k_{\mathrm{phase}}[M]\] 其中异型缓冲仅影响LLPS成熟而非淀粉样形成。64 公式解释 该式区分两条聚集路径:错误折叠形成淀粉样($k_{\mathrm{misfold}}$)和液-液相分离形成液滴($k_{\mathrm{phase}}$)。Rad23B通过“异型缓冲”减缓液滴成熟($k_{\mathrm{buffer}}$ 项),但不影响淀粉样路径。意味着相分离与纤维化可独立调控。 4.4.1 氨基酸调控 PNAS 2024年12月研究发现,三种特异性氨基酸(proline、glutamine、glycine)显著抑制应激颗粒形成,提供了氨基酸调控LLPS的定量框架: 这些氨基酸在应激条件下会整体降低应激颗粒的相分离倾向,可能通过改变局部溶剂化与弱相互作用网络来削弱“粘连”驱动力。 这一结论对 polyQ 特别关键:polyQ 本身富含 glutamine,而 glutamine 在该研究中被归为具有抑制应激颗粒形成效应的氨基酸之一,提示“polyQ 一定促进 LLPS”并非对所有体系都成立。 该工作强调氨基酸组成可以细粒度调控凝聚态相行为,因此在跨蛋白、跨疾病建模时,应把“组成差异”作为与长度同等重要的变量来考虑。 该发现为6.3节相分离机制提供了氨基酸水平的调控维度。来源 4.4.2 细胞拥挤效应 细胞内环境的大分子拥挤(macromolecular crowding)显著影响polyQ聚集动力学,是理解体外与体内差异的关键因素。 拥挤加速聚集的定量效应 实验研究表明,生理水平的拥挤剂(20-30% w/v PEG、Ficoll、dextran)使polyQ聚集加速3-5倍: 参数 稀释缓冲液 拥挤/细胞环境 变化倍数 滞后期 24-48小时 6-12小时 2-4倍缩短 成核速率常数 $k_n$ 基线 4-7倍增强 4-7倍 延伸速率常数 $k_{+}$ 基线 2-3倍增强 2-3倍 二级成核 $k_2$ 基线 5-10倍增强 5-10倍 临界核尺寸 6-7单体 3-4单体 ~50%减少 临界浓度 $C_{\mathrm{crit}}$ 基线 3-5倍降低 降低65-80% 物理机制 排空体积效应(Excluded Volume Effect) 拥挤剂通过减少可用体积来提高“有效浓度”,从而把体系更容易推入过饱和区并加速成核与粗化过程。 该效应可用有效浓度增强因子 $(1-\phi)^{-1}$ 进行近似,其中 $\phi$ 为拥挤剂体积分数。 在典型细胞条件($\phi=0.2$–$0.3$)下,该因子意味着约 $25$–$40$% 的有效浓度提升。 排空吸引力(Depletion Attraction) 由 Asakura–Oosawa 理论描述的“排空吸引力”会给蛋白–蛋白接触额外提供热力学驱动力,从而降低成核所需的有效自由能势垒并促进聚集。 该机制在量级上可表现为把成核相关势垒降低约 $3$–$5\,\mathrm{kJ/mol}$(具体数值高度依赖体系与拥挤剂)。 它会促进蛋白–蛋白形成更稳定的接触并提高有效碰撞成功率,从而缩短滞后期并加速后续增长。 软相互作用 polyQ 与拥挤剂之间还可能存在非特异性的“软相互作用”(例如弱吸附或溶剂化改变),这会让拥挤效应不再只是几何排斥,从而在不同拥挤剂体系间引入显著差异。 修正的成核-延伸模型 拥挤环境下的速率常数可表示为: \[k_{\mathrm{crowded}} = k_0 \times \exp\left(-\dfrac{\Delta G_{\mathrm{crowding}}}{RT}\right)\] 其中$\Delta G_{\mathrm{crowding}} = \Delta G_{\mathrm{excluded}} + \Delta G_{\mathrm{soft}} + \Delta G_{\mathrm{depletion}}$ 或使用修正的Finke-Watzky两步模型: \[\dfrac{\mathrm{d}[P]}{\mathrm{d}t} = k_n \times f_{\mathrm{crowd}}(\phi) \times [M]^n - k_{-} \times [P] \\ \dfrac{\mathrm{d}[A]}{\mathrm{d}t} = k_{+} \times f_{\mathrm{crowd}}(\phi) \times [P] \times [M]\] 其中拥挤增强函数 $f_{\mathrm{crowd}}(\phi) = \exp\left[\dfrac{\alpha\phi}{1-\phi}\right]$。 关键结论 对于Q>40的病理性polyQ,模型预测细胞内聚集比稀释缓冲液中快10-15倍,这与实验观察一致。这解释了为何polyQ疾病在体内的发病速度远超体外实验的预测。 4.4.3 细胞内外聚集差异 定量研究揭示了细胞内与体外polyQ聚集的显著差异。 浓度依赖的滞后相 体外实验 体外动力学常见的经验规律是:滞后时间满足 $t_{\mathrm{lag}}\propto[M]^{-n^{}}$,其中 $[M]$ 为单体浓度,$n^{}$ 为临界核尺寸,因此滞后期对浓度非常敏感。 当把蛋白总浓度提高 $10$ 倍时,滞后期往往会出现“倍数级”缩短,这一现象可用成核项对浓度的幂次依赖来解释。 因此,体外聚集曲线通常可以被归入“成核依赖的聚合动力学”,并可作为把分子尺度速率映射到细胞尺度有效时间的一个入口变量。 细胞内测量使用先进技术: 荧光相关光谱(FCS)用于估算局部 $\mathrm{HTT}$(或 Httex1)的有效浓度,从而把“细胞内真实浓度”与体外配制浓度对齐。 荧光恢复后漂白(FRAP)用于评估液滴或聚集体内部的交换与流动性,进而区分“可逆液态凝聚”与“逐步固化的凝胶/固态状态”。 单分子成像用于捕捉稀有的成核事件与早期寡聚体出现的时间窗,从而约束成核率而不是只拟合宏观终点。 定量对比 细胞内的有效蛋白浓度通常比培养液条件下的体外稀释体系更高,常见量级为约 $2$–$3$ 倍(该比值依细胞类型与定位而变)。 在把拥挤效应纳入后,细胞内聚集速率常表现为比体外更快的时间尺度,典型差异可达约 $3$–$5$ 倍(不同体系差异很大,需逐条核对原始数据)。 从相分离视角看,细胞内的有效饱和浓度或“临界浓度”往往更低,常被报告为可降低约 $3$–$5$ 倍,这会把体系推入更易凝聚的相区。 建立从体外到体内的映射需考虑 需要把“有效浓度增强”显式写入模型(例如用拥挤因子修正 $[M]$ 或修正 $k_{\mathrm{nuc}}$),否则体外速率常数难以外推到细胞内时间尺度。 需要把“局部微环境”作为空间异质性处理(例如细胞器表面或应激颗粒导致的局部浓缩),因为成核对局部浓度的幂次依赖会放大这种异质性。 需要把蛋白质质量控制系统视作与聚集竞争的动力学通路(伴侣蛋白、UPS、自噬等),否则难以解释同一长度在不同细胞类型中的巨大差异。 4.5 LLPS参数速查表 4.5.1 细胞参数速查表 参数类别 参数 数值/关系 参考 细胞拥挤 拥挤加速 聚集加速3-5倍 WebSearch 细胞拥挤 滞后期缩短 从24-48h降至6-12h WebSearch 细胞拥挤 临界浓度降低 3-5倍降低 WebSearch 细胞拥挤 成核速率增强 4-7倍 WebSearch 轴突运输速率 $v_0$ 0.5-2 μm/s(正常) 文献通用值 运输缺陷 $v_{\mathrm{eff}}(Q)$ 随Q指数下降(定量数据缺失) 35 膜表面成核 $k_{\mathrm{nuc}}^{\mathrm{mem}}$ 可能比溶液高数倍(未测量) 推断 注:膜结合热力学参数、膜富集度的Q长度依赖公式见第3.5.3节。 4.5.2 蛋白质降解参数 蛋白质降解与清除的定量参数(UPS速率、aggrephagy速率、Hsp70结合常数等)见第3.5.4节。 4.5.3 细胞毒性参数 参数类别 参数 数值/关系 参考 PC12 Q74 第4天活率 约84% 65 PC12 Q74 第6天活率 约75% 65 PC12 Q74 第8天活率 约60% 65 PC12 Q23 活率维持 >95% 65 Caspase抑制 zVAD-fmk 部分挽救 65 4.5.4 LLPS与相分离参数 参数类别 参数 数值/关系 参考 LLPS 黏度增幅(Q43) 2.9 → 17 Pa·s (约6小时) 56 LLPS 临界浓度趋势 $C_{\mathrm{sat}}$随Q下降(定性) 56 LLPS θ溶剂标度 $\nu_\theta \approx 0.5$ 63 LLPS 氨基酸抑制 Pro, Gln, Gly抑制应激颗粒 WebSearch LLPS理论缺口 $\chi(Q)$ Flory-Huggins参数未测量 - LLPS理论缺口 $C_{\mathrm{sat}}(Q)$ 定量函数关系未知 - LLPS理论缺口 $M(Q)$, $\kappa(Q)$ Cahn-Hilliard参数未测量 - 4.5.5 膜相互作用参数 参数类别 参数 数值/关系 参考 轴突运输速率 $v_0$ 0.5-2 μm/s(正常) 文献通用值 运输缺陷 $v_{\mathrm{eff}}(Q)$ 随Q指数下降(定量数据缺失) 35 膜表面成核 $k_{\mathrm{nuc}}^{\mathrm{mem}}$ 可能比溶液高数倍(未测量) 推断 5. 其他内容(补充内容) 5.1 翻译后修饰与蛋白质量控制 5.1.1 泛素化调控 Qi等(2026)在HD knock-in小鼠模型(Q134)中发现,阻断K6和K9位点特异性泛素化(K>R替换)显著加速mHTT聚集动力学,导致大包涵体形成和专属核定位,运动损伤提前、脑萎缩加速。提示可在模型中添加泛素化修饰项: \[\dfrac{\mathrm{d}[\mathrm{mHTT}_{\mathrm{agg}}]}{\mathrm{d}t} = k_{\mathrm{agg}}^0 \times (1 + \alpha_{\mathrm{K6K9}})\,[\mathrm{mHTT}]\] 其中 $\alpha_{\mathrm{K6K9}} > 0$ 表示K6/K9突变对聚集速率的增强因子。66 公式解释 该式描述泛素化缺失如何加速聚集:$k_{\mathrm{agg}}^0$ 是野生型基础速率,$\alpha_{\mathrm{K6K9}}$ 为K6/K9突变导致的加速倍数。实验显示K>R突变使包涵体更大、更早出现,可将 $\alpha_{\mathrm{K6K9}}$ 设为1.5–3倍估计。 5.1.2 分子伴侣相互作用 系统文献检索确认了5.2.1节提到的数据缺口: 已知 Hsp70 会以 ATP 依赖方式结合并处理 polyQ 相关底物,其结合–释放循环与伴侣蛋白的“分流决策”(折叠或降解)直接相关。 当 Hsp70 完成 ATP 水解并进入 ADP 态时,底物结合会显著更稳定,通常被描述为亲和力提高约10倍(ADP-Hsp70 相对 ATP-Hsp70),因此底物在伴侣体系中的驻留时间会被放大。 现有研究普遍认为 Hsp70 与 polyQ 系统的功能性互作存在长度依赖性,但“长度如何改变 $K_d$ 或 $k_{\mathrm{on}}/k_{\mathrm{off}}$”仍缺少可复用的定量序列。 缺失 polyQ-Hsp70直接结合的 $K_d$ 值在主流文献中未见报道 相关测量仅限于Hsp70与其他底物:多肽结合 $K_d$ 约 50 μM(bacterial DnaK),小分子抑制剂 $K_d$ 约 70 μM 不同长度 polyQ 对 $K_d$、$k_{\mathrm{on}}$ 与 $k_{\mathrm{off}}$ 的影响尚未系统量化,因此目前很难把“伴侣缓冲”写成可跨长度拟合的动力学项。 这证实了蛋白质质量控制(PQC)系统与polyQ相互作用的定量参数是建模的关键缺口。来源, 来源 蛋白质质量控制相互作用缺口:目前公开文献几乎没有报告不同 Q 长度与 Hsp70/Hsp90、泛素连接酶或蛋白酶体亚基之间的结合常数($K_d$)或速率($k_{\mathrm{on}}/k_{\mathrm{off}}$),即便是 HTT 也只停留在定性 Co-IP 规模。为建立多尺度模型,需要系统测定这些参数,才能把“聚集与清除”写成耦合动力学。 已知的伴侣结合实例:现有定量数据表明 Hsc70 以微摩尔亲和力结合同一 Httex1 分子 N17 区域,且结合位点竞争了 Httex1 之间的同型接触并抑制聚集;但尚不清楚 $K_d$ 是否随 Q 长度改变,因而该参数仍需实验补全。39 5.2 小分子与抑制剂 5.2.1 硒纳米颗粒(SeNPs) Torricella等(2025)通过NMR、荧光免疫染色和TEM揭示,SeNPs以纳摩尔亲和力选择性结合到httQ35可延伸末端,亚化学计量地减少纤维形成速率。SeNPs不改变预成核四聚化,而是减少自由可延伸末端池: \[k_{\mathrm{eff}} = k_{+}\,\frac{[E]_{\mathrm{free}}}{[E]_{\mathrm{total}}} = k_{+}\left(1-\frac{[\mathrm{SeNP}]}{K_d+[\mathrm{SeNP}]}\right)\] 其中 $[E]$ 为可延伸末端浓度,$K_d\approx$ nM 量级。为现有统一动力学模型(2.2节)提供了抑制剂调控的量化方案。67 公式解释 SeNPs通过“封端”机制抑制聚集:$k_{+}$ 是自由末端的延伸速率常数($6.4 \times 10^5\ \mathrm{M}^{-1}\mathrm{h}^{-1}$),SeNPs结合后使有效延伸速率 $k_{\mathrm{eff}}$ 按Langmuir吸附式下降。纳摩尔亲和力意味着极低浓度即有效。 5.2.2 姜黄素 Jain等(2025)发现,亚化学计量量的姜黄素影响HttEx1的一级和/或二级成核事件,延长预聚集滞后期。破坏的聚集过程改变了聚集体结构及其细胞代谢特性:当施用于神经元细胞时,“突破”的蛋白聚集体诱导的细胞应激显著低于无抑制剂形成的聚集体。电镜、SAXS和固态NMR分析鉴定了纤维结构变化,探测了fuzzy coat中的flanking结构域和纤维核心,后者变化与polyQ β-hairpin结构的存在或缺失相关。该发现强调小分子抑制剂调控蛋白错误折叠景观的多方面后果,对HD和其他淀粉样疾病治疗策略有潜在意义。68 5.2.3 大分子共溶质 Torricella等(2024)通过NMR监测发现,多糖dextran-20和蛋白lysozyme主要通过改变预成核四聚化平衡影响httQ35聚集动力学,导致“预形成”httQ35四聚体浓度大幅变化。对较短非聚集变体httQ7的类似效应支持该结论。该研究为2.2节四聚体模型提供了共溶质调控的量化参数: \[K_{\mathrm{tetra}}^{\mathrm{eff}} = K_{\mathrm{tetra}}^0 \times f(\text{cosolute}), \quad [T] = K_{\mathrm{tetra}}^{\mathrm{eff}} [M]^4\] 其中 $f(\text{cosolute})$ 为共溶质依赖的修正因子,dextran-20和lysozyme的具体影响可从NMR交叉峰强度拟合获得。69 公式解释 大分子共溶质(如dextran、lysozyme)不直接抑制纤维化,而是改变四聚体平衡常数 $K_{\mathrm{tetra}}$。例如,若 $f < 1$,则四聚体浓度 $[T]$ 降低,延缓成核;若 $f > 1$,则加速。实验显示两者效应类似,可作为调控滞后期的工具。 N17区域因其两亲性螺旋结构可插入脂双层膜,导致线粒体和ER膜富集度与polyQ长度相关,并引起轴突运输缺陷。35 细胞层面的观察 在细胞模型中,Q23 Httex1 往往只呈现较低水平的线粒体与 ER 膜定位,整体更接近弥散分布。 当 polyQ 扩展到 Q73 时,常观察到线粒体与 ER 膜上的富集显著增强,提示定位与聚集倾向可能被同一长度变量共同驱动。 膜富集程度与聚集倾向的正相关关系提示:膜表面可能提供局部二维浓缩与构象引导,从而成为实际成核路径的一部分。 polyQ 扩展还会造成轴突运输缺陷并间接改变线粒体与囊泡运输,从而把“定位–聚集–功能损伤”三者耦合到同一条病理链路中。 定量模型 膜相互作用的详细热力学模型、膜富集度的Q长度依赖公式、以及膜定位与聚集的耦合方程见第3.3节。 轴突运输 扩展polyQ通过多种机制干扰轴突运输: 扩展 polyQ 可能导致马达蛋白(kinesin/dynein)的募集与装配异常,从而降低有效运输复合体的形成概率。 扩展 polyQ 可能引入货物装载与释放的动力学缺陷,使运输过程出现更频繁的停滞或错误卸载。 扩展 polyQ 形成的寡聚体或聚集体还可能造成微管轨道的空间阻塞,进一步降低长距离运输的通量与稳定性。 这些效应导致有效运输速率下降,但具体的Q长度依赖参数尚未测定。 PC12 神经元样毒性曲线:可诱导 PC12 克隆在表达 Q74 Httex1 后,细胞活率于诱导后第 4、6、8 天分别降至约 84%、75%、60%,而 Q23 克隆维持 >95% 生存;终末分化细胞中 Q74 导致 >80% 细胞死亡,广谱 caspase 抑制剂 zVAD-fmk 可部分挽救活率,为累积损伤模型提供现实参数。65 线粒体/氧化应激读数:关于激活 NRF2 或补充 Dopamine-Agapriline 改善 Q111–Q140 神经元凋亡的说法暂无公开实验证据;若在模型中占位,可用 $\mathrm{d}(\mathrm{ROS})/\mathrm{d}t = k_Q - k_{\mathrm{NRF2}}$ 形式暂存,需后续文献验证。 6. 全文总结 本文档系统整理了多聚谷氨酰胺(polyQ)疾病的长度依赖模型,从临床预测到分子机制再到细胞行为,多层次展示了重复长度这一核心参数如何决定疾病表型。以下总结主要结论: 6.1 核心发现 长度阈值效应的普遍性 35-40个谷氨酰胺的临界阈值是polyQ疾病的普适特征。无论是HD、SCA还是其他polyQ疾病,致病性都在此阈值附近显现,且毒性与Q长度呈正相关。这一阈值效应在不同疾病中的具体数值虽有差异(32Q-54Q不等),反映了蛋白背景对polyQ毒性的调制作用。 连续变化与相变跃迁的共存 临床风险随Q长度主要呈连续陡峭的S形曲线,未见绝对突变式跳变;但分子层面的聚集动力学在Q23-26(成核阶数跃迁)和Q36-40(速率陡增)存在类相变特征。两层信息应同时纳入模型:临床层面需用连续函数(Weibull、对数线性),分子层面需考虑阶跃式转变。 多尺度耦合的关键节点 体细胞扩增、HTT1a异常剪接和LLPS是连接分子长度与临床表型的三个关键机制: 体细胞扩增将遗传Q长度转化为组织特异性长度分布($Q_{\mathrm{germline}} \rightarrow Q_{\mathrm{tissue}}$) HTT1a比例随CAG长度线性增加($\mathrm{HTT1a/HTT} \propto \text{CAG}$),提供毒性片段来源 LLPS临界浓度随Q下降($C_{\mathrm{sat}}(Q) \propto Q^{-1}$),促进凝聚体形成 这些机制可以通过耦合方程组纳入定量模型。 6.2 可复用的定量框架 临床预测模型 Weibull分布:$h(Q, t) = (k_{\mathrm{shape}} \cdot k_{\mathrm{scale}} \cdot t^{k_{\mathrm{shape}}-1})$,其中$k_{\mathrm{scale}} = \exp(\beta_0 - \beta_1 Q)$ 对数线性AO模型:$\ln(AO) = \beta_0 - 0.049\,Q_{\mathrm{exp}} + 0.013\,Q_{\mathrm{norm}} + \sigma\,\mathrm{z}(SE)$ 分段模型:$Q\le 44$和$Q\ge 45$采用不同斜率,捕捉阈值附近的陡峭变化 这些公式可直接用于风险预测、个体化随访和遗传咨询。 分子动力学模型 成核-延伸方程组:包含单体浓度、四聚体平衡、纤维延伸速率和自毒性项 相分离热力学:Flory-Huggins相互作用参数$\chi(Q)$和临界体积分数$\phi_c(Q) = N^{-1/2}$的Q依赖性 Lag-time公式:$t_{\mathrm{lag}}(Q) \propto \exp[-\gamma(Q-Q_c)]$,捕捉Q超过阈值时的指数级加速 这些模型为理解polyQ聚集的物理机制提供了定量框架。 6.3 研究缺口与未来方向 定量参数严重缺失 尽管框架已建立,但许多关键参数尚未测定: PQC系统(Hsp70、Hsp90、泛素连接酶、蛋白酶体)与polyQ蛋白的结合常数和速率 膜相互作用的定量热力学参数($K_d$、$\Delta G$、富集系数) LLPS相图的完整测量($\chi(Q)$、$T_c(Q)$、binodal/spinodal曲线) 跨不同Q长度的系统动力学数据 实验到模型的映射挑战 时间尺度跨越:体外lag-time(小时-天)到临床发病年龄(数十年)的映射需引入温度、浓度、细胞环境等校正因子 体内复杂度:细胞拥挤、分子伴侣、膜表面等因素如何定量纳入模型 物种差异:动物模型(酵母、线虫、小鼠)的结果如何外推到人类 跨蛋白统一建模 当前研究主要集中于HTT和少数ataxin蛋白,其他polyQ系统(SBMA、DRPLA、SCA各亚型)的定量数据稀缺。建立跨蛋白的统一模型需要: 标准化实验协议和测量指标 公共数据库和数据共享框架 考虑蛋白背景序列影响的模型扩展 6.4 对建模和实验设计的建议 对建模者的建议 使用分段函数捕捉阈值效应:在Q36-40附近引入平滑过渡函数(如sigmoid或piecewise-linear) 纳入体细胞扩增噪声项:将$Q$视为随年龄和组织演化的随机变量,而非固定参数 耦合聚集与清除动力学:允许出现双稳态和滞后环,反映PQC系统的容量限制 验证跨尺度一致性:确保分子层面参数(如$k_{\mathrm{nuc}}$)与临床层面观测(如AO)的自洽性 对实验设计者的建议 优先测定Q长度依赖的参数:在不同Q长度下系统测量$K_d$、$k_{\mathrm{on/off}}$、$C_{\mathrm{sat}}$等 建立标准化实验条件:统一缓冲液、温度、浓度等参数,便于跨研究比较 同时测量多尺度读数:在同一体系下同时获取分子、细胞和功能层面的数据 关注中等扩增区间:Q23-35的临床意义尚未明确,但可能包含关键的相变信息 通过理论模型与精心实验的结合,我们有望建立更准确的polyQ疾病预测模型,为早期诊断和干预提供定量基础。 7. 参考文献 Chen, S., Ferrone, F. A., & Wetzel, R. (2002). Huntington’s disease age-of-onset linked to polyglutamine aggregation nucleation. Proc Natl Acad Sci USA, 99(Suppl 4), 16483–16488. DOI: 10.1073/pnas.182276099 ↩ ↩2 Thakur, A. K., et al. (2009). Polyglutamine Disruption of the Huntingtin Exon 1 N Terminus Promotes a Conformational Switch to β-Hairpin. PLoS Comput Biol, 5(8), e1000452. DOI: 10.1371/journal.pcbi.0010030 ↩ ↩2 Khare, S. D., Ding, F., & Dokholyan, N. V. (2005). Molecular Origin of Polyglutamine Aggregation in Neurodegenerative Diseases. PLoS Comput Biol, 1(1), e34. DOI: 10.1371/journal.pcbi.0010034 ↩ ↩2 ↩3 Walters, R. H., Jacobson, K. H., Pedersen, J. A., & Murphy, R. M. (2009). Examining Polyglutamine Peptide Length. J Mol Biol, 394(1), 127–135. DOI: 10.1016/j.jmb.2009.09.016 ↩ Swami, M., et al. (2009). Somatic expansion of the Huntington’s disease CAG repeat in the brain is associated with disease progression. Am J Hum Genet. ↩ ↩2 ↩3 Moss, D. J. H., et al. (2023). Somatic instability associates with cortical atrophy in HD. Nat Neurosci. ↩ ↩2 ↩3 FOXP2 polyQ length studies: 多篇文献表明FOXP2的polyQ长度变化影响LLPS和功能,具体文献待查。相关研究 ↩ Lee, J.-M., et al. (2021). Age of onset of Huntington’s disease in carriers of reduced penetrance alleles. J Neurol Neurosurg Psychiatry. ↩ Gruber, D., et al. (2020). Late-onset Huntington’s disease: clinical features and CAG repeat characteristics. Eur J Neurol. ↩ Kar, K., Jayaraman, M., Sahoo, B., Kodali, R., & Wetzel, R. (2011). Critical nucleus size for disease-related polyglutamine aggregation is repeat-length dependent. Nat Struct Mol Biol. DOI: 10.1038/nsmb.1992 ↩ ↩2 Landrum, E., & Wetzel, R. (2014). Biophysical underpinnings of the repeat length dependence of polyglutamine amyloid formation. J Biol Chem. DOI: 10.1074/jbc.C114.552943 ↩ ↩2 ↩3 ↩4 Morley, J. F., et al. (2002). The threshold for polyglutamine-expansion protein aggregation and cellular toxicity is dynamic and influenced by aging in C. elegans. PNAS. ↩ ↩2 Tam, S., et al. (2006). Normal-repeat polyglutamine peptides accelerate aggregation nucleation and cytotoxicity of expanded polyglutamine proteins. PNAS. ↩ Suarez, M., et al. (2022). Concentration-dependent phase transitions of mHttex1 and modulation by profilin. J Mol Biol. ↩ ↩2 Langbehn, D. R., et al. (2004). A new model for prediction of the age of onset and penetrance for Huntington’s disease based on CAG length. Am J Hum Genet. ↩ ↩2 ↩3 ↩4 Lee, J.-M., et al. (2012). Fully dominant modifier model of genetic modifiers in Huntington disease. PLoS Genet. ↩ ↩2 ↩3 Poirier, M. A., et al. (2002). Huntington’s disease age-of-onset linked to polyglutamine aggregation nucleation. Nat Neurosci. ↩ Hoschek, H. A., et al. (2024). HTT1a transcripts accumulate with age and CAG repeat length. Sci Adv. ↩ ↩2 Pu, J., et al. (2024). m6A reader YTHDF reduces mutant HTT exon1 toxicity by suppressing HTT1a. Cell Rep. ↩ Tezenas du Montcel, S., et al. (2014). Modeling age at onset in spinocerebellar ataxias. Brain. ↩ ↩2 Jacobi, H., et al. (2015). Natural history and SARA progression in SCA1/2/3/6. Lancet Neurol. ↩ Jacobi, H., et al. (2016). Determinants of survival in spinocerebellar ataxias. Ann Neurol. ↩ Elert-Dobkowska, E., et al. (2021). Genotype–phenotype correlations in SCA7 families. Orphanet J Rare Dis. ↩ Joncourt, F., et al. (2024). Clinical staging and respiratory decline in SCA7. Front Neurol. ↩ Turon-Viñas, E., et al. (2023). Corneal endothelial cell density loss tracks CAG repeat length in SCA7. Br J Ophthalmol. ↩ Toyoshima, Y., et al. (2023). Clinical spectrum of SCA17 and CAG/CAA repeat size. Mov Disord. ↩ Park, J., et al. (2024). Voxel-based morphometry correlates with TBP repeat size in SCA17. Sci Rep. ↩ Igarashi, S., et al. (1996). Intergenerational instability and phenotypes in DRPLA. Nat Genet. ↩ Akamine, H., et al. (2022). Repeat length correlates with phenotype in DRPLA. Mol Genet Genomic Med. ↩ Abe, Y., et al. (2024). CAG repeat length predicts milestones in DRPLA. Neurology. ↩ Querin, G., et al. (2023). Determinants of disease onset in spinal and bulbar muscular atrophy. J Neurol. ↩ Lee, J.-H., et al. (2015). Clinical features and CAG length correlation in Korean SBMA patients. J Clin Neurol. ↩ Ni, W., et al. (2024). Hormonal and genetic modifiers of muscle strength in SBMA. Neuromuscul Disord. ↩ Elena-Real, C. A., et al. (2023). Structural features of mutant huntingtin correlate with disease severity. Nat Struct Mol Biol. ↩ ↩2 ↩3 Atwal, R. S., et al. (2014). N17 targeting of membranes modulates huntingtin toxicity. Mol Cell. ↩ ↩2 ↩3 ↩4 ↩5 ↩6 ↩7 ↩8 Atwal, R. S., et al. (2014). N17 targeting of membranes modulates huntingtin toxicity. Mol Cell. ↩ ↩2 ↩3 Elena-Real, C. A., et al. (2023). The structure of pathogenic huntingtin exon 1 defines the bases of its aggregation propensity. Nat Struct Mol Biol. DOI: 10.1038/s41594-023-00920-0 ↩ ↩2 Mohanty, P. R., et al. (2025). Transient interdomain interactions control huntingtin self-assembly. Nat Struct Mol Biol. ↩ ↩2 Lakhani, B., et al. (2017). Emerging β-sheet rich conformations in super-compact Huntingtin exon-1 mutant structures. J Am Chem Soc. ↩ ↩2 ↩3 Vitalis, A., Lyle, N., & Pappu, R. V. (2009). Thermodynamics of $\beta$-Sheet Formation in Polyglutamine. Biophys J, 97(1), 303–311. DOI: 10.1016/j.bpj.2009.05.003 ↩ Jakubek, R. S., Workman, R. J., White, S. E., & Asher, S. E. (2019). Polyglutamine Solution-State Structural Propensity Is Repeat-Length-Dependent. J Phys Chem B, 123(19), 4193–4203. DOI: 10.1021/acs.jpcb.9b01433 ↩ Yoo, J.-N., et al. (2025). Concentration-dependent structural transition of huntingtin protein in Huntington’s disease. Biophys Chem. DOI: 10.1016/j.bpc.2025.107473 ↩ ↩2 Bhattacharyya, A. M., Thakur, A. K., & Wetzel, R. (2005). Polyglutamine aggregation nucleation: thermodynamics of a highly unfavorable protein folding reaction. Proc Natl Acad Sci U S A. DOI: 10.1073/pnas.0501651102 ↩ Sugaya, K., & Matsubara, S. (2012). Quantitative connection between polyglutamine aggregation kinetics and neurodegenerative process in patients with Huntington’s disease. Molecular Neurodegeneration, 7, 20. DOI: 10.1186/1750-1326-7-20 ↩ Vieweg, S., et al. (2016). An intein-based strategy for the production of tag-free huntingtin exon 1 proteins enables new insights into the polyglutamine dependence of Httex1 aggregation and fibril formation. J Biol Chem. ↩ ↩2 Sahoo, B., Singer, D., Kodali, R., Züchner, T., & Wetzel, R. (2014). Aggregation Behavior of Chemically Synthesized, Full-Length Huntingtin Exon1. Biochemistry, 53(24), 3897–3907. DOI: 10.1021/bi500300c ↩ ↩2 ↩3 ↩4 Jian, X., et al. (2023). Self-poisoning polymer crystal initiates polyQ aggregation. eLife. ↩ ↩2 Tam, S., et al. (2006). Normal-repeat polyglutamine peptides accelerate aggregation nucleation and cytotoxicity of expanded polyglutamine proteins. PNAS. ↩ Sarkar, S., et al. (2024). Unified kinetic model of huntingtin exon1 aggregation. Adv Sci. ↩ ↩2 ↩3 ↩4 Mishra, R., et al. (2024). A Targetable Self-association Surface of the Huntingtin exon1 Helical Tetramer Required for Assembly of Amyloid Pre-nucleation Oligomers. J Mol Biol. DOI: 10.1016/j.jmb.2024.168607 ↩ Kulshrestha, A., et al. (2025). Multiscale Simulations Elucidate the Mechanism of Polyglutamine Aggregation and the Role of Flanking Domains in Fibril Polymorphism. J Phys Chem B. DOI: 10.1021/acs.jpcb.5c06627 ↩ ↩2 Kulshrestha, A., et al. (2025). Multiscale simulations elucidate the mechanism of polyglutamine aggregation and the role of flanking domains in fibril polymorphism. bioRxiv. DOI: 10.1101/2025.05.19.654960 ↩ ↩2 Dekker, M., et al. (2025). A Coarse-Grained MD Model for Disorder-To-Order Transitions in PolyQ Aggregation. J Chem Theory Comput. DOI: 10.1021/acs.jctc.5c00384 ↩ ↩2 Vaglietti, C., et al. (2023). PolyQ length-based molecular encoding of vocalization frequency in FOXP2. iScience, 26(1), 105720. DOI: 10.1016/j.isci.2022.105720 ↩ ↩2 ↩3 ↩4 ↩5 ↩6 ↩7 ↩8 ↩9 ↩10 ↩11 ↩12 ↩13 ↩14 ↩15 ↩16 Crick, D. C. (2011). Biophysical Underpinnings of the Repeat Length Dependence of Polyglutamine Aggregation. PhD Thesis, University of North Carolina at Chapel Hill.(本地PDF:E:\graduate_study\research\IDP\zgq-length\background\papers\Biophysical Underpinnings of the Repeat Length Dependence of Polyglutamine Aggregation.pdf) ↩ ↩2 ↩3 Peskett, T. R., et al. (2018). A liquid-to-solid phase transition of huntingtin exon1. Mol Cell. ↩ ↩2 ↩3 ↩4 ↩5 ↩6 Starikov, E. B., Lehrach, H., & Wanker, E. E. (1999). Folding of Oligoglutamines: A Theoretical Approach Based Upon Thermodynamics and Molecular Mechanics. J Biomol Struct Dyn, 17(3), 409–427. DOI: 10.1080/07391102.1999.10508374(本地PDF:E:\graduate_study\research\IDP\zgq-length\background\papers\Folding of Oligoglutamines A Theoretical Approach Based Upon Thermodynamics and Molecular Mechanics.pdf)。 ↩ ↩2 ↩3 ↩4 ↩5 ↩6 ↩7 ↩8 polyQ蛋白LLPS定量模型与长度依赖性的研究综述。本地PDF:E:\graduate_study\research\IDP\zgq-length\background\polyQ蛋白LLPS定量模型与长度依赖性的研究综述.pdf ↩ ↩2 ↩3 ↩4 ↩5 Posey, A. E., et al. (2018). Profilin reduces aggregation and phase separation of huntingtin N-terminal fragments. J Biol Chem, 293(15), 5553–5565. DOI: 10.1074/jbc.RA117.001408(本地PDF基于polyQ LLPS综述,完整信息见PDF文献列表) ↩ Liu, Y. H., et al. (2025). PolyQ Expansion Controls Biomolecular Condensation and Aggregation of the N-Terminal Fragments of Ataxin-2. Int J Mol Sci, 26(23), 11538. DOI: 10.3390/ijms262311538 ↩ ↩2 Wijegunawardana, D., et al. (2024). Ataxin-2 polyglutamine expansions aberrantly sequester TDP-43 ribonucleoprotein condensates disrupting mRNA transport and local translation in neurons. Dev Cell. DOI: 10.1016/j.devcel.2024.09.023 ↩ Baidya, L., et al. (2025). Intrinsic stiffness and θ-solvent regime in intrinsically disordered proteins: Implications for liquid-liquid phase separation. PNAS Nexus, 4(2), pgaf039. DOI: 10.1093/pnasnexus/pgaf039 ↩ Baidya, L., et al. (2025). Intrinsic stiffness and θ-solvent regime in intrinsically disordered proteins: Implications for liquid-liquid phase separation. PNAS Nexus. DOI: 10.1093/pnasnexus/pgaf039 ↩ ↩2 Prasad, A., et al. (2025). Rad23B delays ataxin-3 liquid-to-solid phase transition through heterotypic buffering. J Mol Biol. DOI: 10.1016/j.jmb.2025.169351 ↩ Apostol, B. L., et al. (2003). Inducible PC12 cell model of Huntington’s disease shows toxicity and decreased histone acetylation. Hum Mol Genet. ↩ ↩2 ↩3 ↩4 ↩5 ↩6 Qi, P., et al. (2026). Prevention of ubiquitination at K6 and K9 in mutant huntingtin exacerbates disease pathology in a knock-in mouse model. Proc Natl Acad Sci U S A. DOI: 10.1073/pnas.2527258122 ↩ Torricella, F., et al. (2025). Kinetic Mechanism of Substoichiometric Inhibition of Huntingtin Exon-1 Protein Aggregation by Selenium Nanoparticles. Small Sci. DOI: 10.1002/smsc.202500345 ↩ Jain, G., et al. (2025). Inhibitor-based modulation of huntingtin aggregation mechanisms mitigates fibril-induced cellular stress. Nat Commun. DOI: 10.1038/s41467-025-58691-9 ↩ Torricella, F., et al. (2024). Effects of Macromolecular Cosolutes on the Kinetics of Huntingtin Aggregation Monitored by NMR Spectroscopy. J Phys Chem Lett. DOI: 10.1021/acs.jpclett.4c01410 ↩
Field Knowledge
· 2026-01-13
2025年诺贝尔生理学或医学奖:坂口、布伦科和拉姆斯德尔揭示外周免疫耐受之谜
【详细解读】2025年诺贝尔生理学或医学奖:坂口、布伦科和拉姆斯德尔如何揭示外周免疫耐受之谜 引言:免疫系统的悖论性使命 2025年诺贝尔生理学或医学奖授予了三位杰出的科学家:坂口志文(Shimon Sakaguchi)、玛丽·E·布伦科(Mary E. Brunkow)和弗雷德·拉姆斯德尔(Fred Ramsdell),以表彰他们在阐明外周免疫耐受(peripheral immune tolerance)细胞和分子基础方面的融合性发现。这一根本性过程能够主动阻止免疫系统攻击机体自身组织。他们的工作解决了一个长期存在的悖论:一个装备精良、旨在摧毁入侵者的免疫系统,是如何被约束而不引发自我毁灭的。 免疫耐受是免疫系统对自身抗原(即构成机体自身的分子)不产生应答的状态,是维持机体健康的核心原则。这一原则通过两个主要分支来实现。第一个分支是中枢耐受(central tolerance),它在初级淋巴器官——T细胞在胸腺,B细胞在骨髓——中发挥作用,通过克隆删除(clonal deletion)等机制,清除大部分具有强烈自身反应性的淋巴细胞,构成了抵御自身免疫的第一道防线。 然而,中枢耐受并非完美无缺。尽管其筛选过程极为严苛,但仍有相当数量的对自身抗原具有低亲和力反应性的T细胞能够逃逸筛选,进入外周循环系统。这些”漏网之鱼”构成了潜在的威胁,可能在特定条件下被激活,引发自身免疫病。这就引出了对第二个关键检查点的需求,即外周耐受。外周耐受在身体的外周组织和次级淋巴器官中运作,负责控制这些逃逸的自身反应性细胞。其机制包括功能失活(即无能,anergy)、抗原忽视(ignorance)或主动抑制(active suppression)。正是这”主动抑制”的分支,成为了本次诺贝尔奖获得者们研究的核心,他们的工作从根本上定义了这一机制。 本报告将追溯一段从备受争议的历史概念演变为现代免疫学核心支柱的科学历程。报告将详细阐述坂口志文如何鉴定出免疫系统的细胞”卫士”——调节性T细胞(Regulatory T cells, Tregs);布伦科和拉姆斯德尔如何发现其遗传”主开关”——转录因子。最后,本报告将深入探讨这一新认知所开启的深刻且具有双重性的治疗前景,它为自身免疫病和癌症等多种重大疾病带来了革命性的新疗法。 第一章:浴火重生——从”抑制学”的废墟之上 早期假说与”抑制学”的兴衰 在20世纪70年代,免疫学家理查德·格尔雄(Richard Gershon)和近藤(Kondo)首次提出了”抑制性T细胞”(suppressor T cells)的概念,他们假设T细胞不仅能增强免疫应答,还能通过下调某些生物学功能来削弱免疫应答。这一想法在逻辑上极具说服力,因为它暗示了免疫系统必须拥有内在的”刹车”机制,以防止过度反应和自身损伤。 在接下来的十年里,”抑制学”(suppressorology)迅速成为免疫学研究的主流领域。科学家们提出了复杂的细胞相互作用模型,试图解释这些抑制性细胞如何调节免疫反应。然而,到了20世纪80年代中期,这一曾经辉煌的领域却戏剧性地崩塌了。其衰落背后有几个关键原因。首先是”I-J悖论”。当时,研究人员认为一个名为I-J的分子是抑制性T细胞的关键标志物,但分子生物学研究却无法在主要组织相容性复合体(MHC)基因区域内定位到编码I-J的基因,这使得该领域的分子基础受到严重质疑。其次,研究人员始终未能找到稳定且特异的细胞表面标志物来分离和鉴定这些所谓的抑制性细胞。由于无法获得纯化的细胞群体,实验结果往往难以重复,细胞的谱系和功能也充满了不确定性。 这些根本性的问题导致了整个领域的信誉危机。”抑制性T细胞”这一术语几乎从主流科学文献中消失,相关的研究论文难以发表,研究经费也日益枯竭。整个领域被蒙上了一层”污点”,科学界对免疫抑制的概念普遍持怀疑态度。 为突破奠定舞台 理解这段历史对于认识坂口志文工作的开创性至关重要。他所要重新探索的,是一个已经被科学界抛弃和否定的概念。这意味着他必须提供远超常规标准的、无可辩驳的证据,才能克服当时普遍存在的怀疑主义。 第一波抑制性T细胞研究的失败,根源在于缺乏一个可靠的分子”抓手”来识别和纯化目标细胞。分子层面的模糊性(如I-J基因问题)和细胞层面的异质性(无法分离纯净的细胞群体)共同导致了数据的矛盾和领域信誉的崩塌。而坂口志文发起的第二波研究之所以能够成功,恰恰是因为他找到了这样一个决定性的标志物:CD25。他1995年的实验之所以具有里程碑意义,不仅在于观察到了自身免疫的表型,更在于他能将这一表型与一个可通过分子特异性识别的细胞亚群——CD4+CD25+ T细胞——精确地联系起来。 这一标志物的发现,使得对这类细胞进行可重复的分离、鉴定和功能测试成为可能,从而提供了该领域十余年来一直缺失的坚实证据。这段历史深刻地揭示了细胞生物学的一个核心原则:如果没有可靠的方法来识别和分离执行功能的细胞实体,那么其功能就无法被明确界定。坂口志文对标志物的发现,不仅仅是一项观察,更是解开整个领域困局、纠正”抑制学”核心方法论缺陷的关键钥匙。 第二章:坂口志文的突破——鉴定细胞卫士 奠基性兴趣与研究方向 早在医学生时代,坂口志文就对自身免疫的悖论深感兴趣。在当时众多免疫学理论中,他尤其被一个最不受欢迎的假说所吸引:即机体内始终存在自身反应性淋巴细胞,但它们的活性通常受到主动抑制。他坚信,如果能明确定义这些未知的抑制性T细胞,就能揭示自身免疫病的普遍原理。 1995年的关键实验(《免疫学杂志》) 坂口志文的里程碑式研究发表于1995年的《免疫学杂志》(Journal of Immunology),其实验设计严谨而优雅,为调节性T细胞的存在提供了第一个决定性证据。 实验设计
Field Knowledge
· 2025-10-20
解锁T细胞邻近控制:LAG-3与TCR的空间接近如何实现精准自身免疫治疗
【Cell】 解锁免疫“刹车”新用法:强制T细胞“邻近”可精准治疗自身免疫病 一、 本文基本信息 摘要 在自身免疫性疾病中,针对致病性T细胞的治疗一直充满挑战。淋巴细胞活化基因-3(LAG-3)是一种主要在活化T细胞上特异性表达的抑制性检查点受体,已知其可与主要组织相容性复合物II类分子(MHC-II)结合。然而,本研究表明,仅仅与MHC-II相互作用不足以实现LAG-3的最佳功能。相反,由同源肽-MHC-II分子介导的LAG-3与T细胞受体(TCR)的空间邻近性,而非与CD4共受体的邻近性,才是介导CD4+ T细胞抑制的关键。从机制上看,LAG-3通过其胞内FSAL基序与TCR信号组分CD3ε形成凝聚体,从而破坏CD3ε与淋巴细胞特异性蛋白激酶(Lck)的结合。为了利用LAG-3与TCR的邻近性并最大化其依赖的T细胞抑制作用,我们开发了一种Fc功能减弱的LAG-3/TCR抑制性双特异性抗体,以绕过对同源肽-MHC-II的需求。这种方法能够强效抑制CD4$^{+}$和CD8$^{+}$ T细胞,并有效缓解小鼠自身免疫模型的症状。我们的发现揭示了一种复杂且有条件的检查点调控机制,并强调了靶向LAG-3/TCR顺式邻近性对于治疗那些缺乏有效且耐受性良好的免疫疗法的T细胞驱动的自身免疫性疾病具有重要意义。 原文引用信息 Du, J., Chen, H., You, J., Hu, W., Liu, J., Lu, Q.,… & Wang, J. (2025). Proximity between LAG-3 and the T cell receptor guides suppression of T cell activation and autoimmunity. Cell,88,-18. https://doi.org/10.1016/j.cell.2025.06.004 二、 研究背景与科学问题 2.1 免疫系统的“双刃剑”:自身免疫病的困境 免疫系统是人体的精密防御体系,其核心职责是区分“自我”与“非我”,精准清除外来病原体,同时保护自身组织。然而,当这套精密的识别系统出现故障,免疫细胞便会错误地攻击自身健康的组织和器官,导致自身免疫性疾病的发生。这类疾病种类繁多,包括1型糖尿病、类风湿性关节炎、多发性硬化症等,影响着全球数以亿计的人口。在许多自身免疫病中,T淋巴细胞(T细胞)的过度活化被认为是驱动疾病进展的核心元凶。 近年来,自身免疫病的发病率在全球范围内,尤其是在工业化国家,呈现出快速上升的趋势。一个广受关注的理论是“卫生假说”(Hygiene Hypothesis)或其修正版“老朋友假说”(Old Friends Hypothesis)。该假说认为,现代社会过于洁净的环境、抗生素的广泛使用以及生活方式的改变,减少了人们在生命早期接触微生物的机会。这种暴露的缺乏可能导致免疫系统未能得到充分的“训练”和“教育”,使其调控机制发育不全,从而更容易对自身抗原产生过度反应,增加了患过敏性和自身免疫性疾病的风险。这一宏观背景凸显了深入理解免疫调控机制、开发新型精准免疫疗法的迫切性。 2.2 免疫检查点LAG-3:一个熟悉又陌生的“刹车”分子 为了防止免疫反应过度而伤及自身,免疫系统进化出了一系列“刹车”机制,即免疫检查点。这些分子通常是表达在免疫细胞表面的抑制性受体,当它们被激活时,能够抑制免疫细胞的功能,维持免疫稳态。其中,淋巴细胞活化基因-3(LAG-3,又称CD223)是一个关键的抑制性检查点,主要在活化的T细胞上高表达。 LAG-3在结构上是CD4分子的同源物,属于I型跨膜蛋白,其胞外区包含四个免疫球蛋白(Ig)样结构域(D1-D4)。它最主要的配体是主要组织相容性复合物II类分子(MHC-II),同时也能够与纤维蛋白原样蛋白1(FGL1)等其他分子结合。与更为人熟知的PD-1和CTLA-4检查点类似,LAG-3在癌症和慢性感染中标志着T细胞的“耗竭”状态。研究表明,LAG-3与PD-1常常在耗竭的T细胞上共表达,并通过互补的、非冗余的通路协同抑制T细胞功能。因此,靶向LAG-3已成为继PD-1/CTLA-4之后肿瘤免疫治疗的又一重要方向。 2.3 关键科学问题与核心创新点 尽管LAG-3作为免疫“刹车”的重要性已广为人知,但一个核心的谜题长期悬而未决:LAG-3究竟是如何通过与配体结合,将抑制信号传递到细胞内部,从而精确地关闭T细胞的活化引擎的? 其具体的分子作用机制一直不甚明了,被形容为一个“充满谜题的分子”。 本项发表于《细胞》杂志的研究,正是为了解开这个谜题。研究团队通过一系列精巧的实验设计,取得了多项突破性进展,为我们描绘了一幅全新的LAG-3功能图景。其核心创新点可概括为: 揭示了全新的“空间邻近依赖”抑制模型:颠覆了传统的“配体结合即激活”的认知,证明LAG-3的抑制功能需要其与T细胞受体(TCR)在空间上足够靠近才能实现。 阐明了LAG-3在分子水平“釜底抽薪”的物理机制:发现LAG-3的胞内结构域可以直接干扰TCR早期信号通路中关键蛋白复合物的形成。 基于新机制,开创性地设计了“激动型”双特异性抗体BiTS:将基础科学发现迅速转化为一种全新的治疗策略,通过人为强制LAG-3与TCR的邻近来增强其抑制功能。 在多种自身免疫病动物模型中验证了BiTS的强大治疗潜力:证明了这一新策略在1型糖尿病、自身免疫性肝炎和多发性硬化症模型中的有效性,展示了其广阔的临床应用前景。 三、 核心研究内容深度解析 3.1 发现一:LAG-3的抑制功能依赖于其与TCR的“空间邻近” 传统观点认为,抑制性受体的功能主要由其与配体的结合直接触发。然而,本研究的第一个重大发现是,对于LAG-3而言,情况并非如此简单。 研究团队首先构建了一个“简化系统”,即利用不表达天然MHC-II的人工抗原提呈细胞(aAPC)。他们发现,当这些aAPC只表达能够结合LAG-3但不能识别TCR的“非同源”MHC-II分子时,即使LAG-3与MHC-II成功结合,也无法有效抑制T细胞的活化。相反,只有当aAPC表达能够同时被TCR和LAG-3识别的“同源”肽-MHC-II复合物时,LAG-3才表现出强大的抑制功能。 这一现象引出了一个大胆的假设:MHC-II分子可能不仅仅是LAG-3的配体,更像一个“桥梁”,其关键作用是将在细胞膜上原本分离的LAG-3和TCR拉到一起。 只有当LAG-3与TCR在空间上足够“邻近”时,它的抑制“刹车”才能被踩下。 图 1. LAG-3 在同源 pMHC II 存在时介导 CD4+ T 细胞抑制 A. B 细胞与 CD4+ T 细胞界面的配体-受体互作示意图 B. 经 mock 或小鼠 LAG-3 转导的 DO11.10 T 细胞与负载 OVA323-339 同源肽的 LK35.2 B 细胞共培养后的 IL-2 产量(n=3) C. 表达同源 pMHC II 的 aAPC 与 LAG-3+ T 细胞界面互作示意图 D-E. 小鼠 LAG-3 CAR Jurkat NF-κB-GFP 报告 T 细胞与 mock 或 I-Aᵈ OVA323-339 aAPC 共培养的代表性图像及 GFP 信号定量 F-G. 经 mock 或小鼠 LAG-3 转导的 DO11.10 T 细胞被同源 I-Aᵈ OVA323-339 aAPC 激活后的 IL-2 产量及 LAG-3 抑制率(n=3) H-M. 人 LAG-3 CAR 与 mock 或 HLA-DR1 HA306-318 aAPC 共培养的代表性图像及 GFP 信号定量(n=3) N-O. 经 mock 或人 LAG-3 转导、携带 HA1.7 TCR 的 Jurkat NFAT-GFP 报告 T 细胞被 HLA-DR1 HA306-318 aAPC 激活后的 GFP 阳性细胞百分比及 LAG-3 抑制率(n=4) 为了直接验证这一“空间邻近假说”,研究团队设计了一个极为精妙的实验。他们利用一种化学诱导的异源二聚化系统(rapalog系统),在该系统中,可以通过加入一种小分子药物(rapalog)来强制细胞膜上的两种不同蛋白物理性地靠近。他们将TCR的激活分子(抗CD3抗体)与FRB蛋白融合,将非同源MHC-II分子与FKBP12蛋白融合。在没有rapalog时,TCR和LAG-3(通过非同源MHC-II结合)是分离的;而加入rapalog后,FRB和FKBP12结合,从而强制性地将TCR和LAG-3拉到了一起。 结果惊人地清晰:在加入rapalog后,即便没有同源MHC-II的参与,T细胞的活化也受到了显著抑制。这一实验无可辩驳地证明了,LAG-3的抑制功能并不绝对依赖于特定的配体信号,而是由其与TCR的物理邻近性所“授权”的。这是一个全新的、以空间构象为核心的免疫调控模型。 图 2. TCR 空间邻近性对 MHC II/LAG-3 介导的 CD4+ T 细胞抑制至关重要 A. 表达膜锚定抗小鼠 CD3(anti-CD3MT)和非同源 pMHC II(pMHC class IINC)的 aAPC 设计示意图,用于解耦 MHC II 类分子与 CD4+ T 细胞 TCR 的相互作用 B. 流式细胞术检测三种 aAPC 克隆的 anti-CD3MT 表达水平(相对于亲本细胞的 MFI) C-D. 小鼠 LAG-3+ 3A9 T 细胞被上述 aAPC 克隆激活后的 IL-2 产量,有无非同源 I-Aᵇ OVA323-339 的瞬时过表达 E-H. 涉及 anti-CD3MT aAPC 的 IL-2 产量及抑制率分析,包括有无非同源 pMHC II、mock 与 LAG-3 转导的 3A9 T 细胞对比、同型对照与抗 LAG-3(M8)处理对比 I-K. rapalog 诱导异源二聚体实验的示意图、构建设计及流式细胞术检测标记表达(相对于亲本细胞的 MFI) L-N. mock 或小鼠 LAG-3 转导的 3A9 或 DO11.10 T 细胞被表达 rapalog 诱导异源二聚体的 aAPC 激活后的 IL-2 抑制率 3.2 发现二:深入分子内部——LAG-3如何通过破坏“信号凝聚体”来“釜底抽薪” 既然空间邻近是关键,那么下一个问题是:当LAG-3被拉到TCR旁边后,细胞内部究竟发生了什么?为了回答这个问题,研究团队将目光投向了细胞内信号转导的最前沿——生物分子凝聚体(Biomolecular Condensates)。 LAG-3的胞内域(ICD)可以直接与TCR复合物的关键组分CD3ε相互作用。 当LAG-3的ICD被强制带到TCR附近时,它会竞争性地结合CD3e,从而破坏CD3ε与激酶Lck的结合,导致这个至关重要的“信号凝聚体”无法有效形成或维持稳定。 Lck是启动TCR下游信号瀑布的“第一把钥匙”,Lck与CD3ε的分离,相当于直接切断了信号的源头,导致T细胞活化被有效终止。 为了进一步精确定位LAG-3 ICD中的功能区域,研究团队进行了一系列突变分析。结果显示,LAG-3 ICD中两个保守的基序——FSAL基序和EP基序——对于破坏CD3ε/Lck凝聚体以及发挥抑制功能至关重要,而另一个已知的KIEELE基序在此过程中的作用相对较弱。这与以往研究中关于这些基序功能重要性的争论提供了新的见解。 下表总结了关键的突变分析结果: 突变体 突变对象 关键区域/基序 对CD3ε/Lck凝聚体的影响 T细胞抑制功能影响 推断作用 F483A/L486A LAG-3 FSAL Motif 破坏作用丧失 显著减弱 直接参与或稳定与CD3ε的相互作用 ΔEP LAG-3 EP Motif 破坏作用部分减弱 部分减弱 可能通过静电作用协同FSAL基序发挥功能 ΔKIEELE LAG-3 KIEELE Motif 影响较弱 影响较弱 在此邻近模型中的作用非主导 CD3ε BRS基序突变 CD3ε Basic Rich Sequence (BRS) 与LAG-3凝聚能力丧失 - CD3ε上与LAG-3和Lck结合的关键位点 CD3ε mPRS CD3ε Proline Rich Sequence (PRS) 未提及对凝聚体影响(实验显示不影响与LAG-3的凝聚体形成 ) - 非LAG-3与CD3ε结合及破坏CD3ε/Lck凝聚体的关键区域 CD3ε mRK CD3ε RK Motif 未提及对凝聚体影响(实验显示不影响与LAG-3的凝聚体形成 ) - 非LAG-3与CD3ε结合及破坏CD3ε/Lck凝聚体的关键区域 LAG-3 F475A/L478A LAG-3 类似FSAL基序区域(与F483A/L486A类似位置突变) 未明确提及(推测类似F483A/L486A突变体影响 ) 未明确提及(推测抑制功能减弱 ) 可能与FSAL基序功能相似,参与与CD3ε的相互作用 这些发现共同揭示了一个精巧的分子机制:LAG-3通过物理邻近,直接干预TCR信号凝聚体的形成,从源头上关闭了T细胞的活化开关。 图 3. LAG-3 胞内域破坏 CD3ε/Lck 凝聚体 A-F. 溶液中 100 μM 人 CD3ε 与 LAG-3 ICD 的凝聚体形成实验,展示野生型或磷酸化状态的 CD3ε 及野生型或 F483A/L486A 突变的 LAG-3(n=5) G-O. 支持脂质双层(SLB)上 5 μM 磷酸化 CD3ε 与 Lck(开放形式 K273R/Y505F 突变或截短型 LckUD-SH3-SH2)的凝聚体形成实验,有无 5 μM 人或小鼠 LAG-3 ICD 及其他抑制性受体 ICD(n=5/6) P-S. 流式细胞术验证 mock 或人 / 小鼠 LAG-3(野生型或突变体)转导的 Jurkat NFAT-GFP 报告 T 细胞或 3A9 T 细胞中 LAG-3 的表达及抑制功能 3.3 发现三:从机制到疗法——“强制邻近”催生新型激动型双特异性抗体BiTS 这项研究最激动人心的部分在于,它没有停留在基础机制的发现,而是迅速将这一新知识转化为了创新的治疗策略。既然“强制邻近”是激活LAG-3抑制功能的关键,那么是否可以设计一种药物来主动实现这一点呢? 答案是肯定的。研究团队为此设计了一种双特异性抗体(Bispecific Antibody, BsAb)。 一个“手臂”:靶向并结合T细胞表面的LAG-3分子。 另一个“手臂”:靶向并结合同一T细胞表面的TCR复合物。 通过这种设计,BiTS分子就像一根强有力的“分子绳索”,将LAG-3和TCR强行“捆绑”在一起,人为地创造并稳定了发挥抑制功能所必需的“空间邻近”状态。为了避免在体内引发不必要的免疫反应(如抗体依赖的细胞毒性作用,ADCC),研究人员还对其Fc片段进行了N297G突变改造,使其成为一种“功能沉默”的抗体。 体外实验结果证实了这一设计的有效性:BiTS能够以剂量依赖的方式,强效抑制CD4$^{+}$和CD8$^{+}$ T细胞的活化。重要的是,这种抑制作用是严格依赖于LAG-3的——当作用于不表达LAG-3的T细胞时,BiTS的抑制效果大幅减弱,证明其功能确实是通过激活LAG-3通路实现的。 图 4. LAG-3/TCR 双特异性抗体(BiTS)介导 LAG-3 依赖性 T 细胞抑制 A. LAG-3/TCR BiTS 设计及 LAG-3 结合破坏的 BiTS 突变体(BiTSmut)示意图 B. 流式细胞术检测 100 nM BiTS 和 BiTSmut 与 TCR+LAG-3+ 细胞、TCR+LAG-3- 或 TCR-LAG-3+ B3Z 细胞的结合 C. mock 或 LAG-3 转导的 CD4+ TCR+(3A9)或 CD8+ TCR+(B3Z)T 细胞被 anti-CD3MT aAPC 激活后,在有无 10 nM BiTS、BiTSmut 或同型对照存在下的 IL-2 产量(n=3) D-E. 不同浓度 BiTS 或 BiTSmut 对 IL-2 产量的抑制作用及半最大抑制浓度(IC50)(n=3) 图 5. LAG-3/TCR BiTS 强效抑制抗原特异性 CD4+ 和 CD8+ T 细胞反应 A. CD4+ TCR+(3A9)或 CD8+ TCR+(B3Z)T 细胞分别被 I-Aᵏ HEL50-62 CD86+ aAPC 或负载 OVA257-264 的 mutuDC 激活后,在 10 nM BiTS、BiTSmut 或同型对照存在下的 IL-2 产量(n=3) B-D. BiTS、BiTSmut 或 Ly-TCR(无 LAG-3 臂)对 mock 或小鼠 LAG-3 转导的 CD4+ TCR+(3A9)或 CD8+ TCR+(B3Z)T 细胞 IL-2 产量的抑制作用及 IC50(n=3) 3.4 深度对比:BiTS vs. Relatlimab——同一靶点,两种截然相反的疗法 BiTS的出现,使得LAG-3这个靶点呈现出一种迷人的“双重人格”。为了更好地理解BiTS的创新性,有必要将其与已上市的LAG-3靶向药物——Relatlimab——进行对比。 Relatlimab是百时美施贵宝(BMS)开发的全球首个获批的LAG-3抑制剂,它与PD-1抑制剂Nivolumab组成的复方制剂(商品名Opdualag)已被批准用于治疗黑色素瘤。Relatlimab是一种拮抗型(Antagonist)单克隆抗体,其作用机制是阻断LAG-3与其配体(如MHC-II)的结合。在癌症治疗中,肿瘤细胞利用LAG-3通路来抑制T细胞的抗肿瘤活性;而Relatlimab通过切断这条通路,相当于“松开刹车”,重新释放T细胞的杀伤力,从而达到治疗癌症的目的。 而本研究中的BiTS则恰恰相反,它是一种激动型(Agonist)双特异性抗体。它的目标不是“松开刹车”,而是“踩死刹车”。在自身免疫病中,T细胞的过度活化是致病根源,因此需要抑制而非激活它们。BiTS通过强制LAG-3与TCR邻近,主动激活LAG-3的抑制信号,从而精准地“沉默”致病T细胞。 下表清晰地对比了这两种基于同一靶点的、截然不同的治疗策略: 特性 BiTS(本文研究) Relatlimab(已上市抗癌药) 抗体类型 双特异性激动型抗体 (Bispecific Agonist) 单克隆拮抗型抗体 (Monoclonal Antagonist) 作用机制 强制LAG-3与TCR邻近,激活抑制信号 阻断LAG-3与配体结合,解除抑制信号 对T细胞功能的影响 抑制 (Suppression) 激活 (Activation) 治疗领域 自身免疫病 (Autoimmune Diseases) 癌症 (Cancer) 核心科学原理 空间邻近诱导信号调控 配体-受体阻断 这种“一靶两用”的现象,深刻揭示了对靶点生物学机制的深入理解,是如何催生出功能完全相反但同样具有巨大潜力的创新疗法的。 3.5 发现四:BiTS在多种自身免疫病模型中展现强大治疗潜力 基础机制的阐明和创新分子的设计最终都需要在活体模型中得到验证。研究团队在多种经典的自身免疫病小鼠模型中评估了BiTS的治疗效果,结果令人振奋: 1型糖尿病模型(RIP-OVA模型):这是一种由CD8$^{+}$ T细胞攻击胰岛β细胞导致的疾病模型。结果显示,与对照组相比,BiTS治疗显著延缓了糖尿病的发生,并有效保护了小鼠免于发病。组织学分析也证实,BiTS治疗组小鼠胰岛的T细胞浸润(胰岛炎)程度显著降低。 自身免疫性肝炎模型:这是一种由记忆性CD8$^{+}$ T细胞驱动的肝脏炎症模型。BiTS治疗极大地减轻了肝脏的炎症损伤和CD8$^{+}$ T细胞浸润。单细胞测序分析进一步揭示,BiTS不仅减少了致病性T细胞的数量,还深刻改变了它们的活化状态和功能表型。 多发性硬化症模型(EAE模型):这是一种主要由CD4$^{+}$ T细胞介导的中枢神经系统自身免疫病模型。无论是在疾病发生前的预防性治疗,还是在疾病高峰期的治疗性给药,BiTS均能有效降低疾病的临床评分,改善小鼠的神经功能症状。 这些在不同疾病背景、由不同T细胞亚群(CD4$^{+}$或CD8$^{+}$)驱动的模型中取得的一致性阳性结果,强有力地证明了BiTS作为一种平台型疗法的巨大潜力,有望用于治疗多种T细胞介导的自身免疫性疾病。 图 6. LAG-3/TCR BiTS 抑制 CD8+ T 细胞活化并预防糖尿病 A-C. WT 或 Lag3-/- OT-I T 细胞被负载 OVA257-264 肽的 mutuDC 激活后,在 10 nM 同型对照、BiTSmut 或 BiTS 存在下的流式细胞术点图及 IL-2/IFN-γ 产量(n=3) D-E. RIP-OVA 糖尿病模型示意图及糖尿病发病率(每组 n=8,采用 Kaplan-Meier 检验) F-G. 糖尿病小鼠胰岛的代表性组织学图像(无胰岛炎、peri - 胰岛炎、胰岛炎)及组织学评分(每组分析 >10 个胰岛) 图 7. LAG-3/TCR BiTS 改善 CD8+ T 细胞介导的肝炎 A-B. IND-PALF 患者与 Dx-PALF 患者肝组织中免疫相关基因的 RNA-seq 分析及 GO 富集分析 C-D. 抗 - 4-1BB 处理后肝炎小鼠的肝损伤指标(ALT)及肝内 T 细胞亚群分析(n=4-9/5) E-I. 肝炎小鼠肝内免疫细胞的单细胞 RNA-seq 分析,包括 UMAP 降维、细胞类型鉴定及 CD8+ T 细胞亚群的效应 / 记忆基因特征评分 J-O. 肝内 CD8+ T 细胞与髓系细胞的配体 - 受体互作分析及炎症因子水平检测(n=5) 五、 关键结论与批判性总结 5.1 本文关键结论 本研究系统性地揭示了免疫检查点LAG-3的一个全新且精密的调控机制,并基于此开发了一种极具潜力的自身免疫病新疗法。其核心结论如下: 机制重定义:LAG-3的抑制功能并非简单地由配体结合触发,而是“有条件的”,其高效发挥作用的前提是与TCR在细胞膜上形成空间邻近。 分子通路阐明:在空间邻近的条件下,LAG-3的胞内结构域(ICD)通过与CD3ε直接相互作用,破坏了TCR活化所必需的CD3ε/Lck信号凝聚体,从而“釜底抽薪”式地终止了T细胞活化。 创新疗法开发:基于“强制邻近”原理,研究团队设计了激动型双特异性抗体BiTS,它能将LAG-3和TCR强行拉近,从而高效、精准地抑制T细胞功能。 广谱潜力验证:BiTS在1型糖尿病、自身免疫性肝炎和多发性硬化症等多种动物模型中均展现出强大的治疗效果,证明了其作为一种平台型疗法治疗T细胞介导的自身免疫病的广阔前景。 5.2 专家评述与展望 这项工作无疑是免疫学和药物研发领域的一项里程碑式研究,其意义深远,不仅重塑了我们对一个重要免疫检查点的认知,更开辟了一条全新的治疗途径。 意义与优势 LAG-3生物学的范式转移:该研究将我们对LAG-3的理解从一个简单的“配体-受体”模型,提升到了一个复杂的、依赖于“空间构象和信号组织”的全新层面。这种“条件性刹车”机制或许可以解释一个长期存在的现象:为何相比于PD-1,LAG-3在生理条件下的抑制功能显得更为温和。这可能是一种进化上的安全设计,通过增加激活门槛,避免在不适当的情况下过度抑制T细胞,从而为免疫系统提供了更精细的调控层次。 开创“邻近诱导信号调控”新药模式:BiTS的设计理念极具开创性。在药物研发领域,利用小分子诱导蛋白邻近以降解靶蛋白的PROTACs和分子胶技术已是热点。而BiTS则将这一“邻近诱导”概念从“降解”扩展到了“信号调控”领域,代表了一类全新的治疗模式—— 邻近诱导型生物大分子(Proximity-Inducing Biologics)。这为未来针对其他受体对的信号调控药物开发提供了宝贵的范例和理论基础。 严谨优雅的实验设计:从构建简化的人工细胞系统,到利用化学诱导工具精准验证核心假说,再到体外凝聚体实验和多种动物模型的验证,整个研究逻辑链条清晰、证据确凿,体现了极高的科学严谨性和创新性。 局限性与未来方向 尽管这项研究取得了重大突破,但从实验室走向临床,仍有诸多挑战需要克服,同时也为未来的研究指明了方向。 临床转化的挑战: 双特异性抗体的开发壁垒:作为一种复杂的生物大分子,双特异性抗体的生产工艺(CMC)复杂、成本高昂。此外,其独特的结构可能带来不可预测的免疫原性(即诱发抗药抗体),以及独特的毒性谱,如可能引发细胞因子释放综合征(CRS)或神经毒性,这些都是临床开发中必须密切关注和解决的问题。 患者选择生物标志物的缺失:这是BiTS未来临床成功的关键。自身免疫病具有高度异质性,并非所有患者都适合BiTS治疗。未来的临床试验必须建立一套精准的生物标志物(Biomarker)策略来筛选最可能获益的患者群体。基于本研究的机制,潜在的生物标志物可能包括: 致病T细胞的LAG-3表达水平:只有当致病T细胞高表达LAG-3时,BiTS才能有效发挥作用9。 TCR克隆型与疾病驱动因素:识别出由特定T细胞克隆驱动疾病的患者亚群。 基线信号通路状态:评估患者T细胞内源性TCR信号通路的活化状态,可能有助于预测其对BiTS的敏感性。 开发有效的生物标志物是实现精准医疗、提高T细胞疗法成功率的核心0。 在难治性疾病中的定位:对于那些对现有多种疗法(如TNF-α抑制剂、JAK抑制剂等)均无反应的“难治性类风湿关节炎”(D2T RA)患者,他们往往已经耗尽了标准治疗方案。BiTS这种全新机制的疗法,可能为这部分具有巨大未满足医疗需求的患者带来新的希望。 未来研究方向: 机制的进一步深化:LAG-3的胞内域是否还有其他未知的结合蛋白?它是否会影响除Lck之外的其他激酶?这些问题值得进一步探索。 “邻近诱导”概念的拓展:是否可以将这一设计理念应用于其他共刺激或共抑制受体对,开发出更多用于“打开”或“关闭”免疫反应的工具? 小分子药物的探索:是否有可能开发出能够模拟BiTS功能的小分子“分子胶”,通过口服给药,实现更便捷的治疗? 总而言之,这项研究不仅为我们解锁了免疫检查点LAG-3的深层奥秘,更重要的是,它将深刻的机理洞察转化为一种极具前景的创新疗法。BiTS的成功开发,为攻克T细胞介导的自身免疫性疾病这一顽疾点亮了一盏新的明灯,也为整个药物研发领域带来了关于“空间邻近性”的宝贵启示。
Field Knowledge
· 2025-10-08
阿霉素进核机制全景:从被动扩散到纳米载体主动递送的多重路径
阿霉素进核机制综述 阿霉素自由状态的核转运机制 穿膜与胞质分布:阿霉素(DOX)是一种小分子蒽环类抗癌药物(分子量约543 Da),可以通过被动扩散跨越质膜进入细胞。其两亲性(pH 7.4时对数D值约2.4)使其易于溶解于膜脂质层并穿过细胞膜。早在1985年就有研究表明,蒽环药物主要以被动扩散方式进入细胞,这对药物耐受有重要影响。进入胞质后,游离DOX可在胞质中自由扩散,但因其带正电的性质(pKa≈8.2),在酸性细胞器(如溶酶体)中易被质子化并滞留。这意味着部分DOX可能被内吞体/溶酶体隔离,减少进入细胞核的有效浓度。研究显示DOX会在酸性晚期内体和溶酶体中蓄积。因此,提高内体/溶酶体pH值可增加DOX在细胞内的游离量;例如,添加质子泵抑制剂能提升内体pH,从而增强DOX的核内分布和抗肿瘤效力。 与核孔的相互作用:由于DOX分子相对较小,它可以在不存在主动运输信号的情况下通过核孔复合体(NPC)的选择性屏障被动扩散进入细胞核。NPC由多种核孔蛋白(nucleoporins)构成,在核膜上形成通道结构。NPC的亲水通道充满了FG重复序列的核孔蛋白,形成一个允许小分子自由通过而限制大分子扩散的选择性屏障。一般而言,<~40 kDa的分子可通过NPC自由扩散。DOX的分子量远低于这一阈值,因此无需借助主动核进口受体即可进入细胞核。换言之,DOX不含经典核定位序列(NLS),其进核主要依赖浓度驱动的被动扩散。不过,核孔的通透性也受细胞状态调控:核骨架/核膜结构会影响NPC的开放程度。例如,有研究表明细胞骨架与核骨架的力学联系可调节DOX进入细胞核的通量。当细胞通过外界作用使核纤维和细胞骨架瞬时松弛时,NPC对DOX的通透性增加,核内DOX流入加快;而细胞在软基质上长期培养导致机械张力降低时,核对DOX的通透性下降。由此可见,NPC并非静止通道,其通透性会因细胞机械环境变化而改变。 膜转运载体的作用:除了被动扩散外,某些载体蛋白介导的过程也能影响DOX进入细胞的效率。一些溶质转运蛋白被鉴定为DOX的摄取通道。例如,有机阳离子转运蛋白SLC22A16被证明可介导癌细胞对DOX的主动摄取。在表达SLC22A16的Xenopus卵母细胞中,DOX的摄入呈现可饱和的动力学特征(K_m约5.2 μM)。过表达SLC22A16能增强白血病细胞对DOX的敏感性,说明该转运蛋白增加了细胞内累积的DOX水平。同样,有机阴离子转运多肽OATP1A2(SLCO1A2)也被报道可以转运DOX进入细胞(相反,OATP1B1/1B3对DOX摄取的贡献很小)。这些载体在某些组织或癌细胞中高表达,如SLC22A16在急性白血病和部分实体瘤中过度表达,OATP1A2在肝脏、肿瘤血管等部位存在。有研究发现,在弥漫性大B细胞淋巴瘤患者中,SLC22A16基因缺失与化疗(包括DOX)疗效降低相关,提示缺乏该转运体会削弱DOX对肿瘤的作用。因此,在自然状态下,膜转运蛋白的表达差异会导致不同细胞摄取DOX的能力不同:转运体丰富的细胞可主动摄入更多DOX,而缺乏这些载体的细胞主要依赖缓慢的被动扩散。需要注意的是,与摄入方向相反的是外排转运蛋白(如ATP结合盒(ABC)家族)会将DOX泵出胞外,其中多药耐药泵P-糖蛋白(P-gp/MDR1)对DOX的外排尤为重要。P-gp广泛存在于肿瘤细胞膜和正常组织的屏障细胞中,能主动将DOX泵出细胞,从而降低胞内有效浓度。在P-gp高表达的耐药癌细胞中,DOX的累积量仅相当于野生型敏感细胞的20%左右,进入细胞核的净通量也降低近一半。加入P-gp抑制剂(如维拉帕米、环孢素A或PSC833)可将耐药细胞中DOX的积累提高到敏感细胞的70-80%,同时部分恢复其对药物的敏感性。因此,DOX在未经修饰状态下进入细胞核的效率取决于被动扩散、载体摄取和外排泵三者的净效应。 载药系统对DOX核转运路径的影响 为了提高DOX在肿瘤细胞核内的富集、减少对正常组织的毒性,研究者开发了多种药物递送系统(如脂质体、聚合物纳米粒等)。这些载药系统可以显著改变DOX进入细胞和细胞核的途径。总体而言,纳米载体通过内吞途径进入细胞,然后在细胞内释放DOX,从而实现与游离DOX不同的分布动力学。下面按不同载体类型进行讨论: 脂质体与纳米粒载体:脂质体是最成熟的DOX载药系统之一,典型产品如Doxil(长循环PEG化脂质体DOX)。与游离DOX直接透膜不同,脂质体携带的DOX通常通过胞吞作用进入肿瘤细胞。PEG化脂质体在血流中长循环并富集于肿瘤组织(EPR效应),进入肿瘤后被肿瘤细胞内吞,随后在内体/溶酶体内释放DOX。pH响应型脂质体通过设计在酸性环境(如内体pH≈5-6)下膜结构不稳定,从而更快速释放DOX。一项研究比较了酸敏脂质体(SpHL-DOX)与非酸敏脂质体(nSpHL-DOX)在HeLa细胞中的行为:结果发现,SpHL-DOX进入细胞的速度更快,能够更迅速地在内吞途径中释放DOX,并在细胞核内产生更强的药物累积。SpHL-DOX更高的核累积量引发了更显著的细胞凋亡(如激活caspase-3),而非酸敏脂质体则由于释放缓慢,核内DOX水平较低。此外,加入氯喹或蛋白酶体抑制剂E64d以干预溶酶体功能,可进一步增强酸敏脂质体的细胞毒作用,表明加速内吞体逃逸能提高DOX进入细胞核的效率。与之相反,游离DOX由于无需内吞即可扩散进入细胞,通常在更短时间内达到细胞核。实验显示游离DOX处理6小时内就大量出现在细胞核中(符合其被动扩散机制)。但游离DOX也更易被外排泵清除,缺乏肿瘤选择性。相比之下,脂质体等纳米载体通过尺寸排除效应可以一定程度上规避P-gp的作用(P-gp对300–2000 Da的小分子作用显著,对纳米颗粒不能直接泵出)。有报道指出,脂质体包封可逆转部分耐药性:一方面脂质体使DOX隔离于细胞外,P-gp无法立即作用;另一方面脂质体与P-gp的相互作用可能抑制其泵出功能。因此,脂质体载体不仅提高了肿瘤组织中的药物递送,还通过改变入胞路径而在细胞内提高核输送效率。 聚合物纳米粒与杂化载体:聚合物胶束、聚合物-脂质杂化纳米粒等载体同样通过内吞作用进入细胞,并通过促内体逃逸等机制将DOX释放到胞质,从而增强调入细胞核的药物量。例如,一种聚合物-脂质杂化纳米粒(PLN)将DOX与阴离子聚合物复合后包裹。研究发现,在P-gp过表达的耐药乳腺癌细胞中,DOX-PLN处理可显著提高细胞内DOX的摄入量和滞留时间,远超游离DOX。荧光成像显示,经过DOX-PLN递送,更多药物进入细胞核,并且载体的脂质成分也进入了细胞。内吞途径抑制实验进一步揭示,巨胞饮/吞噬作用是这些纳米粒进入细胞的主要方式。机制上,这种纳米载体通过内吞绕过了细胞膜上P-gp的直接外排,部分DOX伴随载体一起被内吞,躲过了P-gp泵出,然后在细胞内释放并保持更高的滞留。因此,该策略使耐药细胞内的DOX净累积增加,从而克服了一定程度的耐药性。另一项研究将DOX通过共价键偶联到大分子上形成 聚合物药物偶联物,如TAT-PEG-多肽-DOX自组装纳米粒(150 nm左右)。HCT8/ADR耐药结肠癌细胞经这种纳米药物处理后,细胞内DOX摄取显著增加,并绕过P-gp介导的外排,细胞内药物保留时间延长。更重要的是,该纳米药物增强了药物在细胞核的分布,在细胞核内累积更多DOX,强烈抑制了DNA/RNA合成。体内实验也表明此纳米制剂对耐药肿瘤有更好的抑瘤效果。总的来说,聚合物纳米载体通过尺寸效应和亚细胞分布差异,使DOX在癌细胞内(特别是细胞核内)的有效浓度提高,减少了多药泵的影响。 核定位和靶向递送:一些载体通过添加核定位序列(NLS)或其他靶向配体来专门增强DOX的核转运。NLS是富含碱性氨基酸的短肽序列,能够被核输入蛋白(importin α/β)识别并将其携带的货物运入细胞核。尽管DOX本身无NLS,但研究人员将NLS肽偶联到纳米颗粒表面或与DOX结合,从而借助细胞自身的核进口机制提高DOX进核效率。Misra等制备了NLS修饰的PLGA纳米粒载DOX(粒径约226 nm),在乳腺癌MCF-7细胞中显示出显著增强的核靶向输送:与未修饰纳米粒或游离DOX相比,NLS-纳米粒处理使更多DOX进入细胞核,细胞毒性大幅提高(MTT测得IC_50下降,以游离DOX的17.6 μM降至2.3 μM)。共聚焦显微镜观察也证实,NLS修饰纳米粒使DOX在细胞核的定位明显增加。同样,将HIV-TAT等穿膜/核导向肽与DOX或载体结合也是有效策略。Pan等的研究将DOX与TAT-PEG-多肽组装成纳米颗粒,发现其细胞摄取率更高且核内药物分布显著加强,能够克服结肠癌耐药细胞的P-gp介导耐药。此外,配体靶向也是改变DOX分布的途径之一。例如,将DOX与铁转运蛋白转铁蛋白(Tf)偶联可利用癌细胞过度表达的Tf受体(TfR/CD71)进行靶向递送。TfR在多种肿瘤(如乳腺癌)中表达上调,在正常细胞中相对较低。Relecka等的研究表明,DOX–Tf偶联物选择性地增加了乳腺癌细胞内的DOX累积,同时对正常内皮细胞的毒性显著降低:与游离DOX相比,DOX–Tf使正常内皮细胞的存活率提高了3.5倍,而在癌细胞中毒性更强。这说明靶向配体介导的递送既增强了肿瘤细胞对DOX的核内积累,又减少了正常细胞的摄取,体现出对肿瘤的选择性。总之,各类载药系统通过改变细胞摄取途径(如利用受体介导胞吞代替被动扩散)以及促进核靶向,提高了DOX在肿瘤细胞核内的浓度。从分子机制上看,载药系统为DOX打开了新的“进核通路”:内吞途径使其绕过了质膜外排泵、靶向序列使其借用了细胞主动核进口机制。这些设计均旨在提高DOX对肿瘤细胞核的有效打击,同时降低对正常细胞核的毒副作用。 肿瘤细胞与正常细胞机制差异 DOX进核机制在不同类型细胞间可能存在显著差异。癌细胞与正常细胞在膜转运蛋白表达、能量代谢、核膜结构及细胞周期状态等方面的不同,都会影响DOX的入核效率: 膜转运与耐药相关差异:癌细胞常出现膜运输蛋白表达重编程。例如,某些肿瘤高表达DOX摄取载体(如SLC22A16、OATP1A2),这可能使DOX更易进入这些细胞的胞质和细胞核。相反,有些正常组织细胞可能几乎不表达这些载体,主要依靠缓慢的被动扩散来吸收DOX。外排泵P-gp在未经药物压力的正常细胞中表达通常较低,而许多经历过化疗选择压力的肿瘤细胞会过度表达P-gp等MDR蛋白。这意味着在相同浓度下,耐药肿瘤细胞往往将大部分DOX泵出,核内累积浓度偏低;而正常细胞由于缺乏高效泵出,反而可能在核内积累较高DOX水平并发生毒性(典型例子是心肌细胞累积DOX导致心脏毒性)。Shen等比较了P-gp过表达的耐药癌细胞与其亲代敏感细胞中的DOX分布,发现耐药细胞核内DOX含量不到敏感细胞的一半。临床上,人们尝试通过联合P-gp抑制剂或使用非P-gp底物型载药系统来提高癌细胞核内DOX水平,同时希望正常组织因P-gp低而不至于过度蓄积药物。 能量依赖性与代谢差异: 游离DOX主要通过被动过程进出细胞,不直接消耗细胞能量;而载药纳米系统进入细胞常需要能量驱动的内吞。因此,细胞代谢水平的差异会影响DOX的入核效率。肿瘤细胞通常具有活跃的代谢和内吞活动,可加速对纳米药物的摄取;相反,某些正常分化细胞(如静止期细胞)代谢率低,内吞途径不活跃,对载药系统的摄入效率较低。实验上,在4℃低温或存在ATP生成抑制剂(如叠氮化钠)时,内吞作用会被抑制,载药DOX进入细胞核的路径几乎被切断,而游离DOX的被动扩散则受温度影响较小(但膜流动性降低仍会减缓扩散)。Wong等的研究用吞噬抑制剂处理细胞,观察到纳米载DOX的内吞被阻断,核内荧光显著下降。这类能量依赖性的验证说明:肿瘤细胞高度活跃的摄取途径对载药DOX核转运至关重要,而在低代谢状态下(如正常细胞或低温条件),这些途径受限导致核内药物减少。 核膜结构与细胞力学差异:癌细胞的核骨架组成和核孔功能常与正常细胞不同。许多癌细胞下调Lamin A/C等核纤层蛋白,使细胞核更加柔软、变形性更强,以利于侵袭和增殖。这种核结构变化可能影响NPC的密度和通透性。核孔复合体数量在不同细胞中的密度也不同:增殖活跃的细胞通常拥有更多的NPC以满足大量核质交换的需求。侵袭性强的癌细胞系中,提取的裸露细胞核对DOX的通透性彼此存在差异;一般来说,恶性程度高的细胞核膜可能更“漏”,允许DOX更快进入核内。同时,DOX本身对核结构也有影响。一些报道指出DOX处理可导致核膜相关蛋白的变化,如核纤层B1蛋白水平下降,可能使核稳定性降低。核膜的完整性和NPC功能的改变都会反馈影响DOX的核内积累。近期有研究从力学角度揭示:增强细胞核-细胞骨架的联结(例如培养在硬质基质上)会减少DOX进入细胞核,而短时间内破坏肌动蛋白、微管或核纤层则瞬时提高DOX核内摄取。这提示肿瘤细胞经常具有异常的细胞骨架和核骨架互动(例如许多肿瘤细胞展现核膜不规则、核纤层变薄),从而可能天然更容易让DOX进入细胞核。相反,正常细胞核骨架完整,机械张力维持,NPC或许更严格地控制分子进出。机械性“松弛”可视为一种促进药物核进口的方式:比如使用低剂量紫杉醇预处理肿瘤细胞可“力学诱导”核膜更通透,实验证实这使随后加入的DOX核摄取量显著提高。因此,肿瘤细胞独特的核结构/力学生态(柔软的核、异常的核孔功能等)可能是其DOX进核机制与正常细胞不同的内在原因之一。 细胞周期因素:细胞周期影响DOX进核有两个方面:(1) 增殖状态: DOX主要作用于DNA,因此对分裂活跃的细胞毒性更强。S期时细胞核DNA展开供复制,可能增加DOX结合机会;G2/M期时若DOX进入,可干扰有丝分裂。而静止期(G0)细胞核膜稳定、代谢慢,DOX进入和作用都相对减少。(2) 有丝分裂核膜破裂:在细胞进入有丝分裂时,核膜暂时解体,胞质与染色体无屏障交流。在这一阶段,无需通过NPC,DOX即可直接接触到染色质。对于游离DOX而言,这提供了一个“机会窗口”使其大量结合染色体DNA,从而在分裂后子核中保持高浓度。这部分解释了为何DOX对高速分裂细胞有更强杀伤力。不过对于正常细胞而言,很多处于分化静止状态,不经历频繁的核膜崩解,因此DOX只能通过NPC缓慢进入核内,对这些细胞的影响相对较小(但如心肌等非分裂细胞由于缺乏外排机制,仍可能积累DOX导致慢性毒性)。值得注意的是,一些抗癌新策略正是利用细胞周期差异,如同步化肿瘤细胞于有丝分裂期以增加药物摄入,或将正常细胞停滞于特定周期保护它们避免药物伤害。这些都凸显了细胞周期对DOX进核效率的影响。总的来说,肿瘤细胞(特别是高度增殖、具耐药表型的)与正常细胞在DOX进核机制上的差异是多因素叠加的结果,包括膜运输体系、能量代谢、核结构和细胞周期等方面。这些差异为我们提供了靶向肿瘤细胞核转运的策略依据,也提醒我们在研究DOX作用时应区分分析细胞种类。 下表总结了DOX进入细胞核的不同路径及其涉及的关键蛋白或影响机制: DOX进核路径 关键蛋白/机制 参考文献 被动扩散跨膜 (自由DOX) - 无需载体,依赖浓度梯度;- 质膜脂质相容性(疏水-亲水平衡);- 分子量小于NPC排除限制,可自由透过核孔复合体。 载体介导跨膜 - 摄入载体:SLC22A16、OATP1A2等膜转运蛋白,将DOX主动转运入胞质;- 外排载体:P-gp、MRP等,将DOX泵出胞外,减少净内流。 内吞/纳米载体路径 - 胞吞相关:網格蛋白介导内吞、巨胞饮等将载药颗粒吞入;需要ATP能量;- 内体逃逸:pH敏感脂质体膜破裂或聚阳离子破膜,实现DOX从内体/溶酶体释放;- 尺寸效应:纳米粒避免P-gp识别,提高细胞内滞留。 NLS介导核进口 - 核定位序列(NLS):如经典的质朴霉素序列PKKKRKV,可偶联于DOX或载体;- 核输入蛋白:Importin-α/β识别NLS并穿行NPC运输复合物进入核;- 核孔蛋白:NPC中的FG核孔蛋白与importin相互作用开辟通道。 有丝分裂核膜开放 - 核膜消失:细胞进入M期时核膜解体,胞质与核物质混合;- 时机依赖:DOX在此阶段可直接接触染色体,大量嵌入DNA;- 非特异过程:不涉及特定载体蛋白,但仅发生于分裂细胞。 (间接体现周期影响) 表:阿霉素进入细胞核的主要路径及其关键机制。 DOX进核机制验证的关键实验步骤 深入研究和验证上述DOX进核机理,可通过多种细胞和分子实验手段相结合来实现。下面列出若干关键实验设计步骤: 实时荧光成像观测DOX定位:利用DOX本身的红色自发荧光或荧光标记的DOX,结合活细胞激光共聚焦显微镜,动态监测DOX从细胞外、胞质到细胞核的时空分布。例如,可在不同时间点拍摄细胞核内DOX荧光积累的图像,以量化进核速度。必要时也可采用高分辨率共聚焦及三维重构,精确定位DOX是弥散在核质中还是结合于染色质。对于超微结构定位,可用电子显微镜(EM)观察经高浓度DOX处理细胞的超薄切片,在细胞核区域辨识具有高电子密度的DOX-DNA复合物沉积。 利用共聚焦显微镜实时观测阿霉素(DOX)在细胞中的分布变化。(A)示意图:实验设置用于测量DOX进入细胞核的过程。(B)代表性图像:红色荧光为DOX与细胞核DNA嵌合,显示随时间推移DOX在细胞核内的积累。(C)量化曲线:DOX核荧光强度随时间上升,每条曲线代表单个细胞核。(D)示意图:DOX跨膜和入核过程,包括质膜通透(受MDR影响)和经核孔进入核内。本图由BioRender绘制。 干扰核转运相关蛋白:采用基因敲低或敲除技术(siRNA或CRISPR-Cas9)靶向细胞的核转运机器,验证其对DOX进核的影响。首先,可敲低核输入受体如importin-β或importin-α亚基,观察游离DOX和NLS修饰DOX输送的核荧光是否有变化。如果DOX主要经被动扩散,则干扰importin可能无明显影响;而对于NLS介导的载药系统,importin缺失应显著降低DOX核进入量。其次,干扰核孔复合体蛋白(如Nup62、Nup98等),破坏NPC选择性屏障,也能提供线索。例如,用siRNA降低关键FG重复核孔蛋白的表达,测定DOX在核/质的分布变化——若DOX被动通透受限于NPC,那么减少FG网架可能增大NPC孔径、提高DOX核内累积。相反,若干扰核孔蛋白导致核屏障功能紊乱,可能出现大分子漏入核,细胞死亡加快等现象,需要结合对照组进行判读。此外,化学抑制剂也可用来瞬时干扰核转运功能,如小分子INI-43抑制importin-β或小麦胚芽凝集素(WGA)封闭核孔,短时间处理细胞并观察DOX荧光变化,以支持基因敲除的结果。通过这些手段确认特定蛋白在DOX核转运中的作用,可进一步佐证相应机制(如验证NLS纳米粒确实经importin通路入核等)。 荧光共定位分析:将DOX的荧光信号与细胞膜、核膜成分进行共定位观察,可直观了解DOX穿膜和入核过程中的空间关系。具体而言,可选用细胞膜荧光探针(如DiI膜染料)以及标记核膜/NPC的抗体(如抗Nup62、抗核孔复合体蛋白的抗体,结合二抗标记染色)进行染色。然后用共聚焦显微镜观察DOX荧光与这些标记的重叠情况。如在早期时间点,DOX荧光与质膜染料局部重叠,提示DOX可能在膜上聚集或通过膜微区进入;与核孔蛋白信号重叠则表明DOX在通过NPC入口处积累。特别对于载药系统,可标记载体本身(例如在纳米粒上标记一种不同颜色的荧光)并追踪:观察载体荧光先集中于胞质囊泡再逐渐消失、同时DOX荧光转而出现在细胞核的过程。如果将溶酶体标记(LysoTracker绿色)一起观察,可发现DOX载体是否陷于溶酶体,以及加入氯喹等是否促使DOX从这些囊泡释放、进入核内。另外,通过FRET技术(若DOX荧光可以作为给体/受体)也可探测DOX与DNA或膜脂的距离改变,验证DOX是否与核内DNA结合定位。 动力学与能量依赖实验:设计一系列实验以解析DOX进入细胞及细胞核的速度及其对能量的依赖程度。首先,可通过时间梯度实验获取动力学参数:如将细胞暴露于DOX并在0、5、15、30、60分钟等不同时间点终止处理,立即固定细胞后测量核内DOX荧光强度,绘制时间-浓度曲线,估算半饱和时间t½等。这个实验对比自由DOX与载药DOX的曲线,可揭示载药系统是否加快或延缓了进核。其次,进行低温与代谢抑制实验:分别在37℃常温和4℃低温条件下处理细胞,同时设置加入ATP抑制剂(如叠氮化钠、2-脱氧葡萄糖组合)的一组,比较不同条件下DOX核内累积量。例如,观察到4℃时游离DOX仍有一定核进入(可能只是较慢),而载药DOX几乎无法进入细胞核,则证明载药途径高度依赖能量的主动过程。再如,在有氧和无氧条件下比较DOX分布,可验证线粒体能量对内吞的影响。最后,可加入特定途径的抑制剂鉴别内吞途径:如氯化铵/巴菲洛霉素A抑制溶酶体酸化,或用胰蛋白酶消化膜受体,或用药物抑制微管(秋水仙碱)和微丝(肌动蛋白,辣根硫蛋白)等,分别针对网格蛋白介导、巨吞饮、微管依赖运输等途径。通过这些抑制剂的组合,可以判定哪种内吞途径在载药DOX进核中占主要地位。举例来说,若使用菲利平(抑制小窝介导)显著降低DOX核荧光,而酵母多糖(抑制巨吞)无影响,则说明主要经网格蛋白途径。只有将动力学数据与能量依赖性结合分析,才能全面了解不同体系下DOX进核的限制步骤。 正常 vs 肿瘤细胞对比实验:选择代表性的正常细胞系与肿瘤细胞系,在相同条件下比较DOX的核转运效率和机制差异。这可以包括:(1) 定量摄取比较:采用高内涵成像或流式细胞术分别测定正常细胞和癌细胞在给药后一段时间内的核内DOX荧光强度平均值,比较其比例。预期可能看到某些正常细胞核内荧光低于肿瘤细胞(如有外排泵缺失的正常细胞,也可能相反);(2) 转运蛋白表达分析:通过qPCR或蛋白质印迹检测两类细胞中SLC22A16、OATP1A2、P-gp等相关载体的表达量,与功能结果相关联;(3) 共定位及超微结构对比:利用上述荧光共定位手段,观察DOX在正常细胞中是否更多滞留于溶酶体(显示与LysoTracker高度重叠)而在肿瘤细胞中更多逃逸进入核。例如,一项针对DOX在肿瘤和正常心内膜细胞中的研究可能发现,肿瘤细胞溶酶体蓄积较少而正常细胞中DOX荧光大部分局限于周边囊泡。这些差异将有助于阐明为什么DOX对某些正常组织毒性大(可能因为缺乏有效外排/酸化机制)而对某些肿瘤反而作用弱(可能因为内吞隔离或外排过强)。(4) 细胞周期与增殖速率:通过EdU掺入或Ki-67免疫染色确定细胞增殖水平,将其与DOX核摄取量相关联。预期增殖指数高的细胞系(通常肿瘤细胞)对DOX核摄取和损伤更敏感,而静止细胞则抗性更高。综上,通过平行对比实验,可以确认肿瘤细胞在结构和功能上哪些特征促成了DOX更高的核内累积,从而为有针对性地改进药物递送提供依据(例如,对正常细胞高毒性的原因也可由此找到缓解策略)。 综而言之,阿霉素进核机制的研究需要多层次手段验证,从宏观成像到分子干预相结合。通过实时观察、特异干扰和跨细胞种比较等实验,我们能够全面描绘DOX穿膜、入核的路径蓝图,并解释不同类型细胞对这种经典化疗药物截然不同的响应。这不仅加深了对药物作用机理的理解,也将为提高阿霉素疗效、降低副作用的策略(如核靶向递药、克服耐药等)提供科学依据。
Field Knowledge
· 2025-10-08
『”别吃我” vs “吃我”』:细胞世界的攻防战,CD47的双重妙用 本文基本信息 标题:Suppressing or Enhancing Macrophage Engulfment through the Use of CD47 and Related Peptides (通过使用CD47及其相关肽抑制或增强巨噬细胞的吞噬作用) 期刊:Bioconjugate Chemistry Citation: Bioconjugate Chem. 2022, 33, 1989-1995 Corresponding Author: Dennis E. Discher Biophysical Engineering Laboratories and Bioengineering Graduate Group, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States orcid.org/0000-0001-6163-2229; Email: discher@seas.upenn.edu 摘要 外来颗粒和微生物在体内会被巨噬细胞迅速清除,尽管许多关键的摄取机制仍不清楚。“自身”细胞表达CD47,它作为巨噬细胞上SIRPα的抗吞噬配体发挥作用,特别是当促吞噬配体(如抗体)同时展示时。在此,我们综述了CD47及相关的“自身”肽作为巨噬细胞摄取的调节剂。与CD47或源自其SIRPα结合位点的肽共轭的纳米颗粒,可以在体外和体内抑制巨噬细胞的吞噬摄取,在展示CD47的病毒上也发现了类似的现象。因此,作为有效载荷的药物、染料和基因对靶细胞的递送效率得以提高。另一方面,癌细胞表达的CD47使其能够逃避巨噬细胞和免疫监视。这推动了针对CD47-SIRPα的可溶性拮抗剂的开发,从临床上的阻断性抗体到临床前模型中的合生肽。因此,**CD47及其肽正在成为具有双重用途的、抗击疾病的吞噬作用调节剂**。 mindmap root(CD47-SIRPα:细胞吞噬的双向调控枢纽) 背景:巨噬细胞的“敌我识别”挑战 **“吃我”信号**<br/>刺激吞噬外来物 **“别吃我”信号**<br/>保护自身细胞 **核心信号轴**<br/>CD47「自身ID卡」<br/>- SIRPα「ID阅读器」 策略一:抑制吞噬「增强递送」 ::icon(fa fa-paper-plane) **核心思想**<br/>为“友军”穿上“自身”的隐身衣 **应用一:纳米药物** **方法**<br/>将CD47或“自身”肽共轭到纳米颗粒表面 **效果**<br/>延长血液循环时间<br/>增强肿瘤靶向递送 **应用二:病毒载体** **方法**<br/>将CD47整合到慢病毒包膜<br/>将“自身”肽展示在AAV衣壳 **效果**<br/>降低免疫清除<br/>提高基因治疗效率 策略二:增强吞噬「癌症免疫治疗」 ::icon(fa fa-crosshairs) **核心思想**<br/>撕掉癌细胞的“ID卡”,使其暴露给免疫系统 **作用机制**<br/>使用可溶性拮抗剂<br/>阻断癌细胞CD47与巨噬细胞SIRPα的结合 **主要方法** **抗体疗法**<br/>针对CD47或SIRPα的单抗<br/>「如Magrolimab」 **重组蛋白**<br/>可溶性CD47/SIRPα作为诱饵 **“纳米自身”肽**<br/>多价肽高效阻断SIRPα **主要挑战**<br/>“抗原库”效应导致的在靶脱瘤毒性<br/>「如贫血、血小板减少」 结论与展望 **双重用途**<br/>同一靶点,可实现抑制与增强两种相反效果 **未来方向**<br/>优化靶向性,减少副作用<br/>进一步研究促吞噬信号的作用机制 引言 吞噬作用是一种古老而基础的细胞过程,指的是对一个目标的吞食行为。对于变形虫而言,细菌和真菌是其吞噬的目标,这个过程几乎不需要或完全不需要辨别。然而,在动物体内,诸如巨噬细胞之类的吞噬细胞必须识别、攻击并优先吞噬“异己”目标,同时避免伤害健康的“自身”细胞。这些先天性免疫吞噬细胞是宿主抵御各种大小入侵微生物的第一道防线。 吞噬作用由“吃我”信号所激发,这些信号会启动肌动蛋白细胞骨架的重塑,从而驱动巨噬细胞伸出突起以包裹——并随后内化和摧毁——一个“异己”目标。相关的驱动因素范围广泛,从高度特异性的生物分子相互作用物(如蛋白质-蛋白质相互作用),到特异性较低的表面效应物(如电荷、吸附的物质、配体模式),再到物理化学特征(如刚度、形状)。与这些“吃我”通路相对的,是能够抑制巨噬细胞摄取的“别吃我”信号分子。这篇简要综述将聚焦于通过调控特定的“别吃我”信号轴——CD47-SIRPα——来调节巨噬细胞对纳米颗粒、病毒和癌细胞清除方面的最新进展。 巨噬细胞检查点,CD47-SIRPα “自身标记”蛋白CD47是一种普遍表达的整合膜蛋白,它通过与巨噬细胞受体SIRPα相互作用来抑制吞噬摄取。尽管CD24与Siglec-10之间的相互作用可能是另一个潜在的巨噬细胞检查点,但CD47-SIRPα相互作用在许多高等动物中的研究更为透彻且更为保守。对巨噬细胞摄取的抑制作用涉及SIRPα胞质区的免疫受体酪氨酸基抑制基序(ITIM)的磷酸化,并激活磷酸酶SHP-1和SHP-2。 CD47与SIRPα之间的结合相互作用倾向于是物种甚至是品系特异性的,但也存在一些显著的交叉相互作用,例如人的CD47可以与NOD小鼠和猪的SIRPα结合,而猪的CD47也能与人的SIRPα结合。因此,这种受体-配体相互作用的抑制效应是由蛋白质的序列和结构决定的。 巨噬细胞高效清除体内循环异物的能力,常常阻碍了基于纳米颗粒的药物递送,无处不在的巨噬细胞的摄取作用使得药物难以到达预期的靶点,例如癌细胞。这催生了利用CD47来使固体颗粒和病毒更具耐受性,并增加纳米药物和基因靶向递送效率的想法。与此相辅相成的一个目标是拮抗CD47-SIRPα轴以增强吞噬作用,其中,基于抗体的阻断正迅速成为癌症免疫疗法中一个具有临床意义的新增手段,而小分子抑制性肽段的设计也带来了新的可能性。 xxxxxxxx 图1:“自身标记”CD47抑制巨噬细胞的吞噬作用 血清蛋白,如血液中的IgG抗体,会吸附到“外来”颗粒表面或与之特异性结合,刺激吞噬作用。然而,如果巨噬细胞表面的SIRPα与其配体CD47(表达在包括红细胞在内的“自身”细胞上)结合,这种吞噬摄取就会被抑制。 展示CD47和“自身”肽的纳米颗粒通过延迟清除来增强递送 静脉注射的纳米颗粒具有在所有组织和病灶部位循环的潜在优势。然而不幸的是,这类被注射的纳米颗粒通常在数分钟到数小时内就被单核吞噬细胞系统(MPS)清除,特别是肝脏和脾脏中的巨噬细胞。作为对比,新鲜的红细胞(RBCs)在输注后可以循环数周甚至更长时间,但最终同样会被巨噬细胞清除,尤其是在脾脏中。 巨噬细胞识别并清除一个纳米颗粒的具体机制尚不完全清楚。已知的是,血液中的血清蛋白会物理吸附并积聚在所有表面,形成一个“蛋白质冠”,这个冠可以与吞噬细胞的受体作用。其中最显著的是免疫球蛋白G(IgG),它可以结合并激活巨噬细胞的Fc受体(FcRs)。这个过程通常被称为调理作用(opsonization),它代表了生物材料领域长期以来的描述:纯净的化学物质在体内几乎总会被“污染”。聚乙二醇化(PEGylation)是延长纳米颗粒循环的经典方法,它倾向于延迟蛋白质在表面的物理吸附,但清除过程仅仅是被推迟了。 调理作用会导致与巨噬细胞的相互作用,但细胞与材料的相互作用还会受到物理性质的进一步调节,例如刚度、尺寸和曲率(形状),这些都已被证明是影响巨噬细胞清除纳米颗粒的因素。由于巨噬细胞的摄取,纳米颗粒的循环半衰期较短,这一局限性为修饰它们以使其更像“自身”提供了机会(如图2所示)。 将CD47的胞外域(约100个氨基酸)通过生物素化连接到涂有抗生物素蛋白的、细胞大小的聚苯乙烯微球上,足以在微球被抗生物素蛋白IgG调理后,抑制其被吞噬。重要的是,CD47对缺少IgG的微球没有影响。微球实验的成功鼓励了后续的纳米微珠研究,并推动了一种相关的、由21个氨基酸组成的“自身”肽的合成。 体内测试表明,CD47和“自身”肽都能通过延迟脾脏巨噬细胞的清除,来增加纳米微珠在小鼠体内的循环半衰期,从而极大地增强了肿瘤成像和药物在肿瘤部位的递送效率。随后,一个独立的实验室将“自身”肽连接到氧化石墨烯纳米片上,报告了相似的结果,并得出结论:“自身”肽比PEGylation更有效。其他实验室的研究也表明,用重组CD47或“自身”肽对纳米材料进行功能化,通常能延长循环时间、抑制清除并改善治疗效果。 另一项应用将“自身”肽连接到纳米脂质体上,与对照组纳米脂质体不同,研究者发现它们能够饱和并“钝化”肝脏巨噬细胞,从而增加了后续注射的其他纳米颗粒的循环时间和功效。在这些不同的纳米颗粒上,血清中起调理作用的IgG的沉积是否在结果中扮演了角色,通常尚不明确。尽管如此,各项研究都凸显了将CD47或更短的“自身”肽偶联到各种纳米材料上在多种应用中的实用价值。 xxxxxxxxxx 图2:CD47肽延长循环并增加靶向递送 载药纳米颗粒(左)和携带基因的慢病毒(右)由于被吞噬细胞摄取,其疗效有限。将CD47及相关的“自身”肽连接到纳米颗粒和病毒表面,它们可以与巨噬细胞上的SIRPα结合,帮助将其识别为“自身”,从而延长循环时间并增强靶向递送。 展示“自身”信号的病毒能够抑制吞噬作用并增强基因递送 基于病毒的基因递送已广泛用于临床,例如疫苗(如新冠病毒的刺突蛋白)和细胞的离体工程改造(如CAR-T细胞)。在静脉注射以实现靶向基因递送的尝试中,慢病毒(Lentivirus)和腺相关病毒(AAV)载体最为常见,但巨噬细胞同样会被激活以清除这些“天然的纳米颗粒”,这可能导致病毒诱导的炎症反应。许多团队尝试通过偶联合成聚合物来抑制单核吞噬细胞系统(MPS)介导的病毒清除,以期最小化调理作用;然而,这类修饰的空间位阻会妨碍病毒与目标靶点结合所必需的关键蛋白质相互作用。 慢病毒通常是在细胞系胞吐后收获的,因此,通过适当改造的细胞系过表达膜蛋白CD47,原则上可以产生展示CD47的慢病毒。已有两项独立研究确实成功生成了CD47-慢病毒,并证实其能减少与巨噬细胞的相互作用并改善基因递送效果。 第一项研究用对照组或CD47-慢病毒递送红色荧光蛋白(RFP)到分化的人类巨噬细胞培养物中,结果显示: CD47-慢病毒的转导效率比对照组低约3倍。 表达SIRPα的A549肺腺癌细胞优先被CD47-慢病毒转导。(不是抑制吗?) 后一个结果表明,SIRPα充当了CD47介导的附着和感染的“停靠受体”。在体内也观察到了相似结果:在A549肿瘤中,使用CD47-慢病毒的转基因表达水平更高,而肝脏和脾脏巨噬细胞中的表达则相对于对照组显著降低。研究还通过基于抗体的SIRPα相互作用抑制实验来验证其特异性。 第二项研究利用人源CD47来提高慢病毒在肝脏基因转移的效率。在确定了肝脏巨噬细胞会清除静脉注射的慢病毒后,研究发现CD47-慢病毒增加了对肝细胞的基因转移,同时减少了对巨噬细胞的转移。这些实验在两种小鼠模型中进行:一种是能表达与人源CD47结合的SIRPα的NOD小鼠,另一种是亲和力较弱的C57BL/6小鼠。结果显示,CD47-慢病毒在C57BL/6小鼠中的清除率更高。其安全性和有效性在与人类CD47和SIRPα序列同源性更高的非人灵长类动物中得到了进一步证实。 这些研究表明,展示CD47能够保护像慢病毒这样的有膜病毒免受巨噬细胞的清除,从而增强基因转移疗法的功效。与慢病毒类似,在临床上具有重要意义的腺相关病毒载体(AAV)上展示“自身”肽,在体外也导致了AAV的吞噬易感性降低。由于AAV没有膜包被,研究人员将“自身”肽直接引入AAV2的衣壳蛋白中,并用甘氨酸-丝氨酸接头(linker)连接以确保衣壳稳定并最小化病毒滴度的损失。这种插入对转导效率几乎没有影响,但与对照组AAV2相比,在人类巨噬细胞中,病毒的摄取量降低了多达10倍。当用抗SIRPα抗体进行阻断后,这种差异再次消失。 AAV的尺寸仅为20纳米,而慢病毒约为100纳米。鉴于CD47-SIRPα是吞噬摄取的特异性抑制剂,而非内吞作用的抑制剂,迄今为止CD47偶联病毒的研究结果,共同凸显了巨噬细胞对纳米颗粒进行吞噬的高效性。在细胞生物学文献中,吞噬作用常被认为仅与较大的实体(颗粒、凋亡细胞或微生物)相关,但早期对尺寸差异巨大的颗粒进行的实验,并未充分考虑到颗粒浮力的差异和其他尺寸效应。如果很少有小颗粒沉降,那么被摄取的自然就少。然而,浮力在体内并不重要。在上述纳米颗粒和病毒的研究中,肝脏和脾脏中的巨噬细胞之所以突出,是因为这些巨噬细胞排列在这些组织的血管壁上,从而能够直接、即时地接触到静脉注射的颗粒。尽管如此,巨噬细胞仍然存在于所有组织中,并且通常是肿瘤或穿刺/损伤部位等病灶处的主要细胞类型。摄取途径对于最终结果也很重要:例如,吞噬体(phagosomes)比内涵体(endosomes)对货物的氧化和破坏性更强。所有这些因素都对数十亿剂作为疫苗注射的病毒(例如,强生或牛津-阿斯利康新冠疫苗中由腺病毒递送的DNA)具有深远影响。 可溶性CD47-SIRPα拮抗剂增强吞噬作用 CD47在细胞表面普遍表达,但早在几十年前,人们不仅记录了CD47在卵巢癌中的过表达,而且在CD47序列被测定之前,用于肿瘤成像的抗体靶向最终也被证明能抑制吞噬作用。随后,针对其他癌症的CD47抗体靶向研究也相继展开,并有证据表明在人类肿瘤异种移植模型中存在治疗窗口,尽管最初尚不清楚这种IgG是抑制了吞噬作用,还是激活了FcR驱动的吞噬作用,或是两者兼有。 此外,一项关于在同系小鼠肿瘤模型中使用抗小鼠CD47治疗的研究,在后来的重复性验证中显示,抗CD47单药治疗没有任何抗肿瘤效果的迹象,反而显示出贫血副作用。这种单药治疗的负面结果在很大程度上也反映在了临床抗癌实践中,并且这与CD47基因敲除小鼠几乎正常(仅有极小缺陷且无明显贫血)的事实似乎是一致的。后一项由免疫学实验室得出的观察结果,一度引起了血液学家们对CD47所谓“自身标记”功能的极大怀疑。 与单药治疗形成对比的是,将拮抗CD47-SIRPα巨噬细胞检查点与一个“吃我”信号相结合的策略展现出了巨大的潜力,并催生了对可溶性拮抗剂研究的爆炸性增长。这些拮抗剂的范围从临床上的各种IgG设计和重组蛋白,到小分子肽,它们都作为可能的药物,对多种血液和实体恶性肿瘤显示出不同程度的疗效。 目前最先进的抗CD47疗法是一种名为magrolimab(或Hu5f9-G4)的人源化IgG4单克隆抗体,它能结合CD47并抑制其与SIRPα的结合,但由于IgG4与巨噬细胞FcR的亲和力较弱,因此不会主动激活巨噬细胞。然而,CD47在体内几乎所有细胞上的表达构成了一个“抗原库”(antigen sink),导致magrolimab等CD47靶向抑制剂的非特异性结合,从而引发了不可避免的在靶毒性,例如贫血和血小板减少症。为解决这一安全问题,正在进行的努力包括开发对CD47具有强结合力但对人红细胞亲和力低的纳米抗体。 靶向SIRPα受体可能更为安全,因为其表达更具限制性,尽管SIRPα的表达不仅限于髓系细胞,也延伸到了如上皮细胞等其他细胞。一些研究确实表明,抗SIRPα阻断与抗CD47同样有效,但能维持安全的血液学指标。一项工程化巨噬细胞的研究进一步证明,将SIRPα阻断与用肿瘤靶向IgG预激活FcR相结合,在缩小已形成的肿瘤方面是有效的。 最近,研究人员基于CD47上与SIRPα结合的β-发夹环结构,开发出了多价的8个氨基酸组成的“纳米自身”(nano-Self)拮抗剂。这些肽的变体能够有效阻断CD47-SIRPα的相互作用,并在低至5 nM的浓度下,增强人类巨噬细胞对经抗体调理的人红白血病细胞的内化。同一项研究中的其他观察结果,进一步证实了巨噬细胞上的CD47能与同一细胞上的SIRPα发生顺式相互作用,传递一种自抑制信号,这与早前的观察结果一致。 然而,并非所有将可溶性CD47多肽添加到培养的巨噬细胞中的研究都显示出吞噬作用的增强。奇怪的是,一项早期的研究报告称,用细菌表达的人源CD47蛋白,在体外反而降低了小鼠巨噬细胞对胶体乳液的吞噬作用。后续其他团队的研究表明,CD47的相互作用能力需要一种细菌所缺乏的翻译后N-末端修饰才能得到改善,并且特定的人-鼠CD47-SIRPα相互作用本身就特别弱,此外,可能还需要对目标进行IgG调理才能揭示出CD47-SIRPα阻断的效果。所有这些仍然是该领域未来发展中需要重点考虑的因素。 图3:用于免疫治疗的SIRPα-CD47可溶性拮抗剂 肿瘤细胞表达巨噬细胞检查点CD47以抑制吞噬作用。单独的IgG抗体调理作用因CD47的“别吃我”信号而不足以引发有效的吞噬。但多种策略可以拮抗这种抑制。目前至少有三种免疫治疗策略正在临床前和临床研究中进行:抗CD47或抗SIRPα的抗体、作为抑制剂的可溶性蛋白版本,以及相关的“自身”肽拮抗剂。小分子(绿色三角)最终可能被开发出来抑制CD47的转录,但仍需要一个“吃我”信号。 结论 SIRPα-CD47轴正成为一个在递送和治疗等多种应用中越来越有吸引力的靶点。展示CD47或相关肽段的纳米颗粒和病毒被巨噬细胞识别为“自身”,从而延迟了这些颗粒的吞噬,延长了循环时间,并增加了染料、药物和基因的靶向递送。未来需要进一步的研究来理解这些纳米颗粒和病毒上的促吞噬信号(即调理作用和蛋白质冠的形成)。 与此同时,旨在增强吞噬作用(特别是对癌细胞的吞噬)的可溶性拮ក抗剂正在持续开发和探索中,这展示了该研究领域发展的双重用途。在全身性注射拮抗剂(如抗CD47 IgG)后,限制其脱靶效应仍然是挑战。至少有一项近期的有趣尝试是使用纳米颗粒来同时阻断CD47并调理癌细胞,但这当然要求纳米颗粒既能躲避巨噬细胞,又能接触到肿瘤细胞。小尺寸有助于渗透到实体瘤中,而最近合成的一种紧凑的环状“纳米自身”肽已被证明能在体外增强原代巨噬细胞对经单抗靶向的黑色素瘤的吞噬,为体内的疗效测试奠定了基础。 CD47-SIRPα信号轴的双向应用策略对比 对比维度 策略一:增强递送 (模拟“别吃我”信号) 策略二:增强吞噬 (阻断“别吃我”信号) 核心思想 为药物/基因载体穿上“自身”隐身衣,使其逃避免疫系统清除。 撕掉癌细胞的“自身”伪装,使其暴露给免疫系统攻击。 关键分子工具 CD47蛋白或其衍生的“自身”肽,直接偶联在纳米载体(如纳米颗粒、病毒)表面。 可溶性拮抗剂,如:抗CD47/SIRPα抗体 (Magrolimab);重组蛋白/“纳米自身”拮抗肽 作用机制 载体表面的“自身”肽激活巨噬细胞SIRPα的抑制性信号通路,从而抑制对载体的吞噬作用。 拮抗剂阻断癌细胞CD47与巨噬细胞SIRPα的结合,从而解除对癌细胞吞噬作用的抑制。 主要应用领域 纳米药物递送(延长循环,增强肿瘤靶向)基因治疗(保护病毒载体,提高转导效率) 癌症免疫治疗(特别是与“吃我”信号药物联用,治疗血液瘤和实体瘤) 主要挑战与副作用 需确保修饰不影响载体自身的功能;蛋白质冠等其他清除机制的影响尚不明确。 严重的在靶脱瘤毒性 (on-target, off-tumor toxicity)因健康细胞(尤其是红细胞)也表达CD47,导致贫血、血小板减少等副作用。
Field Knowledge
<
>
Touch background to close