Dijitalmiyelin:Onarımıhesaplamak

Sorun: biyolojik sinyal var, klinik kanıt yok#

Remiyelinizasyon biyolojik olarak mümkündür. Oligodendrosit öncül hücreleri (OPC) lezyona toplanır, olgun oligodendrositlere farklılaşır, aksonun çevresini sarar ve yeni kılıfı sıkıştırır. Genç kemirgende bu program yaklaşık üç haftada tamamlanır; insanda da çalışır, ancak çoğu hastada farklılaşma aşamasında takılır.

Klinik denemelerin tablosu ise çarpıcı biçimde tek yönlü: klemastin yalnızca elektrofizyolojik bir sinyal üretti, anti-LINGO-1 (opicinumab) birincil sonlanımlarını ıskaladı, beksaroten MTR sonlanımını karşılamadı, klemastin progresif MS'te zarar sinyali verdi. 2025–2026 taramamızda hiçbir pro-remiyelinizan ajanın hakemli faz 2/3 birincil sonlanımını karşıladığı görülmedi.

Başarısızlığın nedenleri arasında sıkça tasarım değişkenleri anılıyor: hasta seçimi, lezyon yaşı, inflamasyon durumu, zamanlama. Bu soruları sistematik biçimde sınayacak, OPC-aracılı onarımı lezyon ölçeğinde VEP ve MTR benzeri ölçülebilir sinyallere bağlayan kalibre bir hesaplamalı model literatürde yok. Simülatör bu boşluğa yerleşiyor.

Dört katmanlı bir dijital miyelin#

Simülatör dört katmandan oluşur ve 200 aksonu 5 tohumda dakikalar içinde tarayabilir. Kod saf numpy / scipy'dir.

A · Yapı ve iletim#

Her aksonun çapı log-normal dağılır (medyan 1,2 µm, σ = 0,35); internod uzunluğu $L_0 = 100\,d$, sağlıklı g-oranı 0,70. İletim hızı:

$$v_i = k \cdot c_{\text{MSS}} \cdot d_i \cdot f(g_i) \cdot h\!\left(\tfrac{L_i}{L_0}\right)$$

Sabitler $k = 6\ \text{m}/(\text{s}\cdot\mu\text{m})$ ve $c_{\text{MSS}} = 0{,}8$'dir. Burada $f(g)$ Rushton profilidir: $f(g) = \sqrt{-\ln g}\,/\,\sqrt{-\ln 0{,}70}$, yani $f(0{,}80) = 0{,}79$ ve $f(0{,}95) = 0{,}38$. İletilen spike oranı dik bir sigmoiddir:

$$s(g) = \frac{1}{1 + e^{(g - 0{,}93)/0{,}012}}$$

Sigmoid o kadar diktir ki $s(0{,}90) = 0{,}92$ iken $s(0{,}95) = 0{,}16$'ya düşer: kılıf biraz daha incelince iletim neredeyse tamamen bloke olur.

Sağlıklı ağda gecikmeler gelişimsel uyarlama ilkesiyle eş zamanlıdır ($\tau_i \approx 6$ ms ± %6). Remiyelinize kılıf $g \to 0{,}80$ platosuna oturur ve internodlar kısalır; lif hızı sağlıklının yaklaşık %70–77'sine döner. Ana sayfadaki canlı akson tam olarak bu iki fonksiyonla çalışıyor: imleçle kılıfı yırttığınızda iletim hızı $f(g)$ ile düşer, blok olasılığı $1 - s(g)$ ile artar.

B · İşlev: 40 Hz'lik gecikmeli bir Kuramoto ağı#

200 kaynak osilatör 40 Hz ortak ritimde faz kilitlidir; her akson 8×5 topografik ızgarada bir hedefe projekte eder:

$$\begin{aligned} \dot\theta_j = {} & \omega_j + \frac{K_{ff}}{n_j}\sum_{i \to j} s_i \sin\!\left(\omega t + \varphi_i - \omega\tau_i - \theta_j\right) \\ & + K_{loc}\,|z|\,\sin(\arg z - \theta_j) + \sqrt{2D}\,\xi \end{aligned}$$

İşlev ölçütü hedef faz senkronisidir: $R = \langle |z| \rangle$, sağlıklı ağda $R = 0{,}85 \pm 0{,}03$. Aşağıdaki deneyde hasarlı aksonların biriktirdiği gecikmenin hedefleri nasıl fazdan düşürdüğünü görebilirsiniz.

L-0240 Hz Senkroni Halkası

Tam ekran ↗

Canlı deney · kaydırınca başlar

C · Hücre dinamiği: onarımın biyolojisi#

Her lezyonlu akson için enkaz $D$, inflamasyon $I(t)$, OPC yoğunluğu $P$, oligodendrosit örtüsü $O$, sarma $W$, kompaksiyon $C$ ve alıcılık $\rho$ izlenir. Çekirdek denklemler:

$$\begin{aligned} \dot O &= k_{\text{diff}}\,A^{0{,}7}\,B(I)\,(1 - D)\,\rho\,P\,(1 - O) \\ \dot W &= k_{\text{wrap}}\,\frac{\max(O - W, 0)}{d/d_{\text{ref}}} \\ \dot C &= \frac{W - C}{\tau_{\text{comp}}} \\ \dot \rho &= -\frac{(1 - O)\,\rho}{\tau_{\text{window}}} \end{aligned}$$

Farklılaşma bloğu $B(I)$ v1'de Hill biçimindedir, $B(I) = 1/(1 + (I/0{,}5)^6)$; v2 buna çarpımsal ($1 - \beta I$) ve kapı tipi (eşik altı 1, üstü 0) mekanizmaları ekler. Yaş skaleri $A$, OPC hızlarını $A^{0{,}7}$ ile, enkaz temizliğini $A$ ile ölçekler. Zaman sabitleri lizolesitin / kuprizon literatüründen gelir: toplanma 3–5 gün, farklılaşma 7–14 gün, sarma 14–21 gün.

L-01Dijital Miyelin: Katman C

Tam ekran ↗

Canlı deney · kaydırınca başlar

D · Gradyan tabanlı onarım: yapay zekâ ne yapardı?#

Biyolojik onarımla aynı gözlenebilirleri optimize eden “işlev-sürümlü” bir tamirci kurduk. Parametre, lezyonlu aksonların g-oranıdır; kayıp:

$$J(g) = \bigl(1 - R(g)\bigr) + \lambda\,\frac{\text{latans}(g)}{\text{latans}_0} + \mu\,\text{hacim}(g)$$

Ağırlıklar $\lambda = 0{,}2$ ve $\mu = 0{,}02$'dir. Optimizerler elle yazılmış SGD, momentum, Nesterov, Adam ve Adam + gürültüdür. v2'de türevlenebilir bir vekil yerine gerçek Kuramoto $R$'si üzerinde ortak rastgele sayılarla sonlu fark gradyanı kullanıldı; vekilin nicel sonuçları yanlış yönlendirdiği görüldü. Bu tamircinin biyolojiden ayrıştığı her nokta, biyolojinin nerede “işlev-kör” olduğunu gösteren bir ipucu.

Bulgular#

Latans hacimden iki kat hızlı düzeliyor#

Odaksal lezyonda (%20 akson, $g \to 0{,}95$) senkroni $R$ sağlıklı 0,852'den 0,683'e düştü ve 120. günde 0,812'ye döndü. VEP-latans vekili 5,99 → 8,31 → 6,53 ms izledi. Latansın T90'ı 22 gün, miyelin hacminin T90'ı 47 gündü. Yorum: ince bir ilk sarma iletimi geri getirir, kompaksiyon MTR sinyalini ancak geç yükseltir. Kalıcı latans artığı ($6{,}53$'e karşı $5{,}99$ ms) ince kılıftan ($g \approx 0{,}80$) ve kısa internodlardan ($L \approx 0{,}54\,L_0$) kaynaklanır. Bu, ölçekten bağımsız bir test öngörüsüne dönüşür: VEP-T90 / MTR-T90 ≈ 0,5.

Yaş süreyi, kronik inflamasyon kapsamı belirliyor#

Düşük ya da geçici inflamasyonda yaş, onarım süresini 1,3–1,9 kat uzattı ama kapsamı değiştirmedi (1,00). Kronik, PRL benzeri profilde ise kapsam yaşla birlikte 0,57 → 0,42 → 0,32'ye indi; lezyon 40'tan 57 aksona genişledi.

Soğuk başlangıç takılır, “kör” ince bir ilk sarma vadiyi aşar#

Blokta başlayan gradyan onarımı ($g_0 = 0{,}95$) yerel bir minimumda takıldı: 40 Hz'de $g = 0{,}95$ iken gecikme 15,7 ms'dir ve bu, sağlıklıya göre ≈140° faz farkı demektir. Bloğu kaldırmak anti-faz bir girdi ekler, $R$ düşer ve gradyan $g$'yi yeniden bloğa iter. İşleve bakmadan yapılan ince bir ilk sarma ($g_0 = 0{,}90$) ise vadiyi aştı. Gerçek-R gradyanıyla bu nitel sonuç korundu (9,4 ± 6,0'a karşı 32,6 ± 7,9 / 40); v1'deki “40/40 tam onarım” iddiası ise vekile özgü bir artefakt olduğu için düşürüldü.

L-03Soğuk Başlangıç, Ilık Başlangıç

Tam ekran ↗

Canlı deney · kaydırınca başlar

Biyoloji ise onarımı işleve göre değil çapa göre sıralar ($\rho_d = 0{,}994$): ince aksonlar önce onarılır. Gradyan sırası biyolojik sırayla ilişkisizdir (Spearman −0,02…+0,03) ve gecikme hatasını izler.

Anti-faz vadisi yalnızca gama bandında#

Onarım yolunda (g 0,95 → 0,80) senkroninin geçici olarak çukura düşmesi taşıyıcı frekansla belirlenir: gecikme hatası (≈9,7 ms) çeyrek periyodu aştığında vadi oluşur.

Model öngörüsü: remiyelinizasyonun ilk haftalarında gama koheransında geçici bir çukur, alfa/beta koheransında ise korunma. Bu, MEG zaman serileriyle doğrudan sınanabilir.

Blok mekanizması ilaç yanıtını belirliyor#

Aynı farklılaşma hızlandırıcısının ($k_{\text{diff}} \times 2$, klemastin-benzeri) kronik lezyondaki etkisi tamamen blok mekanizmasına bağlı çıktı. $I = 0{,}8$'de kapsam Hill blokta 0,39 → 0,63, çarpımsal blokta 0,83 → 0,97 oldu; kapı tipi blokta ise eşik üstünde sıfır kazanç verdi. Bu, faz 2 denemelerindeki yanıt heterojenliğine bir açıklama önerir: hastaların inflamasyon düzeyleri kademeli bir blok eğrisi üzerinde dağılıyorsa kısmi yanıt, kapı modeli geçerliyse iki-modlu bir “tam yanıt / sıfır yanıt” dağılımı beklenir.

Zamanlama: kazanç lezyon yaşıyla sönüyor#

Hızlandırıcı lezyondan $t$ gün sonra başlatıldığında ($\tau_{\text{window}} = 45$ g) kapsam kazancı +0,26'dan (t = 0) +0,06'ya (t = 90 g) tekdüze azaldı; yarılanma ≈45–50 gün. Kazanç yaklaşık $\rho(t) = e^{-t/\tau_{\text{window}}}$ ile aynı ölçekte sönüyor. Tüm kinetik bir insan ölçeği çarpanıyla ölçeklendiğinde kapsam değişmiyor; kapsamı belirleyen boyutsuz grup $\Pi_{\text{window}} = k_{\text{diff}} \cdot \tau_{\text{window}} = 5{,}4$.

Aktivite hızlandırır ama sırayı değiştirmez#

Aktivite kazancı ($\kappa = 3$) toparlanmayı %16 hızlandırdı (T90 22,0 → 18,4 g), ancak onarım sırasını çap belirlemeye devam etti ($\rho_d$ 0,994 → 0,84). Türetilen hipotez: aktivite sırayı değiştirecekse, sarmanın kalınlığına ve internod uzunluğuna etki etmelidir.

Test edilebilir hipotezler#

En güçlü itirazlar#

“Kemirgen zaman sabitleri insana taşınmaz.” Doğru. Bu yüzden v2'de tüm kinetik boyutsuz gruplara indirgendi ($\Pi_{\text{window}} = 5{,}4$, $\Pi_{\text{infl}} = 0{,}36$). Aktarım sorunu “ölçek” değil, bu iki grubun insan değerleridir; longitudinal PET/MTR ile kestirilebilirler.

“Ilık başlangıç sonucu vekil-R artefaktı.” Kısmen. v1'in 40/40'ı vekile özgüydü; gerçek Kuramoto $R$'siyle 32,6 ± 7,9 / 40 (tohumlar arası 21–40) ve yerel minimum 9,4 / 40 doğrulandı. Nitel sonuç korunur; nicel “tam onarım” iddiası düşürüldü. Gerçek-R gradyanının $R \approx 0{,}82$ doyumunda rastgele yürüyüşe dönmesi ayrı bir bulgu olarak raporlandı.

“Bu bir kanıt değil.” Kesinlikle. Simülatörün en verimli çıktısı, gradyan katmanının biyolojik katmandan ayrıştığı noktalardır; her ayrışma bir mekanizma sorusuna ve bir test tasarımına dönüştürülür.

Sırada ne var#

Simülatör v3 için hedefler: kılıf büyümesine bağlı bir aktivite kuralı, longitudinal PET/MTR verisiyle insan $\Pi_{\text{window}}$ ve $\Pi_{\text{infl}}$ kalibrasyonu, eş-şiddetli blok mekanizması eşlemesi, akson kaybı durum değişkeni ve MEG ile karşılaştırılabilir alfa/gama koherans çıktısı.

Çalışmanın ikinci yarısı, aynı gradyan yöntemlerini bir tanı problemine uyguluyor: BOS gerektirmeyen MS tanısını tek bir MR oturumundan hesaplamak.

Tam makale 34 sayfa, 135 doğrulanmış DOI içeriyor: PDF'i indirin. Alan uzmanları için iki sayfalık özet de mevcut.