Elementary effects for models with dimensional inputs of arbitrary type and range: Scaling and trajectory generation.

The Elementary Effects method is a global sensitivity analysis approach for identifying (un)important parameters in a model. However, it has almost exclusively been used where inputs are dimensionless and take values on [0, 1]. Here, we consider models with dimensional inputs, inputs taking values on arbitrary intervals or discrete inputs.

Long-distance hormone transport via the phloem.

Several key plant hormones are synthesised in the shoot and are advected within the phloem to the root tip. In the root tip, these hormones regulate growth and developmental processes, and responses to environmental cues. However, we lack understanding of how environmental factors and biological parameters affect the delivery of hormones to the root tip.

evidence for the utility of parsimonious root phenotypes for improved vegetative growth and carbon sequestration under drought.

Drought is a primary constraint to crop yields and climate change is expected to increase the frequency and severity of drought stress in the future. It has been hypothesized that crops can be made more resistant to drought and better able to sequester atmospheric carbon in the soil by selecting appropriate root phenotypes. We introduce , an upgraded version of the functional-structural plant/soil model , and use it to test the utility of a maize root phenotype with fewer and steeper axial roots, reduced lateral root branching density, and more aerenchyma formation (i.

Uncertainty and error in SARS-CoV-2 epidemiological parameters inferred from population-level epidemic models.

During the SARS-CoV-2 pandemic, epidemic models have been central to policy-making. Public health responses have been shaped by model-based projections and inferences, especially related to the impact of various non-pharmaceutical interventions. Accompanying this has been increased scrutiny over model performance, model assumptions, and the way that uncertainty is incorporated and presented.

Modeling root loss reveals impacts on nutrient uptake and crop development.

Despite the widespread prevalence of root loss in plants, its effects on crop productivity are not fully understood. While root loss reduces the capacity of plants to take up water and nutrients from the soil, it may provide benefits by decreasing the resources required to maintain the root system. Here, we simulated a range of root phenotypes in different soils and root loss scenarios for barley (Hordeum vulgare), common bean (Phaseolus vulgaris), and maize (Zea mays) using and extending the open-source, functional-structural root/soil simulation model OpenSimRoot.

Modelling the emergence of cities and urban patterning using coupled integro-differential equations.

Human residential population distributions show patterns of higher density clustering around local services such as shops and places of employment, displaying characteristic length scales; Fourier transforms and spatial autocorrelation show the length scale between UK cities is around 45 km. We use integro-differential equations to model the spatio-temporal dynamics of population and service density under the assumption that they benefit from spatial proximity, captured via spatial weight kernels. The system tends towards a well-mixed homogeneous state or a spatial pattern.

Simulating Crop Root Systems Using OpenSimRoot.

Functional-structural plant models are valuable modeling tools in analyzing plant development. A functional-structural plant model combines a three-dimensional representation of plant structure with models for physiological functions in order to better understand plant development. We present a guide to simulating crop root systems with OpenSimRoot, a feature-rich, highly cited, and open-source functional-structural root architecture model.

The Virtual Root : Mathematical Modeling of Auxin Transport in the Arabidopsis Root Tip Using the Open-Source Software SimuPlant.

Hormone signals like auxin play a critical role controlling plant growth and development. Determining the mechanisms that regulate auxin distribution in cells and tissues is a vital step in understanding this hormone's role during plant development. Recent mathematical models have enabled us to understand the essential role that auxin influx and efflux carriers play in auxin transport in the Arabidopsis root tip (Band et al.

Differential biosynthesis and cellular permeability explain longitudinal gibberellin gradients in growing roots.

Control over cell growth by mobile regulators underlies much of eukaryotic morphogenesis. In plant roots, cell division and elongation are separated into distinct longitudinal zones and both division and elongation are influenced by the growth regulatory hormone gibberellin (GA). Previously, a multicellular mathematical model predicted a GA maximum at the border of the meristematic and elongation zones.

A Mechanistic Investigation into Ischemia-Driven Distal Recurrence of Glioblastoma.

Glioblastoma (GBM) is the most aggressive primary brain tumor with a short median survival. Tumor recurrence is a clinical expectation of this disease and usually occurs along the resection cavity wall. However, previous clinical observations have suggested that in cases of ischemia following surgery, tumors are more likely to recur distally.

Identifying the spatial and temporal dynamics of molecularly-distinct glioblastoma sub-populations.

Glioblastomas (GBMs) are the most aggressive primary brain tumours and have no known cure. Each individual tumour comprises multiple sub-populations of genetically-distinct cells that may respond differently to targeted therapies and may contribute to disappointing clinical trial results. Image-localized biopsy techniques allow multiple biopsies to be taken during surgery and provide information that identifies regions where particular sub-populations occur within an individual GBM, thus providing insight into their regional genetic variability.

Switching behaviour in vascular smooth muscle cell-matrix adhesion during oscillatory loading.

Integrins regulate mechanotransduction between smooth muscle cells (SMCs) and the extracellular matrix (ECM). SMCs resident in the walls of airways or blood vessels are continuously exposed to dynamic mechanical forces due to breathing or pulsatile blood flow. However, the resulting effects of these forces on integrin dynamics and associated cell-matrix adhesion are not well understood.

Speed Switch in Glioblastoma Growth Rate due to Enhanced Hypoxia-Induced Migration.

We analyze the wave speed of the Proliferation Invasion Hypoxia Necrosis Angiogenesis (PIHNA) model that was previously created and applied to simulate the growth and spread of glioblastoma (GBM), a particularly aggressive primary brain tumor. We extend the PIHNA model by allowing for different hypoxic and normoxic cell migration rates and study the impact of these differences on the wave-speed dynamics. Through this analysis, we find key variables that drive the outward growth of the simulated GBM.

Effect of Loading History on Airway Smooth Muscle Cell-Matrix Adhesions.

Integrin-mediated adhesions between airway smooth muscle (ASM) cells and the extracellular matrix (ECM) regulate how contractile forces generated within the cell are transmitted to its external environment. Environmental cues are known to influence the formation, size, and survival of cell-matrix adhesions, but it is not yet known how they are affected by dynamic fluctuations associated with tidal breathing in the intact airway. Here, we develop two closely related theoretical models to study adhesion dynamics in response to oscillatory loading of the ECM, representing the dynamic environment of ASM cells in vivo.

Capturing the Dynamics of a Hybrid Multiscale Cancer Model with a Continuum Model.

Cancer is a complex disease involving processes at spatial scales from subcellular, like cell signalling, to tissue scale, such as vascular network formation. A number of multiscale models have been developed to study the dynamics that emerge from the coupling between the intracellular, cellular and tissue scales. Here, we develop a continuum partial differential equation model to capture the dynamics of a particular multiscale model (a hybrid cellular automaton with discrete cells, diffusible factors and an explicit vascular network).

OpenSimRoot: widening the scope and application of root architectural models.

OpenSimRoot is an open-source, functional-structural plant model and mathematical description of root growth and function. We describe OpenSimRoot and its functionality to broaden the benefits of root modeling to the plant science community. OpenSimRoot is an extended version of SimRoot, established to simulate root system architecture, nutrient acquisition and plant growth.

Root hydrotropism is controlled via a cortex-specific growth mechanism.

Plants can acclimate by using tropisms to link the direction of growth to environmental conditions. Hydrotropism allows roots to forage for water, a process known to depend on abscisic acid (ABA) but whose molecular and cellular basis remains unclear. Here we show that hydrotropism still occurs in roots after laser ablation removed the meristem and root cap.

Mathematical modelling of cytokines, MMPs and fibronectin fragments in osteoarthritic cartilage.

Osteoarthritis (OA) is a degenerative disease which causes pain and stiffness in joints. OA progresses through excessive degradation of joint cartilage, eventually leading to significant joint degeneration and loss of function. Cytokines, a group of cell signalling proteins, present in raised concentrations in OA joints, can be classified into pro-inflammatory and anti-inflammatory groups.

Dynamic regulation of auxin oxidase and conjugating enzymes AtDAO1 and GH3 modulates auxin homeostasis.

The hormone auxin is a key regulator of plant growth and development, and great progress has been made understanding auxin transport and signaling. Here, we show that auxin metabolism and homeostasis are also regulated in a complex manner. The principal auxin degradation pathways in Arabidopsis include oxidation by Arabidopsis thaliana gene DIOXYGENASE FOR AUXIN OXIDATION 1/2 (AtDAO1/2) and conjugation by Gretchen Hagen3s (GH3s).

The spatiotemporal order of signaling events unveils the logic of development signaling.

Motivation: Animals from worms and insects to birds and mammals show distinct body plans; however, the embryonic development of diverse body plans with tissues and organs within is controlled by a surprisingly few signaling pathways. It is well recognized that combinatorial use of and dynamic interactions among signaling pathways follow specific logic to control complex and accurate developmental signaling and patterning, but it remains elusive what such logic is, or even, what it looks like.

Results: We have developed a computational model for Drosophila eye development with innovated methods to reveal how interactions among multiple pathways control the dynamically generated hexagonal array of R8 cells.

Modelling cell cycle synchronisation in networks of coupled radial glial cells.

Radial glial cells play a crucial role in the embryonic mammalian brain. Their proliferation is thought to be controlled, in part, by ATP mediated calcium signals. It has been hypothesised that these signals act to locally synchronise cell cycles, so that clusters of cells proliferate together, shedding daughter cells in uniform sheets.

Comparing stochastic differential equations and agent-based modelling and simulation for early-stage cancer.

There is great potential to be explored regarding the use of agent-based modelling and simulation as an alternative paradigm to investigate early-stage cancer interactions with the immune system. It does not suffer from some limitations of ordinary differential equation models, such as the lack of stochasticity, representation of individual behaviours rather than aggregates and individual memory. In this paper we investigate the potential contribution of agent-based modelling and simulation when contrasted with stochastic versions of ODE models using early-stage cancer examples.

Modelling the coupling between intracellular calcium release and the cell cycle during cortical brain development.

Most neocortical neurons formed during embryonic brain development arise from radial glial cells which communicate, in part, via ATP mediated calcium signals. Although the intercellular signalling mechanisms that regulate radial glia proliferation are not well understood, it has recently been demonstrated that ATP dependent intracellular calcium release leads to an increase of nearly 100% in overall cellular proliferation. It has been hypothesised that cytoplasmic calcium accelerates entry into S phase of the cell cycle and/or acts to recruit otherwise quiescent cells onto the cell cycle.

On-lattice agent-based simulation of populations of cells within the open-source Chaste framework.

Over the years, agent-based models have been developed that combine cell division and reinforced random walks of cells on a regular lattice, reaction-diffusion equations for nutrients and growth factors; and ordinary differential equations for the subcellular networks regulating the cell cycle. When linked to a vascular layer, this multiple scale model framework has been applied to tumour growth and therapy. Here, we report on the creation of an agent-based multi-scale environment amalgamating the characteristics of these models within a Virtual Physiological Human (VPH) Exemplar Project.

STAT1-cooperative DNA binding distinguishes type 1 from type 2 interferon signaling.

STAT1 is an indispensable component of a heterotrimer (ISGF3) and a STAT1 homodimer (GAF) that function as transcription regulators in type 1 and type 2 interferon signaling, respectively. To investigate the importance of STAT1-cooperative DNA binding, we generated gene-targeted mice expressing cooperativity-deficient STAT1 with alanine substituted for Phe77. Neither ISGF3 nor GAF bound DNA cooperatively in the STAT1F77A mouse strain, but type 1 and type 2 interferon responses were affected differently.