Pith. sign in

REVIEW 2 major objections 4 minor 59 references

Mathematical Modelling of Mechanotransduction via RhoA Signalling Pathways

T0 review · 2 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Pith's one-line read A two-way RhoA–FAK feedback loop can explain why cells deform similarly on soft and stiff surfaces.

desk verdict Solid two-way coupling model, but the homeostasis plateau is computed in a large-deformation regime the linear elasticity assumption cannot support. read the letter →

arxiv 2501.04407 v2 pith:PS7RVG26 submitted 2025-01-08 math.NA cs.NAphysics.bio-phq-bio.CB

classification math.NAcs.NAphysics.bio-phq-bio.CB MSC 65M6092C1035K5774B05
keywords mechanotransductionRhoAsignallingfocaladhesionkinasemechanicalhomeostasisreaction-diffusionequationsbulk-surfacefiniteelementslinearelasticitysubstratestiffness
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper tries to establish that mechanical homeostasis—cells deforming by roughly the same amount whether they sit on a soft or a stiff surface—can emerge from a minimal mechanistic model rather than from a tuned control system. The model couples the RhoA signalling pathway to cell elasticity in two directions: mechanical stress activates focal adhesion kinase (FAK), and active FAK increases the cell's Young's modulus. In simulations, this two-way loop makes the magnitude of cell deformation and cell volume change nearly flat across substrate stiffnesses from 0.1 kPa to 7 GPa, while a model with constant cell stiffness shows no such plateau. A sympathetic reader would take the paper to be claiming that the feedback loop, not any single reaction, is the mechanism behind stiffness homeostasis.

What carries the argument

The central object is a two-way feedback loop. Stress activates FAK through the term $C_1\,\mathrm{tr}(\sigma)_+\,\phi_d$, and the cell's Young's modulus increases with active FAK according to the power law $E_c(\phi_a)=k_7(1+(k_8\phi_a)^p)$, where active FAK is used as a proxy for F-actin. This loop is embedded in a bulk-surface reaction-diffusion system for inactive and active FAK in the cytoplasm, active RhoA on the membrane, and a linear elasticity problem whose boundary traction depends on active RhoA. The numerical machinery is a bulk-surface finite element method with semi-implicit time stepping; the power law is the bridge that converts a biochemical signal into a mechanical property, and removing that bridge removes the homeostasis.

What would settle it

Measure, in a single cell type, both active FAK concentration and deformation on substrates spanning 0.1 kPa to 7 GPa. If deformation follows substrate stiffness monotonically instead of plateauing, or if cell stiffness is unchanged when active FAK is elevated, the claimed two-way mechanism is falsified. Pharmacological FAK inhibition should also destroy the plateau if the mechanism is causal.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is that bidirectional coupling between signalling and cell mechanics is sufficient to reproduce mechanical homeostasis. The authors extend a spatial RhoA signalling model to include linear elastic cell deformation, with the cell's Young's modulus given by $E_c(\phi_a)=k_7(1+(k_8\phi_a)^p)$ and with FAK activation driven by the positive part of the stress trace, $C_1\,\mathrm{tr}(\sigma)_+\,\phi_d$. When both couplings are present, steady-state cell deformation and $\mathrm{div}(u)$ become nearly constant over a wide range of ECM stiffnesses, while active FAK and active RhoA still rise in a threshold-like Hill-function manner. This reproduces the qualitative experimental observations of mechanical homeostasis and of threshold-like stiffness responses. The plateau does not occur when the cell is assigned constant stiffness or when stress feedback onto signalling is switched off, so the paper concludes that the bidirectional coupling itself is the load-bearing mechanism.

Load-bearing premise

The homeostasis prediction rests on the assumed law that cell stiffness rises with active FAK according to $E_c(\phi_a)=k_7(1+(k_8\phi_a)^p)$, with active FAK standing in for F-actin; if real cells stiffen differently, or if FAK is not a good proxy for F-actin, the plateau could disappear.

Editorial extensions

If this is right

  • If the central claim is correct, models that treat cell stiffness as a fixed constant will miss the homeostasis plateau; stiffness must be allowed to respond to signalling.
  • The reduced FAK–RhoA model reproduces the threshold-like dependence of signalling on substrate stiffness seen in experiments, so the threshold response is not an artefact of the full pathway.
  • Cell shape becomes a control variable: axisymmetric and lamellipodium-shaped cells show different concentrations and different threshold responses, implying that geometry is part of the mechanosensing circuit.
  • A linear viscoelastic version of the model produces the same qualitative behaviour, so the homeostasis result does not depend on the purely elastic constitutive law.
  • The framework is designed to be extended to viscoelastic or poroelastic mechanics, nuclear deformation, and spatially varying substrate stiffness.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper, the same feedback loop would predict that a cell pre-exposed to a stiff substrate retains a higher Young's modulus and therefore deforms less on a later soft substrate, offering a concrete mechanism for mechanical memory that could be tested by sequential-culture experiments.
  • Beyond the paper, pharmacological inhibition of FAK or RhoA should erase the plateau: deformation should then scale with substrate stiffness as it does in the constant-stiffness case, providing a direct causal test of the model.
  • Beyond the paper, the fitted parameters $k_8$ and $p$ could be estimated from combined live imaging of active FAK and atomic force microscopy stiffness maps, turning the fitted constitutive law into a measurable prediction.
  • Beyond the paper, because the plateau appears only when both couplings are active, one-way stiffness-sensing models may misestimate how cells respond to matrix gradients, especially in slowly varying environments where the bidirectional term dominates.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The manuscript proposes a reduced bulk-surface reaction-diffusion model for inactive and active FAK and active RhoA, coupled to linear elasticity of the cell. Substrate stiffness activates FAK on the membrane, the cell Young's modulus is made to depend on active FAK concentration through Eq. (2), and the boundary traction depends on active RhoA through Eq. (4). The model is simulated for 2D and 3D stimuli, for axisymmetric and lamellipodium-shaped cells, and the results are used to claim a threshold-like response to substrate stiffness and a robustness of cell deformation to stiffness changes, interpreted as mechanical homeostasis. The paper also presents a bulk-surface finite element method and verifies its convergence on a manufactured solution.

Significance. If the central claim holds, the paper provides a minimal mechanistic explanation of mechanical homeostasis: bidirectional coupling between RhoA signalling and cell mechanics can make cell deformation insensitive to ECM stiffness. The paper has clear strengths: the reduced model is shown to reproduce the qualitative results of the full model of [18]; the numerical scheme is benchmarked on a manufactured solution with the expected orders of convergence; and a parameter sensitivity analysis is provided. However, the main emergent predictions are tied to the assumed power-law stiffness law and to the linear small-strain approximation, so the significance of the homeostasis claim depends on additional validation that the manuscript does not currently provide.

major comments (2)
  1. [Mathematical model for mechanotransduction, Eq. (3); Results, Figs. 1 and 4] The model is formulated under the explicit assumption of small deformations, leading to the linear elastic equation (3). Yet the results report maximum displacement magnitudes of 7 µm for the 2xD stimulus (Fig. 1) and 7.5 µm for the 3D stimulus (Fig. 4) on a cell with |Y|=1193 µm^3, i.e. an equivalent spherical diameter of roughly 13 µm and radius of roughly 6.5 µm. Displacements comparable to the cell radius imply O(1) strains, far outside the small-strain regime in which Eq. (3) is valid. The central homeostasis claim is read off from div(u) and |u| in exactly this regime (Figs. 3, 5, and 7). The Discussion admits that the linear elastic assumption 'is limiting as it assumes small deformations', and the viscoelastic check in S1 Appendix A.7 retains linear kinematics. A finite-strain computation, or a restriction of the homeostasis claim to genuinely small strains, is needed before the plateau can be regarded as a robust prediction of the model rather than a possible artifact of the linearization.
  2. [Eq. (2), Table 1(b), and S1 Appendix A.6] The homeostatic plateau is largely a consequence of the assumed constitutive law E_c(phi_a)=k7(1+(k8 phi_a)^p), with k7, k8, and p fitted to results for phi_a in [18] and F-actin in [38]. The threshold-like dependence on substrate stiffness is similarly inherited from the Hill-type terms and parameters taken from [18]. Because these fitted inputs are used to generate the very curves that are then presented as emergent predictions, the claim that homeostasis 'emerges' from the bidirectional coupling is only as strong as the empirical support for Eq. (2). The sensitivity analysis in S1 Appendix A.6 varies the parameters of Eq. (2) by ±10% and ±20%, but it does not test the form of the law itself. The authors should either test alternative functional forms of E_c(phi_a), or present independent data that constrain Eq. (2), before the homeostasis result is presented as a model prediction rather than a restatement of the fitted stiffness law.
minor comments (4)
  1. [Methods, first paragraph] There is a typo: 'mechanostranduction' should be 'mechanotransduction'.
  2. [Results, Fig. 8 caption] The caption contains a duplicated phrase: 'simulation results showing rho_a, phi_d, phi_a and |u| simulation results showing phi_a and rho_a for the model...' should be cleaned up.
  3. [Introduction, last paragraph] 'mechanisms that underlay the mechanical homeostasis' should be 'mechanisms that underlie mechanical homeostasis'.
  4. [References] Reference [43] is cited before [42] in the text; the citation order in the bibliography could be adjusted if the journal requires strict numerical ordering.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: fitted inputs are external and the homeostasis plateau is not imposed by construction.

full rationale

The derivation chain is self-contained relative to the claims that matter. The biochemical submodel is explicitly a reduction of the Scott et al. model (Eq. 1), with the Michaelis-Menten E/(C+E) and Hill ((gamma*phi_a)^n+1) terms inherited from that external source; the paper frames the resulting threshold-like dependence as reproduction and validation against [18,44], not as a new prediction derived from mechanics. The abstract's phrase "novel emergent features such as ... threshold-like response" overstates novelty because that threshold is structurally present in Eq. (1), but this is a presentation issue, not a circular derivation of the paper's central homeostasis claim. The parameters k7, k8, and p in the cell-stiffness law Eq. (2) are fitted to external F-actin/FAK data ([18,38]), while the homeostasis plateau is not fitted to homeostasis data: it appears only when the bidirectional feedback sigma -> phi_a -> E_c is active, and it is compared qualitatively with [28]. Self-citations ([14,60]) occur only for standard numerical and multiscale tools and are not load-bearing. The paper explicitly flags the linear-elastic small-deformation limitation in the Discussion; that is an accuracy risk, not circularity. No equation is shown to be equivalent by construction to another, and no fitted parameter is renamed as a prediction of the central emergent quantity. The stated but deferred existence proof is a completeness gap, not a circular step.

Assumptions & free parameters 6 free parameters · 6 assumptions · 0 invented entities

The central claim depends on several parameters fitted to prior simulation data and on phenomenological coupling laws. No new physical entities, such as particles, forces, or dimensions, are introduced beyond the mathematical projection P, which is a construction rather than a physical entity.

free parameters (6)
  • C1 = 0.1 (kPa s)^-1; range 0-2 explored
    Coupling strength of cell stress to FAK activation; introduced in this paper, not fitted to specific data, but explored over a range.
  • k6 = 0.1 s^-1
    Traction coefficient in boundary condition (4); fitted to yield a deformation magnitude of 0-10 um as in [18].
  • k7 = 0.2 kPa
    Baseline Young's modulus in Eq (2); fitted to results for active FAK in [18].
  • k8 = 2.4245 dm^3/umol
    FAK sensitivity of cell stiffness in Eq (2); fitted to results for active FAK in [18].
  • p = 2.6
    Exponent in the stiffness law (2); fitted to results for active FAK and F-actin in [18] and [38].
  • nu_c = 0.3
    Poisson ratio of the cell; estimated from the literature range 0.17-0.66 [42].
assumptions (6)
  • domain assumption Linear elasticity with small deformations (Eq. 3)
    The cell is treated as a linear elastic material; small deformation justifies solving signalling on the reference configuration. Stated in the Methods section.
  • ad hoc to paper Cell Young's modulus is a power-law function of active FAK (Eq. 2)
    Proposed based on [18] and [38]; k7, k8, and p are fitted to simulation results. This law is load-bearing for the homeostasis claim.
  • domain assumption Substrate stiffness activates FAK through a Michaelis-Menten term k3 E/(C+E) (Eq. 1)
    Inherited from the reference model [18]; this functional form generates the threshold-like response.
  • domain assumption RhoA activation uses a Hill-type term k5((gamma phi_a)^n + 1)(M_rho/|Y| - rho_a/n_r) (Eq. 1)
    Inherited from [18]; assumes conservation of total RhoA mass.
  • ad hoc to paper FAK is activated by the positive part of the trace of the stress tensor, C1 tr(sigma)+ (Eq. 6)
    Phenomenological proxy for stress-induced FAK activation; no direct molecular derivation is given.
  • ad hoc to paper Boundary traction depends on active RhoA as k6 P(rho_a nu) (Eq. 4)
    Models RhoA-driven contractility; P projects out rigid motions.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Mathematical Modelling of Mechanotransduction via RhoA Signalling Pathways." pith.science (2026). https://pith.science/paper/PS7RVG26

@misc{pith2026250104407,
  author       = {Pith},
  title        = {Pith review of: Mathematical Modelling of Mechanotransduction via RhoA Signalling Pathways},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PS7RVG26}},
  note         = {Machine review of arXiv:2501.04407}
}
read the original abstract

We derive and simulate a mathematical model for mechanotransduction related to the Rho GTPase signalling pathway. The model addresses the bidirectional coupling between signalling processes and cell mechanics. A numerical method based on bulk-surface finite elements is proposed for the approximation of the coupled system of nonlinear reaction-diffusion equations, defined inside the cell and on the cell membrane, and the equations of elasticity. Our simulation results illustrate novel emergent features such as the strong dependence of the dynamics on cell shape, a threshold-like response to changes in substrate stiffness, and the fact that coupling mechanics and signalling can lead to the robustness of cell deformation to larger changes in substrate stiffness, ensuring mechanical homeostasis in agreement with experiments.

Figures

Figures reproduced from arXiv: 2501.04407 by the authors.

Figure 1
Figure 1. Numerical simulation results showing ρa, ϕd, ϕa and |u| for the model in Eqs (3)-(6) for the axisymmetric shape and in the case of 2xD stimulus at a steady state at T = 100 s. Four different scenarios are considered: (A) C1 = 0 (kPa s)−1 (σ ̸→ ϕa) and Ec = 0.6 kPa (ϕa ̸→ Ec); (B) C1 = 0 (kPa s)−1 (σ ̸→ ϕa) and Ec = f(ϕa) (ϕa → Ec); (C) C1 = 0.1 (kPa s)−1 (σ → ϕa) and Ec = 0.6 kPa (ϕa ̸→ Ec); (D) C1 = 0.1 (kPa s)−1 (… view at source ↗
Figure 2
Figure 2. Numerical simulation results showing ρa, ϕd, ϕa and |u| for the model in Eqs (3)-(6) for the lamellipodium shape and in the case of 2xD stimulus at a steady state at T = 100 s. Four different scenarios are considered: (A) C1 = 0 (kPa s)−1 (σ ̸→ ϕa) and Ec = 0.6 kPa (ϕa ̸→ Ec); (B) C1 = 0 (kPa s)−1 (σ ̸→ ϕa) and Ec = f(ϕa) (ϕa → Ec); (C) C1 = 0.1 (kPa s)−1 (σ → ϕa) and Ec = 0.6 kPa (ϕa ̸→ Ec); (D) C1 = 0.1 (kPa s)−1 … view at source ↗
Figure 3
Figure 3. Simulation results showing the mean, 1 |Ω| R Ω · dx, min and max values of f(ϕa), div(u), ϕa and ρa as functions of substrate stiffness E. We consider different couplings with four different values for C1 and two different shapes at T = 100 s by which time the results are at a steady state. All other parameter values as in [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (9 more)
Figure 4
Figure 4. Figure 4: Numerical simulation results showing ρa, ϕd, ϕa and |u| for the model in Eqs (3)-(6) for the axisymmetric shape and in the case of the 3D stimulus at a steady state at T = 100 s. Four different scenarios are considered: (A) C1 = 0 (kPa s)−1 (σ ̸→ ϕa) and Ec = 0.6 kPa (…
Figure 5
Figure 5. Figure 5: Simulation results showing the mean, 1 |Ω| R Ω · dx, min and max values of f(ϕa), div(u), ϕa and ρa as functions of substrate stiffness E. We consider different couplings with four different values for C1 and two different shapes at T = 100 s by which time the results …
Figure 6
Figure 6. Figure 6: Numerical simulation results showing ρa, ϕd, ϕa and |u| for the model in Eqs (3), (4), and (6) for the axisymmetric shape and in the case of the 3D stimulus at a steady state at T = 100 s. Four different scenarios are considered: (A) C1 = 0 (kPa s)−1 (σ ̸→ ϕa) and Ec =…
Figure 7
Figure 7. Figure 7: Simulation results showing the mean, 1 |Ω| R Ω · dx, min and max values of f(ϕa), div(u), ϕa and ρa as functions of substrate stiffness E, in the case of the model in Eqs (3), (4) and (6) and 3D stimulus. We consider different couplings, four different values for C1, a…
Figure 8
Figure 8. Figure 8: Numerical simulation results showing ρa, ϕd, ϕa and |u| simulation results showing ϕa and ρa for the model in Eqs (3), (4), and (6) for the axisymmetric shape and in the case of 2xD stimulus at a steady state at T = 100 s. Four different scenarios are considered: (A) C…
Figure 9
Figure 9. Figure 9: Numerical simulation results showing ρa, ϕd, ϕa and |u| for the model in Eqs (3)-(6) for the lamellipodium shape and in the case of the 3D stimulus at a steady state at T = 100 s. Four different scenarios are considered: (A) C1 = 0 (kPa s)−1 (σ ̸→ ϕa) and Ec = 0.6 kPa …
Figure 10
Figure 10. Figure 10: Numerical simulation results showing ρa, ϕd, ϕa and |u| for the model in Eqs (3), (4), and (6) for the lamellipodium shape and in the case of the 3D stimulus at a steady state at T = 100 s. Four different scenarios are considered: (A) C1 = 0 (kPa s)−1 (σ ̸→ ϕa) and Ec…
Figure 11
Figure 11. Figure 11: Numerical simulation results showing ρa, ϕd, ϕa and |u| for the model in Eqs (3), (4), and (6) for the lamellipodium shape cells and 2xD stimulus at a steady state at T = 100 s. Four different scenarios are considered: (A) C1 = 0 (kPa s)−1 (σ ̸→ ϕa) and Ec = 0.6 kPa (…
Figure 12
Figure 12. Figure 12: Simulation results showing the mean, 1 |Ω| R Ω · dx, min and max values of f(ϕa), div(u), ϕa and ρa as functions of substrate stiffness E, in the case of the model in Eqs (3), (4) and (6) and 2xD stimulus. We consider different couplings, four different values for C1,…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

59 extracted references · 50 canonical work pages

  1. [18]

    A spatial model of YAP/TAZ signaling reveals how stiffness, dimen- sionality, and shape contribute to emergent outcomes

    Scott KE, Fraley SI, Rangamani P. A spatial model of YAP/TAZ signaling reveals how stiffness, dimen- sionality, and shape contribute to emergent outcomes. Proceedings of the National Academy of Sciences. 2021;118(20):e2021571118. doi:10.1073/pnas.2021571118

  2. [38]

    Elastic Behavior of Cross- Linked and Bundled Actin Networks

    Gardel ML, Shin JH, MacKintosh FC, Mahadevan L, Matsudaira P, Weitz DA. Elastic Behavior of Cross- Linked and Bundled Actin Networks. Science. 2004;304(5675):1301–1305. doi:10.1126/science.1095087

  3. [1]

    A F AK-p120RasGAP-p190RhoGAP complex regulates polarity in migrating cells

    Tomar A, Lim ST, Lim Y, Schlaepfer DD. A F AK-p120RasGAP-p190RhoGAP complex regulates polarity in migrating cells. Journal of Cell Science. 2009;122:1852–1862. doi:10.1242/jcs.046870

  4. [2]

    Signalling dynamics, cell decisions, and homeostatic control in health and disease

    Valls PO, Esposito A. Signalling dynamics, cell decisions, and homeostatic control in health and disease. J Cell Sci Curr Opin Cell Biol. 2022;75:102066. doi:10.1016/j.ceb.2022.01.011

  5. [3]

    Crosstalk between mechanotransduction and metabolism

    Romani P, Valcarcel-Jimenez L, Frezza C, Dupont S. Crosstalk between mechanotransduction and metabolism. Nature Reviews Molecular Cell Biology. 2021;22(1):22–38. doi:10.1038/s41580-020-00306-w

  6. [4]

    Mechanoregulation of YAP and TAZ in Cellular Homeostasis and Disease Progres- sion

    Cai X, Wang KC, Meng Z. Mechanoregulation of YAP and TAZ in Cellular Homeostasis and Disease Progres- sion. Front Cell Dev Biol. 2021;9:673599. doi:10.3389/fcell.2021.673599

  7. [5]

    Mechanotransduction and extracellular matrix homeostasis

    Humphrey JD, Dufresne ER, Schwartz MA. Mechanotransduction and extracellular matrix homeostasis. Nature reviews Molecular cell biology. 2014;15(12):802–812. doi:10.1038/nrm3896

  8. [6]

    Cell–extracellular matrix mechanotransduction in 3D

    Saraswathibhatla A, Indana D, Chaudhuri O. Cell–extracellular matrix mechanotransduction in 3D. Nature Reviews Molecular Cell Biology. 2023;24(7):495–516. doi:10.1038/s41580-023-00583-1

Show all 59 references
  1. [7]

    Mechanotransduction: from the cell surface to the nucleus via RhoA

    Burridge K, Monaghan-Benson E, Graham DM. Mechanotransduction: from the cell surface to the nucleus via RhoA. Philosophical Transactions of the Royal Society B: Biological Sciences. 2019;374(1779):20180229. doi:10.1098/rstb.2018.0229

  2. [8]

    Cell response to mechanical microenvironment cues via Rho signaling: From mechanobiology to mechanomedicine

    Xie N, Xiao C, Shu Q, Cheng B, Wang Z, Xue R, et al. Cell response to mechanical microenvironment cues via Rho signaling: From mechanobiology to mechanomedicine. Acta Biomaterialia. 2023;159:1–20. doi:10.1016/j.actbio.2023.01.039

  3. [9]

    Cellular Mechanotransduction: From Tension to Function

    Martino F, Perestrelo AR, Vinarsk´ y V, Pagliari S, Forte G. Cellular Mechanotransduction: From Tension to Function. Frontiers in Physiology. 2018;9

  4. [10]

    Integrin-mediated mechanotransduction

    Sun Z, Guo SS, F¨ assler R. Integrin-mediated mechanotransduction. J Cell Biol. 2016;215(4):445–456. doi:10.1083/jcb.201609037

  5. [11]

    Environmental stiffness restores mechanical homeostasis in vimentin-depleted cells

    Young KM, Reinhart-King CA. Environmental stiffness restores mechanical homeostasis in vimentin-depleted cells. Current Opinion in Cell Biology. 2023;83:102208. doi:10.1016/j.ceb.2023.102208

  6. [12]

    Computational methodologies for modelling, analysis and simulation of signalling networks

    Gilbert D, Fuss H, Gu X, Orton R, Robinson S, Vyshemirsky V, et al. Computational methodologies for modelling, analysis and simulation of signalling networks. Briefings in Bioinformatics. 2006;7(4):339–353. doi:10.1093/bib/bbl043

  7. [13]

    Mathematical modelling and computational study of two- dimensional and three-dimensional dynamics of receptor–ligand interactions in signalling response mechanisms

    Garc ´ ıa-Penarrubia P, J J G´ alvez JJ, G´ alvez J. Mathematical modelling and computational study of two- dimensional and three-dimensional dynamics of receptor–ligand interactions in signalling response mechanisms. J Math Biol. 2014;69:553–582

  8. [14]

    Multiscale Analysis and Simulation of a Signaling Process With Surface Diffusion

    Ptashnyk M, Venkataraman C. Multiscale Analysis and Simulation of a Signaling Process With Surface Diffusion. Multiscale Modeling & Simulation. 2020;18(2):851–886. doi:10.1137/18M1185661

  9. [15]

    Cellular mechanosensing of the biophysical microenvi- ronment: A review of mathematical models of biophysical regulation of cell responses

    Cheng B, Lin M, Huang G, Li Y, Ji B, Genin GM, et al. Cellular mechanosensing of the biophysical microenvi- ronment: A review of mathematical models of biophysical regulation of cell responses. Physics of Life Reviews. 2017;22-23:88–119. doi:10.1016/j.plrev.2017.06.016. 15

  10. [16]

    Coupling biochemistry and mechanics in cell adhesion: a model for inhomogeneous stress fiber contraction

    Besser A, Schwarz US. Coupling biochemistry and mechanics in cell adhesion: a model for inhomogeneous stress fiber contraction. New Journal of Physics. 2007;9(11):425–425. doi:10.1088/1367-2630/9/11/425

  11. [17]

    Spatiotemporal model of cellular mechanotrans- duction via Rho and YAP

    Novev JK, Heltberg ML, Jensen MH, Doostmohammadi A. Spatiotemporal model of cellular mechanotrans- duction via Rho and YAP. Integrative Biology. 2021;13(8):197–209. doi:10.1093/intbio/zyab012

  12. [19]

    Structurally Governed Cell Mechanotransduction through Multiscale Modeling

    Kang J, Puskar KM, Ehrlicher AJ, LeDuc PR, Schwartz RS. Structurally Governed Cell Mechanotransduction through Multiscale Modeling. Scientific Reports. 2015;5(1):8622. doi:10.1038/srep08622

  13. [20]

    Pearling in cells: A clue to understanding cell shape

    Bar-Ziv R, Tlusty T, Moses E, Safran SA, Bershadsky A. Pearling in cells: A clue to understanding cell shape. Proceedings of the National Academy of Sciences. 1999;96(18):10140–10145. doi:10.1073/pnas.96.18.10140

  14. [21]

    Single Cells Spreading on a Pro- tein Lattice Adopt an Energy Minimizing Shape

    Vianay B, K¨ afer J, Planus E, Block M, Graner F, Guillou H. Single Cells Spreading on a Pro- tein Lattice Adopt an Energy Minimizing Shape. Physical Review Letters. 2010;105(12):128101. doi:10.1103/PhysRevLett.105.128101

  15. [22]

    Dynamics of cell shape and forces on micropatterned substrates predicted by a cellular Potts model

    Albert PJ, Schwarz US. Dynamics of cell shape and forces on micropatterned substrates predicted by a cellular Potts model. Biophysical Journal. 2014;106(11):2340–2352. doi:10.1016/j.bpj.2014.04.036

  16. [23]

    A Computational Model of YAP/TAZ Mechanosensing

    Sun M, Spill F, Zaman MH. A Computational Model of YAP/TAZ Mechanosensing. Biophysical Journal. 2016;110(11):2540. doi:10.1016/j.bpj.2016.04.040

  17. [24]

    Exploring the influence of cytosolic and mem- brane F AK activation on YAP/TAZ nuclear translocation

    Eroum´ e KS, Cavill R, Staˇ nkov´ a K, de Boer J, Carlier A. Exploring the influence of cytosolic and mem- brane F AK activation on YAP/TAZ nuclear translocation. Biophysical Journal. 2021;120(20):4360–4377. doi:10.1016/j.bpj.2021.09.009

  18. [25]

    Introduction to Linear Elasticity

    Gould PL. Introduction to Linear Elasticity. Springer New York; 2013. Available from:http://link. springer.com/10.1007/978-1-4614-4833-4

  19. [26]

    Finite element methods for surface PDEs

    Dziuk G, Elliott CM. Finite element methods for surface PDEs. Acta Numerica. 2013;22:289–396. doi:10.1017/S0962492913000056

  20. [28]

    Environmental stiff- ness restores mechanical homeostasis in vimentin-depleted cells

    Grolleman J, van Engeland NCA, Raza M, Azimi S, Conte V, Sahlgren CM, et al. Environmental stiff- ness restores mechanical homeostasis in vimentin-depleted cells. Nature, Scientific Reports. 2023;13:18374. doi:10.1038/s41598-023-44835-8

  21. [29]

    A theoretical model for focal adhesion and cytoskeleton formation in non-motile cells

    McNicol GR, Dalby MJ, Stewart PS. A theoretical model for focal adhesion and cytoskeleton formation in non-motile cells. Journal of Theoretical Biology. 2025;596:111965. doi:10.1016/j.jtbi.2024.111965

  22. [30]

    Cellular mechanotransduction in health and diseases: from molecular mechanism to therapeutic targets

    Di X, Gao X, Peng L, Ai J, Jin X, Qi S, et al. Cellular mechanotransduction in health and diseases: from molecular mechanism to therapeutic targets. Sig Transduct Target Ther. 2023;8:282. doi:10.1038/s41392-023- 01501-9

  23. [31]

    The cell as a material

    Kasza KE, Rowat AC, Liu J, Angelini TE, Brangwynne CP, Koenderink GH, et al. The cell as a material. Current Opinion in Cell Biology. 2007;19(1):101–107. doi:10.1016/j.ceb.2006.12.002

  24. [32]

    The cytoplasm of living cells behaves as a poroelastic material

    Moeendarbary E, Valon L, Fritzsche M, Harris AR, Moulding DA, Thrasher AJ, et al. The cytoplasm of living cells behaves as a poroelastic material. Nature Materials. 2013;12(3):253–261. doi:10.1038/nmat3517

  25. [33]

    Controlling cell–matrix traction forces by extracellular geometry

    Banerjee S, Marchetti MC. Controlling cell–matrix traction forces by extracellular geometry. New Journal of Physics. 2013;15(3):035015. doi:10.1088/1367-2630/15/3/035015

  26. [34]

    Geometry Regulates Traction Stresses in Adherent Cells

    Oakes PW, Banerjee S, Marchetti MC, Gardel ML. Geometry Regulates Traction Stresses in Adherent Cells. Biophysical Journal. 2014;107(4):825–833. doi:10.1016/j.bpj.2014.06.045

  27. [35]

    Reversible elastic phase field approach and application to cell monolayers

    Chojowski R, Schwarz US, Ziebert F. Reversible elastic phase field approach and application to cell monolayers. The European Physical Journal E. 2020;43(10):63. doi:10.1140/epje/i2020-11988-1. 16

  28. [36]

    Predicting YAP/TAZ Nuclear Translocation in Response to ECM Mechanosensing

    Cheng B, Li M, Wan W, Guo H, Genin GM, Lin M, et al. Predicting YAP/TAZ Nuclear Translocation in Response to ECM Mechanosensing. Biophysical Journal. 2023;122(1):43–53. doi:10.1016/j.bpj.2022.11.2943

  29. [37]

    Mechanotransduction and nuclear function

    Graham DM, Burridge K. Mechanotransduction and nuclear function. Current Opinion in Cell Biology. 2016;40:98–105. doi:10.1016/j.ceb.2016.03.006

  30. [39]

    Focal adhesions, stress fibers and mechanical tension

    Burridge K, Guilluy C. Focal adhesions, stress fibers and mechanical tension. Experimental Cell Research. 2016;343:14–20. doi:10.1016/j.yexcr.2015.10.029

  31. [40]

    Local 3D matrix microenvironment regulates cell migration through spatiotemporal dynamics of contractility-dependent adhesions

    Doyle AD, Carvajal N, Jin A, Matsumoto K, Yamada KM. Local 3D matrix microenvironment regulates cell migration through spatiotemporal dynamics of contractility-dependent adhesions. J Cell Sci. 2015;6:8720. doi:10.1038/ncomms9720

  32. [41]

    Force activates smooth muscleα-actin promoter activity through the Rho signaling pathway

    Zhao XH, Laschinger C, Arora P, Sz´ aszi K, Kapus A, McCulloch CA. Force activates smooth muscleα-actin promoter activity through the Rho signaling pathway. J Cell Sci. 2007;120:1801–1809. doi:10.1242/jcs.001586

  33. [42]

    The Poisson Ratio of the Cellular Actin Cortex Is Frequency Dependent

    Mokbel M, Hosseini K, Aland S, Fischer-Friedrich E. The Poisson Ratio of the Cellular Actin Cortex Is Frequency Dependent. Biophysical Journal. 2020;118(8):1968–1976. doi:10.1016/j.bpj.2020.03.002

  34. [43]

    Rho-stimulated contractility drives the formation of stress fibers and focal adhesions

    Chrzanowska-Wodnicka M, Burridge K. Rho-stimulated contractility drives the formation of stress fibers and focal adhesions. Journal of Cell Biology. 1996;133(6):1403–1415. doi:10.1083/jcb.133.6.1403

  35. [44]

    Engineered extracellular matrices with controlled mechanics modulate renal proximal tubular cell epithelialization

    Beamish JA, Chen E, Putnam AJ. Engineered extracellular matrices with controlled mechanics modulate renal proximal tubular cell epithelialization. PLOS ONE. 2017;12(7):e0181085. doi:10.1371/journal.pone.0181085

  36. [45]

    Cell shape provides global control of focal adhesion assembly

    Chen CS, Alonso JL, Ostuni E, Whitesides GM, Ingber DE. Cell shape provides global control of focal adhesion assembly. Biochemical and Biophysical Research Communications. 2003;307(2):355–361. doi:10.1016/S0006- 291X(03)01165-3

  37. [46]

    Cell Shape, Cytoskeletal Tension, and RhoA Regulate Stem Cell Lineage Commitment

    McBeath R, Pirone DM, Nelson CM, Bhadriraju K, Chen CS. Cell Shape, Cytoskeletal Tension, and RhoA Regulate Stem Cell Lineage Commitment. Developmental Cell. 2004;6(4):483–495. doi:10.1016/S1534- 5807(04)00075-9

  38. [47]

    Biomechanical Homeostasis

    Kassab GS. Biomechanical Homeostasis. In: Kassab GS, editor. Coronary Circulation: Mechanobiology, Growth, Remodeling, and Clinical Implications. Cham: Springer International Publishing; 2024. p. 1–43. Available from:https://doi.org/10.1007/978-3-031-62652-4_1

  39. [48]

    Cell size and actin architec- ture determine force generation in optogenetically activated cells

    Andersen T, W¨ orthm¨ uller D, Probst D, Wang I, Moreau P, Fitzpatrick V, et al. Cell size and actin architec- ture determine force generation in optogenetically activated cells. Biophysical Journal. 2023;122(4):684–696. doi:10.1016/j.bpj.2023.01.011

  40. [49]

    Cellular adaptation to biomechanical stress across length scales in tissue homeostasis and disease

    Weaver VM, Gilbert PM. Cellular adaptation to biomechanical stress across length scales in tissue homeostasis and disease. Seminars in cell & developmental biology. 2016;67:141. doi:10.1016/j.semcdb.2016.09.004

  41. [50]

    Cellular stiffness sensing through talin 1 in tissue mechanical homeostasis

    Chanduri M, Kumar A, Weiss D, Emuna N, Barsukov I, Shi M, et al. Cellular stiffness sensing through talin 1 in tissue mechanical homeostasis. Science Advances. 2024;10(34):eadi6286. doi:10.1126/sciadv.adi6286

  42. [51]

    Characterizing and Engineering Biomimetic Materials for Viscoelastic Mechanotransduction Studies

    Cacopardo L, Guazzelli N, Ahluwalia A. Characterizing and Engineering Biomimetic Materials for Viscoelastic Mechanotransduction Studies. Tissue Engineering Part B, Reviews. 2022;28(4):912–925. doi:10.1089/ten.TEB.2021.0151

  43. [52]

    Endothelial actin and cell stiffness is modulated by substrate stiffness in 2D and 3D

    Byfield FJ, Reen RK, Shentu TP, Levitan I, Gooch KJ. Endothelial actin and cell stiffness is modulated by substrate stiffness in 2D and 3D. Journal of Biomechanics. 2009;42(8):1114–1119. doi:10.1016/j.jbiomech.2009.02.012

  44. [53]

    Substrate stress relaxation regulates cell spreading

    Chaudhuri O, Gu L, Darnell M, Klumpers D, Bencherif SA, Weaver JC, et al. Substrate stress relaxation regulates cell spreading. Nature Communications. 2015;6(1):6365. doi:10.1038/ncomms7365

  45. [54]

    Directed Cell Migration towards Softer Environments

    Isomursu A, Park K-Y, Hou J, Cheng B, Mathieu M, Shamsan GA, et al. Directed Cell Migration towards Softer Environments. Nature Materials. 2022;21(9):1081–1090. doi:10.1038/s41563-022-01294-2. 17

  46. [55]

    Physical Aspects of Cell Culture Substrates: Topography, Roughness, and Elasticity

    Ross AM, Jiang Z, Bastmeyer M, Lahann J. Physical Aspects of Cell Culture Substrates: Topography, Roughness, and Elasticity. Small. 2012;8(3):336–355. doi:10.1002/smll.201100934

  47. [56]

    Insights gained from computational modeling of YAP/TAZ signaling for cellular mechanotransduction

    Jafarinia H, Khalilimeybodi A, Barrasa-Fano J, Fraley SI, Rangamani P, Carlier A. Insights gained from computational modeling of YAP/TAZ signaling for cellular mechanotransduction. npj Systems Biology and Applications. 2024;10(1):1–14. doi:10.1038/s41540-024-00414-9

  48. [57]

    The residence time of focal adhesion kinase (F AK) and paxillin at focal adhesions in renal epithelial cells is determined by adhesion size, strength and life cycle status

    Le D´ ev´ edec SE, Geverts B, de Bont H, Yan K, Verbeek FJ, Houtsmuller AB, et al. The residence time of focal adhesion kinase (F AK) and paxillin at focal adhesions in renal epithelial cells is determined by adhesion size, strength and life cycle status. Journal of Cell Scien...

  49. [58]

    Automated Solution of Differential Equations by the Finite Element Method

    Logg A, Mardal KA, Wells G, editors. Automated Solution of Differential Equations by the Finite Element Method. vol. 84 of Lecture Notes in Computational Science and Engineering. Berlin, Heidelberg: Springer Berlin Heidelberg; 2012. Available from:http://link.springer.com/10.1...

  50. [59]

    Gmsh: A 3-D finite element mesh generator with built-in pre- and post- processing facilities

    Geuzaine C, Remacle JF. Gmsh: A 3-D finite element mesh generator with built-in pre- and post- processing facilities. International Journal for Numerical Methods in Engineering. 2009;79(11):1309–1331. doi:10.1002/nme.2579

  51. [60]

    Implicit-explicit timestepping with finite element approximation of reaction-diffusion systems on evolving domains

    Lakkis O, Madzvamuse A, Venkataraman C. Implicit-explicit timestepping with finite element approximation of reaction-diffusion systems on evolving domains. SIAM Journal on Numerical Analysis. 2013;51(4):2309–2330. doi:10.1137/120880112. 18 S1 Fig (A) ρa ϕd ϕa |u| (C) ρa ϕd ϕa ...

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.