Search for: All records

In neuroscience, delayed synaptic activity plays a pivotal and pervasive role in influencing synchronization, oscillation, and information-processing properties of neural networks. In small rhythm-generating networks, such as central pattern generators (CPGs), time-delays may regulate and determine the stability and variability of rhythmic activity, enabling organisms to adapt to environmental changes, and coordinate diverse locomotion patterns in both function and dysfunction. Here, we examine the dynamics of a three-cell CPG model in which time-delays are introduced into reciprocally inhibitory synapses between constituent neurons. We employ computational analysis to investigate the multiplicity and robustness of various rhythms observed in such multi-modal neural networks. Our approach involves deriving exhaustive two-dimensional Poincaré return maps for phase-lags between constituent neurons, where stable fixed points and invariant curves correspond to various phase-locked and phase-slipping/jitter rhythms. These rhythms emerge and disappear through various local (saddle-node, torus) and non-local (homoclinic) bifurcations, highlighting the multi-functionality (modality) observed in such small neural networks with fast inhibitory synapses. 

This paper investigates the origin and onset of chaos in a mathematical model of an individual neuron, arising from the intricate interaction between 3D fast and 2D slow dynamics governing its intrinsic currents. Central to the chaotic dynamics are multiple homoclinic connections and bifurcations of saddle equilibria and periodic orbits. This neural model reveals a rich array of codimension-2 bifurcations, including Shilnikov–Hopf, Belyakov, Bautin, and Bogdanov–Takens points, which play a pivotal role in organizing the complex bifurcation structure of the parameter space. We explore various routes to chaos occurring at the intersections of quiescent, tonic spiking, and bursting activity regimes within this space and provide a thorough bifurcation analysis. Despite the high dimensionality of the model, its fast–slow dynamics allow a reduction to a one-dimensional return map, accurately capturing and explaining the complex dynamics of the neural model. Our approach integrates parameter continuation analysis, newly developed symbolic techniques, and Lyapunov exponents, collectively unveiling the intricate dynamical and bifurcation structures present in the system. 

We present a novel set of quantitative measures for “likeness” (error function) designed to alleviate the time-consuming and subjective nature of manually comparing biological recordings from electrophysiological experiments with the outcomes of their mathematical models. Our innovative “blended” system approach offers an objective, high-throughput, and computationally efficient method for comparing biological and mathematical models. This approach involves using voltage recordings of biological neurons to drive and train mathematical models, facilitating the derivation of the error function for further parameter optimization. Our calibration process incorporates measurements such as action potential (AP) frequency, voltage moving average, voltage envelopes, and the probability of post-synaptic channels. To assess the effectiveness of our method, we utilized the sea slug Melibe leonina swim central pattern generator (CPG) as our model circuit and conducted electrophysiological experiments with TTX to isolate CPG interneurons. During the comparison of biological recordings and mathematically simulated neurons, we performed a grid search of inhibitory and excitatory synapse conductance. Our findings indicate that a weighted sum of simple functions is essential for comprehensively capturing a neuron’s rhythmic activity. Overall, our study suggests that our blended system approach holds promise for enabling objective and high-throughput comparisons between biological and mathematical models, offering significant potential for advancing research in neural circuitry and related fields.