🔬 condensed matter

Cell-induced wrinkling patterns on soft substrates

This paper presents a minimal mechanical framework combining Föppl-von Kármán equations with a phase-field model to explain how cellular contractile forces induce specific wrinkling patterns on soft substrates, successfully predicting pattern transitions and scaling laws that align with experimental observations.

Original authors: Aleksandra Ardaševa, Varun Venkatesh, Daiki Matsunaga, Shinji Deguchi, Amin Doostmohammadi

Published 2026-03-16

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Aleksandra Ardaševa, Varun Venkatesh, Daiki Matsunaga, Shinji Deguchi, Amin Doostmohammadi

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a single cell sitting on a soft, squishy gel, like a tiny jellybean resting on a trampoline. This isn't just a passive rest; the cell is actively pulling on the gel with its internal muscles (made of proteins called actin and myosin).

This paper is about understanding what happens to the trampoline when the jellybean pulls on it.

Here is the breakdown of the research using simple analogies:

1. The Big Picture: The Cell as a "Wrinkle-Maker"

When you step on a soft rug, it bunches up around your feet. Similarly, when a cell pulls on a soft surface, it creates wrinkles. These aren't just random folds; they are organized patterns that can stretch far beyond the cell itself.

The scientists wanted to figure out the "rules of the game":

How far do these wrinkles travel?
Do they look like spokes on a wheel, or lines in a field?
How does the softness of the gel or the shape of the cell change the pattern?

2. The Tool: A "Virtual Lab"

Instead of just watching cells under a microscope (which is hard to measure precisely), the team built a computer simulation. Think of this as a video game where they can control every variable:

The Cell: They modeled it as a drop of "active liquid crystal." Imagine a drop of oil that has its own internal compass and wants to squeeze itself.
The Gel: They modeled it as a thin, stretchy sheet sitting on a thicker, squishy base.
The Physics: They used complex math (the Föppl-von Kármán equations) to predict exactly how the sheet would buckle and fold when the cell pulled on it.

3. The Key Discoveries

A. The "Stress Hotspots" (The Cell's Feet)

The simulation showed that the cell doesn't pull evenly. It pulls hardest at its ends (like the tips of an oval).

Analogy: Imagine holding a rubber band. You don't pull the middle; you pull the ends. Because the cell pulls hard at its tips, that's where the wrinkles start. These wrinkles then shoot out into the surrounding gel, acting like a "mechanical shout" that tells other cells, "I'm here!"

B. The Shape Matters (Round vs. Long)

The shape of the cell changes the pattern of the wrinkles.

Round Cell: If the cell is a perfect circle, the wrinkles spread out evenly in all directions, like ripples in a pond.
Long Cell: If the cell is stretched out (like a football), the wrinkles become more directional.
The Twist: Surprisingly, the simulation found that circular cells actually create wrinkles that travel the farthest. Elongated cells create wrinkles that don't reach as far. It's like a round stone creating bigger ripples in a pond than a long stick does.

C. The "Muscle Tension" (How Hard the Cell Pulls)

The cell can change how hard it pulls (its contractility). The researchers found three distinct "modes" of behavior:

The Calm Mode: The cell pulls steadily. The wrinkles grow smoothly and predictably.
The Wobbly Mode: The cell pulls a bit too hard, and its internal structure starts to wobble (like a jelly shaking). This actually stops the wrinkles from growing as far, even though the cell is pulling harder. It's like trying to walk on a floor that's shaking too much; you can't get a good grip.
The Chaotic Mode: The cell pulls so hard it goes into "turbulence" (like a stormy sea). Surprisingly, the wrinkles start growing again, but they are driven by the chaotic energy of the whole system rather than a steady pull.

D. The "Magic Mirror" (Matching Reality)

Finally, the team tested their computer model against real experiments. They took photos of real cells on real gels and fed the data into their simulation.

The Result: The computer predicted the wrinkles almost perfectly.
The Cool Part: They showed that you don't even need to know the exact forces the cell is using. If you just tell the computer the shape of the cell, it can guess the forces and predict the wrinkles. It's like looking at a shadow and being able to tell exactly what object cast it.

Why Does This Matter?

Cells talk to each other not just with chemicals, but with mechanics. By creating wrinkles, a cell is essentially sending a message across the gel to its neighbors.

The "Sensing Radius": The distance the wrinkles travel is like a cell's "hearing range." If the wrinkles travel far, the cell can "hear" (feel) what its neighbors are doing from a distance.
Medical Applications: This helps us understand how tissues form, how wounds heal, and how cancer cells might spread. If we know how to control these wrinkles, we might be able to guide cells to grow in specific ways or stop them from misbehaving.

In short: This paper built a virtual world to show us that cells are like tiny architects. By pulling on their environment, they create a landscape of wrinkles that serves as a map for communication, and the shape and strength of the cell dictate exactly what that map looks like.

1. Problem Statement

Cells exert traction forces on compliant substrates, inducing surface instabilities that manifest as characteristic wrinkling patterns. While these patterns provide crucial environmental cues for cell sensing, migration, and differentiation, the quantitative relationship between active, cell-generated forces and the emergent substrate folding patterns remains poorly understood.

Existing theoretical models often rely on idealized, uniform stress fields or simplified force dipoles, which fail to capture:

The spatial heterogeneity and temporal dynamics of forces generated by real cells (e.g., actin-myosin contractility).
The non-linear mechanics of substrate folding.
The feedback loop where cell forces create wrinkles, which in turn influence cell behavior and force distribution.

There is a critical need for a computational framework that can self-consistently predict detailed wrinkling patterns based on the underlying cellular mechanics and active stress generation.

2. Methodology

The authors developed a coupled computational modeling framework that integrates two distinct physical models:

A. Active Nematic Phase-Field Model for the Cell

Representation: The cell is modeled as an active nematic droplet using a scalar phase field ( $\phi$ ) to define the cell boundary and a nematic tensor ( $Q$ ) to describe internal cytoskeletal order.
Dynamics: The cell shape is governed by Cahn-Hilliard free energy, while the internal stress is driven by active nematic hydrodynamics (Beris-Edwards equations coupled with Navier-Stokes).
Stress Generation: The total cell stress ( $\sigma_{ij}^{cell}$ ) includes contributions from fluid pressure, the phase-field interface, and dipolar activity ( $\zeta$ ). The activity term drives contractile forces, with topological defects (+1/2 defects) naturally forming at the cell poles, creating localized high-stress regions.

B. Föppl-von Kármán (FvK) Thin Plate Model for the Substrate

Setup: The substrate is modeled as a bilayer system: a thin elastic film (plasma-irradiated PDMS) on a thicker viscoelastic foundation.
Governing Equations: The deformation of the substrate ( $\omega_x, \omega_y, \omega_z$ ) is described using the FvK equations for thin plate deformation.
Coupling: The in-plane contractile stresses generated by the active nematic cell model are applied as boundary conditions/forcing terms to the FvK equations. This allows the simulation of how in-plane cellular forces translate into out-of-plane wrinkling ( $\omega_z$ ).

Numerical Implementation:

The framework was implemented in OpenFOAM using the finite-volume method.
Simulations were run until statistical steady-state was reached, varying parameters such as substrate stiffness, cell aspect ratio, and activity levels.

3. Key Contributions

Unified Framework: First to couple an active nematic phase-field model (capturing internal cell dynamics and defect-driven stress) with FvK thin plate theory to simulate cell-induced wrinkling.
Mechanism Elucidation: Identified that topological defects at the cell poles are the primary drivers of localized stress concentration, which initiates wrinkling.
Validation: Validated the model against experimental Wrinkle-Force Microscopy (WFM) data, demonstrating that the model can reproduce experimental patterns using either measured stress fields or solely cell geometry as input.

4. Key Results

A. Pattern Formation and Propagation

Stress Localization: High stress concentrates at the elongated ends (poles) of the cell due to nematic anchoring and defect formation.
Long-Range Effects: Wrinkles propagate significantly beyond the cell footprint, defining an effective "sensing radius" for mechanical communication.
Amplitude Profile: Wrinkle amplitude peaks near the cell boundary and decays both inward and outward.

B. Influence of Substrate Properties

Stiffness ( $\mu_f$ ) & Poisson's Ratio ( $\nu_f$ ): The extent of wrinkling is highly sensitive to the film's shear modulus and Poisson's ratio.
- Wrinkling extent decreases linearly with increasing film Poisson's ratio.
- Extent shows a shallow maximum with respect to film stiffness before decreasing.
Substrate Independence: The underlying substrate's Poisson ratio ( $\nu_s$ ) has minimal effect on the propagation distance compared to the thin film properties.

C. Influence of Cell Shape

Aspect Ratio: Circular cells produce the maximum wrinkling extent. As cells become more elongated (ellipsoidal), the extent of wrinkles decreases.
Pattern Morphology:
- Circular cells: Produce straight, parallel wrinkles.
- Elongated cells: Produce zig-zag patterns underneath the cell.
Feedback: This suggests cell geometry acts as a powerful regulator of the mechanical environment, creating a feedback loop where substrate mechanics shape cell geometry, which in turn dictates the resulting pattern.

D. Influence of Cell Contractility (Activity $\zeta$ )
The system exhibits three distinct dynamical regimes as contractility increases:

Stationary Nematic Phase: Low activity. Static director field; wrinkle amplitude and extent grow controllably with activity.
Bend Wave Regime: Intermediate activity. The nematic field develops propagating bend waves. This introduces temporal fluctuations that disrupt coherent force transmission, causing wrinkle extent to stagnate despite increasing forces.
Active Turbulence: High activity. Chaotic spatiotemporal dynamics. Wrinkle extent initially drops (due to motile defects) but resumes growth at very high activity levels as time-averaged stresses become sufficient to drive propagation.

E. Experimental Validation

The model successfully reproduced experimental wrinkle patterns from fibroblast cells on PDMS gels.
Predictive Power: When fed only the cell shape (without pre-measured stress fields), the active nematic model generated stress distributions that accurately predicted the orientation and spatial organization of the observed wrinkles.

5. Significance

Mechanical Communication: The study establishes that cells can mechanically communicate over distances much larger than their own size via substrate wrinkling, with the "sensing radius" tunable by cell shape and substrate stiffness.
Minimal Modeling Framework: Provides a robust tool to interpret wrinkling assays, allowing researchers to infer underlying cellular forces from observed surface patterns (inverse problem) or predict patterns from known cell geometries (forward problem).
Biological Implications: Highlights the critical role of cell shape and internal active dynamics (nematic order) in determining how cells remodel their microenvironment. This has implications for understanding tissue morphogenesis, collective cell migration, and mechanotransduction.
Future Directions: The framework lays the groundwork for studying multi-cell interactions, dynamic cell migration, and heterogeneous substrates, bridging the gap between single-cell mechanics and tissue-level organization.