Analysis of a reduced model of epithelial-mesenchymal fate determination in cancer metastasis as a singularly-perturbed monotone system

Tumor metastasis is one of the main factors responsible for the high fatality rate of cancer. Metastasis can occur after malignant cells transition from the epithelial phenotype to the mesenchymal phenotype. This transformation allows cells to migrate via the circulatory system and subsequently settle in distant organs after undergoing the reverse transition from the mesenchymal to the epithelial phenotypes. The core gene regulatory network controlling these transitions consists of a system made up of coupled SNAIL/miRNA-34 and ZEB1/miRNA-200 subsystems. In this work, we formulate a mathematical model of the core regulatory motif and analyze its long-term behavior. We start by developing a detailed reaction network with 24 state variables. Assuming fast promoter and mRNA kinetics, we then show how to reduce our model to a monotone four-dimensional system. For the reduced system, monotone dynamical systems theory can be used to prove generic convergence to the set of equilibria for all bounded trajectories. The theory does not apply to the full model, which is not monotone, but we briefly discuss results for singularly-perturbed monotone systems that provide a tool to extend convergence results from reduced to full systems, under appropriate time separation assumptions.

[1]  Robert J. Plemmons,et al.  Nonnegative Matrices in the Mathematical Sciences , 1979, Classics in Applied Mathematics.

[2]  David Angeli,et al.  New Approach to the Stability of Chemical Reaction Networks: Piecewise Linear in Rates Lyapunov Functions , 2014, IEEE Transactions on Automatic Control.

[3]  Eshel Ben-Jacob,et al.  MicroRNA-based regulation of epithelial–hybrid–mesenchymal fate determination , 2013, Proceedings of the National Academy of Sciences.

[4]  Jeremy Gunawardena,et al.  Time‐scale separation – Michaelis and Menten's old idea, still bearing fruit , 2014, The FEBS journal.

[5]  Eduardo Sontag,et al.  Untangling the wires: A strategy to trace functional interactions in signaling and gene networks , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[6]  Christopher Jones,et al.  Geometric singular perturbation theory , 1995 .

[7]  R. Huang,et al.  Epithelial-Mesenchymal Transitions in Development and Disease , 2009, Cell.

[8]  E. Ben-Jacob,et al.  Tristability in cancer-associated microRNA-TF chimera toggle switch. , 2013, The journal of physical chemistry. B.

[9]  Péter Érdi,et al.  Mathematical Models of Chemical Reactions: Theory and Applications of Deterministic and Stochastic Models , 1989 .

[10]  Liming Wang,et al.  Singularly Perturbed Monotone Systems and an Application to Double Phosphorylation Cycles , 2008, J. Nonlinear Sci..

[11]  Hassan K. Khalil,et al.  Singular perturbation methods in control : analysis and design , 1986 .

[12]  M. Hirsch,et al.  4. Monotone Dynamical Systems , 2005 .

[13]  Gunter Meister,et al.  Argonaute proteins at a glance , 2010, Journal of Cell Science.

[14]  Eduardo D. Sontag,et al.  Monotone and near-monotone biochemical networks , 2007, Systems and Synthetic Biology.

[15]  David Angeli,et al.  A computational framework for a Lyapunov-enabled analysis of biochemical reaction networks , 2020, PLoS computational biology.

[16]  Samy Lamouille,et al.  Molecular mechanisms of epithelial–mesenchymal transition , 2014, Nature Reviews Molecular Cell Biology.

[17]  Neil Fenichel.  Geometric singular perturbation theory for ordinary differential equations , 1979 .

[18]  Alexei A. Sharov,et al.  Database for mRNA Half-Life of 19 977 Genes Obtained by DNA Microarray Analysis of Pluripotent and Differentiating Mouse Embryonic Stem Cells , 2008, DNA research : an international journal for rapid publication of reports on genes and genomes.

[19]  Jun Liu,et al.  Characterization of the mammalian miRNA turnover landscape , 2015, Nucleic acids research.

[20]  Kaspar Nipp,et al.  Smooth attractive invariant manifolds of singularly perturbed ODE's , 1992 .

[21]  Domitilla Del Vecchio,et al.  Multi-modality in gene regulatory networks with slow promoter kinetics , 2017, PLoS Comput. Biol..

[22]  David Angeli,et al.  Graph-theoretic characterizations of monotonicity of chemical networks in reaction coordinates , 2010, Journal of mathematical biology.

[23]  Hal L. Smith,et al.  Monotone Dynamical Systems: An Introduction To The Theory Of Competitive And Cooperative Systems (Mathematical Surveys And Monographs) By Hal L. Smith , 1995 .

[24]  J. Lei.  Monotone Dynamical Systems , 2013 .

[25]  G. Berx,et al.  Regulatory networks defining EMT during cancer initiation and progression , 2013, Nature Reviews Cancer.

[26]  M. Hirsch.  Stability and convergence in strongly monotone dynamical systems. , 1988 .

[27]  B. Kholodenko,et al.  The dynamic control of signal transduction networks in cancer cells , 2015, Nature Reviews Cancer.

[28]  R. Weinberg,et al.  Emerging Biological Principles of Metastasis , 2017, Cell.

[29]  U. Alon.  An introduction to systems biology : design principles of biological circuits , 2019 .

[30]  M. Hirsch.  Systems of Differential Equations that are Competitive or Cooperative II: Convergence Almost Everywhere , 1985 .