我是一位没有科学背景的父亲，从事技术和创意制作工作。以下是我的每一个步骤，详细记录，以便任何面临相同处境的人都能够复现。
Michael的全外显子组测序（WES）报告来自香港儿童医院（实验室编号：23C7500174，2022年12月），列出了两个STRC变异。一个被标注为"致病性"（父系来源的全基因缺失，经MLPA确认）。另一个被标注为"意义不明变异（VUS）"（母系来源的单碱基改变，经Sanger测序确认）：NM_153700.2:c.4976A>C p.(Glu1659Ala)。我需要了解：这第二个变异是否真的有害？
我在UniProt上搜索"STRC"，找到静纤毛蛋白的编号为Q7RTU9。然后在AlphaFold中找到该蛋白的预测三维结构。1659位点的置信度评分（pLDDT）为95.69（满分100），说明该位点的结构预测非常可靠。
AlphaMissense是Google DeepMind开发的工具，用于预测蛋白质突变是否有害。我下载了静纤毛蛋白的预测文件，并搜索"E1659A"（E = 谷氨酸，原始氨基酸；A = 丙氨酸，Michael的变异）。
结果：0.9016/1.0（可能致病）。超过0.564即认为可能有害。我随后检查了1659位点所有其他19种可能的改变，每一种评分均在0.846以上。这意味着1659位点在结构上至关重要：任何改变都会破坏蛋白质功能。
| protein_variant | am_pathogenicity | am_class |
| E1659A | 0.9016 | LPath |
| E1659D | 0.9483 | LPath |
| E1659G | 0.9191 | LPath |
| ... 全部19种替换：LPath（0.846-0.999） | ||
如果某个位点对蛋白质功能重要，那么它在不同物种中应该是相同的氨基酸。我从UniProt下载了9个哺乳动物（人、小鼠、大鼠、牛、猴、猪、狗、蝙蝠、熊）的静纤毛蛋白序列，并在每个物种中搜索1659位点周围的基序。
结果：100%保守。所有9个物种在该位点均含有谷氨酸（E）。周围13个残基的基序（PEIFTEIGTIAAG）在约8000万年的进化历史中完全相同。这是ACMG标准下的PP1支持证据。
通常，遗传学家使用SIFT、PolyPhen-2和CADD来检验变异。我通过Ensembl VEP API尝试了三者，均对该变异无结果返回。
原因：STRC在15号染色体上有一个几乎完全相同的"孪生"基因（称为STRCP1的假基因），会干扰基于序列比对的工具。这正是AlphaMissense对STRC独特重要的原因：它基于蛋白质三维结构而非DNA序列，因此不受假基因影响。
ACMG/AMP指南（Richards et al., 2015）是遗传学家对变异进行分类的标准框架。每项证据均有代码和强度等级。我学习了规则并加以应用：
2项中等证据 + 2项支持证据 = 可能致病。按ACMG联合规则（表5），达到可能致病分类阈值。
我将所有证据整理成正式信函，呈送香港儿童医院化学病理学实验室，请求将该变异从VUS重分类为可能致病性变异。附件包括AlphaMissense数据、保守性分析和ACMG标准分解。我还建立了本网站，使证据透明、可复现，供审阅该病例的任何人查阅。
如果医院接受重分类，Michael的分子诊断将得到确认：双等位基因致病性STRC（DFNB16）。这是参与未来基因治疗临床试验的前提。双AAV基因治疗已在STRC缺陷小鼠中恢复了听力（Iranfar et al., 2026年1月）。人体临床试验预计将在2-3年内启动。届时Michael将年满7-8岁。
一旦我学会了如何阅读蛋白质结构，我意识到自己可以走得更远。我并非遗传学家，但这些工具是免费的，而我儿子的未来或许取决于有人提出正确的问题。于是我继续前行。
与其替换整个基因，能否只修复那一个错误的碱基？我从Ensembl下载了Michael变异周围的基因组序列，检验基因编辑工具能否将其靶向。
碱基编辑（CBE/ABE）：无法修复此变异（C>A颠换超出其编辑范围）。先导编辑：可行。我在距突变位点仅4个碱基对处找到了合适的PAM位点。先导编辑器理论上可以纠正该单碱基改变，但此方案尚未在内耳细胞中得到验证。
目前STRC基因治疗需要两个病毒载体（双AAV），因为基因太长无法装入一个。两个病毒意味着效率更低：两者必须同时进入同一细胞。我分析了AlphaFold结构，注意到前约600个氨基酸的结构置信度极低（pLDDT低于50），提示它们可能无法形成稳定结构，也许是可删除的。
若删除这些区域，剩余的"mini-静纤毛蛋白"（1076个氨基酸，3228 bp）可装入单个AAV载体。这是计算假说，需要实验室验证。但先例已经存在：微型肌营养不良蛋白（删除肌营养不良蛋白的非必需部分）目前已进入肌营养不良症治疗的III期临床试验。
为进一步验证mini-STRC假说，我向AlphaFold 3服务器提交了任务，预测静纤毛蛋白与其相互作用伴侣TMEM145结合的三维结构（TMEM145是2025年《自然通讯》最新发现的静纤毛蛋白功能必需蛋白）。
首批结果已获取（任务1）。ipTM = 0.47，pTM = 0.48。直接结合置信度较低。PAE矩阵分析显示最佳链间接触位于N端残基174-185（但PAE值仍较差，为8.6 A）。
随后我又提交了5个任务，系统性地检验mini-STRC假说：
| # | 实验 | 状态 | 测试内容 |
|---|---|---|---|
| 1 | 完整STRC + TMEM145 | 完成（ipTM 0.47） | 基线相互作用 |
| 2 | Mini-STRC + TMEM145 | 完成（ipTM 0.43） | N端可删除（0.43对0.47基线） |
| 3 | STRC E1659A突变体（单独） | 进行中 | Misha的突变会破坏折叠吗？ |
| 4 | STRC野生型（单独） | 完成（pTM 0.63） | 基线：16%无序（N端拉低分数） |
| 5 | Mini-STRC单独 | 完成（pTM 0.81） | 是！Mini-STRC折叠极佳（7%无序） |
| 6 | 仅N端（1-615） | 完成（pTM 0.27） | 已确认：38%无序，pTM 0.27 |
| 7 | mini-STRC + Piezo2 CED | Submitted | Does mini-STRC interact with mechanosensitive channel? |
| 8 | NFATC1 + Calcineurin A/B | Done (ipTM 0.73) | VALIDATED: CnA-CnB ipTM 0.91, NFAT-CnA ipTM 0.72. Cascade confirmed |
我向美国、法国和中国从事STRC基因治疗研究的顶尖科学家发送了电子邮件，分享了重分类证据、mini-STRC假说以及本网站的链接。
我收到了令人鼓舞的回复，确认计算方法思路正确，相关分析已转达给正在从事STRC基因治疗研究的团队。
On day three, we asked a question that changed the direction of the research: instead of a constitutive promoter that's always on, what if we used a promoter that responds to sound? Hair cells already convert sound to calcium signals. That's a built-in sensor we can hijack.
Hair cells have mechanotransduction (MET) channels that open when sound deflects their stereocilia. Ca²⁺ flows in, activating the calcineurin-NFAT signaling pathway. This pathway is well-characterized in the sonogenetics literature (Wu et al., Nature Communications 2023), where researchers used it to achieve 62-fold gene induction with zero background leakage.
We designed a construct: 6xNFAT promoter + mini-STRC. The 6xNFAT promoter is a synthetic element with 6 copies of the NFAT response element, creating a cooperative digital switch that requires sustained calcium signaling to activate. Combined with mini-STRC (from Day 2), the total construct is 3,893 bp, fitting within the 4,700 bp AAV packaging limit with 807 bp to spare.
To go beyond speculation, we built an ODE (ordinary differential equation) model of the complete signaling cascade: sound level → MET channel open probability → Ca²⁺ influx → calcineurin activation → NFAT nuclear translocation → transcription → translation → protein accumulation. Every parameter comes from peer-reviewed literature.
Result: with a realistic hearing aid schedule (16h on, 8h off), the model predicts therapeutic stereocilin levels (>15,000 molecules per OHC) in just 13 hours. In silence, the system produces only 6.8% of the target (promoter effectively OFF). The full Python code is available on GitHub for anyone to reproduce or extend.
To test the mechanosensitive hypothesis structurally, we submitted two new AlphaFold 3 jobs. Job 7: mini-STRC + Piezo2 CED (does stereocilin physically contact the mechanosensitive channel?). Job 8: NFATC1 + Calcineurin A/B (positive control, validates the cascade model).
Results pending. If Job 7 shows high ipTM (>0.6), it would mean stereocilin directly interacts with Piezo2, implying a feedback loop: broken STRC affects mechanosensation itself, not just structural connections.
Gene therapy currently requires surgical injection into the cochlea. What if ultrasound could push the treatment through the round window membrane non-invasively? I found papers showing that microbubbles + ultrasound create temporary pores in the RWM (Shih et al. 2019: 5.2x permeability increase, full recovery in 24 hours, zero hearing damage).
I built an ODE model of the full delivery chain: ultrasound → pore formation → nanoparticle diffusion → cell uptake → protein production. Honest result: baseline parameters fall short (0.39% per session, 258 sessions needed). Optimized LNPs with ionizable lipids: 78% per session, 2 sessions total. The bottleneck is endosomal escape, not the ultrasound. The most realistic path: one-time AAV surgery, plus non-invasive LNP top-ups years later if expression fades.
所有工具均免费。AlphaFold 3服务器需要Google账号。其他工具无需注册、无需API密钥、无需编程。总费用：0元。