ECB-ART-41248

BMC Syst Biol
2009 Aug 23;3:83. doi: 10.1186/1752-0509-3-83.
Show Gene links
Show Anatomy links

## Monte Carlo analysis of an ODE Model of the Sea Urchin Endomesoderm Network.

Abstract

BACKGROUND: Gene Regulatory Networks (GRNs) control the differentiation, specification and function of cells at the genomic level. The levels of interactions within large GRNs are of enormous depth and complexity. Details about many GRNs are emerging, but in most cases it is unknown to what extent they control a given process, i.e. the grade of completeness is uncertain. This uncertainty stems from limited experimental data, which is the main bottleneck for creating detailed dynamical models of cellular processes. Parameter estimation for each node is often infeasible for very large GRNs. We propose a method, based on random parameter estimations through Monte-Carlo simulations to measure completeness grades of GRNs. RESULTS: We developed a heuristic to assess the completeness of large GRNs, using ODE simulations under different conditions and randomly sampled parameter sets to detect parameter-invariant effects of perturbations. To test this heuristic, we constructed the first ODE model of the whole sea urchin endomesoderm GRN, one of the best studied large GRNs. We find that nearly 48% of the parameter-invariant effects correspond with experimental data, which is 65% of the expected optimal agreement obtained from a submodel for which kinetic parameters were estimated and used for simulations. Randomized versions of the model reproduce only 23.5% of the experimental data. CONCLUSION: The method described in this paper enables an evaluation of network topologies of GRNs without requiring any parameter values. The benefit of this method is exemplified in the first mathematical analysis of the complete Endomesoderm Network Model. The predictions we provide deliver candidate nodes in the network that are likely to be erroneous or miss unknown connections, which may need additional experiments to improve the network topology. This mathematical model can serve as a scaffold for detailed and more realistic models. We propose that our method can be used to assess a completeness grade of any GRN. This could be especially useful for GRNs involved in human diseases, where often the amount of connectivity is unknown and/or many genes/interactions are missing.

PubMed ID: 19698179

PMC ID: PMC2739852

Article link: BMC Syst Biol

Species referenced: Echinodermata

Genes referenced: alx1 arid3a ets1 gata6 hnf6 LOC100887844 LOC100893907 LOC115925415 LOC575170 LOC583082 LOC592057 LOC594353 onecut2 otx2 pmar1 tbr1

Article Images: [+] show captions

References [+] :

Ben-Tabou de-Leon,

*Modeling the dynamics of transcriptional gene regulatory networks for animal development.*2009, Pubmed