Monte Carlo analysis of an ODE Model of the Sea Urchin Endomesoderm Network.
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.


CONCLUSIONS
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.
Kühn
Clemens
C
Max-Planck-Institute for Molecular Genetics, Ihnestr 63-73, 14195 Berlin, Germany. clemens.kuehn@biologie.hu-berlin.de
Wierling
Christoph
C
Kühn
Alexander
A
Klipp
Edda
E
Panopoulou
Georgia
G
Lehrach
Hans
H
Poustka
Albert J
AJ
D000818:Animals; D003198:Computer Simulation; D053263:Gene Regulatory Networks; D007700:Kinetics; D008648:Mesoderm / Q000378:metabolism*; D009010:Monte Carlo Method; D015203:Reproducibility of Results; D012617:Sea Urchins / Q000196:embryology / Q000235:genetics*
Development. 2001 Feb;128(3):365-75
11152635
Nucleic Acids Res. 2006 Jan 1;34(Database issue):D689-91
16381960
Dev Biol. 1997 Jul 15;187(2):253-66
9242422
Nat Struct Mol Biol. 2008 Nov;15(11):1192-8
18849996
Dev Biol. 2006 Dec 1;300(1):74-89
17055477
Nat Rev Mol Cell Biol. 2008 Oct;9(10):770-80
18797474
Brief Funct Genomic Proteomic. 2007 Sep;6(3):240-51
17942476
Bioessays. 2002 Jun;24(6):553-63
12111738
Genome Inform. 2007;18:75-84
18546476
Dev Biol. 2005 May 15;281(2):286-98
15893979
Dev Biol. 2010 Apr 15;340(2):170-8
19523466
Proc Natl Acad Sci U S A. 2007 Jul 24;104(30):12383-8
17636127
Dev Biol. 2006 May 15;293(2):513-25
16581059
Biochim Biophys Acta. 2009 Apr;1789(4):363-74
18757046
Dev Biol. 2002 Jun 1;246(1):162-90
12027441
Science. 2004 Mar 5;303(5663):1538-42
15001784
Cell. 2004 Sep 17;118(6):675-85
15369668
Development. 2001 Mar;128(5):617-29
11171388
Dev Biol. 2009 Jan 15;325(2):317-28
19028486
Gene Expr Patterns. 2004 Dec;5(2):161-5
15567710
J Theor Biol. 2006 Oct 7;242(3):774-89
16765384
Genesis. 2004 Jul;39(3):194-205
15282746
Dev Biol. 2002 Jun 1;246(1):209-28
12027443
Development. 2006 Nov;133(21):4173-81
17038513
Proc Natl Acad Sci U S A. 2008 Apr 22;105(16):5955-62
18413610
Nat Protoc. 2008;3(4):710-8
18388954
Trends Biotechnol. 2003 Jun;21(6):255-62
12788545