BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Bartlomiej Waclaw appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
The problem of antibiotic resistance poses challenges across many disciplines. One such challenge is to understand the fundamental science of how antibiotics work, and how resistance to them can emerge. This is an area where physicists can make important contributions. Here, we highlight cases where this is already happening, and suggest directions for further physics involvement in antimicrobial research.
Stochastic phenotype switching has been suggested to play a beneficial role in microbial populations by leading to the division of labour among cells, or ensuring that at least some of the population survives an unexpected change in environmental conditions. Here we use a computational model to investigate an alternative possible function of stochastic phenotype switching - as a way to adapt more quickly even in a static environment. We show that when a genetic mutation causes a population to become less fit, switching to an alternative phenotype with higher fitness (growth rate) may give the population enough time to develop compensatory mutations that increase the fitness again. The possibility of switching phenotypes can reduce the time to adaptation by orders of magnitude if the "fitness valley" caused by the deleterious mutation is deep enough. Our work has important implications for the emergence of antibiotic-resistant bacteria. In line with recent experimental findings we hypothesise that switching to a slower growing but less sensitive phenotype helps bacteria to develop resistance by exploring a larger set of beneficial mutations while avoiding deleterious ones.
We present a simple non-equilibrium model of mass condensation with Lennard-Jones interactions between particles and the substrate. We show that when some number of particles is deposited onto the surface and the system is left to equilibrate, particles condense into an island if the density of particles becomes higher than some critical density. We illustrate this with numerically obtained phase diagrams for three-dimensional systems. We also solve a two-dimensional counterpart of this model analytically and show that not only the phase diagram but also the shape of the cross-sections of three-dimensional condensates qualitatively matches the two-dimensional predictions. Lastly, we show that when particles are being deposited with a constant rate, the system has two phases: a single condensate for low deposition rates, and multiple condensates for fast deposition. The behaviour of our model is thus similar to that of thin film growth processes, and in particular to Stranski-Krastanov growth.
Most cancers in humans are large, measuring centimeters in diameter, composed of many billions of cells. An equivalent mass of normal cells would be highly heterogeneous as a result of the mutations that occur during each cell division. What is remarkable about cancers is their homogeneity - virtually every neoplastic cell within a large cancer contains the same core set of genetic alterations, with heterogeneity confined to mutations that have emerged after the last clonal expansions. How such clones expand within the spatially-constrained three dimensional architecture of a tumor, and come to dominate a large, pre-existing lesion, has never been explained. We here describe a model for tumor evolution that shows how short-range migration and cell turnover can account for rapid cell mixing inside the tumor. With it, we show that even a small selective advantage of a single cell within a large tumor allows the descendants of that cell to replace the precursor mass in a clinically relevant time frame. We also demonstrate that the same mechanisms can be responsible for the rapid onset of resistance to chemotherapy. Our model not only provides novel insights into spatial and temporal aspects of tumor growth but also suggests that targeting short range cellular migratory activity could have dramatic effects on tumor growth rates.
The universality of many pathways of core metabolism suggests a strong role for evolutionary selection, but it remains unclear whether existing pathways have been selected from a large or small set of biochemical possibilities. To address this question, we construct "in silico" all possible biochemically feasible alternatives to the trunk pathway of glycolysis and gluconeogenesis, one of the most highly conserved pathways in metabolism. We show that, even though a large number of alternative pathways exist, the alternatives carry lower flux than the real pathway under typical physiological conditions. Alternative pathways that could potentially carry higher flux often lead to infeasible intermediate metabolite concentrations. We also find that if physiological conditions were different, different pathways could outperform those found in nature. Our results demonstrate how the rules of biochemistry restrict the alternatives that are open to evolution, and suggest that the existing trunk pathway of glycolysis and gluconeogenesis represents a maximal flux solution.
We introduce a model of a totally asymmetric simple exclusion process (TASEP) on a tree network where the aggregate hopping rate is constant from level to level. With this choice for hopping rates the model shows the same phase diagram as the one-dimensional case. The potential applications of our model are in the area of distribution networks; where a single large source supplies material to a large number of small sinks via a hierarchical network. We show that mean field theory (MFT) for our model is identical to that of the one-dimensional TASEP and that this mean field theory is exact for the TASEP on a tree in the limit of large branching ratio, $b$(or equivalently large coordination number). We then present an exact solution for the two level tree (or star network) that allows the computation of any correlation function and confirm how mean field results are recovered as $b\rightarrow\infty$. As an example we compute the steady-state current as a function of branching ratio. We present simulation results that confirm these results and indicate that the convergence to MFT with large branching ratio is quite rapid.
We study a far-from-equilibrium system of interacting particles, hopping between sites of a 1d lattice with a rate which increases with the number of particles at interacting sites. We find that clusters of particles, which initially spontaneously form in the system, begin to move at increasing speed as they gain particles. Ultimately, they produce a moving condensate which comprises a finite fraction of the mass in the system. We show that, in contrast with previously studied models of condensation, the relaxation time to steady state decreases as an inverse power of ln L with system size L and that condensation is instantenous for L-->infinity.
Drug gradients are believed to play an important role in the evolution of bacteria resistant to antibiotics and tumors resistant to anti-cancer drugs. We use a statistical physics model to study the evolution of a population of malignant cells exposed to drug gradients, where drug resistance emerges via a mutational pathway involving multiple mutations. We show that a non-uniform drug distribution has the potential to accelerate the emergence of resistance when the mutational pathway involves a long sequence of mutants with increasing resistance, but if the pathway is short or crosses a fitness valley, the evolution of resistance may actually be slowed down by drug gradients. These predictions can be verified experimentally, and may help to improve strategies for combatting the emergence of resistance.
Mixed positive and negative feedback loops are often found in biological systems which support oscillations. In this work we consider a prototype of such systems, which has been recently found at the core of many genetic circuits showing oscillatory behaviour. Our model consists of two interacting species A and B, where A activates not only its own production, but also that of its repressor B. While the self-activation of A leads already to a bistable unit, the coupling with a negative feedback loop via B makes the unit frustrated. In the deterministic limit of infinitely many molecules, such a bistable frustrated unit is known to show excitable and oscillatory dynamics, depending on the maximum production rate of A which acts as a control parameter. We study this model in its fully stochastic version and we find oscillations even for parameters which in the deterministic limit are deeply in the fixed-point regime. The deeper we go into this regime, the more irregular these oscillations are, becoming finally random excitations whenever fluctuations allow the system to overcome the barrier for a large excursion in phase space. The fluctuations can no longer be fully treated as a perturbation. The smaller the system size (the number of molecules), the more frequent are these excitations. Therefore, stochasticity caused by demographic noise makes this unit even more flexible with respect to its oscillatory behaviour.
The quasispecies model of biological evolution for asexual organisms such as bacteria and viruses has attracted considerable attention of biological physicists. Many variants of the model have been proposed and subsequently solved using the methods of statistical physics. In this paper I will put forward important but largely overlooked relations between localization theory, random matrices, and the quasispecies model. These relations will help me to study the dynamics of this model. In particular, I will show that the distribution of times between evolutionary jumps in the genotype space follows a power law, in agreement with recent findings in the shell model - a simplified version of the quasispecies model.
Migration between different habitats is ubiquitous among biological populations. In this Letter, we study a simple quasispecies model for evolution in two different habitats, with different fitness landscapes, coupled through one-way migration. Our model applies to asexual, rapidly evolving organisms such as microbes. Our key finding is a dynamical phase transition at a critical value of the migration rate. The time to reach steady state diverges at this critical migration rate. Above the transition, the population is dominated by immigrants from the primary habitat. Below the transition, the genetic composition of the population is highly non-trivial, with multiple coexisting quasispecies which are not native to either habitat. Using results from localization theory, we show that the critical migration rate may be very small --- demonstrating that evolutionary outcomes can be very sensitive to even a small amount of migration.
We present an analytic method for calculating spectral densities of empirical covariance matrices for correlated data. In this approach the data is represented as a rectangular random matrix whose columns correspond to sampled states of the system. The method is applicable to a class of random matrices with radial measures including those with heavy (power-law) tails in the probability distribution. As an example we apply it to a multivariate Student distribution.
We discuss a practical method to determine the eigenvalue spectrum of the empirical correlation matrix. The method is based on the analysis of the behavior of a conformal map at a critical horizon which is defined as a border line of the physical Riemann sheet of this map. The map is a convenient representation of the Marčenko-Pastur equation.