Heterogeneous firing rate response of mouse layer V pyramidal neurons in the fluctuation-driven regime

A mean-field model of gamma-frequency oscillations in networks of excitatory and inhibitory neurons

Gamma oscillations are widely seen in the cerebral cortex in different states of the wake-sleep cycle and are thought to play a role in sensory processing and cognition. Here, we study the emergence of gamma oscillations at two levels, in networks of spiking neurons, and a mean-field model. At the network level, we consider two different mechanisms to generate gamma oscillations and show that they are best seen if one takes into account the synaptic delay between neurons. At the mean-field level, we show that, by introducing delays, the mean-field can also produce gamma oscillations. The mean-field matches the mean activity of excitatory and inhibitory populations of the spiking network, as well as their oscillation frequencies, for both mechanisms. This mean-field model of gamma oscillations should be a useful tool to investigate large-scale interactions through gamma oscillations in the brain.

On the physiological and structural contributors to the overall balance of excitation and inhibition in local cortical networks

Overall balance of excitation and inhibition in cortical networks is central to their functionality and normal operation. Such orchestrated co-evolution of excitation and inhibition is established through convoluted local interactions between neurons, which are organized by specific network connectivity structures and are dynamically controlled by modulating synaptic activities. Therefore, identifying how such structural and physiological factors contribute to establishment of overall balance of excitation and inhibition is crucial in understanding the homeostatic plasticity mechanisms that regulate the balance. We use biologically plausible mathematical models to extensively study the effects of multiple key factors on overall balance of a network. We characterize a network's baseline balanced state by certain functional properties, and demonstrate how variations in physiological and structural parameters of the network deviate this balance and, in particular, result in transitions in spontaneous activity of the network to high-amplitude slow oscillatory regimes. We show that deviations from the reference balanced state can be continuously quantified by measuring the ratio of mean excitatory to mean inhibitory synaptic conductances in the network. Our results suggest that the commonly observed ratio of the number of inhibitory to the number of excitatory neurons in local cortical networks is almost optimal for their stability and excitability. Moreover, the values of inhibitory synaptic decay time constants and density of inhibitory-to-inhibitory network connectivity are critical to overall balance and stability of cortical networks. However, network stability in our results is sufficiently robust against modulations of synaptic quantal conductances, as required by their role in learning and memory. Our study based on extensive bifurcation analyses thus reveal the functional optimality and criticality of structural and physiological parameters in establishing the baseline operating state of local cortical networks.

Diversity-induced trivialization and resilience of neural dynamics

Heterogeneity is omnipresent across all living systems. Diversity enriches the dynamical repertoire of these systems but remains challenging to reconcile with their manifest robustness and dynamical persistence over time, a fundamental feature called resilience. To better understand the mechanism underlying resilience in neural circuits, we considered a nonlinear network model, extracting the relationship between excitability heterogeneity and resilience. To measure resilience, we quantified the number of stationary states of this network, and how they are affected by various control parameters. We analyzed both analytically and numerically gradient and non-gradient systems modeled as non-linear sparse neural networks evolving over long time scales. Our analysis shows that neuronal heterogeneity quenches the number of stationary states while decreasing the susceptibility to bifurcations: a phenomenon known as trivialization. Heterogeneity was found to implement a homeostatic control mechanism enhancing network resilience to changes in network size and connection probability by quenching the system's dynamic volatility.

A multi-layer mean-field model of the cerebellum embedding microstructure and population-specific dynamics

Mean-field (MF) models are computational formalism used to summarize in a few statistical parameters the salient biophysical properties of an inter-wired neuronal network. Their formalism normally incorporates different types of neurons and synapses along with their topological organization. MFs are crucial to efficiently implement the computational modules of large-scale models of brain function, maintaining the specificity of local cortical microcircuits. While MFs have been generated for the isocortex, they are still missing for other parts of the brain. Here we have designed and simulated a multi-layer MF of the cerebellar microcircuit (including Granule Cells, Golgi Cells, Molecular Layer Interneurons, and Purkinje Cells) and validated it against experimental data and the corresponding spiking neural network (SNN) microcircuit model. The cerebellar MF was built using a system of equations, where properties of neuronal populations and topological parameters are embedded in inter-dependent transfer functions. The model time constant was optimised using local field potentials recorded experimentally from acute mouse cerebellar slices as a template. The MF reproduced the average dynamics of different neuronal populations in response to various input patterns and predicted the modulation of the Purkinje Cells firing depending on cortical plasticity, which drives learning in associative tasks, and the level of feedforward inhibition. The cerebellar MF provides a computationally efficient tool for future investigations of the causal relationship between microscopic neuronal properties and ensemble brain activity in virtual brain models addressing both physiological and pathological conditions.

The h-current controls cortical recurrent network activity through modulation of dendrosomatic communication

A fundamental feature of the cerebral cortex is the ability to rapidly turn on and off maintained activity within ensembles of neurons through recurrent excitation balanced by inhibition. Here we demonstrate that reduction of the h-current, which is especially prominent in pyramidal cell dendrites, strongly increases the ability of local cortical networks to generate maintained recurrent activity. Reduction of the h-current resulted in hyperpolarization and increase in input resistance of both the somata and apical dendrites of layer 5 pyramidal cells, while strongly increasing the dendrosomatic transfer of low (<20 Hz) frequencies, causing an increased responsiveness to dynamic clamp-induced recurrent network-like activity injected into the dendrites and substantially increasing the duration of spontaneous Up states. We propose that modulation of the h-current may strongly control the ability of cortical networks to generate recurrent persistent activity and the formation and dissolution of neuronal ensembles.

On the physiological and structural contributors to the overall balance of excitation and inhibition in local cortical networks

Overall balance of excitation and inhibition in cortical networks is central to their functionality and normal operation. Such orchestrated co-evolution of excitation and inhibition is established through convoluted local interactions between neurons, which are organized by specific network connectivity structures and are dynamically controlled by modulating synaptic activities. Therefore, identifying how such structural and physiological factors contribute to establishment of overall balance of excitation and inhibition is crucial in understanding the homeostatic plasticity mechanisms that regulate the balance. We use biologically plausible mathematical models to extensively study the effects of multiple key factors on overall balance of a network. We characterize a network's baseline balanced state by certain functional properties, and demonstrate how variations in physiological and structural parameters of the network deviate this balance and, in particular, result in transitions in spontaneous activity of the network to high-amplitude slow oscillatory regimes. We show that deviations from the reference balanced state can be continuously quantified by measuring the ratio of mean excitatory to mean inhibitory synaptic conductances in the network. Our results suggest that the commonly observed ratio of the number of inhibitory to the number of excitatory neurons in local cortical networks is almost optimal for their stability and excitability. Moreover, the values of inhibitory synaptic decay time constants and density of inhibitory-to-inhibitory network connectivity are critical to overall balance and stability of cortical networks. However, network stability in our results is sufficiently robust against modulations of synaptic quantal conductances, as required by their role in learning and memory.

A comprehensive neural simulation of slow-wave sleep and highly responsive wakefulness dynamics

Hallmarks of neural dynamics during healthy human brain states span spatial scales from neuromodulators acting on microscopic ion channels to macroscopic changes in communication between brain regions. Developing a scale-integrated understanding of neural dynamics has therefore remained challenging. Here, we perform the integration across scales using mean-field modeling of Adaptive Exponential (AdEx) neurons, explicitly incorporating intrinsic properties of excitatory and inhibitory neurons. The model was run using The Virtual Brain (TVB) simulator, and is open-access in EBRAINS. We report that when AdEx mean-field neural populations are connected via structural tracts defined by the human connectome, macroscopic dynamics resembling human brain activity emerge. Importantly, the model can qualitatively and quantitatively account for properties of empirically observed spontaneous and stimulus-evoked dynamics in space, time, phase, and frequency domains. Large-scale properties of cortical dynamics are shown to emerge from both microscopic-scale adaptation that control transitions between wake-like to sleep-like activity, and the organization of the human structural connectome; together, they shape the spatial extent of synchrony and phase coherence across brain regions consistent with the propagation of sleep-like spontaneous traveling waves at intermediate scales. Remarkably, the model also reproduces brain-wide, enhanced responsiveness and capacity to encode information particularly during wake-like states, as quantified using the perturbational complexity index. The model was run using The Virtual Brain (TVB) simulator, and is open-access in EBRAINS. This approach not only provides a scale-integrated understanding of brain states and their underlying mechanisms, but also open access tools to investigate brain responsiveness, toward producing a more unified, formal understanding of experimental data from conscious and unconscious states, as well as their associated pathologies.

Noise Enhancement of Neural Information Processing

Cortical neurons in vivo function in highly fluctuating and seemingly noisy conditions, and the understanding of how information is processed in such complex states is still incomplete. In this perspective article, we first overview that an intense "synaptic noise" was measured first in single neurons, and computational models were built based on such measurements. Recent progress in recording techniques has enabled the measurement of highly complex activity in large numbers of neurons in animals and human subjects, and models were also built to account for these complex dynamics. Here, we attempt to link these two cellular and population aspects, where the complexity of network dynamics in awake cortex seems to link to the synaptic noise seen in single cells. We show that noise in single cells, in networks, or structural noise, all participate to enhance responsiveness and boost the propagation of information. We propose that such noisy states are fundamental to providing favorable conditions for information processing at large-scale levels in the brain, and may be involved in sensory perception.

State-dependent activity dynamics of hypothalamic stress effector neurons

The stress response necessitates an immediate boost in vital physiological functions from their homeostatic operation to an elevated emergency response. However, the neural mechanisms underlying this state-dependent change remain largely unknown. Using a combination of in vivo and ex vivo electrophysiology with computational modeling, we report that corticotropin releasing hormone (CRH) neurons in the paraventricular nucleus of the hypothalamus (PVN), the effector neurons of hormonal stress response, rapidly transition between distinct activity states through recurrent inhibition. Specifically, in vivo optrode recording shows that under non-stress conditions, CRH_PVN neurons often fire with rhythmic brief bursts (RB), which, somewhat counterintuitively, constrains firing rate due to long (~2 s) interburst intervals. Stressful stimuli rapidly switch RB to continuous single spiking (SS), permitting a large increase in firing rate. A spiking network model shows that recurrent inhibition can control this activity-state switch, and more broadly the gain of spiking responses to excitatory inputs. In biological CRH_PVN neurons ex vivo, the injection of whole-cell currents derived from our computational model recreates the in vivo-like switch between RB and SS, providing direct evidence that physiologically relevant network inputs enable state-dependent computation in single neurons. Together, we present a novel mechanism for state-dependent activity dynamics in CRH_PVN neurons.

A framework for preparing a stochastic nonlinear integrate-and-fire model for integrated information theory

This paper presents a framework for spiking neural networks to be prepared for the Integrated Information Theory (IIT) analysis, using a stochastic nonlinear integrate-and-fire model. The model includes the crucial dynamics of the all-or-none law and after-spike refractoriness. The noise is modelled as an additive term in the system's equations. By preparing the model for the IIT analysis, it is meant to determine the length of the analysis time-window and the transition probability distributions required for the IIT 3.0. To this end, a system of differential equations is proposed to estimate the time evolution of the system's mean and covariance. Assuming the binary Fired/Silent activity as the possible states of each neuron, an algorithm is proposed to calculate the required probability distributions. As long as the Fired/Silent probabilities are only concerned, the Gaussian density assumption with the estimated moments is a reasonable estimate. The synaptic inputs are treated as random variables with low variances to avoid the costs of conditioning on the system's past activities. The Monte-Carlo simulation is used to validate the estimation methods. To increase the reliability of the inductive inference behind the Monte-Carlo method, various stimulation protocols are applied to evoke the dynamics of the equations.

Optimal responsiveness and information flow in networks of heterogeneous neurons

Cerebral cortex is characterized by a strong neuron-to-neuron heterogeneity, but it is unclear what consequences this may have for cortical computations, while most computational models consider networks of identical units. Here, we study network models of spiking neurons endowed with heterogeneity, that we treat independently for excitatory and inhibitory neurons. We find that heterogeneous networks are generally more responsive, with an optimal responsiveness occurring for levels of heterogeneity found experimentally in different published datasets, for both excitatory and inhibitory neurons. To investigate the underlying mechanisms, we introduce a mean-field model of heterogeneous networks. This mean-field model captures optimal responsiveness and suggests that it is related to the stability of the spontaneous asynchronous state. The mean-field model also predicts that new dynamical states can emerge from heterogeneity, a prediction which is confirmed by network simulations. Finally we show that heterogeneous networks maximise the information flow in large-scale networks, through recurrent connections. We conclude that neuronal heterogeneity confers different responsiveness to neural networks, which should be taken into account to investigate their information processing capabilities.

The Spectrum of Asynchronous Dynamics in Spiking Networks as a Model for the Diversity of Non-rhythmic Waking States in the Neocortex

The awake cortex exhibits diverse non-rhythmic network states. However, how these states emerge and how each state impacts network function is unclear. Here, we demonstrate that model networks of spiking neurons with moderate recurrent interactions display a spectrum of non-rhythmic asynchronous dynamics based on the level of afferent excitation, from afferent input-dominated (AD) regimes, characterized by unbalanced synaptic currents and sparse firing, to recurrent input-dominated (RD) regimes, characterized by balanced synaptic currents and dense firing. The model predicted regime-specific relationships between different neural biophysical properties, which were all experimentally validated in the somatosensory cortex (S1) of awake mice. Moreover, AD regimes more precisely encoded spatiotemporal patterns of presynaptic activity, while RD regimes better encoded the strength of afferent inputs. These results provide a theoretical foundation for how recurrent neocortical circuits generate non-rhythmic waking states and how these different states modulate the processing of incoming information.

A mean-field approach to the dynamics of networks of complex neurons, from nonlinear Integrate-and-Fire to Hodgkin-Huxley models

We present a mean-field formalism able to predict the collective dynamics of large networks of conductance-based interacting spiking neurons. We apply this formalism to several neuronal models, from the simplest Adaptive Exponential Integrate-and-Fire model to the more complex Hodgkin-Huxley and Morris-Lecar models. We show that the resulting mean-field models are capable of predicting the correct spontaneous activity of both excitatory and inhibitory neurons in asynchronous irregular regimes, typical of cortical dynamics. Moreover, it is possible to quantitatively predict the population response to external stimuli in the form of external spike trains. This mean-field formalism therefore provides a paradigm to bridge the scale between population dynamics and the microscopic complexity of the individual cells physiology.NEW & NOTEWORTHY Population models are a powerful mathematical tool to study the dynamics of neuronal networks and to simulate the brain at macroscopic scales. We present a mean-field model capable of quantitatively predicting the temporal dynamics of a network of complex spiking neuronal models, from Integrate-and-Fire to Hodgkin-Huxley, thus linking population models to neurons electrophysiology. This opens a perspective on generating biologically realistic mean-field models from electrophysiological recordings.

Bridging Single Neuron Dynamics to Global Brain States

Biological neural networks produce information backgrounds of multi-scale spontaneous activity that become more complex in brain states displaying higher capacities for cognition, for instance, attentive awake versus asleep or anesthetized states. Here, we review brain state-dependent mechanisms spanning ion channel currents (microscale) to the dynamics of brain-wide, distributed, transient functional assemblies (macroscale). Not unlike how microscopic interactions between molecules underlie structures formed in macroscopic states of matter, using statistical physics, the dynamics of microscopic neural phenomena can be linked to macroscopic brain dynamics through mesoscopic scales. Beyond spontaneous dynamics, it is observed that stimuli evoke collapses of complexity, most remarkable over high dimensional, asynchronous, irregular background dynamics during consciousness. In contrast, complexity may not be further collapsed beyond synchrony and regularity characteristic of unconscious spontaneous activity. We propose that increased dimensionality of spontaneous dynamics during conscious states supports responsiveness, enhancing neural networks' emergent capacity to robustly encode information over multiple scales.

Biologically Realistic Mean-Field Models of Conductance-Based Networks of Spiking Neurons with Adaptation

Accurate population models are needed to build very large-scale neural models, but their derivation is difficult for realistic networks of neurons, in particular when nonlinear properties are involved, such as conductance-based interactions and spike-frequency adaptation. Here, we consider such models based on networks of adaptive exponential integrate-and-fire excitatory and inhibitory neurons. Using a master equation formalism, we derive a mean-field model of such networks and compare it to the full network dynamics. The mean-field model is capable of correctly predicting the average spontaneous activity levels in asynchronous irregular regimes similar to in vivo activity. It also captures the transient temporal response of the network to complex external inputs. Finally, the mean-field model is also able to quantitatively describe regimes where high- and low-activity states alternate (up-down state dynamics), leading to slow oscillations. We conclude that such mean-field models are biologically realistic in the sense that they can capture both spontaneous and evoked activity, and they naturally appear as candidates to build very large-scale models involving multiple brain areas.

Tune it in: mechanisms and computational significance of neuron-autonomous plasticity

The activity of a neural network is a result of synaptic signals that convey the communication between neurons and neuron-based intrinsic currents that determine the neuron’s input-output transfer function. Ample studies have demonstrated that cell-based excitability, and in particular intrinsic excitability, is modulated by learning and that these modifications play a key role in learning-related behavioral changes. The field of cell-based plasticity is largely growing, and it entails numerous experimental findings that demonstrate a large diversity of currents that are affected by learning. The diverse effect of learning on the neuron’s excitability emphasizes the need for a framework under which cell-based plasticity can be categorized to enable the assessment of the computational roles of the intrinsic modifications. We divide the domain of cell-based plasticity into three main categories, where the first category entails the currents that mediate the passive properties and single-spike generation, the second category entails the currents that mediate spike frequency adaptation, and the third category entails a novel learning-induced mechanism where all excitatory and inhibitory synapses double their strength. Curiously, this elementary division enables a natural categorization of the computational roles of these learning-induced plasticities. The computational roles are diverse and include modification of the neuronal mode of action, such as bursting, prolonged, and fast responsive; attention-like effect where the signal detection is improved; transfer of the network into an active state; biasing the competition for memory allocation; and transforming an environmental cue into a dominant cue and enabling a quicker formation of new memories.

Heterogeneity of the Axon Initial Segment in Interneurons and Pyramidal Cells of Rodent Visual Cortex

The microdomain that orchestrates action potential initiation in neurons is the axon initial segment (AIS). It has long been considered to be a rather homogeneous domain at the very proximal axon hillock with relatively stable length, particularly in cortical pyramidal cells. However, studies in other brain regions paint a different picture. In hippocampal CA1, up to 50% of axons emerge from basal dendrites. Further, in about 30% of thick-tufted layer V pyramidal neurons in rat somatosensory cortex, axons have a dendritic origin. Consequently, the AIS is separated from the soma. Recent in vitro and in vivo studies have shown that cellular excitability is a function of AIS length/position and somatodendritic morphology, undermining a potentially significant impact of AIS heterogeneity for neuronal function. We therefore investigated neocortical axon morphology and AIS composition, hypothesizing that the initial observation of seemingly homogeneous AIS is inadequate and needs to take into account neuronal cell types. Here, we biolistically transfected cortical neurons in organotypic cultures to visualize the entire neuron and classify cell types in combination with immunolabeling against AIS markers. Using confocal microscopy and morphometric analysis, we investigated axon origin, AIS position, length, diameter as well as distance to the soma. We find a substantial AIS heterogeneity in visual cortical neurons, classified into three groups: (I) axons with somatic origin with proximal AIS at the axon hillock; (II) axons with somatic origin with distal AIS, with a discernible gap between the AIS and the soma; and (III) axons with dendritic origin (axon-carrying dendrite cell, AcD cell) and an AIS either starting directly at the axon origin or more distal to that point. Pyramidal cells have significantly longer AIS than interneurons. Interneurons with vertical columnar axonal projections have significantly more distal AIS locations than all other cells with their prevailing phenotype as an AcD cell. In contrast, neurons with perisomatic terminations display most often an axon originating from the soma. Our data contribute to the emerging understanding that AIS morphology is highly variable, and potentially a function of the cell type.

Refractoriness Accounts for Variable Spike Burst Responses in Somatosensory Cortex

Neurons in the primary somatosensory cortex (S1) respond to peripheral stimulation with synchronized bursts of spikes, which lock to the macroscopic 600-Hz EEG waves. The mechanism of burst generation and synchronization in S1 is not yet understood. Using models of single-neuron responses fitted to unit recordings from macaque monkeys, we show that these synchronized bursts are the consequence of correlated synaptic inputs combined with a refractory mechanism. In the presence of noise these models reproduce also the observed trial-to-trial response variability, where individual bursts represent one of many stereotypical temporal spike patterns. When additional slower and global excitability fluctuations are introduced the single-neuron spike patterns are correlated with the population activity, as demonstrated in experimental data. The underlying biophysical mechanism of S1 responses involves thalamic inputs arriving through depressing synapses to cortical neurons in a high-conductance state. Our findings show that a simple feedforward processing of peripheral inputs could give rise to neuronal responses with nontrivial temporal and population statistics. We conclude that neural systems could use refractoriness to encode variable cortical states into stereotypical short-term spike patterns amenable to processing at neuronal time scales (tens of milliseconds).

Sleep regulation of the distribution of cortical firing rates

Sleep is thought to mediate both mnemonic and homeostatic functions. However, the mechanism by which this brain state can simultaneously implement the 'selective' plasticity needed to consolidate novel memory traces and the 'general' plasticity necessary to maintain a well-functioning neuronal system is unclear. Recent findings show that both of these functions differentially affect neurons based on their intrinsic firing rate, a ubiquitous neuronal heterogeneity. Furthermore, they are both implemented by the NREM slow oscillation, which also distinguishes neurons based on firing rate during sequential activity at the DOWN→UP transition. These findings suggest a mechanism by which spiking activity during the slow oscillation acts to maintain network statistics that promote a skewed distribution of neuronal firing rates, and perturbation of that activity by hippocampal replay acts to integrate new memory traces into the existing cortical network.

Heterogeneous firing responses predict diverse couplings to presynaptic activity in mice layer V pyramidal neurons

In this study, we present a theoretical framework combining experimental characterizations and analytical calculus to capture the firing rate input-output properties of single neurons in the fluctuation-driven regime. Our framework consists of a two-step procedure to treat independently how the dendritic input translates into somatic fluctuation variables, and how the latter determine action potential firing. We use this framework to investigate the functional impact of the heterogeneity in firing responses found experimentally in young mice layer V pyramidal cells. We first design and calibrate in vitro a simplified morphological model of layer V pyramidal neurons with a dendritic tree following Rall's branching rule. Then, we propose an analytical derivation for the membrane potential fluctuations at the soma as a function of the properties of the synaptic input in dendrites. This mathematical description allows us to easily emulate various forms of synaptic input: either balanced, unbalanced, synchronized, purely proximal or purely distal synaptic activity. We find that those different forms of dendritic input activity lead to various impact on the somatic membrane potential fluctuations properties, thus raising the possibility that individual neurons will differentially couple to specific forms of activity as a result of their different firing response. We indeed found such a heterogeneous coupling between synaptic input and firing response for all types of presynaptic activity. This heterogeneity can be explained by different levels of cellular excitability in the case of the balanced, unbalanced, synchronized and purely distal activity. A notable exception appears for proximal dendritic inputs: increasing the input level can either promote firing response in some cells, or suppress it in some other cells whatever their individual excitability. This behavior can be explained by different sensitivities to the speed of the fluctuations, which was previously associated to different levels of sodium channel inactivation and density. Because local network connectivity rather targets proximal dendrites, our results suggest that this aspect of biophysical heterogeneity might be relevant to neocortical processing by controlling how individual neurons couple to local network activity.

Speed hysteresis and noise shaping of traveling fronts in neural fields: role of local circuitry and nonlocal connectivity

Neural field models are powerful tools to investigate the richness of spatiotemporal activity patterns like waves and bumps, emerging from the cerebral cortex. Understanding how spontaneous and evoked activity is related to the structure of underlying networks is of central interest to unfold how information is processed by these systems. Here we focus on the interplay between local properties like input-output gain function and recurrent synaptic self-excitation of cortical modules, and nonlocal intermodular synaptic couplings yielding to define a multiscale neural field. In this framework, we work out analytic expressions for the wave speed and the stochastic diffusion of propagating fronts uncovering the existence of an optimal balance between local and nonlocal connectivity which minimizes the fluctuations of the activation front propagation. Incorporating an activity-dependent adaptation of local excitability further highlights the independent role that local and nonlocal connectivity play in modulating the speed of propagation of the activation and silencing wavefronts, respectively. Inhomogeneities in space of local excitability give raise to a novel hysteresis phenomenon such that the speed of waves traveling in opposite directions display different velocities in the same location. Taken together these results provide insights on the multiscale organization of brain slow-waves measured during deep sleep and anesthesia.

Is realistic neuronal modeling realistic?

Scientific models are abstractions that aim to explain natural phenomena. A successful model shows how a complex phenomenon arises from relatively simple principles while preserving major physical or biological rules and predicting novel experiments. A model should not be a facsimile of reality; it is an aid for understanding it. Contrary to this basic premise, with the 21st century has come a surge in computational efforts to model biological processes in great detail. Here we discuss the oxymoronic, realistic modeling of single neurons. This rapidly advancing field is driven by the discovery that some neurons don't merely sum their inputs and fire if the sum exceeds some threshold. Thus researchers have asked what are the computational abilities of single neurons and attempted to give answers using realistic models. We briefly review the state of the art of compartmental modeling highlighting recent progress and intrinsic flaws. We then attempt to address two fundamental questions. Practically, can we realistically model single neurons? Philosophically, should we realistically model single neurons? We use layer 5 neocortical pyramidal neurons as a test case to examine these issues. We subject three publically available models of layer 5 pyramidal neurons to three simple computational challenges. Based on their performance and a partial survey of published models, we conclude that current compartmental models are ad hoc, unrealistic models functioning poorly once they are stretched beyond the specific problems for which they were designed. We then attempt to plot possible paths for generating realistic single neuron models.