Estimating the biophysical properties of neurons with intracellular calcium dynamics

Reduced-Dimension, Biophysical Neuron Models Constructed From Observed Data

Using methods from nonlinear dynamics and interpolation techniques from applied mathematics, we show how to use data alone to construct discrete time dynamical rules that forecast observed neuron properties. These data may come from simulations of a Hodgkin-Huxley (HH) neuron model or from laboratory current clamp experiments. In each case, the reduced-dimension, data-driven forecasting (DDF) models are shown to predict accurately for times after the training period. When the available observations for neuron preparations are, for example, membrane voltage V(t) only, we use the technique of time delay embedding from nonlinear dynamics to generate an appropriate space in which the full dynamics can be realized. The DDF constructions are reduced-dimension models relative to HH models as they are built on and forecast only observables such as V(t). They do not require detailed specification of ion channels, their gating variables, and the many parameters that accompany an HH model for laboratory measurements, yet all of this important information is encoded in the DDF model. As the DDF models use and forecast only voltage data, they can be used in building networks with biophysical connections. Both gap junction connections and ligand gated synaptic connections among neurons involve presynaptic voltages and induce postsynaptic voltage response. Biophysically based DDF neuron models can replace other reduced-dimension neuron models, say, of the integrate-and-fire type, in developing and analyzing large networks of neurons. When one does have detailed HH model neurons for network components, a reduced-dimension DDF realization of the HH voltage dynamics may be used in network computations to achieve computational efficiency and the exploration of larger biological networks.

Identifying the measurements required to estimate rates of COVID-19 transmission, infection, and detection, using variational data assimilation

We demonstrate the ability of statistical data assimilation (SDA) to identify the measurements required for accurate state and parameter estimation in an epidemiological model for the novel coronavirus disease COVID-19. Our context is an effort to inform policy regarding social behavior, to mitigate strain on hospital capacity. The model unknowns are taken to be: the time-varying transmission rate, the fraction of exposed cases that require hospitalization, and the time-varying detection probabilities of new asymptomatic and symptomatic cases. In simulations, we obtain estimates of undetected (that is, unmeasured) infectious populations, by measuring the detected cases together with the recovered and dead - and without assumed knowledge of the detection rates. Given a noiseless measurement of the recovered population, excellent estimates of all quantities are obtained using a temporal baseline of 101 days, with the exception of the time-varying transmission rate at times prior to the implementation of social distancing. With low noise added to the recovered population, accurate state estimates require a lengthening of the temporal baseline of measurements. Estimates of all parameters are sensitive to the contamination, highlighting the need for accurate and uniform methods of reporting. The aim of this paper is to exemplify the power of SDA to determine what properties of measurements will yield estimates of unknown parameters to a desired precision, in a model with the complexity required to capture important features of the COVID-19 pandemic.

Statistical data assimilation for estimating electrophysiology simultaneously with connectivity within a biological neuronal network

A method of data assimilation (DA) is employed to estimate electrophysiological parameters of neurons simultaneously with their synaptic connectivity in a small model biological network. The DA procedure is cast as an optimization, with a cost function consisting of both a measurement error and a model error term. An iterative reweighting of these terms permits a systematic method to identify the lowest minimum, within a local region of state space, on the surface of a nonconvex cost function. In the model, two sets of parameter values are associated with two particular functional modes of network activity: simultaneous firing of all neurons and a pattern-generating mode wherein the neurons burst in sequence. The DA procedure is able to recover these modes if: (i) the stimulating electrical currents have chaotic waveforms and (ii) the measurements consist of the membrane voltages of all neurons in the circuit. Further, this method is able to prune a model of unnecessarily high dimensionality to a representation that contains the maximum dimensionality required to reproduce the provided measurements. This paper offers a proof-of-concept that DA has the potential to inform laboratory designs for estimating properties in small and isolatable functional circuits.

A mathematical and computational model of the calcium dynamics in Caenorhabditis elegans ASH sensory neuron

We propose a mathematical and computational model that captures the stimulus-generated Ca²⁺ transients in the C. elegans ASH sensory neuron. The rationale is to develop a tool that will enable a cross-talk between modeling and experiments, using modeling results to guide targeted experimental efforts. The model is built based on biophysical events and molecular cascades known to unfold as part of neurons' Ca²⁺ homeostasis mechanism, as well as on Ca²⁺ signaling events. The state of ion channels is described by their probability of being activated or inactivated, and the remaining molecular states are based on biochemically defined kinetic equations or known biochemical motifs. We estimate the parameters of the model using experimental data of hyperosmotic stimulus-evoked Ca²⁺ transients detected with a FRET sensor in young and aged worms, unstressed and exposed to oxidative stress. We use a hybrid optimization method composed of a multi-objective genetic algorithm and nonlinear least-squares to estimate the model parameters. We first obtain the model parameters for young unstressed worms. Next, we use these values of the parameters as a starting point to identify the model parameters for stressed and aged worms. We show that the model, in combination with experimental data, corroborates literature results. In addition, we demonstrate that our model can be used to predict ASH response to complex combinations of stimulation pulses. The proposed model includes for the first time the ASH Ca²⁺ dynamics observed during both "on" and "off" responses. This mathematical and computational effort is the first to propose a dynamic model of the Ca²⁺ transients' mechanism in C. elegans neurons, based on biochemical pathways of the cell's Ca²⁺ homeostasis machinery. We believe that the proposed model can be used to further elucidate the Ca²⁺ dynamics of a key C. elegans neuron, to guide future experiments on C. elegans neurobiology, and to pave the way for the development of more mathematical models for neuronal Ca²⁺ dynamics.

Nonlinear statistical data assimilation for HVC[Formula: see text] neurons in the avian song system

With the goal of building a model of the HVC nucleus in the avian song system, we discuss in detail a model of HVC[Formula: see text] projection neurons comprised of a somatic compartment with fast Na[Formula: see text] and K[Formula: see text] currents and a dendritic compartment with slower Ca[Formula: see text] dynamics. We show this model qualitatively exhibits many observed electrophysiological behaviors. We then show in numerical procedures how one can design and analyze feasible laboratory experiments that allow the estimation of all of the many parameters and unmeasured dynamical variables, given observations of the somatic voltage [Formula: see text] alone. A key to this procedure is to initially estimate the slow dynamics associated with Ca, blocking the fast Na and K variations, and then with the Ca parameters fixed estimate the fast Na and K dynamics. This separation of time scales provides a numerically robust method for completing the full neuron model, and the efficacy of the method is tested by prediction when observations are complete. The simulation provides a framework for the slice preparation experiments and illustrates the use of data assimilation methods for the design of those experiments.