The Slow Dynamics of Intracellular Sodium Concentration Increase the Time Window of Neuronal Integration: A Simulation Study

Uncertainty quantification and sensitivity analysis of neuron models with ion concentration dynamics

This paper provides a comprehensive and computationally efficient case study for uncertainty quantification (UQ) and global sensitivity analysis (GSA) in a neuron model incorporating ion concentration dynamics. We address how challenges with UQ and GSA in this context can be approached and solved, including challenges related to computational cost, parameters affecting the system's resting state, and the presence of both fast and slow dynamics. Specifically, we analyze the electrodiffusive neuron-extracellular-glia (edNEG) model, which captures electrical potentials, ion concentrations (Na+, K+, Ca2+, and Cl-), and volume changes across six compartments. Our methodology includes a UQ procedure assessing the model's reliability and susceptibility to input uncertainty and a variance-based GSA identifying the most influential input parameters. To mitigate computational costs, we employ surrogate modeling techniques, optimized using efficient numerical integration methods. We propose a strategy for isolating parameters affecting the resting state and analyze the edNEG model dynamics under both physiological and pathological conditions. The influence of uncertain parameters on model outputs, particularly during spiking dynamics, is systematically explored. Rapid dynamics of membrane potentials necessitate a focus on informative spiking features, while slower variations in ion concentrations allow a meaningful study at each time point. Our study offers valuable guidelines for future UQ and GSA investigations on neuron models with ion concentration dynamics, contributing to the broader application of such models in computational neuroscience.

Stimulus-Induced Theta-Band LFP Oscillations Format Neuronal Representations of Social Chemosignals in the Mouse Accessory Olfactory Bulb

Social communication is crucial for the survival of many species. In most vertebrates, a dedicated chemosensory system, the vomeronasal system (VNS), evolved to process ethologically relevant chemosensory cues. The first central processing stage of the VNS is the accessory olfactory bulb (AOB), which sends information to downstream brain regions via AOB mitral cells (AMCs). Recent studies provided important insights about the functional properties of AMCs, but little is known about the principles that govern their coordinated activity. Here, we recorded local field potentials (LFPs) and single-unit activity in the AOB of adult male and female mice during presentation of natural stimuli. Our recordings reveal prominent LFP theta-band oscillatory episodes with a characteristic spatial pattern across the AOB. Throughout an experiment, the AOB network shows varying degrees of similarity to this pattern, in a manner that depends on the sensory stimulus. Analysis of LFP signal polarity and single-unit activity indicates that oscillatory episodes are generated locally within the AOB, likely representing a reciprocal interaction between AMCs and granule cells. Notably, spike times of many AMCs are constrained to the negative LFP oscillation phase in a manner that can drastically affect integration by downstream processing stages. Based on these observations, we propose that LFP oscillations may gate, bind, and organize outgoing signals from individual AOB neurons to downstream processing stages. Our findings suggest that, as in other neuronal systems and brain regions, population-level oscillations play a key role in organizing and enhancing transmission of socially relevant chemosensory information.SIGNIFICANCE STATEMENT The accessory olfactory bulb (AOB) is the first central stage of the vomeronasal system, a chemosensory system dedicated to processing cues from other organisms. Information from the AOB is conveyed to other brain regions via activity of its principal neurons, AOB mitral cells (AMCs). Here, we show that socially relevant sensory stimulation of the mouse vomeronasal system leads not only to changes in AMC activity, but also to distinct theta-band (∼5 Hz) oscillatory episodes in the local field potential. Notably AMCs favor the negative phase of these oscillatory events. Our findings suggest a novel mechanism for the temporal coordination of distributed patterns of neuronal activity, which can serve to efficiently activate downstream processing stages.

Temporal dynamics of Na/K pump mediated memory traces: insights from conductance-based models of Drosophila neurons

Sodium potassium ATPases (Na/K pumps) mediate long-lasting, dynamic cellular memories that can last tens of seconds. The mechanisms controlling the dynamics of this type of cellular memory are not well understood and can be counterintuitive. Here, we use computational modeling to examine how Na/K pumps and the ion concentration dynamics they influence shape cellular excitability. In a Drosophila larval motor neuron model, we incorporate a Na/K pump, a dynamic intracellular Na⁺ concentration, and a dynamic Na⁺ reversal potential. We probe neuronal excitability with a variety of stimuli, including step currents, ramp currents, and zap currents, then monitor the sub- and suprathreshold voltage responses on a range of time scales. We find that the interactions of a Na⁺-dependent pump current with a dynamic Na⁺ concentration and reversal potential endow the neuron with rich response properties that are absent when the role of the pump is reduced to the maintenance of constant ion concentration gradients. In particular, these dynamic pump-Na⁺ interactions contribute to spike rate adaptation and result in long-lasting excitability changes after spiking and even after sub-threshold voltage fluctuations on multiple time scales. We further show that modulation of pump properties can profoundly alter a neuron's spontaneous activity and response to stimuli by providing a mechanism for bursting oscillations. Our work has implications for experimental studies and computational modeling of the role of Na/K pumps in neuronal activity, information processing in neural circuits, and the neural control of animal behavior.

Identification of Sodium- and Chloride-Sensitive Sites in the Slack Channel

The Slack channel (KCNT1, Slo2.2) is a sodium-activated and chloride-activated potassium channel that regulates heart rate and maintains the normal excitability of the nervous system. Despite intense interest in the sodium gating mechanism, a comprehensive investigation to identify the sodium-sensitive and chloride-sensitive sites has been missing. In the present study, we identified two potential sodium-binding sites in the C-terminal domain of the rat Slack channel by conducting electrophysical recordings and systematic mutagenesis of cytosolic acidic residues in the rat Slack channel C terminus. In particular, by taking advantage of the M335A mutant, which results in the opening of the Slack channel in the absence of cytosolic sodium, we found that among the 92 screened negatively charged amino acids, E373 mutants could completely remove sodium sensitivity of the Slack channel. In contrast, several other mutants showed dramatic decreases in sodium sensitivity but did not abolish it altogether. Furthermore, molecular dynamics (MD) simulations performed at the hundreds of nanoseconds timescale revealed one or two sodium ions at the E373 position or an acidic pocket composed of several negatively charged residues. Moreover, the MD simulations predicted possible chloride interaction sites. By screening predicted positively charged residues, we identified R379 as a chloride interaction site. Thus, we conclude that the E373 site and the D863/E865 pocket are two potential sodium-sensitive sites, while R379 is a chloride interaction site in the Slack channel.SIGNIFICANCE STATEMENT The research presented here identified two distinct sodium and one chloride interaction sites located in the intracellular C-terminal domain of the Slack (Slo2.2, KCNT1) channel. Identification of the sites responsible for the sodium and chloride activation of the Slack channel sets its gating property apart from other potassium channels in the BK channel family. This finding sets the stage for future functional and pharmacological studies of this channel.

Recurrent network interactions explain tectal response variability and experience-dependent behavior

Response variability is an essential and universal feature of sensory processing and behavior. It arises from fluctuations in the internal state of the brain, which modulate how sensory information is represented and transformed to guide behavioral actions. In part, brain state is shaped by recent network activity, fed back through recurrent connections to modulate neuronal excitability. However, the degree to which these interactions influence response variability and the spatial and temporal scales across which they operate, are poorly understood. Here, we combined population recordings and modeling to gain insights into how neuronal activity modulates network state and thereby impacts visually evoked activity and behavior. First, we performed cellular-resolution calcium imaging of the optic tectum to monitor ongoing activity, the pattern of which is both a cause and consequence of changes in network state. We developed a minimal network model incorporating fast, short range, recurrent excitation and long-lasting, activity-dependent suppression that reproduced a hallmark property of tectal activity - intermittent bursting. We next used the model to estimate the excitability state of tectal neurons based on recent activity history and found that this explained a portion of the trial-to-trial variability in visually evoked responses, as well as spatially selective response adaptation. Moreover, these dynamics also predicted behavioral trends such as selective habituation of visually evoked prey-catching. Overall, we demonstrate that a simple recurrent interaction motif can be used to estimate the effect of activity upon the incidental state of a neural network and account for experience-dependent effects on sensory encoding and visually guided behavior.

Inhibition of Voltage-Gated Na+ Currents Exerted by KB-R7943 (2-[2-[4-(4-nitrobenzyloxy)phenyl]ethyl]isothiourea), an Inhibitor of Na+-Ca2+ Exchanging Process

KB-R7943, an isothiourea derivative, has been recognized as an inhibitor in the reverse mode of the Na⁺-Ca²⁺ exchanging process. This compound was demonstrated to prevent intracellular Na⁺-dependent Ca²⁺ uptake in intact cells; however, it is much less effective at preventing extracellular Na⁺-dependent Ca²⁺ efflux. Therefore, whether or how this compound may produce any perturbations on other types of ionic currents, particularly on voltage-gated Na⁺ current (I_Na), needs to be further studied. In this study, the whole-cell current recordings demonstrated that upon abrupt depolarization in pituitary GH₃ cells, the exposure to KB-R7943 concentration-dependently depressed the transient (I_Na(T)) or late component (I_Na(L)) of I_Na with an IC₅₀ value of 11 or 0.9 μM, respectively. Likewise, the dissociation constant for the KB-R7943-mediated block of I_Na on the basis of a minimum reaction scheme was estimated to be 0.97 μM. The presence of benzamil or amiloride could suppress the I_Na(L) magnitude. The instantaneous window Na⁺ current (I_Na(W)) activated by abrupt ascending ramp voltage (V_ramp) was suppressed by adding KB-R7943; however, subsequent addition of deltamethrin or tefluthrin (Tef) effectively reversed KB-R7943-inhibted I_Na(W). With prolonged duration of depolarizing pulses, the I_Na(L) amplitude became exponentially decreased; moreover, KB-R7943 diminished I_Na(L) magnitude. The resurgent Na⁺ current (I_Na(R)) evoked by a repolarizing V_ramp was also suppressed by adding this compound; moreover, subsequent addition of ranolazine or Tef further diminished or reversed, respectively, its reduction in I_Na(R) magnitude. The persistent Na⁺ current (I_Na(P)) activated by sinusoidal voltage waveform became enhanced by Tef; however, subsequent application of KB-R7943 counteracted Tef-stimulated I_Na(P). The docking prediction reflected that there seem to be molecular interactions of this molecule with the hNa_V1.2 or hNa_V1.7 channels. Collectively, this study highlights evidence showing that KB-R7943 has the propensity to perturb the magnitude and gating kinetics of I_Na (e.g., I_Na(T), I_Na(L), I_Na(W), I_Na(R), and I_Na(P)) and that the Na_V channels appear to be important targets for the in vivo actions of KB-R7943 or other relevant compounds.

Phytochemical analysis and toxicological evaluation of the ethanolic Leaves extract of Hypoestes rosea on the morphology and biochemical indices of the Kidneys of albino Wistar Rats

BACKGROUND

Hypoestes rosea (family: Acanthacea), has been harnessed and utilized for treatment of several ailments. However, there is the paucity of available data on nephrotoxicity associated with this herb. Here, we investigated the phytochemical profile and toxicological effect of H. rosea on Wistar Rats.

METHODS

Twenty rats (weight range: 75-100 g) were assigned into five study groups, viz; (a) control (without treatment) (b) treatment group 1, orally administered with 50 mg/kg (c) treatment group 2, orally administered with 100 mg/kg (d) treatment group 3, orally administered with 250 mg/kg, and (e) treatment group 4, orally administered with 300 mg/kg of H. rosea, respectively for 28 days of four rats per group. The rats were made unconscious by using oral administration of chloroform. Cardiac punctures were made, and blood samples collected into 10 ml labeled plain container, allowed to clot and spun to harvest serum for determination of sodium, potassium, chloride, bicarbonate, urea and creatinine using colorimetric, back-titrimetric, Urease-Berthelot and Jaffe's reaction methods respectively. Kidneys of rats were harvested, weighed and immediately fixed in 10% neutral buffered formalin for histological analysis.

RESULT

Mean serum sodium (p = 0.049), potassium (p = 0.007), and urea (p < 0.001) levels were significantly higher among the treatment groups compared to controls. Histopathological findings of kidney sections revealed mild glomerular infiltration in treatment groups 2-4. Additionally, sclerosis was observed in groups 3-4. Phytochemical analysis of H. rosea revealed presence of alkaloids, flavonoids, saponins, tannins, terpenoids, steroids and reducing sugars.

CONCLUSION

From the findings in this study, H. rosea leaf extract causes significant damage to the kidneys of Wistar rats at higher doses. Of which, the damages were dose-dependent in direct proportionality manner. To better determine the safe dosage and ideal duration of consumption, there is the need for further studies on H. rosea.

Interactions among diameter, myelination, and the Na/K pump affect axonal resilience to high-frequency spiking

The reliability of spike propagation in axons is determined by complex interactions among ionic currents, ion pumps, and morphological properties. We use compartment-based modeling to reveal that interactions of diameter, myelination, and the Na/K pump determine the reliability of high-frequency spike propagation. By acting as a “reservoir” of nodal Na⁺ influx, myelinated compartments efficiently increase propagation reliability. Although spike broadening was thought to oppose fast spiking, its effect on spike propagation is complicated, depending on the balance of Na⁺ channel inactivation gate recovery, Na⁺ influx, and axial charge. Our findings suggest that slow Na⁺ removal influences axonal resilience to high-frequency spike propagation and that different strategies may be required to overcome this constraint in different neurons.

Axons reliably conduct action potentials between neurons and/or other targets. Axons have widely variable diameters and can be myelinated or unmyelinated. Although the effect of these factors on propagation speed is well studied, how they constrain axonal resilience to high-frequency spiking is incompletely understood. Maximal firing frequencies range from ∼1 Hz to >300 Hz across neurons, but the process by which Na/K pumps counteract Na⁺ influx is slow, and the extent to which slow Na⁺ removal is compatible with high-frequency spiking is unclear. Modeling the process of Na⁺ removal shows that large-diameter axons are more resilient to high-frequency spikes than are small-diameter axons, because of their slow Na⁺ accumulation. In myelinated axons, the myelinated compartments between nodes of Ranvier act as a “reservoir” to slow Na⁺ accumulation and increase the reliability of axonal propagation. We now find that slowing the activation of K⁺ current can increase the Na⁺ influx rate, and the effect of minimizing the overlap between Na⁺ and K⁺ currents on spike propagation resilience depends on complex interactions among diameter, myelination, and the Na/K pump density. Our results suggest that, in neurons with different channel gating kinetic parameters, different strategies may be required to improve the reliability of axonal propagation.

An electrodiffusive neuron-extracellular-glia model for exploring the genesis of slow potentials in the brain

Within the computational neuroscience community, there has been a focus on simulating the electrical activity of neurons, while other components of brain tissue, such as glia cells and the extracellular space, are often neglected. Standard models of extracellular potentials are based on a combination of multicompartmental models describing neural electrodynamics and volume conductor theory. Such models cannot be used to simulate the slow components of extracellular potentials, which depend on ion concentration dynamics, and the effect that this has on extracellular diffusion potentials and glial buffering currents. We here present the electrodiffusive neuron-extracellular-glia (edNEG) model, which we believe is the first model to combine compartmental neuron modeling with an electrodiffusive framework for intra- and extracellular ion concentration dynamics in a local piece of neuro-glial brain tissue. The edNEG model (i) keeps track of all intraneuronal, intraglial, and extracellular ion concentrations and electrical potentials, (ii) accounts for action potentials and dendritic calcium spikes in neurons, (iii) contains a neuronal and glial homeostatic machinery that gives physiologically realistic ion concentration dynamics, (iv) accounts for electrodiffusive transmembrane, intracellular, and extracellular ionic movements, and (v) accounts for glial and neuronal swelling caused by osmotic transmembrane pressure gradients. The edNEG model accounts for the concentration-dependent effects on ECS potentials that the standard models neglect. Using the edNEG model, we analyze these effects by splitting the extracellular potential into three components: one due to neural sink/source configurations, one due to glial sink/source configurations, and one due to extracellular diffusive currents. Through a series of simulations, we analyze the roles played by the various components and how they interact in generating the total slow potential. We conclude that the three components are of comparable magnitude and that the stimulus conditions determine which of the components that dominate.

An electrodiffusive, ion conserving Pinsky-Rinzel model with homeostatic mechanisms

In most neuronal models, ion concentrations are assumed to be constant, and effects of concentration variations on ionic reversal potentials, or of ionic diffusion on electrical potentials are not accounted for. Here, we present the electrodiffusive Pinsky-Rinzel (edPR) model, which we believe is the first multicompartmental neuron model that accounts for electrodiffusive ion concentration dynamics in a way that ensures a biophysically consistent relationship between ion concentrations, electrical charge, and electrical potentials in both the intra- and extracellular space. The edPR model is an expanded version of the two-compartment Pinsky-Rinzel (PR) model of a hippocampal CA3 neuron. Unlike the PR model, the edPR model includes homeostatic mechanisms and ion-specific leakage currents, and keeps track of all ion concentrations (Na⁺, K⁺, Ca²⁺, and Cl⁻), electrical potentials, and electrical conductivities in the intra- and extracellular space. The edPR model reproduces the membrane potential dynamics of the PR model for moderate firing activity. For higher activity levels, or when homeostatic mechanisms are impaired, the homeostatic mechanisms fail in maintaining ion concentrations close to baseline, and the edPR model diverges from the PR model as it accounts for effects of concentration changes on neuronal firing. We envision that the edPR model will be useful for the field in three main ways. Firstly, as it relaxes commonly made modeling assumptions, the edPR model can be used to test the validity of these assumptions under various firing conditions, as we show here for a few selected cases. Secondly, the edPR model should supplement the PR model when simulating scenarios where ion concentrations are expected to vary over time. Thirdly, being applicable to conditions with failed homeostasis, the edPR model opens up for simulating a range of pathological conditions, such as spreading depression or epilepsy.

Neurons generate their electrical signals by letting ions pass through their membranes. Despite this fact, most models of neurons apply the simplifying assumption that ion concentrations remain effectively constant during neural activity. This assumption is often quite good, as neurons contain a set of homeostatic mechanisms that make sure that ion concentrations vary quite little under normal circumstances. However, under some conditions, these mechanisms can fail, and ion concentrations can vary quite dramatically. Standard models are thus not able to simulate such conditions. Here, we present what to our knowledge is the first multicompartmental neuron model that accounts for ion concentration variations in a way that ensures complete and consistent ion concentration and charge conservation. In this work, we use the model to explore under which activity conditions the ion concentration variations become important for predicting the neurodynamics. We expect the model to be of great value for the field of neuroscience, as it can be used to simulate a range of pathological conditions, such as spreading depression or epilepsy, which are associated with large changes in extracellular ion concentrations.