Developing a hippocampal neural prosthetic to facilitate human memory encoding and recall of stimulus features and categories

Objective

Here, we demonstrate the first successful use of static neural stimulation patterns for specific information content. These static patterns were derived by a model that was applied to a subject's own hippocampal spatiotemporal neural codes for memory.

Approach

We constructed a new model of processes by which the hippocampus encodes specific memory items via spatiotemporal firing of neural ensembles that underlie the successful encoding of targeted content into short-term memory. A memory decoding model (MDM) of hippocampal CA3 and CA1 neural firing was computed which derives a stimulation pattern for CA1 and CA3 neurons to be applied during the encoding (sample) phase of a delayed match-to-sample (DMS) human short-term memory task.

Main results

MDM electrical stimulation delivered to the CA1 and CA3 locations in the hippocampus during the sample phase of DMS trials facilitated memory of images from the DMS task during a delayed recognition (DR) task that also included control images that were not from the DMS task. Across all subjects, the stimulated trials exhibited significant changes in performance in 22.4% of patient and category combinations. Changes in performance were a combination of both increased memory performance and decreased memory performance, with increases in performance occurring at almost 2 to 1 relative to decreases in performance. Across patients with impaired memory that received bilateral stimulation, significant changes in over 37.9% of patient and category combinations was seen with the changes in memory performance show a ratio of increased to decreased performance of over 4 to 1. Modification of memory performance was dependent on whether memory function was intact or impaired, and if stimulation was applied bilaterally or unilaterally, with nearly all increase in performance seen in subjects with impaired memory receiving bilateral stimulation.

Significance

These results demonstrate that memory encoding in patients with impaired memory function can be facilitated for specific memory content, which offers a stimulation method for a future implantable neural prosthetic to improve human memory.

Patterned Hippocampal Stimulation Facilitates Memory in Patients With a History of Head Impact and/or Brain Injury

Rationale: Deep brain stimulation (DBS) of the hippocampus is proposed for enhancement of memory impaired by injury or disease. Many pre-clinical DBS paradigms can be addressed in epilepsy patients undergoing intracranial monitoring for seizure localization, since they already have electrodes implanted in brain areas of interest. Even though epilepsy is usually not a memory disorder targeted by DBS, the studies can nevertheless model other memory-impacting disorders, such as Traumatic Brain Injury (TBI). Methods: Human patients undergoing Phase II invasive monitoring for intractable epilepsy were implanted with depth electrodes capable of recording neurophysiological signals. Subjects performed a delayed-match-to-sample (DMS) memory task while hippocampal ensembles from CA1 and CA3 cell layers were recorded to estimate a multi-input, multi-output (MIMO) model of CA3-to-CA1 neural encoding and a memory decoding model (MDM) to decode memory information from CA3 and CA1 neuronal signals. After model estimation, subjects again performed the DMS task while either MIMO-based or MDM-based patterned stimulation was delivered to CA1 electrode sites during the encoding phase of the DMS trials. Each subject was sorted (post hoc) by prior experience of repeated and/or mild-to-moderate brain injury (RMBI), TBI, or no history (control) and scored for percentage successful delayed recognition (DR) recall on stimulated vs. non-stimulated DMS trials. The subject’s medical history was unknown to the experimenters until after individual subject memory retention results were scored. Results: When examined compared to control subjects, both TBI and RMBI subjects showed increased memory retention in response to both MIMO and MDM-based hippocampal stimulation. Furthermore, effects of stimulation were also greater in subjects who were evaluated as having pre-existing mild-to-moderate memory impairment. Conclusion: These results show that hippocampal stimulation for memory facilitation was more beneficial for subjects who had previously suffered a brain injury (other than epilepsy), compared to control (epilepsy) subjects who had not suffered a brain injury. This study demonstrates that the epilepsy/intracranial recording model can be extended to test the ability of DBS to restore memory function in subjects who previously suffered a brain injury other than epilepsy, and support further investigation into the beneficial effect of DBS in TBI patients.

Altered adult neurogenesis and gliogenesis in patients with mesial temporal lobe epilepsy

The hippocampus is the most common seizure focus in people. In the hippocampus, aberrant neurogenesis plays a critical role in the initiation and progression of epilepsy in rodent models, but it is unknown whether this also holds true in humans. To address this question, we used immunofluorescence on control healthy hippocampus and surgical resections from mesial temporal lobe epilepsy (MTLE), plus neural stem-cell cultures and multi-electrode recordings of ex vivo hippocampal slices. We found that a longer duration of epilepsy is associated with a sharp decline in neuronal production and persistent numbers in astrogenesis. Further, immature neurons in MTLE are mostly inactive, and are not observed in cases with local epileptiform-like activity. However, immature astroglia are present in every MTLE case and their location and activity are dependent on epileptiform-like activity. Immature astroglia, rather than newborn neurons, therefore represent a potential target to continually modulate adult human neuronal hyperactivity.

A Double-Layer Multi-Resolution Classification Model for Decoding Spatiotemporal Patterns of Spikes With Small Sample Size

We build a double-layer, multiple temporal-resolution classification model for decoding single-trial spatiotemporal patterns of spikes. The model takes spiking activities as input signals and binary behavioral or cognitive variables as output signals and represents the input-output mapping with a double-layer ensemble classifier. In the first layer, to solve the underdetermined problem caused by the small sample size and the very high dimensionality of input signals, B-spline functional expansion and L1-regularized logistic classifiers are used to reduce dimensionality and yield sparse model estimations. A wide range of temporal resolutions of neural features is included by using a large number of classifiers with different numbers of B-spline knots. Each classifier serves as a base learner to classify spatiotemporal patterns into the probability of the output label with a single temporal resolution. A bootstrap aggregating strategy is used to reduce the estimation variances of these classifiers. In the second layer, another L1-regularized logistic classifier takes outputs of first-layer classifiers as inputs to generate the final output predictions. This classifier serves as a meta-learner that fuses multiple temporal resolutions to classify spatiotemporal patterns of spikes into binary output labels. We test this decoding model with both synthetic and experimental data recorded from rats and human subjects performing memory-dependent behavioral tasks. Results show that this method can effectively avoid overfitting and yield accurate prediction of output labels with small sample size. The double-layer, multi-resolution classifier consistently outperforms the best single-layer, single-resolution classifier by extracting and utilizing multi-resolution spatiotemporal features of spike patterns in the classification.

Bridging Hierarchies in Multi-Scale Models of Neural Systems: Look-Up Tables Enable Computationally Efficient Simulations of Non-linear Synaptic Dynamics

Synapses are critical actors of neuronal transmission as they form the basis of chemical communication between neurons. Accurate computational models of synaptic dynamics may prove important in elucidating emergent properties across hierarchical scales. Yet, in large-scale neuronal network simulations, synapses are often modeled as highly simplified linear exponential functions due to their small computational footprint. However, these models cannot capture the complex non-linear dynamics that biological synapses exhibit and thus, are insufficient in representing synaptic behavior accurately. Existing detailed mechanistic synapse models can replicate these non-linear dynamics by modeling the underlying kinetics of biological synapses, but their high complexity prevents them from being a suitable option in large-scale models due to long simulation times. This motivates the development of more parsimonious models that can capture the complex non-linear dynamics of synapses accurately while maintaining a minimal computational cost. We propose a look-up table approach that stores precomputed values thereby circumventing most computations at runtime and enabling extremely fast simulations for glutamatergic receptors AMPAr and NMDAr. Our results demonstrate that this methodology is capable of replicating the dynamics of biological synapses as accurately as the mechanistic synapse models while offering up to a 56-fold increase in speed. This powerful approach allows for multi-scale neuronal networks to be simulated at large scales, enabling the investigation of how low-level synaptic activity may lead to changes in high-level phenomena, such as memory and learning.

A sparse multiscale nonlinear autoregressive model for seizure prediction

Objectives.Accurate seizure prediction is highly desirable for medical interventions such as responsive electrical stimulation. We aim to develop a classification model that can predict seizures by identifying preictal states, i.e. the precursor of a seizure, based on multi-channel intracranial electroencephalography (iEEG) signals.Approach.A two-level sparse multiscale classification model was developed to classify interictal and preictal states from iEEG data. In the first level, short time-scale linear dynamical features were extracted as autoregressive (AR) model coefficients; arbitrary (usually long) time-scale linear and nonlinear dynamical features were extracted as Laguerre-Volterra AR model coefficients; root-mean-square error of model prediction was used as a feature representing model unpredictability. In the second level, all features were fed into a sparse classifier to discriminate the iEEG data between interictal and preictal states.Main results. The two-level model can accurately classify seizure states using iEEG data recorded from ten canine and human subjects. Adding arbitrary (usually long) time-scale and nonlinear features significantly improves model performance compared with the conventional AR modeling approach. There is a high degree of variability in the types of features contributing to seizure prediction across different subjects.Significance. This study suggests that seizure generation may involve distinct linear/nonlinear dynamical processes caused by different underlying neurobiological mechanisms. It is necessary to build patient-specific classification models with a wide range of dynamical features.

Topographic Organization of Correlation Along the Longitudinal and Transverse Axes in Rat Hippocampal CA3 Due to Excitatory Afferents

The topographic organization of afferents to the hippocampal CA3 subfield are well-studied, but their role in influencing the spatiotemporal dynamics of population activity is not understood. Using a large-scale, computational neuronal network model of the entorhinal-dentate-CA3 system, the effects of the perforant path, mossy fibers, and associational system on the propagation and transformation of network spiking patterns were investigated. A correlation map was constructed to characterize the spatial structure and temporal evolution of pairwise correlations which underlie the emergent patterns found in the population activity. The topographic organization of the associational system gave rise to changes in the spatial correlation structure along the longitudinal and transverse axes of the CA3. The resulting gradients may provide a basis for the known functional organization observed in hippocampus.

A Computational Model of Mitochondria Motility in Axons

Mitochondria play a critical role in regulating cellular processes including ATP production, intracellular calcium signaling and generation of reactive oxidative species (ROS). Neurons rely on mitochondrial function to perform a range of complex processes, and mitochondrial dysfunctions have been shown to have an impact in pathologies of the nervous system. Yet, neurons contain a finite number of mitochondria, and their location is known to change in response to a number of factors including age and cellular activity, thereby impacting neuronal response. In this paper, we introduce a novel computational model of mitochondria motility that focuses on their movements along the axon. We describe the biological processes involved and the main parameters of the model. We use the model to investigate how some of these parameters affect the ability of mitochondria to position themselves in regions of high energy demand. Finally, we discuss the significance of our work and its downstream applications in further understanding pathologies of the nervous system such as Alzheimer's disease, and help identify potential novel therapeutic targets.

A Computational Model of the Cholinergic Modulation of CA1 Pyramidal Cell Activity

Dysfunction in cholinergic modulation has been linked to a variety of cognitive disorders including Alzheimer's disease. The important role of this neurotransmitter has been explored in a variety of experiments, yet many questions remain unanswered about the contribution of cholinergic modulation to healthy hippocampal function. To address this question, we have developed a model of CA1 pyramidal neuron that takes into consideration muscarinic receptor activation in response to changes in extracellular concentration of acetylcholine and its effects on cellular excitability and downstream intracellular calcium dynamics. This model incorporates a variety of molecular agents to accurately simulate several processes heretofore ignored in computational modeling of CA1 pyramidal neurons. These processes include the inhibition of ionic channels by phospholipid depletion along with the release of calcium from intracellular stores (i.e., the endoplasmic reticulum). This paper describes the model and the methods used to calibrate its behavior to match experimental results. The result of this work is a compartmental model with calibrated mechanisms for simulating the intracellular calcium dynamics of CA1 pyramidal cells with a focus on those related to release from calcium stores in the endoplasmic reticulum. From this model we also make various predictions for how the inhibitory and excitatory responses to cholinergic modulation vary with agonist concentration. This model expands the capabilities of CA1 pyramidal cell models through the explicit modeling of molecular interactions involved in healthy cognitive function and disease. Through this expanded model we come closer to simulating these diseases and gaining the knowledge required to develop novel treatments.

Admittance Method for Estimating Local Field Potentials Generated in a Multi-Scale Neuron Model of the Hippocampus

Significant progress has been made toward model-based prediction of neral tissue activation in response to extracellular electrical stimulation, but challenges remain in the accurate and efficient estimation of distributed local field potentials (LFP). Analytical methods of estimating electric fields are a first-order approximation that may be suitable for model validation, but they are computationally expensive and cannot accurately capture boundary conditions in heterogeneous tissue. While there are many appropriate numerical methods of solving electric fields in neural tissue models, there isn't an established standard for mesh geometry nor a well-known rule for handling any mismatch in spatial resolution. Moreover, the challenge of misalignment between current sources and mesh nodes in a finite-element or resistor-network method volume conduction model needs to be further investigated. Therefore, using a previously published and validated multi-scale model of the hippocampus, the authors have formulated an algorithm for LFP estimation, and by extension, bidirectional communication between discretized and numerically solved volume conduction models and biologically detailed neural circuit models constructed in NEURON. Development of this algorithm required that we assess meshes of (i) unstructured tetrahedral and grid-based hexahedral geometries as well as (ii) differing approaches for managing the spatial misalignment of current sources and mesh nodes. The resulting algorithm is validated through the comparison of Admittance Method predicted evoked potentials with analytically estimated LFPs. Establishing this method is a critical step toward closed-loop integration of volume conductor and NEURON models that could lead to substantial improvement of the predictive power of multi-scale stimulation models of cortical tissue. These models may be used to deepen our understanding of hippocampal pathologies and the identification of efficacious electroceutical treatments.

Accelerating Estimation of a Multi-Input Multi-Output Model of the Hippocampus with a Parallel Computing Strategy

To build hippocampal memory prosthesis for restoring memory functions, we previously developed and implemented a multi-input multi-output (MIMO) nonlinear dynamic model of the hippocampus. This model can successfully predict hippocampal output spike activities based on input spike activities, and thus be used to drive microstimulation to bypass the damaged hippocampal region. Building such a MIMO model involves estimations of a large number of model coefficients, which typically takes hundreds of hours using a single personal computer. In practice, however, due to the requirement of medical care and clinical trials, the modeling processes must be completed within 72 hours after the recording, so that models can be used to drive stimulations. To solve this problem, we utilized a parallelization strategy to divide the whole MIMO model computation involving iterative estimation and optimization into independent computing tasks that can be performed simultaneously in multiple computer nodes. Such a strategy was implemented on the high-performance computing cluster at the University of Southern California. It reduced the model estimation time to tens of hours and thus allowed us to complete the modeling process within the required time frame to further test model-driven electrical stimulation for the hippocampal memory prosthesis.

Generation of Granule Cell Dendritic Morphologies by Estimating the Spatial Heterogeneity of Dendritic Branching

Biological realism of dendritic morphologies is important for simulating electrical stimulation of brain tissue. By adding point process modeling and conditional sampling to existing generation strategies, we provide a novel means of reproducing the nuanced branching behavior that occurs in different layers of granule cell dendritic morphologies. In this study, a heterogeneous Poisson point process was used to simulate branching events. Conditional distributions were then used to select branch angles depending on the orthogonal distance to the somatic plane. The proposed method was compared to an existing generation tool and a control version of the proposed method that used a homogeneous Poisson point process. Morphologies were generated with each method and then compared to a set of digitally reconstructed neurons. The introduction of a conditionally dependent branching rate resulted in the generation of morphologies that more accurately reproduced the emergent properties of dendritic material per layer, Sholl intersections, and proximal passive current flow. Conditional dependence was critically important for the generation of realistic granule cell dendritic morphologies.

Pathogenic Processes Underlying Alzheimer's Disease: Modeling the Effects of Amyloid Beta on Synaptic Transmission

The molecular mechanisms underlying Alzheimer's disease (AD) have been and are still under heavy scrutiny to better understand what leads to the onset and progression of the disease, and to design and develop efficacious therapeutic strategies. These decade-long studies have taught us a lot regarding the various molecular pathways involved in the pathology, but a complete dynamic picture of the underlying pathological mechanisms is still missing.We propose to provide a technological answer to fill this gap by developing and using a computational approach that integrates AD-related experimental findings and their effects on multiple aspects of neuronal function. The present study focuses on implementing one known pathogenic process: the binding of amyloid beta, the hallmark of AD, on NMDA receptors, receptors present in the main type of excitatory synapses in the brain, thereby affecting synaptic transmission and downstream pathways. We describe model implementation and calibration; we then quantify the downstream effects of this disruption both in terms of electrical activity (changes in short-term spiking activity of the postsynaptic neuron), and biochemical pathways activation through changes in calcium dynamics (an important trigger to longer-term changes). The computational approach outlined constitutes an insightful instrument to examine the downstream consequences of multiple pathogenic dysfunctions on higher level observables and sets the path for in-silico discovery and testing of therapeutic agents.

ROOTS: An Algorithm to Generate Biologically Realistic Cortical Axons and an Application to Electroceutical Modeling

Advances in computation and neuronal modeling have enabled the study of entire neural tissue systems with an impressive degree of biological realism. These efforts have focused largely on modeling dendrites and somas while largely neglecting axons. The need for biologically realistic explicit axonal models is particularly clear for applications involving clinical and therapeutic electrical stimulation because axons are generally more excitable than other neuroanatomical subunits. While many modeling efforts can rely on existing repositories of reconstructed dendritic/somatic morphologies to study real cells or to estimate parameters for a generative model, such datasets for axons are scarce and incomplete. Those that do exist may still be insufficient to build accurate models because the increased geometric variability of axons demands a proportional increase in data. To address this need, a Ruled-Optimum Ordered Tree System (ROOTS) was developed that extends the capability of neuronal morphology generative methods to include highly branched cortical axon terminal arbors. Further, this study presents and explores a clear use-case for such models in the prediction of cortical tissue response to externally applied electric fields. The results presented herein comprise (i) a quantitative and qualitative analysis of the generative algorithm proposed, (ii) a comparison of generated fibers with those observed in histological studies, (iii) a study of the requisite spatial and morphological complexity of axonal arbors for accurate prediction of neuronal response to extracellular electrical stimulation, and (iv) an extracellular electrical stimulation strength-duration analysis to explore probable thresholds of excitation of the dentate perforant path under controlled conditions. ROOTS demonstrates a superior ability to capture biological realism in model fibers, allowing improved accuracy in predicting the impact that microscale structures and branching patterns have on spatiotemporal patterns of activity in the presence of extracellular electric fields.

Network Activity Due to Topographic Organization of Schaffer Collaterals in a Large-Scale Model of Rat CA1

Connectivity between neural regions, particularly in the hippocampus, is seldom all-to-all or random, yet it is the predominant method by which connectivity is implemented in most models of neuronal networks. We have been developing a computational platform for simulating the trisynaptic circuit of rat hippocampus with which we have constructed a large-scale, biologically-realistic, spiking neuronal network model of the entorhinal-dentate-CA3 system. Using the model, we had demonstrated a non-trivial effect of topographic connectivity on network dynamics and function. In this work, we detail the introduction of the CA1 subregion to the large-scale model. Using anatomical data, we constrained the distribution of axon collaterals, i.e., Schaffer collaterals, projected from CA3 to CA1 and preserved the topographic organization of the projections. Using a simplified multi-compartmental model of CA1 pyramidal cells and a single compartment model of CA1 parvalbumin basket cells, that were connected with disynaptic feedforward inhibition and feedback inhibition, we demonstrate the network activity of the CA1 network given a topographic organization of Schaffer collaterals. From this introduction of CA1 to the large-scale model, we can then observe the successive transformation of spatio-temporal, spiking neural activity as it propagates through the trisynaptic circuit.

Decoding Position to Analyze Spatial Information Encoding in a Large-Scale Neuronal Network Model of Rat Dentate Gyrus

Spatial information is encoded by the hippocampus, and the factors that contribute to the amount of information that can be encoded and the transformation of spatial information through the trisynaptic circuit remain an important issue. A large-scale neuronal network model of the rat entorhinal-dentate system was developed with multicompartmental representations of the neurons within the dentate gyrus. Spatial information was introduced to the network via grid cell activity, and the spatial information encoding capabilities of the network were assessed using a recursive decoding algorithm to estimate the position of a virtual rat using the dentate activity. To obtain a measure for the information that the network could convey, decoding error was calculated for different decoding population sizes. Decoding error decreased exponentially as a function of population size. Therefore, the time constant and the asymptote of the error curve could be used as metrics to compare the changes in encoding performance. In conjunction with the large-scale model, this paradigm can be used to characterize how neural properties, network composition, and the interactions between different subfields affect spatial information encoding.

Multi-Input, Multi-Output Neuronal Mode Network Approach to Modeling the Encoding Dynamics and Functional Connectivity of Neural Systems

This letter proposes a novel method, multi-input, multi-output neuronal mode network (MIMO-NMN), for modeling encoding dynamics and functional connectivity in neural ensembles such as the hippocampus. Compared with conventional approaches such as the Volterra-Wiener model, linear-nonlinear-cascade (LNC) model, and generalized linear model (GLM), the NMN has several advantages in terms of estimation accuracy, model interpretation, and functional connectivity analysis. We point out the limitations of current neural spike modeling methods, especially the estimation biases caused by the imbalanced class problem when the number of zeros is significantly larger than ones in the spike data. We use synthetic data to test the performance of NMN with a comparison of the traditional methods, and the results indicate the NMN approach could reduce the imbalanced class problem and achieve better predictions. Subsequently, we apply the MIMO-NMN method to analyze data from the human hippocampus. The results indicate that the MIMO-NMN method is a promising approach to modeling neural dynamics and analyzing functional connectivity of multi-neuronal data.

Axonal Anatomy Optimizes Spatial Encoding in the Rat Entorhinal-Dentate System: A Computational Study

OBJECTIVE

The network architecture connecting neural regions is defined by the organization and anatomical properties of the projecting axons, but its contributions to neural encoding and system function are difficult to study experimentally.

METHODS

Using a large-scale, spiking neuronal network model of rat dentate gyrus, the role of the anatomy of the entorhinal-dentate axonal projection was evaluated in the context of spatial encoding by incorporating grid cell activity to provide physiological, spatially-correlated input. The dorso-ventral extents of the entorhinal axon terminal fields were varied to generate different feedforward architectures, and the resulting spatial representations and spatial information scores of the network were evaluated. Position was decoded from the population activity using a point process filter to investigate the contributions of network architecture on spatial encoding.

RESULTS

The model predicted the emergence of anatomical gradients within the dentate gyrus for place field size and spatial information along its dorso-ventral axis, which were dependent on the extents of the entorhinal axon terminal fields. The decoding results revealed an optimal performance at an axon terminal field extent of 2 mm that lies within the biological range.

CONCLUSION

The axonal anatomy mediates a tradeoff between encoding multiple place field sizes or achieving a high spatial information score, and the combination of both properties is necessary to maximize spatial encoding by a network.

SIGNIFICANCE

In total, this paper establishes a mechanistic neuronal network model that, in concert with information-theoretic and statistical methods, can be used to investigate how lower level properties contribute to higher level function.

A Closed-Loop Multi-Channel Asynchronous Neurostimulator to Mimic Neural Code for Cognitive Prosthesis

We describe a novel hardware and embedded system design of a closed loop neurostimulator for generating precise neural code-like, multi-channel, asynchronous electrical stimulation pulses. Such stimulator will be used as the output unit of the cortical prosthesis that aims to restore cognitive functions by reinstating the neural signal transmission.

Cholinergic Modulation of CA1 Pyramidal Cells via M1 Muscarinic Receptor Activation: A Computational Study at Physiological and Supraphysiological Levels

The hippocampus receives extensive cholinergic modulation from the basal forebrain, which has been shown to have a prominent role in attention, learning, and synaptic plasticity. Disruptions of this modulation have been linked to a variety of neural disorders including Alzheimer's Disease. Pyramidal cells of the CA1 region of the hippocampus express several cholinergic receptor types in different locations throughout the cells' morphology. Developing a computational model of these cells and their modulation provides a unique opportunity to explore how each receptor type alters the overall computational role of the cell. To this end we implemented a kinetic model of the most widely distributed receptor type, the M1 muscarinic receptor and examined its role on excitation of a compartmental model of a CA1 pyramidal cell. We demonstrate that the proposed model replicates the increased pyramidal cell excitability seen in experimental results. We then used the model to replicate the effect of organophosphates, a class of pesticides and chemical weapons, whose effects consist in inhibiting the hydrolysis of acetylcholine; we demonstrated the effect of increasing concentrations of acetylcholine on the pyramidal cell's excitability. The cell model we implemented and its associated modulation constitute a basis for exploring the effects of cholinergic modulation in a large scale network model of the hippocampus both under physiological and supraphysiological levels.

Point Process Filtering Estimates of Branching Rate for Neural Dendritic Morphology Generation

Current parametric approaches to dendritic morphology generation are limited in their ability to replicate realistic branching. A non-parametric approach applying a point process filter and the expectation-maximization algorithm offers a data-based solution that estimates the dendritic branching rate based on observations of bifurcation events in real neurons. Point processes can then be simulated using this branching rate estimate to indicate when a generated morphology should branch. Morphologies generated using this technique match both basic and emergent property distributions of the real neurons used as input into the algorithm. Further refinement of branching angles will allow for a flexible tool to generate realistic morphologies of a variety of neuronal stereotypes.

Chronic multi-region recording from the rat hippocampus in vivo with a flexible Parylene-based multi-electrode array

Neural activities of free-moving animals provide valuable insights into behavior, memory formation and cognitive function of the hippocampus. Unitary activities simultaneously recorded from multiple sub-regions of the hippocampus enable detailed study of hippocampal neural circuits, but require high fidelity recordings with high temporal and spatial resolution. In this work, we explored the possibility of using Parylene-C as the structural material for a penetrating, multi-electrode array designed to record from multiple sub-region of the rat hippocampus. A 64-channel Parylene-based flexible electrode array was designed and fabricated. The layout of the electrode array was arranged to conform to the shape of cell body layers of the rat hippocampus. An insertion technique of temporarily reduce the effective length of the probe with polyethylene glycol (PEG) was developed and tested in vivo. The multi-electrode array was implanted into a rat hippocampus for chronic experimentation and unitary activities were collected both during the implantation and after recovery while the animal ran freely in an open field. Unitary activities with an average signal to noise ratios (SNR) of 3 to 4 were recorded with the Parylene probe over the period of one month after implantation.

A Glutamatergic Spine Model to Enable Multi-Scale Modeling of Nonlinear Calcium Dynamics

In synapses, calcium is required for modulating synaptic transmission, plasticity, synaptogenesis, and synaptic pruning. The regulation of calcium dynamics within neurons involves cellular mechanisms such as synaptically activated channels and pumps, calcium buffers, and calcium sequestrating organelles. Many experimental studies tend to focus on only one or a small number of these mechanisms, as technical limitations make it difficult to observe all features at once. Computational modeling enables incorporation of many of these properties together, allowing for more complete and integrated studies. However, the scale of existing detailed models is often limited to synaptic and dendritic compartments as the computational burden rapidly increases when these models are integrated in cellular or network level simulations. In this article we present a computational model of calcium dynamics at the postsynaptic spine of a CA1 pyramidal neuron, as well as a methodology that enables its implementation in multi-scale, large-scale simulations. We first present a mechanistic model that includes individually validated models of various components involved in the regulation of calcium at the spine. We validated our mechanistic model by comparing simulated calcium levels to experimental data found in the literature. We performed additional simulations with the mechanistic model to determine how the simulated calcium activity varies with respect to presynaptic-postsynaptic stimulation intervals and spine distance from the soma. We then developed an input-output (IO) model that complements the mechanistic calcium model and provide a computationally efficient representation for use in larger scale modeling studies; we show the performance of the IO model compared to the mechanistic model in terms of accuracy and speed. The models presented here help achieve two objectives. First, the mechanistic model provides a comprehensive platform to describe spine calcium dynamics based on individual contributing factors. Second, the IO model is trained on the main dynamical features of the mechanistic model and enables nonlinear spine calcium modeling on the cell and network level simulation scales. Utilizing both model representations provide a multi-level perspective on calcium dynamics, originating from the molecular interactions at spines and propagating the effects to higher levels of activity involved in network behavior.

Application of Parylene-Based Flexible Multi-Electrode Array for Recording From Subcortical Brain Regions From Behaving Rats

Obtaining multiple single-unit recordings in particular neural networks from behaving animals is crucial for the understanding of cognitive functions of the brain. Attaining stable, chronic recordings from the brain is also the foundation to develop effective cortical prosthetic devices. However, severe immune response caused by micromotion between stiff implants and surrounding brain tissue often limits the lifetime of penetrating, neural recording devices. To reduce the stiffness mismatch between recording devices and brain tissue, we developed a flexible, polymer based multi-electrode array for recording single neuron activities from the rat hippocampus, a major subcortical structure of the rat brain. Parylene C, a biocompatible polymer, was used as the structural and insulation material of the multi-electrode array. 64 platinum (Pt) recording electrodes were placed in groups along each shank to conform to the anatomical distribution of hippocampal principle neurons. The multi-electrode array was chronically implanted in three animals. After recovery, neural activity together with movement traces were collected from the behaving animals.

Graph-Based Models of Cortical Axons for the Prediction of Neuronal Response to Extracellular Electrical Stimulation

Over the past decade, many important insights to brain function have been obtained through clever application of detailed compartmental model neurons. New computing capabilities brought opportunities to study large networks of model neurons. Certain applications for these models, such as extracellular electrical stimulation, demand a very high degree of biological realism. While dendrites and somatic morphology may be obtained from explicit reconstructions, this approach is less useful for axonal structures, which are more difficult to characterize across a neuronal population. The purpose of this paper is to extend neuronal morphology generative models to highly branched axon terminal arbors as well as to present a clear use-case for such models in the study of cortical tissue response to externally applied electric fields. The results of this work are (i) presentation and quantitative/qualitative description of generated fibers and (ii) an extracellular electrical stimulation strength-duration study.

ASIC Implementation of a Nonlinear Dynamical Model for Hippocampal Prosthesis

A hippocampal prosthesis is a very large scale integration (VLSI) biochip that needs to be implanted in the biological brain to solve a cognitive dysfunction. In this letter, we propose a novel low-complexity, small-area, and low-power programmable hippocampal neural network application-specific integrated circuit (ASIC) for a hippocampal prosthesis. It is based on the nonlinear dynamical model of the hippocampus: namely multi-input, multi-output (MIMO)-generalized Laguerre-Volterra model (GLVM). It can realize the real-time prediction of hippocampal neural activity. New hardware architecture, a storage space configuration scheme, low-power convolution, and gaussian random number generator modules are proposed. The ASIC is fabricated in 40 nm technology with a core area of 0.122 mm[Formula: see text] and test power of 84.4 [Formula: see text]W. Compared with the design based on the traditional architecture, experimental results show that the core area of the chip is reduced by 84.94% and the core power is reduced by 24.30%.

Developing a hippocampal neural prosthetic to facilitate human memory encoding and recall

OBJECTIVE

We demonstrate here the first successful implementation in humans of a proof-of-concept system for restoring and improving memory function via facilitation of memory encoding using the patient's own hippocampal spatiotemporal neural codes for memory. Memory in humans is subject to disruption by drugs, disease and brain injury, yet previous attempts to restore or rescue memory function in humans typically involved only nonspecific, modulation of brain areas and neural systems related to memory retrieval.

APPROACH

We have constructed a model of processes by which the hippocampus encodes memory items via spatiotemporal firing of neural ensembles that underlie the successful encoding of short-term memory. A nonlinear multi-input, multi-output (MIMO) model of hippocampal CA3 and CA1 neural firing is computed that predicts activation patterns of CA1 neurons during the encoding (sample) phase of a delayed match-to-sample (DMS) human short-term memory task.

MAIN RESULTS

MIMO model-derived electrical stimulation delivered to the same CA1 locations during the sample phase of DMS trials facilitated short-term/working memory by 37% during the task. Longer term memory retention was also tested in the same human subjects with a delayed recognition (DR) task that utilized images from the DMS task, along with images that were not from the task. Across the subjects, the stimulated trials exhibited significant improvement (35%) in both short-term and long-term retention of visual information.

SIGNIFICANCE

These results demonstrate the facilitation of memory encoding which is an important feature for the construction of an implantable neural prosthetic to improve human memory.

Designing Patient-Specific Optimal Neurostimulation Patterns for Seizure Suppression

Neurostimulation is a promising therapy for abating epileptic seizures. However, it is extremely difficult to identify optimal stimulation patterns experimentally. In this study, human recordings are used to develop a functional 24 neuron network statistical model of hippocampal connectivity and dynamics. Spontaneous seizure-like activity is induced in silico in this reconstructed neuronal network. The network is then used as a testbed to design and validate a wide range of neurostimulation patterns. Commonly used periodic trains were not able to permanently abate seizures at any frequency. A simulated annealing global optimization algorithm was then used to identify an optimal stimulation pattern, which successfully abated 92% of seizures. Finally, in a fully responsive, or closed-loop, neurostimulation paradigm, the optimal stimulation successfully prevented the network from entering the seizure state. We propose that the framework presented here for algorithmically identifying patient-specific neurostimulation patterns can greatly increase the efficacy of neurostimulation devices for seizures.

Model-Based Analysis of Electrode Placement and Pulse Amplitude for Hippocampal Stimulation

Objective:

The ideal form of a neural-interfacing device is highly dependent upon the anatomy of the region with which it is meant to interface. Multiple-electrode arrays provide a system which can be adapted to various neural geometries. Computational models of stimulating systems have proven useful for evaluating electrode placement and stimulation protocols, but have yet to be adequately adapted to the unique features of the hippocampus.

Methods:

As an approach to understanding potential memory restorative devices, an Admittance Method-NEURON model was constructed to predict the direct and synaptic response of a region of the rat dentate gyrus to electrical stimulation of the perforant path.

Results:

A validation of estimated local field potentials against experimental recordings is performed and results of a bi-linear electrode placement and stimulation amplitude parameter search are presented.

Conclusion:

The parametric analysis presented herein suggests that stimulating electrodes placed between the lateral and medial perforant path, near the crest of the dentate gyrus, yield a larger relative population response to given stimuli.

Significance:

Beyond deepening understanding of the hippocampal tissue system, establishment of this model provides a method to evaluate candidate stimulating devices and protocols.

Evolving Applications, Technological Challenges and Future Opportunities in Neuromodulation: Proceedings of the Fifth Annual Deep Brain Stimulation Think Tank

The annual Deep Brain Stimulation (DBS) Think Tank provides a focal opportunity for a multidisciplinary ensemble of experts in the field of neuromodulation to discuss advancements and forthcoming opportunities and challenges in the field. The proceedings of the fifth Think Tank summarize progress in neuromodulation neurotechnology and techniques for the treatment of a range of neuropsychiatric conditions including Parkinson's disease, dystonia, essential tremor, Tourette syndrome, obsessive compulsive disorder, epilepsy and cognitive, and motor disorders. Each section of this overview of the meeting provides insight to the critical elements of discussion, current challenges, and identified future directions of scientific and technological development and application. The report addresses key issues in developing, and emphasizes major innovations that have occurred during the past year. Specifically, this year's meeting focused on technical developments in DBS, design considerations for DBS electrodes, improved sensors, neuronal signal processing, advancements in development and uses of responsive DBS (closed-loop systems), updates on National Institutes of Health and DARPA DBS programs of the BRAIN initiative, and neuroethical and policy issues arising in and from DBS research and applications in practice.

Multi-resolution multi-trial sparse classification model for decoding visual memories from hippocampal spikes in human

To understand how memories are encoded in the hippocampus, we build memory decoding models to classify visual memories based on hippocampal activities in human. Model inputs are spatio-temporal patterns of spikes recorded in the hippocampal CA3 and CA1 regions of epilepsy patients performing a delayed match-to-sample (DMS) task. Model outputs are binary labels indicating categories and features of sample images. To solve the super high-dimensional estimation problem with short data length, we develop a multi-trial, sparse model estimation method utilizing B-spline basis functions with a large range of temporal resolutions and a regularized logistic classifier. Results show that this model can effectively avoid overfitting and provide significant amount of prediction to memory categories and features using very limited number of data points. Stable estimation of sparse classification function matrices for each label can be obtained with this multi-resolution, multi-trial procedure. These classification models can be used not only to predict memory contents, but also to design optimal spatio-temporal patterns for eliciting specific memories in the hippocampus, and thus have important implications to the development of hippocampal memory prostheses.

Mechanism-Based and Input-Output Modeling of the Key Neuronal Connections and Signal Transformations in the CA3-CA1 Regions of the Hippocampus

This letter examines the results of input-output (nonparametric) modeling based on the analysis of data generated by a mechanism-based (parametric) model of CA3-CA1 neuronal connections in the hippocampus. The motivation is to obtain biological insight into the interpretation of such input-output (Volterra-equivalent) models estimated from synthetic data. The insights obtained may be subsequently used to interpretat input-output models extracted from actual experimental data. Specifically, we found that a simplified parametric model may serve as a useful tool to study the signal transformations in the hippocampal CA3-CA1 regions. Input-output modeling of model-based synthetic data show that GABAergic interneurons are responsible for regulating neuronal excitation, controlling the precision of spike timing, and maintaining network oscillations, in a manner consistent with previous studies. The input-output model obtained from real data exhibits intriguing similarities with its synthetic-data counterpart, demonstrating the importance of a dynamic resonance in the system/model response around 2 Hz to 3 Hz. Using the input-output model from real data as a guide, we may be able to amend the parametric model by incorporating more mechanisms in order to yield better-matching input-output model. The approach we present can also be applied to the study of other neural systems and pathways.

Place field detection using grid-based clustering in a large-scale computational model of the rat dentate gyrus

Place cells are neurons in the hippocampus that are sensitive to location within an environment. Simulations of a large-scale, computational model of the rat dentate gyrus using grid cell input have been performed resulting in granule cells that express multiple place fields. The typical method of detecting place fields using a global threshold on this data is unreliable as the characteristics of the place fields from a single neuron can be highly variable. A grid-based implementation of DENCLUE has been developed to calculate local thresholds to identify each place field. An adaptive binning algorithm used to smooth the rate maps was combined with the DENCLUE implementation to adaptively choose the size of the smoothing kernel and reduce the number of free parameters of the total algorithm. A sensitivity analysis was performed using the threshold parameter to demonstrate the robustness of using local thresholds as opposed to using a single global threshold in detecting the place fields resulting from the large-scale simulation. The analysis supports the use of applying local thresholds for place field detection and will be used to further investigate the role of granule cells in hippocampal function.

Development of a detailed model of calcium dynamics at the postsynaptic spine of an excitatory synapse

Postsynaptic calcium dynamics play a critical role in synaptic plasticity, but are often difficult to measure in experimental protocols due to their relatively fast rise and decay times, and the small spine dimensions. To circumvent these limitations, we propose to develop a computational model of calcium dynamics in the postsynaptic spine. This model integrates the main elements that participate in calcium concentration influx, efflux, diffusion and buffering. These consist of (i) spine geometry; (ii) calcium influx through NMDA receptors and voltage-dependent calcium channels (VDCC); (iii) calcium efflux with plasma membrane calcium pumps (PMCA) and sodium-calcium exchangers (NCX); (iv) intracellular calcium stores; and (v) calcium buffers. We herein present computational results we obtained and compare them with experimentally measured data, thereby validating the proposed model. Overall the development of such postsynaptic calcium model may help us better understand the intricacies of interplay between the different elements that shape calcium dynamics and impact synaptic plasticity in normal functions and pathologies. This model also constitutes a first step in the development of a nonlinear input-output calcium dynamics model for multi-scale, large scale neuronal simulations.

A bi-directional communication paradigm between parallel NEURON and an external non-neuron process

In order to accurately model the pattern of activation due to electrical stimulation of the hippocampus, a multi-scale computational approach is necessary. At the system level, the Admittance Method (ADM) is used to calculate the extracellular voltages created by a stimulating electrode. At the network and cellular levels, a large-scale multi-compartmental neuron network is used to calculate cellular activation. This paper presents a bi-directional communication paradigm between the NEURON model and an external surrogate for the ADM solver, where at each time step, neurons share their membrane currents with the external process, and the external process shares calculated extracellular voltages with the neuronal network. This work constitutes an important first step towards a full multi-scale NEURON-ADM model with bi-directional communication.

Decoding memory features from hippocampal spiking activities using sparse classification models

To understand how memory information is encoded in the hippocampus, we build classification models to decode memory features from hippocampal CA3 and CA1 spatio-temporal patterns of spikes recorded from epilepsy patients performing a memory-dependent delayed match-to-sample task. The classification model consists of a set of B-spline basis functions for extracting memory features from the spike patterns, and a sparse logistic regression classifier for generating binary categorical output of memory features. Results show that classification models can extract significant amount of memory information with respects to types of memory tasks and categories of sample images used in the task, despite the high level of variability in prediction accuracy due to the small sample size. These results support the hypothesis that memories are encoded in the hippocampal activities and have important implication to the development of hippocampal memory prostheses.

Monte Carlo validation of spike-timing-dependent plasticity identification from spiking activity

As increasingly complex non-stationary models become possible to be identified from neural data, rigorous validation approaches must be developed to rule out overfitting and the potential to identify features by chance. Specifically, identification of spike-timing-dependent plasticity (STDP) from recorded spontaneous in vivo spike timing is a potentially powerful tool to quantify activity-dependent plasticity. In previous work, we presented a methodology to perform this STDP identification from spike timing alone and successfully identified a generative model. Validation was straightforward with the generative model because the underlying model was known, but becomes challenging when applied to experimental data. Here, we introduce a set of null hypothesis tests that can be performed with Monte Carlo (MC) simulations of null models to rule out cases of overfitting with experimental data. We demonstrate the identification of these null models and null hypothesis testing on a generative model in two test cases, one with and one without overfitting. Importantly, we show that it is possible to distinguish an identified STDP rule from a null case where there are similar weight fluctuations which are activity-independent. With the development of the null hypothesis tests described here, STDP identification can be effectively applied to experimental data recordings.

Sources of errors and uncertainties in the assessment of forest soil carbon stocks at different scales-review and recommendations

Spatially explicit knowledge of recent and past soil organic carbon (SOC) stocks in forests will improve our understanding of the effect of human- and non-human-induced changes on forest C fluxes. For SOC accounting, a minimum detectable difference must be defined in order to adequately determine temporal changes and spatial differences in SOC. This requires sufficiently detailed data to predict SOC stocks at appropriate scales within the required accuracy so that only significant changes are accounted for. When designing sampling campaigns, taking into account factors influencing SOC spatial and temporal distribution (such as soil type, topography, climate and vegetation) are needed to optimise sampling depths and numbers of samples, thereby ensuring that samples accurately reflect the distribution of SOC at a site. Furthermore, the appropriate scales related to the research question need to be defined: profile, plot, forests, catchment, national or wider. Scaling up SOC stocks from point sample to landscape unit is challenging, and thus requires reliable baseline data. Knowledge of the associated uncertainties related to SOC measures at each particular scale and how to reduce them is crucial for assessing SOC stocks with the highest possible accuracy at each scale. This review identifies where potential sources of errors and uncertainties related to forest SOC stock estimation occur at five different scales-sample, profile, plot, landscape/regional and European. Recommendations are also provided on how to reduce forest SOC uncertainties and increase efficiency of SOC assessment at each scale.

A million-plus neuron model of the hippocampal dentate gyrus: Dependency of spatio-temporal network dynamics on topography

This paper describes a million-plus granule cell compartmental model of the rat hippocampal dentate gyrus, including excitatory, perforant path input from the entorhinal cortex, and feedforward and feedback inhibitory input from dentate interneurons. The model includes experimentally determined morphological and biophysical properties of granule cells, together with glutamatergic AMPA-like EPSP and GABAergic GABAA-like IPSP synaptic excitatory and inhibitory inputs, respectively. Each granule cell was composed of approximately 200 compartments having passive and active conductances distributed throughout the somatic and dendritic regions. Modeling excitatory input from the entorhinal cortex was guided by axonal transport studies documenting the topographical organization of projections from subregions of the medial and lateral entorhinal cortex, plus other important details of the distribution of glutamatergic inputs to the dentate gyrus. Results showed that when medial and lateral entorhinal cortical neurons maintained Poisson random firing, dentate granule cells expressed, throughout the million-cell network, a robust, non-random pattern of spiking best described as spatiotemporal "clustering". To identify the network property or properties responsible for generating such firing "clusters", we progressively eliminated from the model key mechanisms such as feedforward and feedback inhibition, intrinsic membrane properties underlying rhythmic burst firing, and/or topographical organization of entorhinal afferents. Findings conclusively identified topographical organization of inputs as the key element responsible for generating a spatio-temporal distribution of clustered firing. These results uncover a functional organization of perforant path afferents to the dentate gyrus not previously recognized: topography-dependent clusters of granule cell activity as "functional units" that organize the processing of entorhinal signals.

Design of a flexible parylene-based multi-electrode array for multi-region recording from the rat hippocampus

The hippocampus is a critical deep brain structure in several aspects. It is directly related to the formation of new long-term declarative memory. The malfunction of the hippocampus closely relates to various disease and pathological conditions. It is also a model structure for the study of cortical function and synaptic plasticity in general because of its special neuro-anatomical structure and intrinsic connections within the hippocampus formation. Both the understanding of roles that the hippocampus plays in recognition memory and the study of neural plasticity require simultaneously recording of neural activities from multiple sub-regions of the hippocampus from behavioral animals. However the distribution of cells in the hippocampus make the recording from multiple sub-regions a big challenge with the traditional uni-length micro-wire arrays. Well-designed electrode arrays are required to reach multiple regions simultaneously because of the distinctive double C shape of the hippocampus cell body layers. In this work, we designed a multi-shanks electrode which uses Parylene C, a highly biocompatible and flexible polymer, as a base and has multiple recording sites specially positioned along the longitudinal axis to fit the curvy shape of the rat hippocampus.

Sparse Large-Scale Nonlinear Dynamical Modeling of Human Hippocampus for Memory Prostheses

In order to build hippocampal prostheses for restoring memory functions, we build sparse multi-input, multi-output (MIMO) nonlinear dynamical models of the human hippocampus. Spike trains are recorded from hippocampal CA3 and CA1 regions of epileptic patients performing a variety of memory-dependent delayed match-to-sample (DMS) tasks. Using CA3 and CA1 spike trains as inputs and outputs respectively, sparse generalized Laguerre-Volterra models are estimated with group lasso and local coordinate descent methods to capture the nonlinear dynamics underlying the CA3-CA1 spike train transformations. These models can accurately predict the CA1 spike trains based on the ongoing CA3 spike trains during multiple memory events, e.g., sample presentation, sample response, match presentation and match response, of the DMS task, and thus will serve as the computational basis of human hippocampal memory prostheses.

Identification of Stable Spike-Timing-Dependent Plasticity from Spiking Activity with Generalized Multilinear Modeling

Characterization of long-term activity-dependent plasticity from behaviorally driven spiking activity is important for understanding the underlying mechanisms of learning and memory. In this letter, we present a computational framework for quantifying spike-timing-dependent plasticity (STDP) during behavior by identifying a functional plasticity rule solely from spiking activity. First, we formulate a flexible point-process spiking neuron model structure with STDP, which includes functions that characterize the stationary and plastic properties of the neuron. The STDP model includes a novel function for prolonged plasticity induction, as well as a more typical function for synaptic weight change based on the relative timing of input-output spike pairs. Consideration for system stability is incorporated with weight-dependent synaptic modification. Next, we formalize an estimation technique using a generalized multilinear model (GMLM) structure with basis function expansion. The weight-dependent synaptic modification adds a nonlinearity to the model, which is addressed with an iterative unconstrained optimization approach. Finally, we demonstrate successful model estimation on simulated spiking data and show that all model functions can be estimated accurately with this method across a variety of simulation parameters, such as number of inputs, output firing rate, input firing type, and simulation time. Since this approach requires only naturally generated spikes, it can be readily applied to behaving animal studies to characterize the underlying mechanisms of learning and memory.

A large-scale detailed neuronal model of electrical stimulation of the dentate gyrus and perforant path as a platform for electrode design and optimization

Owing to the dramatic rise in treatment of neurological disorders with electrical micro-stimulation it has become apparent that the major technological limitation in deploying effective devices lies in the process of designing efficient, safe, and outcome specific electrode arrays. The time-consuming and low-fidelity nature of gathering test data using experimental means and the immense control and flexibility of computational models, has prompted us and others to build models of electrical stimulation of neural networks that can be simulated in a computer. Because prior work has been focused on single cells, very small networks, or non-biological models of neural tissue, it was expedient that we take advantage of our, 4,040 processor, computing cluster to construct a large-scale 3-dimensional emulation of hippocampal tissue using detailed neuronal models with explicit and unique morphologies. This model, when paired with an equivalent circuit method of estimating voltage signal attenuation throughout anisotropic resistive tissue, can be used to predict tissue response to an exhaustive set of stimulation and tissue conditions: electrode geometry, array geometry, static dielectric properties of tissue, stimulation pulse features, etc. Preliminary experiments demonstrate that this system is capable of yielding neuronal responses with striking similarities to experimental results. This work provides an avenue to qualitative evaluation of electrode arrays, and more meaningful modeling of local field potentials in terms of their contributing sources and sinks.

A flexible parylene probe for in vivo recordings from multiple subregions of the rat hippocampus

The hippocampus is crucial to the formation of long-term memory and declarative memory. It is divided into three sub-fields the CA1, the CA3 and the DG. To understand the neuronal circuitry within the hippocampus and to study the role of the hippocampus in memory function requires the collection of neural activities from multiple subregions of the hippocampus simultaneously. Micro-wire electrode arrays are commonly used as an interface with neural systems. However, recording from multiple deep brain regions with curved anatomical structures such as the thin cell body layers of the hippocampus requires the micro-wires to be arranged into a highly accurate, complex layout that is difficult to fabricated manually. In this work, we designed and developed a flexible parylene-C based neural probe which can be easily micro-machined to the desired dimensions. Sixty-four electrical recording sites are micromachined on to 8 parylene shanks and spaced according to the distribution of hippocampal principal neurons in different hippocampus subregions. Together with our collaborators, we developed and optimized the implantation procedure of the flexible parylene probe and tested the insertion method both in brain tissue phantom and in vivo with a sham device. Immunohistochemistry (IHC) staining post-implantation of the sham probe was used to verify the location of the probe and to evaluate immune responses to the probe. Fully functional devices were fabricated and, in future studies, functional probes will be chronically implanted into the rat hippocampus, and neural activities will be recorded and compared with signals obtained with micro-wire arrays.

A sparse Laguerre-Volterra autoregressive model for seizure prediction in temporal lobe epilepsy

A sparse Laguerre-Volterra autoregressive model has been developed as feature extraction from subdural human EEG data for seizure prediction in temporal lobe epilepsy. The use of Laguerre-Volterra kernel can compactly yield an autoregressive model of longer system memory without increasing the number of the coefficients. In 6 sets of seizure, we used a sparse Laguerre-Volterra autoregressive model with 6 coefficients and the decay parameter of 0.2 and obtained the 10-fold cross-validation prediction results of high Matthews correlation coefficients (0.7-1) and low prediction errors (<;15%). These results demonstrate that the sparse Laguerre-Volterra autoregressive model is effective in the feature extraction for seizure prediction. Finally, this sparse Laguerre-Volterra method can be easily adapted to a potentially more powerful nonlinear autoregressive model as the feature extraction rather than linear autoregressive model that we are currently using.

Estrogen and Neural Plasticity

Converging clinical evidence suggests that postmenopausal estrogen therapy in women is associated with improved cognition and a reduced incidence of Alzheimer's disease. In experimental work, investigators have found estrogen to promote changes in synaptic plasticity within the nervous system. In this article, we review both the clinical and the experimental literature, and consider mechanisms of action of estrogen on neurons and synaptic plasticity, and how they might protect against the cognitive impairments of old age.

Implementing spiking neuron model and spike-timing-dependent plasticity with generalized Laguerre-Volterra models

To perform large-scale simulations of the brain or build biologically-inspired cognitive architectures, it is essential to have a succinct and flexible model of spiking neurons. The model should be able to capture the nonlinear dynamical properties of various types of neurons and the nonstationary properties such as the spike-timing-dependent plasticity (STDP). In this paper, we propose a generalized Laguerre-Volterra modeling approach for such a task. Due to its built-in nonlinear dynamical terms, the generalized Laguerre-Volterra model (GLVM) can capture various biological processes/mechanisms. Using Laguerre expansion of Volterra kernel technique, the model is fully represented with a small set of coefficients. The calculation of the model variables can be expressed recursively based on only the current and the one-step-before values and thus can be performed efficiently. In addition, we show that, using the same methodology, STDP can be implemented as a specific form of second-order Volterra kernel describing the causal relationship between pairs of input-output spikes and the changes of the feedforward kernels in the GLVMs.

A cognitive prosthesis for memory facilitation by closed-loop functional ensemble stimulation of hippocampal neurons in primate brain

Very productive collaborative investigations characterized how multineuron hippocampal ensembles recorded in nonhuman primates (NHPs) encode short-term memory necessary for successful performance in a delayed match to sample (DMS) task and utilized that information to devise a unique nonlinear multi-input multi-output (MIMO) memory prosthesis device to enhance short-term memory in real-time during task performance. Investigations have characterized how the hippocampus in primate brain encodes information in a multi-item, rule-controlled, delayed match to sample (DMS) task. The MIMO model was applied via closed loop feedback micro-current stimulation during the task via conformal electrode arrays and enhanced performance of the complex memory requirements. These findings clearly indicate detection of a means by which the hippocampus encodes information and transmits this information to other brain regions involved in memory processing. By employing the nonlinear dynamic multi-input/multi-output (MIMO) model, developed and adapted to hippocampal neural ensemble firing patterns derived from simultaneous recorded multi-neuron CA1 and CA3 activity, it was possible to extract information encoded in the Sample phase of DMS trials that was necessary for successful performance in the subsequent Match phase of the task. The extension of this MIMO model to online delivery of electrical stimulation patterns to the same recording loci that exhibited successful CA1 firing in the DMS Sample Phase provided the means to increase task performance on a trial-by-trial basis. Increased utility of the MIMO model as a memory prosthesis was exhibited by the demonstration of cumulative increases in DMS task performance with repeated MIMO stimulation over many sessions. These results, reported below in this article, provide the necessary demonstrations to further the feasibility of the MIMO model as a memory prosthesis to recover and/or enhance encoding of cognitive information in humans with memory disruptions resulting from brain injury, disease or aging.

The volterra functional series is a viable alternative to kinetic models for synaptic modeling--calibration and benchmarking

Synaptic transmission is governed by a series of complex and highly nonlinear mechanisms and pathways in which the dynamics have a profound influence on the overall signal sent to the postsynaptic cell. In simulation, these mechanisms are often represented through kinetic models governed by state variables and rate law equations. Calculations of such ordinary differential equations (ODEs) in kinetic models can be computationally intensive, and although algorithms have been optimally developed to handle ODEs efficiently, simulation of numerous, large and complex kinetic models requires a prohibitively large amount of computational power. Here we present an alternative representation of ionotropic glutamatergic receptors AMPAr and NMDAr kinetic models consisting of input-output surrogates of the receptor models which can capture the nonlinear dynamics seen in the kinetic models. We benchmark this Input-Output (IO) synapse model and compare it with kinetic receptor models to evaluate the simulation time required when using either synapse model, as well as the number of time steps each model needs for simulation. While remaining faithful to the original dynamics of the model, our results indicate that the IO synapse model requires less simulation time than the kinetic models under conditions which elicit normal physiological responses, thereby improving computational efficiency while preserving the complex non-linear dynamics of the receptors. These IO surrogates therefore constitute an appealing alternative to kinetic models in large scale networks simulations.