Optogenetic control of intracellular flows and cell migration: a comprehensive mathematical analysis with a minimal active gel model

Oliver M. Drozdowski,Falko Ziebert,Ulrich S. Schwarz

DOI: https://doi.org/10.1103/PhysRevE.104.024406

2021-07-06

Abstract:The actin cytoskeleton of cells is in continuous motion due to both polymerization of new filaments and their contraction by myosin II molecular motors. Through adhesion to the substrate, such intracellular flow can be converted into cell migration. Recently, optogenetics has emerged as a new powerful experimental method to control both actin polymerization and myosin II contraction. While optogenetic control of polymerization can initiate cell migration by generating protrusion, it is less clear if and how optogenetic control of contraction can also affect cell migration. Here we analyze the latter situation using a minimal variant of active gel theory into which we include optogenetic activation as a spatiotemporally constrained perturbation. The model can describe the symmetrical flow of the actomyosin system observed in optogenetic experiments, but not the long-lasting polarization required for cell migration. Motile solutions become possible if cytoskeletal polymerization is included through the boundary conditions. Optogenetic activation of contraction can then initiate locomotion in a symmetrically spreading cell and strengthen motility in an asymmetrically polymerizing one. If designed appropriately, it can also arrest motility even for protrusive boundaries.

Subcellular Processes,Soft Condensed Matter

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **How to use optogenetics to control cell contractility to influence cell migration**.

Specifically, the authors are concerned with whether and how the contractile activity of myosin II can be controlled by optogenetic means to affect cell migration. Although optogenetic control of actin polymerization can trigger cell migration, relatively few studies have been conducted on how optogenetic control of contractility affects cell migration. In this paper, by establishing a minimal active gel model, the influence of activating contractility on intracellular flow and cell migration under specific spatio - temporal conditions was explored.

### Main contributions of the paper:
1. **Theoretical framework**: A simplified active gel model was used. This model does not consider the local polarization of the cytoskeleton but focuses on the velocity field driven by local active stress.
2. **Model analysis**: Through mathematical analysis of the model, the influence of optogenetic activation as a spatio - temporally restricted perturbation on active stress was studied.
3. **Boundary conditions**: Elastic boundary conditions were introduced to simulate the situation where cell size is controlled by elasticity - related processes.
4. **Numerical simulation**: Numerical simulations were carried out by the Finite Volume Method to verify the predictions of the model and show cell behaviors under different optogenetic activation intensities and spatio - temporal distributions.

### Main findings:
- **Symmetric flow**: The model can describe the symmetric flow of the actomyosin system observed in optogenetic experiments, but it cannot describe the long - term polarization required for cell migration.
- **Cell migration**: When cytoskeletal polymerization is introduced into the boundary conditions, optogenetic activation of contractility can initiate the movement of symmetrically expanding cells and enhance the motility of asymmetrically polymerizing cells. With appropriate design, optogenetic activation can even stop the movement of cells with protrusive boundaries.

### Conclusion:
This study theoretically explored the influence of optogenetic control of cell contractility on cell migration, providing a new perspective for understanding the mechanisms of cell migration. Through mathematical modeling and numerical simulation, the authors revealed the potential role of optogenetic activation in cell migration, providing a theoretical basis for further experimental research.

### Formula summary:
- **Relationship between total stress and strain**:
\[
\eta\frac{\partial\epsilon}{\partial t}=\left[1 + \tau\frac{\partial}{\partial t}+\tau v\frac{\partial}{\partial x}\right](\sigma-\sigma_{\text{act}})
\]
where \(\eta\) is the viscosity, \(\tau = \frac{\eta}{E}\) is the Maxwell relaxation time, \(E\) is the elastic modulus, \(\sigma\) is the total stress, and \(\sigma_{\text{act}}\) is the active stress.

- **Relationship between stress and friction**:
\[
\frac{\partial\sigma}{\partial x}=\xi v
\]
where \(\xi\) is the friction coefficient and \(v\) is the flow velocity.

- **Dimensionless equation**:
\[
L^{2}\frac{\partial^{2}\sigma}{\partial x^{2}}-T\frac{\partial\sigma}{\partial t}-T\left(\frac{\partial\sigma}{\partial x}\right)^{2}-\sigma = -\sigma_{\text{act}}
\]
where \(L=\sqrt{\frac{\eta}{\xi L_{0}^{2}}}\) is the relative length scale, \(T = \frac{k\tau}{\xi L_{0}^{2}}=\frac{k\e