Efficient Estimation of the Robustness Region of Biological Models with Oscillatory Behavior

Robustness is an essential feature of biological systems, and any mathematical model that describes such a system should reflect this feature. Especially, persistence of oscillatory behavior is an important issue. A benchmark model for this phenomenon is the Laub-Loomis model, a nonlinear model for cAMP oscillations in Dictyostelium discoideum. This model captures the most important features of biomolecular networks oscillating at constant frequencies. Nevertheless, the robustness of its oscillatory behavior is not yet fully understood. Given a system that exhibits oscillating behavior for some set of parameters, the central question of robustness is how far the parameters may be changed, such that the qualitative behavior does not change. The determination of such a “robustness region” in parameter space is an intricate task. If the number of parameters is high, it may be also time consuming. In the literature, several methods are proposed that partially tackle this problem. For example, some methods only detect particular bifurcations, or only find a relatively small box-shaped estimate for an irregularly shaped robustness region. Here, we present an approach that is much more general, and is especially designed to be efficient for systems with a large number of parameters. As an illustration, we apply the method first to a well understood low-dimensional system, the Rosenzweig-MacArthur model. This is a predator-prey model featuring satiation of the predator. It has only two parameters and its bifurcation diagram is available in the literature. We find a good agreement with the existing knowledge about this model. When we apply the new method to the high dimensional Laub-Loomis model, we obtain a much larger robustness region than reported earlier in the literature. This clearly demonstrates the power of our method. From the results, we conclude that the biological system underlying is much more robust than was realized until now

The Effect of Chemical Information on the Spatial Distribution of Fruit Flies: II Parameterization, Calibration, and Sensitivity

In a companion paper (Lof et al., in Bull. Math. Biol., 2008), we describe a spatio-temporal model for insect behavior. This model includes chemical information for finding resources and conspecifics. As a model species, we used Drosophila melanogaster, because its behavior is documented comparatively well

The Effect of Chemical Information on the Spatial Distribution of Fruit Flies: I Model Results

Animal aggregation is a general phenomenon in ecological systems. Aggregations are generally considered as an evolutionary advantageous state in which members derive the benefits of protection and mate choice, balanced by the costs of limiting resources and competition. In insects, chemical information conveyance plays an important role in finding conspecifics and forming aggregations. In this study, we describe a spatio-temporal simulation model designed to explore and quantify the effects of these infochemicals, i.e., food odors and an aggregation pheromone, on the spatial distribution of a fruit fly (Drosophila melanogaster) population, where the lower and upper limit of local population size are controlled by an Allee effect and competition. We found that during the spatial expansion and strong growth of the population, the use of infochemicals had a positive effect on population size. The positive effects of reduced mortality at low population numbers outweighed the negative effects of increased mortality due to competition. At low resource densities, attraction toward infochemicals also had a positive effect on population size during recolonization of an area after a local population crash, by decreasing the mortality due to the Allee effect. However, when the whole area was colonized and the population was large, the negative effects of competition on population size were larger than the positive effects of the reduction in mortality due to the Allee effect. The use of infochemicals thus has mainly positive effects on population size and population persistence when the population is small and during the colonization of an area

Continuous-time modeling of cell fate determination in Arabidopsis flowers

<p>Abstract</p> <p>Background</p> <p>The genetic control of floral organ specification is currently being investigated by various approaches, both experimentally and through modeling. Models and simulations have mostly involved boolean or related methods, and so far a quantitative, continuous-time approach has not been explored.</p> <p>Results</p> <p>We propose an ordinary differential equation (ODE) model that describes the gene expression dynamics of a gene regulatory network that controls floral organ formation in the model plant <it>Arabidopsis thaliana</it>. In this model, the dimerization of MADS-box transcription factors is incorporated explicitly. The unknown parameters are estimated from (known) experimental expression data. The model is validated by simulation studies of known mutant plants.</p> <p>Conclusions</p> <p>The proposed model gives realistic predictions with respect to independent mutation data. A simulation study is carried out to predict the effects of a new type of mutation that has so far not been made in <it>Arabidopsis</it>, but that could be used as a severe test of the validity of the model. According to our predictions, the role of dimers is surprisingly important. Moreover, the functional loss of any dimer leads to one or more phenotypic alterations.</p

Search for dark matter produced in association with bottom or top quarks in √s = 13 TeV pp collisions with the ATLAS detector

A search for weakly interacting massive particle dark matter produced in association with bottom or top quarks is presented. Final states containing third-generation quarks and miss- ing transverse momentum are considered. The analysis uses 36.1 fb−1 of proton–proton collision data recorded by the ATLAS experiment at √s = 13 TeV in 2015 and 2016. No significant excess of events above the estimated backgrounds is observed. The results are in- terpreted in the framework of simplified models of spin-0 dark-matter mediators. For colour- neutral spin-0 mediators produced in association with top quarks and decaying into a pair of dark-matter particles, mediator masses below 50 GeV are excluded assuming a dark-matter candidate mass of 1 GeV and unitary couplings. For scalar and pseudoscalar mediators produced in association with bottom quarks, the search sets limits on the production cross- section of 300 times the predicted rate for mediators with masses between 10 and 50 GeV and assuming a dark-matter mass of 1 GeV and unitary coupling. Constraints on colour- charged scalar simplified models are also presented. Assuming a dark-matter particle mass of 35 GeV, mediator particles with mass below 1.1 TeV are excluded for couplings yielding a dark-matter relic density consistent with measurements