SNNAP: A Tool for Teaching Neuroscience

Simulation tools aid in learning neuroscience by providing the student with an interactive environment to carry out simulated experiments and test hypotheses. The field of neuroscience is well suited for the use of simulation tools since nerve cell signaling can be described by mathematical equations and solved by computer. Neural signaling entails the propagation of electrical current along nerve membrane and transmission to neighboring neurons through synaptic connections. Action potentials and synaptic transmission can be simulated and results displayed for visualization and analysis. The neurosimulator SNNAP (Simulator for Neural Networks and Action Potentials) is a simulation environment that provides users with editors for model building, simulator engine and visual display editor. This paper presents several modeling examples that illustrate some of the capabilities and features of SNNAP. First, the Hodgkin-Huxley (HH) model is presented and the threshold phenomenon is illustrated. Second, small neural networks are described with HH models using various synaptic connections available with SNNAP. Synaptic connections may be modulated through facilitation or depression with SNNAP. A study of vesicle pool dynamics is presented using an AMPA receptor model. Finally, a central pattern generator model of the Aplysia feeding circuit is illustrated as an example of a complex network that may be studied with SNNAP. Simulation code is provided for each case study described and tasks are suggested for further investigation

GEOMETRY, ACTIVITY-DEPENDENT MECHANISMS, MEMBRANE KINETICS AND CHANNEL DENSITY DISTRIBUTION INTERPLAY IN SINGLE NEURON PLASTICITY: A COMPUTATIONAL STUDY

Conduction of action potentials throughout the complex morphology of neurons may be modulated in an activity-dependent manner. Among modulatory mechanisms, afterhyperpolarization (AHP) plays an important role. To investigate how the AHP modulatory capabilities on transmission were dependent on the axonal geometry as well as on membrane properties such as channel kinetics, channel density distribution and membrane noise, multi-compartment computational neural models were built, using the neurosimulator SNNAP. Two kinetic schema for the sodium and potassium channels were compared. The simulations suggest that channel kinetics profoundly influence the AHP-dependent modulation of action potential conduction through points of impedance mismatch in the highly branched neurites of neurons

Bifurcation and Singularity Analysis of a Molecular Network for the Induction of Long-Term Memory

Withdrawal reflexes of the mollusk Aplysia exhibit sensitization, a simple form of long-term memory (LTM). Sensitization is due, in part, to long-term facilitation (LTF) of sensorimotor neuron synapses. LTF is induced by the modulatory actions of serotonin (5-HT). Pettigrew et al. developed a computational model of the nonlinear intracellular signaling and gene network that underlies the induction of 5-HT-induced LTF. The model simulated empirical observations that repeated applications of 5-HT induce persistent activation of protein kinase A (PKA) and that this persistent activation requires a suprathreshold exposure of 5-HT. This study extends the analysis of the Pettigrew model by applying bifurcation analysis, singularity theory, and numerical simulation. Using singularity theory, classification diagrams of parameter space were constructed, identifying regions with qualitatively different steady-state behaviors. The graphical representation of these regions illustrates the robustness of these regions to changes in model parameters. Because persistent protein kinase A (PKA) activity correlates with Aplysia LTM, the analysis focuses on a positive feedback loop in the model that tends to maintain PKA activity. In this loop, PKA phosphorylates a transcription factor (TF-1), thereby increasing the expression of an ubiquitin hydrolase (Ap-Uch). Ap-Uch then acts to increase PKA activity, closing the loop. This positive feedback loop manifests multiple, coexisting steady states, or multiplicity, which provides a mechanism for a bistable switch in PKA activity. After the removal of 5-HT, the PKA activity either returns to its basal level (reversible switch) or remains at a high level (irreversible switch). Such an irreversible switch might be a mechanism that contributes to the persistence of LTM. The classification diagrams also identify parameters and processes that might be manipulated, perhaps pharmacologically, to enhance the induction of memory. Rational drug design, to affect complex processes such as memory formation, can benefit from this type of analysis

Computational Model of Touch Sensory Cells (T Cells) of the Leech:Role of the Afterhyperpolarization (AHP) in Activity-Dependent Conduction Failure

Abstract. Bursts of spikes in T cells produce an AHP, which results from activation of a Na+/K+ pump and a Ca2+-dependent K+ current. Activity-dependent increases in the AHP are believed to induce conduction block of spikes in several regions of the neuron, which in turn, may decrease presynaptic invasion of spikes and thereby decrease transmitter release. To explore this possibility, we used the neurosimulator SNNAP to develop a multi-compartmental model of the T cell. The model incorporated empirical data that describe the geometry of the cell and activity-dependent changes of the AHP. Simulations indicated that at some branching points, activity-dependent increases of the AHP reduced the number of spikes transmitted from the minor receptive fields to the soma and beyond. More importantly, simulations also suggest that the AHP could modulate, under some circumstances, transmission from the soma to the synaptic terminals, suggesting that the AHP can regulate spike conduction within the presynaptic arborizations of the cell and could in principle contribute to the synaptic depression that is correlated with increases in the AHP

Dynamics of a Minimal Model of Interlocked Positive and Negative Feedback Loops of Transcriptional Regulation by cAMP-Response Element Binding Proteins

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture