How Regional Neurochemistry Shapes Large-Scale Brain Activity Patterns

from work published in PNAS, reframes a long-standing assumption in systems neuroscience: that the cortex speaks with one voice. By layering receptor-density maps onto the brain's structural connections across 68 regions, the authors show that regional neurochemistry shapes — quite literally — the rhythm of large-scale activity.

A cortex with local accents

Built by Leonardo Dalla Porta and Maria V. Sanchez-Vives at IDIBAPS, the model runs on The Virtual Brain (TVB), an open-source whole-brain simulator housed within the EBRAINS research infrastructure, with support from the EBRAINS 2.0 project and The Virtual Brain Twin Project. Its distinguishing move is architectural. Rather than treating muscarinic acetylcholine receptors as a uniform substrate, the team loads detailed maps of their density across 68 cortical regions and lets each area modulate its own dynamics accordingly. The cortex, in effect, inherits a topography of neurochemical accent — a spatial structure of sensitivities that, until now, most large-scale models had quietly smoothed away.

What spatial heterogeneity actually buys you

The team compared this regionally structured model against a uniform baseline while simulating transitions from wakefulness to sleep. Two outcomes stood out: inter-regional coordination rose, and information flow across the network became more efficient. The simulation also spontaneously produced a phenomenon long observed in real brains but rarely captured by models — localized slow waves appearing in one region while the surrounding cortex remains in an awake-like state. This signature has been recorded around brain lesions, during sleep deprivation, and during attentional lapses. The new work suggests it can emerge from receptor-density gradients themselves, without needing ad hoc inputs or lesion-specific tuning.

Acetylcholine is no longer a single knob

Perhaps the more consequential implication lies one layer down. Neuromodulators such as acetylcholine do not paint the cortex with a uniform brush; they land on a receptor landscape whose topography determines local effect. The same chemical broadcast, in other words, is heard differently depending on where the listener sits. For anyone designing experiments or simulations around state transitions — wakefulness, anesthesia, disorders of consciousness — that distinction matters operationally. Receptor maps belong inside the model from the first iteration, not as a post-hoc correction applied once the dynamics fail to match empirical recordings.

A principle for the bench

The takeaway is architectural rather than philosophical: molecular heterogeneity is not a refinement, it is load-bearing. Flatten it, and the very transitions the model is meant to explain — sleep onset, attentional lapses, lesion-adjacent slow waves — quietly disappear into the noise floor. For groups working with whole-brain simulations, the practical step is to treat receptor density as a first-class input alongside the structural connectome. The difference shows up directly in coordination metrics, in the emergence of localized slow-wave states, and in the model's ability to generate testable predictions about how a single neuromodulator can produce several qualitatively different cortical regimes at once.