研究発表プレゼン要約 — 2026-04-24

• S2多くのシナプスモデルはケーブル理論ベースだが、スパインのような微小構造ではイオン濃度の動的変化を扱わないため精度が不足。電位勾配・electric drift を入れたモデルでも、シナプス電位の主要電荷担体である Na+ の動態を正確に再現できないことが多い。本研究は Qian が採用した PNP 方程式を改良し、実験由来の Na+ kinetics に基づく新モデルを開発。
• S3PNP ベース vs ケーブル理論ベースの違い (electric drift なしの EPSP)。
• S4計算の 3 ステップ: ① 細胞内イオン移動 → ② 膜貫通イオン移動 → ③ 膜電位変化。記号: σ=断面積, ρ=占有率 (occupancy), z=区画長, φ=電位勾配。
• S6参照: Arellano et al. (Yuste) 2007、Qian et al. (Sejnowski) 1989。

標準パラメータ (Slide 5)

区画分割: Head 6 / Neck 4 / PSD 1 / Dendrite 50 µm ×2。

大スパインでの Na+ kinetics (Slide 7)

主要な依存性

• S8spine head から dendrite に沿った Na+ 拡散。electric drift は拡散に影響しないが、Cl⁻ の移動には影響する。DNa 値を変えるべきか?
• S9拡散定数は dendrite 全体の Na+ 拡散にほとんど寄与しない。
• S10dendrite 直径が、spine head と親 dendrite の間の Na+ 振幅比を決める支配因子。
• S11electric drift は Na+ 振幅に大きく寄与しない (濃度勾配のみ vs PNP, AMPA only)。
• S12-13Mg²⁺ (1 mM vs 2 mM) の Na+ 電流 decay への影響は有意でない。シミュレーション時間はもっと長くすべき。NMDAR による Na+ influx 寄与は実験から 5.8% と算出。
• S14シナプス電流の大きさは Vm 振幅に影響するが Na+ decay time には影響しない。異なるパラメータで Na+ decay を測る際、Na+ 振幅変化を補償するためのシナプス電流調整は行わない。
• S15DNa が Na+ decay time の支配因子 (Miyazaki et al., 2017 と一致)。
• S16spine neck の DNa を変えても全体効果は無視できるが、関係は指数依存を示す。
• S17spine neck 直径は Na+ decay time に大きく影響しない。DNa が律速で、DNa < 0.2 で初めて大きな影響が出るが、これは spine neck 直径の範囲外。
• S18spine head volume は Na+ decay time に大きく寄与。最大 head volume (head 直径 0.6 µm, volume 0.212 µm³) を使用 — 「Na+ 動態は大きいスパインでのみ検出可能」という実験観察と整合。

• S19PNP と GHK を統合しているため、新たな予測的洞察を生成できる。
• S20新モデルでの GHK 式に基づく Vm 推定: Vm = 66.7 mV, rise 1.8, decay 5.5。spine neck resistance は非常に低い。
• S21新 PNP モデルはより低い spine neck resistance を予測。各イオンの更新式 (d=直径, l=長さ, ρ=体積占有率):

spine neck resistance は ΔVm / I_Na.synapse (= ΔVm / I_synapse) で計算。K+ は electric drift で spine head から出ていくため計算から除外。問題: I_K.synapse が必要 (K+ efflux を除外する処理)。

• S22Na+ 以外のイオン (Cl⁻, K+, Ca²⁺) の動態も予測。未実装/未最適化: VGCC が未実装、NMDAR ファクタ未最適化、Ca²⁺ ポンプパラメータ未最適化。
• S23electric drift の重要性: 濃度勾配のみ (electric drift 除去) にすると Cl⁻ と Vm の挙動が変わる。

• S24average spine での Na+ と Vm kinetics。
• S25より遅い DNa (3.3×10⁻⁶ cm²/s, Mondragão et al., 2016) を使い Miyazaki & Ross (2022) に調整。

• S26MICrONS プロジェクトの正確なデータを使用。Yuste の研究は TEM による平均 neck 径を報告するが、実際の最小径はもっと小さく、spine neck resistance に決定的な影響を与えうる。
• S27-28静止膜電位の影響: 過分極による Na+ 振幅の変化は Vm 振幅の変化より小さい。理由 — Vrest は Na+ 平衡電位から遠いため Na+ 振幅への影響は最小。一方 AMPAR の反転電位は Na+ 平衡電位より Vrest に近いため、変化が大きくなる。
• S29実験的検証: Na+ 振幅は Vrest と EPSP サイズに関して同じ傾向を示す。

現在の問題 (Slide 30)

• spine neck resistance の計算 — 最も影響するパラメータは何か?
• 古典的ケーブル理論の抵抗と異なり、現実の電流はイオンの物理的移動。光速よりはるかに遅いため電位変化に対し電流に伝播遅延が生じ、過渡状態ではオームの法則 (V=IR) が厳密には成り立たない。
• dendrite の Na+ spike 減衰の不一致: 実験では dendrite の Na+ 信号振幅が spine head の 50% 以上を保つが、シミュレーションでは 10–20% に低下。分岐を考慮した sim の方が理論上は現実的なはずなのに、実験は 40–50% 近い高振幅を示す。なぜか?
• bAP を模した全体的 Na+ 増加: 未完了。

検証テスト

• S31bAP の模倣: K+ の機構がないためイオン駆動力に深刻な問題。Tau = 750 ms, baseline 回復 2–4 s。
• S33細い spine neck は Na+ イオンの通過を妨げる。Vm は膜を横切る Na+ と K+ のバランスで決まる。
• S34Cm を 10 倍にすると EPSP の time constant が増加。
• S35MICrONS データからのスパインのランダムサンプリング。

Ca²⁺ / NMDAR の電荷収支計算 (Slide 37)

spine への推定 Ca²⁺ 流入 ≈ 2 µM。APV (NMDAR 拮抗薬) で Ca²⁺ 信号が 60% 減 → NMDAR 経由の Ca²⁺ 流入 = 1.2 µM (2022 年論文では 68% 減と報告)。この流入は全 Ca²⁺ 流入の 5% のみ (95% はバッファ, Sabatini 2002) → 総 Ca²⁺ 流入は約 24 µM。Schneggenburger et al. (1993) によると NMDAR を通る電流の 93% が Na+、7% が Ca²⁺。Ca²⁺ は 2 価なので 24 µM の電荷等価は 48 µM、これが 7% に相当。比 48:X = 7:93 から Na+ 電荷 X ≈ 638 µM。総 NMDAR 経由 Na+ 流入 = 0.638 mM (総電荷 = 686 µM の Na+ + Ca²⁺)。観測される自由 Na+ 変化は 12 mM 超 (Miyazaki et al., 2017) なので総 Na+ 流入はずっと大きく、総 Na+ 流入のうち NMDAR 由来は 5.8% 未満。