Physical biology最新文献

Analyses of structural dynamics of biomolecules hold great promise to deepen the understanding of and ability to construct complex molecular systems. To this end, both experimental and computational means are available, such as fluorescence quenching experiments or molecular dynamics simulations, respectively. We argue that while seemingly disparate, both fields of study have to deal with the same type of data about the same underlying phenomenon of conformational switching. Two central challenges typically arise in both contexts: (i) the amount of obtained data is large, and (ii) it is often unknown how many distinct molecular states underlie these data. In this study, we build on the established idea of Markov state modeling and propose a generative, Bayesian nonparametric hidden Markov state model that addresses these challenges. Utilizing hierarchical Dirichlet processes, we treat different meta-stable molecule conformations as distinct Markov states, the number of which we then do not have to seta priori. In contrast to existing approaches to both experimental as well as simulation data that are based on the same idea, we leverage a mean-field variational inference approach, enabling scalable inference on large amounts of data. Furthermore, we specify the model also for the important case of angular data, which however proves to be computationally intractable. Addressing this issue, we propose a computationally tractable approximation to the angular model. We demonstrate the method on synthetic ground truth data and apply it to known benchmark problems as well as electrophysiological experimental data from a conformation-switching ion channel to highlight its practical utility.

Cellular populations assume an incredible variety of shapes ranging from circular molds to irregular tumors. While we understand many of the mechanisms responsible for these spatial patterns, little is known about how the shape of a population influences its ecology and evolution. Here, we investigate this relationship in the context of microbial colonies grown on hard agar plates. This a well-studied system that exhibits a transition from smooth circular disks to more irregular and rugged shapes as either the nutrient concentration or cellular motility is decreased. Starting from a mechanistic model of colony growth, we identify two dimensionless quantities that determine how morphology and genetic diversity of the population depend on the model parameters. Our simulations further reveal that population dynamics cannot be accurately described by the commonly-used surface growth models. Instead, one has to explicitly account for the emergent growth instabilities and demographic fluctuations. Overall, our work links together environmental conditions, colony morphology, and evolution. This link is essential for a rational design of concrete, biophysical perturbations to steer evolution in the desired direction.

It is widely believed that biological tissues evolved to lower the risks of cancer development. One of the specific ways to minimize the chances of tumor formation comes from proper spatial organization of tissues. However, the microscopic mechanisms of underlying processes remain not fully understood. We present a theoretical investigation on the role of spatial structures in cancer initiation dynamics. In our approach, the dynamics of single mutation fixations are analyzed using analytical calculations and computer simulations by mapping them to Moran processes on graphs with different connectivity that mimic various spatial structures. It is found that while the fixation probability is not affected by modifying the spatial structures of the tissues, the fixation times can change dramatically. The slowest dynamics is observed in 'quasi-one-dimensional' structures, while the fastest dynamics is observed in 'quasi-three-dimensional' structures. Theoretical calculations also suggest that there is a critical value of the degree of graph connectivity, which mimics the spatial dimension of the tissue structure, above which the spatial structure of the tissue has no effect on the mutation fixation dynamics. An effective discrete-state stochastic model of cancer initiation is utilized to explain our theoretical results and predictions. Our theoretical analysis clarifies some important aspects on the role of the tissue spatial structures in the cancer initiation processes.

This paper aims to mathematically model the dynamics of Parkinson's disease with therapeutic strategies. The constructed model consists of five state variables: healthy neurons, infected neurons, extracellularα-syn, active microglia, and resting microglia. The qualitative analysis of the model produced an unstable free equilibrium point and a stable endemic equilibrium point. Moreover, these results are validated by numerical experiments with different initial values. Two therapeutic interventions, reduction of extracellularα-syn and reduction of inflammation induced by activated microglia in the central nervous system, are investigated. It is observed that the latter has no apparent effect in delaying the deterioration of neurons. However, treatment to reduce extracellularα-syn preserves neurons and delays the onset of Parkinson's disease, whether alone or in combination with another treatment.

There have been a number of pharmaceutical and non-pharmaceutical interventions associated with COVID-19 over the past two years. Various non-pharmaceutical interventions were proposed and implemented to control the spread of the COVID-19 pandemic. Most common of these were partial and complete lockdowns that were used in an attempt to minimize the costs associated with mortality, economic losses and social factors, while being subject to constraints such as finite hospital capacity. Here, we use a minimal model posed in terms of optimal control theory to understand the costs and benefits of such strategies. This allows us to determine top-down policies for how to restrict social contact rates given an age-structured model for the dynamics of the disease. Depending on the relative weights allocated to mortality and socioeconomic losses, we see that the optimal strategies range from long-term social-distancing only for the most vulnerable, partial lockdown to ensure not over-running hospitals, and alternating-shifts, all of which lead to significant reduction in mortality and/or socioeconomic losses. Crucially, commonly used strategies that involve long periods of broad lockdown are almost never optimal, as they are highly unstable to reopening and entail high socioeconomic costs. Using parameter estimates from data available for Germany and the USA early in the pandemic, we quantify these policies and use sensitivity analysis in the relevant model parameters and initial conditions to determine the range of robustness of our policies. Finally we also discuss how bottom-up behavioral changes affect the dynamics of the pandemic and show how they can work in tandem with top-down control policies to mitigate pandemic costs even more effectively.

Diatoms are unicellular microalgae with a rigid cell wall, able to glide on surfaces by releasing nanopolymeric fibers through central slits known as raphes. Here we consider the modelNitszchia communisto perform quantitative studies on two complementary aspects involved in diatom gliding. Using video microscopy and automated image analysis, we measure the motion of test beads as they are pulled by extracellular polymeric substances (EPS) fibers at the diatom raphe (particle streaming). A multimodal distribution of particle speed is found, evidencing the appearance of short-time events of high speed and acceleration (known as jerky motion) and suggesting that different mechanisms contribute to set diatom velocity during gliding. Furthermore, we use optical tweezers to obtain force-extension records for extracellular diatom nanofibers; records are well described by the worm-like chain model of polymer elasticity. In contrast to previous studies based on application of denaturing force (in the nN regime), application of low force (up to 6 pN) and using enable us to obtain the persistence length of intact fibers. From these measurements, mechanical parameters of EPS fibers such as radius and elastic constant are estimated. Furthermore, by modeling particle streaming as a spring in parallel with a dashpot, we show that the time involved in the release of mechanical energy after fiber detachment from beads (elastic snapping) agrees with our observations of jerky motion. We conclude that the smooth and jerky motions displayed by gliding diatoms correspond to molecular motors and elastic snapping, respectively, thus providing quantitative elements that incorporate to current models of the mechanics behind diatom locomotion.

We develop a theory for thermodynamic instabilities of complex fluids composed of many interacting chemical species organised in families. This model includes partially structured and partially random interactions and can be solved exactly using tools from random matrix theory. The model exhibits three kinds of fluid instabilities: one in which the species form a condensate with a local density that depends on their family (family condensation); one in which species demix in two phases depending on their family (family demixing); and one in which species demix in a random manner irrespective of their family (random demixing). We determine the critical spinodal density of the three types of instabilities and find that the critical spinodal density is finite for both family condensation and family demixing, while for random demixing the critical spinodal density grows as the square root of the number of species. We use the developed framework to describe phase-separation instability of the cytoplasm induced by a change in pH.