The Hedgehog (Hh) signalling cascade is highly conserved and involved in development and disease throughout evolution. Nevertheless, in comparison with other pathways, our mechanistic understanding ...of Hh signal transduction is remarkably incomplete. In the absence of ligand, the Hh receptor Patched (Ptc) represses the key signal transducer Smoothened (Smo) through an unknown mechanism. Hh binding to Ptc alleviates this repression, causing Smo redistribution to the plasma membrane, phosphorylation and opening of the Smo cytoplasmic tail, and Smo oligomerisation. However, the order and interdependence of these events is as yet poorly understood. We have mathematically modelled and simulated Smo activation for two alternative modes of pathway activation, with Ptc primarily affecting either Smo localisation or phosphorylation. Visualising Smo activation through a novel, fluorescence-based reporter allowed us to test these competing models. Here, we show that Smo localisation to the plasma membrane is sufficient for phosphorylation of the cytoplasmic tail in the presence of Ptc. Using fluorescence cross-correlation spectroscopy (FCCS), we also demonstrate that inactivation of Ptc by Hh induces Smo clustering irrespective of Smo phosphorylation. Our observations therefore support a model of Hh signal transduction whereby Smo subcellular localisation and not phosphorylation is the primary target of Ptc function.
Lattice-gas cellular automata (LGCAs) can serve as stochastic mathematical models for collective behavior (e.g. pattern formation) emerging in populations of interacting cells. In this paper, a ...two-phase optimization algorithm for global parameter estimation in LGCA models is presented. In the first phase, local minima are identified through gradient-based optimization. Algorithmic differentiation is adopted to calculate the necessary gradient information. In the second phase, for global optimization of the parameter set, a multi-level single-linkage method is used. As an example, the parameter estimation algorithm is applied to a LGCA model for early in vitro angiogenic pattern formation.
In many animals, regenerative processes can replace lost body parts. Organ and tissue regeneration consequently also hold great medical promise. The regulation of regenerative processes is achieved ...through concerted actions of multiple organizational levels of the organism, from diffusing molecules and cellular gene expression patterns up to tissue mechanics. Our intuition is usually not adapted well to this degree of complexity and the quantitative aspects of the regulation of regenerative processes remain poorly understood. One way out of this dilemma lies in the combination of experimentation and mathematical modeling within an iterative process of model development/refinement, model predictions for novel experimental conditions, quantitative experiments testing these predictions, and subsequent model refinement. This interdisciplinary approach has already provided key insights into smaller scale processes during embryonic development and a so-far limited number of more complex regeneration processes. This review discusses selected theoretical and interdisciplinary studies and is structured along the three phases of regeneration: (1) initiation of a regeneration response, (2) tissue patterning during regenerate growth, (3) arresting the regeneration response. Moreover, we highlight the opportunities provided by extensions of mathematical models from developmental processes toward the study of related regenerative processes.
Genes and human languages are discrete combinatorial systems (DCSs), in which the basic building blocks are finite sets of elementary units: nucleotides or codons in a DNA sequence, and letters or ...words in a language. Different combinations of these finite units give rise to potentially infinite numbers of genes or sentences. This type of DCSs can be represented as an alphabetic bipartite network (ABN) where there are two kinds of nodes, one type represents the elementary units while the other type represents their combinations. Here, we extend and generalize recent analytical findings for ABNs derived in Peruani, Europhys. Lett. 79, 28001 (2007) and empirically investigate two real world systems in terms of ABNs, the codon gene and the phoneme-language network. The one-mode projections onto the elementary basic units are also studied theoretically as well as in real world ABNs. We propose the use of ABNs as a means for inferring the mechanisms underlying the growth of real world DCSs.
Research software has become a central asset in academic research. It optimizes existing and enables new research methods, implements and embeds research knowledge, and constitutes an essential ...research product in itself. Research software must be sustainable in order to understand, replicate, reproduce, and build upon existing research or conduct new research effectively. In other words, software must be available, discoverable, usable, and adaptable to new needs, both now and in the future. Research software therefore requires an environment that supports sustainability. Hence, a change is needed in the way research software development and maintenance are currently motivated, incentivized, funded, structurally and infrastructurally supported, and legally treated. Failing to do so will threaten the quality and validity of research. In this paper, we identify challenges for research software sustainability in Germany and beyond, in terms of motivation, selection, research software engineering personnel, funding, infrastructure, and legal aspects. Besides researchers, we specifically address political and academic decision-makers to increase awareness of the importance and needs of sustainable research software practices. In particular, we recommend strategies and measures to create an environment for sustainable research software, with the ultimate goal to ensure that software-driven research is valid, reproducible and sustainable, and that software is recognized as a first class citizen in research. This paper is the outcome of two workshops run in Germany in 2019, at deRSE19 - the first International Conference of Research Software Engineers in Germany - and a dedicated DFG-supported follow-up workshop in Berlin.
Spiral and antispiral waves are studied numerically in two examples of oscillatory reaction−diffusion media and analytically in the corresponding complex Ginzburg−Landau equation (CGLE). We argue ...that both these structures are sources of waves in oscillatory media, which are distinguished only by the sign of the phase velocity of the emitted waves. Using known analytical results in the CGLE, we obtain a criterion for the CGLE coefficients that predicts whether antispirals or spirals will occur in the corresponding reaction−diffusion systems. We apply this criterion to the FitzHugh−Nagumo and Brusselator models by deriving the CGLE near the Hopf bifurcations of the respective equations. Numerical simulations of the full reaction−diffusion equations confirm the validity of our simple criterion near the onset of oscillations. They also reveal that antispirals often occur near the onset and turn into spirals further away from it. The transition from antispirals to spirals is characterized by a divergence in the wavelength. A tentative interpretation of recent experimental observations of antispiral waves in the Belousov−Zhabotinsky reaction in a microemulsion is given.
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.
Nonlinear waves emitted from a moving source are studied. A meandering spiral in a reaction-diffusion medium provides an example in which waves originate from a source exhibiting a back-and-forth ...movement in a radial direction. The periodic motion of the source induces a Doppler effect that causes a modulation in wavelength and amplitude of the waves ("superspiral"). Using direct simulations as well as numerical nonlinear analysis within the complex Ginzburg-Landau equation, we show that waves subject to a convective Eckhaus instability can exhibit monotonic growth or decay as well as saturation of these modulations depending on the perturbation frequency. Our findings elucidate recent experimental observations concerning superspirals and their decay to spatiotemporal chaos.