1
|
Ham L, Coomer MA, Öcal K, Grima R, Stumpf MPH. A stochastic vs deterministic perspective on the timing of cellular events. Nat Commun 2024; 15:5286. [PMID: 38902228 PMCID: PMC11190182 DOI: 10.1038/s41467-024-49624-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Accepted: 06/12/2024] [Indexed: 06/22/2024] Open
Abstract
Cells are the fundamental units of life, and like all life forms, they change over time. Changes in cell state are driven by molecular processes; of these many are initiated when molecule numbers reach and exceed specific thresholds, a characteristic that can be described as "digital cellular logic". Here we show how molecular and cellular noise profoundly influence the time to cross a critical threshold-the first-passage time-and map out scenarios in which stochastic dynamics result in shorter or longer average first-passage times compared to noise-less dynamics. We illustrate the dependence of the mean first-passage time on noise for a set of exemplar models of gene expression, auto-regulatory feedback control, and enzyme-mediated catalysis. Our theory provides intuitive insight into the origin of these effects and underscores two important insights: (i) deterministic predictions for cellular event timing can be highly inaccurate when molecule numbers are within the range known for many cells; (ii) molecular noise can significantly shift mean first-passage times, particularly within auto-regulatory genetic feedback circuits.
Collapse
Affiliation(s)
- Lucy Ham
- School of BioSciences, University of Melbourne, Parkville, Australia
- School of Mathematics and Statistics, University of Melbourne, Parkville, Australia
| | - Megan A Coomer
- School of BioSciences, University of Melbourne, Parkville, Australia
- School of Mathematics and Statistics, University of Melbourne, Parkville, Australia
| | - Kaan Öcal
- School of Informatics, University of Edinburgh, Edinburgh, UK
- School of BioSciences, University of Melbourne, Parkville, Australia
| | - Ramon Grima
- School of Biological Sciences, University of Edinburgh, Edinburgh, UK
| | - Michael P H Stumpf
- School of BioSciences, University of Melbourne, Parkville, Australia.
- School of Mathematics and Statistics, University of Melbourne, Parkville, Australia.
| |
Collapse
|
2
|
Martínez N, Deza RR, Montani F. Characterizing the information transmission of inverse stochastic resonance and noise-induced activity amplification in neuronal systems. Phys Rev E 2023; 107:054402. [PMID: 37329070 DOI: 10.1103/physreve.107.054402] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2022] [Accepted: 04/13/2023] [Indexed: 06/18/2023]
Abstract
Purkinje cells exhibit a reduction of the mean firing rate at intermediate-noise intensities, which is somewhat reminiscent of the response enhancement known as "stochastic resonance" (SR). Although the comparison with the stochastic resonance ends here, the current phenomenon has been given the name "inverse stochastic resonance" (ISR). Recent research has demonstrated that the ISR effect, like its close relative "nonstandard SR" [or, more correctly, noise-induced activity amplification (NIAA)], has been shown to stem from the weak-noise quenching of the initial distribution, in bistable regimes where the metastable state has a larger attraction basin than the global minimum. To understand the underlying mechanism of the ISR and NIAA phenomena, we study the probability distribution function of a one-dimensional system subjected to a bistable potential that has the property of symmetry, i.e., if we change the sign of one of its parameters, we can obtain both phenomena with the same properties in the depth of the wells and the width of their basins of attraction subjected to Gaussian white noise with variable intensity. Previous work has shown that one can theoretically determine the probability distribution function using the convex sum between the behavior at small and high noise intensities. To determine the probability distribution function more precisely, we resort to the "weighted ensemble Brownian dynamics simulation" model, which provides an accurate estimate of the probability distribution function for both low and high noise intensities and, most importantly, for the transition of both behaviors. In this way, on the one hand, we show that both phenomena emerge from a metastable system where, in the case of ISR, the global minimum of the system is in a state of lower activity, while in the case of NIAA, the global minimum is in a state of increased activity, the importance of which does not depend on the width of the basins of attraction. On the other hand, we see that quantifiers such as Fisher information, statistical complexity, and especially Shannon entropy fail to distinguish them, but they show the existence of the mentioned phenomena. Thus, noise management may well be a mechanism by which Purkinje cells find an efficient way to transmit information in the cerebral cortex.
Collapse
Affiliation(s)
- Nataniel Martínez
- IFIMAR (CONICET), Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Mar del Plata, B7602AYL Mar del Plata, Argentina
| | - Roberto R Deza
- IFIMAR (CONICET), Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Mar del Plata, B7602AYL Mar del Plata, Argentina
| | - Fernando Montani
- IFLP (CONICET), Facultad de Ciencias Exactas, Universidad Nacional de La Plata, B1900 La Plata, Argentina
| |
Collapse
|
3
|
Ampuero F, Hase MO. First-passage process in degree space for the time-dependent Erdős-Rényi and Watts-Strogatz models. Phys Rev E 2022; 106:034320. [PMID: 36266810 DOI: 10.1103/physreve.106.034320] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2022] [Accepted: 09/07/2022] [Indexed: 06/16/2023]
Abstract
In this work, we investigate the temporal evolution of the degree of a given vertex in a network by mapping the dynamics into a random walk problem in degree space. We analyze when the degree approximates a preestablished value through a parallel with the first-passage problem of random walks. The method is illustrated on the time-dependent versions of the Erdős-Rényi and Watts-Strogatz models, which were originally formulated as static networks. We have succeeded in obtaining an analytic form for the first and the second moments of the first-passage time and showing how they depend on the size of the network. The dominant contribution for large networks with N vertices indicates that these quantities scale on the ratio N/p, where p is the linking probability.
Collapse
Affiliation(s)
- F Ampuero
- Escola de Artes, Ciências e Humanidades, Universidade de São Paulo, Av. Arlindo Béttio 1000, 03828-000 São Paulo, Brazil
| | - M O Hase
- Escola de Artes, Ciências e Humanidades, Universidade de São Paulo, Av. Arlindo Béttio 1000, 03828-000 São Paulo, Brazil
| |
Collapse
|
4
|
Mankin R, Rekker A. Effects of transient subordinators on the firing statistics of a neuron model driven by dichotomous noise. Phys Rev E 2020; 102:012103. [PMID: 32794976 DOI: 10.1103/physreve.102.012103] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2020] [Accepted: 06/15/2020] [Indexed: 06/11/2023]
Abstract
The behavior of a stochastic perfect integrate-and-fire (PIF) model of neurons is considered. The effect of temporally correlated random activity of synaptic inputs is modeled as a combination of an asymmetric dichotomous noise and a random operation time in the form of an inverse strictly increasing Lévy-type subordinator. Using a first-passage-time formulation, we find exact expressions for the output interspike interval (ISI) statistics. Particularly, it is shown that at some parameter regimes the ISI density exhibits a multimodal structure. Moreover, it is demonstrated that the coefficient of variation, the serial correlation coefficient, and the Fano factor display a nonmonotonic dependence on the mean input current μ, i.e., the ISI's regularity is maximized at an intermediate value of μ. The features of spike statistics, analytically revealed in our study, are compared with previously obtained results for a perfect integrate-and-fire neuron model driven by dichotomous noise (without subordination).
Collapse
Affiliation(s)
- Romi Mankin
- School of Natural Sciences and Health, Tallinn University, 29 Narva Road, 10120 Tallinn, Estonia
| | - Astrid Rekker
- School of Natural Sciences and Health, Tallinn University, 29 Narva Road, 10120 Tallinn, Estonia
| |
Collapse
|
5
|
Braun W, Longtin A. Interspike interval correlations in networks of inhibitory integrate-and-fire neurons. Phys Rev E 2019; 99:032402. [PMID: 30999498 DOI: 10.1103/physreve.99.032402] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Indexed: 11/07/2022]
Abstract
We study temporal correlations of interspike intervals, quantified by the network-averaged serial correlation coefficient (SCC), in networks of both current- and conductance-based purely inhibitory integrate-and-fire neurons. Numerical simulations reveal transitions to negative SCCs at intermediate values of bias current drive and network size. As bias drive and network size are increased past these values, the SCC returns to zero. The SCC is maximally negative at an intermediate value of the network oscillation strength. The dependence of the SCC on two canonical schemes for synaptic connectivity is studied, and it is shown that the results occur robustly in both schemes. For conductance-based synapses, the SCC becomes negative at the onset of both a fast and slow coherent network oscillation. We then show by means of offline simulations using prerecorded network activity that a neuron's SCC is highly sensitive to its number of presynaptic inputs. Finally, we devise a noise-reduced diffusion approximation for current-based networks that accounts for the observed temporal correlation transitions.
Collapse
Affiliation(s)
- Wilhelm Braun
- Neural Network Dynamics and Computation, Institut für Genetik, Universität Bonn, Kirschallee 1, 53115 Bonn, Germany.,Department of Physics and Centre for Neural Dynamics, University of Ottawa, 598 King Edward, Ottawa K1N 6N5, Canada
| | - André Longtin
- Department of Physics and Centre for Neural Dynamics, University of Ottawa, 598 King Edward, Ottawa K1N 6N5, Canada
| |
Collapse
|
6
|
Jedynak M, Pons AJ, Garcia-Ojalvo J, Goodfellow M. Temporally correlated fluctuations drive epileptiform dynamics. Neuroimage 2017; 146:188-196. [PMID: 27865920 PMCID: PMC5353705 DOI: 10.1016/j.neuroimage.2016.11.034] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2016] [Revised: 10/17/2016] [Accepted: 11/13/2016] [Indexed: 12/27/2022] Open
Abstract
Macroscopic models of brain networks typically incorporate assumptions regarding the characteristics of afferent noise, which is used to represent input from distal brain regions or ongoing fluctuations in non-modelled parts of the brain. Such inputs are often modelled by Gaussian white noise which has a flat power spectrum. In contrast, macroscopic fluctuations in the brain typically follow a 1/fb spectrum. It is therefore important to understand the effect on brain dynamics of deviations from the assumption of white noise. In particular, we wish to understand the role that noise might play in eliciting aberrant rhythms in the epileptic brain. To address this question we study the response of a neural mass model to driving by stochastic, temporally correlated input. We characterise the model in terms of whether it generates "healthy" or "epileptiform" dynamics and observe which of these dynamics predominate under different choices of temporal correlation and amplitude of an Ornstein-Uhlenbeck process. We find that certain temporal correlations are prone to eliciting epileptiform dynamics, and that these correlations produce noise with maximal power in the δ and θ bands. Crucially, these are rhythms that are found to be enhanced prior to seizures in humans and animal models of epilepsy. In order to understand why these rhythms can generate epileptiform dynamics, we analyse the response of the model to sinusoidal driving and explain how the bifurcation structure of the model gives rise to these findings. Our results provide insight into how ongoing fluctuations in brain dynamics can facilitate the onset and propagation of epileptiform rhythms in brain networks. Furthermore, we highlight the need to combine large-scale models with noise of a variety of different types in order to understand brain (dys-)function.
Collapse
Affiliation(s)
- Maciej Jedynak
- Departament de Física i Enginyeria Nuclear, Universitat Politècnica de Catalunya, Terrassa, Spain; Department of Experimental and Health Sciences, Universitat Pompeu Fabra, Parc de Recerca Biomèdica de Barcelona, Barcelona, Spain.
| | - Antonio J Pons
- Departament de Física i Enginyeria Nuclear, Universitat Politècnica de Catalunya, Terrassa, Spain
| | - Jordi Garcia-Ojalvo
- Department of Experimental and Health Sciences, Universitat Pompeu Fabra, Parc de Recerca Biomèdica de Barcelona, Barcelona, Spain
| | - Marc Goodfellow
- College of Engineering, Mathematics and Physical Sciences, University of Exeter, Exeter, UK; Centre for Biomedical Modelling and Analysis, University of Exeter, Exeter, UK; EPSRC Centre for Predictive Modelling in Healthcare, University of Exeter, Exeter, UK
| |
Collapse
|
7
|
Braun W, Thul R. Sign changes as a universal concept in first-passage-time calculations. Phys Rev E 2017; 95:012114. [PMID: 28208500 DOI: 10.1103/physreve.95.012114] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2016] [Indexed: 06/06/2023]
Abstract
First-passage-time problems are ubiquitous across many fields of study, including transport processes in semiconductors and biological synapses, evolutionary game theory and percolation. Despite their prominence, first-passage-time calculations have proven to be particularly challenging. Analytical results to date have often been obtained under strong conditions, leaving most of the exploration of first-passage-time problems to direct numerical computations. Here we present an analytical approach that allows the derivation of first-passage-time distributions for the wide class of nondifferentiable Gaussian processes. We demonstrate that the concept of sign changes naturally generalizes the common practice of counting crossings to determine first-passage events. Our method works across a wide range of time-dependent boundaries and noise strengths, thus alleviating common hurdles in first-passage-time calculations.
Collapse
Affiliation(s)
- Wilhelm Braun
- Department of Physics and Centre for Neural Dynamics, University of Ottawa, Ottawa, Canada K1N 6N5
| | - Rüdiger Thul
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, NG7 2RD, United Kingdom
| |
Collapse
|
8
|
Thul R, Coombes S, Laing CR. Neural Field Models with Threshold Noise. JOURNAL OF MATHEMATICAL NEUROSCIENCE 2016; 6:3. [PMID: 26936267 PMCID: PMC4775726 DOI: 10.1186/s13408-016-0035-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/26/2015] [Accepted: 02/19/2016] [Indexed: 06/05/2023]
Abstract
The original neural field model of Wilson and Cowan is often interpreted as the averaged behaviour of a network of switch like neural elements with a distribution of switch thresholds, giving rise to the classic sigmoidal population firing-rate function so prevalent in large scale neuronal modelling. In this paper we explore the effects of such threshold noise without recourse to averaging and show that spatial correlations can have a strong effect on the behaviour of waves and patterns in continuum models. Moreover, for a prescribed spatial covariance function we explore the differences in behaviour that can emerge when the underlying stationary distribution is changed from Gaussian to non-Gaussian. For travelling front solutions, in a system with exponentially decaying spatial interactions, we make use of an interface approach to calculate the instantaneous wave speed analytically as a series expansion in the noise strength. From this we find that, for weak noise, the spatially averaged speed depends only on the choice of covariance function and not on the shape of the stationary distribution. For a system with a Mexican-hat spatial connectivity we further find that noise can induce localised bump solutions, and using an interface stability argument show that there can be multiple stable solution branches.
Collapse
Affiliation(s)
- Rüdiger Thul
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, University Park, Nottingham, NG7 2RD, UK.
| | - Stephen Coombes
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, University Park, Nottingham, NG7 2RD, UK.
| | - Carlo R Laing
- Institute of Natural and Mathematical Sciences, Massey University (Albany), Private Bag 102-904, North Shore Mail Centre, Auckland, New Zealand.
| |
Collapse
|