A Theoretical Study on the Role of Astrocytic Activity in Neuronal Hyperexcitability by a Novel Neuron-Glia Mass Model

A neural mass model for disturbance of alpha rhythm in the minimal hepatic encephalopathy

One of the early markers of minimal hepatic encephalopathy (MHE) is the disruption of alpha rhythm observed in electroencephalogram (EEG) signals. However, the underlying mechanisms responsible for this occurrence remain poorly understood. To address this gap, we develop a novel biophysical model MHE-AWD-NCM, encompassing the communication dynamics between a cortical neuron population (CNP) and an astrocyte population (AP), aimed at investigating the relationship between alpha wave disturbance (AWD) and mechanistical principles, specifically concerning astrocyte-neuronal communication in the context of MHE. In addition, we introduce the concepts of peak power density and peak frequency within the alpha band as quantitative measures of AWD. Our model faithfully reproduces the characteristic EEG phenomenology during MHE and shows how impairments of communication between CNP and AP could promote AWD. The results suggest that the disruptions in feedback neurotransmission from AP to CNP, along with the inhibition of GABA uptake by AP from the extracellular space, contribute to the observed AWD. Moreover, we found that the variation of external excitatory stimuli on CNP may play a key role in AWD in MHE. Finally, the sensitivity analysis is also performed to assess the relative significance of above factors in influencing AWD. Our findings align with the physiological observations and provide a more comprehensive understanding of the complex interplay of astrocyte-neuronal communication that underlies the AWD observed in MHE, which potentially may help to explore the targeted therapeutic interventions for the early stage of hepatic encephalopathy.

Investigating the ability of astrocytes to drive neural network synchrony

Recent experimental works have implicated astrocytes as a significant cell type underlying several neuronal processes in the mammalian brain, from encoding sensory information to neurological disorders. Despite this progress, it is still unclear how astrocytes are communicating with and driving their neuronal neighbors. While previous computational modeling works have helped propose mechanisms responsible for driving these interactions, they have primarily focused on interactions at the synaptic level, with microscale models of calcium dynamics and neurotransmitter diffusion. Since it is computationally infeasible to include the intricate microscale details in a network-scale model, little computational work has been done to understand how astrocytes may be influencing spiking patterns and synchronization of large networks. We overcome this issue by first developing an "effective" astrocyte that can be easily implemented to already established network frameworks. We do this by showing that the astrocyte proximity to a synapse makes synaptic transmission faster, weaker, and less reliable. Thus, our "effective" astrocytes can be incorporated by considering heterogeneous synaptic time constants, which are parametrized only by the degree of astrocytic proximity at that synapse. We then apply our framework to large networks of exponential integrate-and-fire neurons with various spatial structures. Depending on key parameters, such as the number of synapses ensheathed and the strength of this ensheathment, we show that astrocytes can push the network to a synchronous state and exhibit spatially correlated patterns.

From phenomenological to biophysical models of seizures

Epilepsy is a complex disease that requires various approaches for its study. In this short review, we discuss the contribution of theoretical and computational models. The review presents theoretical frameworks that underlie the understanding of certain seizure properties and their classification based on their dynamical properties at the onset and offset of seizures. Dynamical system tools are valuable resources in the study of seizures. By analyzing the complex, dynamic behavior of seizures, these tools can provide insights into seizure mechanisms and offer a framework for their classification. Additionally, computational models have high potential for clinical applications, as they can be used to develop more accurate diagnostic and personalized medicine tools. We discuss various modeling approaches that span different scales and levels, while also questioning the neurocentric view, and emphasize the importance of considering glial cells. Finally, we explore the epistemic value provided by this type of approach.

Early deficits in GABA inhibition parallels an increase in L-type Ca2+ currents in the jaw motor neurons of SOD1G93A mouse model for ALS

Amyotrophic Lateral Sclerosis (ALS) involves protracted pre-symptomatic periods of abnormal motor neuron (MN) excitability occurring in parallel with central and peripheral synaptic perturbations. Focusing on inhibitory control of MNs, we first compared longitudinal changes in pre-synaptic terminal proteins for GABA and glycine neurotransmitters around the soma of retrogradely identified trigeminal jaw closer (JC) MNs and ChAT-labeled midbrain extraocular (EO) MNs in the SOD1^G93A mouse model for ALS. Fluorescence immunocytochemistry and confocal imaging were used to quantify GAD67 and GlyT2 synaptic bouton density (SBD) around MN soma at pre-symptomatic ages ∼P12 (postnatal), ∼P50 (adult) and near disease end-stage (∼P135) in SOD1^G93A mice and age-matched wild-type (WT) controls. We noted reduced GAD67 innervation in the SOD1^G93A trigeminal jaw closer MNs around P12, relative to age-matched WT and no significant difference around P50 and P135. In contrast, both GAD67 and GlyT2 innervation were elevated in the SOD1^G93A EO MNs at the pre-symptomatic time points. Considering trigeminal MNs are vulnerable in ALS while EO MNs are spared, we suggest that upregulation of inhibition in the latter might be compensatory. Notable contrast also existed in the innate co-expression patterns of GAD67 and GlyT2 with higher mutual information (co-dependency) in EO MNs compared to JC in both SOD1^G93A and WT mice, especially at adult stages (P50 and P135). Around P12 when GAD67 terminals expression was low in the mutant, we further tested for persistent GABA inhibition in those MNs using in vitro patch-clamp electrophysiology. Our results show that SOD1^G93A JC MNs have reduced persistent GABA inhibition, relative to WT. Pharmacological blocking of an underlying tonically active GABA conductance using the GABA-α5 subunit inverse agonist, L-655-708, disinhibited WT JC MNs and lowered their recruitment threshold, suggesting its role in the control of intrinsic MN excitability. Quantitative RT-PCR in laser dissected JC MNs further supported a reduction in GABA-α5 subunit mRNA expression in the mutant. In light of our previous report that JC MNs forming putative fast motor units have lower input threshold in the SOD1^G93A mice, we suggest that our present result on reduced GABA-α5 tonic inhibition provides for a mechanism contributing to such imbalance. In parallel with reduced GABA inhibition, we noted an increase in voltage-gated L-type Ca²⁺ currents in the mutant JC MNs around P12. Together these results support that, early modifications in intrinsic properties of vulnerable MNs could be an adaptive response to counter synaptic deficits.

Dextromethorphan Dampens Neonatal Astrocyte Activation and Endoplasmic Reticulum Stress Induced by Prenatal Exposure to Buprenorphine

Prenatal exposure to buprenorphine renders offspring vulnerable to cerebral impairments. In this study, our data demonstrate, for the first time, that prenatal exposure to buprenorphine escalates astrocyte activation concurrent with indications of endoplasmic reticulum (ER) stress in the hippocampi of neonates, and this can be prevented by the coadministration of dextromethorphan with buprenorphine. Furthermore, dextromethorphan can inhibit the accumulation of GPR37 in the hippocampus of newborns caused by buprenorphine and is accompanied by the proapoptotic ER stress response that involves the procaspase-3/CHOP pathway. Primary astrocyte cultures derived from the neonates of the buprenorphine group also displayed aberrant ER calcium mobilization and elevated basal levels of cyclooxygenase-2 (COX-2) at 14 days in vitro while showing sensitivity to lipopolysaccharide-activated expression of COX-2. Similarly, these long-lasting defects in the hippocampus and astrocytes were abolished by dextromethorphan. Our findings suggest that prenatal exposure to buprenorphine might instigate long-lasting effects on hippocampal and astrocytic functions. The beneficial effects of prenatal coadministration of dextromethorphan might be, at least in part, attributed to its properties in attenuating astrocyte activation and hippocampal ER stress in neonates.

Computational analysis of NIRS and BOLD signal from neurovascular coupling with three neuron-system feedforward inhibition network

Several neurological disorders occur due to hypoxic condition in brain arising from impairment of cerebral functionality, which can be controlled by neural stimulation driven vasoactive response mediated through biological response in astrocyte, a phenomenon known as neurovascular coupling. Brain can adjust with the problem of hypoxic condition by causing vasodilation with the help of this mechanism. To deduce the mechanism behind vasodilation of blood vessel caused by neuronal stimulus, current study articulates a mathematical model involving neuronal system feedforward inhibition network model (FFI) with two other functional components of neurovascular coupling, i.e. astrocyte and smooth muscle cell lining blood vessel. This study includes the neural inhibition network system where glutamatergic pyramidal neuron and GABAergic interneuron act antagonistically with each other. The proposed model successfully includes the implication of the inhibition system to design mathematical model for neurovascular coupling. Result of the proposed model shows that the increase in neuronal stimulus from 20 to 60 µA/cm² has the ability to increase the vasodilatory activity of blood tissue vasculature. Oxygenation level and hemodynamic response due to input synaptic stimulation has been calculated by regional cerebral oxygenation level (rS0₂) and blood oxygen level dependent (BOLD) imaging signal which supports vasodilation of blood vessel with increase in synaptic input stimulus.

Accuracy and spatial properties of distributed magnetic source imaging techniques in the investigation of focal epilepsy patients

Source localization of interictal epileptiform discharges (IEDs) is clinically useful in the presurgical workup of epilepsy patients. We aimed to compare the performance of four different distributed magnetic source imaging (dMSI) approaches: Minimum norm estimate (MNE), dynamic statistical parametric mapping (dSPM), standardized low-resolution electromagnetic tomography (sLORETA), and coherent maximum entropy on the mean (cMEM). We also evaluated whether a simple average of maps obtained from multiple inverse solutions (Ave) can improve localization accuracy. We analyzed dMSI of 206 IEDs derived from magnetoencephalography recordings in 28 focal epilepsy patients who had a well-defined focus determined through intracranial EEG (iEEG), epileptogenic MRI lesions or surgical resection. dMSI accuracy and spatial properties were quantitatively estimated as: (a) distance from the epilepsy focus, (b) reproducibility, (c) spatial dispersion (SD), (d) map extension, and (e) effect of thresholding on map properties. Clinical performance was excellent for all methods (median distance from the focus MNE = 2.4 mm; sLORETA = 3.5 mm; cMEM = 3.5 mm; dSPM = 6.8 mm, Ave = 0 mm). Ave showed the lowest distance between the map maximum and epilepsy focus (Dmin lower than cMEM, MNE, and dSPM, p = .021, p = .008, p < .001, respectively). cMEM showed the best spatial features, with lowest SD outside the focus (SD lower than all other methods, p < .001 consistently) and high contrast between the generator and surrounding regions. The average map Ave provided the best localization accuracy, whereas cMEM exhibited the lowest amount of spurious distant activity. dMSI techniques have the potential to significantly improve identification of iEEG targets and to guide surgical planning, especially when multiple methods are combined.

Differential temperature sensitivity of synaptic and firing processes in a neural mass model of epileptic discharges explains heterogeneous response of experimental epilepsy to focal brain cooling

Experiments with drug-induced epilepsy in rat brains and epileptic human brain region reveal that focal cooling can suppress epileptic discharges without affecting the brain's normal neurological function. Findings suggest a viable treatment for intractable epilepsy cases via an implantable cooling device. However, precise mechanisms by which cooling suppresses epileptic discharges are still not clearly understood. Cooling experiments in vitro presented evidence of reduction in neurotransmitter release from presynaptic terminals and loss of dendritic spines at post-synaptic terminals offering a possible synaptic mechanism. We show that termination of epileptic discharges is possible by introducing a homogeneous temperature factor in a neural mass model which attenuates the post-synaptic impulse responses of the neuronal populations. This result however may be expected since such attenuation leads to reduced post-synaptic potential and when the effect on inhibitory interneurons is less than on excitatory interneurons, frequency of firing of pyramidal cells is consequently reduced. While this is observed in cooling experiments in vitro, experiments in vivo exhibit persistent discharges during cooling but suppressed in magnitude. This leads us to conjecture that reduction in the frequency of discharges may be compensated through intrinsic excitability mechanisms. Such compensatory mechanism is modelled using a reciprocal temperature factor in the firing response function in the neural mass model. We demonstrate that the complete model can reproduce attenuation of both magnitude and frequency of epileptic discharges during cooling. The compensatory mechanism suggests that cooling lowers the average and the variance of the distribution of threshold potential of firing across the population. Bifurcation study with respect to the temperature parameters of the model reveals how heterogeneous response of epileptic discharges to cooling (termination or suppression only) is exhibited. Possibility of differential temperature effects on post-synaptic potential generation of different populations is also explored.