Showing 1–40 of 40 results for author: Müller, J

Analyzing the Brain's Dynamic Response to Targeted Stimulation using Generative Modeling

Authors: Rishikesan Maran, Eli J. Müller, Ben D. Fulcher

Abstract: Generative models of brain activity have been instrumental in testing hypothesized mechanisms underlying brain dynamics against experimental datasets. Beyond capturing the key mechanisms underlying spontaneous brain dynamics, these models hold an exciting potential for understanding the mechanisms underlying the dynamics evoked by targeted brain-stimulation techniques. This paper delves into this… ▽ More Generative models of brain activity have been instrumental in testing hypothesized mechanisms underlying brain dynamics against experimental datasets. Beyond capturing the key mechanisms underlying spontaneous brain dynamics, these models hold an exciting potential for understanding the mechanisms underlying the dynamics evoked by targeted brain-stimulation techniques. This paper delves into this emerging application, using concepts from dynamical systems theory to argue that the stimulus-evoked dynamics in such experiments may be shaped by new types of mechanisms distinct from those that dominate spontaneous dynamics. We review and discuss: (i) the targeted experimental techniques across spatial scales that can both perturb the brain to novel states and resolve its relaxation trajectory back to spontaneous dynamics; and (ii) how we can understand these dynamics in terms of mechanisms using physiological, phenomenological, and data-driven models. A tight integration of targeted stimulation experiments with generative quantitative modeling provides an important opportunity to uncover novel mechanisms of brain dynamics that are difficult to detect in spontaneous settings. △ Less

Submitted 1 August, 2024; v1 submitted 29 July, 2024; originally announced July 2024.

Unified neural field theory of brain dynamics underlying oscillations in Parkinson's disease and generalized epilepsies

Authors: Eli J Müller, Sacha J van Albada, Jong-Won Kim, Peter A Robinson

Abstract: The mechanisms underlying pathologically synchronized neural oscillations in Parkinson's disease (PD) and generalized epilepsies are jointly explored via a neural field model of the corticothalamic-basal ganglia (CTBG) system. The basal ganglia (BG) are approximated as a single effective population and their roles in modulating oscillatory corticothalamic (CT) dynamics and vice versa are analyzed.… ▽ More The mechanisms underlying pathologically synchronized neural oscillations in Parkinson's disease (PD) and generalized epilepsies are jointly explored via a neural field model of the corticothalamic-basal ganglia (CTBG) system. The basal ganglia (BG) are approximated as a single effective population and their roles in modulating oscillatory corticothalamic (CT) dynamics and vice versa are analyzed. Besides normal EEG rhythms, enhanced activity around 4 Hz and 20 Hz exists in the model, consistent with characteristic frequencies in PD. These rhythms result from resonances in loops between the BG and CT populations, analogous to those underlying epileptic oscillations in a previous CT model. Dopamine depletion is argued to weaken the dampening of these resonances in PD, and network connections explain the significant coherence between BG, thalamic, and cortical activity around 4-8 Hz and 20 Hz. Parallels between the afferent and efferent connection sites of the thalamic reticular nucleus (TRN) and BG predict low dopamine to correspond to a reduced likelihood of tonic-clonic (grand mal) seizures, agreeing with experimental findings. Further, the model predicts an increased likelihood of absence (petit mal) seizure resulting from low dopamine levels matching experimental findings. Suppression of absence seizure activity is shown when afferent and efferent BG connections to the CT system are strengthened, consistent with other CTBG modeling studies. The BG are demonstrated to suppress activity of the CTBG system near tonic-clonic seizure states, providing insight into the reported efficacy of current treatments in BG circuits. Sleep states of the TRN are also found to suppress pathological PD activity matching observations. Overall, the findings demonstrate strong parallels between coherent oscillations in generalized epilepsies and PD, and provide insights into possible comorbidities. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Journal ref: Journal of Theoretical Biology (2017) 428, 132-146

Parameter estimation for contact tracing in graph-based models

Authors: Augustine Okolie, Johannes Müller, Mirjam Kretzschmar

Abstract: We adopt a maximum-likelihood framework to estimate parameters of a stochastic susceptible-infected-recovered (SIR) model with contact tracing on a rooted random tree. Given the number of detectees per index case, our estimator allows to determine the degree distribution of the random tree as well as the tracing probability. Since we do not discover all infectees via contact tracing, this estimati… ▽ More We adopt a maximum-likelihood framework to estimate parameters of a stochastic susceptible-infected-recovered (SIR) model with contact tracing on a rooted random tree. Given the number of detectees per index case, our estimator allows to determine the degree distribution of the random tree as well as the tracing probability. Since we do not discover all infectees via contact tracing, this estimation is non-trivial. To keep things simple and stable, we develop an approximation suited for realistic situations (contract tracing probability small, or the probability for the detection of index cases small). In this approximation, the only epidemiological parameter entering the estimator is $R_0$. The estimator is tested in a simulation study and is furthermore applied to covid-19 contact tracing data from India. The simulation study underlines the efficiency of the method. For the empirical covid-19 data, we compare different degree distributions and perform a sensitivity analysis. We find that particularly a power-law and a negative binomial degree distribution fit the data well and that the tracing probability is rather large. The sensitivity analysis shows no strong dependency of the estimates on the reproduction number. Finally, we discuss the relevance of our findings. △ Less

Submitted 22 November, 2023; v1 submitted 17 July, 2023; originally announced July 2023.

Comments: 24 pages, 8 figures, 3 tables

MSC Class: 92D30

Journal ref: Royal Society Interface 2023

Age structure, replicator equation, and the prisoner's dilemma

Authors: Sona John, Johannes Müller

Abstract: We investigate the evolutionary dynamics of an age-structured population under weak frequency-dependent selection. It turns out that the weak selection is affected in a non-trivial way by the life-history trait. We can disentangle the dynamics, based on the appearance of different time scales. These time scales, which seem to form a universal structure in the interplay of weak selection and life-h… ▽ More We investigate the evolutionary dynamics of an age-structured population under weak frequency-dependent selection. It turns out that the weak selection is affected in a non-trivial way by the life-history trait. We can disentangle the dynamics, based on the appearance of different time scales. These time scales, which seem to form a universal structure in the interplay of weak selection and life-history traits, allow us to reduce the infinite dimensional model to a one-dimensional modified replicator equation. The modified replicator equation is then used to investigate cooperation (the prisoner's dilemma) by means of adaptive dynamics. We identify conditions under which age structure is able to promote cooperation. At the end we discuss the relevance of our findings. △ Less

Submitted 19 December, 2022; originally announced December 2022.

MSC Class: 91A22; 92D25

Quiescence generates moving average in a stochastic epidemiological model with one host and two parasites

Authors: Usman Sanusi, Sona John, Johannes Mueller, Aurélien Tellier

Abstract: Mathematical modelling of epidemiological and coevolutionary dynamics is widely being used to improve disease management strategies of infectious diseases. Many diseases present some form of intra-host quiescent stage, also known as covert infection, while others exhibit dormant stages in the environment. As quiescent/dormant stages can be resistant to drug, antibiotics, fungicide treatments, it i… ▽ More Mathematical modelling of epidemiological and coevolutionary dynamics is widely being used to improve disease management strategies of infectious diseases. Many diseases present some form of intra-host quiescent stage, also known as covert infection, while others exhibit dormant stages in the environment. As quiescent/dormant stages can be resistant to drug, antibiotics, fungicide treatments, it is of practical relevance to study the influence of these two life-history traits on the coevolutionary dynamics. We develop first a deterministic coevolutionary model with two parasite types infecting one host type and study analytically the stability of the dynamical system. We specifically derive a stability condition for a five-by-five system of equations with quiescence. Second, we develop a stochastic version of the model to study the influence of quiescence on stochasticity of the system dynamics. We compute the steady state distribution of the parasite types which follows a multivariate normal distribution. Furthermore, we obtain numerical solutions for the covariance matrix of the system under symmetric and asymmetric quiescence rates between parasite types. When parasite strains are identical, quiescence increases the variance of the number of infected individuals at high transmission rate and vice versa when the transmission rate is low. However, when there is competition between parasite strains with different quiescent rates, quiescence generates a moving average behaviour which dampen off stochasticity and decreases the variance of the number of infected hosts. The strain with the highest rate of entering quiescence determines the strength of the moving average and the magnitude of reduction of stochasticity. Thus, it is worth investigating simple models of multi-strain parasite under quiescence/dormancy to improve disease management strategies. △ Less

Submitted 25 May, 2022; originally announced May 2022.

Seizure prediction with long-term iEEG recordings: What can we learn from data nonstationarity?

Authors: Hongliu Yang, Matthias Eberlein, Jens Müller, Ronald Tetzlaff

Abstract: Repeated epileptic seizures impair around 65 million people worldwide and a successful prediction of seizures could significantly help patients suffering from refractory epilepsy. For two dogs with yearlong intracranial electroencephalography (iEEG) recordings, we studied the influence of time series nonstationarity on the performance of seizure prediction using in-house developed machine learning… ▽ More Repeated epileptic seizures impair around 65 million people worldwide and a successful prediction of seizures could significantly help patients suffering from refractory epilepsy. For two dogs with yearlong intracranial electroencephalography (iEEG) recordings, we studied the influence of time series nonstationarity on the performance of seizure prediction using in-house developed machine learning algorithms. We observed a long-term evolution on the scale of weeks or months in iEEG time series that may be represented as switching between certain meta-states. To better predict impending seizures, retraining of prediction algorithms is therefore necessary and the retraining schedule should be adjusted to the change in meta-states. There is evidence that the nature of seizure-free interictal clips also changes with the transition between meta-states, accwhich has been shown relevant for seizure prediction. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: accepted for MLESP 2021

Journal ref: 2021 IEEE International Conference on Bioinformatics and Biomedicine (BIBM)

A model of opinion dynamics with echo chambers explains the spatial distribution of vaccine hesitancy

Authors: Johannes Müller, Aurelien Tellier, Michael Kurschilgen

Abstract: Vaccination hesitancy is a major obstacle to achieving and maintaining herd immunity. It is therefore of prime importance for public health authorities to understand the dynamics of an anti-vaccine opinion in the population. We introduce a novel mathematical model of opinion dynamics with spatial reinforcement, which can generate echo chambers, i.e. opinion bubbles in which information that is inc… ▽ More Vaccination hesitancy is a major obstacle to achieving and maintaining herd immunity. It is therefore of prime importance for public health authorities to understand the dynamics of an anti-vaccine opinion in the population. We introduce a novel mathematical model of opinion dynamics with spatial reinforcement, which can generate echo chambers, i.e. opinion bubbles in which information that is incompatible with one's entrenched worldview, is likely disregarded. In a first mathematical part, we scale the model both to a deterministic limit and to a weak-effects limit, and obtain bifurcations, phase transitions, and the invariant measure. In a second part, we fit our model to measles and meningococci vaccination coverage across 413 districts in Germany. We reveal that strong echo chambers explain the occurrence and persistence of the anti-vaccination opinion. We predict and compare the effectiveness of different policies aimed at influencing opinion dynamics in order to increase vaccination uptake. △ Less

Submitted 19 December, 2021; originally announced December 2021.

MSC Class: 60J20 (primary) 92D25; 91D30 (secondary)

Life-History traits and the replicator equation

Authors: Johannes Müller, Aurélien Tellier

Abstract: Due to the relevance for conservation biology, there is an increasing interest to extend evolutionary genomics models to plant, animal or microbial species. However, this requires to understand the effect of life-history traits absent in humans on genomic evolution. In this context, it is fundamentally of interest to generalize the replicator equation, which is at the heart of most population geno… ▽ More Due to the relevance for conservation biology, there is an increasing interest to extend evolutionary genomics models to plant, animal or microbial species. However, this requires to understand the effect of life-history traits absent in humans on genomic evolution. In this context, it is fundamentally of interest to generalize the replicator equation, which is at the heart of most population genomics models. However, as the inclusion of life-history traits generates models with a large state space, the analysis becomes involving. We focus, here, on quiescence and seed banks, two features common to many plant, invertebrate and microbial species. We develop a method to obtain a low-dimensional replicator equation in the context of evolutionary game theory, based on two assumptions: (1) the life-history traits are {\it per se} neutral, and (2) frequency-dependent selection is weak. We use the results to investigate the evolution and maintenance of cooperation based on the Prisoner's dilemma. We first consider the generalized replicator equation, and then refine the investigation using adaptive dynamics. It turns out that, depending on the structure and timing of the quiescence/dormancy life-history trait, cooperation in a homogeneous population can be stabilized. We finally discuss and highlight the relevance of these results for plant, invertebrate and microbial communities. △ Less

Submitted 13 November, 2021; originally announced November 2021.

Noradrenergic neuromodulation of nonlinear bursting neurons controls critical dynamics

Authors: Brandon R. Munn, Eli J. Müller, James M. Shine

Abstract: In order to remain adaptable to a dynamic environment, neural activity must be simultaneously both sensitive and stable. To solve this problem, the brain has been hypothesised to sit near a critical boundary. Yet, precisely how criticality and these opposing information processing modes are implemented in the brain remains elusive. A potential solution to this problem involves modulating intrinsic… ▽ More In order to remain adaptable to a dynamic environment, neural activity must be simultaneously both sensitive and stable. To solve this problem, the brain has been hypothesised to sit near a critical boundary. Yet, precisely how criticality and these opposing information processing modes are implemented in the brain remains elusive. A potential solution to this problem involves modulating intrinsically nonlinear neurons within the cerebral cortex with neuromodulatory neurotransmitters such as noradrenaline, a highly-conserved chemical released from the pontine locus coeruleus. Here we confirm that neuronal spiking in mice is poised close to the critical point of a branching process and that time-varying signatures of criticality fluctuate with neuromodulatory tone, as assessed by dynamic alterations in pupil diameter. We explore these results theoretically by creating a dual-compartment model of non-linear pyramidal neurons - capable of both regular spike and bursting modes - that replicates our main empirical findings of slightly subcritical dynamics. We then probe our model at a resolution impossible in vivo to demonstrate that noradrenaline differentially alters spiking- and bursting-criticality to facilitate sensitive and stable dynamics following an inverted-U profile that peaks at intermediate noradrenergic tone. Finally, we demonstrate that this intermediate noradrenergic regime displays burst avalanches with power-law size and duration distributions and scaling relationship belonging to the universality class of self-organized criticality. Our results confirm that the noradrenergic ascending arousal system acts as a control parameter for emergent critical dynamics in the brain. This methodology could be extended to explore other neuromodulators as control parameters of the brain. △ Less

Submitted 6 April, 2023; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: 17 pages, 5 figures

Coupling from the Past for the Stochastic Simulation of Chemical Reaction Networks

Authors: J. N. Mueller, J. N. Corcoran

Abstract: Chemical reaction networks (CRNs) are fundamental computational models used to study the behavior of chemical reactions in well-mixed solutions. They have been used extensively to model a broad range of biological systems, and are primarily used when the more traditional model of deterministic continuous mass action kinetics is invalid due to small molecular counts. We present a perfect sampling a… ▽ More Chemical reaction networks (CRNs) are fundamental computational models used to study the behavior of chemical reactions in well-mixed solutions. They have been used extensively to model a broad range of biological systems, and are primarily used when the more traditional model of deterministic continuous mass action kinetics is invalid due to small molecular counts. We present a perfect sampling algorithm to draw error-free samples from the stationary distributions of stochastic models for coupled, linear chemical reaction networks. The state spaces of such networks are given by all permissible combinations of molecular counts for each chemical species, and thereby grow exponentially with the numbers of species in the network. To avoid simulations involving large numbers of states, we propose a subset of chemical species such that coupling of paths started from these states guarantee coupling of paths started from all states in the state space and we show for the well-known Reversible Michaelis-Menten model that the subset does in fact guarantee perfect draws from the stationary distribution of interest. We compare solutions computed in two ways with this algorithm to those found analytically using the chemical master equation and we compare the distribution of coupling times for the two simulation approaches. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Comments: 27 pages, 30 figures

MSC Class: 60J27; 60J28; 60K30 ACM Class: I.6.3; I.6.8

Reinforcement Learning Control of a Biomechanical Model of the Upper Extremity

Authors: Florian Fischer, Miroslav Bachinski, Markus Klar, Arthur Fleig, Jörg Müller

Abstract: Among the infinite number of possible movements that can be produced, humans are commonly assumed to choose those that optimize criteria such as minimizing movement time, subject to certain movement constraints like signal-dependent and constant motor noise. While so far these assumptions have only been evaluated for simplified point-mass or planar models, we address the question of whether they c… ▽ More Among the infinite number of possible movements that can be produced, humans are commonly assumed to choose those that optimize criteria such as minimizing movement time, subject to certain movement constraints like signal-dependent and constant motor noise. While so far these assumptions have only been evaluated for simplified point-mass or planar models, we address the question of whether they can predict reaching movements in a full skeletal model of the human upper extremity. We learn a control policy using a motor babbling approach as implemented in reinforcement learning, using aimed movements of the tip of the right index finger towards randomly placed 3D targets of varying size. We use a state-of-the-art biomechanical model, which includes seven actuated degrees of freedom. To deal with the curse of dimensionality, we use a simplified second-order muscle model, acting at each degree of freedom instead of individual muscles. The results confirm that the assumptions of signal-dependent and constant motor noise, together with the objective of movement time minimization, are sufficient for a state-of-the-art skeletal model of the human upper extremity to reproduce complex phenomena of human movement, in particular Fitts' Law and the 2/3 Power Law. This result supports the notion that control of the complex human biomechanical system can plausibly be determined by a set of simple assumptions and can easily be learned. △ Less

Submitted 20 April, 2022; v1 submitted 13 November, 2020; originally announced November 2020.

Comments: 19 pages, 7 figures

Journal ref: Fischer, F., Bachinski, M., Klar, M., Fleig, A., & Müller, J. (2021). Reinforcement learning control of a biomechanical model of the upper extremity. Scientific Reports, 11(1), 1-15

Contact Tracing & Super-Spreaders in the Branching-Process Model

Authors: Johannes Müller, Volker Hösel

Abstract: In recent years, it became clear that super-spreader events play an important role, particularly in the spread of airborne infections. We investigate a novel model for super-spreader events, not based on a heterogeneous contact graph but on a random contact rate: Many individuals become infected synchronously in single contact events. We use the branching-process approach for contact tracing to an… ▽ More In recent years, it became clear that super-spreader events play an important role, particularly in the spread of airborne infections. We investigate a novel model for super-spreader events, not based on a heterogeneous contact graph but on a random contact rate: Many individuals become infected synchronously in single contact events. We use the branching-process approach for contact tracing to analyze the impact of super-spreader events on the effect of contact tracing. Here we neglect a tracing delay. Roughly speaking, we find that contact tracing is more efficient in the presence of super-spreaders if the fraction of symptomatics is small, the tracing probability is high, or the latency period is distinctively larger than the incubation period. In other cases, the effect of contact tracing can be decreased by super-spreaders. Numerical analysis with parameters suited for SARS-CoV-2 indicates that super-spreaders do not decrease the effect of contact tracing crucially in case of that infection. △ Less

Submitted 10 April, 2021; v1 submitted 10 October, 2020; originally announced October 2020.

MSC Class: Primary 92D30; Secondary 60J80

Contact Tracing -- Old Models and New Challenges

Authors: Johannes Müller, Mirjam Kretzschmar

Abstract: Contact tracing is an effective method to control emerging diseases. Since the 1980's, mathematical modelers are developing a consistent theory for contact tracing, with the aim to find effective and efficient implementations of contact tracing, and to assess the effects of contact tracing on the spread of an infectious disease. Despite the progress made in the area, there remain important open qu… ▽ More Contact tracing is an effective method to control emerging diseases. Since the 1980's, mathematical modelers are developing a consistent theory for contact tracing, with the aim to find effective and efficient implementations of contact tracing, and to assess the effects of contact tracing on the spread of an infectious disease. Despite the progress made in the area, there remain important open questions. In addition, technological developments, especially in the field of molecular biology (genetic sequencing of pathogens) and modern communication (digital contact tracing), have posed new challenges for the modeling community. In the present paper, we discuss modeling approaches for contact tracing and identify some of the current challenges for the field. △ Less

Submitted 15 August, 2020; v1 submitted 23 July, 2020; originally announced July 2020.

Control of ecological outcomes through deliberate parameter changes in a model of the gut microbiome

Authors: Zipeng Wang, Eric W. Jones, Joshua M. Mueller, Jean M. Carlson

Abstract: The generalized Lotka-Volterra (gLV) equations are a mathematical proxy for ecological dynamics. We focus on a gLV model of the gut microbiome, in which the evolution of the gut microbial state is determined in part by pairwise inter-species interaction parameters that encode environmentally-mediated resource competition between microbes. We develop an in silico method that controls the steady-sta… ▽ More The generalized Lotka-Volterra (gLV) equations are a mathematical proxy for ecological dynamics. We focus on a gLV model of the gut microbiome, in which the evolution of the gut microbial state is determined in part by pairwise inter-species interaction parameters that encode environmentally-mediated resource competition between microbes. We develop an in silico method that controls the steady-state outcome of the system by adjusting these interaction parameters. This approach is confined to a bistable region of the gLV model. The two steady states of interest are idealized as either a "healthy" or "diseased" steady state of the gut microbiome. In this method, a dimensionality reduction technique called steady-state reduction (SSR) is first used to generate a two-dimensional (2D) gLV model that approximates the high-dimensional dynamics on the 2D subspace spanned by the two steady states. Then a bifurcation analysis of the 2D model analytically determines parameter modifications that drive a disease-prone initial condition to the healthy steady state. This parameter modification of the reduced 2D model guides parameter modifications of the original high-dimensional model, resulting in a change of steady-state outcome in the high-dimensional model. This control method, called SPARC (SSR-guided parameter change), bypasses the computational challenge of directly determining parameter modifications in the original high-dimensional system. SPARC could guide the development of indirect bacteriotherapies, which seek to change microbial compositions by deliberately modifying gut environmental variables such as gut acidity or macronutrient availability. △ Less

Submitted 29 March, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Comments: main text 9 pages, 5 figures; supplement 10 pages, 2 figures

Journal ref: Phys. Rev. E 101, 052402 (2020)

Exact and approximate formulas for contact tracing on random trees

Authors: Augustine Okolie, Johannes Müller

Abstract: We consider a stochastic susceptible-infected-recovered (SIR) model with contact tracing on random trees and on the configuration model. On a rooted tree, where initially all individuals are susceptible apart from the root which is infected, we are able to find exact formulas for the distribution of the infectious period. Thereto, we show how to extend the existing theory for contact tracing in ho… ▽ More We consider a stochastic susceptible-infected-recovered (SIR) model with contact tracing on random trees and on the configuration model. On a rooted tree, where initially all individuals are susceptible apart from the root which is infected, we are able to find exact formulas for the distribution of the infectious period. Thereto, we show how to extend the existing theory for contact tracing in homogeneously mixing populations to trees. Based on these formulas, we discuss the influence of randomness in the tree and the basic reproduction. We find the well known results for the homogeneously mixing case as a limit of the present model (tree-shaped contact graph). Furthermore, we develop approximate mean field equations for the dynamics on trees, and - using the message passing method - also for the configuration model. The interpretation and implications of the results are discussed. △ Less

Submitted 22 February, 2020; v1 submitted 15 October, 2019; originally announced October 2019.

Comments: 24 pages, 9 figures

MSC Class: 92Bxx (Primary) 92B05; 05Cxx (Secondary)

Journal ref: Mathematical Biosciences, 108320 (2020)

Homogeneous and Heterogeneous Response of Quorum-Sensing Bacteria in an Evolutionary Context

Authors: Burkhard A. Hense, Matthew McIntosh, Johannes Müller, Martin Schuster

Abstract: To explain the stability of cooperation is a central task of evolutionary theory. We investigate this question in the case of quorum sensing (QS) bacteria, which regulate cooperative traits in response to population density. Cooperation is modeled by the prisoner's dilemma, where individuals produce a costly public good (PG) that equally benefits all members of a community divided into multiple, d… ▽ More To explain the stability of cooperation is a central task of evolutionary theory. We investigate this question in the case of quorum sensing (QS) bacteria, which regulate cooperative traits in response to population density. Cooperation is modeled by the prisoner's dilemma, where individuals produce a costly public good (PG) that equally benefits all members of a community divided into multiple, distinct patches (multilevel selection). Cost and benefit are non-linear functions of the PG production. The analysis of evolutionary stability yields an optimization problem for the expression of PG in dependency on the number of QS individuals within a colony. We find that the optimal total PG production of the QS population mainly depends on the shape of the benefit. A graded and a switch-like response is possible, in accordance with earlier results. Interestingly, at the level of the individual cell, the QS response is determined by the shape of the costs. All QS individuals respond either homogeneously if cost are a convex function of the PG production rate, or they respond heterogeneously with distinct ON/OFF responses if the costs are concave. The latter finding is consistent with recent experimental findings, and contradicts the usual interpretation of QS as a mechanism to establish a uniform, synchronized response of a bacterial population. △ Less

Submitted 1 May, 2019; originally announced May 2019.

MSC Class: 92D25; 91A22

Are the better cooperators dormant or quiescent?

Authors: Thibaut Sellinger, Johannes Müller, Volker Hösel, Aurélien Tellier

Abstract: Despite the wealth of empirical and theoretical studies, the origin and maintenance of cooperation is still an evolutionary riddle. In this context, ecological life-history traits which affect the efficiency of selection may play a role, though these are often ignored. We consider here species such as bacteria, fungi, invertebrates and plants which exhibit resting stages in the form of a quiescent… ▽ More Despite the wealth of empirical and theoretical studies, the origin and maintenance of cooperation is still an evolutionary riddle. In this context, ecological life-history traits which affect the efficiency of selection may play a role, though these are often ignored. We consider here species such as bacteria, fungi, invertebrates and plants which exhibit resting stages in the form of a quiescent state or a seedbank. When quiescent, individuals are inactive and reproduce upon activation, while under seed bank parents produce offspring remaining dormant for different amount of time. We assume weak frequency-dependent selection modeled using game-theory and the prisoners dilemma (cooperation/defect) as payoff matrix. The cooperators and defectors are allowed to evolve different quiescence or dormancy times. By means of singular perturbation theory we reduce the model to a one-dimensional equation resembling the well known replicator equation, where the gain functions are scaled with lumped parameters reflecting the time scale of the resting state of the cooperators and defectors. If both time scales are identical cooperation cannot persist in a homogeneous population. If, however, the time scale of the cooperator is distinctively different from that of the defector, cooperation may become a locally asymptotically stable strategy. Interestingly enough, in the seedbank case the cooperator needs to be faster than the defector, while in the quiescent case the cooperator has to be slower. We use adaptive dynamics to identify situations where cooperation may evolve and form a convergent stable ESS. We conclude by highlighting the relevance fo these results for many non-model species and the maintenance of cooperation in microbial, invertebrate or plant populations. △ Less

Submitted 14 April, 2019; originally announced April 2019.

MSC Class: 92D10; 92D25; 34E20

An extended model of wine fermentation including aromas and acids

Authors: Jan Bartsch, Alfio Borzì, Christina Schenk, Dominik Schmidt, Jonas Müller, Volker Schulz, Kai Velten

Abstract: The art of viticulture and the quest for making wines has a long tradition and it just started recently that mathematicians entered this field with their main contribution of modelling alcoholic fermentation. These models consist of systems of ordinary differential equations that describe the kinetics of the bio-chemical reactions occurring in the fermentation process. The aim of this paper is to… ▽ More The art of viticulture and the quest for making wines has a long tradition and it just started recently that mathematicians entered this field with their main contribution of modelling alcoholic fermentation. These models consist of systems of ordinary differential equations that describe the kinetics of the bio-chemical reactions occurring in the fermentation process. The aim of this paper is to present a new model of wine fermentation that accurately describes the yeast dying component, the presence of glucose transporters, and the formation of aromas and acids. Therefore the new model could become a valuable tool to predict the taste of the wine and provide the starting point for an emerging control technology that aims at improving the quality of the wine by steering a well-behaved fermentation process that is also energetically more efficient. Results of numerical simulations are presented that successfully confirm the validity of the proposed model by comparison with real data. △ Less

Submitted 5 January, 2019; originally announced January 2019.

Convolutional Neural Networks for Epileptic Seizure Prediction

Authors: Matthias Eberlein, Raphael Hildebrand, Ronald Tetzlaff, Nico Hoffmann, Levin Kuhlmann, Benjamin Brinkmann, Jens Müller

Abstract: Epilepsy is the most common neurological disorder and an accurate forecast of seizures would help to overcome the patient's uncertainty and helplessness. In this contribution, we present and discuss a novel methodology for the classification of intracranial electroencephalography (iEEG) for seizure prediction. Contrary to previous approaches, we categorically refrain from an extraction of hand-cra… ▽ More Epilepsy is the most common neurological disorder and an accurate forecast of seizures would help to overcome the patient's uncertainty and helplessness. In this contribution, we present and discuss a novel methodology for the classification of intracranial electroencephalography (iEEG) for seizure prediction. Contrary to previous approaches, we categorically refrain from an extraction of hand-crafted features and use a convolutional neural network (CNN) topology instead for both the determination of suitable signal characteristics and the binary classification of preictal and interictal segments. Three different models have been evaluated on public datasets with long-term recordings from four dogs and three patients. Overall, our findings demonstrate the general applicability. In this work we discuss the strengths and limitations of our methodology. △ Less

Submitted 11 April, 2023; v1 submitted 2 November, 2018; originally announced November 2018.

Comments: accepted for MLESP 2018

Journal ref: 2018 IEEE International Conference on Bioinformatics and Biomedicine (BIBM)

Effects of population- and seed bank noise on neutral evolution and efficacy of natural selection

Authors: Lukas Heinrich, Johannes Müller, Aurélien Tellier, Daniel Zivković

Abstract: Population genetics models typically consider a fixed population size and a unique selection coefficient. However, population dynamics inherently generate noise in numbers of individuals and selection acts on various components of the individuals' fitness. In plant species with seed banks, the size of both the above- and below-ground compartments present noise depending on seed production and the… ▽ More Population genetics models typically consider a fixed population size and a unique selection coefficient. However, population dynamics inherently generate noise in numbers of individuals and selection acts on various components of the individuals' fitness. In plant species with seed banks, the size of both the above- and below-ground compartments present noise depending on seed production and the state of the seed bank. We investigate if this noise has consequences on 1)~the rate of genetic drift, and 2)~the efficacy of selection. We consider four variants of two-allele Moran-type models defined by combinations of presence and absence of noise in above-ground and seed bank compartments. Time scale analysis and dimension reduction methods allow us to reduce the corresponding Fokker-Planck equation to a one-dimensional diffusion approximation of a Moran model. We first show that if the above-ground noise classically affects the rate of genetic drift, below-ground noise reduces the diversity storage effect of the seed bank. Second, we consider that selection can act on four different components of the plant fitness: plant or seed death rate, seed production or seed germination. Our striking result is that the efficacy of selection for seed death rate or germination rate is reduced by seed bank noise, whereas selection occurring on plant death rate or seed production is not affected. We derive the expected site-frequency spectrum reflecting this heterogeneity in selection efficacy between genes underpinning different plant fitness components. Our results highlight the importance to consider the effect of ecological noise to predict the impact of seed banks on neutral and selective evolution. △ Less

Submitted 11 December, 2017; v1 submitted 26 July, 2017; originally announced July 2017.

MSC Class: 92D10; 39A50; 60H10

Global stability properties of renewal epidemic models

Authors: Michael T. Meehan, Daniel G. Cocks, Johannes Müller, Emma S. McBryde

Abstract: We investigate the global dynamics of a general Kermack-McKendrick-type epidemic model formulated in terms of a system of renewal equations. Specifically, we consider a renewal model for which both the force of infection and the infected removal rates are arbitrary functions of the infection age, $τ$, and use the direct Lyapunov method to establish the global asymptotic stability of the equilibriu… ▽ More We investigate the global dynamics of a general Kermack-McKendrick-type epidemic model formulated in terms of a system of renewal equations. Specifically, we consider a renewal model for which both the force of infection and the infected removal rates are arbitrary functions of the infection age, $τ$, and use the direct Lyapunov method to establish the global asymptotic stability of the equilibrium solutions. In particular, we show that the basic reproduction number, $R_0$, represents a sharp threshold parameter such that for $R_0\leq 1$, the infection-free equilibrium is globally asymptotically stable; whereas the endemic equilibrium becomes globally asymptotically stable when $R_0 > 1$, i.e. when it exists. △ Less

Submitted 2 October, 2018; v1 submitted 11 July, 2017; originally announced July 2017.

Comments: 11 pages

Plasmid segregation and accumulation

Authors: Johannes Müller, Karin Münch, Bendix Koopmann, Eva Stadler, Louisa Roselius, Dieter Jahn, Richard Münch

Abstract: The segregation of plasmids in a bacterial population is investigated. Hereby, a dynamical model is formulated in terms of a size-structured population using a hyperbolic partial differential equation incorporating non-local terms (the fragmentation equation). For a large class of parameter functions this PDE can be re-written as an infinite system of ordinary differential equations for the moment… ▽ More The segregation of plasmids in a bacterial population is investigated. Hereby, a dynamical model is formulated in terms of a size-structured population using a hyperbolic partial differential equation incorporating non-local terms (the fragmentation equation). For a large class of parameter functions this PDE can be re-written as an infinite system of ordinary differential equations for the moments of its solution. We investigate the influence of different plasmid production modes, kinetic parameters, and plasmid segregation modes on the equilibrium plasmid distribution. In particular, at small plasmid numbers the distribution is strongly influenced by the production mode, while the kinetic parameters (cell growth rate resp. basic plasmid reproduction rate) influence the distribution mainly at large plasmid numbers. The plasmid transmission characteristics only gradually influence the distribution, but may become of importance for biologically relevant cases. We compare the theoretical findings with experimental results. △ Less

Submitted 12 January, 2017; originally announced January 2017.

MSC Class: 92D25; 35L02

The effect of delay on contact tracing

Authors: Johannes Müller, Bendix Koopmann

Abstract: We consider a model for an infectious disease in the onset of an outbreak. We introduce contact tracing incorporating a tracing delay. The effect of randomness in the delay and the effect of the length of this delay in comparison to the infectious period of the disease respectively to a latency period on the effect of tracing, given e.g. by the change of the reproduction number, is analyzed. We fo… ▽ More We consider a model for an infectious disease in the onset of an outbreak. We introduce contact tracing incorporating a tracing delay. The effect of randomness in the delay and the effect of the length of this delay in comparison to the infectious period of the disease respectively to a latency period on the effect of tracing, given e.g. by the change of the reproduction number, is analyzed. We focus particularly on the effect of randomness in the tracing delay. △ Less

Submitted 22 October, 2016; originally announced October 2016.

MSC Class: 92D30; 60J85

Fisher-Wright model with deterministic seed bank and selection

Authors: Bendix Koopmann, Johannes Müller, Aurélien Tellier, Daniel Živković

Abstract: Seed banks are a common characteristics to many plant species, which allow storage of genetic diversity in the soil as dormant seeds for various periods of time. We investigate an above-ground population following a Fisher-Wright model with selection coupled with a deterministic seed bank assuming the length of the seed bank is kept constant and the number of seeds is large. To assess the combined… ▽ More Seed banks are a common characteristics to many plant species, which allow storage of genetic diversity in the soil as dormant seeds for various periods of time. We investigate an above-ground population following a Fisher-Wright model with selection coupled with a deterministic seed bank assuming the length of the seed bank is kept constant and the number of seeds is large. To assess the combined impact of seed banks and selection on genetic diversity, we derive a general diffusion model. The applied techniques outline a path of approximating a stochastic delay differential equation by an appropriately rescaled stochastic differential equation, which is a common issue in statistical physics. We compute the equilibrium solution of the site-frequency spectrum and derive the times to fixation of an allele with and without selection. Finally, it is demonstrated that seed banks enhance the effect of selection onto the site-frequency spectrum while slowing down the time until the mutation-selection equilibrium is reached. △ Less

Submitted 28 July, 2016; v1 submitted 20 May, 2016; originally announced May 2016.

MSC Class: 92D15

Transcranial Focused Ultrasound for BOLD fMRI Signal Modulation in Humans

Authors: Leo Ai, Jerel K. Mueller, Andrea Grant, Yigitcan Eryaman, Wynn Legon

Abstract: Transcranial focused ultrasound (tFUS) is an emerging form of non-surgical human neuromodulation that confers advantages over existing electro and electromagnetic technologies by providing a superior spatial resolution on the millimeter scale as well as the capability to target sub-cortical structures non-invasively. An examination of the pairing of tFUS and blood oxygen level dependent (BOLD) fun… ▽ More Transcranial focused ultrasound (tFUS) is an emerging form of non-surgical human neuromodulation that confers advantages over existing electro and electromagnetic technologies by providing a superior spatial resolution on the millimeter scale as well as the capability to target sub-cortical structures non-invasively. An examination of the pairing of tFUS and blood oxygen level dependent (BOLD) functional MRI (fMRI) in humans is presented here. △ Less

Submitted 1 March, 2016; originally announced March 2016.

Comments: EMBC 2016

Modeling Persistent Trends in Distributions

Authors: Jonas Mueller, Tommi Jaakkola, David Gifford

Abstract: We present a nonparametric framework to model a short sequence of probability distributions that vary both due to underlying effects of sequential progression and confounding noise. To distinguish between these two types of variation and estimate the sequential-progression effects, our approach leverages an assumption that these effects follow a persistent trend. This work is motivated by the rece… ▽ More We present a nonparametric framework to model a short sequence of probability distributions that vary both due to underlying effects of sequential progression and confounding noise. To distinguish between these two types of variation and estimate the sequential-progression effects, our approach leverages an assumption that these effects follow a persistent trend. This work is motivated by the recent rise of single-cell RNA-sequencing experiments over a brief time course, which aim to identify genes relevant to the progression of a particular biological process across diverse cell populations. While classical statistical tools focus on scalar-response regression or order-agnostic differences between distributions, it is desirable in this setting to consider both the full distributions as well as the structure imposed by their ordering. We introduce a new regression model for ordinal covariates where responses are univariate distributions and the underlying relationship reflects consistent changes in the distributions over increasing levels of the covariate. This concept is formalized as a "trend" in distributions, which we define as an evolution that is linear under the Wasserstein metric. Implemented via a fast alternating projections algorithm, our method exhibits numerous strengths in simulations and analyses of single-cell gene expression data. △ Less

Submitted 24 May, 2017; v1 submitted 13 November, 2015; originally announced November 2015.

Comments: To appear in: Journal of the American Statistical Association

Journal ref: Journal of the American Statistical Association, 113(523):1296-1310, 2018

Analysis of Transcranial Focused Ultrasound Beam Profile Sensitivity for Neuromodulation of the Human Brain

Authors: Jerel K. Mueller, Wynn Legon, William J. Tyler

Abstract: Objective. While ultrasound is largely established for use in diagnostic imaging and heating therapies, its application for neuromodulation is relatively new and not well understood. The objective of the present study was to investigate issues related to interactions between focused acoustic beams and brain tissues to better understand possible limitations of transcranial ultrasound for neuromodul… ▽ More Objective. While ultrasound is largely established for use in diagnostic imaging and heating therapies, its application for neuromodulation is relatively new and not well understood. The objective of the present study was to investigate issues related to interactions between focused acoustic beams and brain tissues to better understand possible limitations of transcranial ultrasound for neuromodulation. Approach. A computational model of transcranial focused ultrasound was constructed and validated against bench top experimental data. The models were then incrementally extended to address and investigate a number of issues related to the use of ultrasound for neuromodulation. These included the effect of variations in skull geometry and gyral anatomy, as well as the effect of transmission across multiple tissue and media layers, such as scalp, skull, CSF, and gray/white matter on ultrasound insertion behavior. In addition, a sensitivity analysis was run to characterize the influence of acoustic properties of intracranial tissues. Finally, the heating associated with ultrasonic stimulation waveforms designed for neuromodulation was modeled. Main results. Depending on factors such as acoustic frequency, the insertion behavior of a transcranial focused ultrasound beam is only subtly influenced by the geometry and acoustic properties of the underlying tissues. Significance. These issues are critical for the refinement of device design and the overall advancement of ultrasound methods for noninvasive neuromodulation. △ Less

Submitted 6 March, 2015; originally announced March 2015.

Comments: 32 pages, 14 figures, and 3 tables

Novel model for wine fermentation including the yeast dying phase

Authors: Alfio Borzì, Juri Merger, Jonas Müller, Achim Rosch, Christina Schenk, Dominik Schmidt, Stephan Schmidt, Volker Schulz, Kai Velten, Christian von Wallbrunn, Michael Zänglein

Abstract: This paper presents a novel model for wine fermentation including a death phase for yeast and the influence of oxygen on the process. A model for the inclusion of the yeast dying phase is derived and compared to a model taken from the literature. The modeling ability of the several models is analyzed by comparing their simulation results. This paper presents a novel model for wine fermentation including a death phase for yeast and the influence of oxygen on the process. A model for the inclusion of the yeast dying phase is derived and compared to a model taken from the literature. The modeling ability of the several models is analyzed by comparing their simulation results. △ Less

Submitted 17 December, 2014; originally announced December 2014.

Comments: 6 pages, 2 figures

Genetic drift opposes mutualism during spatial population expansion

Authors: Melanie JI Muller, Beverly I Neugeboren, David R Nelson, Andrew W Murray

Abstract: Mutualistic interactions benefit both partners, promoting coexistence and genetic diversity. Spatial structure can promote cooperation, but spatial expansions may also make it hard for mutualistic partners to stay together, since genetic drift at the expansion front creates regions of low genetic and species diversity. To explore the antagonism between mutualism and genetic drift, we grew cross-fe… ▽ More Mutualistic interactions benefit both partners, promoting coexistence and genetic diversity. Spatial structure can promote cooperation, but spatial expansions may also make it hard for mutualistic partners to stay together, since genetic drift at the expansion front creates regions of low genetic and species diversity. To explore the antagonism between mutualism and genetic drift, we grew cross-feeding strains of the budding yeast S. cerevisiae on agar surfaces as a model for mutualists undergoing spatial expansions. By supplying varying amounts of the exchanged nutrients, we tuned strength and symmetry of the mutualistic interaction. Strong mutualism suppresses genetic demixing during spatial expansions and thereby maintains diversity, but weak or asymmetric mutualism is overwhelmed by genetic drift even when mutualism is still beneficial, slowing growth and reducing diversity. Theoretical modeling using experimentally measured parameters predicts the size of demixed regions and how strong mutualism must be to survive a spatial expansion. △ Less

Submitted 24 February, 2014; originally announced February 2014.

Comments: 12 pages main text, 20 pages SI

Journal ref: Proc. Natl. Acad. Sci. USA 111, 1037-1042 (2014)

Spatial population expansion promotes the evolution of cooperation in an experimental Prisoner's Dilemma

Authors: J. David Van Dyken, Melanie J. I. Muller, Keenan M. L. Mack, Michael M. Desai

Abstract: Cooperation is ubiquitous in nature, but explaining its existence remains a central interdisciplinary challenge. Cooperation is most difficult to explain in the Prisoner's Dilemma game, where cooperators always lose in direct competition with defectors despite increasing mean fitness. Here we demonstrate how spatial population expansion, a widespread natural phenomenon, promotes the evolution of c… ▽ More Cooperation is ubiquitous in nature, but explaining its existence remains a central interdisciplinary challenge. Cooperation is most difficult to explain in the Prisoner's Dilemma game, where cooperators always lose in direct competition with defectors despite increasing mean fitness. Here we demonstrate how spatial population expansion, a widespread natural phenomenon, promotes the evolution of cooperation. We engineer an experimental Prisoner's Dilemma game in the budding yeast Saccharomyces cerevisiae to show that, despite losing to defectors in nonexpanding conditions, cooperators increase in frequency in spatially expanding populations. Fluorescently labeled colonies show genetic demixing of cooperators and defectors, followed by increase in cooperator frequency as cooperator sectors overtake neighboring defector sectors. Together with lattice-based spatial simulations, our results suggest that spatial population expansion drives the evolution of cooperation by (1) increasing positive genetic assortment at population frontiers and (2) selecting for phenotypes maximizing local deme productivity. Spatial expansion thus creates a selective force whereby cooperator-enriched demes overtake neighboring defector-enriched demes in a "survival of the fastest". We conclude that colony growth alone can promote cooperation and prevent defection in microbes. Our results extend to other species with spatially restricted dispersal undergoing range expansion, including pathogens, invasive species, and humans. △ Less

Submitted 11 November, 2013; originally announced November 2013.

Journal ref: Curr Biol 23, 919-923, 2013 Export Citation Permissions Current Biology, Volume 23, Issue 10, 919-923

General triallelic frequency spectrum under demographic models with variable population size

Authors: Paul A. Jenkins, Jonas W. Mueller, Yun S. Song

Abstract: It is becoming routine to obtain datasets on DNA sequence variation across several thousands of chromosomes, providing unprecedented opportunity to infer the underlying biological and demographic forces. Such data make it vital to study summary statistics which offer enough compression to be tractable, while preserving a great deal of information. One well-studied summary is the site frequency spe… ▽ More It is becoming routine to obtain datasets on DNA sequence variation across several thousands of chromosomes, providing unprecedented opportunity to infer the underlying biological and demographic forces. Such data make it vital to study summary statistics which offer enough compression to be tractable, while preserving a great deal of information. One well-studied summary is the site frequency spectrum---the empirical distribution, across segregating sites, of the sample frequency of the derived allele. However, most previous theoretical work has assumed that each site has experienced at most one mutation event in its genealogical history, which becomes less tenable for very large sample sizes. In this work we obtain, in closed-form, the predicted frequency spectrum of a site that has experienced at most two mutation events, under very general assumptions about the distribution of branch lengths in the underlying coalescent tree. Among other applications, we obtain the frequency spectrum of a triallelic site in a model of historically varying population size. We demonstrate the utility of our formulas in two settings: First, we show that triallelic sites are more sensitive to the parameters of a population that has experienced historical growth, suggesting that they will have use if they can be incorporated into demographic inference. Second, we investigate a recently proposed alternative mechanism of mutation in which the two derived alleles of a triallelic site are created simultaneously within a single individual, and we develop a test to determine whether it is responsible for the excess of triallelic sites in the human genome. △ Less

Submitted 25 November, 2013; v1 submitted 12 October, 2013; originally announced October 2013.

Comments: 29 pages, 11 figures (main text) + 6 pages, 3 figures (Supporting Information)

Journal ref: Genetics, Vol. 196, No. 1 (2014) 295-311

A Quantitative Overview of Biophysical Forces Governing Neural Function

Authors: Jerel Mueller, William J. Tyler

Abstract: The Hodgkin-Huxley (HH) model is the currently accepted formalism of neuronal excitability. However, the HH model does not capture a number of biophysical behaviors associated with action potentials or propagating nerve impulses. Physical mechanisms underlying these processes, such as reversible heat transfer and axonal swelling have been separately investigated and compartmentally modeled to indi… ▽ More The Hodgkin-Huxley (HH) model is the currently accepted formalism of neuronal excitability. However, the HH model does not capture a number of biophysical behaviors associated with action potentials or propagating nerve impulses. Physical mechanisms underlying these processes, such as reversible heat transfer and axonal swelling have been separately investigated and compartmentally modeled to indicate the nervous system is not purely electrical or biochemical. Rather, mechanical forces and principles of thermodynamics also govern neuronal excitability and signaling. To advance our understanding of neural function and dysfunction, compartmentalized analyses of electrical, chemical, and mechanical processes need to revaluated and integrated into more comprehensive theories. The present quantitative perspective is intended to broaden the awareness of known biophysical phenomena, which are often overlooked in neuroscience. By starting to consider the collective influence of the biophysical forces influencing neural function, new paradigms can be applied to the characterization and manipulation of nervous systems. △ Less

Submitted 25 September, 2013; v1 submitted 24 September, 2013; originally announced September 2013.

Comments: 13 pages

MSC Class: 92C20

Supplementary information for "Selective sweeps in growing microbial colonies"

Authors: Kirill S Korolev, Melanie J I Müller, Nilay Karahan, Andrew W Murray, Oskar Hallatschek, David R Nelson

Abstract: This is supplementary information for arXiv:1204.4896, which also appeared in Physical Biology 9, 026008 (2012). Main text abstract: Evolutionary experiments with microbes are a powerful tool to study mutations and natural selection. These experiments, however, are often limited to the well-mixed environments of a test tube or a chemostat. Since spatial organization can significantly affect evolut… ▽ More This is supplementary information for arXiv:1204.4896, which also appeared in Physical Biology 9, 026008 (2012). Main text abstract: Evolutionary experiments with microbes are a powerful tool to study mutations and natural selection. These experiments, however, are often limited to the well-mixed environments of a test tube or a chemostat. Since spatial organization can significantly affect evolutionary dynamics, the need is growing for evolutionary experiments in spatially structured environments. The surface of a Petri dish provides such an environment, but a more detailed understanding of microbial growth on Petri dishes is necessary to interpret such experiments. We formulate a simple deterministic reaction-diffusion model, which successfully predicts the spatial patterns created by two competing species during colony expansion. We also derive the shape of these patterns analytically without relying on microscopic details of the model. In particular, we find that the relative fitness of two microbial strains can be estimated from the logarithmic spirals created by selective sweeps. The theory is tested with strains of the budding yeast Saccharomyces cerevisiae, for spatial competitions with different initial conditions and for a range of relative fitnesses. The reaction-diffusion model also connects the microscopic parameters like growth rates and diffusion constants with macroscopic spatial patterns and predicts the relationship between fitness in liquid cultures and on Petri dishes, which we confirmed experimentally. Spatial sector patterns therefore provide an alternative fitness assay to the commonly used liquid culture fitness assays. △ Less

Submitted 27 April, 2012; originally announced April 2012.

Comments: The main text is available as arXiv:1204.4896

Journal ref: Supplementary information for Physical Biology 9, 026008 (2012)

Selective sweeps in growing microbial colonies

Authors: Kirill S Korolev, Melanie J I Müller, Nilay Karahan, Andrew W Murray, Oskar Hallatschek, David R Nelson

Abstract: Evolutionary experiments with microbes are a powerful tool to study mutations and natural selection. These experiments, however, are often limited to the well-mixed environments of a test tube or a chemostat. Since spatial organization can significantly affect evolutionary dynamics, the need is growing for evolutionary experiments in spatially structured environments. The surface of a Petri dish p… ▽ More Evolutionary experiments with microbes are a powerful tool to study mutations and natural selection. These experiments, however, are often limited to the well-mixed environments of a test tube or a chemostat. Since spatial organization can significantly affect evolutionary dynamics, the need is growing for evolutionary experiments in spatially structured environments. The surface of a Petri dish provides such an environment, but a more detailed understanding of microbial growth on Petri dishes is necessary to interpret such experiments. We formulate a simple deterministic reaction-diffusion model, which successfully predicts the spatial patterns created by two competing species during colony expansion. We also derive the shape of these patterns analytically without relying on microscopic details of the model. In particular, we find that the relative fitness of two microbial strains can be estimated from the logarithmic spirals created by selective sweeps. The theory is tested with strains of the budding yeast Saccharomyces cerevisiae, for spatial competitions with different initial conditions and for a range of relative fitnesses. The reaction-diffusion model also connects the microscopic parameters like growth rates and diffusion constants with macroscopic spatial patterns and predicts the relationship between fitness in liquid cultures and on Petri dishes, which we confirmed experimentally. Spatial sector patterns therefore provide an alternative fitness assay to the commonly used liquid culture fitness assays. △ Less

Submitted 1 May, 2012; v1 submitted 22 April, 2012; originally announced April 2012.

Comments: Supplementary information available at arXiv:1204.6328

Journal ref: Physical Biology 9, 026008 (2012)

Traffic by small teams of molecular motors

Authors: Melanie J. I. Müller, Janina Beeg, Rumiana Dimova, Stefan Klumpp, Reinhard Lipowsky

Abstract: Molecular motors transport various cargos along cytoskeletal filaments, analogous to trucks on roads. In contrast to vehicles, however, molecular motors do not work alone but in small teams. We describe a simple model for the transport of a cargo by one team of motors and by two teams of motors, which walk into opposite directions. The cooperation of one team of motors generates long-range trans… ▽ More Molecular motors transport various cargos along cytoskeletal filaments, analogous to trucks on roads. In contrast to vehicles, however, molecular motors do not work alone but in small teams. We describe a simple model for the transport of a cargo by one team of motors and by two teams of motors, which walk into opposite directions. The cooperation of one team of motors generates long-range transport, which we observed experimentally in vitro. Transport by two teams of motors leads to a variety of bidirectional motility behaviour and to dynamic instabilities reminiscent of spontaneous symmetry breaking. We also discuss how cargo transport by teams of motors allows the cell to generate robust long-range bidirectional transport. △ Less

Submitted 7 July, 2008; originally announced July 2008.

Comments: 6 pages, latex, 1 figure, to appear in the Proceedings of "Traffic and Granular Flow '07" (Springer)

Tug-of-war as a cooperative mechanism for bidirectional cargo transport by molecular motors

Authors: Melanie J. I. Müller, Stefan Klumpp, Reinhard Lipowsky

Abstract: Intracellular transport is based on molecular motors that pull cargos along cytoskeletal filaments. One motor species always moves in one direction, e.g. conventional kinesin moves to the microtubule plus end, while cytoplasmic dynein moves to the microtubule minus end. However, many cellular cargos are observed to move bidirectionally, involving both plus-end and minus-end directed motors. The… ▽ More Intracellular transport is based on molecular motors that pull cargos along cytoskeletal filaments. One motor species always moves in one direction, e.g. conventional kinesin moves to the microtubule plus end, while cytoplasmic dynein moves to the microtubule minus end. However, many cellular cargos are observed to move bidirectionally, involving both plus-end and minus-end directed motors. The presumably simplest mechanism for such bidirectional transport is provided by a tug-of-war between the two motor species. This mechanism is studied theoretically using the load-dependent transport properties of individual motors as measured in single-molecule experiments. In contrast to previous expectations, such a tug-of-war is found to be highly cooperative and to exhibit seven different motility regimes depending on the precise values of the single motor parameters. The sensitivity of the transport process to small parameter changes can be used by the cell to regulate its cargo traffic. △ Less

Submitted 3 July, 2008; originally announced July 2008.

Comments: 17 pages, latex, 11 figures, 4 tables, includes Supporting Information

Journal ref: Proc. Natl. Acad. Sci. USA 105, 4609-4614 (2008)

Cooperative transport by small teams of molecular motors

Authors: Stefan Klumpp, Melanie J. I. Muller, Reinhard Lipowsky

Abstract: Molecular motors power directed transport of cargoes within cells. Even if a single motor is sufficient to transport a cargo, motors often cooperate in small teams. We discuss the cooperative cargo transport by several motors theoretically and explore some of its properties. In particular we emphasize how motor teams can drag cargoes through a viscous environment. Molecular motors power directed transport of cargoes within cells. Even if a single motor is sufficient to transport a cargo, motors often cooperate in small teams. We discuss the cooperative cargo transport by several motors theoretically and explore some of its properties. In particular we emphasize how motor teams can drag cargoes through a viscous environment. △ Less

Submitted 4 December, 2006; originally announced December 2006.

Comments: 9 pages, 1 figure, uses ws-brl.cls, presented at Bio-Systems conference, Berlin, June 2006

Journal ref: Biophysical Reviews and Letters 1, 353-361 (2006)

Traffic of Molecular Motors

Authors: Stefan Klumpp, Melanie J. I. Müller, Reinhard Lipowsky

Abstract: Molecular motors perform active movements along cytoskeletal filaments and drive the traffic of organelles and other cargo particles in cells. In contrast to the macroscopic traffic of cars, however, the traffic of molecular motors is characterized by a finite walking distance (or run length) after which a motor unbinds from the filament along which it moves. Unbound motors perform Brownian moti… ▽ More Molecular motors perform active movements along cytoskeletal filaments and drive the traffic of organelles and other cargo particles in cells. In contrast to the macroscopic traffic of cars, however, the traffic of molecular motors is characterized by a finite walking distance (or run length) after which a motor unbinds from the filament along which it moves. Unbound motors perform Brownian motion in the surrounding aqueous solution until they rebind to a filament. We use variants of driven lattice gas models to describe the interplay of their active movements, the unbound diffusion, and the binding/unbinding dynamics. If the motor concentration is large, motor-motor interactions become important and lead to a variety of cooperative traffic phenomena such as traffic jams on the filaments, boundary-induced phase transitions, and spontaneous symmetry breaking in systems with two species of motors. If the filament is surrounded by a large reservoir of motors, the jam length, i.e., the extension of the traffic jams is of the order of the walking distance. Much longer jams can be found in confined geometries such as tube-like compartments. △ Less

Submitted 6 December, 2005; originally announced December 2005.

Comments: 10 pages, latex, uses Springer styles (included), to appear in the Proceedings of "Traffic and Granular Flow 2005"

Journal ref: Traffic and Granular Flow '05, edited by A. Schadschneider et al. (Springer, Berlin, 2007), pp. 251-261

Thermosynthesis as energy source for the RNA World: a new model for the origin of life

Authors: Anthonie W. J. Muller

Abstract: The thermosynthesis concept, biological free energy gain from thermal cycling, is combined with the concept of the RNA World. The resulting overall origin of life model gives new explanations for the emergence of the genetic code and the ribosome. The first protein named pF1 obtains the energy to support the RNA world by a thermal variation of F1 ATP synthase's binding change mechanism. This pF1… ▽ More The thermosynthesis concept, biological free energy gain from thermal cycling, is combined with the concept of the RNA World. The resulting overall origin of life model gives new explanations for the emergence of the genetic code and the ribosome. The first protein named pF1 obtains the energy to support the RNA world by a thermal variation of F1 ATP synthase's binding change mechanism. This pF1 is the single translation product during the emergence of the genetic machinery. During thermal cycling pF1 condenses many substrates with broad specificity, yielding NTPs and randomly constituted protein and RNA libraries that contain (self)-replicating RNA. The smallness of pF1 permits the emergence of the genetic machinery by selection of RNA that increases the fraction of pF1s in the protein library: (1) a progenitor of rRNA that concatenates amino acids bound to (2) a chain of 'positional tRNAs' linked by mutual recognition, yielding a pF1 (or its main motif); this positional tRNA set gradually evolves to a set of regular tRNAs functioning according to the genetic code, with concomitant emergence of (3) an mRNA coding for pF1. △ Less

Submitted 30 March, 2005; v1 submitted 11 January, 2005; originally announced January 2005.

Comments: 12 pages, 6 figures; changes in text (minor) and 2 figures

Morphological Instability and Dynamics of Fronts in Bacterial Growth Models with Nonlinear Diffusion

Authors: J. Mueller, W. van Saarloos

Abstract: It has been argued that there is biological and modeling evidence that a non-linear diffusion coefficient of the type D(b) = D_0 b^{k} underlies the formation of a number of growth patterns of bacterial colonies. We study a reaction-diffusion system with a non-linear diffusion coefficient introduced by Ben-Jacob et al. Due to the fact that the bacterial diffusion coefficient vanishes when the ba… ▽ More It has been argued that there is biological and modeling evidence that a non-linear diffusion coefficient of the type D(b) = D_0 b^{k} underlies the formation of a number of growth patterns of bacterial colonies. We study a reaction-diffusion system with a non-linear diffusion coefficient introduced by Ben-Jacob et al. Due to the fact that the bacterial diffusion coefficient vanishes when the bacterial density b -> 0, the standard linear stability analysis for fronts cannot be used. We introduce an extension of the stability analysis which can be applied to such singular fronts, map out the region of stability in the D-k-plane and derive an interfacial approximation in some limits. Our linear stability analysis and sharp interface formulation will also be applicable to other examples of interface formation due to nonlinear diffusion, like in porous media or in the problem of vortex motion in superconductors. △ Less

Submitted 21 November, 2002; originally announced November 2002.

Journal ref: Phys.Rev.E65:061111,2002