Brain Simulation Approaches

Levels of Biological Detail

The first decision in any brain simulation is how much biological detail to include. This choice is not merely technical; it reflects a scientific hypothesis about which features of neural computation are essential and which can be safely abstracted away.

At the most detailed level, biophysical models simulate individual ion channels, the protein pores in neuronal membranes that control the flow of sodium, potassium, calcium, and chloride ions. The Hodgkin-Huxley model, published in 1952 and still the gold standard for this level of detail, describes each ion channel as a voltage-dependent conductance with specific activation and inactivation kinetics. A single cortical neuron may express dozens of distinct ion channel types distributed across its soma, axon, and dendritic tree, and a biophysically detailed model must capture all of these to reproduce the neuron electrical behavior accurately.

One step up in abstraction, compartmental models represent the spatial structure of neurons by dividing them into electrical compartments, each described by a set of differential equations for membrane voltage and ionic currents. The NEURON simulator, developed by Michael Hines and John Moore at Yale, is the standard tool for this approach. A typical cortical pyramidal neuron model might include hundreds of compartments to capture the complex interactions between synaptic inputs arriving at different dendritic locations.

Further abstraction produces point neuron models, which collapse the entire neuron into a single electrical compartment. The integrate-and-fire model, the most common variety, simply accumulates incoming currents until the membrane voltage crosses a threshold, at which point the neuron emits a spike and resets. The Izhikevich model adds two differential equations that can reproduce over 20 distinct firing patterns observed in biological neurons while remaining computationally cheap enough to simulate networks of millions of neurons on a single GPU.

At the coarsest level, rate models and mean-field models describe the average activity of populations of neurons rather than individual cells. These models are appropriate for studying large-scale brain dynamics, such as the patterns of activity that underlie sleep, arousal, or cognitive states, where the behavior of individual neurons is less important than the collective dynamics of neural populations.

Major Simulation Platforms

Several software platforms have become standard tools for brain simulation research, each optimized for different levels of detail and different scales of simulation.

NEURON remains the most widely used platform for biophysically detailed simulations. It supports compartmental models with arbitrary morphology and ion channel distributions, and it can simulate networks of detailed neurons connected by chemical and electrical synapses. The CoreNEURON extension enables GPU acceleration, making it possible to simulate networks of hundreds of thousands of biophysically detailed neurons on modern supercomputers.

NEST (Neural Simulation Technology) is optimized for large-scale spiking neural network simulations using point neuron models. NEST can simulate networks of millions to billions of neurons on distributed computing clusters, making it the tool of choice for studying brain-scale dynamics. The NEST simulator was used in the Human Brain Project to simulate cortical circuits at scales approaching the full mouse cortex.

Brian is a Python-based simulator that prioritizes flexibility and ease of use. Researchers define their neuron and synapse models using mathematical equations directly in Python code, and Brian compiles these to efficient C++ or GPU code. This makes it particularly useful for rapid prototyping of novel neuron models and learning rules.

GENESIS (GEneral NEural SImulation System) was one of the earliest neural simulation platforms and is still used for detailed models of subcortical circuits and individual neurons with complex dendritic processing.

The Blue Brain Project Approach

The Blue Brain Project at EPFL represents the most ambitious attempt to build biologically detailed brain simulations constrained by experimental data at every level. The project approach begins with the reconstruction of individual neuron morphologies from microscopy data, continues with the fitting of ion channel models to electrophysiological recordings, and culminates in the assembly of circuit models whose connectivity follows rules derived from paired recording experiments and anatomical tracing studies.

The project first major achievement was a simulation of a cortical microcircuit containing approximately 31,000 neurons and 37 million synapses, representing a small volume of rat somatosensory cortex. Each neuron in the simulation had a detailed three-dimensional morphology with hundreds of compartments and multiple ion channel types. The simulation reproduced several experimentally observed phenomena, including spontaneous activity patterns, the response to simulated sensory input, and the effects of neuromodulators.

The Blue Brain approach has been both celebrated and criticized. Supporters argue that biologically detailed models are necessary to understand how neural circuits actually compute, and that only by matching the biological detail can we identify which features are computationally important. Critics counter that the computational cost is prohibitive, that the models contain too many free parameters to be truly predictive, and that simpler models can often reproduce the same phenomena with less effort and greater interpretability.

Large-Scale Network Simulations

At the opposite end of the detail spectrum, large-scale network simulations sacrifice cellular detail to study the emergent dynamics of brain-sized neural networks. These simulations typically use point neuron models (integrate-and-fire or Izhikevich neurons) connected according to statistical rules derived from neuroanatomy.

The largest spiking neural network simulations to date have reached scales of billions of neurons. A 2018 simulation on the Japanese K computer modeled 1.73 billion neurons with 10.4 trillion synapses, representing roughly 1% of the human brain neural network, though with vastly simplified neuron models. Even at this scale, the simulation ran at roughly 1/1000th of biological real time, meaning that one second of brain activity required roughly 17 minutes of supercomputer time.

More recent work has focused on making large-scale simulations faster and more biologically realistic by incorporating heterogeneous neuron types, distance-dependent connectivity, and realistic patterns of excitation and inhibition. The goal is not to simulate every neuron in the brain but to capture the organizational principles that produce brain-like dynamics, such as the balance between excitation and inhibition, the formation of cell assemblies, and the generation of oscillatory rhythms at different frequencies.

Multi-Scale Simulation

The most promising current direction in brain simulation is multi-scale modeling, which combines detailed biophysical models at the cellular level with more abstract models at larger scales. The idea is that detailed models can identify which cellular properties are functionally important, and these properties can then be incorporated into simplified models that are computationally feasible at larger scales.

Multi-scale simulation requires solving a fundamental mathematical problem: how to derive the macroscopic behavior of a neural population from the microscopic dynamics of individual neurons and synapses. Mean-field theory, borrowed from statistical physics, provides one approach by describing the average behavior of large populations of neurons in terms of a few macroscopic variables like mean firing rate and mean membrane voltage. These mean-field descriptions can then be used as building blocks for whole-brain simulations that capture the dynamics of interacting brain regions.

The Virtual Brain (TVB) project takes this approach to its logical conclusion by simulating the dynamics of the entire human brain using coupled mean-field models for each brain region, connected according to the structural connectivity measured by diffusion MRI tractography. TVB simulations can reproduce patterns of resting-state functional connectivity, the spatiotemporal dynamics of epileptic seizures, and the effects of virtual lesions, making them useful tools for clinical neuroscience as well as basic research.

Validation and the Reality Check

The central challenge in brain simulation is validation: how do you know if your simulation is right? A simulation can produce activity that looks brain-like without actually implementing the same computations that the real brain performs. This problem is particularly acute for large-scale simulations, where the number of adjustable parameters is enormous and the experimental data available for comparison is limited.

Current validation approaches include comparing simulated activity patterns to electrophysiological recordings (do the simulated neurons fire at realistic rates and with realistic temporal patterns?), comparing simulated network dynamics to functional imaging data (does the simulated brain show realistic patterns of correlated activity?), and testing whether simulated circuits can perform the behavioral tasks that the real circuits support (can a simulated hippocampal circuit learn spatial maps?).

Despite decades of progress, no brain simulation has been fully validated against biological data at any scale. Even the C. elegans simulation, with its complete wiring diagram of just 302 neurons, cannot yet reproduce the full behavioral repertoire of the living worm. This humbling fact underscores both the difficulty of the problem and the enormous amount that remains unknown about how neural circuits actually compute.