Decoding the Inflammatory Reflex: A Computational Primer on the VNS CAP Neuroimmune Axis
- Apr 9
- 4 min read
Author & Technical Reporter: SAEHL research
Co-Author: AI
Date Published: 9th April 2026
For decades, immunology and neuroscience were largely treated as non-overlapping magisteria. The central nervous system (CNS) was "immune privileged," and the immune system was an autonomous, self-regulating army. That paradigm fractured with the discovery of the Inflammatory Reflex—a hardwired neural circuit where the brain acts as the ultimate immunomodulator (Tracey, 2002).
Today, we are looking under the hood of a computational "toy model" (the VNS-CAP Neuroimmune Simulation Design) to understand how we can mathematically simulate this reflex.
By translating neurophysiology into code, we can explore the biophysics of bioelectronic medicine e.g. using pacemakers for the vagus nerve to treat systemic inflammation.
The Biological Architecture
To understand the simulation's logic, we must first map the biological entities to computational variables. The Cholinergic Anti-inflammatory Pathway (CAP) operates as a classic closed-loop feedback system:
The Sensor (Afferent Arc): Pathogen-associated molecular patterns (PAMPs), like endotoxin (LPS), trigger macrophages to produce pro-inflammatory cytokines (TNF-𝛼, IL-1𝛽). Afferent vagus nerve fibers express cytokine receptors, detecting this inflammatory spike and transmitting the signal to the nucleus tractus solitarius (NTS) in the brainstem.
The Controller: The brainstem integrates this signal and fires an efferent response via the dorsal motor nucleus of the vagus.
The Effector (Efferent Arc - CAP): Efferent vagal fibers synapse in the celiac ganglion, communicating with the splenic nerve. This triggers a specialised subset of memory T-cells (ChAT+ T-cells) in the spleen to synthesise and release acetylcholine (ACh) (Rosas-Ballina et al., 2011).
The Brake: ACh binds to 𝛼 7 nicotinic acetylcholine receptors (𝛼 7nAChR) on splenic macrophages. This binding inhibits the NF-𝜅B pathway, rapidly halting the transcription and release of TNF-𝛼.

This model explicitly distinguishes between synchronous (Sync) events (e.g., action potentials, direct receptor-ligand binding) and asynchronous (Async) events (e.g., genetic transcription, cytokine synthesis, paracrine diffusion). It also highlights the critical integration points where signals cross distinct biological domains.
Key Architectural Distinctions Mapped:
Synchronous Points (Solid Arrows ->>): Represent immediate, millisecond-to-second physiological reactions. These include the rapid propagation of action potentials down nerve fibers, synaptic neurotransmitter release, and direct ligand-receptor binding events (like ACh hitting the α7nAChR).
Asynchronous Points (Dotted Arrows -->>): Represent molecular synthesis, paracrine diffusion, and genomic responses. These introduce biological "lag" into the system. For instance, the time it takes a macrophage to transcribe and translate TNF-α, or the time required for a T-cell to synthesise non-neuronal ACh, decoupled from the immediate speed of neural firing.
Integration Points: Explicitly marked where the medium of communication changes state. Point 1 converts immune molecules to neural electricity. Point 2 bridges the parasympathetic and sympathetic nervous systems. Point 3 translates a neural transmitter (Norepinephrine) into an immune-cell activation trigger.
The Simulation Design: Logic and Structure
A robust "toy model" reduces this pathway to a system of Ordinary Differential Equations (ODEs) that govern states over time (𝚝). In a Colab environment (typically using Python's scipy.integrate), the model conceptualises the dynamic tension between the inflammatory stimulus and the neural brake.
1. Modeling the Inflammatory Surge
The simulation begins with an inflammatory challenge, usually represented by an injection of LPS. The rate of TNF-𝛼 production relies on the concentration of activated macrophages (M𝛷act) and intrinsic decay kinetics:

(Where K𝑝 is the production rate constant and Kd is the clearance/degradation rate).
"The net rate at which TNF accumulates in the body is equal to the amount being produced by activated macrophages minus the amount that is naturally degrading or being cleared."
2. Conceptualising the Neural Brake (The CAP)
The true elegance of the model lies in how vagus nerve stimulation (VNS) alters macrophage activation. The model treats VNS as an exogenous input (U(𝚝)) that drives ACh concentration. The inhibitory effect of ACh on the macrophage operates via receptor-ligand biophysics, best modelled by a modified Hill equation indicating the fractional occupancy of the 𝛼 7nAChR:

In this logic:
The first term represents standard Michaelis-Menten-style activation by endotoxin.
The second term acts as the inhibitory multiplier. As ACh concentration (driven by VNS) approaches and exceeds the IC50 (the concentration required for 50% inhibition), the multiplier approaches zero, effectively shutting down M𝛷act and, consequently, TNF-𝛼 release.
Why This Toy Model Matters
Why build a simplified model in a Colab notebook? Because it allows us to perturb the system in ways that are ethically or physically impossible in vivo.
Through this structural logic, the simulation conceptualises dose-response non-linearities. Clinically, VNS is not a "more is better" therapy. Overtaxing the vagus nerve can lead to bradycardia or receptor desensitisation. The toy model allows bioengineers to adjust the simulated VNS pulse frequency and amplitude to find the theoretical "Goldilocks zone"—maximising the 𝛼7nAChR inhibitory multiplier while minimising off-target physiological stress.
Furthermore, it conceptualises the therapeutic window for bioelectronic medicine. By tweaking the delay parameters in the ODEs, the model demonstrates that VNS must be applied either preemptively or during the acute upward trajectory of the TNF-𝛼 curve to be effective.
Conclusion
The VNS-CAP Simulation is more than a string of Python logic; it is a biophysical translation of the body's internal peacekeeping mechanism. By mathematical abstraction of the 𝛼7nAChR receptor kinetics and splenic nerve firing rates, these toy models form the foundational blueprints for the next generation of neuromodulation devices aimed at treating rheumatoid arthritis, inflammatory bowel disease, and sepsis.
References:
Pavlov, V. A., & Tracey, K. J. (2012). The vagus nerve and the inflammatory reflex—linking immunity and metabolism. Nature Reviews Endocrinology, 8(12), 743-754. Provides the kinetic framing useful for setting parameters in computational models.
Tracey, K. J. (2002). The inflammatory reflex. Nature, 420(6917), 853-859. This foundational paper defines the macroscopic loop of the VNS-CAP.
Rosas-Ballina, M., et al. (2011). Acetylcholine-synthesising T cells relay neural signals in a vagus nerve circuit. Science, 334(6052), 98-101. Crucial for modeling the intermediate cellular steps in the spleen.


Comments