- Vesicle recycling with time constant τ\_recycle (seconds)
### Complete Model
#### 1. Short-Term Facilitation (STF)
- **Mechanism:** Activity-dependent **increase in release probability (P<sub>r</sub>)** due to accumulation of `residual_Ca²⁺` in the active zone.
- **Model Implementation:**
-`residual_Ca²⁺` is a state variable that acts as a short-term memory buffer.
- Each presynaptic spike injects a fixed amount of Ca²⁺ into `[Ca²⁺]_i`, which then decays exponentially with a fast time constant (τ\_Ca_fast \~ 10-50 ms). A fraction of this is added to `residual_Ca²⁺`, which decays with a slower time constant (τ\_Ca_slow \~ 100-500 ms).
- Where `F_max` is the maximum facilitation strength and `K_d` is the half-saturation constant.
- **Effect:** High-frequency spike trains cause `residual_Ca²⁺` to summate, progressively increasing the `STF_factor` and therefore `P_r` for subsequent spikes.
#### 2. Short-Term Depression (STD)
- **Mechanism:** Activity-dependent **decrease in release** due to depletion of the readily `releasable vesicle_pool`.
- **Model Implementation:**
- Each successful release event (a stochastic outcome with probability `P_r`) reduces the `vesicle_pool` by one quantum.
- Vesicles are replenished from a reserve pool with a **recovery time constant τ\_recycle** (e.g., 0.5 - 10 seconds).
- The **homeostatic_factor** (or depletion factor) is simply the fraction of the pool that is available: `vesicle_pool / vesicle_pool_max`.
- **Effect:** High-frequency firing depletes the pool faster than τ\_recycle can refill it. The `homeostatic_factor` drops, decreasing the *effective* release rate (`P_r * vesicle_pool`), even if `P_r` itself is high.
#### 3. Long-Term Modulation (LTP/LTD of Release)
- **Mechanism:** Retrograde chemical signals (**eCB, NO, BDNF**) induce biochemical cascades that **modify the baseline**`P_r` on timescales of minutes to hours.
- **Model Implementation:**
- These are triggered by specific postsynaptic activity patterns (e.g., postsynaptic Ca²⁺ spikes for eCB, strong depolarization for NO).
- They act as **scaling factors on the baseline**`P_r`:
-`rCa(t)` is the `residual_Ca²⁺` buffer, driven by `[Ca²⁺]_i` spikes.
-`LTD/LTP_factor = α_ECB * β_NO * ...`**(Quasi-static; changes on long timescales)**
-`homeostatic_factor(t) = vesicle_pool(t) / vesicle_pool_max`**(Dynamic; updates per release event)**
**Final Release Decision:** A random number is drawn. If it is `< P_r(t)`, a vesicle fuses, `vesicle_pool` is decremented, and glutamate is released. The `vesicle_pool` then recovers toward its maximum with the time constant `τ_recycle`.
This model captures the core **tension between facilitation (driven by Ca²⁺) and depression (driven by depletion)**, while allowing for slower, homeostatic and Hebbian adjustments, making it a powerful framework for simulating synaptic dynamics.
- **Scaling:** Global adjustment of all `N_AMPA` based on soma firing rate
### Complete Model
#### 1. Voltage-Dependent Gate (NMDA Receptor)
- **Mechanism:** The `NMDA_conductance` (`g_NMDA`) is not constant. It is *gated* by both glutamate binding and the relief of a voltage-dependent Mg²⁺ block.
- **Model Implementation (Instantaneous):**
-`g_NMDA(t) = N_NMDA * γ * B(V_m(t))`
- Where `γ` is the single-channel conductance, and `B(V_m)` is the magnesium unblock fraction.
- **Effect:** This function ensures NMDA receptors are only significant **coincidence detectors**. They pass current only when presynaptic glutamate release (`N_NMDA` is bound) AND postsynaptic depolarization (from `V_m` or bAP) occurs simultaneously.
- **Core Signal:** The key trigger for all plasticity is the postsynaptic **calcium transient**`[Ca²⁺]_i(t)`, which integrates multiple sources:
-`[Ca²⁺]_i(t) = J_NMDA(t) + J_VGCC(t) + J_IP3(t)`
- Where `J_NMDA` is Ca²⁺ influx through NMDA receptors, `J_VGCC` is from voltage-gated channels opened by bAPs, and `J_IP3` is from mGluR/second-messenger pathways.
- A `plasticity_tag` is set to `1` if `[Ca²⁺]_i` crosses a moderate threshold (`θ_tag`) within a short time window (\~100ms). This marks the spine as "eligible."
- The final change in `N_AMPA` is then determined by the peak/amplitude of the calcium signal:
- **LTP:** If peak `[Ca²⁺]_i > θ_LTP` (a high threshold), then `ΔN_AMPA = +A_LTP * plasticity_tag`. This typically requires **strong, coincident** presynaptic glutamate AND a bAP.
- **LTD:** If `θ_LTD < peak [Ca²⁺]_i < θ_LTP` (a moderate, sustained level), then `ΔN_AMPA = -A_LTD`. This can be triggered by presynaptic activity alone or weak pairing.
- The **kinase_activity** and **phosphatase_activity** in your rule are functions of `[Ca²⁺]_i` (e.g., `kinase ∝ ([Ca²⁺]_i - θ_LTP)^2` for LTP).
#### 3. Retrograde Signal Synthesis
- **Mechanism:** The postsynaptic spine functions as a **signal interpreter and broadcaster**. Based on the calcium signal, it synthesizes specific retrograde messengers.
- **Model Implementation (Threshold Logic):**
- **eCB synthesis:** Triggered if `[Ca²⁺]_i > θ_eCB`. eCB is then released, diffusing back to inhibit the presynaptic terminal (lowering `baseline P_r` via your presynaptic `LTD_factor`). This is often a form of heterosynaptic LTD.
- **NO synthesis:** Triggered by a similar high `[Ca²⁺]_i` threshold coupled with activation of specific enzymes (e.g., nNOS). NO diffuses to the presynapse to increase `P_r` (your `LTP_factor`).
- **BDNF synthesis:** Slower, triggered by sustained calcium signals or specific gene activation pathways. BDNF acts both presynaptically and postsynaptically to promote structural changes.
#### 4. Structural & Metaplastic Modulation
- **Metabotropic (mGluR) Pathway:** Activated by sustained or spillover glutamate. It doesn't directly cause plasticity but **modulates the plasticity thresholds** (`θ_LTP`, `θ_LTD`). For example, mGluR activation can lower `θ_LTD`, making LTD easier to induce (metaplasticity).
- **Spine Volume (**`spine_volume`**):** This is a slow variable that couples to receptor counts.
- **Growth:** A sustained increase in `N_AMPA` (from LTP) promotes actin polymerization, increasing `spine_volume`.
- **Shrinkage/Stabilization:** Conversely, a large, stable `spine_volume` creates slots for more `N_AMPA`, stabilizing the potentiation. This creates a positive feedback loop for strong, stable synapses.
#### 5. Homeostatic Scaling
- **Mechanism:** A global, cell-wide feedback mechanism to maintain the soma's average firing rate within a target range.
- **Model Implementation (Multiplicative):**
- Periodically (e.g., every few hours of simulation time), the soma calculates its average firing rate `r_avg`.
- If `r_avg` deviates from a target `r_target`, all `N_AMPA` on all synapses are scaled uniformly:
- If `r_avg < r_target`: `N_AMPA = N_AMPA * β_up` (where `β_up > 1`). This is **up-scaling**.
- If `r_avg > r_target`: `N_AMPA = N_AMPA * β_down` (where `β_down < 1`). This is **down-scaling**.
- This rule is applied **independently** of the Hebbian `plasticity_tag`. It ensures network stability.
Where `B_AMPA(t)` is the fraction of AMPARs bound by glutamate (a transient pulse upon release).
**Simulation Cycle for a Spine:**
1.**Input:** Glutamate binds. bAP may arrive.
2.**Integration:**`V_m` depolarizes locally. Mg²⁺ block is relieved → `g_NMDA(t)` computed.
3.**Calcium:**`[Ca²⁺]_i(t)` is calculated from all sources.
4.**Decision:**
- Set `plasticity_tag`.
- Compute instantaneous `ΔN_AMPA` via the LTP/LTD calcium rule.
- Trigger retrograde signal synthesis if thresholds crossed.
5.**Update:** Apply `ΔN_AMPA`. Slowly update `spine_volume`. Periodically apply global scaling.
This model captures the spine as a **coincidence detector, integrator, biochemical decoder, and structural adaptor**—a fundamental unit of learning and memory.
---
---
---
### **3. Dendrite (Branch - Pattern Detector)**
**Internal State Variables:**
-`V_m(z)`: Space-dependent membrane potential
-`[Ca²⁺]_i(z)`: Local Ca²⁺ concentration
-`NaV_density(z)`: Sodium channel distribution
-`VGCC_density(z)`: Calcium channel distribution
-`branch_excitability`: Global gain factor
**Incoming Signals:**
- **Local:** EPSPs from spines (summed spatially/temporally)
#### 2. Modulation of Branch Excitability & Pattern Separation
Branch-specific modulators (**dopamine, acetylcholine**) reconfigure the branch's **global gain and plasticity thresholds**:
- **Dopamine (D1 receptor):** Increases `branch_excitability` (↑ `NaV_density` sensitivity) and **lowers θ\_Na**. It also gates plasticity: `plasticity_enabled = dopamine_present and coincidence_detected`. This primes the branch for learning salient, reward-predicting patterns.
- **Acetylcholine (muscarinic):** Enhances `VGCC_density(z)` efficacy and NMDA conductance. It promotes **Ca²⁺ spike generation** over Na⁺ spikes, favoring slower, integrative pattern detection over fast propagation.
- **GABA Inhibition:** This is crucial for pattern separation. A GABAergic input onto the branch shunts `V_local(t)`:
- By selectively inhibiting specific branches, interneurons **prevent those branches from reaching spike threshold**, ensuring only the most strongly activated, distinct input combinations generate output. This forces different branches to learn and respond to different patterns.
#### 3. Spatial Computation & Weight Updates
The **pattern separation** emerges from the interaction of localized synaptic inputs and branch-wide thresholds:
1.**Input Pattern:** A set of active spines delivers EPSPs to `V_local(t)`.
2.**Branch Filter:** The combination of `NaV_density(z)`, `VGCC_density(z)`, `branch_excitability`, and local inhibition determines a **unique activation threshold** for that branch.
3.**Pattern Detection:** Only input combinations whose summed `V_local(t)` exceeds this threshold generate a dendritic spike. Slightly different patterns may fail, especially with GABAergic tuning.
4.**Synaptic Tagging (Credit Assignment):** When a dendritic spike coincides with a bAP, it generates a `global_teaching_signal`. This signal is broadcast **retrogradely** but **locally**, tagging all *recently active* spines within a spatial radius `R`. The tag's strength decays with distance from the spike initiation zone.
-`NaV_density(z) += η * (target_activity - dendritic_spike_rate)`
- This allows branches to self-tune their excitability over long timescales.
-`branch_excitability`**:** A slow variable modulated by neuromodulators (↑ by DA) and metaplasticity rules.
#### Summary: The Branch as a Feature Detector
In this model, a dendritic branch is not a passive cable. It is an **active feature detector** with tunable properties:
- **Input:** A vector of synaptic inputs (spatially arranged).
- **Nonlinearity:** A double-threshold operation (Na⁺/Ca²⁺ spike generation) determined by its ion channel makeup and modulatory state.
- **Output:** Either a large, propagating dendritic spike (a **binary feature detection event**) or graded subthreshold voltage.
- **Learning:** Synapses on the branch are updated based on a **three-factor rule**:
1. Presynaptic activity (glutamate release).
2. Postsynaptic dendritic spike (local).
3. Global teaching signal (bAP coincidence, modulated by DA).
- **Function:** Different branches, through their unique channel densities and inhibition, become selective for **different combinations of inputs**, implementing a powerful form of **dendritic pattern separation** that vastly expands the computational capacity of a single neuron.
- Each spike adds to somatic `[Ca²⁺]_i`, which slowly increases `g_sAHP`. This hyperpolarizes the cell, making it harder to reach threshold for subsequent spikes. **Prevents runaway excitation and encodes temporal derivatives.**
- This is the **global synaptic scaling** command sent to all synapses (as referenced in your postsynaptic model). It multiplicatively adjusts `N_AMPA` everywhere, a slow, cell-wide negative feedback loop.
#### 3. Soma as Transcriptional & Metabolic Hub
- **Somatic** `[Ca²⁺]_i`**Integration:** Sustained high `firing_rate_avg` leads to sustained elevated somatic `[Ca²⁺]_i`. This activates transcription factors (e.g., CREB).
- **More ion channels** (long-term excitability changes).
- **Neurotrophic factors** (e.g., BDNF) released to further modify network.
- **Metabolic Gatekeeping:** The `lactate` and `oxygen` signals directly influence ATP production. Low energy → upregulate `I_h` and `g_K_leak` to **reduce metabolic cost** by lowering firing rate—a direct link from metabolism to excitability policy.
#### Summary: The Soma as Central Processor
In this model, the soma is not a simple point neuron. It is a **dynamic policy engine** that:
1.**Integrates** spatially and temporally filtered inputs from dendritic subunits.
2.**Generates** all-or-none output decisions (APs) based on a modifiable threshold.
3.**Broadcasts** teaching signals (bAPs) back to the dendritic computational layers.
4.**Adapts** its own sensitivity on short (sAHP) and long (channel expression) timescales.
5.**Reconfigures** its entire input-output function based on neuromodulatory state (gain, threshold, integration window).
6.**Orchestrates** whole-cell homeostasis via global scaling commands and transcriptional programs.
This transforms the classic "integrate-and-fire" unit into a **biological central processing unit (CPU) with dynamic clock speed, adjustable gain, and multiple feedback control systems**, all dedicated to maintaining stability while allowing for state-dependent, plastic computation.
# Makes threshold higher if cell is too active, lower if too quiet (homeostatic)
```
#### 2. Dynamic Properties: Specific Mechanisms
**A. Threshold Plasticity & Modulation**
This is a key regulatory point. The `threshold_baseline` is not fixed; it's a **homeostatically regulated variable** and a **target for neuromodulation**.
- **Activity-Dependent (Homeostatic):** As shown above, sustained high `recent_spike_rate` increases `threshold_baseline`, making the neuron harder to fire (negative feedback).
- **Phosphorylation-Dependent (Modulatory):** Kinases activated by neuromodulators (PKA, PKC, CK2) phosphorylate specific sites on Nav channels (e.g., Naᵥ1.6).
- **PKA Phosphorylation (e.g., via DA/NE):** Shifts `V_half_inact` to more **depolarized** voltages → increases `h_inf` at resting V_m → effectively **increases NaV_availability** and **lowers effective threshold**. **Policy:***Lower threshold, increase excitability.*
- **CK2 Phosphorylation:** Can shift activation `V_half_act` to more **hyperpolarized** voltages → channels open easier → **lowers threshold**. **Policy:***Increase temporal precision and reliability.*
**B. Reliability (High Safety Factor)**
This is modeled implicitly by the steepness of the Nav activation curve and the high channel density.
- The condition `(V_m_AIS > threshold_dynamic)` is not a linear probability. It's a **step function** because the activation variable `m` of Nav channels is a steep sigmoid:
- With a high density of channels, once `V_m_AIS` crosses `threshold_dynamic` (where `m_inf` becomes significant), the positive feedback of Na⁺ influx is explosive and deterministic. There is no stochastic "maybe" spike.
**C. Timing Precision (Submillisecond Jitter)**
Jitter is minimized by three model features:
1.**Rapid Kinetics:** Very small `τ_m` (activation time constant) for AIS Nav channels (\~0.1 ms).
2.**High dV/dt:** The somatic `V_m` must rise rapidly to cross the AIS threshold. Slow ramps will not trigger a precise spike. This is enforced by the requirement for a strong, synchronous dendritic input to create a fast somatic depolarization.
3.**Refractory State Clarity:** The absolute refractory period (`τ_abs_refractory`) is a hard lockout. The relative refractory period is modeled by the recovery of `NaV_availability` and the elevated `threshold_dynamic` post-spike, which together sharply define the earliest possible next spike time.
#### 3. Role in Backpropagation (bAP) Initiation
The AIS is the **source** of the backpropagating action potential.
- Upon AIS spike generation, the depolarizing current not only propagates down the axon but also **actively back-invades** the soma and dendrites.
-`bAP_strength` in the model can be modulated (e.g., increased by dopamine signaling), affecting the amplitude of this critical teaching signal throughout the dendritic tree.
#### Summary: The AIS as a Programmable Binary Converter
In this model, the Axon Initial Segment is the **final, decisive policy layer**:
- **Processing:** A dynamic threshold function, gated by channel availability and phosphorylation state.
- **Output:** A stereotyped action potential (or not) with high temporal fidelity.
- **Key Modulation:** Its **excitability is tunable** via:
- **Homeostatic Threshold Plasticity:** Keeps average firing rate in check.
- **Phosphorylation States:** Allow neuromodulators (DA, NE, ACh) to directly adjust the "trigger happiness" of the neuron on fast timescales.
- **Refractory Kinetics:** Control maximum firing frequency and temporal precision.
This transforms the AIS from a passive fuse into an **active, tunable decision node** that finalizes the neuron's output based on integrated somatic potential, while itself being subject to meta-level policy controls that set the neuron's overall responsiveness and reliability.
- This effectively **gates synaptic plasticity**: only synapses under active astrocytic "supervision" (high Ca²⁺ in astrocyte) have fully functional NMDARs and can undergo LTP.
- Adenosine acts on neuronal **A1 receptors**, universally inhibiting synaptic release (presynaptic) and excitability.
- **This is a global negative feedback loop:** High network activity → more astrocytic glutamate uptake → more adenosine → stronger network-wide inhibition → enforced **activity quota** leading to sleep.
**E. Waste Manager (Glymphatic Coordination)**
- **State-Dependent Clearance:**
- During sleep/wake cycle, astrocyte **aquaporin-4 (AQP4)** polarization changes.
- Astrocytes dynamically regulate **perivascular space** volume to facilitate convective flow of cerebrospinal fluid during sleep, clearing metabolites like Aβ.
#### 3. Modulation of the Astrocyte
The astrocyte itself is modulated by:
- **Noradrenaline (from locus coeruleus):** ↑ `IP3 production` → ↑ `[Ca²⁺]_i` waves → potentiates gliotransmitter release. **Policy:***Alertness mode*– enhances astrocytic support for heightened neural activity.
- **ATP/Purinergic Signaling:** From active neurons or other astrocytes. Can propagate **Ca²⁺ waves** across the astrocyte syncytium, enabling long-range coordination of homeostatic states.
#### Summary: The Astrocyte as a Multi-Scale Regulator
This model presents the astrocyte as a **biological real-time operating system (RTOS)** for the brain microenvironment:
It introduces **critical non-neuronal constraints** into the neural network model:
- **Energy is limited and dynamically allocated.**
- **Waste products accumulate and must be cleared.**
- **Ionic balance must be maintained to prevent instability (seizures).**
- **Synaptic plasticity is chemically gated by glial oversight.**
In the full system model, the astrocyte is the **homeostatic backdrop** against which the neuron's computational drama plays out—a dynamic environment that both supports and constrains neural activity, ultimately ensuring the system's long-term stability and health.
This model architecture creates a **recursive optimization system** where each component's behavior adjusts based on both immediate inputs and long-term trends, with astrocytes providing the essential metabolic and environmental context that makes sustained neural computation possible.
