Data_Sheet_1_Modeling the Dynamics of Let-7-Coupled Gene Regulatory Networks Linking Cell Proliferation to Malignant Transformation.pdf

Let-7 microRNA controls the expression of proteins that belong to two distinct gene regulatory networks, namely, a cyclin-dependent kinase (Cdk) network driving the cell cycle and a cell transformation network that can undergo an epigenetic switch between a non-transformed and a malignant transformed cell state. Using mathematical modeling and transcriptomic data analysis, we here investigate how Let-7 controls the Cdk-dependent cell cycle network, and how it couples the latter with the transformation network. We also assess the consequence of this coupling on cancer progression. Our analysis shows that the switch from a quiescent to a proliferative state depends on the relative levels of Let-7 and several cell cycle activators. Numerical simulations further indicate that the Let-7-coupled cell cycle and transformation networks mutually control each other, and our model identifies key players for this mutual control. Transcriptomic data analysis from The Cancer Genome Atlas (TCGA) suggests that the two networks are activated in cancer, in particular in gastrointestinal cancers, and that the activation levels vary significantly among patients affected by a same cancer type. Our mathematical model, when applied to a heterogeneous cell population, suggests that heterogeneity among tumors may in part result from stochastic switches between a non-transformed cell state with low proliferative capability and a transformed cell state with high proliferative property. The model further predicts that Let-7 may reduce tumor heterogeneity by decreasing the occurrence of stochastic switches toward a transformed, proliferative cell state. In conclusion, we identified the key components responsible for the qualitative dynamics of two networks interconnected by Let-7. The two networks are heterogeneously activated in several cancers, thereby stressing the need to consider patient’s specific characteristics to optimize therapeutic strategies.