The problem of neuron modeling
Neuron modeling - simulation of neuronal behavior, model based on Hodgkin and Huxley equations, and simplified model using properties of pulse frequency modulated systems
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Neuron modeling - simulation of neuronal behavior, model based on Hodgkin and Huxley equations, and simplified model using properties of pulse frequency modulated systems
Single neuron models are fundamental for computational modeling of the brain's neuronal networks, and understanding how ion channel dynamics mediate neural function. A challenge in defining such models is determining biophysically realistic channel distributions. Here, we present an efficient, highly parallel evolutionary algorithm for developing such models, named NeuroGPU-EA. NeuroGPU-EA uses CPUs and GPUs concurrently to simulate and evaluate neuron membrane potentials with respect to multiple stimuli. We demonstrate a logarithmic cost for scaling the stimuli used in the fitting procedure. NeuroGPU-EA outperforms the typically used CPU based evolutionary algorithm by a factor of 10 on a series of scaling benchmarks. We report observed performance bottlenecks and propose mitigation strategies. Finally, we also discuss the potential of this method for efficient simulation and evaluation of electrophysiological waveforms.
An initial investigation into neuron models aimed at efficient stochastic simulation of turbulent flows.
A method for the rapid learning of nonlinear mappings and topological transformations using a dynamically reconfigurable artificial neural network is presented. This fully-recurrent Adaptive Neuron Model (ANM) network was applied to the highly degenerate inverse kinematics problem in robotics, and its performance evaluation is bench-marked. Once trained, the resulting neuromorphic architecture was implemented in custom analog neural network hardware and the parameters capturing the functional transformation downloaded onto the system. This neuroprocessor, capable of 10(exp 9) ops/sec, was interfaced directly to a three degree of freedom Heathkit robotic manipulator. Calculation of the hardware feed-forward pass for this mapping was benchmarked at approximately 10 microsec.
Melatonin, a hormone primarily produced by the brain’s pineal gland, not only regulates circadian rhythms, but also influences the structural and biophysical properties of neuronal membranes. Its amphiphilic nature enables direct incorporation into lipid bilayers and preferential interactions with cholesterol-rich lipid rafts, critical hubs for cellular signaling and membrane organization. Despite increasing recognition of its membrane activity, the molecular basis of melatonin’s interactions with coexisting liquid-ordered (L o ) and liquid-disordered (L d ) phases remains unclear. Here, in this study, we combine small-angle neutron scattering (SANS) and all-atom molecular dynamics simulations to examine model neuronal membranes composed of DSPC, DOPC, POPC, and cholesterol. Our results show that melatonin preserves domain morphology while adopting distinct orientations within the bilayer and at the membrane interface, allowing both lateral and transmembrane bridging across lipid phases. These findings establish the molecular underpinnings of melatonin’s modulation of membrane heterogeneity and provide strong support for its receptor-independent actions.
One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.
We propose a modification to the Hodgkin-Huxley model to reduce the numerical stiffness of the equations by introducing an explicit voltage threshold. When this threshold is crossed, the voltage and the gating variables are reset to constant values. It is found that, for all of the current stimuli considered, the proposed model accurately reproduces the behavior of the baseline Hodgkin-Huxley model while bypassing the fast dynamics of spiking. Specifically, the model accurately reproduces the spike times and, between spikes, the time courses of the membrane potential and gating variables.
Explore the source record for details and available documents.
The Neuron Ring model employs an avalanche structure with two important distinctions at the neuron level. Each neuron has two memory latches; one traps maximum neuronal activation during pattern presentation, and the other records the time of latch content change. The latches filter short term memory. In the process, they preserve length 1 snapshots of activation theory history. The model finds utility in pattern classification. Its synaptic weights are first conditioned with sample spectra. The model then receives a test or unknown signal. The objective is to identify the sample closest to the test signal. Class decision follows complete presentation of the test data. The decision maker relies exclusively on the latch contents. Presented here is an overview of the Neuron Ring at the seminar level.
Results of a previous metric analysis and an electronic simulation of acceleratory nystagmus are given. On this basis, a tentative mathematical model for describing acceleratory nystagmus is reported. The essential content of the model is Lorente de No's neuron circuit, to which the two-factor theory of excitation has been applied.
Efficient representation of images for human use requires an understanding of how the brain processes and represents visual information. Spatial imagery is represented in the brain in the receptive fields of visual neurons. Models of these neurons lead to models of image fidelity, and to digital implementations of these neural codes. This approach will be illustrated by two example codes. The advantages and difficulties of this approach will be discussed.
The paper presents an artificial neural network concept (the Synchronizable Oscillator Networks) where the instants of individual firings in the form of point processes constitute the only form of information transmitted between joining neurons. In the model, neurons fire spontaneously and regularly in the absence of perturbation. When interaction is present, the scheduled firings are advanced or delayed by the firing of neighboring neurons. Networks of such neurons become global oscillators which exhibit multiple synchronizing attractors. From arbitrary initial states, energy minimization learning procedures can make the network converge to oscillatory modes that satisfy multi-dimensional constraints. Such networks can directly represent routing and scheduling problems that consist of ordering sequences of events.
Most neuromorphic chips implement the single-compartment point neuron model where synapse circuits connect directly to a leaky integrate and fire (LIF) soma circuit. However, when using a biologically plausible soma circuit (e.g., Hodgkin-Huxley neuron model), an interface circuitry, such as a current conveyor circuit, is needed to transmit synaptic current to the soma circuit. This is especially true for ultra-low power neuron circuits, where membrane capacitance is on the order of 20 fF. This need for an interface circuit arises because the parasitic capacitance and leakage current caused by fabrication mismatch and second-order effects of the output transistors in the synapse circuits can disturb the spiking dynamics of the soma circuit if connected without an interface. Using an interface circuit to isolate the soma’s membrane capacitor from synapses resolves this issue. We propose to use a unidirectional resistor (a transconductance circuit) to connect the synapse and soma circuits instead of conventional current conveyor circuits. Using a biologically plausible spike pattern detection model, we show that the on-chip spike-timing-dependent plasticity (STDP) learning performance of the proposed unidirectional two-compartment neuron circuit is similar to a single-compartment circuit (with a current conveyor as an interface) and additionally, it is more power-efficient and biologically plausible. The chip is fabricated in the Taiwan Semiconductor Manufacturing Company (TSMC) 250 nm technology node and comprises a single neuron circuit.
Three adaptive single-neuron models based on neural analogies of behavior modification episodes are proposed, which attempt to bridge the gap between psychology and neurophysiology. The proposed models capture the predictive nature of Pavlovian conditioning, which is essential to the theory of adaptive/learning systems. The models learn to anticipate the occurrence of a conditioned response before the presence of a reinforcing stimulus when training is complete. Furthermore, each model can find the most nonredundant and earliest predictor of reinforcement. The behavior of the models accounts for several aspects of basic animal learning phenomena in Pavlovian conditioning beyond previous related models. Computer simulations show how well the models fit empirical data from various animal learning paradigms.
This work presents a novel method for generating electricity price scenarios from statistical properties of past electricity prices using a hybrid statistical and reduced-form stochastic model. Previous work in applying stochastic differential equations (SDE) to model electricity prices has focused on daily average prices. To extend stochastic price generation methods to hourly or sub-hourly pricing, we address several weaknesses in the state-of-the-art: (1) we replace the mean-reversion component of the SDE with an ARIMA process that is better able to characterize the daily and weekly trends; (2) we extend the price-spike, or jump process to account for conditional probabilities of price spikes occurring in consecutive time steps by replacing the traditional Poisson process for modeling jumps with a generalized point process model inspired by brain neuron models; and (3) we replace the traditional method of estimating spike intensity with empirical variance with a Markov process based on observed price spike intensity transitions. The method is demonstrated with electricity prices from the US ERCOT market and a use-case example is provided for bidding an energy storage unit into the day-ahead and real-time energy markets of ERCOT using stochastic optimization methods. Results show that the the synthetic price model out performs a (naive) persistence forecast model by resulting in 24% to 47% more in profits over 168 simulated days.
Accuracy-optimized convolutional neural networks (CNNs) have emerged as highly effective models at predicting neural responses in brain areas along the primate ventral stream, but it is largely unknown whether they effectively model neurons in the complementary primate dorsal stream. We explored how well CNNs model the optic flow tuning properties of neurons in dorsal area MSTd and we compared our results with the Non-Negative Matrix Factorization (NNMF) model, which successfully models many tuning properties of MSTd neurons. To better understand the role of computational properties in the NNMF model that give rise to optic flow tuning that resembles that of MSTd neurons, we created additional CNN model variants that implement key NNMF constraints – non-negative weights and sparse coding of optic flow. While the CNNs and NNMF models both accurately estimate the observer's self-motion from purely translational or rotational optic flow, NNMF and the CNNs with nonnegative weights yield substantially less accurate estimates than the other CNNs when tested on more complex optic flow that combines observer translation and rotation. Despite its poor accuracy, NNMF gives rise to tuning properties that align more closely with those observed in primate MSTd than any of the accuracy-optimized CNNs. This work offers a step toward a deeper understanding of the computational properties and constraints that describe the optic flow tuning of primate area MSTd.
Highlights: • Establishment of a patient-specific neuronal model for Mucopolysaccharidosis type IIIB using iPSCs. • Lysosomes enlarged and secondary lipid accumulated in Mucopolysaccharidosis type IIIB iPSC-derived neuronal cells. • Recombinant NAGLU, δ-tocopherol, and HPBCD reduced the disease phenotypes in MPS IIIB iPSC-derived neuronal cells. Mucopolysaccharidosis type IIIB (MPS IIIB) is a lysosomal disease caused by mutations in the NAGLU gene encoding α-N-acetylglucosaminidase (NAGLU) which degrades heparan sulfate in lysosomes. Deficiency in NAGLU results in lysosomal accumulation of glycosaminoglycans (GAGs) and neurological symptoms. Currently, there is no effective treatment or cure for this disease. In this study, induced pluripotent stem cell lines were established from two MPS IIIB patient fibroblast lines and differentiated into neural stem cells and neurons. MPS IIIB neural stem cells exhibited NAGLU deficiency accompanied with GAG accumulation, as well as lysosomal enlargement and secondary lipid accumulation. Treatments with recombinant NAGLU, δ-tocopherol, and 2-hydroxypropyl-b-cyclodextrin significantly reduced the disease phenotypes in these cells. These results indicate the MPS IIIB neural stem cells and neurons have the disease relevant phenotype and can be used as a cell-based disease model system for evaluation of drug efficacy and compound screening for drug development.
Representative, reproducible and high-throughput models of human cells and tissues are critical for a meaningful evaluation of host-pathogen interactions and are an essential component of the research developmental pipeline. The most informative infection models - animals, organ explants and human trials - are not suited for extensive evaluation of pathogenesis mechanisms and screening of candidate drugs. At the other extreme, more cost effective and accessible infection models such as conventional cell culture and static co-culture may not capture physiological and three-dimensional aspects of tissue biology that are important in assessing pathogenesis, and effectiveness and cytotoxicity of therapeutics. Our lab has used innovative bioengineering technology to establish biologically meaningful 3-D models of human tissues that recapitulate many aspects of the differentiated structure and function of the parental tissue in vivo, and we have applied these models to study infectious disease. We have established a variety of different 3-D models that are currently being used in infection studies - including small intestine, colon, lung, placenta, bladder, periodontal ligament, and neuronal models. Published work from our lab has shown that our 3-D models respond to infection with bacterial and viral pathogens in ways that reflect the infection process in vivo. By virtue of their physiological relevance, 3-D cell cultures may also hold significant potential as models to provide insight into the neuropathogenesis of HIV infection. Furthermore, the experimental flexibility, reproducibility, cost-efficiency, and high throughput platform afforded by these 3-D models may have important implications for the design and development of drugs with which to effectively treat neurological complications of HIV infection.