Skip to main content

Uncertainty Propagation in Nerve Impulses Through the Action Potential Mechanism

Abstract

We investigate the propagation of probabilistic uncertainty through the action potential mechanism in nerve cells. Using the Hodgkin–Huxley (H-H) model and Stochastic Collocation on Sparse Grids, we obtain an accurate probabilistic interpretation of the deterministic dynamics of the transmembrane potential and gating variables. Using Sobol indices, out of the 11 uncertain parameters in the H-H model, we unravel two main uncertainty sources, which account for more than 90 % of the fluctuations in neuronal responses, and have a direct biophysical interpretation. We discuss how this interesting feature of the H-H model allows one to reduce greatly the probabilistic degrees of freedom in uncertainty quantification analyses, saving CPU time in numerical simulations and opening possibilities for probabilistic generalisation of other deterministic models of great importance in physiology and mathematical neuroscience.

1 Findings

To provide sensible mechanistic models of neuronal circuits in the brain and to facilitate principled interpretations of data from experiments in vitro and in vivo, biophysical models of neurodynamic systems must incorporate variability as a defining feature [16]. Neuronal variability cannot be accounted for by dynamic models formulated as deterministic differential equations. However, such models can be extended to do so by incorporating probabilistic degrees of freedom. For this purpose, the versatile technique of Stochastic Collocation on Sparse Grids (SCSG) [79], from Uncertainty Quantification (UQ) theory [7, 1013], has been recently suggested as a computationally efficient strategy [14].

In this report, we deploy SCSG to investigate uncertainty propagation through the action potential mechanism in the H-H model [15], namely

(1)

with , where is the electrical potential across a neuron’s membrane, , are the gating variables associated with the activation and inactivation of ion currents, respectively, is the gating variable associated with the activation of ions current and where I represents current injected through the cell membrane.

The ODE system (1) exhibits a combination of excitability and highly nonlinear dynamics. We will show how these features can render the model output extremely sensitive to parameter fluctuations under realistic physiological conditions, demonstrating that Uncertainty Quantification might be indispensable to provide sensible biophysical models of nerve impulses. The H-H model [15] includes a vector of 11 parameters, which in state space form reads

(2)

consisting of four initial conditions , , , , three parameters describing the maximum conductances , , , corresponding to , and leakage ion currents respectively, three parameters describing their equilibrium Nernst potentials , , , and the membrane capacitance C. These parameters are uncertain, for they must be determined experimentally and are therefore subject to errors and fluctuations. Their nominal values measured in squid axon preparations are listed in Table 2 [15].

In the probabilistic generalisation of the H-H model, a Banach space describing such 11 probabilistic degrees of freedom is required for uncertainty analysis. Without appropriate discretisation and dimensionality reduction procedures, the computational cost of such analysis depends exponentially on the number of parameters, challenging the capacities of current hardware.

An efficient numerical strategy to deal with such a curse of dimensionality, is to find a selection of points in probability space that yields a good approximation of the model’s response surface. This can be achieved by constructing a Sparse Grid in the multi-dimensional probability space for each time point following Smolyak’s algorithm [8, 9, 1618].

Using Sparse Grid quadrature in conjunction with Sobol indices, we unravelled the parameters whose fluctuation causes most of the observed variability in neuronal responses. Such an analysis hinted at a parsimonious H-H model with two probabilistic degrees of freedom instead of 11, the use of which can achieve tremendous savings in CPU time as discussed below.

The uncertain dynamics of the membrane potential during neuronal discharge is shown in Fig. 1(a). The probability distribution at each time point was obtained by Monte Carlo sampling of a third level Sparse Grid piecewise linear interpolant [8, 19] of the response surface for the H-H model, using 1000 samples per time point from a uniform probability distribution and assuming 20 % of variability in the nominal values of the parameters listed in Table 2. The mean and deterministic solutions are also shown. Notice that they differ when the membrane depolarisation reaches its acme, which is often a signature of high nonlinearity.

Fig. 1
figure 1

a The top panel shows the probabilistic membrane potential, distributions are indicated by box-and-whiskers plots at each time instant. The middle panel shows the first order Sobol indices, revealing the sources of the fluctuations in the membrane potential for each uncertain parameter. The black line is the sum of all indices; deviation from one indicates variance due to parameter interactions. The bottom panel shows the variance difference w.r.t. the model with 11 uncertain parameters for four parsimonious models. Only two parameters ( and ) are required to estimate a variability with less than 10 % error. b Average of the Sobol indices for all uncertain parameters during the time interval shown in a for the membrane potential and gating variables

In Fig. 1, the mean was obtained at the fourth level of a Sparse Grid constructed using the Gauss–Patterson (GP) quadrature rule, which requires 18591 function calls. These deterministic sequences are nested by construction, and for smooth integrands they can achieve the maximum degree of exactness from all nested rules [9, 18, 20].

In Fig. 1 neuronal variability is reported on the top panel by box-and-whisker plots, with one box per time point. On each box, the central mark is the median, the edges of the box are the 25th and 75th percentiles and the whiskers extend to the most extreme membrane potentials observed.

The middle panel of Fig. 1(a) shows the dynamics of the first order Sobol indices [18, 2124], computed by Sparse Grid quadrature. The Sobol indices in Fig. 1(a) allow one to determine the relative contribution of each uncertain parameter to the variability of the model output observed in the top panel at each instant in time. The black line shows the sum of all indices; deviation from one indicates fluctuations in the model output due to parameter interactions not accounted for by the first order indices.

The variance and Sobol indices were obtained at the third level of a GP Sparse Grid in 22 dimensions, requiring 17249 function calls. Notice that only 21 dimensions are needed to obtain the first order Sobol indices via the conditional variance. However, adding one extra dimension allows one to estimate the total variance with the same Sparse Grid, which often increases accuracy [18].

The RMS errors for the mean, variance and Sobol indices shown in Fig. 1 are summarised in Table 1. The most accurate numerical solution is taken as a reference for error analysis, thus the error estimates listed in the table are conservative, since they approximate the numerical error at the last but one quadrature level. Table 1 shows that the numerical errors are sufficiently small, validating the choice of quadrature level for further analysis.

Table 1 RMS error of the Sobol indices show in Fig. 1 (a) for all model outputs

In Fig. 1(b) all uncertainty sources are represented in a colour map, showing the average value of all Sobol indices during the time interval shown in Fig. 1(a). All outputs of the model, namely the membrane potential and the three gating variables are shown.

Interestingly, this colour map reveals that only a small subset of parameters causes most of the uncertainty in the model outputs. Therefore, all other parameters can be fixed to their nominal values preserving the probabilistic interpretation of the model. This is illustrated in the bottom panel of Fig. 1(a), which shows the difference in standard deviation of four effective models with respect to that of the model with the 11 uncertain parameters listed in Table 2. Notice that only two parameters, and , are required to estimate a variability with less than 10 % difference w.r.t. the original model for all time points.

Table 2 Normalised ranking of the average Sobol indices displayed in the colour map in Fig. 1 (b) using SVD

The singular value decomposition (SVD) of the data matrix M displayed in Fig. 1(b) as a colour map, namely , provides a quantitative ranking of the average Sobol indices. The first four columns of U after normalising the patterns so that each column has maximum entry 1 are shown in Table 2. Notice that the first mode (first column of U) contains more than 90 % of the energy, given by the square of the singular values. Notice also that in this maximum energy mode, two probabilistic degrees of freedom are dominant, namely and . The influence of all the other parameters in the variability of the model output is small. Therefore, they can be fixed to their nominal values. This greatly reduces the dimensionality of the model retaining most of its probabilistic features.

This feature of the H-H model is important, since it shows that probabilistic neurodynamic simulations can be performed parsimoniously, in an effective probability subspace of dimension far smaller than that of the original model. This greatly reduces the computational cost and facilitates the analysis of neuronal systems beyond individual cells.

As a non-trivial illustrative application of our parsimonious model Fig. 2 shows the impact of parameter fluctuations on two relevant neurocomputational properties, namely, neural response to current input and membrane refractoriness.

Fig. 2
figure 2

In panel a two current pulses are applied to a deterministic neuron with the nominal values for the parameters listed in Table 2. The first pulse is not strong enough to elicit a spike, the second stronger pulse immediately triggers a spike. Panel c shows the probabilistic counterpart of this experiment assuming 20 % of variability in the nominal parameters and . Large error bars show that the first small current pulse can trigger action potentials in this instance. In panels b and d the stronger current pulse is applied first, immediately eliciting a spike. In the deterministic model shown in panel b, the second small pulse fails to trigger a spike when applied during the refractory period. In the probabilistic model d, large variability in the output indicates a second action potential, which would not be expected from deterministic predictions

In Fig. 2(a) two 1 ms current pulses are applied to a deterministic H-H neuron whose membrane potential is initially at the resting potential. The first pulse of 8.5 μA/cm2 increases the membrane potential, but it does not elicit neuronal discharge. A second stronger pulse of 20 μA/cm2 applied later, immediately triggers a spike. The probabilistic counterpart of this experiment is shown in (c). The large error bars in response to the first pulse show that a small current pulse can indeed trigger action potentials contrary to deterministic predictions.

In panels (b) and (d) the strong current pulse of 20 μA/cm2 is applied first, immediately eliciting a spike. In the deterministic model shown in panel (b), a second small pulse of 10.5 μA/cm2 fails to ignite an action potential due to membrane refractoriness. Under parameter uncertainty, panel (d) shows large variability in the neuronal response, indicating a second action potential, which is impossible on deterministic grounds.

Synaptic input responses and membrane refractoriness are neurocomputational properties of realistic nerve cells, which directly influence collective neuronal behaviour as well as the efficiency of neural codes. Thus, the examples above illustrate how incorporating variability in neurodynamic models might be crucial to our understanding of the nervous system and the behaving brain.

Probabilistic generalisation of other deterministic models of great importance in mathematical neuroscience might be possible using the methods and results in this report, such as those describing neuronal interactions and signal transmission in active neural media, at a feasible computational cost.

References

  1. Fellous JM, Rudolph M, Destexhe A, Sejnowski TJ: Synaptic background noise controls the input/output characteristics of single cells in an in vitro model of in vivo activity. Neuroscience 2003,122(3):811-829.

    Article  Google Scholar 

  2. Destexhe A, Rudolph M, Fellous J-M, Sejnowski TJ: Fluctuating synaptic conductances recreate in vivo-like activity in neocortical neurons. Neuroscience 2001,107(1):13-24. 10.1016/S0306-4522(01)00344-X

    Article  Google Scholar 

  3. Paré D, Shink E, Gaudreau H, Destexhe A, Lang EJ: Impact of spontaneous synaptic activity on the resting properties of cat neocortical pyramidal neurons in vivo. J Neurophysiol 1998, 79: 1450-1460.

    Google Scholar 

  4. Destexhe A, Paré D: Impact of network activity on the integrative properties of neocortical pyramidal neurons in vivo. J Neurophysiol 1999,81(4):1531-1547.

    Google Scholar 

  5. Hô N, Destexhe A: Synaptic background activity enhances the responsiveness of neocortical pyramidal neurons. J Neurophysiol 2000, 84: 1488-1496.

    Google Scholar 

  6. Destexhe A, Rudolph-Lilith M Springer Series in Computational Neuroscience. In Neuronal Noise. Springer, Dordrecht; 2012.

    Chapter  Google Scholar 

  7. Nobile F, Tempone R, Webster CG: A sparse grid stochastic collocation method for partial differential equations with random input data. SIAM J Numer Anal 2008,46(5):2309-2345.

    Article  MathSciNet  MATH  Google Scholar 

  8. Bungartz H-J, Griebel M: Sparse grids. Acta Numer 2004, 13: 147-269. 10.1017/S0962492904000182

    Article  MathSciNet  MATH  Google Scholar 

  9. Gerstner T, Griebel M: Numerical integration using sparse grids. Numer Algorithms 1998, 18: 209-232.

    Article  MathSciNet  MATH  Google Scholar 

  10. Le Maître OP, Knio OM Scientific Computation. In Spectral Methods for Uncertainty Quantification: With Applications to Computational Fluid Dynamics. Springer, Dordrecht; 2010.

    Chapter  Google Scholar 

  11. Grigoriu M Springer Series in Reliability Engineering. In Stochastic Systems: Uncertainty Quantification and Propagation. Springer, New York; 2012.

    Chapter  Google Scholar 

  12. Xiu D: Numerical Methods for Stochastic Computations: A Spectral Method Approach. Princeton University Press, Princeton; 2010.

    MATH  Google Scholar 

  13. Ghanem RG, Spanos PD Civil, Mechanical and Other Engineering Series. In Stochastic Finite Elements: A Spectral Approach. Dover, Mineola; 2003.

    Google Scholar 

  14. Laing C, Zou Y, Smith B, Kevrekidis I: Managing heterogeneity in the study of neural oscillator dynamics. J Math Neurosci 2012., 2(1): Article ID 5 10.1186/2190-8567-2-5

    MATH  Google Scholar 

  15. Hodgkin AL, Huxley AF: A quantitative description of membrane current and its application to conduction and excitation in nerve. J Physiol 1952,117(4):500-544.

    Article  Google Scholar 

  16. Novak E, Ritter K: High dimensional integration of smooth functions over cubes. Numer Math 1996,75(1):79-97. 10.1007/s002110050231

    Article  MathSciNet  MATH  Google Scholar 

  17. Smolyak SA: Quadrature and interpolation formulas for tensor products of certain classes of functions. Dokl Akad Nauk SSSR 1963, 4: 240-243.

    MATH  Google Scholar 

  18. Torres Valderrama A, Blom J: Uncertainty propagation in neuronal dynamical systems. Technical report CWI LS-1304; 2012.

  19. Klimke A, Wohlmuth B: Algorithm 847: Spinterp: piecewise multilinear hierarchical sparse grid interpolation in MATLAB. ACM Trans Math Softw 2005,31(4):561-579. 10.1145/1114268.1114275

    Article  MathSciNet  MATH  Google Scholar 

  20. Patterson TNL: An algorithm for generating interpolatory quadrature rules of the highest degree of precision with preassigned nodes for general weight functions. ACM Trans Math Softw 1989,15(2):123-136. 10.1145/63522.63523

    Article  MathSciNet  MATH  Google Scholar 

  21. Homma T, Saltelli A: Importance measures in global sensitivity analysis of nonlinear models. Reliab Eng Syst Saf 1996, 52: 1-17.

    Article  Google Scholar 

  22. Ishigami T, Homma T: An importance quantification technique in uncertainty analysis for computer models. International Symposium on Uncertainity Modelling and Analysis (ISUMA’90) 1990. December 3–6, University of Maryland

    Google Scholar 

  23. Sobol IM: Sensitivity estimates for nonlinear mathematical models. Math Models 1990, 2: 112-118. (in Russian)

    MATH  Google Scholar 

  24. Saltelli A: Making best use of model evaluations to compute sensitivity indices. Comput Phys Commun 2002,145(2):280-297. 10.1016/S0010-4655(02)00280-1

    Article  MathSciNet  MATH  Google Scholar 

Download references

Acknowledgements

We thank Prof. Stan Gielen for useful discussions and suggestions. This work was supported by European Commission’s 7th Framework Program BioPreDyn grant number 289434. ATV gratefully acknowledges the financial support of the Exact Sciences division of the Netherlands Organisation for Scientific Research NWO grant number: veni 639.071.905.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Aldemar Torres Valderrama, Jeroen Witteveen, Maria Navarro or Joke Blom.

Additional information

Competing Interests

The authors declare that they have no competing interests.

Authors’ Contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0), which permits use, duplication, adaptation, distribution, and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Torres Valderrama, A., Witteveen, J., Navarro, M. et al. Uncertainty Propagation in Nerve Impulses Through the Action Potential Mechanism. J. Math. Neurosc. 5, 3 (2015). https://doi.org/10.1186/2190-8567-5-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/2190-8567-5-3

Keywords