  We generated a computational approach to analyze the biomechanics of
epithelial cell aggregates, either island or stripes or entire monolayers, that
combines both vertex and contact-inhibition-of-locomotion models to include
both cell-cell and cell-substrate adhesion. Examination of the distribution of
cell protrusions (adhesion to the substrate) in the model predicted high order
profiles of cell organization that agree with those previously seen
experimentally. Cells acquired an asymmetric distribution of basal protrusions,
traction forces and apical aspect ratios that decreased when moving from the
edge to the island center. Our in silico analysis also showed that tension on
cell-cell junctions and apical stress is not homogeneous across the island.
Instead, these parameters are higher at the island center and scales up with
island size, which we confirmed experimentally using laser ablation assays and
immunofluorescence. Without formally being a 3-dimensional model, our approach
has the minimal elements necessary to reproduce the distribution of cellular
forces and mechanical crosstalk as well as distribution of principal stress in
cells within epithelial cell aggregates. By making experimental testable
predictions, our approach would benefit the mechanical analysis of epithelial
tissues, especially when local changes in cell-cell and/or cell-substrate
adhesion drive collective cell behavior.
