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