Multiscale Computer Model of the Spinal Dorsal Horn Reveals Changes in Network Processing Associated with Chronic Pain

Large-Scale Mechanistic Models of Brain Circuits with Biophysically and Morphologically Detailed Neurons

Understanding the brain requires studying its multiscale interactions from molecules to networks. The increasing availability of large-scale datasets detailing brain circuit composition, connectivity, and activity is transforming neuroscience. However, integrating and interpreting this data remains challenging. Concurrently, advances in supercomputing and sophisticated modeling tools now enable the development of highly detailed, large-scale biophysical circuit models. These mechanistic multiscale models offer a method to systematically integrate experimental data, facilitating investigations into brain structure, function, and disease. This review, based on a Society for Neuroscience 2024 MiniSymposium, aims to disseminate recent advances in large-scale mechanistic modeling to the broader community. It highlights (1) examples of current models for various brain regions developed through experimental data integration; (2) their predictive capabilities regarding cellular and circuit mechanisms underlying experimental recordings (e.g., membrane voltage, spikes, local-field potential, electroencephalography/magnetoencephalography) and brain function; and (3) their use in simulating biomarkers for brain diseases like epilepsy, depression, schizophrenia, and Parkinson's, aiding in understanding their biophysical underpinnings and developing novel treatments. The review showcases state-of-the-art models covering hippocampus, somatosensory, visual, motor, auditory cortical, and thalamic circuits across species. These models predict neural activity at multiple scales and provide insights into the biophysical mechanisms underlying sensation, motor behavior, brain signals, neural coding, disease, pharmacological interventions, and neural stimulation. Collaboration with experimental neuroscientists and clinicians is essential for the development and validation of these models, particularly as datasets grow. Hence, this review aims to foster interest in detailed brain circuit models, leading to cross-disciplinary collaborations that accelerate brain research.

A Visual and Narrative Timeline Review of Spinal Cord Stimulation Technology and US Food and Drug Administration Milestones

OBJECTIVES

The aim of this study was to present key technologic and regulatory milestones in spinal cord stimulation (SCS) for managing chronic pain on a narrative timeline with visual representation, relying on original sources to the extent possible.

MATERIALS AND METHODS

We identified technical advances in SCS that facilitated and enhanced treatment on the basis of scientific publications and approvals from the United States (US) Food and Drug Administration (FDA). We presented milestones limited to first use in key indications and in the context of new technology validation. We focused primarily on pain management, but other indications (eg, motor disorder in multiple sclerosis) were included when they affected technology development.

RESULTS

We developed a comprehensive visual and narrative timeline of SCS technology and US FDA milestones. Since its conception in the 1960s, the science and technology of SCS neuromodulation have continuously evolved. Advances span lead design (from paddle-type to percutaneous, and increased electrode contacts) and stimulator technology (from wireless power to internally powered and rechargeable, with miniaturized components, and programmable multichannel devices), with expanding stimulation program flexibility (such as burst and kilohertz stimulation frequencies), as well as usage features (such as remote programming and magnetic resonance imaging conditional compatibility).

CONCLUSIONS

This timeline represents the evolution of SCS technology alongside expanding FDA-approved indications for use.

Parvalbumin gates chronic pain through the modulation of firing patterns in inhibitory neurons

Spinal cord dorsal horn inhibition is critical to the processing of sensory inputs, and its impairment leads to mechanical allodynia. How this decreased inhibition occurs and whether its restoration alleviates allodynic pain are poorly understood. Here, we show that a critical step in the loss of inhibitory tone is the change in the firing pattern of inhibitory parvalbumin (PV)-expressing neurons (PVNs). Our results show that PV, a calcium-binding protein, controls the firing activity of PVNs by enabling them to sustain high-frequency tonic firing patterns. Upon nerve injury, PVNs transition to adaptive firing and decrease their PV expression. Interestingly, decreased PV is necessary and sufficient for the development of mechanical allodynia and the transition of PVNs to adaptive firing. This transition of the firing pattern is due to the recruitment of calcium-activated potassium (SK) channels, and blocking them during chronic pain restores normal tonic firing and alleviates chronic pain. Our findings indicate that PV is essential for controlling the firing pattern of PVNs and for preventing allodynia. Developing approaches to manipulate these mechanisms may lead to different strategies for chronic pain relief.

Mechanisms for dysregulation of excitatory-inhibitory balance underlying allodynia in dorsal horn neural subcircuits

Chronic pain is a wide-spread condition that is debilitating and expensive to manage, costing the United States alone around $600 billion in 2010. In a common type of chronic pain called allodynia, non-painful stimuli produce painful responses with highly variable presentations across individuals. While the specific mechanisms remain unclear, allodynia is hypothesized to be caused by the dysregulation of excitatory-inhibitory (E-I) balance in pain-processing neural circuitry in the dorsal horn of the spinal cord. In this work, we analyze biophysically-motivated subcircuit structures that represent common motifs in neural circuits in layers I-II of the dorsal horn. These circuits are hypothesized to be part of the neural pathways that mediate two different types of allodynia: static and dynamic. We use neural firing rate models to describe the activity of populations of excitatory and inhibitory interneurons within each subcircuit. By accounting for experimentally-observed responses under healthy conditions, we specify model parameters defining populations of subcircuits that yield typical behavior under normal conditions. Then, we implement a sensitivity analysis approach to identify the mechanisms most likely to cause allodynia-producing dysregulation of the subcircuit's E-I signaling. We find that disruption of E-I balance generally occurs either due to downregulation of inhibitory signaling so that excitatory neurons are "released" from inhibitory control, or due to upregulation of excitatory neuron responses so that excitatory neurons "escape" their inhibitory control. Which of these mechanisms is most likely to occur, the subcircuit components involved in the mechanism, and the proportion of subcircuits exhibiting the mechanism can vary depending on the subcircuit structure. These results suggest specific hypotheses about diverse mechanisms that may be most likely responsible for allodynia, thus offering predictions for the high interindividual variability observed in allodynia and identifying targets for further experimental studies on the underlying mechanisms of this chronic pain condition.

Mechanism of an Intrinsic Oscillation in Rat Geniculate Interneurons

Depolarizing current injections produced a rhythmic bursting of action potentials - a bursting oscillation - in a set of local interneurons in the lateral geniculate nucleus (LGN) of rats. The current dynamics underlying this firing pattern have not been determined, though this cell type constitutes an important cellular component of thalamocortical circuitry, and contributes to both pathologic and non-pathologic brain states. We thus investigated the source of the bursting oscillation using pharmacological manipulations in LGN slices in vitro and in silico. 1. Selective blockade of calcium channel subtypes revealed that high-threshold calcium currentsI L andI P contributed strongly to the oscillation. 2. Increased extracellular K+ concentration (decreased K+currents) eliminated the oscillation. 3. Selective blockade of K+ channel subtypes demonstrated that the calcium-sensitive potassium current (I A H P ) was of primary importance. A morphologically simplified, multicompartment model of the thalamic interneuron characterized the oscillation as follows: 1. The low-threshold calcium currentI T provided the strong initial burst characteristic of the oscillation. 2. Alternating fluxes through high-threshold calcium channels andI A H P then provided the continuing oscillation's burst and interburst periods respectively. This interplay betweenI L andI A H P contrasts with the current dynamics underlying oscillations in thalamocortical and reticularis neurons, which primarily involveI T andI H , orI T andI A H P respectively. These findings thus point to a novel electrophysiological mechanism for generating intrinsic oscillations in a major thalamic cell type. Because local interneurons can sculpt the behavior of thalamocortical circuits, these results suggest new targets for the manipulation of ascending thalamocortical network activity.

Similar excitability through different sodium channels and implications for the analgesic efficacy of selective drugs

Nociceptive sensory neurons convey pain-related signals to the CNS using action potentials. Loss-of-function mutations in the voltage-gated sodium channel NaV1.7 cause insensitivity to pain (presumably by reducing nociceptor excitability) but clinical trials seeking to treat pain by inhibiting NaV1.7 pharmacologically have struggled. This may reflect the variable contribution of NaV1.7 to nociceptor excitability. Contrary to claims that NaV1.7 is necessary for nociceptors to initiate action potentials, we show that nociceptors can achieve similar excitability using different combinations of NaV1.3, NaV1.7, and NaV1.8. Selectively blocking one of those NaV subtypes reduces nociceptor excitability only if the other subtypes are weakly expressed. For example, excitability relies on NaV1.8 in acutely dissociated nociceptors but responsibility shifts to NaV1.7 and NaV1.3 by the fourth day in culture. A similar shift in NaV dependence occurs in vivo after inflammation, impacting ability of the NaV1.7-selective inhibitor PF-05089771 to reduce pain in behavioral tests. Flexible use of different NaV subtypes exemplifies degeneracy - achieving similar function using different components - and compromises reliable modulation of nociceptor excitability by subtype-selective inhibitors. Identifying the dominant NaV subtype to predict drug efficacy is not trivial. Degeneracy at the cellular level must be considered when choosing drug targets at the molecular level.

Multimodal sensory control of motor performance by glycinergic interneurons of the mouse spinal cord deep dorsal horn

Sensory feedback is integral for contextually appropriate motor output, yet the neural circuits responsible remain elusive. Here, we pinpoint the medial deep dorsal horn of the mouse spinal cord as a convergence point for proprioceptive and cutaneous input. Within this region, we identify a population of tonically active glycinergic inhibitory neurons expressing parvalbumin. Using anatomy and electrophysiology, we demonstrate that deep dorsal horn parvalbumin-expressing interneuron (dPV) activity is shaped by convergent proprioceptive, cutaneous, and descending input. Selectively targeting spinal dPVs, we reveal their widespread ipsilateral inhibition onto pre-motor and motor networks and demonstrate their role in gating sensory-evoked muscle activity using electromyography (EMG) recordings. dPV ablation altered limb kinematics and step-cycle timing during treadmill locomotion and reduced the transitions between sub-movements during spontaneous behavior. These findings reveal a circuit basis by which sensory convergence onto dorsal horn inhibitory neurons modulates motor output to facilitate smooth movement and context-appropriate transitions.

Finding new analgesics: Computational pharmacology faces drug discovery challenges

Despite the worldwide prevalence and huge burden of pain, pain is an undertreated phenomenon. Currently used analgesics have several limitations regarding their efficacy and safety. The discovery of analgesics possessing a novel mechanism of action has faced multiple challenges, including a limited understanding of biological processes underpinning pain and analgesia and poor animal-to-human translation. Computational pharmacology is currently employed to face these challenges. In this review, we discuss the theory, methods, and applications of computational pharmacology in pain research. Computational pharmacology encompasses a wide variety of theoretical concepts and practical methodological approaches, with the overall aim of gaining biological insight through data acquisition and analysis. Data are acquired from patients or animal models with pain or analgesic treatment, at different levels of biological organization (molecular, cellular, physiological, and behavioral). Distinct methodological algorithms can then be used to analyze and integrate data. This helps to facilitate the identification of biological molecules and processes associated with pain phenotype, build quantitative models of pain signaling, and extract translatable features between humans and animals. However, computational pharmacology has several limitations, and its predictions can provide false positive and negative findings. Therefore, computational predictions are required to be validated experimentally before drawing solid conclusions. In this review, we discuss several case study examples of combining and integrating computational tools with experimental pain research tools to meet drug discovery challenges.

A biophysically comprehensive model of urothelial afferent neurons: implications for sensory signalling in urinary bladder

The urothelium is the innermost layer of the bladder wall; it plays a pivotal role in bladder sensory transduction by responding to chemical and mechanical stimuli. The urothelium also acts as a physical barrier between urine and the outer layers of the bladder wall. There is intricate sensory communication between the layers of the bladder wall and the neurons that supply the bladder, which eventually translates into the regulation of mechanical activity. In response to natural stimuli, urothelial cells release substances such as ATP, nitric oxide (NO), substance P, acetylcholine (ACh), and adenosine. These act on adjacent urothelial cells, myofibroblasts, and urothelial afferent neurons (UAN), controlling the contractile activity of the bladder. There is rising evidence on the importance of urothelial sensory signalling, yet a comprehensive understanding of the functioning of the urothelium-afferent neurons and the factors that govern it remains elusive to date. Until now, the biophysical studies done on UAN have been unable to provide adequate information on the ion channel composition of the neuron, which is paramount to understanding the electrical functioning of the UAN and, by extension, afferent signalling. To this end, we have attempted to model UAN to decipher the ionic mechanisms underlying the excitability of the UAN. In contrast to previous models, our model was built and validated using morphological and biophysical properties consistent with experimental findings for the UAN. The model included all the channels thus far known to be expressed in UAN, including; voltage-gated sodium and potassium channels, N, L, T, P/Q, R-type calcium channels, large-conductance calcium-dependent potassium (BK) channels, small conductance calcium-dependent (SK) channels, Hyperpolarisation activated cation (HCN) channels, transient receptor potential melastatin (TRPM8), transient receptor potential vanilloid (TRPV1) channel, calcium-activated chloride(CaCC) channels, and internal calcium dynamics. Our UAN model a) was constrained as far as possible by experimental data from the literature for the channels and the spiking activity, b) was validated by reproducing the experimental responses to current-clamp and voltage-clamp protocols c) was used as a base for modelling the non-urothelial afferent neurons (NUAN). Using our models, we also gained insights into the variations in ion channels between UAN and NUAN neurons.

Network model of nociceptive processing in the superficial spinal dorsal horn reveals mechanisms of hyperalgesia, allodynia, and spinal cord stimulation

The spinal dorsal horn (DH) processes sensory information and plays a key role in transmitting nociception to supraspinal centers. Loss of DH inhibition during neuropathic pain unmasks a pathway from nonnociceptive Aβ-afferent inputs to superficial dorsal horn (SDH) nociceptive-specific (NS) projection neurons, and this change may contribute to hyperalgesia and allodynia. We developed and validated a computational model of SDH neuronal circuitry that links nonnociceptive Aβ-afferent inputs in lamina II/III to a NS projection neuron in lamina I via a network of excitatory interneurons. The excitatory pathway and the NS projection neuron were in turn gated by inhibitory interneurons with connections based on prior patch-clamp recordings. Changing synaptic weights in the computational model to replicate neuropathic pain states unmasked a low-threshold excitatory pathway to NS neurons similar to experimental recordings. Spinal cord stimulation (SCS) is an effective therapy for neuropathic pain, and accumulating experimental evidence indicates that NS neurons in the SDH also respond to SCS. Accounting for these responses may inform therapeutic improvements, and we quantified responses to SCS in the SDH network model and examined the role of different modes of inhibitory control in modulating NS neuron responses to SCS. We combined the SDH network model with a previously published model of the deep dorsal horn (DDH) and identified optimal stimulation frequencies across different neuropathic pain conditions. Finally, we found that SCS-generated inhibition did not completely suppress model NS activity during simulated pinch inputs, providing an explanation of why SCS does not eliminate acute pain.NEW & NOTEWORTHY Chronic pain is a severe public health problem that reduces the quality of life for those affected and exacts an enormous socio-economic burden worldwide. Spinal cord stimulation (SCS) is an effective treatment for chronic pain, but SCS efficacy has not significantly improved over time, in part because the mechanisms of action remain unclear. Most preclinical studies investigating pain and SCS mechanisms have focused on the responses of deep dorsal horn (DDH) neurons, but neural networks in the superficial dorsal horn (SDH) are also important for processing nociceptive information. This work synthesizes heterogeneous experimental recordings from the SDH into a computational model that replicates experimental responses and that can be used to quantify neuronal responses to SCS under neuropathic pain conditions.

Modulation of SK Channels via Calcium Buffering Tunes Intrinsic Excitability of Parvalbumin Interneurons in Neuropathic Pain: A Computational and Experimental Investigation

Parvalbumin-expressing interneurons (PVINs) play a crucial role within the dorsal horn of the spinal cord by preventing touch inputs from activating pain circuits. In both male and female mice, nerve injury decreases PVINs' output via mechanisms that are not fully understood. In this study, we show that PVINs from nerve-injured male mice change their firing pattern from tonic to adaptive. To examine the ionic mechanisms responsible for this decreased output, we used a reparametrized Hodgkin-Huxley type model of PVINs, which predicted (1) the firing pattern transition is because of an increased contribution of small conductance calcium-activated potassium (SK) channels, enabled by (2) impairment in intracellular calcium buffering systems. Analyzing the dynamics of the Hodgkin-Huxley type model further demonstrated that a generalized Hopf bifurcation differentiates the two types of state transitions observed in the transient firing of PVINs. Importantly, this predicted mechanism holds true when we embed the PVIN model within the neuronal circuit model of the spinal dorsal horn. To experimentally validate this hypothesized mechanism, we used pharmacological modulators of SK channels and demonstrated that (1) tonic firing PVINs from naive male mice become adaptive when exposed to an SK channel activator, and (2) adapting PVINs from nerve-injured male mice return to tonic firing on SK channel blockade. Our work provides important insights into the cellular mechanism underlying the decreased output of PVINs in the spinal dorsal horn after nerve injury and highlights potential pharmacological targets for new and effective treatment approaches to neuropathic pain.SIGNIFICANCE STATEMENT Parvalbumin-expressing interneurons (PVINs) exert crucial inhibitory control over Aβ fiber-mediated nociceptive pathways at the spinal dorsal horn. The loss of their inhibitory tone leads to neuropathic symptoms, such as mechanical allodynia, via mechanisms that are not fully understood. This study identifies the reduced intrinsic excitability of PVINs as a potential cause for their decreased inhibitory output in nerve-injured condition. Combining computational and experimental approaches, we predict a calcium-dependent mechanism that modulates PVINs' electrical activity following nerve injury: a depletion of cytosolic calcium buffer allows for the rapid accumulation of intracellular calcium through the active membranes, which in turn potentiates SK channels and impedes spike generation. Our results therefore pinpoint SK channels as potential therapeutic targets for treating neuropathic symptoms.

A systematic review of computational models for the design of spinal cord stimulation therapies: from neural circuits to patient-specific simulations

Seventy years ago, Hodgkin and Huxley published the first mathematical model to describe action potential generation, laying the foundation for modern computational neuroscience. Since then, the field has evolved enormously, with studies spanning from basic neuroscience to clinical applications for neuromodulation. Computer models of neuromodulation have evolved in complexity and personalization, advancing clinical practice and novel neurostimulation therapies, such as spinal cord stimulation. Spinal cord stimulation is a therapy widely used to treat chronic pain, with rapidly expanding indications, such as restoring motor function. In general, simulations contributed dramatically to improve lead designs, stimulation configurations, waveform parameters and programming procedures and provided insight into potential mechanisms of action of electrical stimulation. Although the implementation of neural models are relentlessly increasing in number and complexity, it is reasonable to ask whether this observed increase in complexity is necessary for improved accuracy and, ultimately, for clinical efficacy. With this aim, we performed a systematic literature review and a qualitative meta-synthesis of the evolution of computational models, with a focus on complexity, personalization and the use of medical imaging to capture realistic anatomy. Our review showed that increased model complexity and personalization improved both mechanistic and translational studies. More specifically, the use of medical imaging enabled the development of patient-specific models that can help to transform clinical practice in spinal cord stimulation. Finally, we combined our results to provide clear guidelines for standardization and expansion of computational models for spinal cord stimulation.

Modeling homeostatic and circadian modulation of human pain sensitivity

Introduction

Mathematical modeling has played a significant role in understanding how homeostatic sleep pressure and the circadian rhythm interact to influence sleep-wake behavior. Pain sensitivity is also affected by these processes, and recent experimental results have measured the circadian and homeostatic components of the 24 h rhythm of thermal pain sensitivity in humans. To analyze how rhythms in pain sensitivity are affected by disruptions in sleep behavior and shifts in circadian rhythms, we introduce a dynamic mathematical model for circadian and homeostatic regulation of sleep-wake states and pain intensity.

Methods

The model consists of a biophysically based, sleep-wake regulation network model coupled to data-driven functions for the circadian and homeostatic modulation of pain sensitivity. This coupled sleep-wake-pain sensitivity model is validated by comparison to thermal pain intensities in adult humans measured across a 34 h sleep deprivation protocol.

Results

We use the model to predict dysregulation of pain sensitivity rhythms across different scenarios of sleep deprivation and circadian rhythm shifts, including entrainment to new environmental light and activity timing as occurs with jet lag and chronic sleep restriction. Model results show that increases in pain sensitivity occur under conditions of increased homeostatic sleep drive with nonlinear modulation by the circadian rhythm, leading to unexpected decreased pain sensitivity in some scenarios.

Discussion

This model provides a useful tool for pain management by predicting alterations in pain sensitivity due to varying or disrupted sleep schedules.

Degeneracy in epilepsy: multiple routes to hyperexcitable brain circuits and their repair

Due to its complex and multifaceted nature, developing effective treatments for epilepsy is still a major challenge. To deal with this complexity we introduce the concept of degeneracy to the field of epilepsy research: the ability of disparate elements to cause an analogous function or malfunction. Here, we review examples of epilepsy-related degeneracy at multiple levels of brain organisation, ranging from the cellular to the network and systems level. Based on these insights, we outline new multiscale and population modelling approaches to disentangle the complex web of interactions underlying epilepsy and to design personalised multitarget therapies.

Heterosynaptic long-term potentiation of non-nociceptive synapses requires endocannabinoids, NMDARs, CamKII, and PKCζ

Noxious stimuli or injury can trigger long-lasting sensitization to non-nociceptive stimuli (referred to as allodynia in mammals). Long-term potentiation (LTP) of nociceptive synapses has been shown to contribute to nociceptive sensitization (hyperalgesia) and there is even evidence of heterosynaptic spread of LTP contributing to this type of sensitization. This study will focus on how activation of nociceptors elicits heterosynaptic LTP (hetLTP) in non-nociceptive synapses. Previous studies in the medicinal leech (Hirudo verbana) have demonstrated that high-frequency stimulation (HFS) of nociceptors produces both homosynaptic LTP as well as hetLTP in non-nociceptive afferent synapses. This hetLTP involves endocannabinoid-mediated disinhibition of non-nociceptive synapses at the presynaptic level, but it is not clear if there are additional processes contributing to this synaptic potentiation. In this study, we found evidence for the involvement of postsynaptic level change and observed that postsynaptic N-methyl-d-aspartate (NMDA) receptors (NMDARs) were required for this potentiation. Next, Hirudo orthologs for known LTP signaling proteins, CamKII and PKCζ, were identified based on sequences from humans, mice, and the marine mollusk Aplysia. In electrophysiological experiments, inhibitors of CamKII (AIP) and PKCζ (ZIP) were found to interfere with hetLTP. Interestingly, CamKII was found to be necessary for both induction and maintenance of hetLTP, whereas PKCζ was only necessary for maintenance. These findings show that activation of nociceptors can elicit a potentiation of non-nociceptive synapses through a process that involves both endocannabinoid-mediated disinhibition and NMDAR-initiated signaling pathways.NEW & NOTEWORTHY Pain-related sensitization involves increases in signaling by non-nociceptive sensory neurons. This can allow non-nociceptive afferents to have access to nociceptive circuitry. In this study, we examine a form of synaptic potentiation in which nociceptor activity elicits increases in non-nociceptive synapses. This process involves endocannabinoids, "gating" the activation of NMDA receptors, which in turn activate CamKII and PKCζ. This study provides an important link in how nociceptive stimuli can enhance non-nociceptive signaling related to pain.

Measuring pain and nociception: Through the glasses of a computational scientist. Transdisciplinary overview of methods

In a healthy state, pain plays an important role in natural biofeedback loops and helps to detect and prevent potentially harmful stimuli and situations. However, pain can become chronic and as such a pathological condition, losing its informative and adaptive function. Efficient pain treatment remains a largely unmet clinical need. One promising route to improve the characterization of pain, and with that the potential for more effective pain therapies, is the integration of different data modalities through cutting edge computational methods. Using these methods, multiscale, complex, and network models of pain signaling can be created and utilized for the benefit of patients. Such models require collaborative work of experts from different research domains such as medicine, biology, physiology, psychology as well as mathematics and data science. Efficient work of collaborative teams requires developing of a common language and common level of understanding as a prerequisite. One of ways to meet this need is to provide easy to comprehend overviews of certain topics within the pain research domain. Here, we propose such an overview on the topic of pain assessment in humans for computational researchers. Quantifications related to pain are necessary for building computational models. However, as defined by the International Association of the Study of Pain (IASP), pain is a sensory and emotional experience and thus, it cannot be measured and quantified objectively. This results in a need for clear distinctions between nociception, pain and correlates of pain. Therefore, here we review methods to assess pain as a percept and nociception as a biological basis for this percept in humans, with the goal of creating a roadmap of modelling options.

Pareto optimality, economy-effectiveness trade-offs and ion channel degeneracy: improving population modelling for single neurons

Neurons encounter unavoidable evolutionary trade-offs between multiple tasks. They must consume as little energy as possible while effectively fulfilling their functions. Cells displaying the best performance for such multi-task trade-offs are said to be Pareto optimal, with their ion channel configurations underpinning their functionality. Ion channel degeneracy, however, implies that multiple ion channel configurations can lead to functionally similar behaviour. Therefore, instead of a single model, neuroscientists often use populations of models with distinct combinations of ionic conductances. This approach is called population (database or ensemble) modelling. It remains unclear, which ion channel parameters in the vast population of functional models are more likely to be found in the brain. Here we argue that Pareto optimality can serve as a guiding principle for addressing this issue by helping to identify the subpopulations of conductance-based models that perform best for the trade-off between economy and functionality. In this way, the high-dimensional parameter space of neuronal models might be reduced to geometrically simple low-dimensional manifolds, potentially explaining experimentally observed ion channel correlations. Conversely, Pareto inference might also help deduce neuronal functions from high-dimensional Patch-seq data. In summary, Pareto optimality is a promising framework for improving population modelling of neurons and their circuits.

Modernizing the NEURON Simulator for Sustainability, Portability, and Performance

The need for reproducible, credible, multiscale biological modeling has led to the development of standardized simulation platforms, such as the widely-used NEURON environment for computational neuroscience. Developing and maintaining NEURON over several decades has required attention to the competing needs of backwards compatibility, evolving computer architectures, the addition of new scales and physical processes, accessibility to new users, and efficiency and flexibility for specialists. In order to meet these challenges, we have now substantially modernized NEURON, providing continuous integration, an improved build system and release workflow, and better documentation. With the help of a new source-to-source compiler of the NMODL domain-specific language we have enhanced NEURON's ability to run efficiently, via the CoreNEURON simulation engine, on a variety of hardware platforms, including GPUs. Through the implementation of an optimized in-memory transfer mechanism this performance optimized backend is made easily accessible to users, providing training and model-development paths from laptop to workstation to supercomputer and cloud platform. Similarly, we have been able to accelerate NEURON's reaction-diffusion simulation performance through the use of just-in-time compilation. We show that these efforts have led to a growing developer base, a simpler and more robust software distribution, a wider range of supported computer architectures, a better integration of NEURON with other scientific workflows, and substantially improved performance for the simulation of biophysical and biochemical models.

Pain modulation in the spinal cord

The sensory inflow from the periphery that triggers innocuous and painful sensations is highly complex, capturing key elements of the nature of any stimulus, its location, intensity, and duration, and converting this to dynamic action potential firing across a wide population of afferents. While sensory afferents are highly specialized to detect these features, their input to the spinal cord also triggers active processing and modulation there which determines its output, to drive the sensory percept experienced and behavioral responses. Focus on such active spinal modulation was arguably first introduced by Melzack and Wall in their Spinal Cord Gate Control theory. This theory has had a profound influence on our understanding of pain, and especially its processing, as well as leading directly to the development of clinical interventions, and its historical importance certainly needs to be fully recognized. However, the enormous progress we are making in the understanding of the function of the somatosensory system, means that it is time to incorporate these newly discovered features into a more complex and accurate model of spinal sensory modulation.