We propose a mathematical and computational model that captures the stimulusgenerated Ca 2+ transients in the C. elegans ASH sensory neuron. The model is built based on biophysical events and molecular cascades known to unfold as part of neurons' Ca 2+ homeostasis mechanism, as well as on Ca 2+ 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 with phenomenological adjustments. We estimate the parameters of the model using experimental data of hyperosmotic stimulus-evoked Ca 2+ 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 2+ dynamics observed during both "on" and "off"responses. This mathematical and computational effort is the first to propose a dynamic model of the Ca 2+ transients' mechanism in C. elegans neurons, based on biochemical pathways of the cell's Ca 2+ homeostasis machinery.Key words: C. elegans, calcium dynamics, ASH neuron, aging, oxidative stress
Significance StatementC. elegans is widely used as a model system for monitoring neuronal Ca 2+ transients.The ASH neuron is the subject of several such studies, primarily due to its key importance as a polymodal nociceptor. However, despite its pivotal role in C. elegans biology, and the special characteristics of its stimulus-evoked Ca2 + transients (e.g., the peer-reviewed) is the author/funder. All rights reserved. No reuse allowed without permission.The copyright holder for this preprint (which was not . http://dx.doi.org/10.1101/201962 doi: bioRxiv preprint first posted online Oct. 12, 2017; "off" response), no mathematical or computational model has been developed to include special features of ASH Ca 2+ dynamics, i.e. the "off" response. The model includes for the first time the ASH Ca 2+ dynamics observed during both "on" and "off" responses, and is the first to propose a dynamical model of the C. elegans Ca 2+ transients' mechanism based on biochemical pathways of the cell's Ca 2+ homeostasis machinery.