← All hypotheses

声遗传学STRC假设

DAY 3

用声音治疗听力损失。俄语中有句话叫klin klinom:以楔驱楔。问题本身成为治疗手段。

每种基因治疗载体都需要一个启动子:告诉细胞何时合成蛋白质的开关。标准方法是持续开启。我们的方案:将基因置于一个响应声音的开关后面。毛细胞已经具备这种机制。声音使静纤毛弯曲,钙离子涌入,信号级联启动。我们加以利用。有声音 = 基因开启。安静 = 基因关闭。孩子佩戴他已有的助听器,无需特殊频率,只是日常声音。

蛋白质会在夜间消失吗?不会。立体纤毛蛋白的半衰期约为30天。达到正常水平50%大约需要13小时的助听器使用时间。睡眠不会重置任何东西。

信号级联

Sound
Hearing aid amplifies everyday sounds to 60-80 dB
No special frequencies needed
Stereocilia
MET channels open at tips
~134 channels per cell
15% Ca²⁺ permeability
Ca²⁺ signal
Apical compartment
70 nM → 500-900 nM
V = 0.05 pL (tiny space)
activates
Calcineurin
Ca²⁺-dependent phosphatase
CnA + CnB heterodimer
Kd = 500 nM, Hill n = 4
NFAT
Transcription factor
Dephosphorylated → NLS exposed
Enters nucleus in minutes
enters nucleus, binds promoter
6×NFAT promoter
6 binding sites + minimal promoter
Sharp on/off switch (Hill ~4)
62× induction, zero leakage (Wu 2023)
mini-STRC
1076 aa (residues 700-1775)
Traffics to stereocilia tips
t½ ~30 days. Accumulates over days.
Hearing
Top connectors restored
Bundle moves as unit again
50% in 13h, full in 72h
SLEEP / SILENCE
No sound → no Ca²⁺ → promoter OFF → protein stable (t½ 30 days)
POSITIVE FEEDBACK
As hearing improves → better signal → more protein → better hearing

单AAV构建体

ITR
145bp
Promoter
6xNFAT
~300bp
Transgene
mini-STRC CDS
3228bp
PolyA
bGH
250bp
ITR
145bp
构建体总长度: 4,068 bp AAV包装限制: 4,700 bp 安全余量: 632 bp

Critical Assumptions & Open Questions

HONEST DISCLOSURE

The ODE model is mathematically correct: equations are sound, parameters are from peer-reviewed literature, and the system dynamics are thoroughly analyzed. However, several biological assumptions have not been validated in cochlear outer hair cells. A literature review (April 2026) identified the following critical gaps.

UNVERIFIED
"Ca²⁺ from MET channels activates calcineurin in OHC soma"
Multiple studies show extreme Ca²⁺ compartmentalization in OHC stereocilia. Millimolar calcium buffers, PMCA2 pumps, a mitochondrial barrier below the cuticular plate, and 4-6 mM oncomodulin in the cell body all prevent MET-channel calcium from reaching the soma. Resting somatic Ca²⁺: ~91 nM. Calcineurin Kd: ~500 nM. No published evidence of calcineurin activation during normal hearing.
Bhatt et al. 2000; Ricci et al. 1998; Hackney et al. 2005
NOT OBSERVED
"NFAT translocates to nucleus during normal sound exposure"
Every published observation of NFATc4 nuclear translocation in OHCs is in the context of pathological injury: noise trauma or aminoglycoside exposure. NFATc4 knockout mice develop normal cochlear structure and have normal hearing. NFAT is dispensable for OHC function.
Frontieres Immunol 2019; Front Cell Dev Bio 2021
CONTRADICTED
"NFAT activation in OHCs is safe"
NFATc4 downstream targets in OHCs include TNF, Caspase-8, and Caspase-3, all apoptosis effectors. NFAT inhibitors (11R-VIVIT, pasireotide) protect hair cells from death. Activating NFAT in OHCs may trigger cell death, not therapeutic gene expression.
Pubmed 27787949; Pubmed 31755556
UNTESTED
"6xNFAT promoter functions in hair cells"
The 6xNFAT promoter (Wu et al. 2023) was validated in HEK293T cells with focused ultrasound. It has never been tested in any cochlear cell type. No synthetic NFAT-responsive promoter has been tested in the inner ear.
Wu et al. 2023
The core problem: Ca²⁺ from MET channels stays in stereocilia under physiological conditions. It never reaches the soma at concentrations sufficient to activate calcineurin (~500 nM). The model's V_cell = 0.05 pL (apical compartment) assumes calcium reaches a signaling volume where calcineurin resides. Published data suggests it does not.

Alternative: AC1/cAMP/CREB Pathway

PROPOSED

A more biologically plausible activity-dependent pathway exists. Adenylyl cyclase 1 (AC1/ADCY1) is a calcium-calmodulin-activated enzyme located directly in stereocilia, where MET-channel calcium is abundant. It converts Ca²⁺ signals to cAMP, a small molecule that diffuses freely to the soma, bypassing the calcium compartmentalization barrier entirely.

MET Ca²⁺ CaM AC1 cAMP↑ PKA CREB-P CRE→STRC
NFAT pathway AC1/CREB pathway
Ca²⁺ sensor location Soma (Ca²⁺ doesn't reach) Stereocilia (where Ca²⁺ is)
Second messenger Ca²⁺ itself (heavily buffered) cAMP (freely diffusing)
Pathway in OHC Only activated by damage AC1 mutations cause deafness
Safety Triggers apoptosis (TNF, Caspases) Neuroplasticity pathway
Promoter element 6xNFAT-RE (never tested in ear) CRE (TGACGTCA, well characterized)
Tested in hair cells No No

Why AC1 is credible

ADCY1 in stereocilia: Immunolocalization confirms AC1 protein in hair cell stereocilia, exactly where MET-channel Ca²⁺ is highest.
ADCY1 mutations cause deafness: Recessive hearing impairment in humans and zebrafish. AC1 is functionally essential for hearing, unlike NFATc4.
cAMP diffuses freely: Unlike Ca²⁺ (trapped by millimolar buffers), cAMP passes through the cuticular plate without barrier. Reaches nucleus in seconds.
CREB is safe: Standard neuroplasticity pathway. No link to OHC apoptosis.

What remains unproven

1. Does AC1 produce enough cAMP during normal hearing to activate PKA/CREB?
2. Does cAMP from stereocilia reach the nucleus at sufficient concentration?
3. Would a CRE-responsive promoter fire in OHCs? (Never tested)
4. cAMP/CREB responds to many signals, not just sound. Specificity is a concern.
5. All ODE model parameters would need re-derivation for the AC1/CREB cascade.

Speculative: DFNB16-Specific Calcium Leak

In Strc-/- mice, stereocilia lose horizontal top connectors by P15, bundle stiffness drops 60-74%, and stereocilia progressively disconnect. It is possible that this structural disorganization disrupts calcium compartmentalization, allowing MET-channel Ca²⁺ to leak into the soma, making the NFAT pathway viable only in DFNB16 cells.

If true, this would create an elegant self-dosing mechanism: the disease activates the promoter, the cure silences it. But no one has measured somatic Ca²⁺ in Strc-/- OHCs. This remains pure speculation. It is testable by calcium imaging in cochlear explants from stereocilin-null mice.

支持文献

[1] Wu et al.(2023)。Sonogenetic control of multiplexed genome regulation and base editing. Nature Communications 14:6811。62倍诱导,三周内零泄漏。 doi
[2] Tomida et al.(2003)。NFAT functions as a working memory of Ca²⁺ signals in decoding Ca²⁺ dynamics. EMBO J 22:3825-3832。NFAT核转位动力学。doi
[3] Fettiplace & Kim(2014)。The physiology of mechanoelectrical transduction channels in hearing. Physiological Reviews 94:951-986。MET通道特性。doi
[4] Iranfar et al.(2026)。Dual-vector gene therapy restores cochlear amplification and auditory sensitivity in a mouse model of DFNB16 hearing loss. Science Advances。小鼠STRC基因治疗首例。 PMC
[5] Li, Mu & Yan (2026). Dissecting the Black Box of AlphaFold in Protein–Protein Complex Assembly. bioRxiv. AF3 assembles complexes via monomer geometry + interface pattern matching, not coevolution. Explains why NFAT-Calcineurin (ipTM 0.73) scores well (canonical interface) while STRC partners score low (non-canonical). doi

这是一个计算假设。每个单独组件均已得到验证(NFAT启动子、MET通道生物物理学、mini-STRC包装)。目前尚无人将它们结合用于内耳基因治疗。

ODE模型:钙信号到蛋白质的动力学

DAY 3 模拟结果

空话廉价。我们建立了一个数学模型来验证这一级联是否真的能产生足够的蛋白质。五个微分方程,每个参数均来自已发表文献,测试了四种情景。问题是:佩戴助听器的孩子能否仅通过声音产生治疗性水平的立体纤毛蛋白?

助听器周期(16小时/8小时) 现实情景
立体纤毛蛋白分子数 29,571
目标值: 15,000 197%
Peak Ca²⁺: 689 nM
Peak NFAT(nuclear): 37.3%
达到10%功能的时间: 7.5h
达到50%功能的时间: 13.0h
治疗课程(每天2小时,85 dB) 85 dB
立体纤毛蛋白分子数 29,571
目标值: 15,000 197%
Peak Ca²⁺: 886 nM
Peak NFAT(nuclear): 40.0%
恒定70 dB 70 dB
立体纤毛蛋白分子数 29,733
目标值: 15,000 198%
静默(对照) 对照
立体纤毛蛋白分子数 1,023
目标值: 15,000 6.8%
Peak Ca²⁺: 171 nM
Peak NFAT(nuclear): 1.3%

关键发现:29倍动态范围,13小时内达到治疗水平

在现实的助听器使用时间表下(70 dB开启16小时,睡眠8小时),模型预测72小时后每个OHC中有29,571个立体纤毛蛋白分子(目标:15,000个)。在安静环境中,仅积累1,023个分子(6.8%)。这在声音激活和静默状态之间提供了29倍动态范围。仅需13小时助听器使用即可达到50%治疗阈值。系统自我调节:蛋白质在静纤毛上的可用结合位点处达到饱和,防止过度表达。

Bifurcation Analysis: Self-Dosing Window

FINDING

The system saturates at ~20 dB. Between 22 dB (quiet bedroom) and 70 dB (hearing aid), protein output is identical: ~29,750 molecules. The 29x dynamic range exists only between 0 and 14 dB, a range humans rarely experience. Self-dosing via hearing aid gain modulation does not work with the current 6xNFAT promoter.

0 dB
16%
8 dB
52%
14 dB
123%
22 dB
194%
40 dB
198%
70 dB
198%
100 dB
198%
Steepest transition
14 dB
2,138 molecules/dB
50% therapeutic
8 dB
Near silence
100% therapeutic
14 dB
Whisper at 5m

Eigenvalue Stability

All 5 eigenvalues at 70 dB have negative real parts. System is asymptotically stable. No oscillations, no chaos. Five timescales: Ca²⁺ (seconds), NFAT cycling (minutes), mRNA (hours), protein (days).

Hysteresis

Path-dependence detected due to slow protein degradation (t½ = 30 days). Kinetic memory, not bistability. One steady state per parameter set. System is monostable.

Solution: 9xNFAT-weak Promoter

NEW

Parameter optimization across 245 combinations found that increasing the NFAT activation threshold (0.05→0.30) and cooperativity (Hill 4→8) moves the switch to 40-50 dB. Biologically: use 9 NFAT binding sites with reduced individual affinity (A→G at position 4 of each site). One nucleotide change, repeated 9 times, transforms the therapy from "always ON" to "hearing-aid controlled."

Sound 6xNFAT (current) 9xNFAT-weak
0 dB 11.6% 0.1%
20 dB 179.7% 0.1%
40 dB 198.3% 5.7%
45 dB ~198% 51.2%
50 dB 198.3% 163.3%
70 dB 198.3% 198.1%
// Single nucleotide substitution
Strong: AGGAAAAT Weak:  AGGGAAAT (position 4: A→G)
// Full promoter: 9× [AGGAGAAT] + spacers + TATA = 161 bp
// AAV construct: 4,247 bp total (453 bp margin)
9xNFAT-weak CANDIDATE
Sites:
Sequence: AGGAGAAT
Hill: 8 | Threshold: ~0.30
Switch at 40-50 dB
9xNFAT-mixed FALLBACK
Sites:
Sequence: 5×strong + 4×weak
Hill: 6-7 | Threshold: ~0.15
Switch at 25-35 dB
6xNFAT-strong BASELINE
Sites:
Sequence: AGGAAAAT
Hill: 4 | Threshold: ~0.05
Current model. Saturates at 14 dB

Literature Support

Yamamoto et al. 2018 — Direct 3x/6x/9x NFAT-RE comparison. 9x gave best fold-change (ON/OFF ratio). Supports Hill=8 prediction.
Crabtree 2003, Genes Dev — GGAGA documented as suboptimal NFAT motif with 3-5x reduced affinity vs consensus GGAAA.
CD3γ mutagenesis (JBC 2002) — A→G at position 4 of GGAAAA abrogates binding at individual sites. With 9 cooperative sites, partial binding is retained.

Safety & Stability Summary

6 / 6 PASS
SAFE
Protein over-expression
Hard cap at 30,000 molecules (2× target). Even at 200× fold induction: 199.5% of target. Saturation factor (1 − protein/max) prevents runaway.
SAFE
Promoter strength variation
Fold induction tested from 10× to 200×. All scenarios stay below 200% of target. The saturation cap dominates over promoter strength.
SAFE
Extreme sound exposure
At 120 dB sustained, protein unchanged. Ca²⁺ rises to ~1,155 nM (mild stress). But this is a property of all hair cells, not the construct.
MONITOR
Ca²⁺ toxicity
If buffer ratio drops below 30 (impaired buffering), Ca²⁺ exceeds 1,000 nM at therapeutic levels. Not caused by construct, but pre-existing cochlear health affects safety.
NONE
Oscillations / chaos
All eigenvalues negative. CV < 1% for all variables. No limit cycles. The system behaves like a thermostat, not a bomb.
NONE
Bistability
One steady state per parameter set. Hysteresis exists but is kinetic (slow protein decay), not structural. No dangerous hidden states.

Full bifurcation analysis: Jacobian eigenvalues at 0/70/90 dB, two-parameter sweep (Kd_CaN × NFAT_threshold, 90 points), hysteresis check (ramp up vs down over 72h), and promoter fold-induction sensitivity (10×-200×). All passed. Code: bifurcation_analysis.py

Closed-Loop Model: True Self-Dosing

NEW

Previous ODE models used fixed sound profiles. This model closes the feedback loop: protein accumulates → hearing improves → audiologist reduces hearing aid gain → less MET activation → less STRC production → equilibrium. The hearing aid becomes a true dosing device.

Sound Ca²⁺ NFAT Stereocilin Hearing ↑ HA gain ↓ Sound ↓

Clinical Timeline (9xNFAT-weak, 50 dB ambient)

Day 0.5
21%
Day 1
108%
Day 3
161%
Week 1
177%
Week 2
183%
Month 3
184%

Day 3: Hearing improves enough that HA adds zero gain. System reaches equilibrium at 27,617 molecules (184% of target), 3.1 dB residual hearing loss, maintained by ambient sound alone.

Stress Tests: 6xNFAT-strong vs 9xNFAT-weak

The promoters behave identically under normal conditions. The difference appears when sound is removed. This is the test of real self-dosing.

HA lost 14 days 6x: -0.3% 9x: -26.8%

9x responds to silence. 6x doesn't.

Hospital 30 days (25 dB) 6x: -1.0% 9x: -49.4%

One half-life. Full recovery in 2 weeks.

Min dB for therapy (no HA) 6x: 0 dB 9x: 45 dB

9x needs real sound. 6x is constitutive.

Built-in resilience: Stereocilin t½ = 30 days. Even 2 weeks without any sound loses only ~27% of protein. Recovery: 98% within 2 weeks of HA restoration. Self-dosing is a slow process (weeks), not fast (hours). The long protein half-life buffers against daily sound fluctuations.

模型参数(来自文献)

参数 数值 来源
MET channel conductance 150 pS Beurg et al. 2006
Channels per bundle 134 Fettiplace 2017
Endolymphatic potential +80 mV Standard
Calcineurin Kd (Ca²⁺) 500 nM Stemmer & Klee 1994
NFAT nuclear import t½ ~2 min Tomida et al. 2003
Promoter fold induction up to 62x Wu et al. 2023
Promoter leakage Zero (3 weeks) Wu et al. 2023
Apical compartment volume 0.05 pL Lumpkin & Bhatt 2001

灵敏度分析

结果有多脆弱?我们将每个参数变化±50%,找出最关键的因素。如果一个数字改变时模型失效,该参数就需要最仔细的实验测量。

k_transcription_max 1.50
k_translation 1.50
n_channels 0.46
Kd_CaN -0.47
buffer_ratio -0.46

转录和翻译速率占主导地位(灵敏度1.50)。这意味着最重要的实验是专门表征毛细胞中6xNFAT启动子的强度。通道数量和钙调磷酸酶亲和力的影响较小:系统对这些参数的生物变异具有鲁棒性。

可复现:完整源代码

完整ODE模型以Python脚本形式提供。依赖项:numpy, scipy。自行运行以复现这些结果或修改参数。

在GitHub上查看:ode_model.py
← All hypotheses