Next Article in Journal
Clinical Characteristics and Postoperative Complications in Patients Undergoing Colorectal Cancer Surgery with Perioperative COVID-19 Infection
Next Article in Special Issue
Regulation of Tumor Microenvironment through YAP/TAZ under Tumor Hypoxia
Previous Article in Journal
Recent Pre-Clinical Advancements in Nuclear Medicine: Pioneering the Path to a Limitless Future
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Interaction of Mechanics and the Hippo Pathway in Drosophila melanogaster

1
Department of Mathematics, University of California, Riverside, CA 92507, USA
2
School of Mathematics, University of Minnesota, Minneapolis, MN 55455, USA
*
Author to whom correspondence should be addressed.
Cancers 2023, 15(19), 4840; https://doi.org/10.3390/cancers15194840
Submission received: 17 August 2023 / Revised: 10 September 2023 / Accepted: 15 September 2023 / Published: 3 October 2023
(This article belongs to the Special Issue Hippo Signaling in Cancer)

Abstract

:

Simple Summary

Control of tissue growth is an important component of cancer treatment, but understanding how to do this effectively is still in the early stages. Tissue growth typically involves both complex interacting signal transduction networks as well as mechanical interactions within and between cells, and how they interact to control growth is an open problem. Herein we begin the process of how to understand the interactions by focusing on tissue growth in the fruit fly Drosophila melanogaster. A model of the Hippo pathway, which is one of the primary growth-control pathways, is integrated with a description of the mechanical behavior of cells in wing-disc tissue to predict how mechanics and signaling interact, and a number of significant insights and predictions have emerged from the analysis of the model.

Abstract

Drosophila melanogaster has emerged as an ideal system for studying the networks that control tissue development and homeostasis and, given the similarity of the pathways involved, controlled and uncontrolled growth in mammalian systems. The signaling pathways used in patterning the Drosophila wing disc are well known and result in the emergence of interaction of these pathways with the Hippo signaling pathway, which plays a central role in controlling cell proliferation and apoptosis. Mechanical effects are another major factor in the control of growth, but far less is known about how they exert their control. Herein, we develop a mathematical model that integrates the mechanical interactions between cells, which occur via adherens and tight junctions, with the intracellular actin network and the Hippo pathway so as to better understand cell-autonomous and non-autonomous control of growth in response to mechanical forces.

1. Introduction

A long-standing problem in biology is to understand how the size and shape of multicellular organisms are controlled by the integration of signal transduction pathways that respond to organism-level growth control signals, local extra- and intracellular biochemical signals, and mechanical forces. In particular, understanding cell growth control in a given context is essential for understanding cancer and how particular types may be treated. For example, in the wing disc, ecdysone functions to regulate the disc size via the mTor pathway [1], and mutations or hyperactivation of signal transduction pathways such as the PI3K/Akt/mTOR pathway, the Hippo pathway, and the Ras/MAPK pathway are frequently involved in cancer. A difficulty in identifying type-specific treatments lies in the fact that signal transduction pathways such as these are often highly integrated with positive or negative feedback steps, which often provides a degree of signal transduction redundancy that makes discovering effective treatments difficult. For example, there is evidence that Akt is negatively regulated by Hippo signaling [2] (see Figure 1) and that inactivating Akt prevents entry of the co-transcriptional factor Yorkie into the nucleus [3]. Furthermore, the PI3K and MAPK pathways are strongly interconnected through a number of positive and negative feedback loops, and as a result, targeted inhibition of mTORC1 by treatment with rapamycin can cause MAPK reactivation and can lead to resistance to single mTORC1 inhibition. Consequently, combinatorial targeting of mTOR plus MAPK may induce a better response to therapies [4,5].
Further complicating the establishment of sufficient understanding to develop treatments is the more recently discovered fact that mechanical effects in the form of cell–cell and cell–ECM forces can also play an important role in controlling cell division and tissue size. Recent studies have shown that the nuclear translocation of yes-associated protein (Yap), a transcriptional co-activator which functions as a downstream effector of the Hippo pathway in mammals, is directly controlled by mechanical forces exerted through focal adhesions to the nucleus [6], which establishes a direct connection between substrate stiffness and cell growth. Furthermore, in epithelial tissues, strong cell–cell connections in the form of adherens junctions govern force transmission between cells and are essential in regulating contact-induced growth arrest and tissue integrity and homeostasis [7]. The primary components of the adherens junctions are the cadherin family proteins, and during tumorigenesis, one of the most important processes is the downregulation of E-cadherin, which leads to tumor cell migration and invasion [8]. Knowing the regulatory roles of the mechanical effects on the signaling pathways and the downstream gene expressions will better assist us both in understanding how cancer develops and in designing appropriate treatment strategies. Among the signaling pathways that are affected by cadherin cell–cell adhesion and contribute to tumor progression, the Hippo pathway is perhaps the best understood, due in large part to the high degree of homology between the Hippo pathway in Drosophila, particularly in the wing disc, and the Hippo pathway in mammalian cells (compared in Figure 2). A mathematical model of the Hippo pathway in the wing disc of Drosophila that can explain many of the observed results in mutant transplant experiments was developed earlier [9], and here, we extend this model and integrate it with a model for cell–cell interactions via E-cadherins and intracellular interactions via the actin network.
Figure 1. Interaction of the mTor and Hippo pathways in a mammalian cell. Modified from [10].
Figure 1. Interaction of the mTor and Hippo pathways in a mammalian cell. Modified from [10].
Cancers 15 04840 g001
Figure 2. Comparison of the Hippo pathway in Drosophila and in mammals. From [11] with permission.
Figure 2. Comparison of the Hippo pathway in Drosophila and in mammals. From [11] with permission.
Cancers 15 04840 g002
Our focus is on the development and analysis of a model of the Drosophila wing disc that integrates components of the known biochemical pathways and mechanical effects at the level of growth control so as to determine how normal control functions under various conditions and predict the consequences of dysfunction in pathways. A first step toward understanding the balances between biochemical pathways controlling Hippo activity was successful [9], and a major objective of this paper is to determine whether the current understanding of the biochemical and mechanical effects can explain effects such as disc over- or undergrowth following various interventions.
Several earlier models on integrating mechanics and signaling have been proposed [12,13,14,15]. In [12], the nonuniform growth in a layer of tissue was modeled, and it was proposed that feedback from mechanical stress provides a mechanism for stabilizing uniform growth. In [13], the authors coupled growth factor-induced growth in the disc and stretching-induced proliferation in the periphery of the disc and provided an explanation of the growth cessation in the Drosophila wing disc. Elsewhere, the authors modeled a regulatory network in the disc and coupled interactions between mechanical forces and certain growth factors [15]. A growth cessation mechanism depending on the absolute compression level in the center of the disc and the compression gradient across the disc was proposed. In the model described below, we will include the detailed cell–cell contact through adherens junctions and critical biochemical components in the Hippo pathway. Since disc growth involves the types of interactions found in many other epithelial tissue systems, we anticipate that detailed modeling of the disc will lead to general insights into development in mammalian systems and, in particular, insights into cancer development in epithelial tissues.

2. The Biochemical Regulation of the Hippo Pathway

The Hippo pathway is a highly conserved core kinase cascade that is regulated by many upstream factors and in turn regulates transcription factors involved in cell growth, proliferation, and apoptosis. The key effector in the Hippo pathway is Yorkie (Yki), a co-transcription factor whose nuclear localization is controlled by the kinase Warts (Wts). When Yki is phosphorylated by Warts, it cannot enter the nucleus. In the nucleus, Yki acts by binding to transcription factors such as Scalloped (Sd) to activate the expression of cyclin E, myc, DIAP1, and bantam, which control cell proliferation [16], and it also controls the expression of genes upstream in the Hippo pathway, such as Warts, Expanded, Fat, Ds, Merlin, Kibra, and four-jointed (Fj) (compared in Figure 3). Since Expanded (Ex), Merlin (Mer), Kibra, Fj, and Wts are involved in the regulation of Yki via the Hippo signaling, this establishes a negative feedback loop to regulate its own activity [17,18,19].
Two atypical cadherins, Fat (Ft) and Dachsous (Ds), are key upstream components in the Hippo pathway. Both Ft and Ds are transmembrane molecules that can be localized to the apical cell membrane. The intracellular domains (ICDs) of each can independently mediate signaling through the Hippo pathway within a cell, while Ft and Ds on adjacent cell membranes can bind together to form heterodimers which mediate cell–cell interaction and convey information between neighboring cells.
It was found that much of the Hippo signaling transmitted by Fat depends on the ability of its ICD to alter the subapical localization and activity of Dachs, an unconventional myosin. Dlish/Vamana, an SH3 domain-containing protein, binds with Fat and the Dachs regulator Approximated (App). When palmitoylated by App, Dlish helps localize Dachs to the subapical membrane, where Dachs can promote Warts degradation. Fat reduces the activity of App through direct binding and thus inhibits the subapical localization of Dachs [20]. Localized Dachs can bind the kinase Warts (Wts) and promote its degradation [21]. While overexpression of Dachs increases wing size, loss of Dachs completely suppresses the ability of a Fat mutant to induce overgrowth or target gene overexpression [22,23]. As Fat-Ds signaling also regulates the planar cell polarity (PCP), experiments suggest that while the polarization of Dachs controlled by Fat and Ds is essential for cell polarization, the amount of Dachs localized on the membrane controls cell growth [20].
The ICDs of Ds regulate growth through physical interaction with Riquiqui (Riq), an amino acid protein that contains five WD40 repeats that mediate protein–protein interaction, and Riq binds with Minibrain (Mnb), a DYRK family kinase. Ds is required for localization of Riq at the apical junction of cells in wing discs, and localized Riq does not affect Wts directly but rather potentiates the ability of Mnb to inhibit Wts through phosphorylation, thereby increasing Yki activity. Depletion of Riq or Mnb reduces Yki activity and results in smaller discs, while overexpression of Riq or Mnb increases the wing size [24]. Recent studies also suggest that Dlish/Vamana associates with the ICDs of Ds and Dachs, which provides another link from Ds to the regulation of Hippo signaling [25].
Interaction between Ds and Fat is a component of cell–cell interaction that is regulated by another protein, Four-jointed (Fj), which itself is one of the target genes of the Hippo pathway. Fj functions as a kinase that phosphorylates potential extracellular cadherin domains of Fat and Ds in the Golgi [26], thereby enhancing binding of Fat to Ds and inhibiting the binding of Ds to Fat [27]. However, the observed weaker phenotype of fj mutants compared with ds mutants, as well the high levels of Fat and Ds aggregation in the absence of exogenous Fj, indicates that Fat and Ds binding could occur even without Fj [27,28]. It was also found that either loss of Fj or overexpression of Fj across the disc causes modest reductions in wing size [29,30].
In addition to Fat-Ds signaling, another group of cell cortex-localized proteins, including Ex, Mer, and Kibra, also plays an important role in the regulation of the Hippo pathway through activation of Wts [31,32,33]. The apical-localized transmembrane protein Crumb, which is crucial in establishing and maintaining epithelial apical-basal polarity, directly binds to Ex through its FBM domain, and this localizes Ex to the apical membrane in epithelial cells, where Ex acts like a Wts activator [34,35]. In vivo, most Wts co-localizes with its inhibitor, the LIM protein Ajuba (Jub), at adherens junctions when Hippo signaling is inactive. When activated, Wts shifts from Jub co-localization to Ex association [36]. Ex can also directly bind with Yki to sequester Yki to the membrane and prevent its nuclear activity [33,37]. This effect was found to be independent of the Wts-mediated phosphorylation of Yki [37].
Merlin, a member of the ERM family proteins, and Kibra, another tumor suppressor, were found to activate Hippo signaling at the medial apical cortex in parallel with Ex. Together, they recruit the adapter protein Salvador (Sav) and Hpo, the Drosophila ortholog of the mammalian MST1 and MST2 kinases, to promote the activation of the Hippo pathway through phosphorylation and activation of Wts [32,38,39]. Another study also suggested that Mer could directly bind Wts at the plasma membrane, which promotes the phosphorylation of Wts by an Hpo-Sav complex. Although the details are still unclear, this Mer-Wts binding was shown to be promoted with the disruption of the actin cytoskeleton, which provides another link between the mechanical cues and growth [40]. As Crumbs forms complexes with Ex, Mer, and Kibra at apical junctions, which sequesters Kibra away from the medial apical cortex, it was suggested that Crumbs has a dual effect on Hippo signaling: it promotes the pathway activity through Ex recruitment and represses it by inhibiting the activity of Kibra [38].
The interactions between Ex and Fat in the Hippo pathway regulation remain poorly defined. A reduced apical Ex level was observed in a Fat mutant in spite of the increased ex transcription caused by Yki activity, and Fat is thought to be required for the normal localization and stabilization of Ex at adherens junctions [18,41,42]. Recent studies showed that Dlish/Vamana binds with Fat and Ex, and in vivo, it reduces the subapical accumulation of Ex. Loss of Dlish inhibits the destabilization of Ex caused by loss of Fat [20,43]. Direct binding between Ex and a highly conserved region of Fat ICDs, in parallel to Fat-Ds binding, was also found to recruit Ex to the apical junction and stabilize it [44]. However, another school of thought is that Fat and Ex function in parallel and have additive effects on imaginal disc growth and development. Although a Fat mutant induces partial loss of Ex protein from the membrane, mutation of Fat still promotes growth with overexpressed and accumulated membrane-localized Ex [45].
There are multiple time scales in the processes described above, ranging from extremely fast binding reactions to slower gene expression, and in the following model construction, we assume that gene expression is much slower compared with direct protein-protein binding or relocalization. This implies that the feedback regulation through gene expression controlled by Yki activity has a minor effect on the timescale under consideration, and thus it is not included in the model described below.

Model Description

The complexity of the interactions described above precludes construction of a model that includes all of them, and to include only the essential steps in the regulation of the Hippo pathway, we developed a model that was modified from an earlier one described in [9]. The previous model focused on the collaboration of upstream regulators Fat and Ds and provided explanations to the seemingly contradictory experimental observations from Fat/Ds mutant experiments. The current model extends the earlier one by incorporating both regulatory mechanisms from Fat/Ds and Ex, and it is built around Fat, Ds, Fj, Dachs, Riq, Ex, Wts, and Yki. As the regulatory mechanism from Mer/Kibra was found to parallel that of Ex, their inputs are assumed to be constant for the problem of interest. Similarly, the binding between Crumb and Ex is not modeled explicitly.
The signaling network we considered is shown in Figure 3, and a detailed description of the model and experimental justification of the assumptions can be found in Appendix A.1. As the known regulation of the Wts inhibitor Jub is mainly from biomechanical signals, we first study the rest of the regulatory network for biochemical regulations and consider the integration in later sections. A brief summary of the major assumptions used in the model is given below:
1.
Fat and Ds bind together on adjacent membranes, where the cytosolic Fj acts on Fat to promote the binding and on Ds to inhibit it.
2.
Most protein–protein interactions are described by reversible binding and dissociation steps, except for those between Fat and Dachs, Ds and Riq, and Fj and Fat/Ds and the inhibition of Ex on the active Yki level. For instance, the inhibition by Fat of membrane localization of Dachs is modeled by a reduced binding rate and described by a decreasing Hill function of total membrane-bound Fat. The effect of Ds on facilitating the junctional localization of Riq is modeled in a similar way but with an increasing Hill function.
3.
The activation of Wts by Ex is modeled through direct, reversible binding of the two molecules. The level of the Ex-Wts complex is assumed to be proportional to the active Wts that is involved in Yki phosphorylation. The inhibitory effect of Ex on nuclear Yki activity is modeled as a separate inhibitory term in the governing equation for Yki.
4.
As the interactions between Fat and Ex are still poorly defined, we test the case where Fat stabilizes the junctional localization of Ex and compare it with the case where Fat and Ex act in parallel.
5.
As stated above, the feedback loops through the genetic regulations of Yki on Ex, Wts, and Fj were not incorporated in the current model.
The governing equations for all species can be found in Appendix A.1, and a representative equation that governs the Wts dynamics is shown below:
d [ W ] c i d t = k W A + [ W ] c i j = 1 N m [ A m ] m , j i + k W A j = 1 N m [ W A ¯ ] m , j i k W R + [ W ] c i j = 1 N m [ R ] m , j i + k W R j = 1 N m [ W R ¯ ] m , j i k E W + [ W ] c i j = 1 N m [ E x ] m , j i + k E W j = 1 N m [ E W ¯ ] m , j i ,
Here k W A ± , k W R ± , and k E W ± are the binding/unbinding rates between Wts and Dachs, Riq, and Ex, respectively, and N m represents the number of separate membrane regions in the model. The mechanical regulation of Wts through Jub is not incorporated in the current equation. In accord with the previous model in [9], the current model predicts the inhibitory effect of Fat on the membrane localization of Dachs and the enhancement of membrane localization of Riq in Figure 4A.
Membrane-localized Dachs binds with Wts and promotes its degradation, while membrane Riq potentiates the Mnb-dependent Wts phosphorylation and reduces its activity. Together, they regulate the downstream level of active Wts, and in turn regulate the Yki activity. In addition, as the current model also includes the dynamics of Ex, we also tested how the Fat and Ex interaction affects the downstream level of Yki. In Figure 4B, the Yki level is shown as a function of the Fat production with the stabilization of Ex by Fat. When Fat and Ex function in parallel, Fat can only affect the Yki level through Dachs and Riq, where one is inhibitory and one has a promotional effect. When Fat also acts on membrane Ex, which promotes the activation of Wts, the inhibitory effect of Fat on Yki is strengthened, as shown in Figure 4B. In this case, further complications are added to the regulatory effect of Fat on the Hippo pathway. Furthermore, in the current model, we did not consider the effect of the Crb complex and the impact of Ex on the Hippo pathway through the interactions with Crb, and therefore Ex plays a solely inhibitory role on Yki.
The biochemical model presented here is a simplified version of the upstream regulatory network of the Hippo pathway. As noted earlier, the interactions between Ex and the Crb complex are neglected in the current model, as are their interactions with Hippo/Salvador, which phosphorylates and regulates the activity of Wts. The current model can be further extended by including the regulatory effect of these components and analyzing their impact on the downstream effector. As the Hippo pathway is highly conserved in Drosophila and mammals, a similar model framework can be used to study the regulatory mechanisms of the Hippo pathway in mammalian cells, and this will further help with the development of therapeutic strategies for cancer treatment.

3. Force Transmission within and between Cells

In an epithelial tissue such as the wing disc, tissue-level mechanical interactions are mediated via adherens junctions (AJs) and tight junctions (TJs), utilizing the former to transmit mechanical stress between cells and the latter to allow intercellular exchange of water, ions, and macromolecules and to block mixing of apical and basolateral lipids [46]. These contacts are critical for the establishment of apical-basal polarity, cell fate, and differentiation. In addition, as mentioned earlier, the atypical cadherins Fat and Dachsous form heterophilic bonds of the AJs apically and are involved in planar cell polarity [47]. Cell-level forces that arise in the cytoskeleton (CSK)—which is comprised of actin filaments, actomyosin, and other filament-based structures—are transmitted to the AJs via the cortical actomyosin network. This circumferential cortical bundle produces forces both normal and tangential to the membrane that are transmitted through the AJs, and these cell-to-cell connections in turn mediate tissue-level forces [48,49]. Much of the actomyosin in the wing disc is found in the cortical bundle—also called an actin belt in disc cells—which is located just below the apical membrane, and the apical actomyosin network [50]. The remainder resides in the CSK and in the actomyosin bound to lower portions of the cell membrane (see Figure 5). During growth, many regulators of actomyosin are upregulated, and the process of cell division involves precise regulation of actomyosin contractility [51]. Actomyosin levels are spatially variable throughout the disc, being lower in the center and elevated in the periphery and in cells along compartment boundaries [52,53,54]. As will be seen in the following, the role of actomyosin in exerting tensile force on AJs is crucial in the activation of many biochemical pathways, including the Hippo pathway [55].
The AJs provide an interface between the CSK and adjacent cells via E-cadherins (E-cads) [48,56] (compared in Figure 6). E-cads on opposing membranes form homophilic bonds with their extracellular domains, and these bonds have catch-slip properties, which means that their stability increases under tension. Adjacent E-cads on the same cell can dimerize via cis bonds by “swapping” their terminal regions, and the connection can be terminated by undoing this process [57]. Both tangential and normal forces acting on the AJs are transmitted to the contiguous cell, and these cell-to-cell connections in turn mediate tissue-level contractile forces.
The intracellular domains of the E-cads act as mechanotransducers that employ adapter proteins to connect them to the cortex and the CSK. The adapters include p120-catenin, α - and β -catenin (Armadillo (Arm) in Drosophila), and Vinculin (Vinc), which link the E-cads to the actin network (compared in Figure 6), and p120-catenin binds with E-cadherin at the inner surface of the membrane and interacts with myosin VI and RhoGTPase [59]. Arm binds directly to E-cadherin [60], and this occurs prior to E-cadherin reaching the membrane [59]. At the AJs, Arm forms part of the mechanotransduction link via its binding to both E-cadherin and α -catenin [49,61].
The tangential forces arise from actomyosin filaments in the actin belt aligned parallel to the membrane, while the normal forces arise from the CSK and from branched networks produced by Arp2/3 in the cortex [62,63]. The latter can generate protrusive forces that are important for maintaining stable cell-cell contacts at the AJs. Although the cortical belt to which the AJs are attached appears as a relatively homogeneous band of an actomyosin network along the cell junctions, there is additional internal organization. At a tricellular junction (TCJ) where three cells meet, specialized proteins are involved in cell division and other processes [64]. In particular, the protein Mud is involved in division orientation control [65], and loss of Mud randomizes division orientation in the wing disc, thereby inhibiting oriented division [66]. A distinct type of puncta forms adjacent to Jub/Vinculin puncta. Zyxin is recruited to junctions under tension, and its primary function is to repair actin in conjunction with Ena/Vasp [67,68,69]. In addition, Hippo pathway proteins have been found to be concentrated near the TCJs, and thus TCJs may serve as important signaling centers.
The response of the AJs to applied forces is anisotropic [70]; those in the direction tangential to the membrane tend to disrupt E-cad binding, whereas those in the normal direction can strengthen cell–cell adhesion via mechanical feedback. The strengthening of junctions under normal force and weakening under tangential force is a mechanism that serves to maintain epithelial integrity while allowing enough plasticity for division and cell rearrangements.
The actin nucleating and branching protein complex Arp2/3 is also found at AJs. It is believed to be necessary for AJ stability and formation by creating protrusive actin networks that push the AJ complexes on opposite cells together [62,71]. It was also found that the apical Crumb/Ex/Kibra/Mer complex responds to mechanical strains. When the apical area of the Drosophila follicular epithelium was stretched, reduction of the concentration of these complexes was observed, which in turn activated Yki activity [72]. But because the dynamics of Crumb, Kibra, and Mer, as well as their interactions with Ex are not considered in the current model, we assume constant inputs from these components in response to mechanical stretches and study the mechano-effect of adherens junctions separately.

3.1. The Structure of α -Catenin

The most thoroughly-studied mechanosensitive protein in an AJ is α -catenin [73,74,75]. It can bind both F-actin and β -catenin, and under tension it can undergo an “unfurling” change in conformation which exposes cryptic binding sites for Vinculin, actin, and Jub. As shown in Figure 7, α -catenin contains three domains:
  • An N domain which can bind β -catenin and participate in the homo-dimerization of non-junctional α -catenin;
  • An M domain (central modulatory domain) which contains the Vinculin binding site (VBS);
  • An ABD domain, namely the C-terminal actin-binding domain.
A schematic diagram of these sites and the proteins that bind to them is shown in Figure 7. The M and ABD domains can undergo force-dependent conformational changes (the unfurling mentioned earlier) that can lead to force-dependent binding of various proteins. Several AJ-associated proteins can bind with α -catenin under different levels of tension [74], including the F-actin binding protein Vinculin [63] and the upstream Hippo pathway regulator Jub [76].
In Drosophila, binding of Jub to the N-terminal domain of α -catenin is tension-dependent [75,76]. The recruitment of Jub is negatively regulated by the M1 domain of α -catenin, but the binding rate does not depend on the presence of Vinculin [74]. Once recruited, junctional Jub sequesters Wts [77], preventing it from binding with the Wts-activating proteins Ex [36] and Crb [78] located apically to the AJs in the subapical region of disc cells. The effect of Jub on wing disc growth is known to be dynamic and asymmetric with regard to elevated or diminshed tension; younger discs are more sensitive to Jub underexpression, whereas older discs are more sensitive to Jub overexpression [52], and it is believed that most Wts is bound to Jub in late third-instar discs [36]. Interestingly, Jub also appears to be an essential component for AJ stability, possibly via a Rac-mediated mechanism that enhances actin association with E-cadherin complexes [79]. Other results show that Jub also interacts directly with F-actin [74,80]. This system exhibits multiple examples of feedback, as Jub mutants also cause Hippo-independent changes in myosin activity and may be involved in a negative feedback loop that downregulates tension via the activation of the cytohesin Step [81]. Finally, Jub may play a parallel role in the recruitment of Yki to the cell membrane, where it can activate myosin [82].
Although Vinculin does not affect the binding rate between Jub and the N domain of α -catenin, it does inhibit the folding process of the M1 domain of α -catenin. In addition, Vinculin can also exhibit a force-dependent enhancement of α -catenin binding to F-actin in that association of Vinculin with the ternary complex of E-cadherin, β -catenin, and α E-catenin increases the bound lifetime of individual complexes on F-actin as a function of the number of load-bearing complexes bound but only when force is directed toward the pointed end of the filament [83].
A summary of the force-dependent binding of Jub to α -catenin is given below [74]:
  • M1 limits Jub recruitment, and loss of M1 causes Jub hyper-recruitment to AJs, promoting tissue tension-independent overgrowth. Although M1 binds Vinculin, Vinculin is not responsible for this effect.
  • M1 normally limits junctional localization of Jub, while M2 and M3 normally appear to enhance junctional localization of Jub.
  • The N domain is essential for Jub recruitment to AJs, but it was also found that the α -catenin-mediated recuitment of Jub to AJs does not seem to have a linear relationship with growth regulation in the wing disc epithelium, as normal growth is compatible with normal or low concentrations of Jub at AJs.
A summary of the effect of tension sensing is shown in Figure 8, where there is a Jub-dependent mechanism ① and an unknown Jub-independent mechanism for growth regulation ②.
The authors show that stability of the cadherin-catenin complex (CCC) is controlled by the M2 domain and the α -catenin actin-binding domain (ABD) and that lower CCC levels reduce Jub recruitment.

3.2. A Three-State Model for Binding and Unfurling of the AJs

Because there are different energy wells associated with the different states of α -catenin, the response to tension involves discontinuous jumps between these states under the force-driven transitions of α -catenin. An early model to describe cell–cell adhesion was developed by Bell [84], who used an exponential function to describe the lifetime τ of molecular bonds. This takes the form
τ = τ 0 e ( E 0 γ F ) / k B T ,
where τ 0 is the reciprocal of a natural frequency, E 0 is the bond energy, F is the applied force per bond, and γ is a parameter that must be determined empirically to account for the structure of the material. This model has been used to describe protein folding and unfolding under force (reviewed in [85]) and leads to the following expression for the rate constant of protein folding and unfolding:
k ( F ) = k m k 0 e F x / k B T ,
Here k ( f ) is the rate constant under a force f, k m includes the contributions of the components of the system to the observed rates, k 0 is the intrinsic rate constant in the absence of force, F is the applied force, and x is the distance to the transition state. If the forces for transitions are known, then one can approximate the value of k * 0 k m k 0 [86] (see Appendix A.2).
We can define three distinct states of α -catenin: (1) only the N region is exposed, and all M regions are folded; (2) the N and M1 regions are exposed; and (3) the N and all the M regions are unfolded. Here, we assume that the M2 and M3 regions fold and unfold together, since current experiments have not identified a clear transition state between these two states. Based on this, we can classify the E-cad- β -cat- α -cat complex into three different states: C C C 1 , C C C 2 , and C C C 3 . This leads to the system
C C C 1 k f 1 ( f ) k u 1 ( f ) C C C 2 k f 2 ( f ) k u 2 ( f ) C C C 3 ,
where
k f i ( f ) = k f i 0 e f x f i / k B T k u i ( f ) = k u i 0 e f x u i / k B T , i = 1 , 2 .
and C C C 1 + C C C 2 + C C C 3 = 1 . Hereafter, the states are labeled by dropping the leading C C . Further details on the development of the equations can be found in Appendix A.2.
In the following computations, we hold all components in Figure 9 but α -catenin constant and vary the forces to determine the steady state levels of the three states. At zero force, C 1 was open to Jub binding, but as the force increased, α -catenin started unfolding from its initial (fully folded) state (Figure 10A), and exposure of M1 led to inhibition of Jub binding. A crossover between the decreasing C 1 state and the increasing C 2 state occurred at a force of ≈4.5 pN. A similar crossover from M1 into the M2 and M3 unfolded states occurred at ≈10 pN, which accords with the experimental result [87]. Furthermore, the time for this simple module to reach a steady state was around 2∼4 s, which also accords with the experimental result [87]. Figure 10B shows how the bound and unbound Jub depended on the force. The level of bound Jub was largest when the force was large and at its second largest when the force was small. In the range from 5 pN to 9 pN, where the M1 region was unfolding, the amount of bound Jub was lowest because M1 inhibited Jub binding.This quantifies more precisely the known result that the amount of bound Jub varies with the level of force [74,88].
We show the crossover points for the state transitions as a function of Vinculin in Figure 11, where one sees that Vinculin enhanced the C 1 -to- C 2 transition by reducing the force required. Although Vinculin does not affect the binding between Jub and α -catenin, it can inhibit the folding from C 2 to C 1 , which agrees with our expectation.

3.3. A Mechanical Model for Tissue-Level Interactions

Next, we introduce a model for a tissue that incorporates both the mechanical interactions within cells that occur via the actomyosin network at the apical membrane and the interactions between cells via the cortical belt that connects TCJs. Because the quantities involved in growth control via Yorkie are predominantly localized at or near the apical membrane [36,89], we focus on a two-dimensional section of a cell that includes the apical membrane and a small vertical component of the lateral membrane where the tight junctions and AJs are found. We furthermore assume that the cytosolic fluid is quiescent, and therefore the disc mechanics are determined by the mechanical properties of the actin belt and the cyctosolic component of the actomyosin network. It is known that on a longer time scale of days, the disc mechanics change with age, with younger third-instar wing discs held in a higher tension state, and that central compression occurs as the discs age due to rapid growth in the disc’s center [52,76]. The time scale used here is in the range of seconds to minutes (for example, Jub binds to α -catenin in seconds [90]), and the results indicate the short-term changes in the network that resulted from changes in various factors rather than the final state.
We studied a small network of cells—a central cell and six neighbors, as shown in Figure 12A—to understand how stress and imposed growth interact at the single-cell level. All cells were equal-sized hexagons initially, and we prescribed the elements in each cell as shown. The nodes received forces from the belt connections in the network and transverse forces from other nodes in their cell that reflected the actin network attached to the apical membrane. How the forces were computed is described in Appendix A.3.
A model that incorporates details of the actin network growth and its interaction with myosin and numerous control molecules such as capping proteins would be far too complex to implement here [91,92,93], and thus we describe the passive components of the network with a Kelvin–Voigt (KV) viscoelastic model. We combined this with an active element for the contractile actomyosin network as shown in Figure 12B, and we refer to the combined module as the KVA model. The form of the active component is given in Appendix A.3.
The KV model represents elastic interactions and viscous damping within and between cells, and it was chosen because it describes the stress response of Drosophila embryonic wing cells well [94], and because it is the simplest of the standard models in which the tissue relaxes to its original configuration when the stress is released. Cell-ECM interactions are not considered here, but their importance for tissue pattern formation has been investigated [95], and it was found that the stress-strain constitutive equations used for the ECM play an important role in that context.
Since tissue rheology has been studied extensively, there are good estimates of the properties of the cytosol [96], the cortex (including cortical tension) [97,98], and the effects of perturbations of actomyosin activity [99]. This level of detail in the model is necessary for several reasons. Vertex-based models have been widely used for modeling epithelial sheets [100,101,102], but they are inadequate here because they ignore active stress [103,104], since cell surface tensions—represented here by the cytosolic connections—are comparable to cell-cell interface tensions in embryonic Drosophila tissue [105], and because such models assume instantaneous relaxation of mechanical interactions, which may not hold in wing disc tissue in general [106].
To develop the governing equations, we defined two matrices to describe the connectivity in the tissue: one ( Δ N ) for the actin belt connections between nodes and another for the transverse connections that represent the actomyosin network underlying the apical membrane, which we call Δ T . Formally, these represent the Laplacian operator on the underlying graph of the network, and we introduced two matrices because different constants may appear in the node-to-node connections and the transverse connections. Under the earlier assumptions, the structure of the governing dynamical equations for the network represents a force balance given by the following momentum equation:
M d 2 x d t 2 = K N Δ N ( x l ) + K T Δ N ( x l ) + Ξ d x d t + F A ( d l d t ) + F I m p .
Here, M is the diagonal matrix whose entries are the masses of the nodes, x is the vector of the spatial positions of the nodes in R 2 , l is the rest length of the springs, x l is the vector of spring extensions, K N ( K T ) represents the diagonal matrices whose entries are the elasticity coefficients for the belt connections (transverse connections), and Ξ is the matrix of damping coefficents for the viscous response. When all nodes have the same properties, the matricies M , K N , K T , and Ξ are replaced by constants. The term F A represents the active force, and F I m p represents an externally imposed force. Details for these and other terms in Equation (4) are given in Appendix A.3.
In addition to the mechanical dynamics, we allowed for growth of the nodal and transverse connections, and this led to the equation
d l d t = F G
where F G represents the growth function. In reality, this can depend on numerous factors such as capping proteins, cofiln, and others that affect the growth of the actomyosin network [91,92], but we used the simplified form given in Appendix A.3.
As an illustration of the effect of abnormal growth of one cell under purely mechanical conditions—the biochemical network is turned off—we began with a configuration of equal hexagonal cells in which the active and passive forces were balanced. We then artificially increased the size of the central cell, which was intended to model a volume change of a cell that had taken up water due to osmotic changes, and computed the new stationary state when the system relaxed. The results are shown in Figure 13, where here and hereafter, forces applied at the exterior nodes were directed along a ray connecting the node and the center of the network. One can see that in each case, the swelling of the central cell changed the surrounding cells into irregular hexagons (i.e., hexagons of unequal sides). In Figure 13b, the contractile forces shrank the network. (We did not impose volume conservation to reflect the fact that cells can adjust their volume by controlling the water content via tension-dependent aquaporins.) Application of a stretching force of 800 pN per node in Figure 13c had a dramatic effect on the size, whereas a compressive force of 100 pN in Figure 13d had little effect.

4. The Interaction of Biochemical Signaling and Mechanics

In this section, we combine the biochemical and biomechanical models developed earlier, the result of which we call the biochemical signaling and mechanics (BSM), and we compare simulation results from the BSM with experimental observations for different conditions. We have seen that α -catenin plays a central role in transmitting signals between the mechanical and biochemical components, which in turn respond with direct or indirect feedback. We have focused on α -catenin, but it should be noted that the experimentally observed junctional complex includes other components such as E-cad and β -catenin, and the formation of the complex is regulated not only by force [48,58,107,108] but also by biochemical factors such as Vinculin [109,110,111]. Here, we ignore the formation of junctions and only study their role in regulating junctional Jub localization, which in turn regulates Wts and Yki. Using the BSM, we aim to understand the experimental observations related to regulation of the Hippo pathway by mechanical inputs and to obtain insight into how their interactions can control cell and tissue growth. Possible extensions of the current model are discussed in Section 5.
Based on the studies of α -catenin described in Section 3, we assume it has three possible states: C 1 , C 2 , and C 3 , which correspond to the degree to which it is unfolded. In the simulations, we assumed that the total amount of α -catenin was maintained at a constant level, and in light of the fact that the expression level of E-cad in epithelial cells is approximately 26–46 molecules/ μ m 2 [112], we assumed that the membrane-localized α -catenin was maintained at a similar level. In addition, as most of the E-cad was localized at AJs in a region which was around 1 μ m below the apical membrane [113], and the averaged edge length of the cell edges in the Drosophila wing was around 1 μ m, we assumed that the total amount of the different states of α -catenin was about 26–46 molecules per cell edge. In our model, the mass is carried by vertices as shown in Figure 12, and therefore the total number of α -catenins was set to 26–46 molecules per node. Furthermore, despite the fact that the level of membrane-localized α -catenin can change dramatically at different stages of cell development [113], we will see that the qualitative model predictions will not change with changes in α -catenin.
In the diagram shown in the left panel of Figure 14, we illustrate qualitatively the Yki activity level as a function of the tension forces acting on an α -catenin molecule. This shape can be understood by comparison with the level of bound Jub in Figure 10B, because when Jub is low, Wts will be high, and Yki will be low. The dip in Jub at intermediate forces stems from the fact that when M1 is open, it inhibits binding of Jub to the N-domain of α -catenin (see Section 3), and this inhibition is suppressed at sufficiently high tensions. Thus, the system can be in one of three states denoted by 1, 2, and 3, in which the relative level of Yki first decreases from stage 1 to 2 and then increases from stage 2 to 3. While this qualitative description is robust, details may change as parameters are changed. For example, the level of Vinculin controls the location of the transition from stage 1 to 2, as shown in Figure 11. Increasing Vinculin decreases the transition force needed for a transition from C 1 to C 2 , which will shorten the 1 stage in Figure 14 (Left) as a result. Furthermore, cells may experience different cytoskeletal tensions and possess different levels of α -catenin in different developmental stages [110,113], and their states may progress along a non-monotonic path in the x direction of Figure 14. This will be elaborated upon later.
In addition to analyzing the BSM under different cytoskeletal forces, we also explored how the system responds to changes in α -catenin. For this purpose, we constructed a growth model for cells shown in the right panel of Figure 14 and assumed that the cell area grew following a Yki-dependent equation described in Appendix A.3. As an initial condition, we assumed that the volume of the center cell enlarged from the region enclosed by black dashed lines to that enclosed by red dashed lines due to, for example, changes in osmotic conditions. Tension and compression forces can also be applied to the boundary nodes, as in the pure mechanical example shown in Figure 13. In any case, one must now solve the system comprising Equations (4) and (5).
According to the three-state model of α -catenin described earlier, increasing a tension force from zero should first produce a decrease in Yki, as shown in Figure 14, followed by an increase at larger forces. However, experimental observations show that the Yki level increased with increasing cytoskeletal tension [76]. Based on the model, we suggest that this discrepancy between the molecular- and cell-level observations may be due to the magnitude of the imposed force. Indeed, when a single α -catenin is bound with E-cad and β -catenin, a force of around 5 pN is needed to unfold the M1 region (although when vinculin exists, this value may decrease) [87], which converts the α -catenin to state II. Given that the tension force generated by a single Myo-II molecule is around 3–5 pN [110], it is possible that in cell-level experiments, α -catenin unfolds the M1 region directly without passing through state I in the left panel of Figure 14. To test this hypothesis, we plotted the Yki activity as a function of the cytoskeletal tension and the level of membrane-localized α -catenin in Figure 15. When the α -catenin level was low (20–50 molecules/node, which was the normal range we set in our system), the Yki level increased in response to increasing force (right lower panel of Figure 15). However, if the α -catenin level increased to >50 molecules/node, then the transitions were those predicted by the molecular-level mechanism, as shown in the upper right panel of Figure 15 [74,75,88]. This is due to the fact that as the level of α -catenin increases with the same magnitude of tension applied to each node, the force experienced by individual α -catenin molecules decreases. Further insight can be obtained from a vertical slice in the heat map shown in Figure 15. When the magnitude of the applied tension was fixed at, for example, 600 pN/node, Yki would first increase and then decrease as α -catenin increased. Since α -catenin modulates the Hippo pathway through membrane localization of Jub, which is a negative regulator of Warts, an elevated level of α -catenin promotes Yki activity at intermediate levels. It was also observed that a loss of α -catenin can cause a decrease in the Yki level [74,76], which leads to tissue undergrowth. Beginning with α -catenin at 60–70 and the tension fixed at 300 pN/node, one sees in Figure 15 that Yki is essentially monotone decreasing as α -catenin decreases, which comports with the observations. It was observed that the E-cad level increased as the wing disc cells aged [113], and the localized α -catenin level may also increase. If the tension does not change proportionally the Yki activity will be inhibited by the increase of α -catenin and possibly lead to cell death.
We further investigated the effect of the α -catenin level on Yki activity with imposed compression and tension forces. For the results shown in Figure 16, we varied the amount of α -catenin in the center cell while keeping the α -catenin in the boundary cells constant. One can see there that by fixing the boundary tension, increasing α -catenin in the center cell can induce higher Yki activitiy there, which again agrees with the experimental observation described earlier. One can also see that at low α -catenin levels, Yki increased more slowly at the center under tension than compression, but the difference disappeared at sufficiently large α -catenin levels. Interestingly, increasing α -catenin in the center cell had a negative effect on the Yki level in the boundary cells under both tension and compression. We could not solve the governing system analytically, but this effect was caused by the mechanical feedback among cells in the tissue, which Yki would transduce into an effect on the growth. This observation provides a possible explanation for the observation of cell-autonomous decreases and non-cell-autonomous increases in Yki activity in tissues with E-cad or α -catenin knockdown [114].
Next, we investigated the effect of α -catenin on the discrepancies in the Yki level and cell growth in the center and boundary cells. As shown in Figure 17, with a higher α -catenin level, the difference in the Yki level and the growth between the center and boundary cells increased. The effect of the magnitude of the applied environmental force only had a significant effect when the α -catenin level was relatively low. To understand this result, we should first notice that a lower boundary tension will decrease the cytoskeletal tension, especially in this case, where the boundary tension varied from 800 pN to 100 pN. Furthermore, as shown in Figure 15, if α -catenin is in the normal range (from 20 to 50), then the system will generate an increasing Yki level as the cytoskeletal tension increases, but if the α -catenin amount is higher (>50), then the system may produce the same high Yki level at lower cytoskeleton tension as a system with high cytoskeletal tension.
We also analyzed the model under a similar setting with and without cell growth, as shown in Figure 18. After the same evolution time for each case, with cell growth, the Yki level was always lower, although slightly, compared with that without cell growth for both magnitudes of the applied environmental forces. We think this was caused by the unbalanced growth; the center cell had a higher growth rate than the boundary cell, which could also be reflected by the Yki level. This unbalanced growth will increase the compression between cells in the tissue and, as a result, inhibit the Yki activity compared with the no-growth cases. This observation can be interpreted as the feedback of the biochemical module on the mechanical changes. Since the growth of the center cell caused expansion of the cell area, it compressed the surrounding cells. In return, the Yki level in the center cell decreased as a responses to mechanical cues to avoid tissue overgrowth and overcompression. Such feedback has been widely explored both theoretically and experimentally [7,115,116,117,118,119], and our model provides a platform for further investigation of this mechanism. However, as shown in Figure 18, this effect was not all that significant in the model results, and this was caused by the relatively simple description of growth. The real growth process is more complicated, with dynamics of the cell size, actin, etc. such as that shown in [120]. A more comprehensive model with the incorporation of a detailed growth model could amplify this feedback effect and provide a better understanding of the details.
As discussed earlier, α -catenin plays a central role in bridging the mechanical and biochemical signals in growth control by regulating the membrane localization of Jub [76,80]. As shown in Figure 10, we used a single-molecule model to show that under higher cytoskeletal tension, Jub accumulated at the AJ, which accords with the experimental observations [76,121].
We could also analyze the difference in Yki level and growth between the center and boundary cells as functions of the α -catenin level when different Jub expressions were imposed (Figure 19). This setting could be used to mimic the case where a patch of Jub overexpressed or underexpressed cells surrounded by WT cells. As shown in Figure 19, by fixing the environmental force at 800 pN, we observed that the Yki level and growth effect both increased as the Jub production rate increased, which accords with the experimental observations that increased Jub recruitment to α -catenin is associated with increased Yki activity and wing growth [75].

5. Conclusions

5.1. Overview

Studies on vertebrate systems showed that signal transduction and growth control pathways are highly conserved across species, and thus much of what is learned about Drosophila applies in higher organisms. Tissue overgrowth in Drosophila is similar to tumor growth in mammalian systems, and a better understanding of how molecular signals and mechanics interact in development will shed light on how these signaling processes interact in the microenvironment of a tumor to affect its growth [122,123]. The ubiquity of the pathways across species and the balances of signals between them within an organ highlight the need for mathematical models to understand these complex systems.
Growth control in the Drosophila disc involves both disc-level and extrinsic control [124], and since Yki-dependent growth requires dTOR activity, the extrinsic control structure may depend on a hypothetical mechanism in which dTOR affects the state of nuclear Yki [125]. The dTOR and Hippo pathways may act in parallel in that Hippo assesses local growth suitability, and dTOR ensures that organismal factors such as nutrition are sufficient to support the increased growth in the early stage of the cell cycle. In [126], we suggested a description called the mechanical feedback model (MFM), in which mechanics and signaling are integrated for describing disc growth, but while the MFM incorporates part of the Hippo pathway, a description of growth also requires the integration of extrinsic growth controls transduced via the dTOR pathway [127].
Mechanical effects due to an influx or efflux of water may also be important. As we showed in the example in Section 3.3, enlargement of a cell can have a significant effect on the cells connected to it, and this example suggests that osmotic effects can have an important impact, as noted in several recent reviews [128,129,130]. In particular, recent research has shown that volume and shape control can play a significant role in cancer cell migration, and this involves a variety of control factors [131,132,133,134]. Ion channels and aquaporins are central components of cell volume regulation, and water flow driven by osmotic gradients generated by ion transport contributes to the driving force for cell migration. In a simple experimental example, it was shown that tumor cells confined in a narrow channel can translocate by establishing an asymmetric end-to-end distribution of Na+/H+ pumps and aquaporins, which creates a net inflow of water and ions at the cell leading edge and a net outflow of water and ions at the trailing edge [135]. Moreover, it has been reported that metastatic cancer cells have higher expression of these transport proteins than nonmetastatic cancer cells [136].

5.2. Current Results

Herein, we developed a model for tissue development in the Drosophila wing disc that integrates an earlier model of the Hippo pathway with a model of the tissue mechanics. In light of the complexity of the actomyosin network and the cell-cell connections that govern the mechanics, it was necessary to use the high-level Kelvin–Voigt description of the cortical and cytoskeletal dynamics and a simplified model of cell-cell connections. Nonetheless, a number of significant insights have emerged and a number of important results have been obtained despite the simplifications.
In Section 3, we developed a model for the unfolding of α -catenin under increasing force levels, and in particular, we showed how the transition between the C 1 and C 2 states depends on the level of Vinculin. This and other results reported herein involving the unfolding process provide a starting point for further investigations of the dynamical response of the cell-to-cell connection machinery based on E-cads, α -catenin, and other components. The universality of the components involved in cell–cell adhesion suggests that what is learned from Drosophila will be applicable to other systems and may shed light on alterations in the interactions that arise in cancer.
In Section 4, we performed numerical experiments that validate our model by comparison with experimental results and enabled us to make testable predictions. It has been widely reported that α -catenin is a critical component that connects the mechanical cues with the Hippo pathway, but the detailed mechanism is not well understood [75]. Recently, with the analysis of phenotypic defects resulting from the differential reduction of gene function [74,88], it was found that different binding regions of α -catenin have different effects on the binding between α -catenin and Jub. Furthermore, the folding and unfolding behavior of α -catenin was found to be force-dependent [87,137]. Based on these findings, we constructed a model that describes the conversion between different states of α -catenin and employed Bell’s model to describe the force-dependent transition rates. The simulation results were compared with experimental observations, which showed that the model successfully reproduced the force-regulated folding and unfolding of α -catenin.
Next, we integrated this model with the model for the Hippo pathway, and from this, were able to reproduce the effect of elevated α -catenin and the cytoskeletal tension on Yki activity. In addition, we can explain a critical discord between the molecule-level and cell-level observations: Although the molecule-level experiments predicted three distinct stages of force-dependent Yki activities, in the cell-level experiments, only the latter two were observed [76]. We suggest that this discrepancy may be due to the relatively high magnitude of forces exerted by Myo-II on single α -catenin molecules, which implies rare observation of the first stage. We can also speculate, based on this analysis, that in older cells, a higher E-cad level may reduce Yki activity, which could relate to cell death.
Finally, integration of the mechanical and biochemical model with a model for growth enabled us to reproduce many growth-related observations, including the role of α -catenin level-dependent tissue growth and the effect of Jub expression on growth. The results show that the integrated response of the mechanical, biochemical, and growth modules to variations in the α -catenin level across a tissue emphasized its role as a critical pivot in tissue growth regulation [7,115,116,117,118,119].
In summary, the interaction between the Hippo pathway and the mechanical forces can lead to control of growth in the Drosophila wing disc under various conditions, although a more detailed model is needed to investigate the mechanical effects further.

5.3. Future Directions

In light of the fact that this is the first attempt known to us to model this complex system, we have indicated a number of potential directions for future investigations.
  • As stated earlier, in the current model, we ignored the feedback loop in which Yki activity activates the expression of upstream regulators, such as Four-jointed, Ex, Mer, and Kibra [138], and a similar feedback loop exists in that of mammalian cells [139]. On a longer time scale, this feedback may play a significant role in controlling growth.
  • A more detailed model of the actomyosin dynamics is needed to more accurately describe the tissue mechanics. While much is known about individual components, a cell-level model that describes the details of the actin belt dynamics and its interaction with the E-cad- α -catenin system is not at hand. For example, a more detailed model would incorporate spectrin, a contractile protein that forms a membrane-attached skeleton beneath the plasma membrane by crosslinking short F-actin and binding-integral membrane proteins [140]. Spectrin is required for the formation of epithelia, and unlike other regulators such as Crumbs and Merlin, it regulates Hippo signaling by modulating cortical actomyosin activity through non-muscle myosin II [141].
  • A more detailed and dynamic model of the mechanics of cell–cell interactions is needed, including the binding of E-cadherin between cells and α -catenin binding and production or mobilization, is needed. This is particularly important for understanding relative movement between cells and T1 transitions.
  • Because the copy numbers of signaling molecules and other key components are frequently low, stochastic effects should be considered with a view toward understanding if and how the network structure plays a role in adapting to noisy signals. Noise can affect the precision of gene expression in simple networks [142], and key components of the Dpp pathway, for example, are present at nanomolar concentrations [143] in the disc, yet the disc patterning and size are remarkably reproducible. Though it has been argued that even such low concentrations are sufficient to mitigate stochastic noise [144], the network structure may play a role in mitigation. How the functions of signal transduction networks and mechanical regulation are maintained in the presence of fluctuations is still a major question in cellular biology.
  • The current model involves regulation of both the Hippo pathway and the mechanical pathway, and in both components many of the parameters used in the model are unknown or difficult to measure in experiments. Therefore we selected the parameters within biologically-meaningful ranges, but did not do a detailed sensitivity analysis. In the previous version of the biochemical model [9], of which the current biochemical model is a small variation, both local and global sensitivity analysis was performed to analyze the impact of the parameters in the results, and it was found that a small number (6–8) of parameters were very important, and the majority were much less important. However, the Hippo signaling network in both cases has a top-down structure and feedback loops are not considered, which may affect the results. For instance, it was found earlier that the variations of parameters close to the downstream output – the cytosolic Yki concentration – have a greater effect on the variations of the output than do the upstream parameters.
    A similar sensitivity issue applies to the mechanical component of the model, and when it has been done, the sensitivity analysis for the combined model has to be done. This has the potential to identify the key steps in the combined model and lead to a reduced model that can be used in tissue-level computations, but this is a major project in its own right, in part because how the parameters affect the qualitative responses, for instance, the Yki profile as a function of Fat production as shown in Figure 3, is difficult to analyze and required a method that differs from the traditional global and local sensitivity analysis.

Author Contributions

J.G. and T.Z. contributed equally to this work. J.G. and T.Z. performed the numerical simulations. J.G., T.Z. and H.G.O. analyzed the data and wrote the article. All authors have read and agreed to the published version of the manuscript.

Funding

J.G. is funded by start-up research funds from the University of California, Riverside. This work was supported in part by NSF Grant DMS 1311974 and by NIH Grant 54-CA-210190 to HGO. HGO would like to thank the Isaac Newton Institute for Mathematical Sciences, Cambridge, for support and hospitality during the programme Collective Behaviour where work on this paper was undertaken. This work was supported by EPSRC grant no EP/R014604/1. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the National Science Foundation or of their current employers, nor do they necessarily represent the official views of the National Institutes of Health.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. The Biochemical Model

In the tissue model, all of the membrane species are located on the six nodes of a hexagon, and the concentration of each species in the cytosol is assumed to be uniform. In the reactions given below, the cell index is i, the index j represents the index of one of the six sides, the α values are the production rates, and the β values are the degradation rates.
Table A1. Parameters and their values used in the model.
Table A1. Parameters and their values used in the model.
ProteinNotation
Membrane-Bound Species
Fat [ F ] m , j i
Dacsous(Ds) [ S ] m , j i
Fat-Ds complex [ F S ¯ ] m , j i
Dachs [ A ] m , j i
Warts-Dachs complex [ W A ¯ ] m , j i
Expanded [ E r ] m , j i
E-Cad- β -cat- α -cat complex (N, M1, M2+M3) [ C 1 ] m , j i , [ C 2 ] m , j i , [ C 3 ] m , j i
E-Cad- β -cat- α -cat-Vinculin complex [ C 2 V ¯ ] m , j i , [ C 3 V ¯ ] m , j i
E-Cad- β -cat- α -cat-Jub complex [ C 1 B ¯ ] m , j i , [ C 2 B ¯ ] m , j i , [ C 3 B ¯ ]
E-Cad- β -cat- α -cat-Jub-Vinculin complex [ C 2 B V ¯ ] m , j i , [ C 3 B V ¯ ] m , j i
E-Cad- β -cat- α -cat-Jub-Warts complex [ C 1 B W ¯ ] m , j i , [ C 2 B W ¯ ] m , j i , [ C 3 B W ¯ ] m , j i
E-Cad- β -cat- α -cat-Jub-Vinculin-Warts complex [ C 2 B V W ¯ ] m , j i , [ C 3 B V W ¯ ] m , j i
Cytosolic Species
Fat [ F ] c i
Ds [ S ] c i
Dachs [ A ] c i
Warts [ W ] c i
Yki [ Y ] c i
Four-jointed [ J ] c i
Jub [ B ] c i
Vinculin [ V ] c i

Appendix A.1.1. The Module for Membrane-Bound Species

Fat and Dacsous

d [ F ] m , j i d t = k F + [ F ] c i k F [ F ] m , j i k F S + ( [ J ] c i ) [ F ] m , j i n ( i ) , m ( i , j ) k S F + ( [ J ] c n ( i ) ) [ S ] m , m ( i , j ) n ( i ) + k F S n ( i ) , m ( i , j ) [ F S ¯ ] m , m ( i , j ) n ( i ) , d [ S ] m , j i d t = k S + [ S ] c i k S [ S ] m , j i k S F + ( [ J ] c i ) [ S ] m , j i n ( i ) , m ( i , j ) k F S + ( [ J ] c n ( i ) ) [ F ] m , m ( i , j ) n ( i ) + k F S n ( i ) , m ( i , j ) [ F S ¯ ] m , m ( i , j ) i , d [ F S ¯ ] m , j i d t = k S F + ( [ J ] c i ) [ S ] m , j i n ( i ) , m ( i , j ) k F S + ( [ J ] c n ( i ) ) [ F ] m , m ( i , j ) n ( i ) k F S [ F S ¯ ] m , j i ,
where
k S F + ( [ J ] c i ) = k F S + k F S + [ J ] c i , k F S + ( [ J ] c i ) = k F S + [ J ] c i k F S + [ J ] c i .

Dachs

d [ A ] m , j i d t = k A + ( [ F ] m , j i , n ( i ) , m ( i , j ) [ F S ¯ ] m , m ( i , j ) n ( i ) ) [ A ] c i k A [ A ] m , j i k W A + [ W ] c i [ A ] m i + k W A [ W A ¯ ] m i ,
where
k A + ( [ F ] m , j i , n ( i ) , m ( i , j ) [ F S ¯ ] m , m ( i , j ) n ( i ) ) = k A + ( 1 + ( [ F ] m , j i K 1 ) h 1 ) ( 1 + ( n ( i ) , m ( i , j ) [ F S ¯ ] m , m ( i , j ) n ( i ) K 2 ) h 2 )
d [ W A ¯ ] m i d t = k W A + [ W ] c i [ A ] m i k W A [ W A ¯ ] m i ,

Riq

d [ R ] m , j i d t = k R + ( [ S ] m , j i , n ( i ) , m ( i , j ) [ F S ¯ ] m , m ( i , j ) n ( i ) ) [ R ] c i k R [ R ] m , j i k W R + [ W ] c i [ R ] m i + k W R [ W R ¯ ] m i ,
where
k R + ( [ S ] m , j i , n ( i ) , m ( i , j ) [ F S ¯ ] m , m ( i , j ) n ( i ) ) = k R + ( [ S ] m , j i ) h 3 ( n ( i ) , m ( i , j ) [ F S ¯ ] m , m ( i , j ) n ( i ) ) h 4 ( 1 + ( [ S ] m , j i K 3 ) h 3 ) ( 1 + ( n ( i ) , m ( i , j ) [ F S ¯ ] m , m ( i , j ) n ( i ) K 4 ) h 4 )
d [ W R ¯ ] m i d t = k W R + [ W ] c i [ R ] m i k W R [ W R ¯ ] m i ,

Expanded and Expanded Binding with Warts

d [ E r ] m , j i d t = k E W + [ E r ] m , j i [ W ] c + k E W [ E W ¯ ] m , j i d [ E W ¯ ] m , j i d t = k E W + [ E r ] m , j i [ W ] c k E W [ E W ¯ ] m , j i ,

E-Cad-β-Cat-α-Cat Complex

d [ C 1 ] m , j i d t = k u 1 ( f ) [ C 1 ] m , j i + k f 1 ( f ) [ C 2 ] m , j i k 1 + [ C 1 ] m , j i [ B ] c + k 1 [ C 1 B ¯ ] m , j i d [ C 2 ] m , j i d t = k u 1 ( f ) [ C 1 ] m , j i k f 1 ( f ) [ C 2 ] m , j i k u 2 ( f ) [ C 2 ] m , j i + k f 2 ( f ) [ C 3 ] m , j i k C V + [ C 2 ] m , j i [ V ] c + k C V [ C 2 V ¯ ] m , j i k 2 + [ C 2 ] m , j i [ B ] c + k 2 [ C 2 B ¯ ] m , j i d [ C 3 ] m , j i d t = k u 2 ( f ) [ C 2 ] m , j i k f 2 ( f ) [ C 3 ] m , j i k C V + [ C 3 ] m , j i [ V ] c + k C V [ C 3 V ¯ ] m , j i k 3 + [ C 3 ] m , j i [ B ] c + k 3 [ C 3 B ¯ ] m , j i
where f is the force and
k f 1 ( f ) = k f 1 0 e f x f 1 / ( k B T ) , k u 1 ( f ) = k u 1 0 e f x u 1 / ( k B T ) , k f 2 ( f ) = k f 2 0 e f x f 2 / ( k B T ) , k u 2 ( f ) = k u 2 0 e f x u 2 / ( k B T ) .

E-Cad-β-Cat-α-Cat-Vinculin Complex

d [ C 2 V ¯ ] m , j i d t = k u 2 ( f ) [ C 2 V ¯ ] m , j i + k f 2 ( f ) [ C 3 V ¯ ] m , j i + k C V + [ C 2 ] m , j i [ V ] c k C V [ C 2 V ¯ ] m , j i k 2 + [ C 2 V ¯ ] m , j i [ B ] c + k 2 [ C 2 B V ¯ ] m , j i d [ C 3 V ¯ ] m , j i d t = k u 2 ( f ) [ C 2 V ¯ ] m , j i k f 2 ( f ) [ C 3 V ¯ ] m , j i + k C V + [ C 3 ] m , j i [ V ] c k C V [ C 3 V ¯ ] m , j i k 3 + [ C 3 V ¯ ] m , j i [ B ] c + k 3 [ C 3 B V ¯ ] m , j i

E-Cad-β-Cat-α-Cat-Jub Complex

d [ C 1 B ¯ ] m , j i d t = k u 1 ( f ) [ C 1 B ¯ ] m , j i + k f 1 ( f ) [ C 2 B ¯ ] m , j i + k 1 + [ C 1 ] m , j i [ B ] c k 1 [ C 1 B ¯ ] m , j i d [ C 2 B ¯ ] m , j i d t = k u 1 ( f ) [ C 1 B ¯ ] m , j i k f 1 ( f ) [ C 2 B ¯ ] m , j i k u 2 ( f ) [ C 2 B ¯ ] m , j i + k f 2 ( f ) [ C 3 B ¯ ] m , j i + k 2 + [ C 2 ] m , j i [ B ] c k 2 [ C 2 B ¯ ] m , j i k C V + [ C 2 B ¯ ] m , j i [ V ] c + k C V [ C 2 B V ¯ ] m , j i d [ C 3 B ¯ ] m , j i d t = k u 2 ( f ) [ C 2 B ¯ ] m , j i k f 2 ( f ) [ C 3 B ¯ ] m , j i + k 3 + [ C 3 ] m , j i [ B ] c k 3 [ C 3 B ¯ ] m , j i k C V + [ C 3 B ¯ ] m , j i [ V ] c + k C V [ C 3 B V ¯ ] m , j i

E-Cad-β-Cat-α-Cat-Jub-Vinculin Complex

d [ C 2 B V ¯ ] m , j i d t = k u 2 ( f ) [ C 2 B V ¯ ] m , j i + k f 2 ( f ) [ C 3 B V ¯ ] m , j i + k 2 + [ C 2 V ¯ ] m , j i [ B ] c k 2 [ C 2 B V ¯ ] m , j i + k C V + [ C 2 B ¯ ] m , j i [ V ] c k C V [ C 2 B V ¯ ] m , j i d [ C 3 B V ¯ ] m , j i d t = k u 2 ( f ) [ C 2 B V ¯ ] m , j i k f 2 ( f ) [ C 3 B V ¯ ] m , j i + k 3 + [ C 3 V ¯ ] m , j i [ B ] c k 3 [ C 3 B V ¯ ] m , j i + k C V + [ C 3 B ¯ ] m , j i [ V ] c k C V [ C 3 B V ¯ ] m , j i

E-Cad-β-Cat-α-Cat-Jub-Warts Complex

d [ C 1 B W ¯ ] m , j i d t = k u 1 ( f ) [ C 1 B W ¯ ] m , j i + k f 1 ( f ) [ C 2 B W ¯ ] m , j i + k W B + [ W ] c [ C 1 B ¯ ] m , j i k W B [ C 1 B W ¯ ] m , j i d [ C 2 B W ¯ ] m , j i d t = k u 1 ( f ) [ C 1 B W ¯ ] m , j i k f 1 ( f ) [ C 2 B W ¯ ] m , j i k u 2 ( f ) [ C 2 B W ¯ ] m , j i + k f 2 ( f ) [ C 3 B W ¯ ] m , j i + k W B + [ W ] c [ C 2 B ¯ ] m , j i k W B [ C 2 B W ¯ ] m , j i k C V + [ C 2 B W ¯ ] m , j i [ V ] c + k C V [ C 2 B W V ¯ ] m , j i d [ C 3 B W ¯ ] m , j i d t = k u 2 ( f ) [ C 2 B W ¯ ] m , j i k f 2 ( f ) [ C 3 B W ¯ ] m , j i + k W B + [ W ] c [ C 3 B ¯ ] m , j i k W B [ C 3 B W ¯ ] m , j i k C V + [ C 3 B W ¯ ] m , j i [ V ] c + k C V [ C 3 B W V ¯ ] m , j i

E-Cad-β-Cat-α-Cat-Jub-Vinculin-Warts Complex

d [ C 2 B V W ¯ ] m , j i d t = k u 2 ( f ) [ C 2 B V W ¯ ] m , j i + k f 2 ( f ) [ C 3 B V W ¯ ] m , j i + k W B + [ W ] c [ C 2 B V ¯ ] m , j i k W B [ C 2 B V W ¯ ] m , j i + k C V + [ C 2 B W ¯ ] m , j i [ V ] c k C V [ C 2 B W V ¯ ] m , j i d [ C 3 B V W ¯ ] m , j i d t = k u 2 ( f ) [ C 2 B V W ¯ ] m , j i k f 2 ( f ) [ C 3 B V W ¯ ] m , j i + k W B + [ W ] c [ C 3 B V ¯ ] m , j i k W B [ C 3 B V W ¯ ] m , j i + k C V + [ C 3 B W ¯ ] m , j i [ V ] c k C V [ C 3 B W V ¯ ] m , j i

Appendix A.1.2. The Cytosolic Module

Fat and Dacsous

d [ F ] c i d t = α F β F [ F ] c i N m k F + [ F ] c i + k F j = 1 N m [ F ] m , j i , d [ S ] c i d t = α S β S [ S ] c i N m k S + [ S ] c i + k S j = 1 N m [ S ] m , j i ,

Dachs

d [ A ] c i d t = α A β A [ A ] c i [ A ] c i j = 1 N m k A + ( [ F ] m , j i , [ F S ¯ ] m , j i ) + k A j = 1 N m [ A ] m , j i ,

Riq

d [ R ] c i d t = α R β R [ R ] c i [ R ] c i j = 1 N m k A + ( [ S ] m , j i , [ F S ¯ ] m , j i ) + k A j = 1 N m [ R ] m , j i ,

Warts

d [ W ] c i d t = k W R + [ W ] c i j = 1 N m [ R ] m i + k W R j = 1 N m [ W R ¯ ] m i k W A + [ W ] c i j = 1 N m [ A ] m i + k W R j = 1 N m [ W A ¯ ] m i k E W + [ W ] c i j = 1 N m [ E r ] m , j i + k E W j = 1 N m [ E W ¯ ] m , j i k W B + [ W ] c i j = 1 N m ( [ C 1 B ¯ ] m , j i + [ C 2 B ¯ ] m , j i + [ C 3 B ¯ ] m , j i + [ C 2 B V ¯ ] m , j i + [ C 3 B V ¯ ] m , j i ) + k W B j = 1 N m ( [ C 1 B W ¯ ] m , j i + [ C 2 B W ¯ ] m , j i + [ C 3 B W ¯ ] m , j i + [ C 2 B V W ¯ ] m , j i + [ C 3 B V W ¯ ] m , j i )

Yorki

d [ Y ] c i d t = α Y β Y [ Y ] c i [ Y ] c i j = 1 N m k E W Y + [ E W ¯ ] m , j i k E W Y + [ E W ¯ ] m , j i [ Y ] c i j = 1 N m k E r Y + [ E r ] m , j i k E r Y + [ E r ] m , j i ,

Four-Jointed

This is a constant level in our system.

Jub

d [ B ] c i d t = α B β B [ B ] c i [ B ] c i j = 1 N m ( k 1 + [ C 1 ] m , j i + k 2 + [ C 2 ] m , j i + k 3 + [ C 3 ] m , j i + k 2 + [ C 2 V ¯ ] m , j i + k 3 + [ C 3 V ¯ ] m , j i ) + j = 1 N m ( k 1 [ C 1 B ¯ ] m , j i + k 2 [ C 2 B ¯ ] m , j i + k 3 [ C 3 B ¯ ] m , j i + k 2 [ C 2 B V ¯ ] m , j i + k 3 [ C 3 B V ¯ ] m , j i ) .

Vinculin

d [ V ] c i d t = α V β V [ V ] c i k C V + [ V ] c i j = 1 N m ( [ C 2 ] m , j i + [ C 3 ] m , j i + [ C 2 B ¯ ] m , j i + [ C 3 B ¯ ] m , j i + [ C 2 B W ¯ ] m , j i + [ C 3 B W ¯ ] m , j i ) + k C V j = 1 N m ( [ C 2 V ¯ ] m , j i + [ C 3 V ¯ ] m , j i + [ C 2 B V ¯ ] m , j i + [ C 3 B V ¯ ] m , j i + [ C 2 B V W ¯ ] m , j i + [ C 3 B V W ¯ ] m , j i )
Table A2. Parameters and their values used in the model.
Table A2. Parameters and their values used in the model.
NameDescriptionValuesUnitsReference
Basic parameters
Parameters in C 1 , C 2 , and C 3 folding and unfolding module
k f 1 0 Transition rate of C 2 C 1 with no load3.30 × 10 4 -[87]
x f 1 Transition distance from C 2 to C 1 8.20 nm[87]
k u 1 0 Transition rate of C 1 C 2 with no load7.00 × 10 4 s 1 [87]
x u 1 Transition distance from C 1 to C 2 7.40nm[87]
k f 2 0 Transition rate of C 3 C 2 with no load9.40 × 10 5 -[87]
x f 2 Transition distance from C 3 to C 2 13 nm[87]
k u 2 0 Transition rate of C 2 C 3 with no load3.00 × 10 7 s 1 [87]
x u 2 Transition distance from C 2 to C 3 5.60nm[87]
Production rates
α F Fat production rate0.4 μ M·min 1 [9]
α S Dacsous production rate0.2 μ M·min 1 [9]
α A Dachs production rate0.25 μ M·min 1 [9]
α R Riq production rate0.1 μ M·min 1 [9]
α Y Yki production rate0.5 μ M·min 1 [9]
α B Jub production rate0.25 μ M·min 1 This paper
α V Vinculin production rate0.1 μ M·min 1 This paper
Parameters in Jub binding with CCC complex
k 1 + , k 2 + , k 3 + Binding rate of Jub and α -catenin10 μ M·min 1 This paper
k 1 Dissociation rate of Jub and C 1 4.5 μ M·min 1 This paper
k 2 Dissociation rate of Jub and C 2 500 μ M·min 1 This paper
k 3 Dissociation rate of Jub and C 3 4 μ M·min 1 This paper
Parameters in Vh binding with M1 region
k C V + Binding rate of Vh and M1 region(2.88 ± 0.42) × 10 2 μ M 1 · min 1 [145,146]
k C V Dissociation rate of Vh and M1 region ( 4.8 ± 3 ) × 10 1 min 1 [145,146]
Degradation rates
β F , β S , β A , β R Ft, Ds, Dachs, Riq degradation rate0.1min 1 [9]
β Y , β B Yki and Jub degradation rate0.01min 1 This paper
β V Vinculin degradation rate0.05min 1 This paper
Parameters in Fat-Ds binding
k F + , k S + Forward mem localization rate ofcytosolic Ft and Ds4min 1 [9]
k F , k S Backward mem localization rate ofcytosolic Ft and Ds3min 1 [9]
k F S + Binding rate of Ft-Ds4 μ M·min 1 [9]
k F S Unbinding rate of Ft-Ds3min 1 [9]
Parameters in k A + ( · ) function
k A + Max localization rate of Dachs to mem0.21min 1 This paper
k A Unbinding rate of Dachs from mem0.2min 1 This paper
k R + Max localization rate of Riq to mem0.3min 1 This paper
k R Unbinding rate of Riq from mem0.5min 1 This paper
K 1 , K 2 , h 1 , h 2 Scaling parameters1, 2, 1, 1-This paper
K 3 , K 4 , h 3 , h 4 Scaling parameters1, 2, 1, 1-This paper
Miscellaneous
k W A + Binding rates between Dachs and Warts0.1min 1 This paper
k W A Release rates from Dachs and Warts0.05min 1 This paper
k E r Y + Maximum inhibition rate of Expanded on Yki0.1min 1 This paper
k E r Y Scaling parameter0.05min 1 This paper
k E W + Binding rate between Expanded and Warts0.1min 1 This paper
k E W Unbinding rate between Expanded and Warts0.05min 1 This paper
k W B + Binding rate between Wts and mem localized Jub10min 1 This paper
k W B Unbinding rate between Wts and mem localized Jub1min 1 This paper
k E W Y + Maximum inhibition rate of Expanded-Warts on Yki1min 1 This paper
k E W Y Scaling parameter0.05min 1 This paper

Appendix A.2. The Three States of α-Catenin

The binding strength between the different states of CCC and J u b is different. We denote the on and off binding rates that correspond to different states of the CCC as k on i and k off i , respectively, where i = 1 , 2 , 3 . For different states of CCC, based on the effect of different M regions on the J u b and α -catenin binding, we can assume
k 1 + = k 2 + = k 3 + ,
and the off binding rates have the following relation:
k 2 > k 1 > ( = ) k 3 .
Moreover, once V h (Vinculin) binds with the M1 region, the first folded process in Equation (2) cannot happen again. Therefore, we can use another binding complex to represent this behavior:
C 2 + V h k C V k C V + C 2 - V h ,
where k C V + = ( 4.8 ± 0.7 ) × 10 6 M 1 s 1 and k C V = ( 8 ± 5 ) × 10 3 s 1 and hence K D = ( 2 ± 1 ) × 10 9 M [145]. Furthermore, ref. [146] gives a more comprehensive (under a different situation) rates for this binding process.
Moreover, since this binding complex can also be unfolded, we have
C 3 + V h k C 2 V k C 2 V + C 3 - V h .
Aside from that, for the dissociation process of C 3 - V h , we assume it has the same binding rate as C 2 - V h , e.g., k C 2 V + = k C V + and k C 2 V = k C V . However, this process may also be force-related [87], and as the force keeps increasing (>30 pN), there will be a conformation change in α -catenin that causes the dissociation of Vh, but for now, we simply assume that this binding process is force-independent.
As explained in [74,88], the bindings of J u b and V h are independent of each other, and thus they can share the same binding parameter (e.g., the binding process of C 2 and C 2 - J u b with V h shares the same binding parameter).

Appendix A.3. The Detailed Model of the Tissue Mechanics

The basic structure is the hexagonal network shown in Figure 12A. (We only sketched the derivation of the equations and how the parameters were estimated. Further details can be obtained by contacting the authors. ) At each vertex, there is an element of mass m i , and the force balance for each vertex can be written as
m i x ¨ i ( t ) = i F i ( x , t ) ,
where x denotes a vector of all the points x i . There are connections between vertices via the edges and via the transverse connections in the network. The forces between vertices are represented by the Kelvin–Voigt module, with the active forces shown in Figure 12B. Each spring in a network of KVA elements has a resting length l , and the modules can have different parameters on the edges than on the transverse connections. We use the linearized force–velocity relation for the motor activity given in [147,148] and write it in the form
F a c t i v e = F s ( 1 v v 0 )
where F s is the stall force, the free velocity is v 0 , and v = | d l d t | is the speed.
The evolution equation for the point mass m i connected to the point masses m j via a KVA connection is given by
m i x ¨ i ( t ) = F i e l a s t i c ( x ) + F i d a m p i n g ( x ) + F i a c t i v e ( x ) = j N ( v i ) ( k i , j ( x i , j l i , j ) + ξ i , j d x i , j d t F i , j s ( 1 1 v i , j 0 | d x i , j d t | ) ) x j x i x j x i + F i , i m p
where N ( v i ) denotes the set of vertices adjacent x i , j is the length of boundary element e i , j connecting i and j, d x i , j d t is its time derivative, l i , j is the resting length related to the spring at boundary element e i , j , x j x i x j x i is the normalized vector accounting for the orientation and direction of the element e i , j , and F i a c t i v e is the active force described earlier.
To ease the notation, the above equations are written for the forces along the belt connections between nodes, but an identical set of equations, albeit with different parameters, can be written for the transverse connections that represent the connections in the apical network. Taken together, they give a complete description of the evolution equation for the network. The parameters used for computations based on the KVA model are given in the following table.
Table A3. Parameters for the mechanical model.
Table A3. Parameters for the mechanical model.
NameDescriptionValuesUnitsRef.
Parameter
ρ 0 Density of the cell1g/cm 3
A 0 Area of the cell 5.85 μ m 2
k 0 Elastic stiffness coefficient of a single stress fiber 4.14 × 10 5 N/m[149]
τ s Timescale constant ( ξ 0 / k 0 )2s[150]
v 0 Maximum value of polymerization velocity 6.92 × 10 7 m/s[151]
F 0 s Stall force of a single polymerizing actin filament1 (1–10)pN[152,153]

Appendix A.3.1. Interaction with the Biochemical Pathway

Force Applied on the α-Catenin and Expanded

Next, we outline how to compute the components of forces exerted on the α -catenin. In Figure A1 (left), we note that the forces we computed stemmed from the transverse connections across the cell, and here we ignored the forces in the actin belts (the lateral connections between nodes) as far as unfurling was concerned.
Figure A1. The transverse components of the forces on α -catenin in the undeformed and deformed states.
Figure A1. The transverse components of the forces on α -catenin in the undeformed and deformed states.
Cancers 15 04840 g0a1
If we consider the node labeled 5, then in the hexagonal configuration, the force is in the diagonal direction x 2 x 5 , and if we ignore damping and denote the cytoskeleton force between x 1 and x 5 , x 2 and x 5 , and x 3 and x 5 as F 1 , 5 , F 2 , 5 , and F 3 , 5 , respectively, then the force is given by
x 5 x 2 x 5 x 2 | F 1 , 5 | x 5 x 1 x 5 x 1 + | F 2 , 5 | x 5 x 2 x 5 x 2 + | F 3 , 5 | x 5 x 3 x 5 x 3 .
When the cell is deformed, as in Figure A1 (right), the calculation is similar, but the neighbors to 2 contribute different amounts to the force at 5 since the forces between 1 and 5 will generally be different than those between 3 and 5.

Yki-Controlled Cell Growth

Note that active Yki not only acts in the nucleus to stimulate growth but is also recruited to apical junctions where it enhances myosin II activity, which in turn promotes tissue growth through the Hippo/Yki pathway [82,116]. Because the details of how all the components are involved interact are not known, we used the following simple quadratic function to model the growth rate:
r ( [ Y k i ] ) = γ 1 ( [ Y k i ] [ Y k i ] 0 ) ( [ Y k i ] 1 [ Y k i ] ) [ Y k i ] < [ Y k i ] 1 , 0 otherwise .
This rate is negative when [Yki] = 0, which allows for cell death at insufficient Yki levels, and when [ Y k i ] > [ Y k i ] 1 to reflect the fact that other components may be limiting. We simply set γ 1 = 1 /h.

References

  1. Strassburger, K.; Lutz, M.; Müller, S.; Teleman, A.A. Ecdysone regulates Drosophila wing disc size via a TORC1 dependent mechanism. Nat. Commun. 2021, 12, 6684. [Google Scholar] [CrossRef]
  2. Ye, X.; Deng, Y.; Lai, Z.C. Akt is negatively regulated by Hippo signaling for growth inhibition in Drosophila. Dev. Biol. 2012, 369, 115–123. [Google Scholar] [CrossRef] [PubMed]
  3. Borreguero-Muñoz, N.; Fletcher, G.C.; Aguilar-Aragon, M.; Elbediwy, A.; Vincent-Mistiaen, Z.I.; Thompson, B.J. The Hippo pathway integrates PI3K–Akt signals with mechanical and polarity cues to control tissue growth. PLoS Biol. 2019, 17, e3000509. [Google Scholar] [CrossRef] [PubMed]
  4. Hua, H.; Kong, Q.; Zhang, H.; Wang, J.; Luo, T.; Jiang, Y. Targeting mTOR for cancer therapy. J. Hematol. Oncol. 2019, 12, 1–19. [Google Scholar] [CrossRef]
  5. Paton, E.L.; Turner, J.A.; Schlaepfer, I.R. Overcoming resistance to therapies targeting the MAPK pathway in BRAF-mutated tumours. J. Oncol. 2020, 2020. [Google Scholar] [CrossRef] [PubMed]
  6. Elosegui-Artola, A.; Andreu, I.; Beedle, A.E.; Lezamiz, A.; Uroz, M.; Kosmalska, A.J.; Oria, R.; Kechagia, J.Z.; Rico-Lastres, P.; Le Roux, A.L.; et al. Force triggers YAP nuclear entry by regulating transport across nuclear pores. Cell 2017, 171, 1397–1410. [Google Scholar] [CrossRef] [PubMed]
  7. Dasgupta, I.; McCollum, D. Control of cellular responses to mechanical cues through YAP/TAZ regulation. J. Biol. Chem. 2019, 294, 17693–17706. [Google Scholar] [CrossRef]
  8. Kaszak, I.; Witkowska-Piłaszewicz, O.; Niewiadomska, Z.; Dworecka-Kaszak, B.; Ngosa Toka, F.; Jurka, P. Role of cadherins in cancer—A review. Int. J. Mol. Sci. 2020, 21, 7624. [Google Scholar] [CrossRef]
  9. Gou, J.; Lin, L.; Othmer, H.G. A Model for the Hippo Pathway in the Drosophila Wing Disc. Biophys. J. 2018, 115, 737–747, PMID: 30041810. [Google Scholar] [CrossRef]
  10. Csibi, A.; Blenis, J. Hippo–YAP and mTOR pathways collaborate to regulate organ size. Nat. Cell Biol. 2012, 14, 1244–1245. [Google Scholar] [CrossRef]
  11. Mo, J.S.; Park, H.W.; Guan, K.L. The Hippo signaling pathway in stem cell biology and cancer. EMBO Rep. 2014, 15, 642–656. [Google Scholar] [CrossRef] [PubMed]
  12. Shraiman, B.I. Mechanicall feedback as a possible regulator of tissue growth. Proc. Nat. Acad. Sci. USA 2005, 102, 3318–3323. [Google Scholar] [CrossRef] [PubMed]
  13. Aegerter-Wilmsen, T.; Aegerter, C.M.; Bisseling, T. Model for the robust establishment of precise proportions in the early Drosophila embryo. J. Theor. Biol. 2005, 234, 13–19. [Google Scholar] [CrossRef] [PubMed]
  14. Irvine, K.D.; Shraiman, B.I. Mechanical control of growth: Ideas, facts and challenges. Development 2017, 144, 4238–4248. [Google Scholar] [CrossRef]
  15. Aegerter-Wilmsen, T.; Heimlicher, M.B.; Smith, A.C.; de Reuille, P.B.; Smith, R.S.; Aegerter, C.M.; Basler, K. Integrating force-sensing and signaling pathways in a model for the regulation of wing imaginal disc size. Development 2012, 139, 3221–3231. [Google Scholar] [CrossRef] [PubMed]
  16. Pan, D. Hippo signaling in organ size control. Genes Dev. 2007, 21, 886–897. [Google Scholar] [CrossRef]
  17. Ikmi, A.; Gaertner, B.; Seidel, C.; Srivastava, M.; Zeitlinger, J.; Gibson, M.C. Molecular evolution of the Yap/Yorkie proto-oncogene and elucidation of its core transcriptional program. Mol. Biol. Evol. 2014, 31, 1375–1390. [Google Scholar] [CrossRef]
  18. Willecke, M.; Hamaratoglu, F.; Kango-Singh, M.; Udan, R.; Chen, C.l.; Tao, C.; Zhang, X.; Halder, G. The fat cadherin acts through the hippo tumor-suppressor pathway to regulate tissue size. Curr. Biol. 2006, 16, 2090–2100. [Google Scholar] [CrossRef]
  19. Andl, T. miRNAs: Miracle or Mirage? The Limes Against the Barbaric Floods of Leaky and Undesired Transcripts. Organogenesis 2007, 3, 25–33. [Google Scholar] [CrossRef]
  20. Zhang, Y.; Wang, X.; Matakatsu, H.; Fehon, R.; Blair, S.S. The novel SH3 domain protein Dlish/CG10933 mediates fat signaling in Drosophila by binding and regulating Dachs. Elife 2016, 5, e16624. [Google Scholar] [CrossRef] [PubMed]
  21. Cho, E.; Feng, Y.; Rauskolb, C.; Maitra, S.; Fehon, R.; Irvine, K.D. Delineation of a Fat tumor suppressor pathway. Nat. Genet. 2006, 38, 1142. [Google Scholar] [CrossRef]
  22. Cho, E.; Irvine, K.D. Action of fat, four-jointed, dachsous and dachs in distal-to-proximal wing signaling. Development 2004, 131, 4489–4500. [Google Scholar] [CrossRef] [PubMed]
  23. Mao, Y.; Rauskolb, C.; Cho, E.; Hu, W.L.; Hayter, H.; Minihan, G.; Katz, F.N.; Irvine, K.D. Dachs: An unconventional myosin that functions downstream of Fat to regulate growth, affinity and gene expression in Drosophila. Development 2006, 133, 2539–2551. [Google Scholar] [CrossRef] [PubMed]
  24. Degoutin, J.L.; Milton, C.C.; Yu, E.; Tipping, M.; Bosveld, F.; Yang, L.; Bellaiche, Y.; Veraksa, A.; Harvey, K.F. Riquiqui and minibrain are regulators of the hippo pathway downstream of Dachsous. Nat. Cell Biol. 2013, 15, 1176. [Google Scholar] [CrossRef] [PubMed]
  25. Misra, J.R.; Irvine, K.D. Vamana couples fat signaling to the hippo pathway. Developmental Cell 2016, 39, 254–266. [Google Scholar] [CrossRef]
  26. Ishikawa, H.O.; Takeuchi, H.; Haltiwanger, R.S.; Irvine, K.D. Four-jointed is a Golgi kinase that phosphorylates a subset of cadherin domains. Science 2008, 321, 401–404. [Google Scholar] [CrossRef]
  27. Simon, M.A.; Xu, A.; Ishikawa, H.O.; Irvine, K.D. Modulation of fat: Dachsous binding by the cadherin domain kinase four-jointed. Curr. Biol. 2010, 20, 811–817. [Google Scholar] [CrossRef]
  28. Matakatsu, H.; Blair, S.S. Interactions between Fat and Dachsous and the regulation of planar cell polarity in the Drosophila wing. Development 2004, 131, 3785–3794. [Google Scholar] [CrossRef] [PubMed]
  29. Feng, Y.; Irvine, K.D. Processing and phosphorylation of the Fat receptor. Proc. Natl. Acad. Sci. USA 2009, 106, 11989–11994. [Google Scholar] [CrossRef]
  30. Rogulja, D.; Rauskolb, C.; Irvine, K.D. Morphogen control of wing growth through the Fat signaling pathway. Dev. Cell 2008, 15, 309–321. [Google Scholar] [CrossRef] [PubMed]
  31. Hamaratoglu, F.; Willecke, M.; Kango-Singh, M.; Nolo, R.; Hyun, E.; Tao, C.; Jafar-Nejad, H.; Halder, G. The tumour-suppressor genes NF2/Merlin and Expanded act through Hippo signalling to regulate cell proliferation and apoptosis. Nat. Cell Biol. 2006, 8, 27. [Google Scholar] [CrossRef] [PubMed]
  32. Yu, J.; Zheng, Y.; Dong, J.; Klusza, S.; Deng, W.M.; Pan, D. Kibra functions as a tumor suppressor protein that regulates Hippo signaling in conjunction with Merlin and Expanded. Dev. Cell 2010, 18, 288–299. [Google Scholar] [CrossRef] [PubMed]
  33. Oh, H.; Reddy, B.; Irvine, K.D. Phosphorylation-independent repression of Yorkie in Fat-Hippo signaling. Dev. Biol. 2009, 335, 188–197. [Google Scholar] [CrossRef]
  34. Ling, C.; Zheng, Y.; Yin, F.; Yu, J.; Huang, J.; Hong, Y.; Wu, S.; Pan, D. The apical transmembrane protein Crumbs functions as a tumor suppressor that regulates Hippo signaling by binding to Expanded. Proc. Natl. Acad. Sci. USA 2010, 107, 10532–10537. [Google Scholar] [CrossRef]
  35. Ribeiro, P.; Holder, M.; Frith, D.; Snijders, A.P.; Tapon, N. Crumbs promotes expanded recognition and degradation by the SCFSlimb/β-TrCP ubiquitin ligase. Proc. Natl. Acad. Sci. USA 2014, 111, E1980–E1989. [Google Scholar] [CrossRef]
  36. Sun, S.; Reddy, B.; Irvine, K.D. Localization of Hippo signalling complexes and Warts activation in vivo. Nat. Commun. 2015, 6, 8402. [Google Scholar] [CrossRef] [PubMed]
  37. Badouel, C.; Gardano, L.; Amin, N.; Garg, A.; Rosenfeld, R.; Le Bihan, T.; McNeill, H. The FERM-domain protein Expanded regulates Hippo pathway activity via direct interactions with the transcriptional activator Yorkie. Dev. Cell 2009, 16, 411–420. [Google Scholar] [CrossRef] [PubMed]
  38. Su, T.; Ludwig, M.Z.; Xu, J.; Fehon, R.G. Kibra and merlin activate the hippo pathway spatially distinct from and independent of expanded. Dev. Cell 2017, 40, 478–490. [Google Scholar] [CrossRef]
  39. Baumgartner, R.; Poernbacher, I.; Buser, N.; Hafen, E.; Stocker, H. The WW domain protein Kibra acts upstream of Hippo in Drosophila. Dev. Cell 2010, 18, 309–316. [Google Scholar] [CrossRef]
  40. Yin, F.; Yu, J.; Zheng, Y.; Chen, Q.; Zhang, N.; Pan, D. Spatial organization of Hippo signaling at the plasma membrane mediated by the tumor suppressor Merlin/NF2. Cell 2013, 154, 1342–1355. [Google Scholar] [CrossRef]
  41. Bennett, F.C.; Harvey, K.F. Fat cadherin modulates organ size in Drosophila via the Salvador/Warts/Hippo signaling pathway. Curr. Biol. 2006, 16, 2101–2110. [Google Scholar] [CrossRef] [PubMed]
  42. Silva, E.; Tsatskis, Y.; Gardano, L.; Tapon, N.; McNeill, H. The tumor-suppressor gene fat controls tissue growth upstream of expanded in the hippo signaling pathway. Curr. Biol. 2006, 16, 2081–2089. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, X.; Zhang, Y.; Blair, S.S. Fat-regulated adaptor protein Dlish binds the growth suppressor Expanded and controls its stability and ubiquitination. Proc. Natl. Acad. Sci. USA 2019, 116, 1319–1324. [Google Scholar] [CrossRef] [PubMed]
  44. Fulford, A.D.; Enderle, L.; Rusch, J.; Hodzic, D.; Holder, M.V.; Earl, A.; Oh, R.H.; Tapon, N.; McNeill, H. Expanded directly binds conserved regions of Fat to restrain growth via the Hippo pathway. J. Cell Biol. 2023, 222, e202204059. [Google Scholar] [CrossRef]
  45. Feng, Y.; Irvine, K.D. Fat and expanded act in parallel to regulate growth through warts. Proc. Natl. Acad. Sci. USA 2007, 104, 20362–20367. [Google Scholar] [CrossRef]
  46. Campbell, H.K.; Maiers, J.L.; DeMali, K.A. Interplay between tight junctions & adherens junctions. Exp. Cell Res. 2017, 358, 39–44. [Google Scholar]
  47. Thomas, C.; Strutt, D. The roles of the cadherins Fat and Dachsous in planar polarity specification in Drosophila. Dev. Dyn. 2012, 241, 27–39. [Google Scholar] [CrossRef]
  48. Leckband, D.E.; de Rooij, J. Cadherin Adhesion and Mechanotransduction. Annu. Rev. Cell Dev. Biol. 2014, 30, 291–315. [Google Scholar] [CrossRef]
  49. Pinheiro, D.; Bellaiche, Y. Mechanical force-driven adherens junction remodeling and epithelial dynamics. Dev. Cell 2018, 47, 3–19. [Google Scholar] [CrossRef]
  50. Tang, V.W.; Brieher, W.M. α-Actinin-4/FSGS1 is required for Arp2/3-dependent actin assembly at the adherens junction. J. Cell Biol. 2012, 196, 115–130. [Google Scholar] [CrossRef]
  51. Pinheiro, D.; Hannezo, E.; Herszterg, S.; Bosveld, F.; Gaugue, I.; Balakireva, M.; Wang, Z.; Cristo, I.; Rigaud, S.U.; Markova, O.; et al. Transmission of cytokinesis forces via E-cadherin dilution and actomyosin flows. Nature 2017, 545, 103. [Google Scholar] [CrossRef] [PubMed]
  52. Pan, Y.; Alégot, H.; Rauskolb, C.; Irvine, K.D. The dynamics of hippo signaling during Drosophila wing development. Development 2018, 145, dev165712. [Google Scholar] [CrossRef] [PubMed]
  53. Major, R.J.; Irvine, K.D. Localization and requirement for Myosin II at the dorsal-ventral compartment boundary of the Drosophila wing. Dev. Dyn. 2006, 235, 3051–3058. [Google Scholar] [CrossRef]
  54. Aliee, M.; Röper, J.C.; Landsberg, K.P.; Pentzold, C.; Widmann, T.J.; Jülicher, F.; Dahmann, C. Physical mechanisms shaping the Drosophila dorsoventral compartment boundary. Curr. Biol. 2012, 22, 967–976. [Google Scholar] [CrossRef] [PubMed]
  55. Sun, S.; Irvine, K.D. Cellular organization and cytoskeletal regulation of the Hippo signaling network. Trends Cell Biol. 2016, 26, 694–704. [Google Scholar] [CrossRef] [PubMed]
  56. Lecuit, T.; Yap, A.S. E-cadherin junctions as active mechanical integrators in tissue dynamics. Nat. Cell Biol. 2015, 17, 533. [Google Scholar] [CrossRef]
  57. Hong, S.; Troyanovsky, R.B.; Troyanovsky, S.M. Cadherin exits the junction by switching its adhesive bond. J. Cell Biol. 2011, 192, 1073–1083. [Google Scholar] [CrossRef]
  58. Angulo-Urarte, A.; van der Wal, T.; Huveneers, S. Cell-cell junctions as sensors and transducers of mechanical forces. Biochim. Biophys. Acta-Bba-Biomembr. 2020, 1862, 183316. [Google Scholar] [CrossRef] [PubMed]
  59. Bryant, D.M.; Stow, J.L. The ins and outs of E-cadherin trafficking. Trends in Cell Biology 2004, 14, 427–434. [Google Scholar] [CrossRef] [PubMed]
  60. Wirtz-Peitz, F.; Zallen, J.A. Junctional trafficking and epithelial morphogenesis. Curr. Opin. Genet. Dev. 2009, 19, 350–356. [Google Scholar] [CrossRef]
  61. Röper, J.C.; Mitrossilis, D.; Stirnemann, G.; Waharte, F.; Brito, I.; Fernandez-Sanchez, M.E.; Baaden, M.; Salamero, J.; Farge, E. The major β-catenin/E-cadherin junctional binding site is a primary molecular mechano-transductor of differentiation in vivo. Elife 2018, 7, e33381. [Google Scholar] [CrossRef] [PubMed]
  62. Efimova, N.; Svitkina, T.M. Branched actin networks push against each other at adherens junctions to maintain cell–cell adhesion. J. Cell Biol. 2018, 217, 1827–1845. [Google Scholar] [CrossRef]
  63. Röper, K. Integration of cell–cell adhesion and contractile actomyosin activity during morphogenesis. Curr. Top. Dev. Biol. 2015, 112, 103–127. [Google Scholar]
  64. Bosveld, F.; Wang, Z.; Bellaïche, Y. Tricellular junctions: A hot corner of epithelial biology. Curr. Opin. Cell Biol. 2018, 54, 80–88. [Google Scholar] [CrossRef] [PubMed]
  65. Bosveld, F.; Markova, O.; Guirao, B.; Martin, C.; Wang, Z.; Pierre, A.; Balakireva, M.; Gaugue, I.; Ainslie, A.; Christophorou, N.; et al. Epithelial tricellular junctions act as interphase cell shape sensors to orient mitosis. Nature 2016, 530, 495. [Google Scholar] [CrossRef]
  66. Zhou, Z.; Alégot, H.; Irvine, K.D. Oriented Cell Divisions Are Not Required for Drosophila Wing Shape. Curr. Biol. 2019, 29, 856–864. [Google Scholar] [CrossRef]
  67. Hirata, H.; Tatsumi, H.; Sokabe, M. Zyxin emerges as a key player in the mechanotransduction at cell adhesive structures. Commun. Integr. Biol. 2008, 1, 192–195. [Google Scholar] [CrossRef] [PubMed]
  68. Oldenburg, J.; Van Der Krogt, G.; Twiss, F.; Bongaarts, A.; Habani, Y.; Slotman, J.A.; Houtsmuller, A.; Huveneers, S.; De Rooij, J. VASP, zyxin and TES are tension-dependent members of Focal Adherens Junctions independent of the α-catenin-vinculin module. Sci. Rep. 2015, 5, 17225. [Google Scholar] [CrossRef]
  69. Sala, S.; Van Troys, M.; Medves, S.; Catillon, M.; Timmerman, E.; Staes, A.; Schaffner-Reckinger, E.; Gevaert, K.; Ampe, C. Expanding the Interactome of TES by Exploiting TES Modules with Different Subcellular Localizations. J. Proteome Res. 2017, 16, 2054–2071. [Google Scholar] [CrossRef]
  70. Kale, G.R.; Yang, X.; Philippe, J.M.; Mani, M.; Lenne, P.F.; Lecuit, T. Distinct contributions of tensile and shear stress on E-cadherin levels during morphogenesis. Nat. Commun. 2018, 9, 5021. [Google Scholar] [CrossRef] [PubMed]
  71. Baum, B.; Georgiou, M. Dynamics of adherens junctions in epithelial establishment, maintenance, and remodeling. J. Cell Biol. 2011, 192, 907–917. [Google Scholar] [CrossRef]
  72. Fletcher, G.C.; Borreguero-Muñoz, N.; Holder, M.; Aguilar-Aragon, M.; Thompson, B.J. Mechanical strain regulates the Hippo pathway in Drosophila. Development 2018, 145, dev–159467. [Google Scholar] [CrossRef]
  73. Huveneers, S.; de Rooij, J. Mechanosensitive systems at the cadherin–F-actin interface. J. Cell Sci. 2013, 126, 403–413. [Google Scholar] [CrossRef] [PubMed]
  74. Sarpal, R.; Yan, V.; Kazakova, L.; Sheppard, L.; Yu, J.C.; Fernandez-Gonzalez, R.; Tepass, U. Role of α-Catenin and its mechanosensing properties in regulating Hippo/YAP-dependent tissue growth. PLOS Genet. 2019, 15, e1008454. [Google Scholar] [CrossRef] [PubMed]
  75. Alégot, H.; Markosian, C.; Rauskolb, C.; Yang, J.; Kirichenko, E.; Wang, Y.C.; Irvine, K.D. Recruitment of Jub by α-catenin promotes Yki activity and Drosophila wing growth. J. Cell Sci. 2019, 132, jcs222018. [Google Scholar] [CrossRef]
  76. Rauskolb, C.; Sun, S.; Sun, G.; Pan, Y.; Irvine, K.D. Cytoskeletal Tension Inhibits Hippo Signaling through an Ajuba-Warts Complex. Cell 2014, 158, 143–156. [Google Scholar] [CrossRef] [PubMed]
  77. Duszyc, K.; Viasnoff, V.; Yap, A.S. Mechanosensing and mechanotransduction at cell–cell junctions. CSH Perspect. Biol. 2018, 10, a028761. [Google Scholar]
  78. Fletcher, G.C.; Elbediwy, A.; Khanal, I.; Ribeiro, P.S.; Tapon, N.; Thompson, B.J. The Spectrin cytoskeleton regulates the Hippo signalling pathway. EMBO J. 2015, 34, 940–954. [Google Scholar] [CrossRef]
  79. Nola, S.; Daigaku, R.; Smolarczyk, K.; Carstens, M.; Martin-Martin, B.; Longmore, G.; Bailly, M.; Braga, V.M. Ajuba is required for Rac activation and maintenance of E-cadherin adhesion. J. Cell Biol. 2011, 195, 855–871. [Google Scholar] [CrossRef]
  80. Marie, H.; Pratt, S.J.; Betson, M.; Epple, H.; Kittler, J.T.; Meek, L.; Moss, S.J.; Troyanovsky, S.; Attwell, D.; Longmore, G.D.; et al. The LIM Protein Ajuba Is Recruited to Cadherin-dependent Cell Junctions through an Association with α-Catenin*. J. Biol. Chem. 2003, 278, 1220–1228. [Google Scholar] [CrossRef]
  81. Rauskolb, C.; Cervantes, E.; Madere, F.; Irvine, K.D. Organization and function of tension-dependent complexes at adherens junctions. J. Cell Sci. 2019, 132, jcs224063. [Google Scholar] [CrossRef] [PubMed]
  82. Xu, J.; Vanderzalm, P.J.; Ludwig, M.; Su, T.; Tokamov, S.A.; Fehon, R.G. Yorkie Functions at the Cell Cortex to Promote Myosin Activation in a Non-transcriptional Manner. Dev. Cell 2018, 46, 271–284.e5. [Google Scholar] [CrossRef]
  83. Bax, N.A.; Wang, A.; Huang, D.L.; Pokutta, S.; Weis, W.I.; Dunn, A.R. Multi-level Force-dependent Allosteric Enhancement of αE-catenin Binding to F-actin by Vinculin. J. Mol. Biol. 2023, 435, 167969. [Google Scholar] [CrossRef]
  84. Bell, G.I. Models for the specific adhesion of cells to cells. Science 1978, 200, 618–627. [Google Scholar] [CrossRef]
  85. Jagannathan, B.; Marqusee, S. Protein folding and unfolding under force. Biopolymers 2013, 99, 860–869. [Google Scholar] [CrossRef]
  86. Strunz, T.; Oroszlan, K.; Schäfer, R.; Güntherodt, H.J. Dynamic force spectroscopy of single DNA molecules. Proc. Natl. Acad. Sci. USA 1999, 96, 11277–11282. [Google Scholar] [CrossRef]
  87. Yao, M.; Qiu, W.; Liu, R.; Efremov, A.K.; Cong, P.; Seddiki, R.; Payre, M.; Lim, C.T.; Ladoux, B.; Mège, R.M.; et al. Force-dependent conformational switch of α-catenin controls vinculin binding. Nat. Commun. 2014, 5, 4525. [Google Scholar] [CrossRef]
  88. Sheppard, L.; Green, D.G.; Lerchbaumer, G.; Rothenberg, K.E.; Fernandez-Gonzalez, R.; Tepass, U. The α-Catenin mechanosensing M region is required for cell adhesion during tissue morphogenesis. J. Cell Biol. 2022, 222, e202108091. [Google Scholar] [CrossRef] [PubMed]
  89. Rauskolb, C.; Irvine, K.D. Localization of Hippo signaling components in drosophila by fluorescence and immunofluorescence. In The Hippo Pathway; Springer: Berlin/Heidelberg, Germany, 2019; pp. 61–73. [Google Scholar]
  90. Razzell, W.; Bustillo, M.E.; Zallen, J.A. The force-sensitive protein Ajuba regulates cell adhesion during epithelial morphogenesis. J. Cell Biol. 2018, 217, 3715–3730. [Google Scholar] [CrossRef] [PubMed]
  91. Hu, J.; Matzavinos, A.; Othmer, H.G. A theoretical approach to actin filament dynamics. J. Stat. Phys. 2007, 128, 111–138. [Google Scholar] [CrossRef]
  92. Matzavinos, A.; Othmer, H.G. A stochastic analysis of actin polymerization in the presence of twinfilin and gelsolin. J. Theor. Biol. 2007, 249, 723–736. [Google Scholar] [CrossRef] [PubMed]
  93. Blanchoin, L.; Boujemaa-Paterski, R.; Sykes, C.; Plastino, J. Actin dynamics, architecture, and mechanics in cell motility. Physiol. Rev. 2014, 94, 235–263. [Google Scholar] [CrossRef]
  94. Iyer, K.V.; Piscitello-Gómez, R.; Paijmans, J.; Jülicher, F.; Eaton, S. Epithelial viscoelasticity is regulated by mechanosensitive E-cadherin turnover. Curr. Biol. 2019, 29, 578–591. [Google Scholar] [CrossRef]
  95. Villa, C.; Chaplain, M.A.; Gerisch, A.; Lorenzi, T. Mechanical models of pattern and form in biological tissues: The role of stress–strain constitutive equations. Bull. Math. Biol. 2021, 83, 80. [Google Scholar] [CrossRef]
  96. Kollmannsberger, P.; Fabry, B. Linear and nonlinear rheology of living cells. Annu. Rev. Mater. Res. 2011, 41, 75–97. [Google Scholar] [CrossRef]
  97. Salbreux, G.; Barthel, L.K.; Raymond, P.A.; Lubensky, D.K. Coupling Mechanical Deformations and Planar Cell Polarity to Create Regular Patterns in the Zebrafish Retina. PLoS Comput. Biol. 2012, 8, e1002618. [Google Scholar] [CrossRef]
  98. Charras, G.; Yap, A.S. Tensile forces and mechanotransduction at cell–cell junctions. Curr. Biol. 2018, 28, R445–R457. [Google Scholar] [CrossRef]
  99. Fischer, S.C.; Blanchard, G.B.; Duque, J.; Adams, R.J.; Arias, A.M.; Guest, S.D.; Gorfinkiel, N. Contractile and Mechanical Properties of Epithelia with Perturbed Actomyosin Dynamics. PloS ONE 2014, 9, e95695. [Google Scholar] [CrossRef] [PubMed]
  100. Farhadifar, R.; Röper, J.C.; Aigouy, B.; Eaton, S.; Jülicher, F. The influence of cell mechanics, cell-cell interactions, and proliferation on epithelial packing. Curr. Biol. 2007, 17, 2095–2104. [Google Scholar] [CrossRef]
  101. Hannezo, E.; Prost, J.; Joanny, J.F. Theory of epithelial sheet morphology in three dimensions. Proc. Nat. Acad. Sci. USA 2014, 111, 27–32. [Google Scholar] [CrossRef]
  102. Fletcher, A.G.; Osborne, J.M.; Maini, P.K.; Gavaghan, D.J. Implementing vertex dynamics models of cell populations in biology within a consistent computational framework. Prog. Biophys. Mol. Biol. 2013, 113, 299–326. [Google Scholar] [CrossRef] [PubMed]
  103. Rauzi, M.; Verant, P.; Lecuit, T.; Lenne, P.F. Nature and anisotropy of cortical forces orienting Drosophila tissue morphogenesis. Nat. Cell Biol. 2008, 10, 1401–1410. [Google Scholar] [CrossRef] [PubMed]
  104. Käfer, J.; Hayashi, T.; Marée, A.F.; Carthew, R.W.; Graner, F. Cell adhesion and cortex contractility determine cell patterning in the Drosophila retina. Proc. Nat. Acad. Sci. USA 2007, 104, 18549–18554. [Google Scholar] [CrossRef]
  105. Ma, X.; Lynch, H.E.; Scully, P.C.; Hutson, M.S. Probing embryonic tissue mechanics with laser hole drilling. Phys. Biol. 2009, 6, 036004. [Google Scholar] [CrossRef]
  106. Mao, Y.; Tournier, A.L.; Hoppe, A.; Kester, L.; Thompson, B.J.; Tapon, N. Differential proliferation rates generate patterns of mechanical tension that orient tissue growth. EMBO J. 2013, 32, 2790–2803. [Google Scholar] [CrossRef] [PubMed]
  107. Yonemura, S. Cadherin–actin interactions at adherens junctions. Curr. Opin. Cell Biol. 2011, 23, 515–522. [Google Scholar] [CrossRef] [PubMed]
  108. Borghi, N.; Sorokina, M.; Shcherbakova, O.G.; Weis, W.I.; Pruitt, B.L.; Nelson, W.J.; Dunn, A.R. E-cadherin is under constitutive actomyosin-generated tension that is increased at cell–cell contacts upon externally applied stretch. Proc. Natl. Acad. Sci. USA 2012, 109, 12568–12573. [Google Scholar] [CrossRef]
  109. Li, X.; Zhai, W.; Guo, Z.; Chen, H. Force-dependent unfolding and folding dynamics of protein alpha-catenin modulation domains. J. Innov. Opt. Health Sci. 2018, 12, 1841001. [Google Scholar] [CrossRef]
  110. Clarke, D.N.; Martin, A.C. Actin-based force generation and cell adhesion in tissue morphogenesis. Curr. Biol. 2021, 31, R667–R680. [Google Scholar] [CrossRef] [PubMed]
  111. Wang, A.; Dunn, A.R.; Weis, W.I.; Hill, C.P.; Pfeffer, S.R.; Yan, J. Mechanism of the cadherin–catenin F-actin catch bond interaction. eLife 2022, 11, e80130. [Google Scholar] [CrossRef]
  112. Silvestre, J.; Kenis, P.J.A.; Leckband, D.E. Cadherin and Integrin Regulation of Epithelial Cell Migration. Langmuir 2009, 25, 10092–10099. [Google Scholar] [CrossRef] [PubMed]
  113. Truong Quang, B.A.; Mani, M.; Markova, O.; Lecuit, T.; Lenne, P.F. Principles of E-Cadherin Supramolecular Organization In Vivo. Curr. Biol. 2013, 23, 2197–2207. [Google Scholar] [CrossRef]
  114. Yang, C.C.; Graves, H.K.; Moya, I.M.; Tao, C.; Hamaratoglu, F.; Gladden, A.B.; Halder, G. Differential regulation of the Hippo pathway by adherens junctions and apical–basal cell polarity modules. Proc. Natl. Acad. Sci. USA 2015, 112, 1785–1790. [Google Scholar] [CrossRef]
  115. Dupont, S.; Morsut, L.; Aragona, M.; Enzo, E.; Giulitti, S.; Cordenonsi, M.; Zanconato, F.; Le Digabel, J.; Forcato, M.; Bicciato, S.; et al. Role of YAP/TAZ in mechanotransduction. Nature 2011, 474, 179–183. [Google Scholar] [CrossRef] [PubMed]
  116. Koontz, L.M.; Liu-Chittenden, Y.; Yin, F.; Zheng, Y.; Yu, J.; Huang, B.; Chen, Q.; Wu, S.; Pan, D. The Hippo Effector Yorkie Controls Normal Tissue Growth by Antagonizing Scalloped-Mediated Default Repression. Dev. Cell 2013, 25, 388–401. [Google Scholar] [CrossRef] [PubMed]
  117. Irvine, K.D.; Harvey, K.F. Control of Organ Growth by Patterning and Hippo Signaling in Drosophila. Cold Spring Harb. Perspect. Biol. 2015, 7. [Google Scholar] [CrossRef] [PubMed]
  118. Yu, F.X.; Zhao, B.; Guan, K.L. Hippo Pathway in Organ Size Control, Tissue Homeostasis, and Cancer. Cell 2015, 163, 811–828. [Google Scholar] [CrossRef]
  119. Eder, D.; Aegerter, C.; Basler, K. Forces controlling organ growth and size. Mech. Dev. 2017, 144, 53–61. [Google Scholar] [CrossRef]
  120. Sknepnek, R.; Djafer-Cherif, I.; Chuai, M.; Weijer, C.; Henkes, S.; Kruse, K.; Walczak, A.M. Generating active T1 transitions through mechanochemical feedback. eLife 2023, 12, e79862. [Google Scholar] [CrossRef]
  121. Pan, Y.; Heemskerk, I.; Ibar, C.; Shraiman, B.I.; Irvine, K.D. Differential growth triggers mechanical feedback that elevates Hippo signaling. Proc. Natl. Acad. Sci. USA 2016, 113, E6974–E6983. [Google Scholar] [CrossRef]
  122. Herranz, H.; Weng, R.; Cohen, S.M. Crosstalk between epithelial and mesenchymal tissues in tumorigenesis and imaginal disc development. Curr. Biol. 2014, 24, 1476–1484. [Google Scholar] [CrossRef] [PubMed]
  123. Milán, M. Tumor Models: Tumor–Stroma Interactions Drive Neoplastic Transformation in Drosophila. Curr. Biol. 2014, 24, R658–R659. [Google Scholar] [CrossRef] [PubMed]
  124. Buchmann, A.; Alber, M.; Zartman, J.J. Sizing it up: The mechanical feedback hypothesis of organ growth regulation. Semin. Cell Dev. Biol. 2014, 35, 73–81. [Google Scholar] [CrossRef] [PubMed]
  125. Parker, J.; Struhl, G. Scaling the Drosophila wing: TOR-dependent target gene access by the Hippo pathway transducer Yorkie. PLoS Biol. 2015, 13, e1002274. [Google Scholar] [CrossRef]
  126. Gou, J.; Stotsky, J.A.; Othmer, H.G. Growth control in the Drosophila wing disk. Wiley Interdiscip. Rev. Syst. Biol. Med. 2020, 12, e1478. [Google Scholar] [CrossRef] [PubMed]
  127. Vollmer, J.; Casares, F.; Iber, D. Growth and size control during development. Open Biol. 2017, 7, 170190. [Google Scholar] [CrossRef] [PubMed]
  128. Ma, S.; Meng, Z.; Chen, R.; Guan, K.L. The Hippo pathway: Biology and pathophysiology. Annu. Rev. Biochem. 2019, 88, 577–604. [Google Scholar] [CrossRef]
  129. van Soldt, B.J.; Cardoso, W.V. Hippo-Yap/Taz signaling: Complex network interactions and impact in epithelial cell behavior. Wiley Interdiscip. Rev. Dev. Biol. 2020, 9, e371. [Google Scholar] [CrossRef]
  130. Fu, M.; Hu, Y.; Lan, T.; Guan, K.L.; Luo, T.; Luo, M. The Hippo signalling pathway and its implications in human health and diseases. Signal Transduct. Target. Ther. 2022, 7, 376. [Google Scholar] [CrossRef] [PubMed]
  131. Morishita, K.; Watanabe, K.; Ichijo, H. Cell volume regulation in cancer cell migration driven by osmotic water flow. Cancer Sci. 2019, 110, 2337–2347. [Google Scholar] [CrossRef]
  132. Catacuzzeno, L.; Sforna, L.; Esposito, V.; Limatola, C.; Franciolini, F. Ion channels in glioma malignancy. Transp. Malfunct. Cancer Spectr. 2020; 181, 223–267. [Google Scholar]
  133. Stock, C.; Schwab, A. Ion channels and transporters in metastasis. Biochim. Biophys. Acta-Bba-Biomembr. 2015, 1848, 2638–2646. [Google Scholar] [CrossRef]
  134. So, J.S.; Kim, H.; Han, K.S. Mechanisms of Invasion in Glioblastoma: Extracellular Matrix, Ca2+ Signaling, and Glutamate. Front. Cell. Neurosci. 2021, 15, 663092. [Google Scholar] [CrossRef]
  135. Stroka, K.M.; Jiang, H.; Chen, S.H.; Tong, Z.; Wirtz, D.; Sun, S.X.; Konstantopoulos, K. Water permeation drives tumor cell migration in confined microenvironments. Cell 2014, 157, 611–623. [Google Scholar] [CrossRef] [PubMed]
  136. Papadopoulos, M.C.; Saadoun, S. Key roles of aquaporins in tumor biology. Biochim. Biophys. Acta-BBA-Biomembr. 2015, 1848, 2576–2583. [Google Scholar] [CrossRef]
  137. Maki, K.; Han, S.W.; Hirano, Y.; Yonemura, S.; Hakoshima, T.; Adachi, T. Mechano-adaptive sensory mechanism of α-catenin under tension. Sci. Rep. 2016, 6, 24878. [Google Scholar] [CrossRef] [PubMed]
  138. Jukam, D.; Xie, B.; Rister, J.; Terrell, D.; Charlton-Perkins, M.; Pistillo, D.; Gebelein, B.; Desplan, C.; Cook, T. Opposite feedbacks in the Hippo pathway for growth control and neural fate. Science 2013, 342, 1238016. [Google Scholar] [CrossRef] [PubMed]
  139. Moroishi, T.; Park, H.W.; Qin, B.; Chen, Q.; Meng, Z.; Plouffe, S.W.; Taniguchi, K.; Yu, F.X.; Karin, M.; Pan, D.; et al. A YAP/TAZ-induced feedback mechanism regulates Hippo pathway homeostasis. Genes Dev. 2015, 29, 1271–1284. [Google Scholar] [CrossRef] [PubMed]
  140. Baines, A.J. Evolution of spectrin function in cytoskeletal and membrane networks. Biochem. Soc. Trans. 2009, 37, 796–803. [Google Scholar] [CrossRef]
  141. Deng, H.; Wang, W.; Yu, J.; Zheng, Y.; Qing, Y.; Pan, D. Spectrin regulates Hippo signaling by modulating cortical actomyosin activity. Elife 2015, 4, e06567. [Google Scholar] [CrossRef]
  142. England, J.L.; Cardy, J. Morphogen gradient from a noisy source. Phys. Rev. Lett. 2005, 94, 078101. [Google Scholar] [CrossRef]
  143. Shimmi, O.; O’Connor, M.B. Physical properties of Tld, Sog, Tsg and Dpp protein interactions are predicted to help create a sharp boundary in Bmp signals during dorsoventral patterning of the Drosophila embryo. Development 2003, 130, 4673–4682. [Google Scholar] [CrossRef] [PubMed]
  144. Kicheva, A.; Pantazis, P.; Bollenbach, T.; Kalaidzidis, Y.; Bittig, T.; Jülicher, F.; Gonzalez-Gaitan, M. Kinetics of morphogen gradient formation. Science 2007, 315, 521–525. [Google Scholar] [CrossRef]
  145. Pang, S.M.; Le, S.; Kwiatkowski, A.V.; Yan, J. Mechanical stability of αT-catenin and its activation by force for vinculin binding. Mol. Biol. Cell 2019, 30, 1930–1937. [Google Scholar] [CrossRef] [PubMed]
  146. Choi, H.J.; Pokutta, S.; Cadwell, G.W.; Bobkov, A.A.; Bankston, L.A.; Liddington, R.C.; Weis, W.I. αE-catenin is an autoinhibited molecule that coactivates vinculin. Proc. Natl. Acad. Sci. USA 2012, 109, 8576–8581. [Google Scholar] [CrossRef]
  147. Mayer, M.; Depken, M.; Bois, J.S.; Jülicher, F.; Grill, S.W. Anisotropies in cortical tension reveal the physical basis of polarizing cortical flows. Nature 2010, 467, 617–621. [Google Scholar] [CrossRef] [PubMed]
  148. Kassianidou, E.; Brand, C.A.; Schwarz, U.S.; Kumar, S. Geometry and network connectivity govern the mechanics of stress fibers. Proc. Natl. Acad. Sci. USA 2017, 114, 2622–2627. [Google Scholar] [CrossRef]
  149. De La Cruz, E.M.; Gardel, M.L. Actin Mechanics and Fragmentation*. J. Biol. Chem. 2015, 290, 17137–17144. [Google Scholar] [CrossRef]
  150. Hwang, Y.; Barakat, A.I. Dynamics of Mechanical Signal Transmission through Prestressed Stress Fibers. PLoS ONE 2012, 7, e35343. [Google Scholar] [CrossRef]
  151. Zimmermann, J.; Brunner, C.; Enculescu, M.; Goegler, M.; Ehrlicher, A.; Käs, J.; Falcke, M. Actin Filament Elasticity and Retrograde Flow Shape the Force-Velocity Relation of Motile Cells. Biophys. J. 2012, 102, 287–295. [Google Scholar] [CrossRef]
  152. Footer, M.J.; Kerssemakers, J.W.J.; Theriot, J.A.; Dogterom, M. Direct measurement of force generation by actin filament polymerization using an optical trap. Proc. Natl. Acad. Sci. USA 2007, 104, 2181–2186. [Google Scholar] [CrossRef]
  153. Jasnin, M.; Hervy, J.; Balor, S.; Bouissou, A.; Proag, A.; Voituriez, R.; Schneider, J.; Mangeat, T.; Maridonneau-Parini, I.; Baumeister, W.; et al. Elasticity of podosome actin networks produces nanonewton protrusive forces. Nat. Commun. 2022, 13, 3842. [Google Scholar] [CrossRef] [PubMed]
Figure 3. Basic components of the Hippo pathway described herein. Positive feedback steps occur on a slower time scale.
Figure 3. Basic components of the Hippo pathway described herein. Positive feedback steps occur on a slower time scale.
Cancers 15 04840 g003
Figure 4. (A) The membrane-localized Dachs and Riq as a function of Fat production. (B) The active Yki level as a function of Fat production with Fat-dependent stabilization of membrane Ex.
Figure 4. (A) The membrane-localized Dachs and Riq as a function of Fat production. (B) The active Yki level as a function of Fat production with Fat-dependent stabilization of membrane Ex.
Cancers 15 04840 g004
Figure 5. The mechanical structure of a wing disc cell. Green denotes the actin network.
Figure 5. The mechanical structure of a wing disc cell. Green denotes the actin network.
Cancers 15 04840 g005
Figure 6. Components of the junctional complex. PM = plasma membrane. Taken from [58] with permission.
Figure 6. Components of the junctional complex. PM = plasma membrane. Taken from [58] with permission.
Cancers 15 04840 g006
Figure 7. The domains of α -catenin. VBS = Vinculin binding site; ABD = actin binding domain. Taken from [58] with permission.
Figure 7. The domains of α -catenin. VBS = Vinculin binding site; ABD = actin binding domain. Taken from [58] with permission.
Cancers 15 04840 g007
Figure 8. A model of how binding of Jub to α -catenin controls tissue growth in the Drosophila wing disc. Taken from [74] with permission.
Figure 8. A model of how binding of Jub to α -catenin controls tissue growth in the Drosophila wing disc. Taken from [74] with permission.
Cancers 15 04840 g008
Figure 9. A simple model used to simulate the AJ under different forces.
Figure 9. A simple model used to simulate the AJ under different forces.
Cancers 15 04840 g009
Figure 10. (A) The progression of the states as the force is increased. (B) The variation between bound and unbound Jub as a function of force.
Figure 10. (A) The progression of the states as the force is increased. (B) The variation between bound and unbound Jub as a function of force.
Cancers 15 04840 g010
Figure 11. The force for the crossover from C 1 to C 2 as a function of Vinculin. The crossover from C 2 to C 3 is Vinculin-independent.
Figure 11. The force for the crossover from C 1 to C 2 as a function of Vinculin. The crossover from C 2 to C 3 is Vinculin-independent.
Cancers 15 04840 g011
Figure 12. The basic network (A), the KVA model (B), and the details of the notation (C), where v j denotes a node and m j its mass.
Figure 12. The basic network (A), the KVA model (B), and the details of the notation (C), where v j denotes a node and m j its mass.
Cancers 15 04840 g012
Figure 13. The dynamics under an initial increase in the size of the central cell. Black lines denote the initial configuration in each panel, red lines denote the deformed configuration. The table gives the conditions in the four panels, where a “+” force corresponds to stretching and a “-” force corresponds to compression. Growth is not incorporated here; size changes are due to rebalancing of the forces under the different conditions and result from elastic changes in the lengths of the edges.
Figure 13. The dynamics under an initial increase in the size of the central cell. Black lines denote the initial configuration in each panel, red lines denote the deformed configuration. The table gives the conditions in the four panels, where a “+” force corresponds to stretching and a “-” force corresponds to compression. Growth is not incorporated here; size changes are due to rebalancing of the forces under the different conditions and result from elastic changes in the lengths of the edges.
Cancers 15 04840 g013
Figure 14. (A) A qualitative representation of the Yki level in response to tension exerted on α -catenin. (B) A simple framework of the growth model.
Figure 14. (A) A qualitative representation of the Yki level in response to tension exerted on α -catenin. (B) A simple framework of the growth model.
Cancers 15 04840 g014
Figure 15. The heat map shows how the Yki level varies with different tension force applied on the CCC complex at each node and with the level of α -catenin at each node. The two plots on the right correspond to the two horizontal slices indicated by the white dashed lines shown in the heat map, where the level of α -catenin takes the value of 80 molecules/node (upper panel) and 20 molecules/node (lower panel).
Figure 15. The heat map shows how the Yki level varies with different tension force applied on the CCC complex at each node and with the level of α -catenin at each node. The two plots on the right correspond to the two horizontal slices indicated by the white dashed lines shown in the heat map, where the level of α -catenin takes the value of 80 molecules/node (upper panel) and 20 molecules/node (lower panel).
Cancers 15 04840 g015
Figure 16. The Yki level (nM/cell) in the boundary cell ([boundary]) and center cell ([center]) as a response to a change in the α -catenin level (molecules/node) in the center cell, while the α -catenin level in the boundary cells is constant at 30 molecules/node. Shown are the results for the boundary forces of −100 pN (compression) and 100 pN (tension). The dark- and light-blue curves coincide.
Figure 16. The Yki level (nM/cell) in the boundary cell ([boundary]) and center cell ([center]) as a response to a change in the α -catenin level (molecules/node) in the center cell, while the α -catenin level in the boundary cells is constant at 30 molecules/node. Shown are the results for the boundary forces of −100 pN (compression) and 100 pN (tension). The dark- and light-blue curves coincide.
Cancers 15 04840 g016
Figure 17. The difference in the Yki level (nM/cell) (left) and the ratio of change in the rest length (%) (right) between the boundary cells and the center cell as a response to a change in the α -catenin level (molecules/node) in the center cell, while the α -catenin level of the boundary cell is constant at 30 molecules/node. Here, the subscripts “c" and “b" indicate the center cell and a boundary cell, respectively, and R denotes the ratio of the final resting length to the initial resting length.
Figure 17. The difference in the Yki level (nM/cell) (left) and the ratio of change in the rest length (%) (right) between the boundary cells and the center cell as a response to a change in the α -catenin level (molecules/node) in the center cell, while the α -catenin level of the boundary cell is constant at 30 molecules/node. Here, the subscripts “c" and “b" indicate the center cell and a boundary cell, respectively, and R denotes the ratio of the final resting length to the initial resting length.
Cancers 15 04840 g017
Figure 18. Comparison of the Yki levels as a function of the α -catenin levels with and without growth of the center cell under different environmental forces. Other settings are the same as those shown in Figure 17.
Figure 18. Comparison of the Yki levels as a function of the α -catenin levels with and without growth of the center cell under different environmental forces. Other settings are the same as those shown in Figure 17.
Cancers 15 04840 g018
Figure 19. The difference in the Yki level (nM/cell) (left) and the ratio of change in the rest length (%) (right) between the boundary cells and the center cell as a response to a change in the production rate of Jub (nM/min) in the center cell, while the production rate of Jub for the boundary cell was set at a constant value of 250 nM/min.
Figure 19. The difference in the Yki level (nM/cell) (left) and the ratio of change in the rest length (%) (right) between the boundary cells and the center cell as a response to a change in the production rate of Jub (nM/min) in the center cell, while the production rate of Jub for the boundary cell was set at a constant value of 250 nM/min.
Cancers 15 04840 g019
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gou, J.; Zhang, T.; Othmer, H.G. The Interaction of Mechanics and the Hippo Pathway in Drosophila melanogaster. Cancers 2023, 15, 4840. https://doi.org/10.3390/cancers15194840

AMA Style

Gou J, Zhang T, Othmer HG. The Interaction of Mechanics and the Hippo Pathway in Drosophila melanogaster. Cancers. 2023; 15(19):4840. https://doi.org/10.3390/cancers15194840

Chicago/Turabian Style

Gou, Jia, Tianhao Zhang, and Hans G. Othmer. 2023. "The Interaction of Mechanics and the Hippo Pathway in Drosophila melanogaster" Cancers 15, no. 19: 4840. https://doi.org/10.3390/cancers15194840

APA Style

Gou, J., Zhang, T., & Othmer, H. G. (2023). The Interaction of Mechanics and the Hippo Pathway in Drosophila melanogaster. Cancers, 15(19), 4840. https://doi.org/10.3390/cancers15194840

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop