Stochasticity and robustness analysis of microRNA-mediated ERK signaling network

MicroRNAs (miRNAs) play a critical role in regulating the signaling network such as the extracellular signal regulated kinase (ERK) pathway. However, the mechanisms of miRNA-mediated posttranscriptional regulations with regard to their impacts on signaling require further investigation from a systematic level. Therefore, we develop a mathematical model by analyzing the precise dynamic mechanisms with a bevy of miRNAs involved in the Ras/Raf/MEK/ERK pathway. A systems-based analysis approach is introduced into this model and the dynamics have been implemented deterministically and stochastically. Our analysis reveals that miRNAs are key participants regulating the gene expression of ERK network, and the cooperative actions of miRNAs are important in keeping the normal biological characteristics and amplification effects of the system. Meanwhile, the appearance of system disorder in the absence of miRNAs suggests that miRNAs may play a role in the pathological processes, such as tumor and inflammation. The sensitivity analysis and the kinetic parameter perturbation show that the binding of receptor (EGFR) and adaptor protein (Shc, Grb2 and Sos) under normal physiological conditions is crucial for the robustness of the whole pathway. In addition, the stochastic dynamic patterns are in a good agreement with the deterministic results, further demonstrating that the variability of the system due to the presence of some stochastic noise is low. All these will be helpful for a deeper understanding of the dynamic mechanism of miRNA-mediated ERK signal network, which might present a rich area of future research with the relevant regulatory roles of miRNAs in cell signaling pathway.

[1]  C. Widmann,et al.  Mitogen-activated protein kinase: conservation of a three-kinase module from yeast to human. , 1999, Physiological reviews.

[2]  T K Christopoulos,et al.  Signal amplification system for DNA hybridization assays based on in vitro expression of a DNA label encoding apoaequorin. , 1999, Nucleic acids research.

[3]  E. O’Shea,et al.  Global analysis of protein expression in yeast , 2003, Nature.

[4]  Yunqing Li,et al.  microRNA-7 inhibits the epidermal growth factor receptor and the Akt pathway and is down-regulated in glioblastoma. , 2008, Cancer research.

[5]  Xiaodong Cai,et al.  Exact stochastic simulation of coupled chemical reactions with delays. , 2007, The Journal of chemical physics.

[6]  Oliver E. Sturm,et al.  Computational modelling of the receptor-tyrosine-kinase-activated MAPK pathway. , 2005, The Biochemical journal.

[7]  B. Kholodenko,et al.  Quantification of Short Term Signaling by the Epidermal Growth Factor Receptor* , 1999, The Journal of Biological Chemistry.

[8]  F. Slack,et al.  RAS Is Regulated by the let-7 MicroRNA Family , 2005, Cell.

[9]  D. Fell,et al.  Differential feedback regulation of the MAPK cascade underlies the quantitative differences in EGF and NGF signalling in PC12 cells , 2000, FEBS letters.

[10]  B. Kholodenko,et al.  Negative feedback and ultrasensitivity can bring about oscillations in the mitogen-activated protein kinase cascades. , 2000, European journal of biochemistry.

[11]  Danish Sayed,et al.  MicroRNA-21 targets Sprouty2 and promotes cellular outgrowths. , 2008, Molecular biology of the cell.

[12]  Tony Pawson,et al.  The Grb2-mSos1 Complex Binds Phosphopeptides with Higher Affinity than Grb2* , 1996, The Journal of Biological Chemistry.

[13]  Stefano Piccolo,et al.  MicroRNA control of signal transduction , 2010, Nature Reviews Molecular Cell Biology.

[14]  A. Arkin,et al.  Stochastic mechanisms in gene expression. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[15]  Isil Aksan,et al.  A Computer‐Based Model for the Regulation of Mitogen Activated Protein Kinase (MAPK) Activation , 2003, Journal of receptor and signal transduction research.

[16]  M. Cobb,et al.  Mitogen-activated protein (MAP) kinase pathways: regulation and physiological functions. , 2001, Endocrine reviews.

[17]  Olli Yli-Harja,et al.  Simulation tools for biochemical networks: evaluation of performance and usability , 2005, Bioinform..

[18]  J Godovac-Zimmermann,et al.  Phosphotyrosine 1173 Mediates Binding of the Protein-tyrosine Phosphatase SHP-1 to the Epidermal Growth Factor Receptor and Attenuation of Receptor Signaling* , 1998, The Journal of Biological Chemistry.

[19]  Xia Wang,et al.  MiR-206-mediated dynamic mechanism of the mammalian circadian clock , 2011, BMC Systems Biology.

[20]  Michael A. Gibson,et al.  Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels , 2000 .

[21]  Hiroaki Kitano,et al.  Biological robustness , 2008, Nature Reviews Genetics.

[22]  J B Hoek,et al.  Differential inhibition of epidermal growth factor signaling pathways in rat hepatocytes by long-term ethanol treatment. , 1997, Gastroenterology.

[23]  Yan Li,et al.  Dynamic mechanisms for pre-miRNA binding and export by Exportin-5. , 2011, RNA.

[24]  E. Gilles,et al.  Computational modeling of the dynamics of the MAP kinase cascade activated by surface and internalized EGF receptors , 2002, Nature Biotechnology.

[25]  Yan Li,et al.  Stochastic simulations of the cytochrome P450 catalytic cycle. , 2007, The journal of physical chemistry. B.

[26]  S. Kimura,et al.  A computational model on the modulation of mitogen-activated protein kinase (MAPK) and Akt pathways in heregulin-induced ErbB signalling. , 2003, The Biochemical journal.

[27]  Steffen Hauptmann,et al.  Activation of mitogen-activated protein kinase is required for migration and invasion of placental site trophoblastic tumor. , 2005, The American journal of pathology.

[28]  W. Rottbauer,et al.  MicroRNA-21 contributes to myocardial disease by stimulating MAP kinase signalling in fibroblasts , 2008, Nature.

[29]  D. Engelman,et al.  Two EGF molecules contribute additively to stabilization of the EGFR dimer , 1997, The EMBO journal.

[30]  G. Hannon,et al.  Control of translation and mRNA degradation by miRNAs and siRNAs. , 2006, Genes & development.

[31]  Chi-Ying F. Huang,et al.  Ultrasensitivity in the mitogen-activated protein kinase cascade. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[32]  D. Gillespie A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[33]  Gozoh Tsujimoto,et al.  MicroRNA-34a Inhibits Cell Proliferation by Repressing Mitogen-Activated Protein Kinase Kinase 1 during Megakaryocytic Differentiation of K562 Cells , 2010, Molecular Pharmacology.

[34]  D. Mcnamara,et al.  Determinants of Substrate Recognition in the Protein-tyrosine Phosphatase, PTP1 (*) , 1996, The Journal of Biological Chemistry.

[35]  Q. Cui,et al.  Principles of microRNA regulation of a human cellular signaling network , 2006, Molecular systems biology.

[36]  Sonia Duprey,et al.  Global sensitivity analysis of the joint kinematics during gait to the parameters of a lower limb multi-body model , 2015, Medical & Biological Engineering & Computing.

[37]  H. Schaeffer,et al.  Mitogen-Activated Protein Kinases: Specific Messages from Ubiquitous Messengers , 1999, Molecular and Cellular Biology.