Boolean modeling of mechanosensitive epithelial to mesenchymal transition and its reversal.

The significance of biophysical modulators of the epithelial to mesenchymal transition (EMT) is demonstrated by experiments that document full EMT on stiff, nano-patterned substrates in the absence of biochemical induction. Yet, current models focus on biochemical triggers of EMT without addressing its mechanosensitive nature. Here, we built a Boolean model of EMT triggered by mechanosensing - mitogen crosstalk. Our model reproduces epithelial, hybrid E/M and mesenchymal phenotypes, the role of autocrine TGFß signaling in maintaining mesenchymal cells in the absence of external drivers, inhibition of proliferation by TGFß, and its apoptotic effects on soft ECM. We offer testable predictions on the density-dependence of partial EMT, its molecular drivers, and the conflict between mitosis and hybrid E/M stability. Our model opens the door to modeling the effects of the biomechanical environment on cancer cell stemness linked to the hybrid E/M state, as well as the mutually inhibitory crosstalk between EMT and senescence.

Boolean model of anchorage dependence and contact inhibition points to coordinated inhibition but semi-independent induction of proliferation and migration.

Epithelial cells respond to their physical neighborhood with mechano-sensitive behaviors required for development and tissue maintenance. These include anchorage dependence, matrix stiffness-dependent proliferation, contact inhibition of proliferation and migration, and collective migration that balances cell crawling with the maintenance of cell junctions. While required for development and tissue repair, these coordinated responses to the microenvironment also contribute to cancer metastasis. Predictive models of the signaling networks that coordinate these behaviors are critical in controlling cell behavior to halt disease. Here we propose a Boolean regulatory network model that synthesizes mechanosensitive signaling that links anchorage to a matrix of varying stiffness and cell density sensing to contact inhibition, proliferation, migration, and apoptosis. Our model can reproduce anchorage dependence and anoikis, detachment-induced cytokinesis errors, the effect of matrix stiffness on proliferation, and contact inhibition of proliferation and migration by two mechanisms that converge on the YAP transcription factor. In addition, we offer testable predictions related to cell cycle-dependent anoikis sensitivity, the molecular requirements for abolishing contact inhibition, and substrate stiffness dependent expression of the catalytic subunit of PI3K. Moreover, our model predicts heterogeneity in migratory vs. non-migratory phenotypes in sub-confluent monolayers, and co-inhibition but semi-independent induction of proliferation vs. migration as a function of cell density and mitogenic stimulation. Our model serves as a stepping-stone towards modeling mechanosensitive routes to the epithelial to mesenchymal transition, capturing the effects of the mesenchymal state on anoikis resistance, and understanding the balance between migration versus proliferation at each stage of the epithelial to mesenchymal transition.