Spike-Conducting Integrate-and-Fire Model

Computational Modeling of the Prefrontal-Cingulate Cortex to Investigate the Role of Coupling Relationships for Balancing Emotion and Cognition

Within the prefrontal-cingulate cortex, abnormalities in coupling between neuronal networks can disturb the emotion-cognition interactions, contributing to the development of mental disorders such as depression. Despite this understanding, the neural circuit mechanisms underlying this phenomenon remain elusive. In this study, we present a biophysical computational model encompassing three crucial regions, including the dorsolateral prefrontal cortex, subgenual anterior cingulate cortex, and ventromedial prefrontal cortex. The objective is to investigate the role of coupling relationships within the prefrontal-cingulate cortex networks in balancing emotions and cognitive processes. The numerical results confirm that coupled weights play a crucial role in the balance of emotional cognitive networks. Furthermore, our model predicts the pathogenic mechanism of depression resulting from abnormalities in the subgenual cortex, and network functionality was restored through intervention in the dorsolateral prefrontal cortex. This study utilizes computational modeling techniques to provide an insight explanation for the diagnosis and treatment of depression.

Modeling and simulation for prediction of multiple sclerosis progression

In light of extensive work that has created a wide range of techniques for predicting the course of multiple sclerosis (MS) disease, this paper attempts to provide an overview of these approaches and put forth an alternative way to predict the disease progression. For this purpose, the existing methods for estimating and predicting the course of the disease have been categorized into clinical, radiological, biological, and computational or artificial intelligence-based markers. Weighing the weaknesses and strengths of these prognostic groups is a profound method that is yet in need and works directly at the level of diseased connectivity. Therefore, we propose using the computational models in combination with established connectomes as a predictive tool for MS disease trajectories. The fundamental conduction-based Hodgkin-Huxley model emerged as promising from examining these studies. The advantage of the Hodgkin-Huxley model is that certain properties of connectomes, such as neuronal connection weights, spatial distances, and adjustments of signal transmission rates, can be taken into account. It is precisely these properties that are particularly altered in MS and that have strong implications for processing, transmission, and interactions of neuronal signaling patterns. The Hodgkin-Huxley (HH) equations as a point-neuron model are used for signal propagation inside a small network. The objective is to change the conduction parameter of the neuron model, replicate the changes in myelin properties in MS and observe the dynamics of the signal propagation across the network. The model is initially validated for different lengths, conduction values, and connection weights through three nodal connections. Later, these individual factors are incorporated into a small network and simulated to mimic the condition of MS. The signal propagation pattern is observed after inducing changes in conduction parameters at certain nodes in the network and compared against a control model pattern obtained before the changes are applied to the network. The signal propagation pattern varies as expected by adapting to the input conditions. Similarly, when the model is applied to a connectome, the pattern changes could give an insight into disease progression. This approach has opened up a new path to explore the progression of the disease in MS. The work is in its preliminary state, but with a future vision to apply this method in a connectome, providing a better clinical tool.

A dynamical computational model of theta generation in hippocampal circuits to study theta-gamma oscillations during neurostimulation

Neurostimulation of the hippocampal formation has shown promising results for modulating memory but the underlying mechanisms remain unclear. In particular, the effects on hippocampal theta-nested gamma oscillations and theta phase reset, which are both crucial for memory processes, are unknown. Moreover, these effects cannot be investigated using current computational models, which consider theta oscillations with a fixed amplitude and phase velocity. Here, we developed a novel computational model that includes the medial septum, represented as a set of abstract Kuramoto oscillators producing a dynamical theta rhythm with phase reset, and the hippocampal formation, composed of biophysically realistic neurons and able to generate theta-nested gamma oscillations under theta drive. We showed that, for theta inputs just below the threshold to induce self-sustained theta-nested gamma oscillations, a single stimulation pulse could switch the network behavior from non-oscillatory to a state producing sustained oscillations. Next, we demonstrated that, for a weaker theta input, pulse train stimulation at the theta frequency could transiently restore seemingly physiological oscillations. Importantly, the presence of phase reset influenced whether these two effects depended on the phase at which stimulation onset was delivered, which has practical implications for designing neurostimulation protocols that are triggered by the phase of ongoing theta oscillations. This novel model opens new avenues for studying the effects of neurostimulation on the hippocampal formation. Furthermore, our hybrid approach that combines different levels of abstraction could be extended in future work to other neural circuits that produce dynamical brain rhythms.

A computational model to simulate spectral modulation and speech perception experiments of cochlear implant users

Speech understanding in cochlear implant (CI) users presents large intersubject variability that may be related to different aspects of the peripheral auditory system, such as the electrode–nerve interface and neural health conditions. This variability makes it more challenging to proof differences in performance between different CI sound coding strategies in regular clinical studies, nevertheless, computational models can be helpful to assess the speech performance of CI users in an environment where all these physiological aspects can be controlled. In this study, differences in performance between three variants of the HiRes Fidelity 120 (F120) sound coding strategy are studied with a computational model. The computational model consists of (i) a processing stage with the sound coding strategy, (ii) a three-dimensional electrode-nerve interface that accounts for auditory nerve fiber (ANF) degeneration, (iii) a population of phenomenological ANF models, and (iv) a feature extractor algorithm to obtain the internal representation (IR) of the neural activity. As the back-end, the simulation framework for auditory discrimination experiments (FADE) was chosen. Two experiments relevant to speech understanding were performed: one related to spectral modulation threshold (SMT), and the other one related to speech reception threshold (SRT). These experiments included three different neural health conditions (healthy ANFs, and moderate and severe ANF degeneration). The F120 was configured to use sequential stimulation (F120-S), and simultaneous stimulation with two (F120-P) and three (F120-T) simultaneously active channels. Simultaneous stimulation causes electric interaction that smears the spectrotemporal information transmitted to the ANFs, and it has been hypothesized to lead to even worse information transmission in poor neural health conditions. In general, worse neural health conditions led to worse predicted performance; nevertheless, the detriment was small compared to clinical data. Results in SRT experiments indicated that performance with simultaneous stimulation, especially F120-T, were more affected by neural degeneration than with sequential stimulation. Results in SMT experiments showed no significant difference in performance. Although the proposed model in its current state is able to perform SMT and SRT experiments, it is not reliable to predict real CI users' performance yet. Nevertheless, improvements related to the ANF model, feature extraction, and predictor algorithm are discussed.

A phenomenological computational model of the evoked action potential fitted to human cochlear implant responses

There is a growing interest in biomedical engineering in developing procedures that provide accurate simulations of the neural response to electrical stimulus produced by implants. Moreover, recent research focuses on models that take into account individual patient characteristics.

We present a phenomenological computational model that is customized with the patient’s data provided by the electrically evoked compound action potential (ECAP) for simulating the neural response to electrical stimulus produced by the electrodes of cochlear implants (CIs). The model links the input currents of the electrodes to the simulated ECAP.

Potentials and currents are calculated by solving the quasi-static approximation of the Maxwell equations with the finite element method (FEM). In ECAPs recording, an active electrode generates a current that elicits action potentials in the surrounding auditory nerve fibers (ANFs). The sum of these action potentials is registered by other nearby electrode. Our computational model emulates this phenomenon introducing a set of line current sources replacing the ANFs by a set of virtual neurons (VNs). To fit the ECAP amplitudes we assign a suitable weight to each VN related with the probability of an ANF to be excited. This probability is expressed by a cumulative beta distribution parameterized by two shape parameters that are calculated by means of a differential evolution algorithm (DE). Being the weights function of the current density, any change in the design of the CI affecting the current density produces changes in the weights and, therefore, in the simulated ECAP, which confers to our model a predictive capacity.

The results of the validation with ECAP data from two patients are presented, achieving a satisfactory fit of the experimental data with those provided by the proposed computational model.

The cochlea, found in the inner ear, is the organ where the sound is transformed into an electrical pulse to be transmitted by the neurons to the auditory cortex. Hearing loss can be caused by damage to the hair cells, in which case neuronal excitation is impaired. CIs are devices that replace the normal function of the impaired/damaged Organ of Corti. Computational models allow a better understanding of the mechanisms involved in the electrical stimulation of the auditory nerve. These models can help biomedical engineers to develop new CIs with improved auditory performance. One important aspect of our model is its customization with the patient’s data provided by the recording of the evoked compound action potential (the synchronous firing of a population of electrically stimulated auditory nerve fibers). This phenomenological model allows us to predict the registers of neural stimulation produced when the auditory nerve is stimulated with the CIs. We have validated the proposed model with real data obtained from two patients with CIs.

Recording Saltatory Conduction Along Sensory Axons Using a High-Density Microelectrode Array

Myelinated fibers are specialized neurological structures used for conducting action potentials quickly and reliably, thus assisting neural functions. Although demyelination leads to serious functional impairments, little is known the relationship between myelin structural change and increase in conduction velocity during myelination and demyelination processes. There are no appropriate methods for the long-term evaluation of spatial characteristics of saltatory conduction along myelinated axons. Herein, we aimed to detect saltatory conduction from the peripheral nervous system neurons using a high-density microelectrode array. Rat sensory neurons and intrinsic Schwann cells were cultured. Immunofluorescence and ultrastructure examination showed that the myelinating Schwann cells appeared at 1 month, and compact myelin was formed by 10 weeks in vitro. Activity of rat sensory neurons was evoked with optogenetic stimulation, and axon conduction was detected with high-density microelectrode arrays. Some conductions included high-speed segments with low signal amplitude. The same segment could be detected with electrical recording and immunofluorescent imaging for a myelin-related protein. The spatiotemporal analysis showed that some segments show a velocity of more than 2 m/s and that ends of the segments show a higher electrical sink, suggesting that saltatory conduction occurred in myelinated axons. Moreover, mathematical modeling supported that the recorded signal was in the appropriate range for axon and electrode sizes. Overall, our method could be a feasible tool for evaluating spatial characteristics of axon conduction including saltatory conduction, which is applicable for studying demyelination and remyelination.

Impedance Spectroscopy Dynamics of Biological Neural Elements: From Memristors to Neurons and Synapses

Understanding the operation of neurons and synapses is essential to reproducing biological computation. Building artificial neuromorphic networks opens the door to a new generation of faster and low-energy-consuming electronic circuits for computation. The main candidates to imitate the natural biocomputation processes, such as the generation of action potentials and spiking, are memristors. Generally, the study of the performance of material neuromorphic elements is done by the analysis of time transient signals. Here, we present an analysis of neural systems in the frequency domain by small-amplitude ac impedance spectroscopy. We start from the constitutive equations for the conductance and memory effect, and we derive and classify the impedance spectroscopy spectra. We first provide a general analysis of a memristor and demonstrate that this element can be expressed as a combination of simple parts. In particular, we derive a basic equivalent circuit where the memory effect is represented by an RL branch. We show that this ac model is quite general and describes the inductive/negative capacitance response in many systems such as halide perovskites and organic LEDs. Thereafter, we derive the impedance response of the integrate-and-fire exponential adaptative neuron model that introduces a negative differential resistance and a richer set of spectra. On the basis of these insights, we provide an interpretation of the varied spectra that appear in the more general Hodgkin-Huxley neuron model. Our work provides important criteria to determine the properties that must be found in material realizations of neuronal elements. This approach has the great advantage that the analysis of highly complex phenomena can be based purely on the shape of experimental impedance spectra, avoiding the need for specific modeling of rather involved material processes that produce the required response.

Effect of Channel Interaction on Vocal Cue Perception in Cochlear Implant Users

Speech intelligibility in multitalker settings is challenging for most cochlear implant (CI) users. One possibility for this limitation is the suboptimal representation of vocal cues in implant processing, such as the fundamental frequency (F0), and the vocal tract length (VTL). Previous studies suggested that while F0 perception depends on spectrotemporal cues, VTL perception relies largely on spectral cues. To investigate how spectral smearing in CIs affects vocal cue perception in speech-on-speech (SoS) settings, adjacent electrodes were simultaneously stimulated using current steering in 12 Advanced Bionics users to simulate channel interaction. In current steering, two adjacent electrodes are simultaneously stimulated forming a channel of parallel stimulation. Three such stimulation patterns were used: Sequential (one current steering channel), Paired (two channels), and Triplet stimulation (three channels). F0 and VTL just-noticeable differences (JNDs; Task 1), in addition to SoS intelligibility (Task 2) and comprehension (Task 3), were measured for each stimulation strategy. In Tasks 2 and 3, four maskers were used: the same female talker, a male voice obtained by manipulating both F0 and VTL (F0+VTL) of the original female speaker, a voice where only F0 was manipulated, and a voice where only VTL was manipulated. JNDs were measured relative to the original voice for the F0, VTL, and F0+VTL manipulations. When spectral smearing was increased from Sequential to Triplet, a significant deterioration in performance was observed for Tasks 1 and 2, with no differences between Sequential and Paired stimulation. Data from Task 3 were inconclusive. These results imply that CI users may tolerate certain amounts of channel interaction without significant reduction in performance on tasks relying on voice perception. This points to possibilities for using parallel stimulation in CIs for reducing power consumption.

Protonic conductor: better understanding neural resting and action potential

With the employment of the transmembrane electrostatic proton localization theory with a new membrane potential equation, neural resting and action potential is now much better understood as the voltage contributed by the localized protons/cations at a neural liquid- membrane interface. Accordingly, the neural resting/action potential is essentially a protonic/cationic membrane capacitor behavior. It is now understood with a newly formulated action potential equation: when action potential is <0 (negative number), the localized protons/cations charge density at the liquid-membrane interface along the periplasmic side is >0 (positive number); when the action potential is >0, the concentration of the localized protons and localized nonproton cations is <0, indicating a "depolarization" state. The nonlinear curve of the localized protons/cations charge density in the real-time domain of an action potential spike appears as an inverse mirror image to the action potential. The newly formulated action potential equation provides biophysical insights for neuron electrophysiology, which may represent a complementary development to the classic Goldman-Hodgkin-Katz equation. With the use of the action potential equation, the biological significance of axon myelination is now also elucidated as to provide protonic insulation and prevent any ions both inside and outside of the neuron from interfering with the action potential signal, so that the action potential can quickly propagate along the axon with minimal (e.g., 40 times less) energy requirement.NEW & NOTEWORTHY The newly formulated action potential equation provides biophysical insights for neuron electrophysiology, which may represent a complementary development to the classic Goldman-Hodgkin-Katz equation. The nonlinear curve of the localized protons/cations charge density in the real-time domain of an action potential spike appears as an inverse mirror image to the action potential. The biological significance of axon myelination is now elucidated as to provide protonic insulation and prevent any ions from interfering with action potential signal.