Approximate Bayesian Computation (ABC) is a widely applicable and popular approach to estimating unknown parameters of mechanistic models. As ABC analyses are computationally expensive, ...parallelization on high-performance infrastructure is often necessary. However, the existing parallelization strategies leave computing resources unused at times and thus do not optimally leverage them yet. We present look-ahead scheduling, a wall-time minimizing parallelization strategy for ABC Sequential Monte Carlo algorithms, which avoids idle times of computing units by preemptive sampling of subsequent generations. This allows to utilize all available resources. The strategy can be integrated with e.g. adaptive distance function and summary statistic selection schemes, which is essential in practice. Our key contribution is the theoretical assessment of the strategy of preemptive sampling and the proof of unbiasedness. Complementary, we provide an implementation and evaluate the strategy on different problems and numbers of parallel cores, showing speed-ups of typically 10-20% and up to 50% compared to the best established approach, with some variability. Thus, the proposed strategy allows to improve the cost and run-time efficiency of ABC methods on high-performance infrastructure.
Abstract We introduce a novel lattice-gas cellular automaton (LGCA) for compressible vectorial active matter with polar and nematic velocity alignment. Interactions are, by construction, zero-range. ...For polar alignment, we show the system undergoes a phase transition that promotes aggregation with strong resemblance to the classic zero-range process. We find that above a critical point, the states of a macroscopic fraction of the particles in the system coalesce into the same state, sharing the same position and momentum (polar condensate). For nematic alignment, the system also exhibits condensation, but there exist fundamental differences: a macroscopic fraction of the particles in the system collapses into a filament, where particles possess only two possible momenta. Furthermore, we derive hydrodynamic equations for the active LGCA model to understand the phase transitions and condensation that undergoes the system. We also show that generically the discrete lattice symmetries—e.g. of a square or hexagonal lattice—affect drastically the emergent large-scale properties of on-lattice active systems. The study puts in evidence that aligning active matter on the lattice displays new behavior, including phase transitions to states that share similarities to condensation models.
The secretion of osmolytes into a lumen and thereby caused osmotic water inflow can drive fluid flows in organs without a mechanical pump. Such fluids include saliva, sweat, pancreatic juice and ...bile. The effects of elevated fluid pressure and the associated mechanical limitations of organ function remain largely unknown since fluid pressure is difficult to measure inside tiny secretory channels in vivo. We consider the pressure profile of the coupled osmolyte-flow problem in a secretory channel with a closed tip and an open outlet. Importantly, the entire lateral boundary acts as a dynamic fluid source, the strength of which self-organizes through feedback from the emergent pressure solution itself. We derive analytical solutions and compare them to numerical simulations of the problem in three-dimensional space. The theoretical results reveal a phase boundary in a four-dimensional parameter space separating the commonly considered regime with steady flow all along the channel, here termed "wet-tip" regime, from a "dry-tip" regime suffering ceased flow downstream from the closed tip. We propose a relation between the predicted phase boundary and the onset of cholestasis, a pathological liver condition with reduced bile outflow. The phase boundary also sets an intrinsic length scale for the channel which could act as a length sensor during organ growth.
Early disease diagnosis is key to the effective treatment of diseases. Histopathological analysis of human biopsies is the gold standard to diagnose tissue alterations. However, this approach has low ...resolution and overlooks 3D (three-dimensional) structural changes resulting from functional alterations. Here, we applied multiphoton imaging, 3D digital reconstructions and computational simulations to generate spatially resolved geometrical and functional models of human liver tissue at different stages of non-alcoholic fatty liver disease (NAFLD). We identified a set of morphometric cellular and tissue parameters correlated with disease progression, and discover profound topological defects in the 3D bile canalicular (BC) network. Personalized biliary fluid dynamic simulations predicted an increased pericentral biliary pressure and micro-cholestasis, consistent with elevated cholestatic biomarkers in patients' sera. Our spatially resolved models of human liver tissue can contribute to high-definition medicine by identifying quantitative multiparametric cellular and tissue signatures to define disease progression and provide new insights into NAFLD pathophysiology.
Cell fate reprogramming, such as the generation of insulin-producing β cells from other pancreas cells, can be achieved by external modulation of key transcription factors. However, the known gene ...regulatory interactions that form a complex network with multiple feedback loops make it increasingly difficult to design the cell reprogramming scheme because the linear regulatory pathways as schemes of causal influences upon cell lineages are inadequate for predicting the effect of transcriptional perturbation. However, sufficient information on regulatory networks is usually not available for detailed formal models. Here we demonstrate that by using the qualitatively described regulatory interactions as the basis for a coarse-grained dynamical ODE (ordinary differential equation) based model, it is possible to recapitulate the observed attractors of the exocrine and β, δ, α endocrine cells and to predict which gene perturbation can result in desired lineage reprogramming. Our model indicates that the constraints imposed by the incompletely elucidated regulatory network architecture suffice to build a predictive model for making informed decisions in choosing the set of transcription factors that need to be modulated for fate reprogramming.
Simulations of tissue-specific effects of primary acute viral infections like COVID-19 are essential for understanding disease outcomes and optimizing therapies. Such simulations need to support ...continuous updating in response to rapid advances in understanding of infection mechanisms, and parallel development of components by multiple groups. We present an open-source platform for multiscale spatiotemporal simulation of an epithelial tissue, viral infection, cellular immune response and tissue damage, specifically designed to be modular and extensible to support continuous updating and parallel development. The base simulation of a simplified patch of epithelial tissue and immune response exhibits distinct patterns of infection dynamics from widespread infection, to recurrence, to clearance. Slower viral internalization and faster immune-cell recruitment slow infection and promote containment. Because antiviral drugs can have side effects and show reduced clinical effectiveness when given later during infection, we studied the effects on progression of treatment potency and time-of-first treatment after infection. In simulations, even a low potency therapy with a drug which reduces the replication rate of viral RNA greatly decreases the total tissue damage and virus burden when given near the beginning of infection. Many combinations of dosage and treatment time lead to stochastic outcomes, with some simulation replicas showing clearance or control (treatment success), while others show rapid infection of all epithelial cells (treatment failure). Thus, while a high potency therapy usually is less effective when given later, treatments at late times are occasionally effective. We illustrate how to extend the platform to model specific virus types (e.g., hepatitis C) and add additional cellular mechanisms (tissue recovery and variable cell susceptibility to infection), using our software modules and publicly-available software repository.
Functional tissue architecture originates by self-assembly of distinct cell types, following tissue-specific rules of cell-cell interactions. In the liver, a structural model of the lobule was ...pioneered by Elias in 1949. This model, however, is in contrast with the apparent random 3D arrangement of hepatocytes. Since then, no significant progress has been made to derive the organizing principles of liver tissue. To solve this outstanding problem, we computationally reconstructed 3D tissue geometry from microscopy images of mouse liver tissue and analyzed it applying soft-condensed-matter-physics concepts. Surprisingly, analysis of the spatial organization of cell polarity revealed that hepatocytes are not randomly oriented but follow a long-range liquid-crystal order. This does not depend exclusively on hepatocytes receiving instructive signals by endothelial cells, since silencing Integrin-β1 disrupted both liquid-crystal order and organization of the sinusoidal network. Our results suggest that bi-directional communication between hepatocytes and sinusoids underlies the self-organization of liver tissue.
Activation of ErbB receptors by epidermal growth factor (EGF) or heregulin (HRG) determines distinct cell-fate decisions, although signals propagate through shared pathways. Using mathematical ...modeling and experimental approaches, we unravel how HRG and EGF generate distinct, all-or-none responses of the phosphorylated transcription factor c-Fos. In the cytosol, EGF induces transient and HRG induces sustained ERK activation. In the nucleus, however, ERK activity and
c-fos mRNA expression are transient for both ligands. Knockdown of dual-specificity phosphatases extends HRG-stimulated nuclear ERK activation, but not
c-fos mRNA expression, implying the existence of a HRG-induced repressor of
c-fos transcription. Further experiments confirmed that this repressor is mainly induced by HRG, but not EGF, and requires new protein synthesis. We show how a spatially distributed, signaling-transcription cascade robustly discriminates between transient and sustained ERK activities at the c-Fos system level. The proposed control mechanisms are general and operate in different cell types, stimulated by various ligands.
Display omitted
► ErbB ligands EGF and HRG induce transient versus sustained cytoplasmic ERK activity ► Nuclear ERK activity and
c-fos mRNA expression are transient for both ligands ► HRG induces increased expression of
c-fos, ERK phosphatases, and a
c-fos repressor ► Signaling-transcription cascades robustly discriminate different ERK dynamics
The mechanisms of organ size control remain poorly understood. A key question is how cells collectively sense the overall status of a tissue. We addressed this problem focusing on mouse liver ...regeneration. Using digital tissue reconstruction and quantitative image analysis, we found that the apical surface of hepatocytes forming the bile canalicular network expands concomitant with an increase in F‐actin and phospho‐myosin, to compensate an overload of bile acids. These changes are sensed by the Hippo transcriptional co‐activator YAP, which localizes to apical F‐actin‐rich regions and translocates to the nucleus in dependence of the integrity of the actin cytoskeleton. This mechanism tolerates moderate bile acid fluctuations under tissue homeostasis, but activates YAP in response to sustained bile acid overload. Using an integrated biophysical–biochemical model of bile pressure and Hippo signaling, we explained this behavior by the existence of a mechano‐sensory mechanism that activates YAP in a switch‐like manner. We propose that the apical surface of hepatocytes acts as a self‐regulatory mechano‐sensory system that responds to critical levels of bile acids as readout of tissue status.
Synopsis
This study shows that the apical surface of hepatocytes acts as a self‐regulatory mechano‐sensory system that responds to critical levels of bile acids as readout of tissue status and activates YAP, by a mechanism dependent on the actin cytoskeleton.
During regeneration, the bile canalicular network expands concomitant with an increase of F‐actin and phospho‐Myosin, to compensate an overload of bile acids.
The bile canaliculi expansion is sensed by the transcriptional co‐activator YAP, which localizes to apical F‐actin‐rich regions of hepatocytes and translocates to the nucleus depending on the actin cytoskeleton.
An integrated biophysical‐biochemical model of bile pressure and Hippo signalling suggests a mechano‐sensory mechanism that activates YAP in a switch‐like manner upon bile acid overload.
This study shows that the apical surface of hepatocytes acts as a self‐regulatory mechano‐sensory system that responds to critical levels of bile acids as readout of tissue status and activates YAP, by a mechanism dependent on the actin cytoskeleton.
Chemical reaction networks are ubiquitous in biology, and their dynamics is fundamentally stochastic. Here, we present the software library pSSAlib, which provides a complete and concise ...implementation of the most efficient partial-propensity methods for simulating exact stochastic chemical kinetics. pSSAlib can import models encoded in Systems Biology Markup Language, supports time delays in chemical reactions, and stochastic spatiotemporal reaction-diffusion systems. It also provides tools for statistical analysis of simulation results and supports multiple output formats. It has previously been used for studies of biochemical reaction pathways and to benchmark other stochastic simulation methods. Here, we describe pSSAlib in detail and apply it to a new model of the endocytic pathway in eukaryotic cells, leading to the discovery of a stochastic counterpart of the cut-out switch motif underlying early-to-late endosome conversion. pSSAlib is provided as a stand-alone command-line tool and as a developer API. We also provide a plug-in for the SBMLToolbox. The open-source code and pre-packaged installers are freely available from http://mosaic.mpi-cbg.de.