Source author record

Namiko Mitarai

Namiko Mitarai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2020arXiv

Dynamics of bacterial populations under the feast-famine cycles

Bacterial populations in natural conditions are expected to experience stochastic environmental fluctuations, and in addition, environments are affected by bacterial activities since they consume substrates and excrete various chemicals. We here study possible outcomes of population dynamics and evolution under the repeated cycle of substrate-rich conditions and starvation, called the "feast-famine cycle", by a simple stochastic model with the trade-off relationship between the growth rate and the growth yield or the death rate. In the model, the feast (substrate-rich) period is led by a stochastic substrate addition event, while the famine (starvation) period is evoked because bacteria use the supplied substrate. Under the repeated feast-famine cycle, the bacterial population tends to increase the growth rate, even though that tends to decrease the total population size due to the trade-off. Analysis of the model shows that the ratio between the growth rate and the death rate becomes the effective fitness of the population. Hence, the functional form of the trade-off between the growth and death rate determines if the bacterial population eventually goes extinct as an evolutionary consequence. We then show that the increase of the added substrate in the feast period can drive the extinction faster. Overall, the model sheds light on non-trivial possible outcomes under repeated feast-famine cycles.

preprint2020arXiv

When to wake up? The optimal waking-up strategies for starvation-induced persistence

Prolonged lag time can be induced by starvation contributing to the antibiotic tolerance of bacteria. We analyze the optimal lag time to survive and grow the iterative and stochastic application of antibiotics. A simple model shows that the optimal lag time exhibits a discontinuous transition when the severeness of the antibiotic is increased. This suggests the possibility of reducing tolerant bacteria by controlled usage of antibiotics application. When the bacterial populations are able to have two phenotypes with different lag times, the fraction of the second phenotype that has different lag time shows a continuous transition. We then present a generic framework to investigate the optimal lag time distribution for total population fitness for a given distribution of the antibiotic application duration. The obtained optimal distributions have multiple peaks for a wide range of the antibiotic application duration distributions, including the case where the latter is monotonically decreasing. The analysis supports the advantage in evolving multiple, possibly discrete phenotypes in lag time for bacterial long-term fitness.

preprint2016arXiv

Bifurcation of Transition Paths Induced by Coupled Bistable Systems

We discuss the transition paths in a coupled bistable system consisting of interacting multiple identical bistable motifs. We propose a simple model of coupled bistable gene circuits as an example, and show that its transition paths are bifurcating. We then derive a criterion to predict the bifurcation of transition paths in a generalized coupled bistable system. We confirm the validity of the theory for the example system by numerical simulation. We also demonstrate in the example system that, if the steady states of individual gene circuits are not changed by the coupling, the bifurcation pattern is not dependent on the number of gene circuits. We further show that the transition rate exponentially decreases with the number of gene circuits when the transition path does not bifurcate, while a bifurcation facilitates the transition by lowering the quasi-potential energy barrier.

preprint2016arXiv

Three is much more than two in coarsening dynamics of cyclic competitions

The classical game of rock-paper-scissors have inspired experiments and spatial model systems that address robustness of biological diversity. In particular the game nicely illustrates that cyclic interactions allow multiple strategies to coexist for long time intervals. When formulated in terms of a one-dimensional cellular automata, the spatial distribution of strategies exhibits coarsening with algebraically growing domain size over time, while the two-dimensional version allows domains to break and thereby opens for long-time coexistence. We here consider a quasi-one-dimensional implementation of the cyclic competition, and study the long-term dynamics as a function of rare invasions between parallel linear ecosystems. We find that increasing the complexity from two to three parallel subsystems allows a transition from complete coarsening to an active steady state where the domain size stays finite. We further find that this transition happens irrespective of whether the update is done in parallel for all sites simultaneously, or done randomly in sequential order. In both cases the active state is characterized by localized bursts of dislocations, followed by longer periods of coarsening. In the case of the parallel dynamics, we find that there is another phase transition between the active steady state and the coarsening state within the three-line system when the invasion rate between the subsystems is varied. We identify the critical parameter for this transition, and show that the density of active boundaries have critical exponents that are consistent with the directed percolation universality class. On the other hand, numerical simulations with the random sequential dynamics suggest that the system may exhibit an active steady state as long as the invasion rate is finite.

preprint2015arXiv

Dynamic competition between transcription initiation and repression: Role of nonequilibrium steps in cell-to-cell heterogeneity

Transcriptional repression may cause transcriptional noise by a competition between repressor and RNA polymerase binding. Although promoter activity is often governed by a single limiting step, we argue here that the size of the noise strongly depends on whether this step is the initial equilibrium binding or one of the subsequent unidirectional steps. Overall, we show that nonequilibrium steps of transcription initiation systematically increase the cell-to-cell heterogeneity in bacterial populations. In particular, this allows also weak promoters to give substantial transcriptional noise.

preprint2015arXiv

Food web assembly rules

In food webs, many interacting species coexist despite the restrictions imposed by the competitive exclusion principle and apparent competition. For the generalized Lotka-Volterra equations, sustainable coexistence necessitates nonzero determinant of the interaction matrix. Here we show that this requirement is equivalent to demanding that each species be part of a non-overlapping pairing, which substantially constrains the food web structure. We demonstrate that a stable food web can always be obtained if a non-overlapping pairing exists. If it does not, the matrix rank can be used to quantify the lack of niches, corresponding to unpaired species. For the species richness at each trophic level, we derive the food web assembly rules, which specify sustainable combinations. In neighboring levels, these rules allow the higher level to avert competitive exclusion at the lower, thereby incorporating apparent competition. In agreement with data, the assembly rules predict high species numbers at intermediate levels and thinning at the top and bottom. Using comprehensive food web data, we demonstrate how omnivores or parasites with hosts at multiple trophic levels can loosen the constraints and help obtain coexistence in food webs. Hence, omnivory may be the glue that keeps communities intact even under extinction or ecological release of species.

preprint2014arXiv

Disturbance accelerates the transition from low- to high- diversity state in a model ecosystem

The effect of disturbance on a model ecosystem of sessile and mutually competitive species [Mathiesen et al. Phys. Rev. Lett. 107, 188101 (2011); Mitarai et al. Phys. Rev. E 86, 011929 (2012) ] is studied. The disturbance stochastically removes individuals from the system, and the created empty sites are re-colonized by neighbouring species. We show that the stable high-diversity state, maintained by occasional cyclic species interactions that create isolated patches of meta-populations, is robust against small disturbance. We further demonstrate that finite disturbance can accelerate the transition from the low- to high-diversity state by helping creation of small patches through diffusion of boundaries between species with stand-off relation.

preprint2014arXiv

Fixation in cyclically competing species on a directed graph with quenched disorder

A simple model of cyclically competing species on a directed graph with quenched disorder is proposed as an extension of the rock-paper-scissors model. By assuming that the effects of loops in a directed random graph can be ignored in the thermodynamic limit, it is proved for any finite disorder that the system fixates to a frozen configuration when the species number $s$ is larger than the spatial connectivity $c$, and otherwise stays active. Nontrivial lower and upper bounds for the persistence probability of a site never changing its state are also analytically computed. The obtained bounds and numerical simulations support the existence of a phase transition as a function of disorder for $1<c_l(s)\le c <s$, with a $s$-dependent threshold of the connectivity $c_l(s)$.

preprint2013arXiv

Control of ribosome traffic by position-dependent choice of synonymous codons

Messenger RNA encodes a sequence of amino acids by using codons. For most amino acids there are multiple synonymous codons that can encode the amino acid. The translation speed can vary from one codon to another, thus there is room for changing the ribosome speed while keeping the amino acid sequence and hence the resulting protein. Recently, it has been noticed that the choice of the synonymous codon, via the resulting distribution of slow- and fast-translated codons, affects not only on the average speed of one ribosome translating the messenger RNA (mRNA) but also might have an effect on nearby ribosomes by affecting the appearance of "traffic jams" where multiple ribosomes collide and form queues. To test this "context effect" further, we here investigate the effect of the sequence of synonymous codons on the ribosome traffic by using a ribosome traffic model with codon-dependent rates, estimated from experiments. We compare the ribosome traffic on wild type sequences and sequences where the synonymous codons were swapped randomly. By simulating translation of 87 genes, we demonstrate that the wild type sequences, especially those with a high bias in codon usage, tend to have the ability to reduce ribosome collisions, hence optimizing the cellular investment in the translation apparatus. The magnitude of such reduction of the translation time might have a significant impact on the cellular growth rate and thereby have importance for the survival of the species.

preprint2013arXiv

Entrainment of noise-induced and limit cycle oscillators under weak noise

Theoretical models that describe oscillations in biological systems are often either a limit cycle oscillator, where the deterministic nonlinear dynamics gives sustained periodic oscillations, or a noise-induced oscillator, where a fixed point is linearly stable with complex eigenvalues and addition of noise gives oscillations around the fixed point with fluctuating amplitude. We investigate how each class of model behaves under the external periodic forcing, taking the well-studied van der Pol equation as an example. We find that, when the forcing is additive, the noise-induced oscillator can show only one-to-one entrainment to the external frequency, in contrast to the limit cycle oscillator which is known to entrain to any ratio. When the external forcing is multiplicative, on the other hand, the noise-induced oscillator can show entrainment to a few ratios other than one-to-one, while the limit cycle oscillator shows entrain to any ratio. The noise blurs the entrainment in general, but clear entrainment regions for limit cycles can be identified as long as the noise is not too strong.

preprint2013arXiv

Experimental observation of shear thickening oscillation

We report experimental observation of the shear thickening oscillation, i.e. the spontaneous macroscopic oscillation in the shear flow of severe shear thickening fluid. The shear thickening oscillation is caused by the interplay between the fluid dynamics and the shear thickening, and has been predicted theoretically by the present authors using a phenomenological fluid dynamics model for the dilatant fluid, but never been reported experimentally. Using a density-matched starch-water mixture, in the cylindrical shear flow of a few centimeters flow width, we observed strong vibrations of the frequency around 20 Hz, which is consistent with our theoretical prediction.

preprint2012arXiv

Circuit architecture explains functional similarity of bacterial heat shock responses

Heat shock response is a stress response to temperature changes and a consecutive increase in amounts of unfolded proteins. To restore homeostasis, cells upregulate chaperones facilitating protein folding by means of transcription factors (TF). We here investigate two heat shock systems: one characteristic to gram negative bacteria, mediated by transcriptional activator sigma32 in E. coli, and another characteristic to gram positive bacteria, mediated by transcriptional repressor HrcA in L. lactis. We construct simple mathematical model of the two systems focusing on the negative feedbacks, where free chaperons suppress sigma32 activation in the former, while they activate HrcA repression in the latter. We demonstrate that both systems, in spite of the difference at the TF regulation level, are capable of showing very similar heat shock dynamics. We find that differences in regulation impose distinct constrains on chaperone-TF binding affinities: the binding constant of free sigma32 to chaperon DnaK, known to be in 100 nM range, set the lower limit of amount of free chaperon that the system can sense the change at the heat shock, while the binding affinity of HrcA to chaperon GroE set the upper limit and have to be rather large extending into the micromolar range.

preprint2012arXiv

Emergence of diversity in a model ecosystem

The biological requirements for an ecosystem to develop and maintain species diversity are in general unknown. Here we consider a model ecosystem of sessile and mutually excluding organisms competing for space [Mathiesen et al. Phys. Rev. Lett. 107, 188101 (2011)]. The competition is controlled by an interaction network with fixed links chosen by a Bernoulli process. New species are introduced in the system at a predefined rate. In the limit of small introduction rates, the system becomes bistable and can undergo a phase transition from a state of low diversity to high diversity. We suggest that patches of isolated meta-populations formed by the collapse of cyclic relations are essential for the transition to the state of high diversity.

preprint2012arXiv

Information spreading and development of cultural centers

The historical interplay between societies are governed by many factors, including in particular spreading of languages, religion and other symbolic traits. Cultural development, in turn, is coupled to emergence and maintenance of information spreading. Strong centralized cultures exist thanks to attention from their members, which faithfulness in turn relies on supply of information. Here, we discuss a culture evolution model on a planar geometry that takes into account aspects of the feedback between information spreading and its maintenance. Features of model are highlighted by comparing it to cultural spreading in ancient and medieval Europe, where it in particular suggests that long lived centers should be located in geographically remote regions.

preprint2012arXiv

Switching between oscillations and homeostasis in competing negative and positive feedback motifs

We analyze a class of network motifs in which a short, two-node positive feed- back motif is inserted in a three-node negative feedback loop. We demonstrate that such networks can undergo a bifurcation to a state where a stable fixed point and a stable limit cycle coexist. At the bifurcation point the period of the oscillations diverges. Further, intrinsic noise can make the system switch between oscillatory state and the stationary state spontaneously. We find that this switching also occurs in previous models of circadian clocks that use this combination of positive and negative feedback. Our results suggest that real- life circadian systems may need specific regulation to prevent or minimize such switching events.

preprint2011arXiv

Ecosystems with mutually exclusive interactions self-organize to a state of high diversity

Ecological systems comprise an astonishing diversity of species that cooperate or compete with each other forming complex mutual dependencies. The minimum requirements to maintain a large species diversity on long time scales are in general unknown. Using lichen communities as an example, we propose a model for the evolution of mutually excluding organisms that compete for space. We suggest that chain-like or cyclic invasions involving three or more species open for creation of spatially separated sub-populations that subsequently can lead to increased diversity. In contrast to its non-spatial counterpart, our model predicts robust co-existence of a large number of species, in accordance with observations on lichen growth. It is demonstrated that large species diversity can be obtained on evolutionary timescales, provided that interactions between species have spatial constraints. In particular, a phase transition to a sustainable state of high diversity is identified.

preprint2011arXiv

Fluid dynamics of dilatant fluid

Dense mixture of granules and liquid often shows a sever shear thickening and is called a dilatant fluid. We construct a fluid dynamics model for the dilatant fluid by introducing a phenomenological state variable for a local state of dispersed particles. With simple assumptions for an equation of the state variable, we demonstrate that the model can describe basic features of the dilatant fluid such as the stress-shear rate curve that represents discontinuous severe shear thickening, hysteresis upon changing shear rate, instantaneous hardening upon external impact. Analysis of the model reveals that the shear thickening fluid shows an instability in a shear flow for some regime and exhibits {\it the shear thickening oscillation}, i.e. the oscillatory shear flow alternating between the thickened and the relaxed states. Results of numerical simulations are presented for one and two-dimensional systems.

preprint2011arXiv

Modelling the Spatial Dynamics of Culture Spreading in the Presence of Cultural Strongholds

Cultural competition has throughout our history shaped and reshaped the geography of boundaries between humans. Language and culture are intimately connected and linguists often use distinctive keywords to quantify the dynamics of information spreading in societies harbouring strong culture centres. One prominent example, which is addressed here, is Kyoto's historical impact on Japanese culture. We construct a first minimal model, based on shared properties of linguistic maps, to address the interplay between information flow and geography. In particular, we show that spreading of information over Japan in the pre-modern time can be described as a Eden growth process, with noise levels corresponding to coherent spatial patches of sizes given by a single days walk, and with patch-to-patch communication time comparable to the time between human generations.

preprint2010arXiv

Fluid dynamics and jamming in a dilatant fluid

We present a phenomenological fluid dynamics model for a dilatant fluid, i.e. a severe shear thickening fluid, by introducing a state variable. The Navier-Stokes equation is coupled with the state variable field, which evolves in response to the local shear stress as the fluid is sheared. The viscosity is assumed to depend upon the state variable and to diverge at a certain value due to jamming. We demonstrate that the coupling of the fluid dynamics with the shear thickening leads to an oscillatory instability in the shear flow. The model also shows a peculiar response of the fluid to a strong external impact.

preprint2010arXiv

Theory for stability and regulation of epigenetic landscapes

Cells can often choose among several stably heritable phenotypes. Examples are the expression of genes in eukaryotic cells where long chromosomal regions can adopt persistent and heritable silenced or active states, that may be associated with positive feedback in dynamic modification of nucleosomes. We generalize this mechanism in terms of bistability associated with valleys in an epigenetic landscape. A transfer matrix method was used to rigorously follow the system through the disruptive process of cell division. This combined treatment of noisy dynamics both between and during cell division provides an efficient way to calculate the stability of alternative states in a broad range of epigenetic systems.

preprint2009arXiv

Non-ideal behavior of intramolecular structure factor of dilute polymers in a theta solvent

We study the configurational properties of single polymers in a theta solvent by Monte Carlo simulation of the bond fluctuation model. The intramolecular structure factor at the theta point is found to be distinctively different from that of the ideal chain. The structure factor shows a hump around $q\sim 5/R_g$ and a dip around $q\sim 10/R_g$ in the Kratky plot with $R_g$ being the radius of gyration. This feature is apparently similar to that in a melt. The theoretical expression by the simple perturbation expansion to the first order in terms of the Mayer function can be fitted to the obtained structure factor quite well, but the second virial coefficient cannot be set to zero.

preprint2009arXiv

Simple Model for Wet Granular Materials with Liquid Clusters

We propose a simple phenomenological model for wet granular media to take into account many particle interaction through liquid in the funicular state as well as two-body cohesive force by a liquid bridge in the pendular state. In the wet granular media with small liquid content, liquid forms a bridge at each contact point, which induces two-body cohesive force due to the surface tension. As the liquid content increases, some liquid bridges merge, and more than two grains interact through a single liquid cluster. In our model, the cohesive force acts between the grains connected by a liquid-gas interface. As the liquid content increases, the number of grains that interact through the liquid increases, but the liquid-gas interface may decrease when liquid clusters are formed. Due to this competition, our model shows that the shear stress has a maximum as a function of the liquid-content.

preprint2008arXiv

The Generation of Promoter-Mediated Transcriptional Noise in Bacteria

Noise in the expression of a gene produces fluctuations in the concentration of the gene product. These fluctuations can interfere with optimal function or can be exploited to generate beneficial diversity between cells; gene expression noise is therefore expected to be subject to evolutionary pressure. Shifts between modes of high and low rates of transcription initiation at a promoter appear to contribute to this noise both in eukaryotes and prokaryotes. However, models invoked for eukaryotic promoter noise such as stable activation scaffolds or persistent nucleosome alterations seem unlikely to apply to prokaryotic promoters. We consider the relative importance of the steps required for transcription initiation. The 3-step transcription initiation model of McClure is extended into a mathematical model that can be used to predict consequences of additional promoter properties. We show in principle that the transcriptional bursting observed at an E. coli promoter by Golding et al. (2005) can be explained by stimulation of initiation by the negative supercoiling behind a transcribing RNA polymerase (RNAP) or by the formation of moribund or dead-end RNAP-promoter complexes. Both mechanisms are tunable by the alteration of promoter kinetics and therefore allow the optimization of promoter mediated noise.

preprint2007arXiv

Velocity correlations in the dense granular shear flows: Effects on energy dissipation and normal stress

We study the effect of pre-collisional velocity correlations on granular shear flow by molecular dynamics simulations of the inelastic hard sphere system. Comparison of the simulations with the kinetic theory reveals that the theory overestimates both the energy dissipation rate and the normal stress in the dense flow region. We find that the relative normal velocity of colliding particles is smaller than that expected from random collisions, and the discrepancies in the dissipation and the normal stress can be adjusted by introducing the idea of the collisional temperature, from which we conclude that the velocity correlation neglected in the kinetic theory is responsible for the discrepancies. Our analysis of the distributions of the pre-collisional velocity suggests that the correlation grows through multiple inelastic collisions during the time scale of the inverse of the shear rate. As for the shear stress, the discrepancy is also found in the dense region, but it depends strongly on the particle inelasticity.