Stochastic Hybrid Systems in Cellular Neuroscience

Efficient stochastic simulation of piecewise-deterministic Markov processes and its application to the Morris-Lecar model of neural dynamics

Piecewise-deterministic Markov processes combine continuous in time dynamics with jump events, the rates of which generally depend on the continuous variables and thus are not constants. This leads to a problem in a Monte-Carlo simulation of such a system, where, at each step, one must find the time instant of the next event. The latter is determined by an integral equation and usually is rather slow in numerical implementation. We suggest a reformulation of the next event problem as an ordinary differential equation where the independent variable is not the time but the cumulative rate. This reformulation is similar to the Hénon approach to efficiently constructing the Poincaré map in deterministic dynamics. The problem is then reduced to a standard numerical task of solving a system of ordinary differential equations with given initial conditions on a prescribed interval. We illustrate the method with a stochastic Morris-Lecar model of neuron spiking with stochasticity in the opening and closing of voltage-gated ion channels.

Competition of two time scales determines the performance of a voltage-gated potassium channel

The dynamics of a voltage-gated potassium channel in real environments represents a crucial bridge between its molecular structure and functions. However, it is still missing due to the mathematical difficulty that arises from the high dimensionality and nonlinear interregulation. Here we present a method for solving the stationary distribution of a hybrid process that contains two negatively interregulating kinetics: channel gating and voltage decay. The results can be summarized as follows: first, the voltage distribution is determined by the competition of their time scales; second, the fluctuation structures in parameter space illustrate that, to perform the voltage-controlling task, the channel gating is elastic while the membrane produces the stabilizing function; third, the power dissipated by the capacitive currents and the internal battery current are calculated and explained. Based on these findings, we examine the manner in which macroscopic functions of potassium channels are manifested. Our methodology provides an accurate characterization of hybrid processes that are pervasive in the life sciences.

Modular architecture facilitates noise-driven control of synchrony in neuronal networks

High-level information processing in the mammalian cortex requires both segregated processing in specialized circuits and integration across multiple circuits. One possible way to implement these seemingly opposing demands is by flexibly switching between states with different levels of synchrony. However, the mechanisms behind the control of complex synchronization patterns in neuronal networks remain elusive. Here, we use precision neuroengineering to manipulate and stimulate networks of cortical neurons in vitro, in combination with an in silico model of spiking neurons and a mesoscopic model of stochastically coupled modules to show that (i) a modular architecture enhances the sensitivity of the network to noise delivered as external asynchronous stimulation and that (ii) the persistent depletion of synaptic resources in stimulated neurons is the underlying mechanism for this effect. Together, our results demonstrate that the inherent dynamical state in structured networks of excitable units is determined by both its modular architecture and the properties of the external inputs.

Multiscale molecular kinetics by coupling Markov state models and reaction-diffusion dynamics

A novel approach to simulate simple protein-ligand systems at large time and length scales is to couple Markov state models (MSMs) of molecular kinetics with particle-based reaction-diffusion (RD) simulations, MSM/RD. Currently, MSM/RD lacks a mathematical framework to derive coupling schemes, is limited to isotropic ligands in a single conformational state, and lacks multiparticle extensions. In this work, we address these needs by developing a general MSM/RD framework by coarse-graining molecular dynamics into hybrid switching diffusion processes. Given enough data to parameterize the model, it is capable of modeling protein-protein interactions over large time and length scales, and it can be extended to handle multiple molecules. We derive the MSM/RD framework, and we implement and verify it for two protein-protein benchmark systems and one multiparticle implementation to model the formation of pentameric ring molecules. To enable reproducibility, we have published our code in the MSM/RD software package.

Modelling the motion of organelles in an elongated cell via the coordination of heterogeneous drift-diffusion and long-range transport

 Cellular distribution of organelles in living cells is achieved via a variety of transport mechanisms, including directed motion, mediated by molecular motors along microtubules (MTs), and diffusion which is predominantly heterogeneous in space. In this paper, we introduce a model for particle transport in elongated cells that couples poleward drift, long-range bidirectional transport and diffusion with spatial heterogeneity in a three-dimensional space. Using stochastic simulations and analysis of a related population model, we find parameter regions where the three-dimensional model can be reduced to a coupled one-dimensional model or even a one-dimensional scalar model. We explore the efficiency with which individual model components can overcome drift towards one of the cell poles to reach an approximately even distribution. In particular, we find that if lateral movement is well mixed, then increasing the binding ability of particles to MTs is an efficient way to overcome a poleward drift, whereas if lateral motion is not well mixed, then increasing the axial diffusivity away from MTs becomes an efficient way to overcome the poleward drift. Our three-dimensional model provides a new tool that will help to understand the mechanisms by which eukaryotic cells organize their organelles in an elongated cell, and in particular when the one-dimensional models are applicable.

Special Issue from the 2017 International Conference on Mathematical Neuroscience

The ongoing acquisition of large and multifaceted data sets in neuroscience requires new mathematical tools for quantitatively grounding these experimental findings. Since 2015, the International Conference on Mathematical Neuroscience (ICMNS) has provided a forum for researchers to discuss current mathematical innovations emerging in neuroscience. This special issue assembles current research and tutorials that were presented at the 2017 ICMNS held in Boulder, Colorado from May 30 to June 2. Topics discussed at the meeting include correlation analysis of network activity, information theory for plastic synapses, combinatorics for attractor neural networks, and novel data assimilation methods for neuroscience-all of which are represented in this special issue.