← All hypotheses

Sonogenetic STRC Hypothesis

DAY 3

Treating hearing loss with sound. In Russian we say klin klinom: driving a wedge out with a wedge. The problem itself becomes the cure.

Every gene therapy vector needs a promoter: a switch that tells the cell when to make the protein. Standard approach is always-on. Ours: put the gene behind a switch that responds to sound. Hair cells already have this machinery. Sound bends stereocilia, calcium floods in, a signaling cascade fires. We hijack it. Sound in = gene on. Silence = gene off. The child wears the hearing aid he already has. No special frequencies. Just everyday sounds.

Does the protein vanish at night? No. Stereocilin has a half-life of ~30 days. It takes ~13 hours of hearing aid use to hit 50% of normal levels. Sleep doesn't reset anything.

The cascade

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

Single-AAV construct

ITR
145bp
Promoter
6xNFAT
~300bp
Transgene
mini-STRC CDS
3228bp
PolyA
bGH
250bp
ITR
145bp
Total construct: 4,068 bp AAV packaging limit: 4,700 bp Safety margin: 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.

Supporting literature

[1] Wu et al. (2023). Sonogenetic control of multiplexed genome regulation and base editing. Nature Communications 14:6811. 62-fold induction, zero leakage over 3 weeks. 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 nuclear translocation kinetics. doi
[3] Fettiplace & Kim (2014). The physiology of mechanoelectrical transduction channels in hearing. Physiological Reviews 94:951-986. MET channel properties. 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. First STRC gene therapy in mice. 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

This is a computational hypothesis. Each individual component is proven (NFAT promoters, MET channel biophysics, mini-STRC packaging). Nobody has combined them for inner ear gene therapy yet.

ODE Model: Calcium-to-Protein Dynamics

DAY 3 SIMULATED

Words are cheap. We built a mathematical model to check whether this cascade actually produces enough protein. Five differential equations, every parameter from published literature, four scenarios tested. The question: does a hearing-aid-wearing child produce therapeutic levels of stereocilin through sound alone?

Hearing Aid Cycle (16h/8h) REALISTIC
Stereocilin molecules 29,571
Target: 15,000 197%
Peak Ca²⁺: 689 nM
Peak NFAT(nuclear): 37.3%
Time to 10% function: 7.5h
Time to 50% function: 13.0h
Therapy Sessions (2h/day 85 dB) 85 dB
Stereocilin molecules 29,571
Target: 15,000 197%
Peak Ca²⁺: 886 nM
Peak NFAT(nuclear): 40.0%
Constant 70 dB 70 dB
Stereocilin molecules 29,733
Target: 15,000 198%
Silence (control) CONTROL
Stereocilin molecules 1,023
Target: 15,000 6.8%
Peak Ca²⁺: 171 nM
Peak NFAT(nuclear): 1.3%

Key finding: 29x dynamic range, therapeutic levels in 13 hours

With a realistic hearing aid schedule (16 hours ON at 70 dB, 8 hours sleep), the model predicts 29,571 stereocilin molecules per OHC after 72 hours (target: 15,000). In silence, only 1,023 molecules accumulate (6.8%). This gives a 29-fold dynamic range between sound-activated and silent states. The 50% therapeutic threshold is reached in just 13 hours of hearing aid use. The system self-regulates: protein saturates at the available binding sites on stereocilia, preventing overexpression.

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.

Model parameters (from literature)

Parameter Value Source
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

Sensitivity analysis

How fragile is the result? We varied each parameter by ±50% to find what matters most. If the model breaks when one number changes, that parameter needs the most careful experimental measurement.

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

Transcription and translation rates dominate (sensitivity 1.50). This means the most important experiment is characterizing the 6xNFAT promoter strength specifically in hair cells. Channel count and calcineurin affinity matter less: the system is robust to biological variation in these parameters.

Reproducible: full source code

The complete ODE model is available as a Python script. Dependencies: numpy, scipy. Run it yourself to reproduce these results or modify parameters.

View on GitHub: ode_model.py
← All hypotheses