The Effect of Thermal Fluctuation on the Receptor-Mediated Adhesion of a Cell Membrane to an Elastic Substrate

Abstract

Mechanics of the bilayer membrane play an important role in many biological and bioengineering problems such as cell–substrate and cell–nanomaterial interactions. In this work, we study the effect of thermal fluctuation and the substrate elasticity on the cell membrane–substrate adhesion. We model the adhesion of a fluctuating membrane on an elastic substrate as a two-step reaction comprised of the out-of-plane membrane fluctuation and the receptor–ligand binding. The equilibrium closed bond ratio as a function of substrate rigidity was computed by developing a coupled Fourier space Brownian dynamics and Monte Carlo method. The simulation results show that there exists a crossover value of the substrate rigidity at which the closed bond ratio is maximal.

1. Introduction

Mechanics of the bilayer membrane play an important role in many biological and bioengineering problems such as membrane remodeling during cell mechanoadaptation [1], cell–substrate and cell–nanomaterial interactions [2], and supported bilayer membranes [3,4,5,6]. In this work, we focus on the effect of thermal fluctuation and the substrate elasticity on the cell membrane–substrate adhesion. Cell–substrate adhesion has attracted considerable interest in past decades due to its importance in mechanobiology as well as in biotechnology where living cells interact with a variety of substrates [7,8,9]. Molecular and cellular mechanisms by which cells “sense” and respond to the mechanical properties of the extracellular matrices can, in general, fall into two categories. In the first category, mechanosensing is actively executed through molecular mechanoenzymatics by which mechanical forces is converted into biochemical signals [10]. The other category is through “passive” physical laws that are derived from the kinetics or energetics of the processes in cell–matrix interactions [11]. Due to the complexity of the biological system, it is very likely that principles from both of these two categories contribute to the cell–matrix interactions. When focusing on the cell–substrate adhesion, it is generally considered that different phases can be distinguished [12]. In the early phase (e.g., the first a few minutes of cell–substrate adhesion), cells come in contact with the substrate to form nascent adhesion patch mainly via receptor–ligand binding, and in the later phases, much more complicated processes are involved such as receptor clustering, cytoskeleton remodeling, and focal adhesion formation [13].

While considerable mechanics modeling efforts have been devoted to the later phase (see a review [11] for more details), the present work is mainly focused on the early phase. The early phase, despite being relatively simpler than the later phase, has also attracted great interest. For example, dedicated vesicle adhesion experiments [14,15] where binding molecules are present in the vesicle membrane and on the substrate have been carried out to mimic cell adhesion. In these experiments, out-of-plane membrane fluctuation was used to detect receptor–ligand binding. In particular, it was found [15] that equilibrium contact zones by ligand–receptor binding were formed and suppressing membrane fluctuation leads to increase of adhesion strength. Furthermore, it has been postulated [16] that cells utilize membrane fluctuation to probe surface properties in the first minutes of cell–substrate contact.

Membrane fluctuation was thought to play an essential role in the early stage of cell adhesion, which is also the central theme of the present work. Besides the membrane–substrate adhesion, membrane fluctuation has also been thought to play important roles in many other cellular processes, such as being utilized to avoid non-specific adhesion [17,18], creating transient gap for actin monomer intercalation during cell migration [19], inducing long-range interactions between membrane proteins or domains [20], and affecting protein mobility on cell membranes [21,22]. In this paper, we aim at developing a quantitative simulation model for the adhesion of the fluctuating membrane and to ascertain whether the membrane fluctuation can be “utilize” to probe the substrate rigidity and receptor density.

2. The Model

A pioneering work on theoretical modeling of cell adhesion mediated by reversible receptor–ligand bonds was presented by Bell [23], in which the cell adhesion was modeled as a two-step reaction. In the first step, binding proteins diffuse in membranes to form receptor–ligand encounters, and in the second step, receptor–ligand encounters react to form closed bonds. In our case of membrane–substrate, we assume receptors in the membrane are dilute and ligand density on the substrate is saturated, so the effect of the in-plane receptor diffusion on forming receptor–ligand encounters is neglected and receptors are assumed to be immobile and uniformly distributed in the membrane (as illustrated in Figure 1). On the other hand, the cell membrane undergoes out-of-plane undulation as a form of Brownian motion, which brings receptors close to the substrate to form receptor–ligand encounters or inversely pull closed bonds apart. Qualitatively, we can take the same idea as Bell to model the adhesion of a fluctuating membrane as a two-step reaction by simply replacing the in-plane binder diffusion of Bell’s model with the out-of-plane membrane fluctuation,

$$\mathrm{R}+\mathrm{L}\rightleftarrows \mathrm{R}--\mathrm{L}\underset{k_{off}}{\overset{k_{on}^0}{\rightleftharpoons }}\mathrm{RL}$$

(1)

where the first step concerns the formation of encounter complex $\mathrm{R}--\mathrm{L}$ via membrane fluctuation and the second step is the formation of receptor–ligand bonds with single-molecule reaction rate $k_{on}^0$. The forward reaction of the second step can only occur when the ligand–receptor distance is within the encounter distance R_RL. The unbinding rate constant k_off takes Bell’s formula [23] as $k_{off}=k_{off}^0\exp \left(f{x}_b/k_{B}T\right)$, where f is the force acting on a ligand–receptor bond by the membrane pulling, x_b describes how strongly the reaction rate change with force, $k_{off}^0$ is the rate at zero force, k_B is Boltzmann constant, and T is the temperature. As showed by Bell [23], for the case of classical particle diffusion, rates of the diffusion step are simple functions of particle diffusivities and encounter distance. In contrast, the two-step reaction in this model is complicated by the membrane fluctuation and we resort to computer simulations. Throughout this paper, we use the closed bond ratio ϕ, which is simply defined as the number of closed ligand–receptor bonds divided by the total number of receptors in the membrane, to characterize the adhesion strength. The objective of this paper is to carry out numerical simulations of the chemical reaction in Equation (1) to calculate the closed bond ratio ϕ, and to study how it depends on substrate rigidity and receptor density by involving the membrane fluctuation step. It shall be mentioned that previous cell adhesion models [23,24,25,26,27] mainly focused on the rupture strength of a cluster of bonds under an external pulling force, while in our present work, we are concerned with the equilibrium closed bond ratio in the case of no external pulling force.

The fluctuation or Brownian motion of the membrane in potential fields will be simulated using a Fourier space Brownian dynamics (FSBD) scheme developed by Lin and Brown [22,28]. Single ligand–receptor binding-unbinding kinetics will be modeled via Monte Carlo (MC) method and hereafter we refer to this coupled scheme as the FSBD-MC (Fourier space Brownian dynamics- Monte Carlo) simulation. Membrane fluctuation mediated cooperative behavior [29] among different ligand–receptor pairs will be naturally taken into account in the FSBD-MC simulations. On the other hand, particle-based membrane simulation models [30,31,32] are more suitable and straightforward if one needs to simulate receptor diffusion and membrane fluctuation simultaneously. However, in the present work, since receptor diffusion is neglected and only single-valued out-of-plane deformation of a planar membrane is examined, adopting such a Fourier mode-based continuum model can be efficient for the large membranes. The overdamped Langevin equation in Fourier space, with the amplitudes of the undulation modes as the degrees of freedom, reads [28]

$$\frac{\partial h_q\left(t\right)}{\partial t}={\Lambda }_q\left[F_q\left(t\right)+{\zeta }_q\left(t\right)\right]$$

(2)

where $h_q={\displaystyle {\int }_{A}h\left(x\right)\exp \left(-iq\cdot x\right)dx}$ and inversely $h\left(x\right)=L^{-2}{\displaystyle \sum _q{h}_q\exp \left(iq\cdot x\right)}$. As illustrated in Figure 1, $x=\left(x,y\right)$, and $h\left(x\right)$ is the membrane height, A = L², L is the membrane side length, ${\Lambda }_q=1/\left(4\eta q\right)$, η is the viscosity of the surrounding fluid, and $q=\left|q\right|$, $q=\left(2\pi \alpha /L,2\pi \beta /L\right)$, α and β are integers from 0 to α_max. Note that, in the numerical simulations, two wavelength cutoffs are defined: λ_min = L/α_max and λ_max = L. The stochastic force ${\zeta }_q\left(t\right)$ has a Gaussian distribution with $\langle \left|{\zeta }_q\left(t\right)\right|\rangle =0$ and $\langle \left|{\zeta }_q\left(t\right){\zeta }_{q^{\prime }}\left(t^{\prime }\right)\right|\rangle =2k_{\mathrm{B}}T{L}^2{\Lambda }_q^{-1}{\delta }_{q{q}^{\prime }}\delta \left(t-t^{\prime }\right)$, and $F_q\left(t\right)$ is the deterministic force derived from potential energies which will be described in detail later. Numerical integration is used to update $h_q$ in Equation (2) from time t to t + Δt, and $h\left(x\right)$ is obtained by inverse Fourier transform from $h_q$. At the end of each time step, $k_{on}^0$ is calculated for each open bond if ligand–receptor distance is within the encounter distance R_RL, and $k_{off}=k_{off}^0\exp \left(f{x}_b/k_{\mathrm{B}}T\right)$ is computed for each closed bond with given bond force f calculated based on the membrane position. Then, $k_{on}^0\Delta t$ and $k_{off}\Delta t$ are compared with a uniformly distributed random number in [0,1] to determine whether binding or unbinding occurs.

The deterministic force $F_q$ in Fourier space is computed from $F_q={\displaystyle {\int }_{A}P\left(x\right)\exp \left(-iq\cdot x\right)dx}$, where $P\left(x\right)=-\delta \widehat{E}/\delta h$ is the force per unit area acting on the membrane and $\widehat{E}$ is the potential or interaction energy per unit area. Next we define all the potential energies involved in the membrane–substrate system illustrated in Figure 1 as follows. First, the membrane bending energy is described by the Canham–Helfrich fluid membrane model [33,34]: $E_b$ = $(1/2){\displaystyle {\int }_{A}B{\left({\nabla }^{2}h\right)}^{2}dx}$, where B denotes bending rigidity, and ${\nabla }^2$ is the Laplace operator in 2D. Thus, the force derived from $E_b$ is $F_q^{\mathrm{b}}=-B{q}^4{h}_q$. Second, we define the equivalent spring constant k of the membrane–substrate linkage (MSL) that takes three components into account: shear deformation resistance in the thickness direction of the lipid bilayer membrane, receptor–ligand complex stiffness, and the substrate elasticity. In this paper, we model the linkage as a linear spring, which is a simplified case for biomaterials that are usually nonlinear. Thus, the potential energy stored in MSL can be written as $E_L$ = ${\displaystyle \sum _{m=0}^{N_b-1}{\displaystyle \sum _{n=0}^{N_b-1}{\displaystyle {\int }_A\frac{1}{2}k{\left(h-h_L\right)}^2\delta \left(x_{mn}-x\right){\varphi }_{mn}dx}}}$, where h_L denotes the vertical position of rest unbound ligands, $x_{mn}=\left(m\Delta ,n\Delta \right)$ are uniform locations of ligand–receptor pairs where Δ denotes receptor spacing as shown in Figure 1. The receptor density is then characterized by 1/Δ². The total number of receptors in the domain of interest is denoted by ${\left(N_b\right)}^2$. The binary function φ_mn describes ligand–receptor binding status with φ_mn = 1 for closed bonds and φ_mn = 0 for open bonds. The force derived from $E_L$ is $F_q^L=-{\displaystyle \sum _{m=0}^{N_b-1}{\displaystyle \sum _{n=0}^{N_b-1}k\left[h\left(x_{mn}\right)-h_L\right]\exp \left(-iq\cdot x_{mn}\right)}}$. Since the substrate is a solid bulk material, the thermal fluctuation of the substrate shape is neglected in this model. Therefore, unbound ligands embedded in the substrate are considered to be at rest with the substrate. The elevated ligands elastically restore to their rest positions when unbinding from receptors and such retracting process takes no time by assuming the time scale of elastic restoring is small compared to membrane fluctuation time scale. The substrate is modeled using a soft-wall repulsive interaction $E_w$ = ${\epsilon }_w{\displaystyle {\int }_A{\left({\sigma }_w/(h+{\sigma }_w-h_{sub})\right)}^{8}dx}$, and E_w is truncated so that $E_w\equiv 0$ when h > h_sub where h_sub is the upper bound of the repulsive interaction range, h_sub = 0 is assumed in this paper. Here, σ_w and ε_w are parameters determining the repulsive interaction strength. Note that other types of enthalpic repulsive interactions [35], if considered, can also be absorbed into E_w. The force acting on the membrane derived from E_w is

$$F_q^w={\displaystyle {\int }_A\left[8\left({\epsilon }_w/{\sigma }_w\right){\left({\sigma }_w/(h\left(x\right)+{\sigma }_w-h_{sub})\right)}^9\right]\exp \left(-iq\cdot x\right)dx}$$

(3)

3. Results

To validate the implementation of the FSBD simulation method, we first carried out free membrane simulations (integrating Equation (2) when $F_q\left(t\right)$ = $F_q^{\mathrm{b}}$), for which a simple analytical expression (the equation inset in Figure 2) exists for the fluctuation spectrum. As shown in Figure 2a, the simulation result (square dots) matches with the theoretical curve (solid line) very well. The equation $\sqrt{<h^2(x)>}\approx \sqrt{A{k}_{\mathrm{B}}T/(4{\pi }^{3}B)}$ [36] theoretically predicts the mean fluctuation in physics space as ~33 nm for B = 20 pN·nm and $\sqrt{A}=L$ = 800 nm, which is comparable to the fluctuation range shown in the inset. The relative mean fluctuation amplitude $\sqrt{<h^2(x)>}/L\approx \sqrt{k_{\mathrm{B}}T/(4{\pi }^{3}B)}$ ~4%. The integration time step is chosen so that it is two orders of magnitude smaller than the relaxation time of the deformation mode of the smallest wavelength. For λ_min = 20 nm, B = 20 pN·nm, and η = 0.06 Poise, the relaxation time is estimated as [37] $\tau =4\eta {\lambda }^3/B{\left(2\pi \right)}^3$~40 ns.

Membrane fluctuation can induce an entropic repulsive force between the membrane and the substrate due to the confinement of the fluctuation [38]. Thus, the membrane is pushed away from the substrate by the entropic pressure. Consequently, the membrane is pinned closer to the substrate on a stiffer substrate by the MSLs for the same amount of closed bonds. We will in the later section show that it is this entropic pressure that gives rise to different closed bond ratios when substrate rigidity is varied. Therefore, it is important to revisit the classical problem of how the entropic pressure is quantitatively dependent on the membrane–substrate distance d using the FSBD simulations. Here, instead of considering a fluctuating membrane between two hard walls, we simulate a fluctuating membrane between a soft wall and an externally applied pressure p that pushes the membrane against the soft wall. A soft wall is used to make the time–space FSBD numerical simulations applicable. Upon equilibrium, the magnitude of the applied pressure p will be equal to that of the entropic pressure. Therefore, we use p to represent the entropic pressure in this section. If the membrane is infinitely large and the hard wall is assumed, the only length scale involved is the distance d. A simple dimensional analysis yields the following relation between p and d,

$$p=c\frac{{\left(k_{\mathrm{B}}T\right)}^2}{B{d}^3}$$

(4)

where c is the coefficient to be determined [36,39]. Substantial work has been devoted to determine the coefficient c (see ref. [40] and also refs cited therein). Here, we carry out the FSBD simulations by integrating Equation (2) with deterministic forces $F_q\left(t\right)$ = $F_q^{\mathrm{b}}$ + $F_q^w$+$F_q^p$, where $F_q^p=-L^{2}p{\delta }_{q,0}$ to determine the coefficient c. Three more length scales are introduced in our simulations compared to the ideal case: the upper cutoff of the Fourier mode wavelength, which is side length L, the lower cutoff wavelength λ_min, and σ_w from the soft wall interaction.

We investigated how the simulation result deviates from Equation (4) when these three length scales are varied. We found that the p – d relation deviates from Equation (4) drastically when λ_min is varied. This is rather counterintuitive, because small wavelength modes with relatively small fluctuation seem to be less affected by the hard wall confinement. Varying membrane periodic size from 800 nm to 240 nm only results in negligible changes in the p – d relation as indicated by the overlapping three curves, as shown in Figure 2b. In addition, increasing the length scale parameter σ_w renders the wall potential softer and thus causes the p – d relation to deviate from Equation (4). By fitting the simulation data that is closest to Equation (4), we obtain c ~ 0.078, which is in good agreement with other simulation results (see the references cited in [40]). Note that, due to the entropic repulsive force, a weak pressure p₀ pushing the membrane towards the substrate is necessary in our FSBD-MC simulations below to prevent the membrane from diffusing far away from the substrate in dynamic processes.

Next, the FSBD-MC simulations were carried out to ascertain the effect of the substrate rigidity and the receptor density on the membrane adhesion. In this case, $F_q\left(t\right)$ = $F_q^{\mathrm{L}}$ + $F_q^{\mathrm{b}}$ + $F_q^p$ + $F_q^w$ is used in Equation (2). We first briefly summarize the value of the parameters used in the simulations. The second step in Equation (1) in the case of zero bond force is a two-state reaction and the equilibrium constant $k_{on}^0/k_{off}^0$ is approximated to be $\exp \left({\epsilon }_b/k_{\mathrm{B}}T\right)$ [41], where ε_b is the ligand–receptor binding energy. Assuming ε_b to be 10 k_BT yields $k_{on}^0/k_{off}^0$~10⁴. In the literature, $k_{on}^0$ was estimated to be from 10⁻³ to 1 ns⁻¹ for hapten-antibody [23], and $k_{off}^0$ was measured to be from 10⁻¹¹ to 10⁻⁶ ns⁻¹ for integrins [42,43]. Since varying the substrate Young’s modulus is equivalent to changing the MSL spring constant k if other segments of the linkage remain unchanged, the lumped spring stiffness k is used to represent the change of substrate stiffness. Typical morphologies of the fluctuating membrane in our FSBD-MC simulations are shown in Figure 3.

The closed bond ratio ϕ as a function of the spring constant k for different $k_{on}^0$ and $k_{off}^0$ is plotted in Figure 4a. In all cases, ϕ first increases with k and then decreases. Such a general behavior can be interpreted as a result of the competition between the change of the rebinding rate and the unbinding rate. Qualitatively, the rebinding rate decreases with the average membrane fluctuation distance, which is characterized here by the average distance d between the membrane and the substrate. The distance d, as shown on the right Y-axis of Figure 4a, decreases with the linkage stiffness k, which implies that the membrane is pinned more closely on the stiffer substrate. Therefore, the rebinding rate is larger on a more rigid substrate, which accounts for the increasing region of the ϕ-k curves. On the other hand, the linkage stiffness k also affects the receptor–ligand bond force. Assuming the bond force is induced by the thermal energy, which is constant, the bond force f scales with k as $f\propto \sqrt{k_{\mathrm{B}}T\times k}$. Thus, the bond force increases with k, as does the unbinding rate $k_{off}=k_{off}^0\exp \left(f{x}_b/k_{\mathrm{B}}T\right)$. When the increase in the unbinding rate outperforms the rebinding rate, the equilibrium closed bond ratio decreases, which corresponds to the decreasing region of the ϕ-k curves. Note that such a substrate rigidity-dependent adhesion scenario only exists in a certain parameter range. As for diffusion-limited reactions, in order for the membrane fluctuation to affect the overall reaction in Equation (1), the rates of the second step need to be much faster than the first step. Therefore, in order to see the substrate rigidity effect on the membrane adhesion, we picked the values of $k_{on}^0$ and $k_{off}^0$ at their upper bound. As shown in Figure 4a for smaller $k_{on}^0$ and $k_{off}^0$, the rigidity dependence is weakened and even disappear based on our other simulation data that is not shown here. Simpler than the rigidity dependence, the closed bond ratio ϕ monotonically decreases with the receptor distance Δ, as shown in Figure 4b. This is obvious because lowering the receptor density directly reduces the maximum closed bonds per unit area and thus results in a greater membrane–substrate distance subject to the entropic repulsive pressure. Therefore, with other parameters remaining constant, a solely decreasing receptor density leads to a smaller rebinding rate and thus a lower closed bond ratio.

4. Conclusions

The simulation results presented in this paper show that when thermally excited membrane fluctuation is taken into account in the specific membrane–substrate adhesion, the adhesion strength characterized by the closed bond ratio becomes dependent on the substrate rigidity and receptor density. The fluctuation induced entropic repulsive force provides the means for the system to “probe” these variations. Such dependence is most evident when the “rate” of the membrane fluctuation step dictated by the viscosity of the surrounding fluid and the membrane bending rigidity is much slower than that of the single receptor–ligand binding and unbinding process. For the classical particle diffusion, the diffusion step only changes the forward and reverse reaction rates but not the overall equilibrium constant. However, it was shown here that the membrane diffusion does affect the overall equilibrium and thus the closed bond ratio. For example, there exists a crossover value of the substrate rigidity at which the closed bond ratio is maximal, and the closed bond ratio monotonically increases with the receptor density. In addition, it is interesting to point out that the membrane fluctuation always works to reduce the closed bond ratio compared to a zero-temperature membrane because, in the absence of membrane fluctuation, the equilibrium constant in Equation (1) $K$ = $k_{on}^0/k_{off}^0$, which is maximal. Therefore, as shown in our additional simulations that are not presented in this paper, when the membrane fluctuating is suppressed or reduced by, for example, increasing the membrane bending rigidity or increasing the membrane tension [16], the adhesion strength increases.