Simon van Mourik

Simulation of Organ Patterning on the Floral Meristem Using a Polar Auxin Transport Model

An intriguing phenomenon in plant development is the timing and positioning of lateral organ initiation, which is a fundamental aspect of plant architecture. Although important progress has been made in elucidating the role of auxin transport in the vegetative shoot to explain the phyllotaxis of leaf formation in a spiral fashion, a model study of the role of auxin transport in whorled organ patterning in the expanding floral meristem is not available yet. We present an initial simulation approach to study the mechanisms that are expected to play an important role. Starting point is a confocal imaging study of Arabidopsis floral meristems at consecutive time points during flower development. These images reveal auxin accumulation patterns at the positions of the organs, which strongly suggests that the role of auxin in the floral meristem is similar to the role it plays in the shoot apical meristem. This is the basis for a simulation study of auxin transport through a growing floral meristem, which may answer the question whether auxin transport can in itself be responsible for the typical whorled floral pattern. We combined a cellular growth model for the meristem with a polar auxin transport model. The model predicts that sepals are initiated by auxin maxima arising early during meristem outgrowth. These form a pre-pattern relative to which a series of smaller auxin maxima are positioned, which partially overlap with the anlagen of petals, stamens, and carpels. We adjusted the model parameters corresponding to properties of floral mutants and found that the model predictions agree with the observed mutant patterns. The predicted timing of the primordia outgrowth and the timing and positioning of the sepal primordia show remarkable similarities with a developing flower in nature.

A mathematical model for BRI1 mediated signaling in root growth and hypocotyl elongation

Brassinosteroid (BR) signaling is essential for plant growth and development. In Arabidopsis (Arabidopsis thaliana), BRs are perceived by the BRASSINOSTEROID INSENSITIVE1 (BRI1) receptor. Root growth and hypocotyl elongation are convenient downstream physiological outputs of BR signaling. A computational approach was employed to predict root growth solely on the basis of BRI1 receptor activity. The developed mathematical model predicts that during normal root growth, few receptors are occupied with ligand. The model faithfully predicts root growth, as observed in bri1 loss-of-function mutants. For roots, it incorporates one stimulatory and two inhibitory modules, while for hypocotyls, a single inhibitory module is sufficient. Root growth as observed when BRI1 is overexpressed can only be predicted assuming that a decrease occurred in the BRI1 half-maximum response values. Root growth appears highly sensitive to variation in BR concentration and much less to reduction in BRI1 receptor level, suggesting that regulation occurs primarily by ligand availability and biochemical activity.

A quantitative and dynamic model of the Arabidopsis flowering time gene regulatory network.

Various environmental signals integrate into a network of floral regulatory genes leading to the final decision on when to flower. Although a wealth of qualitative knowledge is available on how flowering time genes regulate each other, only a few studies incorporated this knowledge into predictive models. Such models are invaluable as they enable to investigate how various types of inputs are combined to give a quantitative readout. To investigate the effect of gene expression disturbances on flowering time, we developed a dynamic model for the regulation of flowering time in Arabidopsis thaliana. Model parameters were estimated based on expression time-courses for relevant genes, and a consistent set of flowering times for plants of various genetic backgrounds. Validation was performed by predicting changes in expression level in mutant backgrounds and comparing these predictions with independent expression data, and by comparison of predicted and experimental flowering times for several double mutants. Remarkably, the model predicts that a disturbance in a particular gene has not necessarily the largest impact on directly connected genes. For example, the model predicts that SUPPRESSOR OF OVEREXPRESSION OF CONSTANS (SOC1) mutation has a larger impact on APETALA1 (AP1), which is not directly regulated by SOC1, compared to its effect on LEAFY (LFY) which is under direct control of SOC1. This was confirmed by expression data. Another model prediction involves the importance of cooperativity in the regulation of APETALA1 (AP1) by LFY, a prediction supported by experimental evidence. Concluding, our model for flowering time gene regulation enables to address how different quantitative inputs are combined into one quantitative output, flowering time.

Continuous-time modeling of cell fate determination in Arabidopsis flowers

BackgroundThe 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.ResultsWe 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 Arabidopsis thaliana. 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.ConclusionsThe 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 Arabidopsis, 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.

Mutational Robustness of Gene Regulatory Networks

Mutational robustness of gene regulatory networks refers to their ability to generate constant biological output upon mutations that change network structure. Such networks contain regulatory interactions (transcription factor – target gene interactions) but often also protein-protein interactions between transcription factors. Using computational modeling, we study factors that influence robustness and we infer several network properties governing it. These include the type of mutation, i.e. whether a regulatory interaction or a protein-protein interaction is mutated, and in the case of mutation of a regulatory interaction, the sign of the interaction (activating vs. repressive). In addition, we analyze the effect of combinations of mutations and we compare networks containing monomeric with those containing dimeric transcription factors. Our results are consistent with available data on biological networks, for example based on evolutionary conservation of network features. As a novel and remarkable property, we predict that networks are more robust against mutations in monomer than in dimer transcription factors, a prediction for which analysis of conservation of DNA binding residues in monomeric vs. dimeric transcription factors provides indirect evidence.

Prediction uncertainty assessment of a systems biology model requires a sample of the full probability distribution of its parameters

Multi-parameter models in systems biology are typically ‘sloppy’: some parameters or combinations of parameters may be hard to estimate from data, whereas others are not. One might expect that parameter uncertainty automatically leads to uncertain predictions, but this is not the case. We illustrate this by showing that the prediction uncertainty of each of six sloppy models varies enormously among different predictions. Statistical approximations of parameter uncertainty may lead to dramatic errors in prediction uncertainty estimation. We argue that prediction uncertainty assessment must therefore be performed on a per-prediction basis using a full computational uncertainty analysis. In practice this is feasible by providing a model with a sample or ensemble representing the distribution of its parameters. Within a Bayesian framework, such a sample may be generated by a Markov Chain Monte Carlo (MCMC) algorithm that infers the parameter distribution based on experimental data. Matlab code for generating the sample (with the Differential Evolution Markov Chain sampler) and the subsequent uncertainty analysis using such a sample, is supplied as Supplemental Information.

Switching Control for a Class of Non-Linear Systems with an Application to Post-Harvest Food Storage

For a class of scalar nonlinear systems with switching input a controller is designed using design theory for linear systems. A stability criterion is derived that contains all the physical system parameters, allowing a stability analysis without the need for numerical simulation. The results are motivated by and applied to a model of a bulk storage room for food products. It is shown that for this model a controller with excellent robustness and performance properties can be designed.

A Mathematical Model for the Coreceptors SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE1 and SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE3 in BRASSINOSTEROID INSENSITIVE1-Mediated Signaling

Complex physiological responses of the brassinosteroid pathway can be explained by a biochemical model based on kinase loss-of-function alleles with brassinosteroid receptor complex occupancy as the central estimated parameter. Brassinosteroids (BRs) are key regulators in plant growth and development. The main BR-perceiving receptor in Arabidopsis (Arabidopsis thaliana) is BRASSINOSTEROID INSENSITIVE1 (BRI1). Seedling root growth and hypocotyl elongation can be accurately predicted using a model for BRI1 receptor activity. Genetic evidence shows that non-ligand-binding coreceptors of the SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE (SERK) family are essential for BRI1 signal transduction. A relatively simple biochemical model based on the properties of SERK loss-of-function alleles explains complex physiological responses of the BRI1-mediated BR pathway. The model uses BRI1-BR occupancy as the central estimated parameter and includes BRI1-SERK interaction based on mass action kinetics and accurately describes wild-type root growth and hypocotyl elongation. Simulation studies suggest that the SERK coreceptors primarily act to increase the magnitude of the BRI1 signal. The model predicts that only a small number of active BRI1-SERK complexes are required to carry out BR signaling at physiological ligand concentration. Finally, when calibrated with single mutants, the model predicts that roots of the serk1serk3 double mutant are almost completely brassinolide (BL) insensitive, while the double mutant hypocotyls remain sensitive. This points to residual BRI1 signaling or to a different coreceptor requirement in shoots.

Tracing the molecular basis of transcriptional dynamics in noisy data by using an experiment-based mathematical model

Changes in transcription factor levels, epigenetic status, splicing kinetics and mRNA degradation can each contribute to changes in the mRNA dynamics of a gene. We present a novel method to identify which of these processes is changed in cells in response to external signals or as a result of a diseased state. The method employs a mathematical model, for which the kinetics of gene regulation, splicing, elongation and mRNA degradation were estimated from experimental data of transcriptional dynamics. The time-dependent dynamics of several species of adipose differentiation-related protein (ADRP) mRNA were measured in response to ligand activation of the transcription factor peroxisome proliferator-activated receptor δ (PPARδ). We validated the method by monitoring the mRNA dynamics upon gene activation in the presence of a splicing inhibitor. Our mathematical model correctly identifies splicing as the inhibitor target, despite the noise in the data.

Identifying Optimal Models to Represent Biochemical Systems

Biochemical systems involving a high number of components with intricate interactions often lead to complex models containing a large number of parameters. Although a large model could describe in detail the mechanisms that underlie the system, its very large size may hinder us in understanding the key elements of the system. Also in terms of parameter identification, large models are often problematic. Therefore, a reduced model may be preferred to represent the system. Yet, in order to efficaciously replace the large model, the reduced model should have the same ability as the large model to produce reliable predictions for a broad set of testable experimental conditions. We present a novel method to extract an “optimal” reduced model from a large model to represent biochemical systems by combining a reduction method and a model discrimination method. The former assures that the reduced model contains only those components that are important to produce the dynamics observed in given experiments, whereas the latter ensures that the reduced model gives a good prediction for any feasible experimental conditions that are relevant to answer questions at hand. These two techniques are applied iteratively. The method reveals the biological core of a model mathematically, indicating the processes that are likely to be responsible for certain behavior. We demonstrate the algorithm on two realistic model examples. We show that in both cases the core is substantially smaller than the full model.