State and location dependence of action potential metabolic cost in cortical pyramidal neurons

Aging-associated weakening of the action potential in fast-spiking interneurons in the human neocortex

Aging is associated with the slowdown of neuronal processing and cognitive performance in the brain; however, the exact cellular mechanisms behind this deterioration in humans are poorly elucidated. Recordings in human acute brain slices prepared from tissue resected during brain surgery enable the investigation of neuronal changes with age. Although neocortical fast-spiking cells are widely implicated in neuronal network activities underlying cognitive processes, they are vulnerable to neurodegeneration. Herein, we analyzed the electrical properties of 147 fast-spiking interneurons in neocortex samples resected in brain surgery from 106 patients aged 11-84 years. By studying the electrophysiological features of action potentials and passive membrane properties, we report that action potential overshoot significantly decreases and spike half-width increases with age. Moreover, the action potential maximum-rise speed (but not the repolarization speed or the afterhyperpolarization amplitude) significantly changed with age, suggesting a particular weakening of the sodium channel current generated in the soma. Cell passive membrane properties measured as the input resistance, membrane time constant, and cell capacitance remained unaffected by senescence. Thus, we conclude that the action potential in fast-spiking interneurons shows a significant weakening in the human neocortex with age. This may contribute to the deterioration of cortical functions by aging.

A Multimodal Fitting Approach to Construct Single-Neuron Models With Patch Clamp and High-Density Microelectrode Arrays

In computational neuroscience, multicompartment models are among the most biophysically realistic representations of single neurons. Constructing such models usually involves the use of the patch-clamp technique to record somatic voltage signals under different experimental conditions. The experimental data are then used to fit the many parameters of the model. While patching of the soma is currently the gold-standard approach to build multicompartment models, several studies have also evidenced a richness of dynamics in dendritic and axonal sections. Recording from the soma alone makes it hard to observe and correctly parameterize the activity of nonsomatic compartments. In order to provide a richer set of data as input to multicompartment models, we here investigate the combination of somatic patch-clamp recordings with recordings of high-density microelectrode arrays (HD-MEAs). HD-MEAs enable the observation of extracellular potentials and neural activity of neuronal compartments at subcellular resolution. In this work, we introduce a novel framework to combine patch-clamp and HD-MEA data to construct multicompartment models. We first validate our method on a ground-truth model with known parameters and show that the use of features extracted from extracellular signals, in addition to intracellular ones, yields models enabling better fits than using intracellular features alone. We also demonstrate our procedure using experimental data by constructing cell models from in vitro cell cultures. The proposed multimodal fitting procedure has the potential to augment the modeling efforts of the computational neuroscience community and provide the field with neuronal models that are more realistic and can be better validated.

Complex biophysical changes and reduced neuronal firing in an SCN8A variant associated with developmental delay and epilepsy

Mutations in the SCN8A gene, encoding the voltage-gated sodium channel NaV1.6, are associated with a range of neurodevelopmental syndromes. The p.(Gly1625Arg) (G1625R) mutation was identified in a patient diagnosed with developmental epileptic encephalopathy (DEE). While most of the characterized DEE-associated SCN8A mutations were shown to cause a gain-of-channel function, we show that the G1625R variant, positioned within the S4 segment of domain IV, results in complex effects. Voltage-clamp analyses of NaV1.6G1625R demonstrated a mixture of gain- and loss-of-function properties, including reduced current amplitudes, increased time constant of fast voltage-dependent inactivation, a depolarizing shift in the voltage dependence of activation and inactivation, and increased channel availability with high-frequency repeated depolarization. Current-clamp analyses in transfected cultured neurons revealed that these biophysical properties caused a marked reduction in the number of action potentials when firing was driven by the transfected mutant NaV1.6. Accordingly, computational modeling of mature cortical neurons demonstrated a mild decrease in neuronal firing when mimicking the patients' heterozygous SCN8A expression. Structural modeling of NaV1.6G1625R suggested the formation of a cation-π interaction between R1625 and F1588 within domain IV. Double-mutant cycle analysis revealed that this interaction affects the voltage dependence of inactivation in NaV1.6G1625R. Together, our studies demonstrate that the G1625R variant leads to a complex combination of gain and loss of function biophysical changes that result in an overall mild reduction in neuronal firing, related to the perturbed interaction network within the voltage sensor domain, necessitating personalized multi-tiered analysis for SCN8A mutations for optimal treatment selection.

Metabolic response of auditory brainstem neurons to their broad physiological activity range

Neurons exhibit a high energetic need, and the question arises as how they metabolically adapt to changing activity states. This is relevant for interpreting functional neuroimaging in different brain areas. Particularly, neurons with a broad firing range might exhibit metabolic adaptations. Therefore, we studied MNTB (medial nucleus of the trapezoid body) principal neurons, which generate action potentials (APs) at frequencies up to several hundred hertz. We performed the experiments in acute brainstem slices of the Mongolian gerbil (Meriones unguiculatus) at 22.5-24.5°C. Upon electrical stimulation of afferent MNTB fibres with 400 stimuli at varying frequencies, we monitored autofluorescence levels of NAD(P)H and FAD and determined the extremum amplitudes of their biphasic response. Additionally, we compared these data with alterations in O2 concentrations measured with an electrochemical sensor. These O2 changes are prominent since MNTB neurons rely on oxidative phosphorylation as shown by our pharmacological experiments. We calculated the O2 consumption rate as change in O2 concentration divided by stimulus durations, because these periods varied inversely with stimulus frequency as a result of the constant number of 400 stimuli applied. The O2 consumption rate increased with stimulation frequency up to a constant value at 600 Hz; that is, energy demand depends on temporal characteristics of activity despite the same number of stimuli. The rates showed no correlation with peak amplitudes of NAD(P)H or FAD, whilst the values of the two molecules were linearly correlated. This points at the complexity of analysing autofluorescence imaging for quantitative metabolic studies, because these values report only relative net changes of many superimposed oxidative and reductive processes. Monitoring O2 concentration rates is, thus, an important tool to improve the interpretation of NAD(P)H/FAD autofluorescence data, as they do not under all conditions and in all systems appropriately reflect the metabolic activity or energy demand.

A critical review of ethanol effects on neuronal firing: A metabolic perspective

Ethanol metabolism is relatively understudied in neurons, even though changes in neuronal metabolism are known to affect their activity. Recent work demonstrates that ethanol is preferentially metabolized over glucose as a source of carbon and energy, and it reprograms neurons to a state of reduced energy potential and diminished capacity to utilize glucose once ethanol is exhausted. Ethanol intake has been associated with changes in neuronal firing and specific brain activity (EEG) patterns have been linked with risk for alcohol use disorder (AUD). Furthermore, a haplotype of the inwardly rectifying potassium channel subunit, GIRK2, which plays a critical role in regulating excitability of neurons, has been linked with AUD and shown to be directly regulated by ethanol. At the same time, overexpression of GIRK2 prevents ethanol-induced metabolic changes. Based on the available evidence, we conclude that the mechanisms underlying the effects of ethanol on neuronal metabolism are a novel target for developing therapies for AUD.

Kinetics and functional consequences of BK channels activation by N-type Ca2+ channels in the dendrite of mouse neocortical layer-5 pyramidal neurons

The back-propagation of an action potential (AP) from the axon/soma to the dendrites plays a central role in dendritic integration. This process involves an intricate orchestration of various ion channels, but a comprehensive understanding of the contribution of each channel type remains elusive. In this study, we leverage ultrafast membrane potential recordings (Vm) and Ca2+ imaging techniques to shed light on the involvement of N-type voltage-gated Ca2+ channels (VGCCs) in layer-5 neocortical pyramidal neurons' apical dendrites. We found a selective interaction between N-type VGCCs and large-conductance Ca2+-activated K+ channels (BK CAKCs). Remarkably, we observe that BK CAKCs are activated within a mere 500 μs after the AP peak, preceding the peak of the Ca2+ current triggered by the AP. Consequently, when N-type VGCCs are inhibited, the early broadening of the AP shape amplifies the activity of other VGCCs, leading to an augmented total Ca2+ influx. A NEURON model, constructed to replicate and support these experimental results, reveals the critical coupling between N-type and BK channels. This study not only redefines the conventional role of N-type VGCCs as primarily involved in presynaptic neurotransmitter release but also establishes their distinct and essential function as activators of BK CAKCs in neuronal dendrites. Furthermore, our results provide original functional validation of a physical interaction between Ca2+ and K+ channels, elucidated through ultrafast kinetic reconstruction. This insight enhances our understanding of the intricate mechanisms governing neuronal signaling and may have far-reaching implications in the field.

Laminar organization of neocortical activities during systemic anoxia

The neocortex is highly susceptible to metabolic dysfunction. When exposed to global ischemia or anoxia, it suffers a slowly propagating wave of collective neuronal depolarization that ultimately impairs its structure and function. While the molecular signature of anoxic depolarization (AD) is well documented, little is known about the brain states that precede and follow AD onset. Here, by means of multisite extracellular local field potentials and intracellular recordings from identified pyramidal cells, we investigated the laminar expression of cortical activities induced by transient anoxia in rat primary somatosensory cortex. Soon after the interruption of brain oxygenation, we observed a well-organized sequence of stereotyped activity patterns across all cortical layers. This sequence included an initial period of beta-gamma activity, rapidly replaced by delta-theta oscillations followed by a decline in all spontaneous activites, marking the entry into a sustained period of electrical silence. Intracellular recordings revealed that cortical pyramidal neurons were depolarized and highly active during high-frequency activity, became inactive and devoid of synaptic potentials during the isoelectric state, and showed subthreshold composite synaptic depolarizations during the low-frequency period. Contrasting with the strong temporal coherence of pre-AD activities along the vertical axis of the cortical column, the onset of AD was not uniform across layers. AD initially occurred in layer 5 or 6 and then propagated bidirectionally in the upward and downward direction. Conversely, the post-anoxic waves that indicated the repolarization of cortical neurons upon brain reoxygenation did not exhibit a specific spatio-temporal profile. Pyramidal neurons from AD initiation site had a more depolarized resting potential and higher spontaneous firing rate compared to superficial cortical cells. We also found that the propagation pattern of AD was reliably reproduced by focal injection of an inhibitor of sodium‑potassium ATPases, suggesting that cortical AD dynamics could reflect layer-dependent variations in cellular metabolic regulations.

Human voltage-gated Na+ and K+ channel properties underlie sustained fast AP signaling

Human cortical pyramidal neurons are large, have extensive dendritic trees, and yet have unexpectedly fast input-output properties: Rapid subthreshold synaptic membrane potential changes are reliably encoded in timing of action potentials (APs). Here, we tested whether biophysical properties of voltage-gated sodium (Na+) and potassium (K+) currents in human pyramidal neurons can explain their fast input-output properties. Human Na+ and K+ currents exhibited more depolarized voltage dependence, slower inactivation, and faster recovery from inactivation compared with their mouse counterparts. Computational modeling showed that despite lower Na+ channel densities in human neurons, the biophysical properties of Na+ channels resulted in higher channel availability and contributed to fast AP kinetics stability. Last, human Na+ channel properties also resulted in a larger dynamic range for encoding of subthreshold membrane potential changes. Thus, biophysical adaptations of voltage-gated Na+ and K+ channels enable fast input-output properties of large human pyramidal neurons.

Establishing brain states in neuroimaging data

The definition of a brain state remains elusive, with varying interpretations across different sub-fields of neuroscience-from the level of wakefulness in anaesthesia, to activity of individual neurons, voltage in EEG, and blood flow in fMRI. This lack of consensus presents a significant challenge to the development of accurate models of neural dynamics. However, at the foundation of dynamical systems theory lies a definition of what constitutes the 'state' of a system-i.e., a specification of the system's future. Here, we propose to adopt this definition to establish brain states in neuroimaging timeseries by applying Dynamic Causal Modelling (DCM) to low-dimensional embedding of resting and task condition fMRI data. We find that ~90% of subjects in resting conditions are better described by first-order models, whereas ~55% of subjects in task conditions are better described by second-order models. Our work calls into question the status quo of using first-order equations almost exclusively within computational neuroscience and provides a new way of establishing brain states, as well as their associated phase space representations, in neuroimaging datasets.

Pathogenic SCN2A variants cause early-stage dysfunction in patient-derived neurons

Pathogenic heterozygous variants in SCN2A, which encodes the neuronal sodium channel NaV1.2, cause different types of epilepsy or intellectual disability (ID)/autism without seizures. Previous studies using mouse models or heterologous systems suggest that NaV1.2 channel gain-of-function typically causes epilepsy, whereas loss-of-function leads to ID/autism. How altered channel biophysics translate into patient neurons remains unknown. Here, we investigated iPSC-derived early-stage cortical neurons from ID patients harboring diverse pathogenic SCN2A variants [p.(Leu611Valfs*35); p.(Arg937Cys); p.(Trp1716*)] and compared them with neurons from an epileptic encephalopathy (EE) patient [p.(Glu1803Gly)] and controls. ID neurons consistently expressed lower NaV1.2 protein levels. In neurons with the frameshift variant, NaV1.2 mRNA and protein levels were reduced by ~ 50%, suggesting nonsense-mediated decay and haploinsufficiency. In other ID neurons, only protein levels were reduced implying NaV1.2 instability. Electrophysiological analysis revealed decreased sodium current density and impaired action potential (AP) firing in ID neurons, consistent with reduced NaV1.2 levels. In contrast, epilepsy neurons displayed no change in NaV1.2 levels or sodium current density, but impaired sodium channel inactivation. Single-cell transcriptomics identified dysregulation of distinct molecular pathways including inhibition of oxidative phosphorylation in neurons with SCN2A haploinsufficiency and activation of calcium signaling and neurotransmission in epilepsy neurons. Together, our patient iPSC-derived neurons reveal characteristic sodium channel dysfunction consistent with biophysical changes previously observed in heterologous systems. Additionally, our model links the channel dysfunction in ID to reduced NaV1.2 levels and uncovers impaired AP firing in early-stage neurons. The altered molecular pathways may reflect a homeostatic response to NaV1.2 dysfunction and can guide further investigations.

Evidence of Metabolic Dysfunction in Amyotrophic Lateral Sclerosis (ALS) Patients and Animal Models

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease that affects motor neurons, leading to muscle weakness, paralysis, and eventual death. Research from the past few decades has appreciated that ALS is not only a disease of the motor neurons but also a disease that involves systemic metabolic dysfunction. This review will examine the foundational research of understanding metabolic dysfunction in ALS and provide an overview of past and current studies in ALS patients and animal models, spanning from full systems to various metabolic organs. While ALS-affected muscle tissue exhibits elevated energy demand and a fuel preference switch from glycolysis to fatty acid oxidation, adipose tissue in ALS undergoes increased lipolysis. Dysfunctions in the liver and pancreas contribute to impaired glucose homeostasis and insulin secretion. The central nervous system (CNS) displays abnormal glucose regulation, mitochondrial dysfunction, and increased oxidative stress. Importantly, the hypothalamus, a brain region that controls whole-body metabolism, undergoes atrophy associated with pathological aggregates of TDP-43. This review will also cover past and present treatment options that target metabolic dysfunction in ALS and provide insights into the future of metabolism research in ALS.

Somatic and terminal CB1 receptors are differentially coupled to voltage-gated sodium channels in neocortical neurons

Endogenous cannabinoid signaling is vital for important brain functions, and the same pathways can be modified pharmacologically to treat pain, epilepsy, and posttraumatic stress disorder. Endocannabinoid-mediated changes to excitability are predominantly attributed to 2-arachidonoylglycerol (2-AG) acting presynaptically via the canonical cannabinoid receptor, CB1. Here, we identify a mechanism in the neocortex by which anandamide (AEA), another major endocannabinoid, but not 2-AG, powerfully inhibits somatically recorded voltage-gated sodium channel (VGSC) currents in the majority of neurons. This pathway involves intracellular CB1 that, when activated by anandamide, decreases the likelihood of recurrent action potential generation. WIN 55,212-2 similarly activates CB1 and inhibits VGSC currents, indicating that this pathway is also positioned to mediate the actions of exogenous cannabinoids on neuronal excitability. The coupling between CB1 and VGSCs is absent at nerve terminals, and 2-AG does not block somatic VGSC currents, indicating functional compartmentalization of the actions of two endocannabinoids.

Integrated transcriptomic and metabolomic profiling reveals dysregulation of purine metabolism during the acute phase of spinal cord injury in rats

Introduction

Spinal cord injury (SCI) results in drastic dysregulation of microenvironmental metabolism during the acute phase, which greatly affects neural recovery. A better insight into the potential molecular pathways of metabolic dysregulation by multi-omics analysis could help to reveal targets that promote nerve repair and regeneration in the future.

Materials and methods

We established the SCI model and rats were randomly divided into two groups: the acute-phase SCI (ASCI) group (n = 14, 3 days post-SCI) and the sham group with day-matched periods (n = 14, without SCI). In each group, rats were sacrificed at 3 days post-surgery for histology study (n = 3), metabolome sequencing (n = 5), transcriptome sequencing (n = 3), and quantitative real-time polymerase chain reaction (n = 3). The motor function of rats was evaluated by double-blind Basso, Beattie, and Bresnahan (BBB) Locomotor Scores at 0, 1, 2, 3 days post-SCI in an open field area. Then the transcriptomic and metabolomic data were integrated in SCI model of rat to reveal the underlying molecular pathways of microenvironmental metabolic dysregulation.

Results

The histology of the microenvironment was significantly altered in ASCI and the locomotor function was significantly reduced in rats. Metabolomics analysis showed that 360 metabolites were highly altered during the acute phase of SCI, of which 310 were up-regulated and 50 were down-regulated, and bioinformatics analysis revealed that these differential metabolites were mainly enriched in arginine and proline metabolism, D-glutamine and D-glutamate metabolism, purine metabolism, biosynthesis of unsaturated fatty acids. Transcriptomics results showed that 5,963 genes were clearly altered, of which 2,848 genes were up-regulated and 3,115 genes were down-regulated, and these differentially expressed genes were mainly involved in response to stimulus, metabolic process, immune system process. Surprisingly, the Integrative analysis revealed significant dysregulation of purine metabolism at both transcriptome and metabolome levels in the acute phase of SCI, with 48 differential genes and 16 differential metabolites involved. Further analysis indicated that dysregulation of purine metabolism could seriously affect the energy metabolism of the injured microenvironment and increase oxidative stress as well as other responses detrimental to nerve repair and regeneration.

Discussion

On the whole, we have for the first time combined transcriptomics and metabolomics to systematically analyze the potential molecular pathways of metabolic dysregulation in the acute phase of SCI, which will contribute to broaden our understanding of the sophisticated molecular mechanisms of SCI, in parallel with serving as a foundation for future studies of neural repair and regeneration after SCI.

Ultrafast simulation of large-scale neocortical microcircuitry with biophysically realistic neurons

Understanding the activity of the mammalian brain requires an integrative knowledge of circuits at distinct scales, ranging from ion channel gating to circuit connectomics. Computational models are regularly employed to understand how multiple parameters contribute synergistically to circuit behavior. However, traditional models of anatomically and biophysically realistic neurons are computationally demanding, especially when scaled to model local circuits. To overcome this limitation, we trained several artificial neural network (ANN) architectures to model the activity of realistic multicompartmental cortical neurons. We identified an ANN architecture that accurately predicted subthreshold activity and action potential firing. The ANN could correctly generalize to previously unobserved synaptic input, including in models containing nonlinear dendritic properties. When scaled, processing times were orders of magnitude faster compared with traditional approaches, allowing for rapid parameter-space mapping in a circuit model of Rett syndrome. Thus, we present a novel ANN approach allowing for rapid, detailed network experiments using inexpensive and commonly available computational resources.

Tracking axon initial segment plasticity using high-density microelectrode arrays: A computational study

Despite being composed of highly plastic neurons with extensive positive feedback, the nervous system maintains stable overall function. To keep activity within bounds, it relies on a set of negative feedback mechanisms that can induce stabilizing adjustments and that are collectively termed “homeostatic plasticity.” Recently, a highly excitable microdomain, located at the proximal end of the axon—the axon initial segment (AIS)—was found to exhibit structural modifications in response to activity perturbations. Though AIS plasticity appears to serve a homeostatic purpose, many aspects governing its expression and its functional role in regulating neuronal excitability remain elusive. A central challenge in studying the phenomenon is the rich heterogeneity of its expression (distal/proximal relocation, shortening, lengthening) and the variability of its functional role. A potential solution is to track AISs of a large number of neurons over time and attempt to induce structural plasticity in them. To this end, a promising approach is to use extracellular electrophysiological readouts to track a large number of neurons at high spatiotemporal resolution by means of high-density microelectrode arrays (HD-MEAs). However, an analysis framework that reliably identifies specific activity signatures that uniquely map on to underlying microstructural changes is missing. In this study, we assessed the feasibility of such a task and used the distal relocation of the AIS as an exemplary problem. We used sophisticated computational models to systematically explore the relationship between incremental changes in AIS positions and the specific consequences observed in simulated extracellular field potentials. An ensemble of feature changes in the extracellular fields that reliably characterize AIS plasticity was identified. We trained models that could detect these signatures with remarkable accuracy. Based on these findings, we propose a hybrid analysis framework that could potentially enable high-throughput experimental studies of activity-dependent AIS plasticity using HD-MEAs.

Associations of Human Milk Oligosaccharides with Infant Brain Tissue Organization and Regional Blood Flow at 1 Month of Age

Animal studies have shown that human milk oligosaccharides (HMOs) are important in early brain development, yet their roles have not been assessed in humans. The purpose of this study was to determine the associations of HMOs with MRI indices of tissue microstructure and regional cerebral blood flow (rCBF) in infants. Mother–infant pairs (N = 20) were recruited at 1 month postpartum. Milk was assayed for the concentrations of the HMOs 2′-fucosyllactose (2′FL), 3-fucosyllactose (3FL), 3′-sialyllactose (3′SL), and 6′-sialyllactose (6′SL). Diffusion and arterial spin labeling measures were acquired using a 3.0-Tesla MRI scanner. Multiple linear regression was used to assess the voxel-wise associations of HMOs with fractional anisotropy (FA), mean diffusivity (MD), and rCBF values across the brain. After adjusting for pre-pregnancy BMI, sex, birthweight, and postmenstrual age at time of scan, a higher 2′FL concentration was associated with reduced FA, increased MD, and reduced rCBF in similar locations within the cortical mantle. Higher 3FL and 3′SL concentrations were associated with increased FA, reduced MD, and increased rCBF in similar regions within the developing white matter. The concentration of 6′SL was not associated with MRI indices. Our data reveal that fucosylated and sialylated HMOs differentially associate with indices of tissue microstructure and rCBF, suggesting specific roles for 2′FL, 3FL, and 3′SL in early brain maturation.

Using tissue microstructure and multimodal MRI to parse the phenotypic heterogeneity and cellular basis of autism spectrum disorder

BACKGROUND

Identifying the brain bases for phenotypic heterogeneity in Autism Spectrum Disorder (ASD) will advance understanding of its pathogenesis and improve its clinical management.

METHODS

We compared Diffusion Tensor Imaging (DTI) indices and connectome measures between 77 ASD and 88 Typically Developing (TD) control participants. We also assessed voxel-wise associations of DTI indices with measures of regional cerebral blood flow (rCBF) and N-acetylaspartate (NAA) to understand how tissue microstructure associates with cellular metabolism and neuronal density, respectively.

RESULTS

Autism Spectrum Disorder participants had significantly lower fractional anisotropy (FA) and higher diffusivity values in deep white matter tracts, likely representing ether reduced myelination by oligodendrocytes or a reduced density of myelinated axons. Greater abnormalities in these measures and regions were associated with higher ASD symptom scores. Participant age, sex and IQ significantly moderated these group differences. Path analyses showed that reduced NAA levels accounted significantly for higher diffusivity and higher rCBF values in ASD compared with TD participants.

CONCLUSIONS

Reduced neuronal density (reduced NAA) likely underlies abnormalities in DTI indices of white matter microstructure in ASD, which in turn are major determinants of elevated blood flow. Together, these findings suggest the presence of reduced axonal density and axonal pathology in ASD white matter. Greater pathology in turn accounts for more severe symptoms, lower intellectual ability, and reduced global efficiency for measures of white matter connectivity in ASD.

An automated method for precise axon reconstruction from recordings of high-density micro-electrode arrays

Objective

Neurons communicate with each other by sending action potentials through their axons. The velocity of axonal signal propagation describes how fast electrical action potentials can travel. This velocity can be affected in a human brain by several pathologies, including multiple sclerosis, traumatic brain injury and channelopathies. High-density microelectrode arrays (HD-MEAs) provide unprecedented spatio-temporal resolution to extracellularly record neural electrical activity. The high density of the recording electrodes enables to image the activity of individual neurons down to subcellular resolution, which includes the propagation of axonal signals. However, axon reconstruction, to date, mainly relies on manual approaches to select the electrodes and channels that seemingly record the signals along a specific axon, while an automated approach to track multiple axonal branches in extracellular action-potential recordings is still missing.

Approach

In this article, we propose a fully automated approach to reconstruct axons from extracellular electrical-potential landscapes, so-called “electrical footprints” of neurons. After an initial electrode and channel selection, the proposed method first constructs a graph based on the voltage signal amplitudes and latencies. Then, the graph is interrogated to extract possible axonal branches. Finally, the axonal branches are pruned, and axonal action-potential propagation velocities are computed.

Main results

We first validate our method using simulated data from detailed reconstructions of neurons, showing that our approach is capable of accurately reconstructing axonal branches. We then apply the reconstruction algorithm to experimental recordings of HD-MEAs and show that it can be used to determine axonal morphologies and signal-propagation velocities at high throughput.

Significance

We introduce a fully automated method to reconstruct axonal branches and estimate axonal action-potential propagation velocities using HD-MEA recordings. Our method yields highly reliable and reproducible velocity estimations, which constitute an important electrophysiological feature of neuronal preparations.

Therapeutic Induction of Energy Metabolism Reduces Neural Tissue Damage and Increases Microglia Activation in Severe Spinal Cord Injury

Neural tissue has high metabolic requirements. Following spinal cord injury (SCI), the damaged, tissue suffers from a severe metabolic impairment, which aggravates axonal degeneration and, neuronal loss. Impaired cellular energetic, tricarboxylic acid (TCA) cycle and oxidative, phosphorylation metabolism in neuronal cells has been demonstrated to be a major cause of neural tissue death and regeneration failure following SCI. Therefore, rewiring the spinal cord cell metabolism may be an innovative therapeutic strategy for the treatment of SCI. In this study, we evaluated the therapeutic effect of the recovery of oxidative metabolism in a mouse model of severe contusive SCI. Oral administration of TCA cycle intermediates, co-factors, essential amino acids, and branched-chain amino acids was started 3 days post-injury and continued until the end of the experimental procedures. Metabolomic, immunohistological, and biochemical analyses were performed on the injured spinal cord sections. Administration of metabolic precursors enhanced spinal cord oxidative metabolism. In line with this metabolic shift, we observed the activation of the mTORC1 anabolic pathway, the increase in mitochondrial mass, and ROS defense which effectively prevented the injury-induced neural cell apoptosis in treated animals. Consistently, we found more choline acetyltransferase (ChAT)-expressing motor neurons and increased neurofilament positive corticospinal axons in the spinal cord parenchyma of the treated mice. Interestingly, oral administration of the metabolic precursors increased the number of activated microglia expressing the CD206 marker suggestive of a, pro-resolutive, M2-like phenotype. These molecular and histological modifications observed in treated animals ultimately led to a significant, although partial, improvement of the motor functions. Our data demonstrate that rewiring the cellular metabolism can represent an effective strategy to treat SCI.

NeuroGPU: Accelerating multi-compartment, biophysically detailed neuron simulations on GPUs

BACKGROUND

The membrane potential of individual neurons depends on a large number of interacting biophysical processes operating on spatial-temporal scales spanning several orders of magnitude. The multi-scale nature of these processes dictates that accurate prediction of membrane potentials in specific neurons requires the utilization of detailed simulations. Unfortunately, constraining parameters within biologically detailed neuron models can be difficult, leading to poor model fits. This obstacle can be overcome partially by numerical optimization or detailed exploration of parameter space. However, these processes, which currently rely on central processing unit (CPU) computation, often incur orders of magnitude increases in computing time for marginal improvements in model behavior. As a result, model quality is often compromised to accommodate compute resources.

NEW METHOD

Here, we present a simulation environment, NeuroGPU, that takes advantage of the inherent parallelized structure of the graphics processing unit (GPU) to accelerate neuronal simulation.

RESULTS & COMPARISON WITH EXISTING METHODS

NeuroGPU can simulate most biologically detailed models 10-200 times faster than NEURON simulation running on a single core and 5 times faster than GPU simulators (CoreNEURON). NeuroGPU is designed for model parameter tuning and best performs when the GPU is fully utilized by running multiple (> 100) instances of the same model with different parameters. When using multiple GPUs, NeuroGPU can reach to a speed-up of 800 fold compared to single core simulations, especially when simulating the same model morphology with different parameters. We demonstrate the power of NeuoGPU through large-scale parameter exploration to reveal the response landscape of a neuron. Finally, we accelerate numerical optimization of biophysically detailed neuron models to achieve highly accurate fitting of models to simulation and experimental data.

CONCLUSIONS

Thus, NeuroGPU is the fastest available platform that enables rapid simulation of multi-compartment, biophysically detailed neuron models on commonly used computing systems accessible by many scientists.

A multimodal 3D neuro-microphysiological system with neurite-trapping microelectrodes

Three-dimensional cell technologies as pre-clinical models are emerging tools for mimicking the structural and functional complexity of the nervous system. The accurate exploration of phenotypes in engineered 3D neuronal cultures, however, demands morphological, molecular and especially functional measurements. Particularly crucial is measurement of electrical activity of individual neurons with millisecond resolution. Current techniques rely on customized electrophysiological recording set-ups, characterized by limited throughput and poor integration with other readout modalities. Here we describe a novel approach, using multiwell glass microfluidic microelectrode arrays, allowing non-invasive electrical recording from engineered 3D neural tissues. We demonstrate parallelized studies with reference compounds, calcium imaging and optogenetic stimulation. Additionally, we show how microplate compatibility allows automated handling and high-content analysis of human induced pluripotent stem cell-derived neurons. This microphysiological platform opens up new avenues for high-throughput studies on the functional, morphological and molecular details of neurological diseases and their potential treatment by therapeutic compounds.

Synaptic vesicle pools are a major hidden resting metabolic burden of nerve terminals

The brain is a metabolically fragile organ as compromises in fuel availability rapidly degrade cognitive function. Nerve terminals are likely loci of this vulnerability as they do not store sufficient ATP molecules, needing to synthesize them during activity or suffer acute degradation in performance. The ability of on-demand ATP synthesis to satisfy activity-driven ATP hydrolysis will depend additionally on the magnitude of local resting metabolic processes. We show here that synaptic vesicle (SV) pools are a major source of presynaptic basal energy consumption. This basal metabolic processes arises from SV-resident V-ATPases compensating for a hidden resting H+ efflux from the SV lumen. We show that this steady-state H+ efflux (i) is mediated by vesicular neurotransmitter transporters, (ii) is independent of the SV cycle, (iii) accounts for up to 44% of the resting synaptic energy consumption, and (iv) contributes substantially to nerve terminal intolerance of fuel deprivation.

Paradoxical hyperexcitability from NaV1.2 sodium channel loss in neocortical pyramidal cells

Loss-of-function variants in the gene SCN2A, which encodes the sodium channel NaV1.2, are strongly associated with autism spectrum disorder and intellectual disability. An estimated 20%-30% of children with these variants also suffer from epilepsy, with altered neuronal activity originating in neocortex, a region where NaV1.2 channels are expressed predominantly in excitatory pyramidal cells. This is paradoxical, as sodium channel loss in excitatory cells would be expected to dampen neocortical activity rather than promote seizure. Here, we examined pyramidal neurons lacking NaV1.2 channels and found that they were intrinsically hyperexcitable, firing high-frequency bursts of action potentials (APs) despite decrements in AP size and speed. Compartmental modeling and dynamic-clamp recordings revealed that NaV1.2 loss prevented potassium channels from properly repolarizing neurons between APs, increasing overall excitability by allowing neurons to reach threshold for subsequent APs more rapidly. This cell-intrinsic mechanism may, therefore, account for why SCN2A loss-of-function can paradoxically promote seizure.

Active dendrites enable strong but sparse inputs to determine orientation selectivity

The dendrites of neocortical pyramidal neurons are excitable. However, it is unknown how synaptic inputs engage nonlinear dendritic mechanisms during sensory processing in vivo, and how they in turn influence action potential output. Here, we provide a quantitative account of the relationship between synaptic inputs, nonlinear dendritic events, and action potential output. We developed a detailed pyramidal neuron model constrained by in vivo dendritic recordings. We drive this model with realistic input patterns constrained by sensory responses measured in vivo and connectivity measured in vitro. We show mechanistically that under realistic conditions, dendritic Na+ and NMDA spikes are the major determinants of neuronal output in vivo. We demonstrate that these dendritic spikes can be triggered by a surprisingly small number of strong synaptic inputs, in some cases even by single synapses. We predict that dendritic excitability allows the 1% strongest synaptic inputs of a neuron to control the tuning of its output. Active dendrites therefore allow smaller subcircuits consisting of only a few strongly connected neurons to achieve selectivity for specific sensory features.

Optimising the energetic cost of the glutamatergic synapse

As for electronic computation, neural information processing is energetically expensive. This is because information is coded in the brain as membrane voltage changes, which are generated largely by passive ion movements down electrochemical gradients, and these ion movements later need to be reversed by active ATP-dependent ion pumping. This article will review how much of the energetic cost of the brain reflects the activity of glutamatergic synapses, consider the relative amount of energy used pre- and postsynaptically, outline how evolution has energetically optimised synapse function by adjusting the presynaptic release probability and the postsynaptic number of glutamate receptors, and speculate on how energy use by synapses may be sensed and adjusted.

Subcellular Distribution of Persistent Sodium Conductance in Cortical Pyramidal Neurons

Cortical pyramidal neurons possess a persistent Na+ current (I NaP), which, in contrast to the larger transient current, does not undergo rapid inactivation. Although relatively quite small, I NaP is active at subthreshold voltages and therefore plays an important role in neuronal input-output processing. The subcellular distribution of channels responsible for I NaP and the mechanisms that render them persistent are not known. Using high-speed fluorescence Na+ imaging and whole-cell recordings in brain slices obtained from mice of either sex, we reconstructed the I NaP elicited by slow voltage ramps in soma and processes of cortical pyramidal neurons. We found that in all neuronal compartments, the relationship between persistent Na+ conductance and membrane voltage has the shape of a Boltzmann function. Although the density of channels underlying I NaP was about twofold lower in the axon initial segment (AIS) than in the soma, the axonal channels were activated by ∼10 mV less depolarization than were somatic channels. This difference in voltage dependence explains why, at functionally critical subthreshold voltages, most I NaP originates in the AIS. Finally, we show that endogenous polyamines constrain I NaP availability in both somatodendritic and axonal compartments of nondialyzed cortical neurons.SIGNIFICANCE STATEMENT The most salient characteristic of neuronal sodium channels is fast inactivation. However, a fraction of the sodium current does not inactivate. In cortical neurons, persistent current (I NaP) plays a prominent role in many important functions. Its subcellular distribution and generation mechanisms are, however, elusive. Using high-speed fluorescence Na+ imaging and electrical recordings, we reconstructed the I NaP in soma and processes of cortical pyramidal neurons. We found that at near-threshold voltages I NaP originates predominately from the axon, because of the distinctive voltage dependence of the underlying channels and not because of their high density. Finally, we show that the presence of endogenous polyamines significantly constrains I NaP availability in all compartments of nondialyzed cortical neurons.

Small Size of Recorded Neuronal Structures Confines the Accuracy in Direct Axonal Voltage Measurements

Patch-clamp instruments including amplifier circuits and pipettes affect the recorded voltage signals. We hypothesized that realistic and complete in silico representation of recording instruments together with detailed morphology and biophysics of small recorded structures will reveal signal distortions and provide a tool that predicts native, instrument-free electrical signals from distorted voltage recordings. Therefore, we built a model that was verified by small axonal recordings. The model accurately recreated actual action potential (AP) measurements with typical recording artefacts and predicted the native electrical behavior. The simulations verified that recording instruments substantially filter voltage recordings. Moreover, we revealed that instrumentation directly interferes with local signal generation depending on the size of the recorded structures, which complicates the interpretation of recordings from smaller structures, such as axons. However, our model offers a straightforward approach that predicts the native waveforms of fast voltage signals and the underlying conductances even from the smallest neuronal structures.

Frequency-Dependent Energy Demand of Dendritic Responses to Deep Brain Stimulation in Thalamic Neurons: A Model-Based Study

Thalamic deep brain stimulation (DBS) generates excitatory postsynaptic currents and action potentials (APs) by triggering large numbers of synaptic inputs to local cells, which also activates axonal spikes to antidromically invade the soma and dendrites. To maintain signaling, the evoked dendritic responses require metabolic energy to restore ion gradients in each dendrite. The objective of this study is to estimate the energy demand associated with dendritic responses to thalamic DBS. We use a morphologically realistic computational model to simulate dendritic activity in thalamocortical (TC) relay neurons with axonal intracellular stimulation or DBS-like extracellular stimulation. We determine the metabolic cost by calculating the number of adenosine triphosphate (ATP) expended to pump Na⁺ and Ca²⁺ ions out of each dendrite. The ATP demand of dendritic activity exhibits frequency dependence, which is determined by the number of spikes in the dendrites. Each backpropagating AP from the soma activates a spike in the dendrites, and the dendritic firing is dominated by antidromic activation of the soma. High stimulus frequencies decrease dendritic ATP cost by reducing the fidelity of antidromic activation. Synaptic inputs and stimulus-induced polarization govern the ATP cost of dendritic responses by facilitating/suppressing antidromic activation, which also influences the ATP cost by depolarizing/hyperpolarizing each dendrite. These findings are important for understanding the synaptic signaling energy in TC relay neurons and metabolism-dependent functional imaging data of thalamic DBS.

Human larynx motor cortices coordinate respiration for vocal-motor control

Vocal flexibility is a hallmark of the human species, most particularly the capacity to speak and sing. This ability is supported in part by the evolution of a direct neural pathway linking the motor cortex to the brainstem nucleus that controls the larynx the primary sound source for communication. Early brain imaging studies demonstrated that larynx motor cortex at the dorsal end of the orofacial division of motor cortex (dLMC) integrated laryngeal and respiratory control, thereby coordinating two major muscular systems that are necessary for vocalization. Neurosurgical studies have since demonstrated the existence of a second larynx motor area at the ventral extent of the orofacial motor division (vLMC) of motor cortex. The vLMC has been presumed to be less relevant to speech motor control, but its functional role remains unknown. We employed a novel ultra-high field (7T) magnetic resonance imaging paradigm that combined singing and whistling simple melodies to localise the larynx motor cortices and test their involvement in respiratory motor control. Surprisingly, whistling activated both 'larynx areas' more strongly than singing despite the reduced involvement of the larynx during whistling. We provide further evidence for the existence of two larynx motor areas in the human brain, and the first evidence that laryngeal-respiratory integration is a shared property of both larynx motor areas. We outline explicit predictions about the descending motor pathways that give these cortical areas access to both the laryngeal and respiratory systems and discuss the implications for the evolution of speech.

Multiscale representation of very large environments in the hippocampus of flying bats

Hippocampal place cells encode the animal's location. Place cells were traditionally studied in small environments, and nothing is known about large ethologically relevant spatial scales. We wirelessly recorded from hippocampal dorsal CA1 neurons of wild-born bats flying in a long tunnel (200 meters). The size of place fields ranged from 0.6 to 32 meters. Individual place cells exhibited multiple fields and a multiscale representation: Place fields of the same neuron differed up to 20-fold in size. This multiscale coding was observed from the first day of exposure to the environment, and also in laboratory-born bats that never experienced large environments. Theoretical decoding analysis showed that the multiscale code allows representation of very large environments with much higher precision than that of other codes. Together, by increasing the spatial scale, we discovered a neural code that is radically different from classical place codes.

Alternative splicing potentiates dysfunction of early-onset epileptic encephalopathy SCN2A variants

Epileptic encephalopathies are severe forms of infantile-onset epilepsy often complicated by severe neurodevelopmental impairments. Some forms of early-onset epileptic encephalopathy (EOEE) have been associated with variants in SCN2A, which encodes the brain voltage-gated sodium channel NaV1.2. Many voltage-gated sodium channel genes, including SCN2A, undergo developmentally regulated mRNA splicing. The early onset of these disorders suggests that developmentally regulated alternative splicing of NaV1.2 may be an important consideration when elucidating the pathophysiological consequences of epilepsy-associated variants. We hypothesized that EOEE-associated NaV1.2 variants would exhibit greater dysfunction in a splice isoform that is prominently expressed during early development. We engineered five EOEE-associated NaV1.2 variants (T236S, E999K, S1336Y, T1623N, and R1882Q) into the adult and neonatal splice isoforms of NaV1.2 and performed whole-cell voltage clamp to elucidate their functional properties. All variants exhibited functional defects that could enhance neuronal excitability. Three of the five variants (T236S, E999K, and S1336Y) exhibited greater dysfunction in the neonatal isoform compared with those observed in the adult isoform. Computational modeling of a developing cortical pyramidal neuron indicated that T236S, E999K, S1336Y, and R1882Q showed hyperexcitability preferentially in immature neurons. These results suggest that both splice isoform and neuronal developmental stage influence how EOEE-associated NaV1.2 variants affect neuronal excitability.

Communication consumes 35 times more energy than computation in the human cortex, but both costs are needed to predict synapse number

Darwinian evolution tends to produce energy-efficient outcomes. On the other hand, energy limits computation, be it neural and probabilistic or digital and logical. Taking a particular energy-efficient viewpoint, we define neural computation and make use of an energy-constrained computational function. This function can be optimized over a variable that is proportional to the number of synapses per neuron. This function also implies a specific distinction between adenosine triphosphate (ATP)-consuming processes, especially computation per se vs. the communication processes of action potentials and transmitter release. Thus, to apply this mathematical function requires an energy audit with a particular partitioning of energy consumption that differs from earlier work. The audit points out that, rather than the oft-quoted 20 W of glucose available to the human brain, the fraction partitioned to cortical computation is only 0.1 W of ATP [L. Sokoloff, Handb. Physiol. Sect. I Neurophysiol. 3, 1843-1864 (1960)] and [J. Sawada, D. S. Modha, "Synapse: Scalable energy-efficient neurosynaptic computing" in Application of Concurrency to System Design (ACSD) (2013), pp. 14-15]. On the other hand, long-distance communication costs are 35-fold greater, 3.5 W. Other findings include 1) a [Formula: see text]-fold discrepancy between biological and lowest possible values of a neuron's computational efficiency and 2) two predictions of N, the number of synaptic transmissions needed to fire a neuron (2,500 vs. 2,000).

A Metabolic Landscape for Maintaining Retina Integrity and Function

Neurons have high metabolic demands that are almost exclusively met by glucose supplied from the bloodstream. Glucose is utilized in complex metabolic interactions between neurons and glia cells, described by the astrocyte-neuron lactate shuttle (ANLS) hypothesis. The neural retina faces similar energy demands to the rest of the brain, with additional high anabolic needs to support continuous renewal of photoreceptor outer segments. This demand is met by a fascinating variation of the ANLS in which photoreceptors are the central part of a metabolic landscape, using glucose and supplying surrounding cells with metabolic intermediates. In this review we summarize recent evidence on how neurons, in particular photoreceptors, meet their energy and biosynthetic requirements by comprising a metabolic landscape of interdependent cells.

Memory retention in pyramidal neurons: a unified model of energy-based homo and heterosynaptic plasticity with homeostasis

The brain can learn new tasks without forgetting old ones. This memory retention is closely associated with the long-term stability of synaptic strength. To understand the capacity of pyramidal neurons to preserve memory under different tasks, we established a plasticity model based on the postsynaptic membrane energy state, in which the change in synaptic strength depends on the difference between the energy state after stimulation and the resting energy state. If the post-stimulation energy state is higher than the resting energy state, then synaptic depression occurs. On the contrary, the synapse is strengthened. Our model unifies homo- and heterosynaptic plasticity and can reproduce synaptic plasticity observed in multiple experiments, such as spike-timing-dependent plasticity, and cooperative plasticity with few and common parameters. Based on the proposed plasticity model, we conducted a simulation study on how the activation patterns of dendritic branches by different tasks affect the synaptic connection strength of pyramidal neurons. We further investigate the formation mechanism by which different tasks activate different dendritic branches. Simulation results show that compare to the classic plasticity model, the plasticity model we proposed can achieve a better spatial separation of different branches activated by different tasks in pyramidal neurons, which deepens our insight into the memory retention mechanism of brains.

Optical measurement of physiological sodium currents in the axon initial segment

KEY POINTS

Τhe axonal Na⁺ fluorescence underlying an action potential in the axon initial segment was optically measured at unprecedented temporal resolution. The measurement allowed resolution of the kinetics of the Na⁺ current at different axonal locations. The distinct components of the Na⁺ current were correlated with the kinetics of the action potential. NEURON simulations from a modified published model qualitatively predicted the experimentally measured Na⁺ current. The present method permits the direct investigation of the kinetic behaviour of native Na⁺ channels under physiological and pathological conditions.

ABSTRACT

In most neurons of the mammalian central nervous system, the action potential (AP) is generated in the axon initial segment (AIS) by a fast Na⁺ current mediated by voltage-gated Na⁺ channels. While the axonal Na⁺ signal associated with the AP has been measured using fluorescent Na⁺ indicators, the insufficient resolution of these recordings has not allowed tracking the Na⁺ current kinetics underlying this fundamental event. In this article, we report the first optical measurement of Na⁺ currents in the AIS of pyramidal neurons of layer 5 of the somatosensory cortex from brain slices of the mouse. This measurement was obtained by achieving a temporal resolution of 100 μs in the Na⁺ imaging technique, with a pixel resolution of 0.5 μm, and by calculating the time-derivative of the Na⁺ change corrected for longitudinal diffusion. We identified a subthreshold current before the AP, a fast-inactivating current peaking during the rise of the AP and a non-inactivating current during the AP repolarization. We established a correlation between the kinetics of the non-inactivating current at different distances from the soma and the kinetics of the somatic AP. We quantitatively compared the experimentally measured Na⁺ current with the current obtained by computer simulation of published NEURON models, demonstrating how the present approach can lead to the correct estimate of the native behaviour of Na⁺ channels. Finally, we discuss how the present approach can be used to investigate the physiological or pathological function of different channel types during AP initiation and propagation.

Axonal mechanisms mediating γ-aminobutyric acid receptor type A (GABA-A) inhibition of striatal dopamine release

Axons of dopaminergic neurons innervate the striatum where they contribute to movement and reinforcement learning. Past work has shown that striatal GABA tonically inhibits dopamine release, but whether GABA-A receptors directly modulate transmission or act indirectly through circuit elements is unresolved. Here, we use whole-cell and perforated-patch recordings to test for GABA-A receptors on the main dopaminergic neuron axons and branching processes within the striatum of adult mice. Application of GABA depolarized axons, but also decreased the amplitude of axonal spikes, limited propagation and reduced striatal dopamine release. The mechanism of inhibition involved sodium channel inactivation and shunting. Lastly, we show the positive allosteric modulator diazepam enhanced GABA-A currents on dopaminergic axons and directly inhibited release, but also likely acts by reducing excitation from cholinergic interneurons. Thus, we reveal the mechanisms of GABA-A receptor modulation of dopamine release and provide new insights into the actions of benzodiazepines within the striatum.

Relationship between oxygen consumption and neuronal activity in a defined neural circuit

BACKGROUND

Neuronal computations related to sensory and motor activity along with the maintenance of spike discharge, synaptic transmission, and associated housekeeping are energetically demanding. The most efficient metabolic process to provide large amounts of energy equivalents is oxidative phosphorylation and thus dependent on O2 consumption. Therefore, O2 levels in the brain are a critical parameter that influences neuronal function. Measurements of O2 consumption have been used to estimate the cost of neuronal activity; however, exploring these metabolic relationships in vivo and under defined experimental conditions has been limited by technical challenges.

RESULTS

We used isolated preparations of Xenopus laevis tadpoles to perform a quantitative analysis of O2 levels in the brain under in vivo-like conditions. We measured O2 concentrations in the hindbrain in relation to the spike discharge of the superior oblique eye muscle-innervating trochlear nerve as proxy for central nervous activity. In air-saturated bath Ringer solution, O2 levels in the fourth ventricle and adjacent, functionally intact hindbrain were close to zero. Inhibition of mitochondrial activity with potassium cyanide or fixation of the tissue with ethanol raised the ventricular O2 concentration to bath levels, indicating that the brain tissue consumed the available O2. Gradually increasing oxygenation of the Ringer solution caused a concurrent increase of ventricular O2 concentrations. Blocking spike discharge with the local anesthetics tricaine methanesulfonate diminished the O2 consumption by ~ 50%, illustrating the substantial O2 amount related to neuronal activity. In contrast, episodes of spontaneous trochlear nerve spike bursts were accompanied by transient increases of the O2 consumption with parameters that correlated with burst magnitude and duration.

CONCLUSIONS

Controlled experimental manipulations of both the O2 level as well as the neuronal activity under in vivo-like conditions allowed to quantitatively relate spike discharge magnitudes in a particular neuronal circuitry with the O2 consumption in this area. Moreover, the possibility to distinctly manipulate various functional parameters will yield more insight in the coupling between metabolic and neuronal activity. Thus, apart from providing quantitative empiric evidence for the link between physiologically relevant spontaneous spike discharge in the brain and O2-dependent metabolism, isolated amphibian preparations are promising model systems to further dissociate the O2 dynamics in relation to neuronal computations.

Ca2+ entry through NaV channels generates submillisecond axonal Ca2+ signaling

Calcium ions (Ca²⁺) are essential for many cellular signaling mechanisms and enter the cytosol mostly through voltage-gated calcium channels. Here, using high-speed Ca²⁺ imaging up to 20 kHz in the rat layer five pyramidal neuron axon we found that activity-dependent intracellular calcium concentration ([Ca²⁺]_i) in the axonal initial segment was only partially dependent on voltage-gated calcium channels. Instead, [Ca²⁺]_i changes were sensitive to the specific voltage-gated sodium (Na_V) channel blocker tetrodotoxin. Consistent with the conjecture that Ca²⁺ enters through the Na_V channel pore, the optically resolved I_Ca in the axon initial segment overlapped with the activation kinetics of Na_V channels and heterologous expression of Na_V1.2 in HEK-293 cells revealed a tetrodotoxin-sensitive [Ca²⁺]_i rise. Finally, computational simulations predicted that axonal [Ca²⁺]_i transients reflect a 0.4% Ca²⁺ conductivity of Na_V channels. The findings indicate that Ca²⁺ permeation through Na_V channels provides a submillisecond rapid entry route in Na_V-enriched domains of mammalian axons.

Nerve cells communicate using tiny electrical impulses called action potentials. Special proteins termed ion channels produce these electric signals by allowing specific charged particles, or ions, to pass in or out of cells across its membrane. When a nerve cell ‘fires’ an action potential, specific ion channels briefly open to let in a surge of positively charged ions which electrify the cell. Action potentials begin in the same place in each nerve cell, at an area called the axon initial segment. The large number of sodium channels at this site kick-start the influx of positively charged sodium ions ensuring that every action potential starts from the same place.

Previous research has shown that, when action potentials begin, the concentration of calcium ions at the axon initial segment also increases, but it was not clear which ion channels were responsible for this entry of calcium. Channels that are selective for calcium ions are the prime candidates for this process. However, research in squid nerve cells gave rise to an unexpected idea by suggesting that sodium channels may not exclusively let in sodium but also allow some calcium ions to pass through. Hanemaaijer, Popovic et al. therefore wanted to test the routes that calcium ions take and see whether the sodium channels in mammalian nerve cells are also permeable to calcium.

Experiments using fluorescent dyes to track the concentration of calcium in rat and human nerve cells showed that calcium ions accumulated at the axon initial segment when action potentials fired. Most of this increase in calcium could be stopped by treating the neurons with a toxin that prevents sodium channels from opening. Electrical manipulations of the cells revealed that, in this context, the calcium ions were effectively behaving like sodium ions. Human kidney cells were then engineered to produce the sodium channel protein. This confirmed that calcium and sodium ions were indeed both passing through the same channel.

These results shed new light on the relationship between calcium ions and sodium channels within the mammalian nervous system and that this interplay occurs at the axon initial segment of the cell. Genetic mutations that ‘nudge’ sodium channels towards favoring calcium entry are also found in patients with autism spectrum disorders, and so this new finding may contribute to our understanding of these conditions.

Complete neural stem cell (NSC) neuronal differentiation requires a branched chain amino acids-induced persistent metabolic shift towards energy metabolism

Neural stem cell (NSC) neuronal differentiation requires a metabolic shift towards oxidative phosphorylation. We now show that a branched-chain amino acids-driven, persistent metabolic shift toward energy metabolism is required for full neuronal maturation. We increased energy metabolism of differentiating neurons derived both from murine NSCs and human induced pluripotent stem cells (iPSCs) by supplementing the cell culture medium with a mixture composed of branched-chain amino acids, essential amino acids, TCA cycle precursors and co-factors. We found that treated differentiating neuronal cells with enhanced energy metabolism increased: i) total dendritic length; ii) the mean number of branches and iii) the number and maturation of the dendritic spines. Furthermore, neuronal spines in treated neurons appeared more stable with stubby and mushroom phenotype and with increased expression of molecules involved in synapse formation. Treated neurons modified their mitochondrial dynamics increasing the mitochondrial fusion and, consistently with the increase of cellular ATP content, they activated cellular mTORC1 dependent p70S6 K1 anabolism. Global transcriptomic analysis further revealed that treated neurons induce Nrf2 mediated gene expression. This was correlated with a functional increase in the Reactive Oxygen Species (ROS) scavenging mechanisms. In conclusion, persistent branched-chain amino acids-driven metabolic shift toward energy metabolism enhanced neuronal differentiation and antioxidant defences. These findings offer new opportunities to pharmacologically modulate NSC neuronal differentiation and to develop effective strategies for treating neurodegenerative diseases.

Theoretical relation between axon initial segment geometry and excitability

In most vertebrate neurons, action potentials are triggered at the distal end of the axon initial segment (AIS). Both position and length of the AIS vary across and within neuron types, with activity, development and pathology. What is the impact of AIS geometry on excitability? Direct empirical assessment has proven difficult because of the many potential confounding factors. Here, we carried a principled theoretical analysis to answer this question. We provide a simple formula relating AIS geometry and sodium conductance density to the somatic voltage threshold. A distal shift of the AIS normally produces a (modest) increase in excitability, but we explain how this pattern can reverse if a hyperpolarizing current is present at the AIS, due to resistive coupling with the soma. This work provides a theoretical tool to assess the significance of structural AIS plasticity for electrical function.

Saltatory Conduction along Myelinated Axons Involves a Periaxonal Nanocircuit

The propagation of electrical impulses along axons is highly accelerated by the myelin sheath and produces saltating or "jumping" action potentials across internodes, from one node of Ranvier to the next. The underlying electrical circuit, as well as the existence and role of submyelin conduction in saltatory conduction remain, however, elusive. Here, we made patch-clamp and high-speed voltage-calibrated optical recordings of potentials across the nodal and internodal axolemma of myelinated neocortical pyramidal axons combined with electron microscopy and experimentally constrained cable modeling. Our results reveal a nanoscale yet conductive periaxonal space, incompletely sealed at the paranodes, which separates the potentials across the low-capacitance myelin sheath and internodal axolemma. The emerging double-cable model reproduces the recorded evolution of voltage waveforms across nodes and internodes, including rapid nodal potentials traveling in advance of attenuated waves in the internodal axolemma, revealing a mechanism for saltation across time and space.

Intermittent Failure of Spike Propagation in Primary Afferent Neurons during Tactile Stimulation

Primary afferent neurons convey somatosensory information to the CNS. Low-threshold mechanoreceptors are classified as slow-adapting (SA) or rapid-adapting (RA) based on whether or not they spike repetitively during sustained tactile stimulation; the former are subclassified as Type 1 or 2 based on the regularity of their spiking. Recording in vivo from DRGs of mice, we observed irregular- and regular-spiking units consistent with SA1 and SA2 low-threshold mechanoreceptors, but some units, which we labeled "semiregular," did not fit cleanly into the existing classification scheme. Analysis of their spiking revealed integer-multiple patterning in which spike trains comprised a fundamental interspike interval and multiples thereof. Integer-multiple-patterned spiking was reproduced by randomly removing spikes from an otherwise regular spike train, suggesting that semiregular units represent SA2 units in which some spikes are "missing." We hypothesized that missing spikes arose from intermittent failure of spikes to initiate or to propagate. Intermittent failure of spike initiation was ruled out by several observations: integer-multiple-patterned spiking was not induced by intradermal lidocaine, was independent of stimulus modality (mechanical vs optogenetic), and could not be reproduced in a conductance-based model neuron given constant input. On the other hand, integer-multiple-patterned spiking was induced by application of lidocaine to the DRG, thus pinpointing intermittent failure of spike propagation as the basis for integer-multiple-patterned spiking. Indeed, half of all SA2 units exhibited some missing spikes, mostly at low rate (<5%), which suggests that axons are efficient in using the lowest safety factor capable of producing near-perfect propagation reliability.SIGNIFICANCE STATEMENT The impedance mismatch at axon branch points can impede spike propagation. Reliability of spike propagation across branch points remains an open question and is especially important for primary afferents whose spikes must cross a T-junction to reach the CNS. Past research on propagation reliability has relied almost entirely on simulations and in vitro experiments. Here, recording in vivo, we linked a distinctive pattern of spiking to the intermittent failure of spike propagation at the T-junction. The rarity of failures argues that safety factor is high under physiological conditions, yet the occurrence of such failures argues that safety factor is just high enough to ensure near-perfect reliability, consistent with a good balance between propagation reliability and energy efficiency.

Metabolic Cost of Dendritic Ca2+ Action Potentials in Layer 5 Pyramidal Neurons

Pyramidal neurons consume most signaling-related energy to generate action potentials (APs) and perform synaptic integration. Dendritic Ca²⁺ spike is an important integration mechanism for coupling inputs from different cortical layers. Our objective was to quantify the metabolic energy associated with the generation of Ca²⁺ APs in the dendrites. We used morphology-based computational models to simulate the dendritic Ca²⁺ spikes in layer 5 pyramidal neurons. We calculated the energy cost by converting Ca²⁺ influx into the number of ATP required to restore and maintain the homeostasis of intracellular Ca²⁺ concentrations. We quantified the effects of synaptic inputs, dendritic voltage, back-propagating Na⁺ spikes, and Ca²⁺ inactivation on Ca²⁺ spike cost. We showed that much more ATP molecules were required for reversing Ca²⁺ influx in the dendrites than for Na⁺ ion pumping in the soma during a Ca²⁺ AP. Increasing synaptic input increased the rate of dendritic depolarization and underlying Ca²⁺ influx, resulting in higher ATP consumption. Depolarizing dendritic voltage resulted in the inactivation of Ca²⁺ channels and reduced the ATP cost, while dendritic hyperpolarization increased the spike cost by de-inactivating Ca²⁺ channels. A back-propagating Na⁺ AP initiated in the soma increased Ca²⁺ spike cost in the apical dendrite when it coincided with a synaptic input within a time window of several milliseconds. Increasing Ca²⁺ inactivation rate reduced Ca²⁺ spike cost, while slowing Ca²⁺ inactivation increased the spike cost. The results revealed that the energy demand of a Ca²⁺ AP was dynamically dependent on the state of dendritic activity. These findings were important for predicting the energy budget for signaling in pyramidal cells, interpreting functional imaging data, and designing energy-efficient neuromorphic devices.

The energetic brain - A review from students to students

The past 20 years have resulted in unprecedented progress in understanding brain energy metabolism and its role in health and disease. In this review, which was initiated at the 14th International Society for Neurochemistry Advanced School, we address the basic concepts of brain energy metabolism and approach the question of why the brain has high energy expenditure. Our review illustrates that the vertebrate brain has a high need for energy because of the high number of neurons and the need to maintain a delicate interplay between energy metabolism, neurotransmission, and plasticity. Disturbances to the energetic balance, to mitochondria quality control or to glia-neuron metabolic interaction may lead to brain circuit malfunction or even severe disorders of the CNS. We cover neuronal energy consumption in neural transmission and basic ('housekeeping') cellular processes. Additionally, we describe the most common (glucose) and alternative sources of energy namely glutamate, lactate, ketone bodies, and medium chain fatty acids. We discuss the multifaceted role of non-neuronal cells in the transport of energy substrates from circulation (pericytes and astrocytes) and in the supply (astrocytes and microglia) and usage of different energy fuels. Finally, we address pathological consequences of disrupted energy homeostasis in the CNS.

HCN channel-mediated neuromodulation can control action potential velocity and fidelity in central axons

Hyperpolarization-activated cyclic-nucleotide-gated (HCN) channels control electrical rhythmicity and excitability in the heart and brain, but the function of HCN channels at the subcellular level in axons remains poorly understood. Here, we show that the action potential conduction velocity in both myelinated and unmyelinated central axons can be bidirectionally modulated by a HCN channel blocker, cyclic adenosine monophosphate (cAMP), and neuromodulators. Recordings from mouse cerebellar mossy fiber boutons show that HCN channels ensure reliable high-frequency firing and are strongly modulated by cAMP (EC50 40 µM; estimated endogenous cAMP concentration 13 µM). In addition, immunogold-electron microscopy revealed HCN2 as the dominating subunit in cerebellar mossy fibers. Computational modeling indicated that HCN2 channels control conduction velocity primarily by altering the resting membrane potential and are associated with significant metabolic costs. These results suggest that the cAMP-HCN pathway provides neuromodulators with an opportunity to finely tune energy consumption and temporal delays across axons in the brain.

Mechanisms for the maintenance and regulation of axonal energy supply

The unique polarization and high-energy demand of neurons necessitates specialized mechanisms to maintain energy homeostasis throughout the cell, particularly in the distal axon. Mitochondria play a key role in meeting axonal energy demand by generating adenosine triphosphate through oxidative phosphorylation. Recent evidence demonstrates how axonal mitochondrial trafficking and anchoring are coordinated to sense and respond to altered energy requirements. If and when these mechanisms are impacted in pathological conditions, such as injury and neurodegenerative disease, is an emerging research frontier. Recent evidence also suggests that axonal energy demand may be supplemented by local glial cells, including astrocytes and oligodendrocytes. In this review, we provide an updated discussion of how oxidative phosphorylation, aerobic glycolysis, and oligodendrocyte-derived metabolic support contribute to the maintenance of axonal energy homeostasis.

Complementary Tuning of Na+ and K+ Channel Gating Underlies Fast and Energy-Efficient Action Potentials in GABAergic Interneuron Axons

Fast-spiking, parvalbumin-expressing GABAergic interneurons (PV⁺-BCs) express a complex machinery of rapid signaling mechanisms, including specialized voltage-gated ion channels to generate brief action potentials (APs). However, short APs are associated with overlapping Na⁺ and K⁺ fluxes and are therefore energetically expensive. How the potentially vicious combination of high AP frequency and inefficient spike generation can be reconciled with limited energy supply is presently unclear. To address this question, we performed direct recordings from the PV⁺-BC axon, the subcellular structure where active conductances for AP initiation and propagation are located. Surprisingly, the energy required for the AP was, on average, only ∼1.6 times the theoretical minimum. High energy efficiency emerged from the combination of fast inactivation of Na⁺ channels and delayed activation of Kv3-type K⁺ channels, which minimized ion flux overlap during APs. Thus, the complementary tuning of axonal Na⁺ and K⁺ channel gating optimizes both fast signaling properties and metabolic efficiency.