Publications

A critical cell cycle checkpoint for most bacteria is the onset of constriction when the septal peptidoglycan synthesis starts. According to the current understanding, the arrival of FtsN to midcell triggers this checkpoint in Escherichia coli. Recent structural and in vitro data suggests that recruitment of FtsN to the Z-ring leads to a conformational switch in actin-like FtsA, which links FtsZ protofilaments to the cell membrane and acts as a hub for the late divisome proteins. Here, we investigate this putative pathway using in vivo measurements and stochastic cell cycle modeling at moderately fast growth rates. Quantitatively upregulating protein concentrations and determining the resulting division timings shows that FtsN and FtsA numbers are not rate-limiting for the division in E. coli. However, at higher overexpression levels, they affect divisions: FtsN by accelerating and FtsA by inhibiting them. At the same time, we find that the FtsZ numbers in the cell are one of the rate-limiting factors for cell divisions in E. coli. Altogether, these findings suggest that instead of FtsN, accumulation of FtsZ in the Z-ring is one of the main drivers of the onset of constriction in E. coli at faster growth rates.

Mechanical energy, specifically in the form of ultrasound, can induce pressure variations and temperature fluctuations when applied to an aqueous media. These conditions can both positively and negatively affect protein complexes, consequently altering their stability, folding patterns, and self-assembling behavior. Despite much scientific progress, our current understanding of the effects of ultrasound on the self-assembly of amyloidogenic proteins remains limited. In the present study, we demonstrate that when the amplitude of the delivered ultrasonic energy is sufficiently low, it can induce refolding of specific motifs in protein monomers, which is sufficient for primary nucleation; this has been revealed by MD. These ultrasound-induced structural changes are initiated by pressure perturbations and are accelerated by a temperature factor. Furthermore, the prolonged action of low-amplitude ultrasound enables the elongation of amyloid protein nanofibrils directly from natively folded monomeric lysozyme protein, in a controlled manner, until it reaches a critical length. Using solution X-ray scattering, we determined that nanofibrillar assemblies, formed either under the action of sound or from natively fibrillated lysozyme, share identical structural characteristics. Thus, these results provide insights into the effects of ultrasound on fibrillar protein self-assembly and lay the foundation for the potential use of sound energy in protein chemistry.

Error correction is central to many biological systems and is critical for protein function and cell health. During mitosis, error correction is required for the faithful inheritance of genetic material. When functioning properly, the mitotic spindle segregates an equal number of chromosomes to daughter cells with high fidelity. Over the course of spindle assembly, many initially erroneous attachments between kinetochores and microtubules are fixed through the process of error correction. Despite the importance of chromosome segregation errors in cancer and other diseases, there is a lack of methods to characterize the dynamics of error correction and how it can go wrong. Here, we present an experimental method and analysis framework to quantify chromosome segregation error correction in human tissue culture cells with live cell confocal imaging, timed premature anaphase, and automated counting of kinetochores after cell division. We find that errors decrease exponentially over time during spindle assembly. A coarse-grained model, in which errors are corrected in a chromosome-autonomous manner at a constant rate, can quantitatively explain both the measured error correction dynamics and the distribution of anaphase onset times. We further validated our model using perturbations that destabilized microtubules and changed the initial configuration of chromosomal attachments. Taken together, this work provides a quantitative framework for understanding the dynamics of mitotic error correction.

Cell growth and gene expression, essential elements of all living systems, have long been the focus of biophysical interrogation. Advances in single-cell methods have invigorated theoretical studies into these processes. However, until recently, there was little dialog between the two areas of study. Most theoretical models for gene regulation assumed gene activity to be oblivious to the progression of the cell cycle between birth and division. But there are numerous ways in which the periodic character of all cellular observables can modulate gene expression. The molecular factors required for transcription and translation increase in number during the cell cycle, but are also diluted due to the continuous increase in cell volume. The replication of the genome changes the dosage of those same cellular players but also provides competing targets for regulatory binding. Finally, cell division reduces their number again, and so forth. Stochasticity is inherent to all these biological processes, manifested in fluctuations in the synthesis and degradation of new cellular components as well as the random partitioning of molecules at each cell division. The notion of gene expression as stationary is thus hard to justify. In this review, we survey the emerging paradigm of cell-cycle regulated gene expression, with an emphasis on the global expression patterns rather than gene-specific regulation. We discuss recent experimental reports where cell growth and gene expression were simultaneously measured in individual cells, providing first glimpses into the coupling between the two. While the experimental findings, not surprisingly, differ among genes and organisms, several theoretical models have emerged that attempt to reconcile these differences and form a unifying framework for understanding gene expression in growing cells.

In the past decade, great strides have been made to quantify the dynamics of single-cell growth and division in microbes. In order to make sense of the evolutionary history of these organisms, we must understand how features of single-cell growth and division influence evolutionary dynamics. This requires us to connect processes on the single-cell scale to population dynamics. Here, we consider a model of microbial growth in finite populations which explicitly incorporates the single-cell dynamics. We study the behavior of a mutant population in such a model and ask: can the evolutionary dynamics be coarse-grained so that the forces of natural selection and genetic drift can be expressed in terms of the long-term fitness We show that it is in fact not possible, as there is no way to define a single fitness parameter (or reproductive rate) that defines the fate of an organism even in a constant environment. This is due to fluctuations in the population averaged division rate. As a result, various details of the single-cell dynamics affect the fate of a new mutant independently from how they affect the long-term growth rate of the mutant population. In particular, we show that in the case of neutral mutations, variability in generation times increases the rate of genetic drift, and in the case of beneficial mutations, variability decreases its fixation probability. Furthermore, we explain the source of the persistent division rate fluctuations and provide analytic solutions for the fixation probability as a multispecies generalization of the Euler-Lotka equation.

The surge in interest in individual differences has coincided with the latest replication crisis centered around brain-wide association studies of brain-behavior correlations. Yet the reliability of the measures we use in cognitive neuroscience, a crucial component of this brain-behavior relationship, is often assumed but not directly tested. Here, we evaluate the reliability of different cognitive tasks on a large dataset of over 250 participants, who each completed a multi-day task battery. We show how reliability improves as a function of number of trials, and describe the convergence of the reliability curves for the different tasks, allowing us to score tasks according to their suitability for studies of individual differences. To improve the accessibility of these findings, we designed a simple web-based tool that implements this function to calculate the convergence factor and predict the expected reliability for any given number of trials and participants, even based on limited pilot data.Competing Interest StatementThe authors have declared no competing interest.

The statistics of noise emitted by ultrathin crumpled sheets is measured while they exhibit logarithmic relaxations under load. We find that the logarithmic relaxation advanced via a series of discrete, audible, micromechanical events that are log-Poisson distributed (i.e., the process becomes a Poisson process when time stamps are replaced by their logarithms). The analysis places constraints on the possible mechanisms underlying the glasslike slow relaxation and memory retention in these systems.

How cells regulate their cell cycles is a central question for cell biology. Models of cell size homeostasis have been proposed for bacteria, archaea, yeast, plant, and mammalian cells. New experiments bring forth high volumes of data suitable for testing existing models of cell size regulation and proposing new mechanisms. In this paper, we use conditional independence tests in conjunction with data of cell size at key cell cycle events (birth, initiation of DNA replication, and constriction) in the model bacterium Escherichia coli to select between the competing cell cycle models. We find that in all growth conditions that we study, the division event is controlled by the onset of constriction at midcell. In slow growth, we corroborate a model where replicationrelated processes control the onset of constriction at midcell. In faster growth, we find that the onset of constriction is affected by additional cues beyond DNA replication. Finally, we also find evidence for the presence of additional cues triggering initiations of DNA replication apart from the conventional notion where the mother cells solely determine the initiation event in the daughter cells via an adder per origin model. The use of conditional independence tests is a different approach in the context of understanding cell cycle regulation and it can be used in future studies to further explore the causal links between cell events.

Dynamic instability - the growth, catastrophe, and shrinkage of quasi-one-dimensional filaments - has been observed in multiple biopolymers. Scientists have long understood the catastrophic cessation of growth and subsequent depolymerization as arising from the interplay of hydrolysis and polymerization at the tip of the polymer. Here we show that for a broad class of catastrophe models, the expected catastrophe time distribution is exponential. We show that the distribution shape is insensitive to noise, but that depletion of monomers from a finite pool can dramatically change the distribution shape by reducing the polymerization rate. We derive a form for this finite-pool catastrophe time distribution and show that finite-pool effects can be important even when the depletion of monomers does not greatly alter the polymerization rate.

Adaptation dynamics on fitness landscapes is often studied theoretically in the strong-selection, weak-mutation regime. However, in a large population, multiple beneficial mutants can emerge before any of them fixes in the population. Competition between mutants is known as clonal interference, and while it is known to slow down the rate of adaptation (when compared to the strong-selection, weak-mutation model with the same parameters), how it affects the shape of long-term fitness trajectories in the presence of epistasis is an open question. Here, by considering how changes in fixation probabilities arising from weak clonal interference affect the dynamics of adaptation on fitness-parameterized landscapes, we find that the change in the shape of fitness trajectory arises only through changes in the supply of beneficial mutations (or equivalently, the beneficial mutation rate). Furthermore, a depletion of beneficial mutations as a population climbs up the fitness landscape can speed up the rescaled fitness trajectory (where adaptation speed is measured relative to its value at the start of the experiment), while an enhancement of the beneficial mutation rate does the opposite of slowing it down. Our findings suggest that by carrying out evolution experiments in both regimes (with and without clonal interference), one could potentially distinguish the different sources of macroscopic epistasis (fitness effect of mutations vs change in fraction of beneficial mutations).

The observation that phenotypic variability is ubiquitous in isogenic populations has led to a multitude of experimental and theoretical studies seeking to probe the causes and consequences of this variability. Whether it be in the context of antibiotic treatments or exponential growth in constant environments, non-genetic variability has significant effects on population dynamics. Here, we review research that elucidates the relationship between cell-to-cell variability and population dynamics. After summarizing the relevant experimental observations, we discuss models of bet-hedging and phenotypic switching. In the context of these models, we discuss how switching between phenotypes at the single-cell level can help populations survive in uncertain environments. Next, we review more fine-grained models of phenotypic variability where the relationship between single-cell growth rates, generation times and cell sizes is explicitly considered. Variability in these traits can have significant effects on the population dynamics, even in a constant environment. We show how these effects can be highly sensitive to the underlying model assumptions. We close by discussing a number of open questions, such as how environmental and intrinsic variability interact and what the role of non-genetic variability in evolutionary dynamics is.

Scientists have observed and studied diffusive waves in contexts as disparate as population genetics and cell signaling. Often, these waves are propagated by discrete entities or agents, such as individual cells in the case of cell signaling. For a broad class of diffusive waves, we characterize the transition between the collective propagation of diffusive waves, in which the wave speed is well described by continuum theory, and the propagation of diffusive waves by individual agents. We show that this transition depends heavily on the dimensionality of the system in which the wave propagates and that disordered systems yield dynamics largely consistent with lattice systems. In some system dimensionalities, the intuition that closely packed sources more accurately mimic a continuum can be grossly violated.

Mechanical rupture, or lysis, of the cytoplasmic membrane is a common cell death pathway in bacteria occurring in response to β-lactam antibiotics. A better understanding of the cellular design principles governing the susceptibility and response of individual cells to lysis could indicate methods of potentiating β-lactam antibiotics and clarify relevant aspects of cellular physiology. Here, we take a single-cell approach to bacterial cell lysis to examine three cellular featuresturgor pressure, mechanosensitive channels, and cell shape changesthat are expected to modulate lysis. We develop a mechanical model of bacterial cell lysis and experimentally analyze the dynamics of lysis in hundreds of single Escherichia coli cells. We find that turgor pressure is the only factor, of these three cellular features, which robustly modulates lysis. We show that mechanosensitive channels do not modulate lysis due to insufficiently fast solute outflow, and that cell shape changes result in more severe cellular lesions but do not influence the dynamics of lysis. These results inform a single-cell view of bacterial cell lysis and underscore approaches of combatting antibiotic tolerance to β-lactams aimed at targeting cellular turgor.

Gene expression is a stochastic process. Despite the increase of protein numbers in growing cells, the protein concentrations are often found to be confined within small ranges throughout the cell cycle. Generally, the noise in protein concentration can be decomposed into an intrinsic and an extrinsic component, where the former vanishes for high expression levels. Considering the time trajectory of protein concentration as a random walker in the concentration space, an effective restoring force (with a corresponding "spring constant") must exist to prevent the divergence of concentration due to random fluctuations. In this work, we prove that the magnitude of the effective spring constant is directly related to the fraction of intrinsic noise in the total protein concentration noise. We show that one can infer the magnitude of intrinsic, extrinsic, and measurement noises of gene expression solely based on time-resolved data of protein concentration, without any a priori knowledge of the underlying gene expression dynamics. We apply this method to experimental data of single-cell bacterial gene expression. The results allow us to estimate the average copy numbers and the translation burst parameters of the studied proteins.

In biological contexts as diverse as development, apoptosis, and synthetic microbial consortia, collections of cells or subcellular components have been shown to overcome the slow signaling speed of simple diffusion by utilizing diffusive relays, in which the presence of one type of diffusible signaling molecule triggers participation in the emission of the same type of molecule. This collective effect gives rise to fast-traveling diffusive waves. Here, in the context of cell signaling, we show that system dimensionality the shape of the extracellular medium and the distribution of cells within it can dramatically affect the wave dynamics, but that these dynamics are insensitive to details of cellular activation. As an example, we show that neutrophil swarming experiments exhibit dynamical signatures consistent with the proposed signaling motif. We further show that cell signaling relays generate much steeper concentration profiles than does simple diffusion, which may facilitate neutrophil chemotaxis.

Selection of mutants in a microbial population depends on multiple cellular traits. In serial-dilution evolution experiments, three key traits are the lag time when transitioning from starvation to growth, the exponential growth rate, and the yield (number of cells per unit resource). Here, we investigate how these traits evolve in laboratory evolution experiments using a minimal model of population dynamics, where the only interaction between cells is competition for a single limiting resource. We find that the fixation probability of a beneficial mutation depends on a linear combination of its growth rate and lag time relative to its immediate ancestor, even under clonal interference. The relative selective pressure on growth rate and lag time is set by the dilution factor; a larger dilution factor favors the adaptation of growth rate over the adaptation of lag time. The model shows that yield, however, is under no direct selection. We also show how the adaptation speeds of growth and lag depend on experimental parameters and the underlying supply of mutations. Finally, we investigate the evolution of covariation between these traits across populations, which reveals that the population growth rate and lag time can evolve a nonzero correlation even if mutations have uncorrelated effects on the two traits. Altogether these results provide useful guidance to future experiments on microbial evolution.

In isogenic microbial populations, phenotypic variability is generated by a combination of stochastic mechanisms, such as gene expression, and deterministic factors, such as asymmetric segregation of cell volume. Here we address the question: how does phenotypic variability of a microbial population affect its fitness? While this question has previously been studied for exponentially growing populations, the situation when the population size is kept fixed has received much less attention, despite its relevance to many natural scenarios. We show that the outcome of competition between multiple microbial species can be determined from the distribution of phenotypes in the culture using a generalization of the well-known Euler-Lotka equation, which relates the steady-state distribution of phenotypes to the population growth rate. We derive a generalization of the Euler-Lotka equation for finite cultures, which relates the distribution of phenotypes among cells in the culture to the exponential growth rate. Our analysis reveals that in order to predict fitness from phenotypes, it is important to understand how distributions of phenotypes obtained from different subsets of the genealogical history of a population are related. To this end, we derive a mapping between the various ways of sampling phenotypes in a finite population and show how to obtain the equivalent distributions from an exponentially growing culture. Finally, we use this mapping to show that species with higher growth rates in exponential growth conditions will have a competitive advantage in the finite culture.

The generalized central limit theorem is a remarkable generalization of the central limit theorem, showing that the sum of a large number of independent, identically-distributed (i.i.d) random variables with infinite variance may converge under appropriate scaling to a distribution belonging to a special family known as Levy stable distributions. Similarly, the maximum of i.i.d. variables may converge to a distribution belonging to one of three universality classes (Gumbel, Weibull and Frechet). Here, we rederive these known results following a mathematically non-rigorous yet highly transparent renormalization-group-inspired approach that captures both of these universal results following a nearly identical procedure.

We study the spread of a quantum-mechanical wave packet in a noisy environment, modeled using a tight-binding Hamiltonian. Despite the coherent dynamics, the fluctuating environment may give rise to diffusive behavior. When correlations between different level-crossing events can be neglected, we use the solution of the Landau-Zener problem to find how the diffusion constant depends on the noise. We also show that when an electric field or external disordered potential is applied to the system, the diffusion constant is suppressed with no drift term arising. The results are relevant to various quantum systems, including exciton diffusion in photosynthesis and electronic transport in solid-state physics.

Membrane lysis, or rupture, is a cell death pathway in bacteria frequently caused by cell wall-targeting antibiotics. Although previous studies have clarified the biochemical mechanisms of antibiotic action, a physical understanding of the processes leading to lysis remains lacking. Here, we analyze the dynamics of membrane bulging and lysis in Escherichia coli, in which the formation of an initial, partially subtended spherical bulge (\u201cbulging\u201d) after cell wall digestion occurs on a characteristic timescale of 1 s and the growth of the bulge (\u201cswelling\u201d) occurs on a slower characteristic timescale of 100 s. We show that bulging can be energetically favorable due to the relaxation of the entropic and stretching energies of the inner membrane, cell wall, and outer membrane and that the experimentally observed timescales are consistent with model predictions. We then show that swelling is mediated by the enlargement of wall defects, after which cell lysis is consistent with both the inner and outer membranes exceeding characteristic estimates of the yield areal strains of biological membranes. These results contrast biological membrane physics and the physics of thin, rigid shells. They also have implications for cellular morphogenesis and antibiotic discovery across different species of bacteria.

It is well known that the contribution of harmonic phonons to the thermal conductivity of 1D systems diverges with the harmonic chain length L. Within various one-dimensional models containing disorder it was shown that the thermal conductivity scales as under certain boundary conditions. Here we show that when the chain is weakly coupled to the heat reservoirs and there is strong disorder this scaling can be violated. We find a weaker power-law dependence on L, and show that for sufficiently strong disorder the thermal conductivity ceases to be anomalous - it does not depend on L and hence obeys Fourier's law. This is despite both density of states and the diverging localization length scaling anomalously. Surprisingly, in this strong disorder regime two anomalously scaling quantities cancel each other to recover Fourier's law of heat transport.

Discriminating between correct and incorrect substrates is a core process in biology, but how is energy apportioned between the conflicting demands of accuracy (μ), speed (σ), and total entropy production rate (P)? Previous studies have focused on biochemical networks with simple structure or relied on simplifying kinetic assumptions. Here, we use the linear framework for timescale separation to analytically examine steady-state probabilities away from thermodynamic equilibrium for networks of arbitrary complexity. We also introduce a method of scaling parameters that is inspired by Hopfield's treatment of kinetic proofreading. Scaling allows asymptotic exploration of high-dimensional parameter spaces. We identify in this way a broad class of complex networks and scalings for which the quantity σln(μ)/P remains asymptotically finite whenever accuracy improves from equilibrium, so that μeq/μ→0. Scalings exist, however, even for Hopfield's original network, for which σln(μ)/P is asymptotically infinite, illustrating the parametric complexity. Outside the asymptotic regime, numerical calculations suggest that, under more restrictive parametric assumptions, networks satisfy the bound, σln(μ/μeq)/P

For many decades, the wedding of quantitative data with mathematical modeling has been fruitful, leading to important biological insights. Here, we review some of the ongoing efforts to gain insights into problems in microbiology and, in particular, cell-cycle progression and its regulation through observation and quantitative analysis of the natural fluctuations in the system. We first illustrate this idea by reviewing a classic example in microbiology the LuriaDelbrück experiment and discussing how, in that case, useful information was obtained by looking beyond the mean outcome of the experiment, but instead paying attention to the variability between replicates of the experiment. We then switch gears to the contemporary problem of cell cycle progression and discuss in more detail how insights into cell size regulation and, when relevant, coupling between the cell cycle and the circadian clock, can be gained by studying the natural fluctuations in the system and their statistical properties. We end with a more general discussion of how (in this context) the correct level of phenomenological model should be chosen, as well as some of the pitfalls associated with this type of analysis. Throughout this review the emphasis is not on providing details of the experimental setups or technical details of the models used, but rather, in fleshing out the conceptual structure of this particular approach to the problem. For this reason, we choose to illustrate the framework on a rather broad range of problems, and on organisms from all domains of life, to emphasize the commonality of the ideas and analysis used (as well as their differences).

In model bacteria, such as E. coli and B. subtilis, regulation of cell-cycle progression and cellular organization achieves consistency in cell size, replication dynamics, and chromosome positioning [13]. Mycobacteria elongate and divide asymmetrically, giving rise to significant variation in cell size and elongation rate among closely related cells [4, 5]. Given the physical asymmetry of mycobacteria, the models that describe coordination of cellular organization and cell-cycle progression in model bacteria are not directly translatable [1, 2, 68]. Here, we used time-lapse microscopy and fluorescent reporters of DNA replication and chromosome positioning to examine the coordination of growth, division, and chromosome dynamics at a single-cell level in Mycobacterium smegmatis (M. smegmatis) and Mycobacterium bovis Bacillus Calmette-Guérin (BCG). By analyzing chromosome and replisome localization, we demonstrated that chromosome positioning is asymmetric and proportional to cell size. Furthermore, we found that cellular asymmetry is maintained throughout the cell cycle and is not established at division. Using measurements and stochastic modeling of mycobacterial cell size and cell-cycle timing in both slow and fast growth conditions, we found that well-studied models of cell-size control are insufficient to explain the mycobacterial cell cycle. Instead, we showed that mycobacterial cell-cycle progression is regulated by an unprecedented mechanism involving parallel adders (i.e., constant growth increments) that start at replication initiation. Together, these adders enable mycobacterial populations to regulate cell size, growth, and heterogeneity in the face of varying environments. Logsdon et al. study cell-cycle dynamics of asymmetrically growing mycobacteria and show how cell cycle, chromosome organization, and division are coordinated in single cells. A parallel adder model, where cells add a constant length between initiations and initiation to division, describes cell-size control in M. smegmatis and BCG.

Establishing a quantitative connection between the population growth rate and the generation times of single cells is a prerequisite for understanding evolutionary dynamics of microbes. However, existing theories fail to account for the experimentally observed correlations between mother-daughter generation times that are unavoidable when cell size is controlled for, which is essentially always the case. Here, we study population-level growth in the presence of cell size control and corroborate our theory using experimental measurements of single-cell growth rates. We derive a closed formula for the population growth rate and demonstrate that it only depends on the single-cell growth rate variability, not other sources of stochasticity. Our work provides an evolutionary rationale for the narrow growth rate distributions often observed in nature: when single-cell growth rates are less variable but have a fixed mean, the population will exhibit an enhanced population growth rate as long as the correlations between the mother and daughter cells growth rates are not too strong. The effects of single-cell variability on the population growth are studied, and it is found that the only stochastic term of consequence is the growth rate variability, which typically decreases the population growth rate.

At low temperatures the dynamical degrees of freedom in amorphous solids are tunneling two-level systems (TLSs). Concentrating on these degrees of freedom, and taking into account disorder and TLS-TLS interactions, we obtain a "TLS glass," described by the random-field Ising model with random 1/r3 interactions. In this paper we perform a self-consistent mean-field calculation, previously used to study the electron-glass (EG) model [A. Amir, Phys. Rev. B 77, 165207 (2008)PRBMDO1098-012110.1103/PhysRevB.77.165207]. Similarly to the electron glass, we find a 1λ distribution of relaxation rates λ, leading to logarithmic slow relaxation. However, with increased interactions the EG model shows slower dynamics whereas the TLS-glass model shows faster dynamics. This suggests that given system-specific properties, glass dynamics can be slowed down or sped up by the interactions.

Photonic crystals with slowly varying lattice period are responsible for broadband light reflectance in many biological contexts, ranging from the coatings of shiny beetles to the eye shine of butterflies. We utilize a quantum scattering analogy to obtain an approximate but closed-form expression for the reflectance of these adiabatically chirped photonic crystals (ACPCs). This expression allows us to devise a method for designing ACPCs with tailored reflectance spectra and obtain an estimate for the minimal number of layers required to exceed a given reflectance threshold. Comparison against the number of layers found in ACPCs throughout nature gives a quantitative measure of the optimality of chirped biological reflectors. Together, these results elucidate the design of chirped reflectors in nature and their possible application to future optical technologies.

We consider a class of biologically motivated stochastic processes in which a unicellular organism divides its resources (volume or damaged proteins, in particular) symmetrically or asymmetrically between its progeny. Assuming the final amount of the resource is controlled by a growth policy and subject to additive and multiplicative noise, we derive the recursive integral equation describing the evolution of the resource distribution over subsequent generations and use it to study the properties of stable resource distributions. We find conditions under which a unique stable resource distribution exists and calculate its moments for the class of affine linear growth policies. Moreover, we apply an asymptotic analysis to elucidate the conditions under which the stable distribution (when it exists) has a power-law tail. Finally, we use the results of this asymptotic analysis along with the moment equations to draw a stability phase diagram for the system that reveals the counterintuitive result that asymmetry serves to increase stability while at the same time widening the stable distribution. We also briefly discuss how cells can divide damaged proteins asymmetrically between their progeny as a form of damage control. In the appendixes, motivated by the asymmetric division of cell volume in Saccharomyces cerevisiae, we extend our results to the case wherein mother and daughter cells follow different growth policies.

To maintain a constant cell size, dividing cells have to coordinate cell-cycle events with cell growth. This coordination has long been supposed to rely on the existence of size thresholds determining cell-cycle progression [1]. In budding yeast, size is controlled at the G1/S transition [2]. In agreement with this hypothesis, the size at birth influences the time spent in G1: smaller cells have a longer G1 period [3]. Nevertheless, even though cells born smaller have a longer G1, the compensation is imperfect and they still bud at smaller cell sizes. In bacteria, several recent studies have shown that the incremental model of size control, in which size is controlled by addition of a constant volume (in contrast to a size threshold), is able to quantitatively explain the experimental data on four different bacterial species [4-7]. Here, we report on experimental results for the budding yeast Saccharomyces cerevisiae, finding, surprisingly, that cell size control in this organism is very well described by the incremental model, suggesting a common strategy for cell size control with bacteria. Additionally, we argue that for S. cerevisiae the "volume increment" is not added from birth to division, but rather between two budding events.

Characterizing the frequency-dependent response of amorphous systems and glasses can provide important insights into their physics. Here, we study the response of an electron glass, where Coulomb interactions are important and have previously been shown to significantly modify the conductance and lead to memory effects and aging. We propose a model which allows us to take the interactions into account in a self-consistent way, and explore the frequency-dependent conduction at all frequencies. At low frequencies conduction occurs on the percolation backbone, and the model captures the variable-range-hopping behavior. At high frequencies conduction is dominated by localized clusters. Despite the difference in physical mechanisms at low and high frequency, we are able to approximately scale all numerical data onto a single curve, using two parameters: the DC conduction and the DC dielectric constant. The behavior follows the universal scaling that is experimentally observed for a large class of amorphous solids.

Cell walls define a cell's shape in bacteria. The walls are rigid to resist large internal pressures, but remarkably plastic to adapt to a wide range of external forces and geometric constraints. Currently, it is unknown how bacteria maintain their shape. In this paper, we develop experimental and theoretical approaches and show that mechanical stresses regulate bacterial cell wall growth. By applying a precisely controllable hydrodynamic force to growing rod-shaped Escherichia coli and Bacillus subtilis cells, we demonstrate that the cells can exhibit two fundamentally different modes of deformation. The cells behave like elastic rods when subjected to transient forces, but deform plastically when significant cell wall synthesis occurs while the force is applied. The deformed cells always recover their shape. The experimental results are in quantitative agreementwith the predictions of the theory of dislocation-mediated growth. In particular, we find that a single dimensionless parameter, which depends on a combination of independently measured physical properties of the cell, can describe the cell's responses under various experimental conditions. These findings provide insight into how living cells robustly maintain their shape under varying physical environments.

Interference is the source of some of the spectacular colors of animals and plants in nature. In some of these systems, the physical structure consists of an ordered array of layers with alternating high and low refractive indices. This periodicity leads to an optical band structure that is analogous to the electronic band structure encountered in semiconductor physics: specific bands of wavelengths (the stop bands) are perfectly reflected. Here, we present a minimal model for optical band structure in a periodic multilayer structure and solve it using recursion relations. The stop bands emerge in the limit of an infinite number of layers by finding the fixed point of the recursion. We compare to experimental data for various beetles, whose optical structure resembles the proposed model. Thus, using only the phenomenon of interference and the idea of recursion, we are able to elucidate the concept of band structure in the context of the experimentally observed high reflectance and iridescent appearance of structurally colored beetles.

Slow relaxation occurs in many physical and biological systems. "Creep" is an example from everyday life. When stretching a rubber band, for example, the recovery to its equilibrium length is not, as one might think, exponential: The relaxation is slow, in many cases logarithmic, and can still be observed after many hours. The form of the relaxation also depends on the duration of the stretching, the "waiting time." This ubiquitous phenomenon is called aging, and is abundant both in natural and technological applications. Here, we suggest a general mechanism for slow relaxations and aging, which predicts logarithmic relaxations, and a particular aging dependence on the waiting time. We demonstrate the generality of the approach by comparing our predictions to experimental data on a diverse range of physical phenomena, from conductance in granular metals to disordered insulators and dirty semiconductors, to the low temperature dielectric properties of glasses.

We study the spectral properties of a class of random matrices where the matrix elements depend exponentially on the distance between uniformly and randomly distributed points. This model arises naturally in various physical contexts, such as the diffusion of particles, slow relaxations in glasses, and scalar phonon localization. Using a combination of a renormalization group procedure and a direct moment calculation, we find the eigenvalue distribution density (i.e., the spectrum), for low densities, and the localization properties of the eigenmodes, for arbitrary dimension. Finally, we discuss the physical implications of the results.

Glassy systems are ubiquitous in nature. They are characterized by slow relaxations to equilibrium without a typical time scale, aging, and memory effects. Understanding this has been a long-standing problem in physics. We study the aging of the electron glass, a system showing remarkable slow relaxations of the conductance. We find that the appropriate broad distribution of relaxation rates leads to a universal relaxation of the form log(1+tw/t) for the common aging protocol, where tw is the length of time the perturbation driving the system out of equilibrium was on, and t the time of measurement. These results agree well with several experiments performed on different glassy systems, and examining different physical observables, for times ranging from seconds to several hours. The suggested theoretical framework appears to offer a paradigm for aging in a broad class of glassy materials.