Antibiotic responses in bacteria are highly dynamic and heterogeneous, with sudden exposure of bacterial colonies to high drug doses resulting in the coexistence of recovered and arrested cells. The dynamics of the response is determined by regulatory circuits controlling the expression of resistance genes, which are in turn modulated by the drug's action on cell growth and metabolism. Despite advances in understanding gene regulation at the molecular level, we still lack a framework to describe how feedback mechanisms resulting from the interdependence between expression of resistance and cell metabolism can amplify naturally occurring noise and create heterogeneity at the population level. To understand how this interplay affects cell survival upon exposure, we constructed a mathematical model of the dynamics of antibiotic responses that links metabolism and regulation of gene expression, based on the tetracycline resistancetetoperon inE. coli. We use this model to interpret measurements of growth and expression of resistance in microfluidic experiments, both in single cells and in biofilms. We also implemented a stochastic model of the drug response, to show that exposure to high drug levels results in large variations of recovery times and heterogeneity at the population level. We show that stochasticity is important to determine how nutrient quality affects cell survival during exposure to high drug concentrations. A quantitative description of how microbes respond to antibiotics in dynamical environments is crucial to understand population-level behaviors such as biofilms and pathogenesis.
Understanding the structural and functional development of human-induced pluripotent stem-cell-derived cardiomyocytes (hiPSC-CMs) is essential to engineering cardiac tissue that enables pharmaceutical testing, modeling diseases, and designing therapies. Here we use a method not commonly applied to biological materials, small angle x-ray scattering, to characterize the structural development of hiPSC-CMs within three-dimensional engineered tissues during their preliminary stages of maturation. An x-ray scattering experimental method enables the reliable characterization of the cardiomyocyte myofilament spacing with maturation time. The myofilament lattice spacing monotonically decreases as the tissue matures from its initial post-seeding state over the span of 10 days. Visualization of the spacing at a grid of positions in the tissue provides an approach to characterizing the maturation and organization of cardiomyocyte myofilaments and has the potential to help elucidate mechanisms of pathophysiology, and disease progression, thereby stimulating new biological hypotheses in stem cell engineering.
Dopaminergic neurons are specialized cells in the substantia nigra, tasked with dopamine secretion. This secretion relies on intracellular calcium signaling coupled to neuronal electrical activity. These neurons are known to display spontaneous calcium oscillationsin-vitroandin-vivo, even in synaptic isolation, controlling the basal dopamine levels. Here we outline a kinetic model for the ion exchange across the neuronal plasma membrane. Crucially, we relax the assumption of constant, cytoplasmic sodium and potassium concentration. We show that sodium-potassium dynamics are strongly coupled to calcium dynamics and are essential for the robustness of spontaneous firing frequency. The model predicts several regimes of electrical activity, including tonic and 'burst' oscillations, and predicts the switch between those in response to perturbations. 'Bursting' correlates with increased calcium amplitudes, while maintaining constant average, allowing for a vast change in the calcium signal responsible for dopamine secretion. All the above traits provide the flexibility to create rich action potential dynamics that are crucial for cellular function.
Computational modeling of cancer can help unveil dynamics and interactions that are hard to replicate experimentally. Thanks to the advancement in cancer databases and data analysis technologies, these models have become more robust than ever. There are many mathematical models which investigate cancer through different approaches, from sub-cellular to tissue scale, and from treatment to diagnostic points of view. In this study, we lay out a step-by-step methodology for a data-driven mechanistic model of the tumor microenvironment. We discuss data acquisition strategies, data preparation, parameter estimation, and sensitivity analysis techniques. Furthermore, we propose a possible approach to extend mechanistic ordinary differential equation models to PDE models coupled with mechanical growth. The workflow discussed in this article can help understand the complex temporal and spatial interactions between cells and cytokines in the tumor microenvironment and their effect on tumor growth.
A fundamental question in complex systems is how to relate interactions between individual components ('microscopic description') to the global properties of the system ('macroscopic description'). Furthermore, it is unclear whether such a macroscopic description exists and if such a description can capture large-scale properties. Here, we address the validity of a macroscopic description of a complex biological system using the collective motion of desert locusts as a canonical example. One of the world's most devastating insect plagues begins when flightless juvenile locusts form 'marching bands'. These bands display remarkable coordinated motion, moving through semiarid habitats in search of food. We investigated how well macroscopic physical models can describe the flow of locusts within a band. For this, we filmed locusts within marching bands during an outbreak in Kenya and automatically tracked all individuals passing through the camera frame. We first analyzed the spatial topology of nearest neighbors and found individuals to be isotropically distributed. Despite this apparent randomness, a local order was observed in regions of high density in the radial distribution function, akin to an ordered fluid. Furthermore, reconstructing individual locust trajectories revealed a highly aligned movement, consistent with the one-dimensional version of the Toner-Tu equations, a generalization of the Navier-Stokes equations for fluids, used to describe the equivalent macroscopic fluid properties of active particles. Using this effective Toner-Tu equation, which relates the gradient of the pressure to the acceleration, we show that the effective 'pressure' of locusts increases as a linear function of density in segments with the highest polarization (for which the one-dimensional approximation is most appropriate). Our study thus demonstrates an effective hydrodynamic description of flow dynamics in plague locust swarms.
In recentin vitroexperiments on co-culture between breast tumour spheroids and activated immune cells, it was observed that the introduction of the stress hormone cortisol resulted in a decreased immune cell infiltration into the spheroids. Moreover, the presence of cortisol deregulated the normal levels of the pro- and anti-inflammatory cytokines IFN-γand IL-10. We present an individual-based model to explore the interaction dynamics between tumour and immune cells under psychological stress conditions. With our model, we explore the processes underlying the emergence of different levels of immune infiltration, with particular focus on the biological mechanisms regulated by IFN-γand IL-10. The set-up of numerical simulations is defined to mimic the scenarios considered in the experimental study. Similarly to the experimental quantitative analysis, we compute a score that quantifies the level of immune cell infiltration into the tumour. The results of numerical simulations indicate that the motility of immune cells, their capability to infiltrate through tumour cells, their growth rate and the interplay between these cell parameters can affect the level of immune cell infiltration in different ways. Ultimately, numerical simulations of this model support a deeper understanding of the impact of biological stress-induced mechanisms on immune infiltration.
Proteins populate a manifold in the high-dimensional sequence space whose geometrical structure guides their natural evolution. Leveraging recently-developed structure prediction tools based on transformer models, we first examine the protein sequence landscape as defined by an effective energy that is a proxy of sequence foldability. This landscape shares characteristics with optimization challenges encountered in machine learning and constraint satisfaction problems. Our analysis reveals that natural proteins predominantly reside in wide, flat minima within this energy landscape. To investigate further, we employ statistical mechanics algorithms specifically designed to explore regions with high local entropy in relatively flat landscapes. Our findings indicate that these specialized algorithms can identify valleys with higher entropy compared to those found using traditional methods such as Monte Carlo Markov Chains. In a proof-of-concept case, we find that these highly entropic minima exhibit significant similarities to natural sequences, especially in critical key sites and local entropy. Additionally, evaluations through Molecular Dynamics suggests that the stability of these sequences closely resembles that of natural proteins. Our tool combines advancements in machine learning and statistical physics, providing new insights into the exploration of sequence landscapes where wide, flat minima coexist alongside a majority of narrower minima.
Fungi expand in space and time to form complex multicellular communities. The mechanisms by which they do so can vary dramatically and determine the life-history and dispersal traits of expanding populations. These traits influence deterministic and stochastic components of evolution, resulting in complex eco-evolutionary dynamics during colony expansion. We perform experiments on budding yeast strains genetically engineered to display rough-surface and smooth-surface phenotypes in colony-like structures called 'mats'. Previously, it was shown that the rough-surface strain has a competitive advantage over the smooth-surface strain when grown on semi-solid media. We experimentally observe the emergence and expansion of segments with a distinct smooth-surface phenotype during rough-surface mat development. We propose a trade-off between dispersal and local carrying capacity to explain the relative fitness of these two phenotypes. Using a modified stepping-stone model, we demonstrate that this trade-off gives the high-dispersing, rough-surface phenotype a competitive advantage from standing variation, but that it inhibits this phenotype's ability to invade a resident smooth-surface population via mutation. However, the trade-off improves the ability of the smooth-surface phenotype to invade in rough-surface mats, replicating the frequent emergence of smooth-surface segments in experiments. Together, these computational and experimental findings advance our understanding of the complex eco-evolutionary dynamics of fungal mat expansion.