The variational modulus density theory explains mechanical responses of cell membranes and membrane crosslinkers

TL;DR

This paper introduces a new theory that combines classical and quantum mechanics to better simulate cell membranes and their interactions.

Contribution

A novel continuum theory that integrates classical mechanics with quantum wavefunction compatibility for modeling cell membranes.

Findings

The theory explains confined Brownian motion of cell membrane inclusions via elastic crosslinkers.

Calculations using the model matched experimental data on adhesion molecular machineries and cell membranes.

The model offers improved simulation capabilities by interpreting the quantum wavefunction physically and mathematically.

Abstract

Both classical mechanics and quantum mechanics explain Brownian motion. However, it remains unclear whether they are compatible with each other, as the physical and mathematical identity of the wavefunction in quantum mechanics has been elusive. In this work, a continuum theory using grammars in classical mechanics modeling, but potentially compatible with the quantum wavefunction, is introduced. The theory explains the confined Brownian motion of cell membrane inclusions interacting with extracellular matrices or cytoskeletons via elastic molecular crosslinkers. This crosslinker theory is integrated into the Canham-Helfrich-Evans model for fluid membranes. Calculations, based on a finite element method for the combined theory, reproduced measured data from adhesion molecular machineries and cell membranes. Overall, by providing physical and mathematical interpretations of the quantum wavefunction, the presented theoretical model provides improved capabilities for the realistic simulation of cell membranes and membrane linker proteins.

Figures

Descriptions of the membrane-crosslinker model. (A) Two different views in mathematical modeling for the confined Brownian motion of crosslinkers based on conventional elasticity (left) and the modulus density theory proposed in this work (right). (A, left) The conventional viewpoint. Similar to an elastic spring on the solid membrane, the crosslinker force at a fixed time point is calculated using a constant elastic modulus. (A, right) The alternative viewpoint (in this work). As a modeling treatment, the average distribution of the mobile crosslinker was considered to calculate the crosslink

Lipid sorting and formation of lipid nanodomains. (A) Three-dimensional deformation of a planar membrane with the surface map for sorting of eight different lipid classes with varying bending moduli. Crosslinkers are not shown. Left: Angled view. Right: Planar view. (B) Lipid number strain α profile for the membrane in (A). Crosslinkers are shown here. (C‒J) The pattern of lipid nanodomains, shown separately for each lipid class. (K‒N) Expanded views of patterns in (D‒G). (O) Plots of the mean and standard deviation (SD) of the area of the nanodomains. Mean values are 20,183 $\documentclass[12

Calculations of the interaction among integrins, talin, and membranes in the formation of nascent adhesions. (A) Force vs. extension responses with respect to inputs in Fig. S3. The extension of the crosslinker was defined from its fixed point (red bead) to the maximum $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document}$ point on the membrane (black bead). Expanded views of indicated regions (a

The effect of the number and distribution of crosslinkers and lipid composition on nanomechanical responses of membrane-crosslinker complexes. (A) The deformed configuration with the surface strain map for a $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.944\times 1.944 \; \mathrm{\mu {m}^{2}}$$\end{document}$ membrane ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfon

Application of the modulus density theory to the crosslinker with two fixed ends. Calculations using Eq. 7 with $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_\mathrm{fixed-talin}=1\times {10}^{-7}\sqrt{\mathrm{N}}$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage

A combined analysis for the force transfer at the membrane-crosslinker interface using the modulus density theory. (A) The $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.944\times 1.944 \; \mathrm{{\mu m}^{2}}$$\end{document}$ membrane was stimulated by using four different complexes of crosslinkers and membrane inclusions ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{am

Full text

Introduction

Understanding interactions between cell membranes and crosslinkers that connect the membrane to extracellular matrices (ECMs) or cytoskeletons is crucial for studying cellular signaling across cell membranes. A variety of adhesion receptors and cytoplasmic adaptors that act as crosslinkers have been identified. For example, activated integrins inserted into the membrane mechanically interact with ECMs to facilitate migration, gene expression, and homeostasis^1^. The talin adaptor, which interacts with membrane-inserted integrins, forms physical connections between membranes and cytoskeletons to support various physiological processes^2,3^. Similarly, CD44 membrane receptors are linked to underlying cytoskeletons to mediate cellular signaling^4,5^. Despite rapid advances in experimental technologies and the accumulation of a vast amount of measured data on membrane-crosslinker interaction, current theoretical continuum models do not fully incorporate these experimental developments.

The membrane-crosslinker interaction is complex. One end of the crosslinker is linked to a membrane-inclusion molecule (or domain) inserted into the membrane characterized by a certain lipid environment, while the other end is fixed on ECMs or cytoskeletons. Therefore, one end of the crosslinker exhibits confined Brownian motion on the surface of a deforming fluid membrane. At the same time, the spring-like crosslinker can be stretched and relaxed due to applied forces. In addition to these physical complexities, there exists a viewpoint that considers the Brownian motion as a quantum phenomenon^6–8^, while cell membranes have been explained by classical continuum mechanics^9–11^. A unified theoretical and mathematical framework for classical mechanics and quantum mechanics is thus required for realistic calculations of the membrane-crosslinker complex.

In this work, two simple yet novel principles are proposed to incorporate all these complexities and discrepancies into a single variational framework. First, the average of multiple configurations of the crosslinker under confined diffusion can be considered to calculate the crosslinker force (Fig. 1A, right). It is assumed that, during confined diffusion, forces at the fixed end of the crosslinker are constant if all the other boundaries in the membrane-crosslinker complex are also fixed (even though the mobile end is under confined diffusion). This assumption aligns with the notion that membrane curvature change and lateral stretching are negligible during the confined diffusion, because a mechanical change of the membrane can alter the crosslinker force. In this case, the constant crosslinker force can be calculated by finding the average configuration of the crosslinker under the confined diffusion, because the single constant force representing all individual crosslinker configurations is the same as the averaged force, even though the new configuration after the averaging process may differ from those individual configurations. In continuum theory, finding the average configuration of the mobile crosslinker is equivalent to finding the modulus density profile on the membrane that is used to determine the averaged distribution of internal elastic energies from the individual configurations of the crosslinker under the confined diffusion (Fig. 1B). The approach here is, therefore, statistical and independent of the time scale of confined diffusion, even though confined diffusion itself is a time-dependent process.Fig. 1. Descriptions of the membrane-crosslinker model. (A) Two different views in mathematical modeling for the confined Brownian motion of crosslinkers based on conventional elasticity (left) and the modulus density theory proposed in this work (right). (A, left) The conventional viewpoint. Similar to an elastic spring on the solid membrane, the crosslinker force at a fixed time point is calculated using a constant elastic modulus. (A, right) The alternative viewpoint (in this work). As a modeling treatment, the average distribution of the mobile crosslinker was considered to calculate the crosslinker force. (B‒D) An example configuration calculated using the finite element model of membrane-crosslinker complexes. (B) The elastic modulus density is shown in the surface map. (C) Lateral lipid number strain of the membrane is shown in the surface map. (D) The bending moduli of different lipid classes in the membrane are shown in the surface map. (E) Description of how to index the group of crosslinkers in this model. (F) Description of how the constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} can be split into real numbers assigned over a parametric surface using two degrees of freedom (DOFs) at each point. The functions \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} determine values for the two DOFs. With two DOFs, both “positive” and “negative” values can be independently considered at each point during the decomposition of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} . Such a mathematical strategy may maximize the variability of the shape of the modulus density function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} made from the given \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} value. The illustration here depicts the region where the modulus density is non-zero (except the boundary).

Second, another assumption is made that any confined Brownian motion of a mobile membrane inclusion linked to a crosslinker and embedded within a certain lipid environment (i.e., all possible elastic modulus density profiles) is a physical output that depends on a certain intrinsic physical constant of the inclusion-crosslinker complex. Therefore, a combined mathematical framework is introduced to generate various functions with values greater than or equal to zero (i.e., various configurations for the modulus density profile) on a two-dimensional surface as many as possible, by using a single constant value. This is achieved by defining the modulus density value at each point with two degrees of freedom (DOFs) (using functions \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} , see Fig. 1F). The sum of squares of these two \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} values is the modulus density at each point ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}={{K}_{1c}}^{2}+{{K}_{2c}}^{2}$$\end{document} , see Fig. 1F and Eq. 6 in Results). To correlate the two-dimensional functions \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} with the single constant value, it is necessary to employ the integral operator. The area integral of the functions serves as the intrinsic constant for the inclusion-crosslinker complex embedded in a certain lipid environment ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\int \left({K}_{1c}+{K}_{2c}\right)dA={\varphi }_{c}$$\end{document} , see Fig. 1F, the relation is implemented in Eq. 1).

There must be many configurations for the function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} that can yield the same constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} when the function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} (without \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} ) is integrated over the surface. However, “plus” and “minus” signs are two fundamental operators in integration. To enhance the variability of the shape of the function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} made from the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} constant, therefore, two DOFs using \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} are required so that both the “plus” and “minus” operations can be independently considered at each point on the surface in the integration process leading to the constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} . In other words, the two DOFs framework allows the most dynamic decomposition of the single constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} into real numbers distributed over the continuous membrane surface. Having both the “positive” and “negative” components at each point still guarantees the modulus density value that is greater than or equal to zero. Additional DOFs beyond these two at each point might be redundant in considering the two operations, and thus may not affect the final results. With the foregoing treatments using \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} , therefore, the maximum variability of the shape of the modulus density function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} , with respect to the given constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} , might be achieved. Finally, we can generalize that all physically possible confined diffusion configurations (i.e., all physically possible \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} functions) for a certain membrane inclusion, linked to a certain crosslinker and interacting with a certain lipid environment, might be generated using this mathematical framework. Among the infinitely many possibilities, however, the principle of stationary total potential energy provides one unique \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} profile.

This crosslinker theory is coupled with continuum models for lipid membranes^11–15^. Calculations, performed using a finite element method for the combined model (Fig. 1B‒D), directly reproduced numerous experimental data on adhesion molecular machineries and cell membranes. Furthermore, the model provides biological insight into the following: how cell membranes can be compartmentalized with a minimal number of membrane-cytoskeleton crosslinkers; how the presence of different lipid classes and their sorting affect the size of lipid reservoirs during mechanical responses of membranes; and how the new theory offers a more general framework for explaining one-dimensional crosslinker protein compared to conventional elasticity. Overall, the presented theoretical model provides invaluable biological insights into cell membranes and their interactions with crosslinkers.

Results

A combined continuum theory for cell membranes and membrane crosslinkers

The total functional \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Pi }^\mathrm{tot}$$\end{document} , which describes mechanical responses of a cell membrane and elastic crosslinkers connected to mobile membrane inclusions, is defined in Eq. 1.

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Pi }^\mathrm{tot}=\int {E}^\mathrm{curv}dA+\int {E}^\mathrm{st}dA+{\sum }_{c=1}^{{n}_\mathrm{crs}}\int {E}_{c}^\mathrm{crs}dA+{\uplambda }^\mathrm{lipid}\left(\int \frac{{\phi }_\mathrm{0}}{\alpha +1}dA-{n}_\mathrm{lipid}\right)+{\sum }_{c=1}^{{n}_\mathrm{crs}}{\uplambda }_{c}^\mathrm{crs}\left(\int \left({K}_{1c}+{K}_{2c}\right)dA-{\varphi }_{c}\right)$$\end{document}

where

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${E}^\mathrm{curv}=2{\kappa}{H}^{2}+{k}_\mathrm{g}G\quad\text{ for }\quad0\le \alpha \le {\alpha }_\mathrm{cross}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${E}^\mathrm{st}=\frac{{\upsigma }_{0}}{{c}_{1}}\mathrm{exp}\left({c}_{1}\alpha \right)-\frac{{\upsigma }_{0}}{{c}_{1}}\quad\text{ for }\quad0\le \alpha \le {\alpha }_\mathrm{cross}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${E}^\mathrm{st}=\frac{{K}_\mathrm{app}{\alpha }^{2}}{2}-{K}_\mathrm{app}{\alpha }_\mathrm{cut}\alpha +{c}_{2}\quad\text{ for }\quad\alpha>{\alpha }_\mathrm{cross}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${E}_{c}^\mathrm{crs}=\frac{1}{3}{M}_{c}{\sqrt{{\left(X-{\overline{X} }_{c}\right)}^{2}+{\left(Y-{\overline{Y} }_{c}\right)}^{2}+{\left(Z-{\overline{Z} }_{c}\right)}^{2}}}^{3}+{\sum }_{p=1}^{{n}_{c}^\mathrm{crs-add}}\frac{1}{3}{{G}_{cp}^\mathrm{M}M}_{c}{\sqrt{{\left(X-{\overline{X} }_{cp}\right)}^{2}+{\left(Y-{\overline{Y} }_{cp}\right)}^{2}+{\left(Z-{\overline{Z} }_{cp}\right)}^{2}}}^{3}$$\end{document}

and

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}={{K}_{1c}}^{2}+{{K}_{2c}}^{2}$$\end{document}

In Eq. 1, the first term represents the membrane energy generated from the mean curvature \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H$$\end{document} and the Gaussian curvature \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G$$\end{document} that can be expressed with displacement functions \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z$$\end{document} in parametric coordinates \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$u$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v$$\end{document} ^12,13,16^ (Fig. 1C and Fig. S1 in Supplementary Material (SM)). The second term accounts for energy due to the lateral strain \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} of membranes^11,14,15^ (Fig. 1C). Details on membrane modeling are described in the APPENDIX and SM. Briefly, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\alpha }_\mathrm{cross}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\alpha }_\mathrm{cut}$$\end{document} are the crossover and cut-off strains, respectively. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\kappa}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_\mathrm{g}$$\end{document} are the bending modulus and the Gaussian curvature modulus. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_\mathrm{app}$$\end{document} is the apparent area stretching modulus. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${c}_{1}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${c}_{2}$$\end{document} are constants (see APPENDIX). \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$dA$$\end{document} is the area element. Throughout this article, all calculations were performed in the range \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0\le \alpha \le {\alpha }_\mathrm{cross}$$\end{document} . The Gaussian curvature energy contribution in Eq. 2 was omitted in this work based on the Gauss-Bonnet theorem. A summary of the membrane parameters used in this work is provided in Table S1 in SM.

The expression in the third term of Eq. 1 denotes elastic energies from the crosslinkers that are mobile with the membrane inclusion at one end and fixed in the three-dimensional space at the other end (Fig. 1B and C). The energy density for a group of crosslinkers is written as given in Eq. 5. It gives tensile forces in a quadratic form when both ends of a single crosslinker are fixed (see Eq. 7). \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c$$\end{document} indexes the group of crosslinkers that share the membrane-targeting point i.e., the membrane inclusion (Fig. 1E). \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c$$\end{document} also globally indexes the first crosslinker in each group. It is possible that the group is composed of a single crosslinker (Fig. 1B and C). \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p$$\end{document} indexes additional crosslinkers within the group indexed by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c$$\end{document} . \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_\mathrm{crs}$$\end{document} is the number of crosslinker groups (Fig. 1E). \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{c}^\mathrm{crs-add}$$\end{document} is the number of additional crosslinkers in the group indexed with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c$$\end{document} . \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\overline{X} }_{c}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\overline{Y} }_{c}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\overline{Z} }_{c}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\overline{X} }_{cp}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\overline{Y} }_{cp}$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\overline{Z} }_{cp}$$\end{document} are fixed point values for the crosslinkers (red beads in Fig. 1B and C). \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{G}_{cp}^\mathrm{M}M}_{c}$$\end{document} in Eq. 5 are the modulus density functions used to determine the averaged internal elastic energy distribution of the mobile crosslinkers indexed by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left(c,p\right)$$\end{document} , respectively. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${G}_{cp}^\mathrm{M}$$\end{document} is defined as a unitless gain in this work (see APPENDIX). \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} are functions to provide two orthogonal DOFs per each point on the membrane surface. As illustrated in Fig. 1F, by introducing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} , the function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}={{K}_{1c}}^{2}+{{K}_{2c}}^{2}$$\end{document} ) along the membrane surface can be defined. Using two  \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} values as an orthogonal basis,  \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{K}_{1c}}^{2}+{{K}_{2c}}^{2}$$\end{document} defines their intrinsic value greater thanor equal to zero (i.e., squared length). The  \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} value at a certain point on the membrane may provide a combined quantity for the relative residence frequency and the elastic property of the mobile crosslinker. Here, it is not clear whether the elastic properties at different diffusion points are the same or different from each other. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} has the unit of Newton per fourth power of meter i.e., \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{N/{m}^{4}}$$\end{document} . Note that for all calculated solutions in this article, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} was almost equal to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} i.e., \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}\approx {K}_{2c}$$\end{document} . Finally, it is important to note that, in principle, the membrane area with a defined lipid environment needs to be large enough to fully cover the region with modulus density values greater than zero, because discontinuation of the lipid environment contradicts the second principle of this model.

The fourth term of Eq. 1 defines a constraint to preserve the total number of lipids \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_\mathrm{lipid}$$\end{document} by introducing a Lagrange multiplier \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\lambda }^\mathrm{lipid}$$\end{document} . Another Lagrange multiplier \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\lambda }_{c}^\mathrm{crs}$$\end{document} in the fifth term of Eq. 1 constrains the total value of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} i.e., \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\int \left({K}_{1c}+{K}_{2c}\right)dA$$\end{document} , to be an intrinsic constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} (unit: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{\mathrm{Newton}}$$\end{document} ). \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} determines the distribution and magnitude of the elastic modulus density, which reflect the confined diffusion area, residence frequency, and elastic property of the mobile crosslinker. Therefore, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} is the intrinsic constant that collectively characterizes the elastic property of the crosslinker, as well as the interfacial viscosity between the fluid membrane characterized by a certain lipid environment and the membrane inclusion connected to the elastic crosslinker. Function values for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} can be any real number. Infinitely many possibilities for the shape of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} function with values greater than or equal to zero on the membrane can be made by satisfying Eq. 6 and the condition \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\int \left({K}_{1c}+{K}_{2c}\right)dA={\varphi }_{c}$$\end{document} . Among the possibilities, however, the variational analysis gives us one solution profile satisfying the stationary condition of the total energy functional \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Pi }^\mathrm{tot}$$\end{document} .

It is worth noting that the expressions for the crosslinker are similar to quantum mechanical descriptions. Since the modulus density in this model and the probability density in quantum mechanics are conceptually similar, scaled versions of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} i.e., \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{{K}_{1c}}{\sqrt{{G}_\mathrm{scaling}}}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{{K}_{2c}}{\sqrt{{G}_\mathrm{scaling}}}$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{{M}_{c}}{{G}_\mathrm{scaling}}$$\end{document} where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${G}_\mathrm{scaling}$$\end{document} is the positive scaling factor with the unit of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{N/{m}^{2}}$$\end{document} , may correspond to the real part of the wavefunction \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Psi$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{Re}\left[\Psi \right]$$\end{document} ) and the imaginary part of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Psi$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{Im}\left[\Psi \right]$$\end{document} ), and the probability density i.e., the square of the modulus of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Psi$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\left|\Psi \right|}^{2}$$\end{document} ), respectively. This mathematical analogy may confirm our initial modeling objective to account for quantum mechanics for confined Brownian motion. According to the given notions, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left|\Psi \right|=\sqrt{\frac{{M}_{c}}{{G}_\mathrm{scaling}}}$$\end{document} can be held.

In this model paradigm, therefore, the wavefunction in quantum mechanics might be alternatively determined by asking the modulus density profile \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} and the corresponding \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1c}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2c}$$\end{document} functions. Here, the calculated \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} defines the combined quantity at each point for the relative residence frequency of a mobile object (the mobile end of a crosslinker in this work) and its energetic attraction property with respect to its conjugate (the fixed end of the crosslinker in this work) with given physical constraints (the membrane deformability in this work) and intrinsic physical constants. Overall, it is required to apply the theory to other quantum mechanical systems such as the particle-in-a-box problem^17,18^ and atomic orbitals to demonstrate its general applicability (see Discussion). It is also worth investigating how solutions of the Schrödinger equation and the presented variational theory of the modulus density are similar to or different from each other. Finally, we might conclude as follows regarding Brownian motion―there is no requirement for the energetic modulus to be a constant value, as in classical mechanics, according to this potentially quantum-compatible modulus density theory. Therefore, Brownian motion in classical mechanics and quantum mechanics are not fully compatible.

Using Eqs. 1‒6, a finite element method was developed, with a full description provided in the APPENDIX and SM. As a quasi-static problem, other time-dependent phenomena beyond confined diffusion—such as membrane viscoelasticity due to prescribed boundaries—are not considered in this model. An analysis of strain gradient generation around membrane-inserted proteins is presented in Fig. S2 (also see SM for discussion on this analysis). Finally, the method was used to investigate: (1) lipid sorting in the deformed membrane; (2) the interaction between the cell membrane and the adhesion proteins integrin and talin; (3) nanomechanical responses of cell membranes; and (4) the correlation of the variational modulus density theory to current theories in classical elasticity.

Lipid sorting and the formation of lipid nanodomains

In Fig. 2, a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${1.547\times 1.547 \; \mathrm{\mu{m}^{2}}}$$\end{document} planar membrane and nine regularly defined crosslinkers ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi =2\times {10}^{-5}\sqrt{\mathrm{N}}$$\end{document} ) were prepared to investigate lipid sorting due to membrane curvature^19^ and stretching (Fig. 2A and B). The distribution of the crosslinkers was determined by simplifying recent experimental data on CD44^4,20^. Eight lipid classes were assumed in the square membrane (see SM, Lipid sorting simulation). For each lipid class, the total number of lipids was the same, while the bending modulus \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\kappa}$$\end{document} was different. With systematic mechanical inputs, the eight different lipid classes were sorted, and lipid nanodomains were generated as shown in Fig. 2A and Movie S1 (see SM for details on lipid sorting and segmentation algorithms). The calculated membrane configuration is reminiscent of a mosaic made of different lipid nanodomains as described previously^21^. The lateral strain profile shown in Fig. 2B demonstrates that the lipid domain with a lower bending modulus is more stretched. According to Movie S1, regions with nonzero \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} values belong to the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\kappa}={10k}_\mathrm{b}T$$\end{document} membrane domain.Fig. 2. Lipid sorting and formation of lipid nanodomains. (A) Three-dimensional deformation of a planar membrane with the surface map for sorting of eight different lipid classes with varying bending moduli. Crosslinkers are not shown. Left: Angled view. Right: Planar view. (B) Lipid number strain α profile for the membrane in (A). Crosslinkers are shown here. (C‒J) The pattern of lipid nanodomains, shown separately for each lipid class. (K‒N) Expanded views of patterns in (D‒G). (O) Plots of the mean and standard deviation (SD) of the area of the nanodomains. Mean values are 20,183 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${10k}_\mathrm{b}T$$\end{document} , SD: 11 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} , n: 9), 19,785 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} (15.7 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_\mathrm{b}T$$\end{document} , SD: 39 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} , n: 9), 12,897 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} (21.4 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_\mathrm{b}T$$\end{document} , SD: 8476 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} , n: 13), 1786 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} (27.1 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_\mathrm{b}T$$\end{document} , SD: 1975 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} , n: 146), 1209 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} (32.9 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_\mathrm{b}T$$\end{document} , SD: 1117 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} , n: 222), 1060 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} (38.6 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_\mathrm{b}T$$\end{document} , SD: 1194 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} , n: 262), 1480 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} (44.3 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_\mathrm{b}T$$\end{document} , SD: 1886 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} , n: 194), and 8571 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${50k}_\mathrm{b}T$$\end{document} , SD: 9642 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{nm}}^{2}$$\end{document} , n: 35). Individual domain areas are indicated by colored dots. Blue: Multi-element domain. Red: Single-element domain. See Fig. S14 for the degrees of freedom of each element. See Movie S1 for all membrane-crosslinker configurations used in the analyses in this figure.

How the lipid sorting calculation can be correlated to measured compartmentalized diffusion of lipids^22–24^ was also investigated. To this end, the lipid domains with different bending moduli are separately plotted in Fig. 2C‒J. Figure 2C shows the disk-shaped domain, which is similar to the area formed by the confined diffusion of lipids (as an example see Fig. 4B in reference^23^). In Fig. 2F and G (and Fig. 2M and N), donut patterns for the lipid nanodomain are demonstrated. Here, multiple lipid nanodomains collectively form a single donut shape. This donut shape may explain measured hop diffusion trajectories in a closed form (see an example in Fig. 4C of reference^23^). The donut shape with a single lipid domain in Fig. 2D and E (and Fig. 2K and L) might be consistent with the orbital trajectory of the lipid observed from supported lipid bilayers^25^ and living cell membranes^26^. The other patterns in Fig. 2H‒J may explain more randomized hop diffusion trajectories reported in numerous previous experiments with an assumption that a certain lipid can only diffuse to energetically favored lipid nanodomains i.e., the domains with similar lipid parameter values. For a more direct comparison, areas of the nanodomains were calculated and plotted in Fig. 2O. These areas in the length scale of 32.6–142 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{nm}$$\end{document} were similar to measured values for the hop diffusion compartment^27^.

Mechanical responses of talin and integrins and their interaction with cell membranes in the formation of nascent adhesions

Integrins are focal adhesion receptors that interact with the cytoplasmic adaptor talin^1–3^. To investigate the formation of nascent adhesions mediated by integrins, a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2.32\times 2.32 \; \mathrm{\mu {m}^{2}}$$\end{document} membrane and two crosslinker groups with a total of three crosslinkers around the center region of the membrane were prepared as a basis complex for adhesion molecular machineries. Here, the crosslinker with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=1$$\end{document} represents a talin monomer with an inactivated integrin heterodimer (see inset diagram in Fig. 3, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{1}=9\times {10}^{-6} \sqrt{\mathrm{N}}$$\end{document} ). The crosslinker with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=2$$\end{document} represents a complex of a single integrin heterodimer and its extracellular ligand ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{2}=2\times {10}^{-5} \sqrt{\mathrm{N}}$$\end{document} ). The crosslinker with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left(c,p\right)=\left(\mathrm{2,1}\right)$$\end{document} represents another talin monomer connected to the integrin-ligand complex ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${G}_{21}^\mathrm{M}=0.03$$\end{document} ). For this additional crosslinker with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left(c,p\right)=\left(\mathrm{2,1}\right)$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\overline{X} }_{cp}={\overline{X} }_{c=1}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\overline{Y} }_{cp}={\overline{Y} }_{c=1}$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\overline{Z} }_{cp}={\overline{Z} }_{c=1}$$\end{document} were initially defined to simulate a talin dimer whose two monomers share a fixed point on actin cytoskeletons.Fig. 3. Calculations of the interaction among integrins, talin, and membranes in the formation of nascent adhesions. (A) Force vs. extension responses with respect to inputs in Fig. S3. The extension of the crosslinker was defined from its fixed point (red bead) to the maximum \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} point on the membrane (black bead). Expanded views of indicated regions (arrows) are provided in Fig. S5. (B, C) Half-max \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} area vs. extension responses for the crosslinkers. (D‒F) Deformed shapes of the membrane-crosslinker complex. The data point index is indicated at the top. Surface maps indicate the membrane strain \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} and the sum of normalized elastic modulus densities at the top and bottom, respectively. See Fig. S4 for the membrane-crosslinker complex from a different viewpoint. See Fig. S6 for the surface map of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=1$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=2$$\end{document} . See Fig. S15 for the degrees of freedom. See Movie S2 for all membrane-crosslinker configurations used in the analyses in this figure.

Investigations of the integrin-membrane-talin complex consisted of five parts, each with different inputs (Fig. S3 and Fig. 3). The first simulation (sim1, data point index 1‒28), shown with black traces, was performed by pulling the crosslinker with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=2$$\end{document} . This simulation corresponds to the integrin-ligand interaction during the formation of initial adhesion. The lateral size of membrane tenting was microscopic (Fig. 3D and Fig. S4B). Forces applied to the crosslinker with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=2$$\end{document} are plotted in Fig. 3A and Fig. S5A (black circles). Here, the final force, approximately 6.5 pN, was within one significant force population in 5.5–7 pN estimated for the integrin-ligand interaction by using fluorescence resonance energy transfer (FRET) sensors^28^. It is important to note that agreements between FRET force measurements and model calculations here and in the remainder of this paper cannot be used to support the validity of the proposed theory, since the FRET force sensor is based on the conventional mechanics approach (see Discussion). However, the agreements may indicate that molecular systems with fluid membranes exist where both the conventional approach and the new theory work well in estimating applied forces (see Discussion). By using the given \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{1}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{2}$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${G}_{21}^\mathrm{M}$$\end{document} values, forces applied on the talin crosslinkers were 2.3 pN (for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=1$$\end{document} ) and 2.2 pN (for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left(c,p\right)=\left(\mathrm{2,1}\right)$$\end{document} ) at the end of sim1.

The second simulation (sim2, data point index 29‒32), shown with green traces in Fig. S3 and Fig. 3A‒C, was performed by pulling the talin dimer ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=1$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left(c,p\right)=\left(\mathrm{2,1}\right)$$\end{document} ) in the cytoplasmic direction. This pulling of talin dimer during the initial phase of the adhesion may happen with the contraction of actin networks before the formation of stress fibers^29^. Forces applied to talin, ranging from 2.39 to 2.94 pN, were similar to the minimum estimated talin force of about 2.65 pN in a previous FRET experiment^30^. Forces applied to the activated integrin crosslinker remained within the measured FRET force population in 5.5–7 pN (Fig. 3A and Fig. S5A, green circles)^28^.

In the third simulation (sim3, red traces in Fig. S3 and Fig. 3, data point index 33–52), the fixed point for the talin dimer was randomly displaced within an area radially defined by 100 nm from the center. In this part, the average separation between two membrane-targeting ends of the talin dimer was 162.38 nm when projected onto the X–Y plane. This value is consistent with a previous measurement using super-resolution microscopy^31^. The separation was 72.9 nm when the point was randomly displaced within an area defined by 20 nm from the center. The result suggests that the motility of the underlying actin networks is critical for the talin dimer separation. Finally, together with the surface map plots for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} (bottom panel of Fig. 3D, E and Fig. S6A, B), the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\int {M}_{c}dA$$\end{document} plot in Fig. S7 demonstrates that the model finds the solution i.e., the deformed configuration of the membrane-crosslinker complex, by varying both the spatial distribution and the integrated total amount of the elastic modulus density \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} with given fixed \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${G}_{cp}^\mathrm{M}$$\end{document} values.

The fourth simulation (sim4, data point index 53‒54), shown with cyan traces, was performed by increasing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{1}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${G}_{21}^\mathrm{M}$$\end{document} . Since \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} is the intrinsic constant, any change in its value may correspond to structural modification of the crosslinker, including its significant refolding, unfolding, and overstretching when the membrane interface remains unchanged. For talin, this can also occur with the change of force-acting sites due to the binding of vinculin along the talin rod^32–34^. The change in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${G}_{cp}^\mathrm{M}$$\end{document} modifies the relative magnitude of the elastic modulus density of the corresponding crosslinker within its group. By increasing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{1}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${G}_{21}^\mathrm{M}$$\end{document} , as shown with cyan traces in Fig. S3D and E to mimic the effect of talin fastening with vinculin, the initially tented membrane was contracted under the X–Y plane (Fig. 3F and Fig. S4D).

In the fifth simulation (sim5, data point index 55–74) shown with blue traces, the shared point for the talin dimer was separated, and each monomer was randomly displaced, as shown in Fig. S3A–C. The integrin forces ranged from 11.2 pN to 11.9 pN. The calculations were consistent with previous FRET experimental results—these experiments indicated that vinculin is required for the generation of integrin forces greater than 7 pN^28^. With significant membrane deformation here, the forces on the two talin rods, ranging from 5.5 to 6.2 pN, were also within measured values using the FRET sensors^30,35^. In sim2‒5, the extension values of the two talin monomers were within the experimentally identified range from living cells^36^.

According to a previous study, integrins on living cell surfaces can undergo confined diffusion on the length scale of a hundred or hundreds of nanometers^37^. In the experiment, the radii of the confinement ranged from about 119 nm (inside of the focal adhesion, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4.4488\times {10}^{4} \; \mathrm{{nm}^{2}}$$\end{document} ) to 236 nm (outside of the focal adhesion, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.7497\times {10}^{5} \; \mathrm{{nm}^{2}}$$\end{document} ) (see Fig. 1i in the reference paper^37^). Furthermore, a recent experiment directly visualized that the initially freely mobile integrin was confined in a similar length scale with the talin colocalization^38^. This confined diffusion area was about \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.102\times {10}^{5} \; \mathrm{{nm}^{2}}$$\end{document} (Fig. 2A right in the reference paper^38^; see SM for the area calculation)^38^. These observations support the crosslinker mechanism for the confined integrin diffusion. Remarkably, areas defined by half of the maximum \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} value in Fig. 3B and C were consistent with the confined integrin diffusion area observed in the experiments^37,38^ (see SM for the calculation of the half-max \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} area).

Stretching and sliding of crosslinkers with membrane deformation are visualized in Movie S2 for all data in Fig. 3. In Figs. S8–10 and Movie S3, molecular machineries with four adhesion sites are investigated by using the same membrane and parameter values. The results were largely consistent with those from Fig. 3. Overall, the presented analyses demonstrate the remarkable ability of the model to predict biological responses in the formation of nascent adhesion at the single-molecule level.

The effect of the number and distribution of crosslinkers, as well as the composition of lipids, on nanomechanical responses of membrane-crosslinker complexes

To reproduce recently identified nanomechanical responses of cell membranes^11^, five crosslinkers were defined in a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.944\times 1.944 \; \mathrm{\mu {m}^{2}}$$\end{document} membrane where the distance from one at the center to the other four crosslinkers is about 488.8 nm (Fig. 4A, Fig. S11A, and Movie S4). The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} value for the crosslinker located at the center ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=3$$\end{document} ) was \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2.6\times {10}^{-5}\sqrt{\mathrm{N}}$$\end{document} . \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} for the other crosslinkers ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{c}=1, 2, 4, 5$$\end{document} ) that represent the membrane-cytoskeleton linkage was \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$6.5\times {10}^{-6}\sqrt{\mathrm{N}}$$\end{document} . This value was selected to investigate whether \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document} smaller than that of the full-length talin monomer (with the inactivated integrin) can generate the nanomechanical responses. To calculate the force vs. displacement response, the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\overline{Z} }_{c}$$\end{document} value for the center crosslinker ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=3$$\end{document} ) was increased. The force vs. displacement curve shape in Fig. 4C (blue circles) showed good agreement with that of a measured nanomechanical response of cell membranes^11^.Fig. 4. The effect of the number and distribution of crosslinkers and lipid composition on nanomechanical responses of membrane-crosslinker complexes. (A) The deformed configuration with the surface strain map for a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.944\times 1.944 \; \mathrm{\mu {m}^{2}}$$\end{document} membrane ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\kappa}=20{k}_\mathrm{b}T$$\end{document} ) and five crosslinkers when the center crosslinker ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=3$$\end{document} ) was pulled with a 249 nm displacement (left). An expanded view of the membrane-crosslinker complex with the surface map for the sum of normalized \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} values (right). (B) A deformed configuration with the surface strain map for the membrane-crosslinker complex in (A) by considering 75% \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$50{k}_\mathrm{b}T$$\end{document} lipids and 25% \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$20{k}_\mathrm{b}T$$\end{document} lipids. The corresponding lipid sorting map is shown on the right. (C) The force vs. displacement of the red bead at the center ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=3$$\end{document} ) in (A) is shown with blue circles, and compared to experimental data^11^. The force vs. displacement for the center red bead in (B) is shown with red diamonds. (D) The deformed configuration with the surface strain map for the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.944\times 1.944 \; \mathrm{\mu {m}^{2}}$$\end{document} membrane ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\kappa}=20{k}_\mathrm{b}T$$\end{document} ) and twenty-one crosslinkers when the center crosslinker ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=11$$\end{document} ) was pulled with a 281 nm displacement (left). An expanded view of the membrane-crosslinker complex with the surface map for the sum of normalized \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} values (right). (E) A deformed configuration with the surface strain map for the membrane-crosslinker complex in (D) by considering 75% \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$50{k}_\mathrm{b}T$$\end{document} lipids and 25% \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$20{k}_\mathrm{b}T$$\end{document} lipids. The corresponding lipid sorting map is shown on the right. (F) The force vs. displacement of the red bead at the center ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c=11$$\end{document} ) in (D) is shown with blue circles and compared to experimental data^11^. The force vs. displacement for the center red bead in (E) is shown with red diamonds. See Figs. S16 and S17 for the degrees of freedom. Movies S4‒S7 visualize all membrane-crosslinker complexes used in the analyses in this figure.

To reproduce another type of nanomechanical response, twenty-one crosslinkers were compactly inserted in a region spanning 305.5 nm from the center of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.944\times 1.944 \; \mathrm{{\mu {m}}^{2}}$$\end{document} membrane (see Fig. 4D, Fig. S11B, and Movie S6). The shortest distance from the center to a crosslinker was 61.1 nm (Fig. S11B). \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}=1.6\times {10}^{-5}\sqrt{\mathrm{N}}$$\end{document} for the centered crosslinker ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{c}=11$$\end{document} ), while \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}=6.5\times {10}^{-6}\sqrt{\mathrm{N}}$$\end{document} for the others were used. Increasing the number of crosslinkers within a reduced area resulted in a more localized deformation shape of the membrane. The crescent shape of the area for the elastic modulus density \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${M}_{c}$$\end{document} around the sharp membrane curvature ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{c}=7, 8, 14, 15$$\end{document} ) demonstrates the mechanistic and hydrostatic interaction between the fluid membrane and the crosslinkers connected to mobile membrane inclusions (Fig. S12). With the given twenty-one crosslinkers, the force vs. displacement response showed zero stiffness at an intermediate displacement (Fig. 4F, blue circles). A comparison between the calculated zero stiffness response and experimental data from a previous work also showed good agreement (Fig. 4F, blue circles)^11^. According to Fig. S13 and Movie S8, reducing the distance to the membrane-cytoskeleton linker from the center, without a sufficient increase in the number of crosslinkers, can result in instabilities that break the symmetry of the deformation of the membrane-crosslinker complex.

The effect of lipid composition on the generation of the nanomechanical response was also investigated. The force vs. displacement responses with five and twenty-one crosslinkers were reexamined using two different lipids ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\kappa}=50{k}_\mathrm{b}T$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\kappa}=20{k}_\mathrm{b}T$$\end{document} ) in Fig. 4C and F, respectively (red diamonds). For simplicity, the same \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{c}$$\end{document} values in Fig. 4A and D were used for the five and twenty-one crosslinkers, respectively, in the calculations with the two lipid classes, even though the lipid environment here is not exactly the same as that of Fig. 4A and D. With mechanical pulling of the center protein, two lipid classes were sorted and patterns were generated in the membrane in both cases, as shown in Movies S5 and S7. The central region with high membrane curvatures was filled with the lipid of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\kappa}=20{k}_\mathrm{b}T$$\end{document} for both cases. Unlike the sorting investigation in Fig. 2, the collapse of the boundary between two lipid domains was identified as the lateral membrane stretching became significant with larger displacements (Fig. 4B and E; Movies S5 and S7). Notably, for both cases, the response in the high displacement region shifted when the lipid of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\kappa}=50{k}_\mathrm{b}T$$\end{document} was considered alongside the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$20{k}_\mathrm{b}T$$\end{document} lipid (Fig. 4C and F). This shifting effect is consistent with the effect of reducing the size of lipid reservoirs in mechanical stretching of membranes, as analyzed in the previous work^11^. According to the presented analysis, therefore, the reduced size of lipid reservoirs in membrane stretching is a result of the presence of different lipid classes with higher bending moduli and sorting of the lipids in the deformed membrane.

The correlation of the variational modulus density theory to current theories in classical elasticity

How the two principles in the modulus density theory can be applied to the crosslinker with two fixed ends was investigated by using previous magnetic tweezer data on a single monomeric talin^31^. In Fig. 5, the experimental tensile force vs. extension data, obtained from single full-length talin monomers fixed at both ends (one at the magnetic bead and the other at the substrate), and interacting with vinculin molecules, were compared to a calculated curve using the following Eq. 7:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{f}}_\mathrm{fixed-talin}=\frac{{{\varphi }_\mathrm{fixed-talin}}^{2}{{\Delta {L}}_{\mathrm{talin}}}^{2}}{2{{A}}_{\mathrm{talin}}}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Delta L}_\mathrm{talin}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${A}_\mathrm{talin}$$\end{document} are the talin extension and the coarse-grained cross-section area for talin, respectively. The two principles, described in the Introduction section with respect to the mobile crosslinker, are applied in deriving Eq. 7 as follows. First, the force for talin in a fixed configuration can be calculated by finding the modulus density value used to calculate the internal elastic energy of talin. Here, the statistical description in the first principle for the mobile crosslinker is conceptually simplified for the fixed crosslinker. Second, by assuming that the response of the fixed talin is also governed by its intrinsic physical constant termed \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_\mathrm{fixed-talin}$$\end{document} , all possible values for the modulus density are derived from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_\mathrm{fixed-talin}$$\end{document} , where the single modulus density value that generates the minimum energy is selected to calculate the talin force. In mathematics, which is employed to describe natural phenomena, decomposition (or addition) is the most fundamental operation. Accordingly, as the best way,  \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_\mathrm{fixed-talin}$$\end{document} is decomposed to generate multiple variations of the modulus density value. As similarly done for the mobile crosslinker, using two variables to decompose \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_\mathrm{fixed-talin}$$\end{document}  might be enough to maximize the size of the set of possible values (see APPENDIX for the derivation of Eq. 7). The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_\mathrm{fixed-talin}$$\end{document} value was \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\times {10}^{-7 } \sqrt{\mathrm{N}}$$\end{document} when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${A}_{\mathrm{talin}}=10.2 \; \mathrm{{nm}^{2}}$$\end{document} in an extension range of about 0–800 nm (Fig. 5). It is important to note that the nature of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_\mathrm{fixed-talin}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document} for the mobile crosslinker are different from each other, as the concept of the interface viscosity is absent in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_\mathrm{fixed-talin}$$\end{document} .Fig. 5. Application of the modulus density theory to the crosslinker with two fixed ends. Calculations using Eq. 7 with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_\mathrm{fixed-talin}=1\times {10}^{-7}\sqrt{\mathrm{N}}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${A}_{\mathrm{talin}}=10.2 \; \mathrm{{nm}^{2}}$$\end{document} (circles). Both ends of the crosslinker are assumed fixed. Calculations are compared with magnetic tweezer measurements for single full-length talin monomers interacting with vinculin molecules in solution^31^. The calculations and measurements show reasonable agreement.

Classical mechanics laws such as Hooke’s law and nonlinear elasticity state that the force on a one-dimensional elastic spring is proportional to its extension (or a function of the extension), through a modulus such as Hooke’s spring constant or Young’s modulus. The formalism of the variational modulus density theory for the two-fixed-ends crosslinker is aligned with the notion of those conventional theories. According to the derivation of Eq. 7, the intrinsic constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_\mathrm{fixed-talin}$$\end{document} yields the constant elastic modulus value under the minimum energy assumption (see APPENDIX for the derivation of Eq. 7). Overall, the analysis in Fig. 5 suggests that the principles of the modulus density theory are identically applicable to the elastic linker with two fixed ends.

In Fig. 6, a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.944\times 1.944 \; \mathrm{{\mu m}^{2}}$$\end{document} membrane with a single crosslinker at its center was used to further investigate the characteristics of the crosslinker interacting with the fluid membrane. With four different \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document} values, the responses of the crosslinker varied significantly (Fig. 6A). In this model framework, a single \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document} value collectively defines both the elastic property of the crosslinker and the viscosity of the membrane inclusion connected to it. Therefore, such non-negligible variations may be generated among identical crosslinkers if their inclusion properties are significantly different. This suggests that estimating the crosslinker force (or energy) from its extension using conventional elasticity might not always be valid when the crosslinker interacts with the fluid membrane.Fig. 6A combined analysis for the force transfer at the membrane-crosslinker interface using the modulus density theory. (A) The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.944\times 1.944 \; \mathrm{{\mu m}^{2}}$$\end{document} membrane was stimulated by using four different complexes of crosslinkers and membrane inclusions ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi = 2.6\times {10}^{-5}\sqrt{\mathrm{N}}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi = 1.6\times {10}^{-5}\sqrt{\mathrm{N}}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi = 0.6\times {10}^{-5}\sqrt{\mathrm{N}}$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi = 0.2\times {10}^{-5}\sqrt{\mathrm{N}}$$\end{document} ). Forces applied to the crosslinkers were plotted against the extension of the crosslinkers. The four responses vary significantly. (B) The forces on the crosslinkers in (A) were plotted against the displacement of the membrane. Here, the response of pulling the membrane without the center crosslinker was additionally plotted (black empty squares). The responses overlap with small variations. See Fig. S18 for the degrees of freedom.

According to Fig. 6B, however, the crosslinker force vs. membrane displacement responses overlapped for all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document} values used. They were also consistent with the membrane force response directly obtained by taking the derivative of the membrane energy without considering the pulling crosslinker (black empty squares in Fig. 6B). While the validity of Hooke’s law for the crosslinker under confined Brownian diffusion is elusive, the results in Fig. 6B suggest that Newton’s third law is still valid in this new paradigm, specifically at the membrane-crosslinker interface. Therefore, the membrane deformation can be used to approximate the force on the crosslinker, regardless of its elastic property and membrane-inserted interface, as similarly conducted in previous works^10,11,39,40^. Small variations among the five curves in Fig. 6B might be generated numerically or due to differences in modulus density profiles that can minimally influence the curvature of membranes at their tips. Overall, the results in Figs. 5 and 6 suggest that the presented variational modulus density theory provides a more general framework for explaining one-dimensional elastic linkers not only interacting with the fluid membrane but also in standalone configurations. Finally, the results emphasize the importance of characterizing the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document} value for various complexes of crosslinkers and their membrane-inserted molecules.

Discussion

There is a viewpoint that considers Brownian diffusion as a quantum phenomenon^6–8^, while membrane deformation has been explained by classical continuum mechanics^9–11^. To explain confined Brownian diffusion in the membrane, the main modeling objective in this work was to develop a new theory in which quantum mechanical descriptions can be meaningfully reflected in the classical mechanics framework. To this end, the variational modulus density theory, incorporating two novel principles, was developed. The principles can be summarized as follows: (1) the shared value for a certain physical characteristic across multiple bodies is equal to the value of the characteristic in their average configuration; and (2) all ensemble averages of confined physical motions of an object might be mathematically predictable with an assumption that the motions are governed by an intrinsic physical constant. This new modeling approach, for the confined Brownian motion and mechanical deformation of crosslinkers interacting with deforming cell membranes, establishes a mathematical analogy between the proposed theory and the quantum mechanical wavefunction.

For a crosslinker interacting with a solid membrane, a one-dimensional force based on conventional elasticity can be used to represent the crosslinker force at a specific point on the membrane (corresponding to Fig. 1A, left). However, this approach may not be directly applicable to a fluid membrane, which differs fundamentally from a solid one. As a modeling strategy in this study, the crosslinker force was calculated by averaging over multiple configurations of the mobile crosslinker (Fig. 1A, right). In this new paradigm, it is no longer necessary to assume a constant elastic modulus as in the conventional theory. Instead, the force calculation relies on the averaged configuration of the mobile crosslinker, rather than individual configurations. Assuming that the crosslinker force remains constant during confined diffusion, it is reasonable to equate the force applied by each individual configuration with that of the averaged configuration. To represent this averaged behavior, the concept of elastic modulus density was introduced.

A mathematical framework to maximize the variability of the shape of the modulus density function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M$$\end{document} made from the intrinsic constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document} is also proposed. In this framework, it is necessary to use the integral summation to correlate the constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document} with the two-dimensional functions \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2}$$\end{document} . At each point on the membrane surface, two function values from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2}$$\end{document} define the modulus density value. By doing so, the “plus” and “minus” operations (i.e., two fundamental operators in summation) can be independently considered at each point in the integration process leading to the constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document} . In other words, employing the two degrees of freedom (DOFs) allows the most dynamic decomposition of the single constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document} into real numbers distributed across a continuous area. Consequently, this maximizes the variability of the resulting modulus density profile derived from the constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document} . More precisely, all physically possible confined diffusion configurations for a certain membrane inclusion, linked to a certain crosslinker and interacting with a certain lipid environment, might be generated using this mathematical framework, with an assumption that the physical motion of the inclusion-crosslinker complex is governed by its intrinsic constant. The functions \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{1}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${K}_{2}$$\end{document} may correspond to the real and imaginary parts of the wavefunction in quantum mechanics. How the mathematical identification for the “positive” and “negative” components of the modulus density value at each point can be interpreted in a more physical manner is an additional aspect to be explored in the future.

The consistency of the calculated results with experimental data suggests that the mobile and force-bearing crosslinker may be better explained using this new approach. Most notably, the area defined by half of the maximum \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M$$\end{document} value for the force-bearing talin was similar to the measured area for the integrin under confined diffusion^37,38^. Additionally, single-molecule magnetic tweezer responses for talin were reproduced^31^. The separation distance between two membrane-targeting ends of the talin dimer was consistent with a previous experiment^31^. Also, the model directly reproduced measured nonlinear nanomechanical responses of living cell membranes^11^. These agreements collectively support the validity of the model.

The possibility of generating different force values for the same crosslinker with the same extension but with different membrane interfaces may question the validity of FRET force sensors interacting with the fluid membrane, because the technology parameterizes the extension of sensor molecules to estimate the force, as in the conventional mechanics theory^28,30,35^. In this work, calculated forces for talin were similar to FRET measurements from living cells^30,35^. Calculated forces for integrin-ligand interactions at adhesions were also consistent with experimental FRET data^28^. Assuming that the new theory is reasonable, these agreements suggest that the conventional mechanics approach and the corresponding FRET sensor might be effective in determining forces on various membrane-associated proteins such as talin and the integrin. Nevertheless, it is still unclear whether the conventional approach can generate the same confined diffusion area. Overall, more direct experimental investigations into the validity of the FRET force sensor with respect to many other membrane protein systems and lipid environments might be required.

One particular focus in this work was the sorting of lipids and the generation of lipid nanodomains due to membrane curvatures and stretching. The size and shape of the nanodomains identified in this work were similar to those of membrane compartments identified by tracking the diffusion of lipids in living cell membranes. The picket-fence model suggests that anchored transmembrane proteins serve as physical barriers against the diffusion of lipids^22,41^. However, recently identified CD44 proteins do not fully fence a closed area in the membrane while they are responsible for membrane compartmentalization^4,20^. The result may ask an alternative explanation in which the compartment can be made with few picket proteins. The presented analyses suggest that lipids can diffuse only into energetically favored nanodomains i.e., the domains with the same lipid parameter values, and that might be responsible for the generation of the observed compartmentalized diffusion. In this alternative paradigm, CD44 simply serves as the membrane-cytoskeleton crosslinker to generate membrane curvatures and stretching. The presented interpretation can also explain the hop diffusion observed in the supported lipid bilayer where picket-fence proteins are absent^25,42,43^. The lipid sorting study in this work further explains how the presence of different lipid classes reduces the size of lipid reservoirs in mechanical deformation of cell membranes^11,44–46^. Overall, the work demonstrated that the effect of molecular composition in lipid membranes can be analyzed in detail within the framework of continuum mechanics, as actively investigated recently for various problems including vesicle budding with lipid demixing^47^ and the interaction between amphipathic peptides and lipid membranes^48,49^.

As mentioned earlier, the demonstration of the general applicability of the modulus density theory to other physical systems might be important. In this regard, investigating atomic orbitals using the theory would be intriguing. The electron is mobile in three-dimensional space and is energetically attracted to the nucleus. Additionally, multiple electrons can physically constrain each other’s motion due to mutual repulsion. Moreover, the kinetic energy of electrons can push them away from the nucleus. Therefore, the physical and mathematical characteristics of the atomic orbital may share similarities with those of the membrane-crosslinker complex examined in this work (though the kinetic energy is not considered in this work; see SM for an estimate for the kinetic energy of confined diffusion). The electron cloud in the conventional model could be viewed as a volume characterized by the modulus density value for the electron-nucleus interaction in this model paradigm. Likewise, the presented theory may be applicable to modeling various flow phenomena, such as those involving liquids, information, or currency. Additionally, conventional theories based on spring-like, one-dimensional representations for the interaction between two fluid particles—such as molecular fluid dynamic models—may need to be reformulated based on this work. Finally, investigating planetary trajectories and artificial intelligence through the lens of this theory and its underlying principles could open up fascinating extensions of the presented research.

In summary, the area of the confined Brownian motion for force-bearing crosslinkers was explored by introducing a continuum theory that is potentially compatible with the quantum wavefunction. The variational modulus density theory for crosslinkers was developed and integrated into the continuum model for cell membranes. A novel finite element method for the combined model was developed, which not only reproduced numerous experimental data but also provided valuable predictions regarding the lateral organization of the cell membrane and the membrane-crosslinker interaction. The model also suggested potential answers for the physical and mathematical identity of the quantum wavefunction. The theoretical model presented in this work is expected to provide enhanced capabilities for the realistic simulation of cell membranes and membrane crosslinkers.

Footnote

All figure data, as well as the initial version of the manuscript and the Supplementary Material, were generated at INTEGRITY Co.,Ltd. (9, Gangnamseo-ro, Giheung-gu, Yongin-si, Gyeonggi-do, Republic of Korea, 16977), which is no longer active.

Electronic supplementary material

Below is the link to the electronic supplementary material.Supplementary Information 1.Supplementary Video 1.Supplementary Video 2.Supplementary Video 3.Supplementary Video 4.Supplementary Video 5.Supplementary Video 6.Supplementary Video 7.Supplementary Video 8.