@article{7490, abstract = {In plants, clathrin mediated endocytosis (CME) represents the major route for cargo internalisation from the cell surface. It has been assumed to operate in an evolutionary conserved manner as in yeast and animals. Here we report characterisation of ultrastructure, dynamics and mechanisms of plant CME as allowed by our advancement in electron microscopy and quantitative live imaging techniques. Arabidopsis CME appears to follow the constant curvature model and the bona fide CME population generates vesicles of a predominantly hexagonal-basket type; larger and with faster kinetics than in other models. Contrary to the existing paradigm, actin is dispensable for CME events at the plasma membrane but plays a unique role in collecting endocytic vesicles, sorting of internalised cargos and directional endosome movement that itself actively promote CME events. Internalized vesicles display a strongly delayed and sequential uncoating. These unique features highlight the independent evolution of the plant CME mechanism during the autonomous rise of multicellularity in eukaryotes.}, author = {Narasimhan, Madhumitha and Johnson, Alexander J and Prizak, Roshan and Kaufmann, Walter and Tan, Shutang and Casillas Perez, Barbara E and Friml, Jiří}, issn = {2050-084X}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{Evolutionarily unique mechanistic framework of clathrin-mediated endocytosis in plants}}, doi = {10.7554/eLife.52067}, volume = {9}, year = {2020}, } @article{7569, abstract = {Genes differ in the frequency at which they are expressed and in the form of regulation used to control their activity. In particular, positive or negative regulation can lead to activation of a gene in response to an external signal. Previous works proposed that the form of regulation of a gene correlates with its frequency of usage: positive regulation when the gene is frequently expressed and negative regulation when infrequently expressed. Such network design means that, in the absence of their regulators, the genes are found in their least required activity state, hence regulatory intervention is often necessary. Due to the multitude of genes and regulators, spurious binding and unbinding events, called “crosstalk”, could occur. To determine how the form of regulation affects the global crosstalk in the network, we used a mathematical model that includes multiple regulators and multiple target genes. We found that crosstalk depends non-monotonically on the availability of regulators. Our analysis showed that excess use of regulation entailed by the formerly suggested network design caused high crosstalk levels in a large part of the parameter space. We therefore considered the opposite ‘idle’ design, where the default unregulated state of genes is their frequently required activity state. We found, that ‘idle’ design minimized the use of regulation and thus minimized crosstalk. In addition, we estimated global crosstalk of S. cerevisiae using transcription factors binding data. We demonstrated that even partial network data could suffice to estimate its global crosstalk, suggesting its applicability to additional organisms. We found that S. cerevisiae estimated crosstalk is lower than that of a random network, suggesting that natural selection reduces crosstalk. In summary, our study highlights a new type of protein production cost which is typically overlooked: that of regulatory interference caused by the presence of excess regulators in the cell. It demonstrates the importance of whole-network descriptions, which could show effects missed by single-gene models.}, author = {Grah, Rok and Friedlander, Tamar}, issn = {1553-7358}, journal = {PLOS Computational Biology}, number = {2}, publisher = {Public Library of Science}, title = {{The relation between crosstalk and gene regulation form revisited}}, doi = {10.1371/journal.pcbi.1007642}, volume = {16}, year = {2020}, } @phdthesis{8155, abstract = {In the thesis we focus on the interplay of the biophysics and evolution of gene regulation. We start by addressing how the type of prokaryotic gene regulation – activation and repression – affects spurious binding to DNA, also known as transcriptional crosstalk. We propose that regulatory interference caused by excess regulatory proteins in the dense cellular medium – global crosstalk – could be a factor in determining which type of gene regulatory network is evolutionarily preferred. Next,we use a normative approach in eukaryotic gene regulation to describe minimal non-equilibrium enhancer models that optimize so-called regulatory phenotypes. We find a class of models that differ from standard thermodynamic equilibrium models by a single parameter that notably increases the regulatory performance. Next chapter addresses the question of genotype-phenotype-fitness maps of higher dimensional phenotypes. We show that our biophysically realistic approach allows us to understand how the mechanisms of promoter function constrain genotypephenotype maps, and how they affect the evolutionary trajectories of promoters. In the last chapter we ask whether the intrinsic instability of gene duplication and amplification provides a generic alternative to canonical gene regulation. Using mathematical modeling, we show that amplifications can tune gene expression in many environments, including those where transcription factor-based schemes are hard to evolve or maintain. }, author = {Grah, Rok}, issn = {2663-337X}, pages = {310}, publisher = {Institute of Science and Technology Austria}, title = {{Gene regulation across scales – how biophysical constraints shape evolution}}, doi = {10.15479/AT:ISTA:8155}, year = {2020}, } @unpublished{7675, abstract = {In prokaryotes, thermodynamic models of gene regulation provide a highly quantitative mapping from promoter sequences to gene expression levels that is compatible with in vivo and in vitro bio-physical measurements. Such concordance has not been achieved for models of enhancer function in eukaryotes. In equilibrium models, it is difficult to reconcile the reported short transcription factor (TF) residence times on the DNA with the high specificity of regulation. In non-equilibrium models, progress is difficult due to an explosion in the number of parameters. Here, we navigate this complexity by looking for minimal non-equilibrium enhancer models that yield desired regulatory phenotypes: low TF residence time, high specificity and tunable cooperativity. We find that a single extra parameter, interpretable as the “linking rate” by which bound TFs interact with Mediator components, enables our models to escape equilibrium bounds and access optimal regulatory phenotypes, while remaining consistent with the reported phenomenology and simple enough to be inferred from upcoming experiments. We further find that high specificity in non-equilibrium models is in a tradeoff with gene expression noise, predicting bursty dynamics — an experimentally-observed hallmark of eukaryotic transcription. By drastically reducing the vast parameter space to a much smaller subspace that optimally realizes biological function prior to inference from data, our normative approach holds promise for mathematical models in systems biology.}, author = {Grah, Rok and Zoller, Benjamin and Tkačik, Gašper}, booktitle = {bioRxiv}, publisher = {Cold Spring Harbor Laboratory}, title = {{Normative models of enhancer function}}, doi = {10.1101/2020.04.08.029405}, year = {2020}, } @article{7656, abstract = {We propose that correlations among neurons are generically strong enough to organize neural activity patterns into a discrete set of clusters, which can each be viewed as a population codeword. Our reasoning starts with the analysis of retinal ganglion cell data using maximum entropy models, showing that the population is robustly in a frustrated, marginally sub-critical, or glassy, state. This leads to an argument that neural populations in many other brain areas might share this structure. Next, we use latent variable models to show that this glassy state possesses well-defined clusters of neural activity. Clusters have three appealing properties: (i) clusters exhibit error correction, i.e., they are reproducibly elicited by the same stimulus despite variability at the level of constituent neurons; (ii) clusters encode qualitatively different visual features than their constituent neurons; and (iii) clusters can be learned by downstream neural circuits in an unsupervised fashion. We hypothesize that these properties give rise to a “learnable” neural code which the cortical hierarchy uses to extract increasingly complex features without supervision or reinforcement.}, author = {Berry, Michael J. and Tkačik, Gašper}, issn = {1662-5188}, journal = {Frontiers in Computational Neuroscience}, publisher = {Frontiers}, title = {{Clustering of neural activity: A design principle for population codes}}, doi = {10.3389/fncom.2020.00020}, volume = {14}, year = {2020}, } @phdthesis{8657, abstract = {Synthesis of proteins – translation – is a fundamental process of life. Quantitative studies anchor translation into the context of bacterial physiology and reveal several mathematical relationships, called “growth laws,” which capture physiological feedbacks between protein synthesis and cell growth. Growth laws describe the dependency of the ribosome abundance as a function of growth rate, which can change depending on the growth conditions. Perturbations of translation reveal that bacteria employ a compensatory strategy in which the reduced translation capability results in increased expression of the translation machinery. Perturbations of translation are achieved in various ways; clinically interesting is the application of translation-targeting antibiotics – translation inhibitors. The antibiotic effects on bacterial physiology are often poorly understood. Bacterial responses to two or more simultaneously applied antibiotics are even more puzzling. The combined antibiotic effect determines the type of drug interaction, which ranges from synergy (the effect is stronger than expected) to antagonism (the effect is weaker) and suppression (one of the drugs loses its potency). In the first part of this work, we systematically measure the pairwise interaction network for translation inhibitors that interfere with different steps in translation. We find that the interactions are surprisingly diverse and tend to be more antagonistic. To explore the underlying mechanisms, we begin with a minimal biophysical model of combined antibiotic action. We base this model on the kinetics of antibiotic uptake and binding together with the physiological response described by the growth laws. The biophysical model explains some drug interactions, but not all; it specifically fails to predict suppression. In the second part of this work, we hypothesize that elusive suppressive drug interactions result from the interplay between ribosomes halted in different stages of translation. To elucidate this putative mechanism of drug interactions between translation inhibitors, we generate translation bottlenecks genetically using in- ducible control of translation factors that regulate well-defined translation cycle steps. These perturbations accurately mimic antibiotic action and drug interactions, supporting that the interplay of different translation bottlenecks partially causes these interactions. We extend this approach by varying two translation bottlenecks simultaneously. This approach reveals the suppression of translocation inhibition by inhibited translation. We rationalize this effect by modeling dense traffic of ribosomes that move on transcripts in a translation factor-mediated manner. This model predicts a dissolution of traffic jams caused by inhibited translocation when the density of ribosome traffic is reduced by lowered initiation. We base this model on the growth laws and quantitative relationships between different translation and growth parameters. In the final part of this work, we describe a set of tools aimed at quantification of physiological and translation parameters. We further develop a simple model that directly connects the abundance of a translation factor with the growth rate, which allows us to extract physiological parameters describing initiation. We demonstrate the development of tools for measuring translation rate. This thesis showcases how a combination of high-throughput growth rate mea- surements, genetics, and modeling can reveal mechanisms of drug interactions. Furthermore, by a gradual transition from combinations of antibiotics to precise genetic interventions, we demonstrated the equivalency between genetic and chemi- cal perturbations of translation. These findings tile the path for quantitative studies of antibiotic combinations and illustrate future approaches towards the quantitative description of translation.}, author = {Kavcic, Bor}, isbn = {978-3-99078-011-4}, issn = {2663-337X}, pages = {271}, publisher = {Institute of Science and Technology Austria}, title = {{Perturbations of protein synthesis: from antibiotics to genetics and physiology}}, doi = {10.15479/AT:ISTA:8657}, year = {2020}, } @article{8250, abstract = {Antibiotics that interfere with translation, when combined, interact in diverse and difficult-to-predict ways. Here, we explain these interactions by “translation bottlenecks”: points in the translation cycle where antibiotics block ribosomal progression. To elucidate the underlying mechanisms of drug interactions between translation inhibitors, we generate translation bottlenecks genetically using inducible control of translation factors that regulate well-defined translation cycle steps. These perturbations accurately mimic antibiotic action and drug interactions, supporting that the interplay of different translation bottlenecks causes these interactions. We further show that growth laws, combined with drug uptake and binding kinetics, enable the direct prediction of a large fraction of observed interactions, yet fail to predict suppression. However, varying two translation bottlenecks simultaneously supports that dense traffic of ribosomes and competition for translation factors account for the previously unexplained suppression. These results highlight the importance of “continuous epistasis” in bacterial physiology.}, author = {Kavcic, Bor and Tkačik, Gašper and Bollenbach, Tobias}, issn = {2041-1723}, journal = {Nature Communications}, publisher = {Springer Nature}, title = {{Mechanisms of drug interactions between translation-inhibiting antibiotics}}, doi = {10.1038/s41467-020-17734-z}, volume = {11}, year = {2020}, } @unpublished{7673, abstract = {Combining drugs can improve the efficacy of treatments. However, predicting the effect of drug combinations is still challenging. The combined potency of drugs determines the drug interaction, which is classified as synergistic, additive, antagonistic, or suppressive. While probabilistic, non-mechanistic models exist, there is currently no biophysical model that can predict antibiotic interactions. Here, we present a physiologically relevant model of the combined action of antibiotics that inhibit protein synthesis by targeting the ribosome. This model captures the kinetics of antibiotic binding and transport, and uses bacterial growth laws to predict growth in the presence of antibiotic combinations. We find that this biophysical model can produce all drug interaction types except suppression. We show analytically that antibiotics which cannot bind to the ribosome simultaneously generally act as substitutes for one another, leading to additive drug interactions. Previously proposed null expectations for higher-order drug interactions follow as a limiting case of our model. We further extend the model to include the effects of direct physical or allosteric interactions between individual drugs on the ribosome. Notably, such direct interactions profoundly change the combined drug effect, depending on the kinetic parameters of the drugs used. The model makes additional predictions for the effects of resistance genes on drug interactions and for interactions between ribosome-targeting antibiotics and antibiotics with other targets. These findings enhance our understanding of the interplay between drug action and cell physiology and are a key step toward a general framework for predicting drug interactions.}, author = {Kavcic, Bor and Tkačik, Gašper and Bollenbach, Tobias}, booktitle = {bioRxiv}, publisher = {Cold Spring Harbor Laboratory}, title = {{A minimal biophysical model of combined antibiotic action}}, doi = {10.1101/2020.04.18.047886}, year = {2020}, } @article{7652, abstract = {Organisms cope with change by taking advantage of transcriptional regulators. However, when faced with rare environments, the evolution of transcriptional regulators and their promoters may be too slow. Here, we investigate whether the intrinsic instability of gene duplication and amplification provides a generic alternative to canonical gene regulation. Using real-time monitoring of gene-copy-number mutations in Escherichia coli, we show that gene duplications and amplifications enable adaptation to fluctuating environments by rapidly generating copy-number and, therefore, expression-level polymorphisms. This amplification-mediated gene expression tuning (AMGET) occurs on timescales that are similar to canonical gene regulation and can respond to rapid environmental changes. Mathematical modelling shows that amplifications also tune gene expression in stochastic environments in which transcription-factor-based schemes are hard to evolve or maintain. The fleeting nature of gene amplifications gives rise to a generic population-level mechanism that relies on genetic heterogeneity to rapidly tune the expression of any gene, without leaving any genomic signature.}, author = {Tomanek, Isabella and Grah, Rok and Lagator, M. and Andersson, A. M. C. and Bollback, Jonathan P and Tkačik, Gašper and Guet, Calin C}, issn = {2397-334X}, journal = {Nature Ecology & Evolution}, number = {4}, pages = {612--625}, publisher = {Springer Nature}, title = {{Gene amplification as a form of population-level gene expression regulation}}, doi = {10.1038/s41559-020-1132-7}, volume = {4}, year = {2020}, } @inbook{19988, abstract = {Quantitative studies of cell metabolism are often based on large chemical reaction network models. A steady-state approach is suited to analyze phenomena on the timescale of cell growth and circumvents the problem of incomplete experimental knowledge on kinetic laws and parameters, but it should be supported by a correct implementation of thermodynamic constraints. In this chapter, we review the latter aspect, highlighting its computational challenges and physical insights. The simple introduction of Gibbs inequalities avoids the presence of unfeasible loops allowing for correct timescale analysis, but leads to possibly non-convex feasible flux spaces whose exploration needs efficient algorithms. We briefly review the implementation of thermodynamics through variational principles in constraint-based models of metabolic networks.}, author = {De Martino, A and De Martino, Daniele and Marinari, E}, booktitle = {Chemical Kinetics}, isbn = {9781786347008}, pages = {455--471}, publisher = {World Scientific Publishing}, title = {{The Essential Role of Thermodynamics in Metabolic Network Modeling: Physical Insights and Computational Challenges}}, doi = {10.1142/9781786347015_0018}, year = {2019}, } @article{196, abstract = {The abelian sandpile serves as a model to study self-organized criticality, a phenomenon occurring in biological, physical and social processes. The identity of the abelian group is a fractal composed of self-similar patches, and its limit is subject of extensive collaborative research. Here, we analyze the evolution of the sandpile identity under harmonic fields of different orders. We show that this evolution corresponds to periodic cycles through the abelian group characterized by the smooth transformation and apparent conservation of the patches constituting the identity. The dynamics induced by second and third order harmonics resemble smooth stretchings, respectively translations, of the identity, while the ones induced by fourth order harmonics resemble magnifications and rotations. Starting with order three, the dynamics pass through extended regions of seemingly random configurations which spontaneously reassemble into accentuated patterns. We show that the space of harmonic functions projects to the extended analogue of the sandpile group, thus providing a set of universal coordinates identifying configurations between different domains. Since the original sandpile group is a subgroup of the extended one, this directly implies that it admits a natural renormalization. Furthermore, we show that the harmonic fields can be induced by simple Markov processes, and that the corresponding stochastic dynamics show remarkable robustness over hundreds of periods. Finally, we encode information into seemingly random configurations, and decode this information with an algorithm requiring minimal prior knowledge. Our results suggest that harmonic fields might split the sandpile group into sub-sets showing different critical coefficients, and that it might be possible to extend the fractal structure of the identity beyond the boundaries of its domain. }, author = {Lang, Moritz and Shkolnikov, Mikhail}, issn = {1091-6490}, journal = {Proceedings of the National Academy of Sciences of the United States of America}, number = {8}, pages = {2821--2830}, publisher = {National Academy of Sciences}, title = {{Harmonic dynamics of the Abelian sandpile}}, doi = {10.1073/pnas.1812015116}, volume = {116}, year = {2019}, } @misc{9786, author = {Ruess, Jakob and Pleska, Maros and Guet, Calin C and Tkačik, Gašper}, publisher = {Public Library of Science}, title = {{Supporting text and results}}, doi = {10.1371/journal.pcbi.1007168.s001}, year = {2019}, } @article{5817, abstract = {We theoretically study the shapes of lipid vesicles confined to a spherical cavity, elaborating a framework based on the so-called limiting shapes constructed from geometrically simple structural elements such as double-membrane walls and edges. Partly inspired by numerical results, the proposed non-compartmentalized and compartmentalized limiting shapes are arranged in the bilayer-couple phase diagram which is then compared to its free-vesicle counterpart. We also compute the area-difference-elasticity phase diagram of the limiting shapes and we use it to interpret shape transitions experimentally observed in vesicles confined within another vesicle. The limiting-shape framework may be generalized to theoretically investigate the structure of certain cell organelles such as the mitochondrion.}, author = {Kavcic, Bor and Sakashita, A. and Noguchi, H. and Ziherl, P.}, issn = {1744-6848}, journal = {Soft Matter}, number = {4}, pages = {602--614}, publisher = {Royal Society of Chemistry}, title = {{Limiting shapes of confined lipid vesicles}}, doi = {10.1039/c8sm01956h}, volume = {15}, year = {2019}, } @article{5945, abstract = {In developing organisms, spatially prescribed cell identities are thought to be determined by the expression levels of multiple genes. Quantitative tests of this idea, however, require a theoretical framework capable of exposing the rules and precision of cell specification over developmental time. We use the gap gene network in the early fly embryo as an example to show how expression levels of the four gap genes can be jointly decoded into an optimal specification of position with 1% accuracy. The decoder correctly predicts, with no free parameters, the dynamics of pair-rule expression patterns at different developmental time points and in various mutant backgrounds. Precise cellular identities are thus available at the earliest stages of development, contrasting the prevailing view of positional information being slowly refined across successive layers of the patterning network. Our results suggest that developmental enhancers closely approximate a mathematically optimal decoding strategy.}, author = {Petkova, Mariela D. and Tkacik, Gasper and Bialek, William and Wieschaus, Eric F. and Gregor, Thomas}, journal = {Cell}, number = {4}, pages = {844--855.e15}, publisher = {Cell Press}, title = {{Optimal decoding of cellular identities in a genetic network}}, doi = {10.1016/j.cell.2019.01.007}, volume = {176}, year = {2019}, } @article{6046, abstract = {Sudden stress often triggers diverse, temporally structured gene expression responses in microbes, but it is largely unknown how variable in time such responses are and if genes respond in the same temporal order in every single cell. Here, we quantified timing variability of individual promoters responding to sublethal antibiotic stress using fluorescent reporters, microfluidics, and time‐lapse microscopy. We identified lower and upper bounds that put definite constraints on timing variability, which varies strongly among promoters and conditions. Timing variability can be interpreted using results from statistical kinetics, which enable us to estimate the number of rate‐limiting molecular steps underlying different responses. We found that just a few critical steps control some responses while others rely on dozens of steps. To probe connections between different stress responses, we then tracked the temporal order and response time correlations of promoter pairs in individual cells. Our results support that, when bacteria are exposed to the antibiotic nitrofurantoin, the ensuing oxidative stress and SOS responses are part of the same causal chain of molecular events. In contrast, under trimethoprim, the acid stress response and the SOS response are part of different chains of events running in parallel. Our approach reveals fundamental constraints on gene expression timing and provides new insights into the molecular events that underlie the timing of stress responses.}, author = {Mitosch, Karin and Rieckh, Georg and Bollenbach, Mark Tobias}, journal = {Molecular systems biology}, number = {2}, publisher = {Embo Press}, title = {{Temporal order and precision of complex stress responses in individual bacteria}}, doi = {10.15252/msb.20188470}, volume = {15}, year = {2019}, } @article{6049, abstract = {In this article it is shown that large systems with many interacting units endowing multiple phases display self-oscillations in the presence of linear feedback between the control and order parameters, where an Andronov–Hopf bifurcation takes over the phase transition. This is simply illustrated through the mean field Landau theory whose feedback dynamics turn out to be described by the Van der Pol equation and it is then validated for the fully connected Ising model following heat bath dynamics. Despite its simplicity, this theory accounts potentially for a rich range of phenomena: here it is applied to describe in a stylized way (i) excess demand-price cycles due to strong herding in a simple agent-based market model; (ii) congestion waves in queuing networks triggered by user feedback to delays in overloaded conditions; and (iii) metabolic network oscillations resulting from cell growth control in a bistable phenotypic landscape.}, author = {De Martino, Daniele}, journal = {Journal of Physics A: Mathematical and Theoretical}, number = {4}, publisher = {IOP Publishing}, title = {{Feedback-induced self-oscillations in large interacting systems subjected to phase transitions}}, doi = {10.1088/1751-8121/aaf2dd}, volume = {52}, year = {2019}, } @article{6090, abstract = {Cells need to reliably sense external ligand concentrations to achieve various biological functions such as chemotaxis or signaling. The molecular recognition of ligands by surface receptors is degenerate in many systems, leading to crosstalk between ligand-receptor pairs. Crosstalk is often thought of as a deviation from optimal specific recognition, as the binding of noncognate ligands can interfere with the detection of the receptor's cognate ligand, possibly leading to a false triggering of a downstream signaling pathway. Here we quantify the optimal precision of sensing the concentrations of multiple ligands by a collection of promiscuous receptors. We demonstrate that crosstalk can improve precision in concentration sensing and discrimination tasks. To achieve superior precision, the additional information about ligand concentrations contained in short binding events of the noncognate ligand should be exploited. We present a proofreading scheme to realize an approximate estimation of multiple ligand concentrations that reaches a precision close to the derived optimal bounds. Our results help rationalize the observed ubiquity of receptor crosstalk in molecular sensing.}, author = {Carballo-Pacheco, Martín and Desponds, Jonathan and Gavrilchenko, Tatyana and Mayer, Andreas and Prizak, Roshan and Reddy, Gautam and Nemenman, Ilya and Mora, Thierry}, journal = {Physical Review E}, number = {2}, publisher = {American Physical Society}, title = {{Receptor crosstalk improves concentration sensing of multiple ligands}}, doi = {10.1103/PhysRevE.99.022423}, volume = {99}, year = {2019}, } @article{6784, abstract = {Mathematical models have been used successfully at diverse scales of biological organization, ranging from ecology and population dynamics to stochastic reaction events occurring between individual molecules in single cells. Generally, many biological processes unfold across multiple scales, with mutations being the best studied example of how stochasticity at the molecular scale can influence outcomes at the population scale. In many other contexts, however, an analogous link between micro- and macro-scale remains elusive, primarily due to the challenges involved in setting up and analyzing multi-scale models. Here, we employ such a model to investigate how stochasticity propagates from individual biochemical reaction events in the bacterial innate immune system to the ecology of bacteria and bacterial viruses. We show analytically how the dynamics of bacterial populations are shaped by the activities of immunity-conferring enzymes in single cells and how the ecological consequences imply optimal bacterial defense strategies against viruses. Our results suggest that bacterial populations in the presence of viruses can either optimize their initial growth rate or their population size, with the first strategy favoring simple immunity featuring a single restriction modification system and the second strategy favoring complex bacterial innate immunity featuring several simultaneously active restriction modification systems.}, author = {Ruess, Jakob and Pleska, Maros and Guet, Calin C and Tkačik, Gašper}, issn = {1553-7358}, journal = {PLoS Computational Biology}, number = {7}, publisher = {Public Library of Science}, title = {{Molecular noise of innate immunity shapes bacteria-phage ecologies}}, doi = {10.1371/journal.pcbi.1007168}, volume = {15}, year = {2019}, } @article{7103, abstract = {Origin and functions of intermittent transitions among sleep stages, including short awakenings and arousals, constitute a challenge to the current homeostatic framework for sleep regulation, focusing on factors modulating sleep over large time scales. Here we propose that the complex micro-architecture characterizing the sleep-wake cycle results from an underlying non-equilibrium critical dynamics, bridging collective behaviors across spatio-temporal scales. We investigate θ and δ wave dynamics in control rats and in rats with lesions of sleep-promoting neurons in the parafacial zone. We demonstrate that intermittent bursts in θ and δ rhythms exhibit a complex temporal organization, with long-range power-law correlations and a robust duality of power law (θ-bursts, active phase) and exponential-like (δ-bursts, quiescent phase) duration distributions, typical features of non-equilibrium systems self-organizing at criticality. Crucially, such temporal organization relates to anti-correlated coupling between θ- and δ-bursts, and is independent of the dominant physiologic state and lesions, a solid indication of a basic principle in sleep dynamics.}, author = {Wang, Jilin W. J. L. and Lombardi, Fabrizio and Zhang, Xiyun and Anaclet, Christelle and Ivanov, Plamen Ch.}, issn = {1553-7358}, journal = {PLoS Computational Biology}, number = {11}, publisher = {Public Library of Science}, title = {{Non-equilibrium critical dynamics of bursts in θ and δ rhythms as fundamental characteristic of sleep and wake micro-architecture}}, doi = {10.1371/journal.pcbi.1007268}, volume = {15}, year = {2019}, } @article{7422, abstract = {Biochemical reactions often occur at low copy numbers but at once in crowded and diverse environments. Space and stochasticity therefore play an essential role in biochemical networks. Spatial-stochastic simulations have become a prominent tool for understanding how stochasticity at the microscopic level influences the macroscopic behavior of such systems. While particle-based models guarantee the level of detail necessary to accurately describe the microscopic dynamics at very low copy numbers, the algorithms used to simulate them typically imply trade-offs between computational efficiency and biochemical accuracy. eGFRD (enhanced Green’s Function Reaction Dynamics) is an exact algorithm that evades such trade-offs by partitioning the N-particle system into M ≤ N analytically tractable one- and two-particle systems; the analytical solutions (Green’s functions) then are used to implement an event-driven particle-based scheme that allows particles to make large jumps in time and space while retaining access to their state variables at arbitrary simulation times. Here we present “eGFRD2,” a new eGFRD version that implements the principle of eGFRD in all dimensions, thus enabling efficient particle-based simulation of biochemical reaction-diffusion processes in the 3D cytoplasm, on 2D planes representing membranes, and on 1D elongated cylinders representative of, e.g., cytoskeletal tracks or DNA; in 1D, it also incorporates convective motion used to model active transport. We find that, for low particle densities, eGFRD2 is up to 6 orders of magnitude faster than conventional Brownian dynamics. We exemplify the capabilities of eGFRD2 by simulating an idealized model of Pom1 gradient formation, which involves 3D diffusion, active transport on microtubules, and autophosphorylation on the membrane, confirming recent experimental and theoretical results on this system to hold under genuinely stochastic conditions.}, author = {Sokolowski, Thomas R and Paijmans, Joris and Bossen, Laurens and Miedema, Thomas and Wehrens, Martijn and Becker, Nils B. and Kaizu, Kazunari and Takahashi, Koichi and Dogterom, Marileen and ten Wolde, Pieter Rein}, issn = {1089-7690}, journal = {The Journal of Chemical Physics}, number = {5}, publisher = {AIP Publishing}, title = {{eGFRD in all dimensions}}, doi = {10.1063/1.5064867}, volume = {150}, year = {2019}, }