@article{14901, abstract = {Global services like navigation, communication, and Earth observation have increased dramatically in the 21st century due to advances in outer space industries. But as orbits become increasingly crowded with both satellites and inevitable space debris pollution, continued operations become endangered by the heightened risks of debris collisions in orbit. Kessler Syndrome is the term for when a critical threshold of orbiting debris triggers a runaway positive feedback loop of debris collisions, creating debris congestion that can render orbits unusable. As this potential tipping point becomes more widely recognized, there have been renewed calls for debris mitigation and removal. Here, we combine complex systems and social-ecological systems approaches to study how these efforts may affect space debris accumulation and the likelihood of reaching Kessler Syndrome. Specifically, we model how debris levels are affected by future launch rates, cleanup activities, and collisions between extant debris. We contextualize and interpret our dynamic model within a discussion of existing space debris governance and other social, economic, and geopolitical factors that may influence effective collective management of the orbital commons. In line with previous studies, our model finds that debris congestion may be reached in less than 200 years, though a holistic management strategy combining removal and mitigation actions can avoid such outcomes while continuing space activities. Moreover, although active debris removal may be particularly effective, the current lack of market and governance support may impede its implementation. Research into these critical dynamics and the multi-faceted variables that influence debris outcomes can support policymakers in curating impactful governance strategies and realistic transition pathways to sustaining debris-free orbits. Overall, our study is useful for communicating about space debris sustainability in policy and education settings by providing an exploration of policy portfolio options supported by a simple and clear social-ecological modeling approach.}, author = {Nomura, Keiko and Rella, Simon and Merritt, Haily and Baltussen, Mathieu and Bird, Darcy and Tjuka, Annika and Falk, Dan}, issn = {1875-0281}, journal = {International Journal of the Commons}, keywords = {Sociology and Political Science}, number = {1}, publisher = {Ubiquity Press}, title = {{Tipping points of space debris in low earth orbit}}, doi = {10.5334/ijc.1275}, volume = {18}, year = {2024}, } @phdthesis{15020, abstract = {This thesis consists of four distinct pieces of work within theoretical biology, with two themes in common: the concept of optimization in biological systems, and the use of information-theoretic tools to quantify biological stochasticity and statistical uncertainty. Chapter 2 develops a statistical framework for studying biological systems which we believe to be optimized for a particular utility function, such as retinal neurons conveying information about visual stimuli. We formalize such beliefs as maximum-entropy Bayesian priors, constrained by the expected utility. We explore how such priors aid inference of system parameters with limited data and enable optimality hypothesis testing: is the utility higher than by chance? Chapter 3 examines the ultimate biological optimization process: evolution by natural selection. As some individuals survive and reproduce more successfully than others, populations evolve towards fitter genotypes and phenotypes. We formalize this as accumulation of genetic information, and use population genetics theory to study how much such information can be accumulated per generation and maintained in the face of random mutation and genetic drift. We identify the population size and fitness variance as the key quantities that control information accumulation and maintenance. Chapter 4 reuses the concept of genetic information from Chapter 3, but from a different perspective: we ask how much genetic information organisms actually need, in particular in the context of gene regulation. For example, how much information is needed to bind transcription factors at correct locations within the genome? Population genetics provides us with a refined answer: with an increasing population size, populations achieve higher fitness by maintaining more genetic information. Moreover, regulatory parameters experience selection pressure to optimize the fitness-information trade-off, i.e. minimize the information needed for a given fitness. This provides an evolutionary derivation of the optimization priors introduced in Chapter 2. Chapter 5 proves an upper bound on mutual information between a signal and a communication channel output (such as neural activity). Mutual information is an important utility measure for biological systems, but its practical use can be difficult due to the large dimensionality of many biological channels. Sometimes, a lower bound on mutual information is computed by replacing the high-dimensional channel outputs with decodes (signal estimates). Our result provides a corresponding upper bound, provided that the decodes are the maximum posterior estimates of the signal.}, author = {Hledik, Michal}, issn = {2663 - 337X}, keywords = {Theoretical biology, Optimality, Evolution, Information}, pages = {158}, publisher = {Institute of Science and Technology Austria}, title = {{Genetic information and biological optimization}}, doi = {10.15479/at:ista:15020}, year = {2024}, } @article{13127, abstract = {Cooperative disease defense emerges as group-level collective behavior, yet how group members make the underlying individual decisions is poorly understood. Using garden ants and fungal pathogens as an experimental model, we derive the rules governing individual ant grooming choices and show how they produce colony-level hygiene. Time-resolved behavioral analysis, pathogen quantification, and probabilistic modeling reveal that ants increase grooming and preferentially target highly-infectious individuals when perceiving high pathogen load, but transiently suppress grooming after having been groomed by nestmates. Ants thus react to both, the infectivity of others and the social feedback they receive on their own contagiousness. While inferred solely from momentary ant decisions, these behavioral rules quantitatively predict hour-long experimental dynamics, and synergistically combine into efficient colony-wide pathogen removal. Our analyses show that noisy individual decisions based on only local, incomplete, yet dynamically-updated information on pathogen threat and social feedback can lead to potent collective disease defense.}, author = {Casillas Perez, Barbara E and Bod'Ová, Katarína and Grasse, Anna V and Tkačik, Gašper and Cremer, Sylvia}, issn = {2041-1723}, journal = {Nature Communications}, publisher = {Springer Nature}, title = {{Dynamic pathogen detection and social feedback shape collective hygiene in ants}}, doi = {10.1038/s41467-023-38947-y}, volume = {14}, year = {2023}, } @article{12762, abstract = {Neurons in the brain are wired into adaptive networks that exhibit collective dynamics as diverse as scale-specific oscillations and scale-free neuronal avalanches. Although existing models account for oscillations and avalanches separately, they typically do not explain both phenomena, are too complex to analyze analytically or intractable to infer from data rigorously. Here we propose a feedback-driven Ising-like class of neural networks that captures avalanches and oscillations simultaneously and quantitatively. In the simplest yet fully microscopic model version, we can analytically compute the phase diagram and make direct contact with human brain resting-state activity recordings via tractable inference of the model’s two essential parameters. The inferred model quantitatively captures the dynamics over a broad range of scales, from single sensor oscillations to collective behaviors of extreme events and neuronal avalanches. Importantly, the inferred parameters indicate that the co-existence of scale-specific (oscillations) and scale-free (avalanches) dynamics occurs close to a non-equilibrium critical point at the onset of self-sustained oscillations.}, author = {Lombardi, Fabrizio and Pepic, Selver and Shriki, Oren and Tkačik, Gašper and De Martino, Daniele}, issn = {2662-8457}, journal = {Nature Computational Science}, pages = {254--263}, publisher = {Springer Nature}, title = {{Statistical modeling of adaptive neural networks explains co-existence of avalanches and oscillations in resting human brain}}, doi = {10.1038/s43588-023-00410-9}, volume = {3}, year = {2023}, } @article{14515, abstract = {Most natural and engineered information-processing systems transmit information via signals that vary in time. Computing the information transmission rate or the information encoded in the temporal characteristics of these signals requires the mutual information between the input and output signals as a function of time, i.e., between the input and output trajectories. Yet, this is notoriously difficult because of the high-dimensional nature of the trajectory space, and all existing techniques require approximations. We present an exact Monte Carlo technique called path weight sampling (PWS) that, for the first time, makes it possible to compute the mutual information between input and output trajectories for any stochastic system that is described by a master equation. The principal idea is to use the master equation to evaluate the exact conditional probability of an individual output trajectory for a given input trajectory and average this via Monte Carlo sampling in trajectory space to obtain the mutual information. We present three variants of PWS, which all generate the trajectories using the standard stochastic simulation algorithm. While direct PWS is a brute-force method, Rosenbluth-Rosenbluth PWS exploits the analogy between signal trajectory sampling and polymer sampling, and thermodynamic integration PWS is based on a reversible work calculation in trajectory space. PWS also makes it possible to compute the mutual information between input and output trajectories for systems with hidden internal states as well as systems with feedback from output to input. Applying PWS to the bacterial chemotaxis system, consisting of 182 coupled chemical reactions, demonstrates not only that the scheme is highly efficient but also that the number of receptor clusters is much smaller than hitherto believed, while their size is much larger.}, author = {Reinhardt, Manuel and Tkačik, Gašper and Ten Wolde, Pieter Rein}, issn = {2160-3308}, journal = {Physical Review X}, number = {4}, publisher = {American Physical Society}, title = {{Path weight sampling: Exact Monte Carlo computation of the mutual information between stochastic trajectories}}, doi = {10.1103/PhysRevX.13.041017}, volume = {13}, year = {2023}, } @article{14656, abstract = {Although much is known about how single neurons in the hippocampus represent an animal's position, how circuit interactions contribute to spatial coding is less well understood. Using a novel statistical estimator and theoretical modeling, both developed in the framework of maximum entropy models, we reveal highly structured CA1 cell-cell interactions in male rats during open field exploration. The statistics of these interactions depend on whether the animal is in a familiar or novel environment. In both conditions the circuit interactions optimize the encoding of spatial information, but for regimes that differ in the informativeness of their spatial inputs. This structure facilitates linear decodability, making the information easy to read out by downstream circuits. Overall, our findings suggest that the efficient coding hypothesis is not only applicable to individual neuron properties in the sensory periphery, but also to neural interactions in the central brain.}, author = {Nardin, Michele and Csicsvari, Jozsef L and Tkačik, Gašper and Savin, Cristina}, issn = {1529-2401}, journal = {The Journal of Neuroscience}, number = {48}, pages = {8140--8156}, publisher = {Society of Neuroscience}, title = {{The structure of hippocampal CA1 interactions optimizes spatial coding across experience}}, doi = {10.1523/JNEUROSCI.0194-23.2023}, volume = {43}, year = {2023}, } @article{12487, abstract = {Sleep plays a key role in preserving brain function, keeping the brain network in a state that ensures optimal computational capabilities. Empirical evidence indicates that such a state is consistent with criticality, where scale-free neuronal avalanches emerge. However, the relationship between sleep, emergent avalanches, and criticality remains poorly understood. Here we fully characterize the critical behavior of avalanches during sleep, and study their relationship with the sleep macro- and micro-architecture, in particular the cyclic alternating pattern (CAP). We show that avalanche size and duration distributions exhibit robust power laws with exponents approximately equal to −3/2 e −2, respectively. Importantly, we find that sizes scale as a power law of the durations, and that all critical exponents for neuronal avalanches obey robust scaling relations, which are consistent with the mean-field directed percolation universality class. Our analysis demonstrates that avalanche dynamics depends on the position within the NREM-REM cycles, with the avalanche density increasing in the descending phases and decreasing in the ascending phases of sleep cycles. Moreover, we show that, within NREM sleep, avalanche occurrence correlates with CAP activation phases, particularly A1, which are the expression of slow wave sleep propensity and have been proposed to be beneficial for cognitive processes. The results suggest that neuronal avalanches, and thus tuning to criticality, actively contribute to sleep development and play a role in preserving network function. Such findings, alongside characterization of the universality class for avalanches, open new avenues to the investigation of functional role of criticality during sleep with potential clinical application.Significance statementWe fully characterize the critical behavior of neuronal avalanches during sleep, and show that avalanches follow precise scaling laws that are consistent with the mean-field directed percolation universality class. The analysis provides first evidence of a functional relationship between avalanche occurrence, slow-wave sleep dynamics, sleep stage transitions and occurrence of CAP phase A during NREM sleep. Because CAP is considered one of the major guardians of NREM sleep that allows the brain to dynamically react to external perturbation and contributes to the cognitive consolidation processes occurring in sleep, our observations suggest that neuronal avalanches at criticality are associated with flexible response to external inputs and to cognitive processes, a key assumption of the critical brain hypothesis.}, author = {Scarpetta, Silvia and Morrisi, Niccolò and Mutti, Carlotta and Azzi, Nicoletta and Trippi, Irene and Ciliento, Rosario and Apicella, Ilenia and Messuti, Giovanni and Angiolelli, Marianna and Lombardi, Fabrizio and Parrino, Liborio and Vaudano, Anna Elisabetta}, issn = {2589-0042}, journal = {iScience}, number = {10}, pages = {107840}, publisher = {Elsevier}, title = {{Criticality of neuronal avalanches in human sleep and their relationship with sleep macro- and micro-architecture}}, doi = {10.1016/j.isci.2023.107840}, volume = {26}, year = {2023}, } @inproceedings{14862, author = {Rella, Simon and Kulikova, Y and Minnegalieva, Aygul and Kondrashov, Fyodor}, booktitle = {European Journal of Public Health}, issn = {1464-360X}, keywords = {Public Health, Environmental and Occupational Health}, number = {Supplement_2}, publisher = {Oxford University Press}, title = {{Complex vaccination strategies prevent the emergence of vaccine resistance}}, doi = {10.1093/eurpub/ckad160.597}, volume = {33}, year = {2023}, } @article{14402, abstract = {Alpha oscillations are a distinctive feature of the awake resting state of the human brain. However, their functional role in resting-state neuronal dynamics remains poorly understood. Here we show that, during resting wakefulness, alpha oscillations drive an alternation of attenuation and amplification bouts in neural activity. Our analysis indicates that inhibition is activated in pulses that last for a single alpha cycle and gradually suppress neural activity, while excitation is successively enhanced over a few alpha cycles to amplify neural activity. Furthermore, we show that long-term alpha amplitude fluctuations—the “waxing and waning” phenomenon—are an attenuation-amplification mechanism described by a power-law decay of the activity rate in the “waning” phase. Importantly, we do not observe such dynamics during non-rapid eye movement (NREM) sleep with marginal alpha oscillations. The results suggest that alpha oscillations modulate neural activity not only through pulses of inhibition (pulsed inhibition hypothesis) but also by timely enhancement of excitation (or disinhibition).}, author = {Lombardi, Fabrizio and Herrmann, Hans J. and Parrino, Liborio and Plenz, Dietmar and Scarpetta, Silvia and Vaudano, Anna Elisabetta and De Arcangelis, Lucilla and Shriki, Oren}, issn = {2211-1247}, journal = {Cell Reports}, number = {10}, publisher = {Elsevier}, title = {{Beyond pulsed inhibition: Alpha oscillations modulate attenuation and amplification of neural activity in the awake resting state}}, doi = {10.1016/j.celrep.2023.113162}, volume = {42}, year = {2023}, } @unpublished{10821, abstract = {Rhythmical cortical activity has long been recognized as a pillar in the architecture of brain functions. Yet, the dynamic organization of its underlying neuronal population activity remains elusive. Here we uncover a unique organizational principle regulating collective neural dynamics associated with the alpha rhythm in the awake resting-state. We demonstrate that cascades of neural activity obey attenuation-amplification dynamics (AAD), with a transition from the attenuation regime—within alpha cycles—to the amplification regime—across a few alpha cycles—that correlates with the characteristic frequency of the alpha rhythm. We find that this short-term AAD is part of a large-scale, size-dependent temporal structure of neural cascades that obeys the Omori law: Following large cascades, smaller cascades occur at a rate that decays as a power-law of the time elapsed from such events—a long-term AAD regulating brain activity over the timescale of seconds. We show that such an organization corresponds to the "waxing and waning" of the alpha rhythm. Importantly, we observe that short- and long-term AAD are unique to the awake resting-state, being absent during NREM sleep. These results provide a quantitative, dynamical description of the so-far-qualitative notion of the "waxing and waning" phenomenon, and suggest the AAD as a key principle governing resting-state dynamics across timescales.}, author = {Lombardi, Fabrizio and Herrmann, Hans J. and Parrino, Liborio and Plenz, Dietmar and Scarpetta, Silvia and Vaudano, Anna Elisabetta and de Arcangelis, Lucilla and Shriki, Oren}, booktitle = {bioRxiv}, pages = {25}, publisher = {Cold Spring Harbor Laboratory}, title = {{Alpha rhythm induces attenuation-amplification dynamics in neural activity cascades}}, doi = {10.1101/2022.03.03.482657}, year = {2022}, } @article{11638, abstract = {Statistical inference is central to many scientific endeavors, yet how it works remains unresolved. Answering this requires a quantitative understanding of the intrinsic interplay between statistical models, inference methods, and the structure in the data. To this end, we characterize the efficacy of direct coupling analysis (DCA)—a highly successful method for analyzing amino acid sequence data—in inferring pairwise interactions from samples of ferromagnetic Ising models on random graphs. Our approach allows for physically motivated exploration of qualitatively distinct data regimes separated by phase transitions. We show that inference quality depends strongly on the nature of data-generating distributions: optimal accuracy occurs at an intermediate temperature where the detrimental effects from macroscopic order and thermal noise are minimal. Importantly our results indicate that DCA does not always outperform its local-statistics-based predecessors; while DCA excels at low temperatures, it becomes inferior to simple correlation thresholding at virtually all temperatures when data are limited. Our findings offer insights into the regime in which DCA operates so successfully, and more broadly, how inference interacts with the structure in the data.}, author = {Ngampruetikorn, Vudtiwat and Sachdeva, Vedant and Torrence, Johanna and Humplik, Jan and Schwab, David J. and Palmer, Stephanie E.}, issn = {2643-1564}, journal = {Physical Review Research}, number = {2}, publisher = {American Physical Society}, title = {{Inferring couplings in networks across order-disorder phase transitions}}, doi = {10.1103/PhysRevResearch.4.023240}, volume = {4}, year = {2022}, } @article{12156, abstract = {Models of transcriptional regulation that assume equilibrium binding of transcription factors have been less successful at predicting gene expression from sequence in eukaryotes than in bacteria. This could be due to the non-equilibrium nature of eukaryotic regulation. Unfortunately, the space of possible non-equilibrium mechanisms is vast and predominantly uninteresting. The key question is therefore how this space can be navigated efficiently, to focus on mechanisms and models that are biologically relevant. In this review, we advocate for the normative role of theory—theory that prescribes rather than just describes—in providing such a focus. Theory should expand its remit beyond inferring mechanistic models from data, towards identifying non-equilibrium gene regulatory schemes that may have been evolutionarily selected, despite their energy consumption, because they are precise, reliable, fast, or otherwise outperform regulation at equilibrium. We illustrate our reasoning by toy examples for which we provide simulation code.}, author = {Zoller, Benjamin and Gregor, Thomas and Tkačik, Gašper}, issn = {2452-3100}, journal = {Current Opinion in Systems Biology}, keywords = {Applied Mathematics, Computer Science Applications, Drug Discovery, General Biochemistry, Genetics and Molecular Biology, Modeling and Simulation}, number = {9}, publisher = {Elsevier}, title = {{Eukaryotic gene regulation at equilibrium, or non?}}, doi = {10.1016/j.coisb.2022.100435}, volume = {31}, year = {2022}, } @article{10530, abstract = {Cell dispersion from a confined area is fundamental in a number of biological processes, including cancer metastasis. To date, a quantitative understanding of the interplay of single cell motility, cell proliferation, and intercellular contacts remains elusive. In particular, the role of E- and N-Cadherin junctions, central components of intercellular contacts, is still controversial. Combining theoretical modeling with in vitro observations, we investigate the collective spreading behavior of colonies of human cancer cells (T24). The spreading of these colonies is driven by stochastic single-cell migration with frequent transient cell-cell contacts. We find that inhibition of E- and N-Cadherin junctions decreases colony spreading and average spreading velocities, without affecting the strength of correlations in spreading velocities of neighboring cells. Based on a biophysical simulation model for cell migration, we show that the behavioral changes upon disruption of these junctions can be explained by reduced repulsive excluded volume interactions between cells. This suggests that in cancer cell migration, cadherin-based intercellular contacts sharpen cell boundaries leading to repulsive rather than cohesive interactions between cells, thereby promoting efficient cell spreading during collective migration. }, author = {Zisis, Themistoklis and Brückner, David and Brandstätter, Tom and Siow, Wei Xiong and d’Alessandro, Joseph and Vollmar, Angelika M. and Broedersz, Chase P. and Zahler, Stefan}, issn = {0006-3495}, journal = {Biophysical Journal}, keywords = {Biophysics}, number = {1}, pages = {P44--60}, publisher = {Elsevier}, title = {{Disentangling cadherin-mediated cell-cell interactions in collective cancer cell migration}}, doi = {10.1016/j.bpj.2021.12.006}, volume = {121}, year = {2022}, } @article{10736, abstract = {Predicting function from sequence is a central problem of biology. Currently, this is possible only locally in a narrow mutational neighborhood around a wildtype sequence rather than globally from any sequence. Using random mutant libraries, we developed a biophysical model that accounts for multiple features of σ70 binding bacterial promoters to predict constitutive gene expression levels from any sequence. We experimentally and theoretically estimated that 10–20% of random sequences lead to expression and ~80% of non-expressing sequences are one mutation away from a functional promoter. The potential for generating expression from random sequences is so pervasive that selection acts against σ70-RNA polymerase binding sites even within inter-genic, promoter-containing regions. This pervasiveness of σ70-binding sites implies that emergence of promoters is not the limiting step in gene regulatory evolution. Ultimately, the inclusion of novel features of promoter function into a mechanistic model enabled not only more accurate predictions of gene expression levels, but also identified that promoters evolve more rapidly than previously thought.}, author = {Lagator, Mato and Sarikas, Srdjan and Steinrueck, Magdalena and Toledo-Aparicio, David and Bollback, Jonathan P and Guet, Calin C and Tkačik, Gašper}, issn = {2050-084X}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{Predicting bacterial promoter function and evolution from random sequences}}, doi = {10.7554/eLife.64543}, volume = {11}, year = {2022}, } @article{12332, abstract = {Activity of sensory neurons is driven not only by external stimuli but also by feedback signals from higher brain areas. Attention is one particularly important internal signal whose presumed role is to modulate sensory representations such that they only encode information currently relevant to the organism at minimal cost. This hypothesis has, however, not yet been expressed in a normative computational framework. Here, by building on normative principles of probabilistic inference and efficient coding, we developed a model of dynamic population coding in the visual cortex. By continuously adapting the sensory code to changing demands of the perceptual observer, an attention-like modulation emerges. This modulation can dramatically reduce the amount of neural activity without deteriorating the accuracy of task-specific inferences. Our results suggest that a range of seemingly disparate cortical phenomena such as intrinsic gain modulation, attention-related tuning modulation, and response variability could be manifestations of the same underlying principles, which combine efficient sensory coding with optimal probabilistic inference in dynamic environments.}, author = {Mlynarski, Wiktor F and Tkačik, Gašper}, issn = {1545-7885}, journal = {PLoS Biology}, number = {12}, pages = {e3001889}, publisher = {Public Library of Science}, title = {{Efficient coding theory of dynamic attentional modulation}}, doi = {10.1371/journal.pbio.3001889}, volume = {20}, year = {2022}, } @article{12081, abstract = {Selection accumulates information in the genome—it guides stochastically evolving populations toward states (genotype frequencies) that would be unlikely under neutrality. This can be quantified as the Kullback–Leibler (KL) divergence between the actual distribution of genotype frequencies and the corresponding neutral distribution. First, we show that this population-level information sets an upper bound on the information at the level of genotype and phenotype, limiting how precisely they can be specified by selection. Next, we study how the accumulation and maintenance of information is limited by the cost of selection, measured as the genetic load or the relative fitness variance, both of which we connect to the control-theoretic KL cost of control. The information accumulation rate is upper bounded by the population size times the cost of selection. This bound is very general, and applies across models (Wright–Fisher, Moran, diffusion) and to arbitrary forms of selection, mutation, and recombination. Finally, the cost of maintaining information depends on how it is encoded: Specifying a single allele out of two is expensive, but one bit encoded among many weakly specified loci (as in a polygenic trait) is cheap.}, author = {Hledik, Michal and Barton, Nicholas H and Tkačik, Gašper}, issn = {1091-6490}, journal = {Proceedings of the National Academy of Sciences}, number = {36}, publisher = {Proceedings of the National Academy of Sciences}, title = {{Accumulation and maintenance of information in evolution}}, doi = {10.1073/pnas.2123152119}, volume = {119}, year = {2022}, } @article{10535, abstract = {Realistic models of biological processes typically involve interacting components on multiple scales, driven by changing environment and inherent stochasticity. Such models are often analytically and numerically intractable. We revisit a dynamic maximum entropy method that combines a static maximum entropy with a quasi-stationary approximation. This allows us to reduce stochastic non-equilibrium dynamics expressed by the Fokker-Planck equation to a simpler low-dimensional deterministic dynamics, without the need to track microscopic details. Although the method has been previously applied to a few (rather complicated) applications in population genetics, our main goal here is to explain and to better understand how the method works. We demonstrate the usefulness of the method for two widely studied stochastic problems, highlighting its accuracy in capturing important macroscopic quantities even in rapidly changing non-stationary conditions. For the Ornstein-Uhlenbeck process, the method recovers the exact dynamics whilst for a stochastic island model with migration from other habitats, the approximation retains high macroscopic accuracy under a wide range of scenarios in a dynamic environment.}, author = {Bod'ová, Katarína and Szep, Eniko and Barton, Nicholas H}, issn = {1553-7358}, journal = {PLoS Computational Biology}, number = {12}, publisher = {Public Library of Science}, title = {{Dynamic maximum entropy provides accurate approximation of structured population dynamics}}, doi = {10.1371/journal.pcbi.1009661}, volume = {17}, year = {2021}, } @unpublished{10912, abstract = {Brain dynamics display collective phenomena as diverse as neuronal oscillations and avalanches. Oscillations are rhythmic, with fluctuations occurring at a characteristic scale, whereas avalanches are scale-free cascades of neural activity. Here we show that such antithetic features can coexist in a very generic class of adaptive neural networks. In the most simple yet fully microscopic model from this class we make direct contact with human brain resting-state activity recordings via tractable inference of the model's two essential parameters. The inferred model quantitatively captures the dynamics over a broad range of scales, from single sensor fluctuations, collective behaviors of nearly-synchronous extreme events on multiple sensors, to neuronal avalanches unfolding over multiple sensors across multiple time-bins. Importantly, the inferred parameters correlate with model-independent signatures of "closeness to criticality", suggesting that the coexistence of scale-specific (neural oscillations) and scale-free (neuronal avalanches) dynamics in brain activity occurs close to a non-equilibrium critical point at the onset of self-sustained oscillations.}, author = {Lombardi, Fabrizio and Pepic, Selver and Shriki, Oren and Tkačik, Gašper and De Martino, Daniele}, pages = {37}, publisher = {arXiv}, title = {{Quantifying the coexistence of neuronal oscillations and avalanches}}, doi = {10.48550/ARXIV.2108.06686}, year = {2021}, } @unpublished{10579, abstract = {We consider a totally asymmetric simple exclusion process (TASEP) consisting of particles on a lattice that require binding by a "token" to move. Using a combination of theory and simulations, we address the following questions: (i) How token binding kinetics affects the current-density relation; (ii) How the current-density relation depends on the scarcity of tokens; (iii) How tokens propagate the effects of the locally-imposed disorder (such a slow site) over the entire lattice; (iv) How a shared pool of tokens couples concurrent TASEPs running on multiple lattices; (v) How our results translate to TASEPs with open boundaries that exchange particles with the reservoir. Since real particle motion (including in systems that inspired the standard TASEP model, e.g., protein synthesis or movement of molecular motors) is often catalyzed, regulated, actuated, or otherwise mediated, the token-driven TASEP dynamics analyzed in this paper should allow for a better understanding of real systems and enable a closer match between TASEP theory and experimental observations.}, author = {Kavcic, Bor and Tkačik, Gašper}, booktitle = {arXiv}, title = {{Token-driven totally asymmetric simple exclusion process}}, doi = {10.48550/arXiv.2112.13558}, year = {2021}, } @article{7463, abstract = {Resting-state brain activity is characterized by the presence of neuronal avalanches showing absence of characteristic size. Such evidence has been interpreted in the context of criticality and associated with the normal functioning of the brain. A distinctive attribute of systems at criticality is the presence of long-range correlations. Thus, to verify the hypothesis that the brain operates close to a critical point and consequently assess deviations from criticality for diagnostic purposes, it is of primary importance to robustly and reliably characterize correlations in resting-state brain activity. Recent works focused on the analysis of narrow-band electroencephalography (EEG) and magnetoencephalography (MEG) signal amplitude envelope, showing evidence of long-range temporal correlations (LRTC) in neural oscillations. However, brain activity is a broadband phenomenon, and a significant piece of information useful to precisely discriminate between normal (critical) and pathological behavior (non-critical), may be encoded in the broadband spatio-temporal cortical dynamics. Here we propose to characterize the temporal correlations in the broadband brain activity through the lens of neuronal avalanches. To this end, we consider resting-state EEG and long-term MEG recordings, extract the corresponding neuronal avalanche sequences, and study their temporal correlations. We demonstrate that the broadband resting-state brain activity consistently exhibits long-range power-law correlations in both EEG and MEG recordings, with similar values of the scaling exponents. Importantly, although we observe that the avalanche size distribution depends on scale parameters, scaling exponents characterizing long-range correlations are quite robust. In particular, they are independent of the temporal binning (scale of analysis), indicating that our analysis captures intrinsic characteristics of the underlying dynamics. Because neuronal avalanches constitute a fundamental feature of neural systems with universal characteristics, the proposed approach may serve as a general, systems- and experiment-independent procedure to infer the existence of underlying long-range correlations in extended neural systems, and identify pathological behaviors in the complex spatio-temporal interplay of cortical rhythms.}, author = {Lombardi, Fabrizio and Shriki, Oren and Herrmann, Hans J and de Arcangelis, Lucilla}, issn = {1872-8286}, journal = {Neurocomputing}, pages = {657--666}, publisher = {Elsevier}, title = {{Long-range temporal correlations in the broadband resting state activity of the human brain revealed by neuronal avalanches}}, doi = {10.1016/j.neucom.2020.05.126}, volume = {461}, year = {2021}, }