6.**Result:** The synapse stops firing to save itself from excitotoxicity.
**CASCADE 1 — Vesicle Depletion** appears at `Max_RRP`, `Max_RP`, `p_release_base`, `k_rec_fast/slow`, `stochastic_release`, `map_trace_to_speed`, and the recruitment block in Loop 1. The key annotation explains the asymmetry: `p_release_base * Ca_micro` makes each spike draw *more* vesicles as Ca_micro rises early in a burst — a positive feedback that accelerates the collapse before recruitment can respond.
**CASCADE 2 — ATP Depletion** is anchored at `Glucose_level` (the root input) and at `compute_astrocyte_metabolic_health` in Loop 3, where it explains that `ATP_level` is the bridge variable that carries minute-scale metabolic state into the millisecond Ca²⁺ world.
**CASCADE 3 — Pump Failure** is annotated at `k_PMCA`, `k_NCX`, `k_SERCA`, `ATP_half`, and `compute_pump_atp_factor`. The NCX comment explicitly notes its role as a floor-not-rescue — it keeps clearing during failure and enables the auto-reset, but cannot prevent accumulation alone.
**CASCADE 4 — Residual Ca²⁺** appears at `B_total`, `tau_buffer_rebind`, the `capture_fraction` block, and the buffer recharge lines. The buffer saturation note explains the two-phase dynamic: buffer is protective early but becomes invisible once `B_free → 0`.
**CASCADE 5 — CDI Lock-out** is annotated at `k_CDI_rise`, `Ca_micro_saturation`, `k_CDI_rec`, and both the rise and recovery lines in Loop 1. The recovery comment spells out the self-locking logic explicitly as a chain.
**CASCADE 6 — Silence** sits at `effective_conductance` with a timing note showing that mGluR fires first, eCB second, and CDI last but irreversibly. The eCB and mGluR annotations in Loop 2 explain their roles as early partial brakes versus the terminal lock.
## 5. Model Summary Checklist
- \[x\] **Timing:** Spans 0.1 ms (AP) to 300,000 ms (Metabolism).
- \[x\] **Conservation:** Vesicles and Neurotransmitters are conserved through the Gln→RP→RRP→Cleft→Astro loop.
- \[x\] **Stability:** CDI and mGluR/eCB provide three layers of negative feedback to prevent runaway excitation.
- \[x\] **Resource Constraints:** ATP and Pool guards prevent physically impossible negative values or infinite accumulation.
---
This pseudocode serves as a comprehensive blueprint for a biologically consistent tripartite synapse. It is structured into three nested temporal loops to handle the transition from microseconds to minutes.
Gln_pool*=0.9# Account for metabolic overhead/loss
t+=dt# Increment simulation time
```
---
### 3. Biological Consistency Summary
1.**Metabolic Coupling:** The `atp_efficiency` variable creates a physical link between the 5-minute astrocyte clock and the 0.1ms calcium clock. If the astrocyte is exhausted, the pumps fail, and the `CDI_fac` locks the synapse into **silence**.
2.**Double-Lock Inhibition:** You have two distinct brakes. The `mGluR_pre` is "homosynaptic" (local self-check), while the `eCB_level` is "heterosynaptic" (feedback from the receiver). This prevents runaway excitation and excitotoxicity.
3.**Conservation of Mass:** Neurotransmitters follow a closed loop: `Gln_pool` -> `N_RP` -> `N_RRP` -> `Glu_cleft` -> `Astrocyte` -> `Gln_pool`.
4.**Buffer Buffer:** The `B_free` mechanism prevents numerical instability during high-frequency bursts by absorbing excess calcium influx instantly.
2.**Calcium buffer proteins** (calbindin, calmodulin) — fast buffering before pumps clear
3.**VGCC inactivation recovery** — CDI is mentioned but the recovery (de-inactivation) is missing
4.**RRP hard cap** — the refill logic can overshoot without a ceiling
5.**Reserve Pool (RP) depletion guard** — refill can go negative
6.**Presynaptic autoreceptors** (mGluR/CB1R feedback closing the loop from NT in cleft → presynaptic suppression) — separate from eCB which comes from postsynapse
7.**Spike refractory / AP waveform duration** — the effective window for Ca²⁺ entry is finite
Let me build an annotated code diagram and then write the enhanced model.Here is a full analysis of every missing loop, followed by the updated code for each section.
---
## Missing behaviors and the reasoning behind each
**Why Ca²⁺ clearance is the most critical gap.**`Ca_micro` currently has no exit route — it only grows. Without PMCA pumps, NCX exchangers, and SERCA (ER uptake), every spike leaves residual Ca²⁺ that accumulates across the simulation and eventually locks the synapse in a permanently over-activated (or CDI-locked) state. The clearance mechanisms also run at different speeds: NCX is fast (tens of ms), PMCA is slower but higher-capacity, SERCA is slowest and stores calcium for later use as an internal buffer.
**Why Ca²⁺ buffer proteins must precede clearance.** Calbindin and calmodulin bind free Ca²⁺ within microseconds and act as a fast, temporary "sponge". They blunt the initial `Ca_micro` peak, protecting against excess CDI. They also slowly release Ca²⁺ back into the cytosol, which feeds the trace integrator more smoothly. Without buffers, the microdomain pulse is unrealistically sharp.
**Why CDI recovery closes a loop without itself.** The model already writes `CDI_factor` but never resets it. A VGCC that inactivated on spike N stays inactivated on spike N+1. CDI recovery is simply a decay back toward zero, with a time constant of \~100 ms, driven by Ca²⁺ falling (i.e., it depends on clearance — another reason clearance comes first).
**Why mGluR autoreceptors are needed.** The eCB pathway is a *retrograde* signal synthesized by the *postsynapse*. But the presynapse also has its own direct cleft-sensing system: presynaptic mGluR2/3 receptors bind glutamate in the cleft and suppress VGCC conductance and cAMP. This is a homosynaptic feedback loop that is entirely local to the presynapse and missing from the current model.
**Why pool guards matter.** The `refill_amount` calculation can produce `N_RP < 0` if `current_recruitment_rate * N_RP > N_RP`. The `N_RRP` overshoot is subtler but also real: if two slow-loop updates stack before the fast loop consumes RRP, you can exceed `Max_RRP`.
---
## Updated code, section by section
### Loop 1A — Ca²⁺ dynamics (replaces the current `Ca_micro +=` block)
# This is the realistic metabolic-silence cascade.
defcompute_pump_atp_factor(ATP_level):
# Hill function: half-maximal at ATP_half = 0.3
returnATP_level**2/(ATP_level**2+ATP_half**2)
# The slowest metabolic loop already writes ATP_level via compute_astrocyte_metabolic_health().
# No further wiring needed — pump_scale above picks it up automatically.
```
---
## The closed loop, stated plainly
Every Ca²⁺ that enters now has exactly one exit: PMCA, NCX, or SERCA. Buffers slow the peak. CDI rises with Ca²⁺ and falls as Ca²⁺ falls — it can no longer lock permanently. The cleft-sensing mGluR autoreceptor gives the presynapse its own, faster brake independent of the postsynapse. Pool arithmetic is bounded on both ends. And ATP depletion now cascades naturally: less ATP → slower pumps → higher residual Ca²⁺ → more CDI → fewer effective VGCCs → silence — which is precisely the metabolic-fatigue endpoint the deep loop was trying to express but couldn't reach without the pump link.
