# Undulation of a moving fluid membrane pushed by filament growth

Hiroshi Noguchi, Olivier Pierre-Louis

<figure><img src="image_1.png" /></figure>

## To cite this version:

Hiroshi Noguchi, Olivier Pierre-Louis. Undulation of a moving fluid membrane pushed by filament growth. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-87073-6. hal-03234485

HAL Id: hal-03234485

https://hal.science/hal-03234485

Submitted on 25 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN

# Undulation of a moving fluid membrane pushed by filament growth

Hiroshi Noguchi<fnref n="1" />,²✉ & Olivier Pierre-Louis²

Biomembranes experience out-of-equilibrium conditions in living cells. Their undulation spectra are different from those in thermal equilibrium. Here, we report on the undulation of a fluid membrane pushed by the stepwise growth of filaments as in the leading edge of migrating cells, using three-dimensional Monte Carlo simulations. The undulations are largely modified from equilibrium behavior. When the tension is constrained, the low-wave-number modes are suppressed or enhanced at small or large growth step sizes, respectively, for high membrane surface tensions. In contrast, they are always suppressed for the tensionless membrane, wherein the wave-number range of the suppression depends on the step size. When the membrane area is constrained, in addition to these features, a specific mode is excited for zero and low surface tensions. The reduction of the undulation first induces membrane buckling at the lowest wave-number, and subsequently, other modes are excited, leading to a steady state.

In nonequilibrium, surfaces and interfaces often exhibit different fluctuations from thermal equilibrium. The fluctuations of growing surfaces have been studied by the Kardar-Parisi-Zhang equation and other models<fnref n="1" />⁻⁷. Scaling laws for the spatial and temporal evolution of surface roughness and correlation functions have been intensively discussed in various physical systems. On the other hand, larger undulations of fluid membranes have been observed out of equilibrium than in equilibrium⁸⁻¹¹. For example, red blood cells exhibit non-thermal fluctuations¹⁰. However, nonequilibrium membrane undulations have only been investigated under limited conditions.

In living cells, membranes often interact with protein filaments (e.g., actin¹²,¹³ and microtubules¹⁴,¹⁵), which push and/or pull the membranes. For example, actin filaments grow at the front side of crawling cells and push the membrane forward. This membrane motion has been intensively studied by Brownian ratchet theories and simulations¹⁵⁻²⁸. However, membrane fluctuations have not yet been studied. In most studies reported in the literature¹⁶⁻²⁴, the membrane is modeled as a flat rigid surface; thus, no fluctuations are accounted for. Although the membrane fluctuations are assumed to be a Gaussian distribution of the membrane height in Ref.²⁵, the spatial correlation was not accounted for. Membrane fluctuations were also considered with a one-dimensional lattice model²⁶⁻²⁸, but the surface tension was treated by a rough approximation of the rectangular contour length in the absence of the bending energy. Thus, the effects of the bending rigidity have not yet been investigated.

In this study, we report on a minimal model that describes the fluctuations of a fluid membrane with the bending rigidity in interaction with growing filaments. We examine the fluctuations of the membrane pushed by filament growth using Monte Carlo (MC) simulations. The filaments grow by a stepwise random walk under the membrane and have an excluded-volume interaction with the membrane. The analysis of the fluctuation spectrum reveals that membranes deviate from the well-known equilibrium behavior under either condition in which the bending energy or surface tension is dominant. Fluctuations are enhanced or suppressed depending on the conditions. Moreover, enhancement of a specific wave-number is obtained when the membrane area is constrained.

## Results

The filaments are modeled as a square array of columns growing vertically in the *z* direction by the addition of new growth units of $\Delta z_{fil}$. The filaments grow only upward in the +*z* direction with a probability *p*<sub>fil</sub>, as shown in Fig. 1. The filament retraction and interactions between filaments are not considered for simplicity. The fluid

<footnote n="1">¹Institute for Solid State Physics, University of Tokyo, Kashiwa, Chiba 277-8581, Japan. ²Institut Lumière Matière, UMR5306 Université Lyon 1-CNRS, Université de Lyon, 69622 Villeurbanne, France. ✉email: noguchi@issp.u-tokyo.ac.jp</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Motion of the membrane and filaments. (a) Schematic of the membrane and filaments. According to the Metropolis MC method, the membrane moves with a vertical (z) step taken from a uniform random number in [-Δz<sub>mb</sub>, Δz<sub>mb</sub>]. The filament grows stepwise with a probability p<sub>fil</sub>. (b) Snapshot of a part of the membrane and filaments at Δz<sub>fil</sub> = 0.4, p<sub>fil</sub> = 0.5, γ = 0, and N = 4096 in case I (tension constraint). The red mesh and blue bars represent the membrane and filaments, respectively. The origin of the z coordinate is taken at the minimum value of the filament position z<sub>fil</sub>.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Spectra of undulation modes $\langle |h(q)|^2 \rangle$. (a) Membrane under thermal equilibrium at $\gamma = 0, 1$, and 10. The symbols (×, ∇, and +) and (◇, △, and □) represent the spectra in cases I (tension constraint) and II (area constraint), respectively. Solid lines show $1/(\kappa q^4 + \gamma q^2)$ with $\kappa = 10$. (b, c) The undulation spectra of the membrane (red diamonds) and filament tips (blue cross marks) of tensionless membranes ($\gamma = 0$) at $\Delta z_{\text{fil}} = 0.004$ for (b) case I and (c) case II. The error bars are smaller than the symbol size.</figcaption></figure>

membrane is described by a continuous height z, defined on a square mesh above each filament, and is moved stochastically upward and downward by the Metropolis MC method. Unless otherwise specified, the number N = 64² = 4096 of filaments (as well as membrane vertices) and the filament growth probability p<sub>fil</sub> = 0.5 are used (see “Methods” for more details). The results are displayed with the thermal energy k<sub>B</sub>T and lateral distance between neighboring filaments l<sub>fil</sub> for the energy and length units, respectively.

**Thermal equilibrium.** In thermal equilibrium (no interactions with the filaments), the height z = h(x, y) of the membrane exhibits a static spectrum that is controlled by the bending energy and tension<sup>29–32</sup>:

$$ \langle |h(q)|^2 \rangle = \frac{1}{\kappa q^4 + \gamma q^2}, \qquad (1) $$

where q is the norm of the wave-vector in the x, y plane, κ is the bending rigidity, and γ is the mechanical surface tension<sup>33</sup> conjugated to the projected membrane area A<sub>xy</sub> = N, respectively (see Fig. 2a). For low wave-numbers (q ≲ √γ/κ), the tension is the dominant factor in determining membrane undulation, while the bending rigidity is at high wave-numbers (q ≳ √γ/κ). In this study, κ = 10 is used. Thus, for γ = 10, surface tension is dominant for most of the wave-number range; for γ = 1, the tension and bending energies are dominant for low and high wave-numbers, respectively. To maintain the surface tension or membrane area, we employ two types of constraints: case I (tension constraint), in which the intrinsic surface tension γ<sub>it</sub> conjugated to the real membrane area A is imposed by the potential U<sub>it</sub> = γ<sub>it</sub>A, and case II (area constraint), in which area A is
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Mean growth velocity ($v_z$) and mean distance between the membrane and filament tips ($z_{gap}$). In (a) and (b), these quantities are plotted as a function of $\Delta z_{fil}$ for $\gamma = 0, 1,$ and 10. The solid and dashed lines represent the data in cases I (tension constraint) and II (area constraint), respectively. In (c) and (d), the data around the maximum velocity and minimum distance are enlarged, and the mean filament growth distance $p_{fil} \Delta z_{fil}$ normalized by the mean distance $\Delta z_{m0}$ of the membrane motion at one MC step. The error bars are smaller than the line thickness.</figcaption></figure>

constrained by a harmonic potential $U_{ar} = (K_A/2A_0)(A - A_0)^2$, where $K_A$ is the area compression modulus. In these cases, $\gamma_{it}$ or $A_0$ is adjusted to obtain the spectrum given by Eq. (1). Note that $\gamma_{it}$ is slightly greater than $\gamma$ because of the membrane undulation as discussed in the literature$^{33-36}$. In case I, the membrane is a part of a cell membrane or liposome and the rest part of the membrane can act as a membrane reservoir to change the area, while no reservoir exist in case II. Both types of constraints lead to the same static spectrum as shown in Fig. 2a.

**Steady-state velocity.** As the filaments push the membrane upward, the membrane and filament tips relax to a steady-state, in which they move at the same speed ($v_z$). The growth speed ($v_z$) exhibits a non-monotonic dependence on the size of the filament growth units $\Delta z_{fil}$. For intermediate values of $\Delta z_{fil}$, a maximum speed is reached concomitantly with a minimum in the average of the gap distance $z_{gap} = z_{mb} - z_{fil}$ between the membrane and filament tips, as shown in Fig. 3. Hence, the averaged motion is not sensitive to the surface tension $\gamma$ or to the constraint types (case I or II). In all cases a similar dependence on $\Delta z_{fil}$ is found.

The two asymptotic regimes of large and small $\Delta z_{fil}$ present distinct behaviors. For a large step of $\Delta z_{fil}$, the mean growth velocity ($v_z$) and the gap distance ($z_{gap}$) are independent of $p_{fil}$ and determined by $\Delta z_{fil}$, as shown in Fig. S1, because frequent growth trials of the filaments are rejected until the membrane moves by a sufficiently large distance for the filament growth units to be inserted. This corresponds to the diffusion limit in the Brownian ratchet model$^{16,17}$. On the other hand, the behavior in response to small steps of $\Delta z_{fil}$ is predominantly determined by the ratio of the filament growth rate $p_{fil} \Delta z_{fil}$ to the mean membrane step size $\Delta z_{m0} = p_{mb} \Delta z_{mb}/2$, where $p_{mb}$ is the mean acceptance ratio of the membrane motion in the absence of filaments ($p_{mb} \approx 0.5$ in our simulations, see “Methods”). The velocity ($v_z$) for $p_{fil} = 0.25, 0.5,$ and 1 merge into one curve when they are plotted with the horizontal axis of $p_{fil} \Delta z_{fil} / \Delta z_{m0}$ (see Fig. S1d), whereas a small difference remains in ($z_{gap}$) (see Fig. S1e). The velocity approaches the growth velocity of free filaments, $p_{fil} \Delta z_{fil}$, for $\Delta z_{fil} \to 0$, so that this corresponds to the reaction limit in the Brownian ratchet model$^{16,17}$. The maximum velocity and minimum distance are obtained at $p_{fil} \Delta z_{fil} \sim \Delta z_{m0}$, i.e., when the motion of the membrane opens gaps corresponding to the size of filament growth units at a frequency which is similar to that of the insertion of new growth units (see Fig. 3c and d as well as Fig. S1d and e). Furthermore, irrespective of the average value of the gap, the probability distribution $P(z_{gap})$ of the gap $z_{gap}$ decreases monotonically with increasing $z_{gap}$ and exhibits a stepwise discontinuity at $z_{gap} = \Delta z_{fil}$, because the filament growth is rejected at a smaller distance (see Fig. S2).

**Undulation spectrum.** The pushed membrane exhibits a different undulation spectrum from that of equilibrium. Figure 2b and c show the spectra of the tensionless membrane ($\gamma = 0$) at $\Delta z_{fil} = 0.004$. At a low wave-number $q$ (i.e., long wavelength), the membrane and filament surface connecting their tips have identical spectra. On the other hand, at a high $q$ (i.e., short wavelength), the membrane spectra are not modified from the equilibrium spectrum in Fig. 2a, and the filament surface exhibits a flat spectrum of white noise, as in the
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Undulation spectra normalized by the equilibrium values $\langle|h(q)|^2\rangle/\langle|h_{eq}(q)|^2\rangle$ for $\Delta z_{fil} = 0.002, 0.01, 0.06$ (or 0.065), 0.4, and 0.8. (a–c) Case I (tension constraint) and (d–f) case II (area constraint) for (a, d) $\gamma = 0$; (b, e) $\gamma = 1$; and (c, f) $\gamma = 10$. The blue triangles represent the data at $\Delta z_{fil} = 0.06$ or 0.065, where the distance $\langle z_{gap} \rangle$ has the minimum value. The black pluses represent the data of the membrane with the totally asymmetric diffusion (tamb). The error bars are smaller than the symbol size except for low wave-numbers at small $\Delta z_{fil}$ in case II.</figcaption></figure>

absence of membrane-filament interactions. The threshold $q_{sep}$ of separation of the spectra of the membrane and filament surface is correlated to the gap distance. With increasing $\Delta z_{fil}$, $q_{sep}$ increases, reaches a maximum when $\langle z_{gap} \rangle$ is minimum, and then decreases (see Fig. S3). Up to the maximum of $q_{sep}$, the filament spectrum at high $q$ is flat. However, above the maximum, the filament spectrum is not flat at high $q$; thus, the spectrum is modified by the membrane-filament interactions (see Fig. S3c). These high $q$ behaviors are not qualitatively changed by the choice of constraints (cases I or II). However, a distinct difference is found in the spectra of low $q$. In case I (tension constraint), the membrane undulation are suppressed at low $q$. In contrast, the spectrum is enhanced at $q \approx 0.07\pi$ in case II (area constraint; compare Fig. 2b and c). Hence, the membrane buckles to maintain the membrane area A. A similar nonequilibrium buckling scenario resulting from the suppression of membrane fluctuations has been discussed as a mechanism that induces instability of the lamellar phase in shear flow³⁷. Note that the membrane-filament interactions cannot be interpreted by the effective change of the surface tension, since the spectrum shapes in Fig. 2b and c are different from that for any value of $\gamma$.

To further examine the spectrum changes, Fig. 4 shows the membrane spectrum normalized by the equilibrium spectrum as $\langle|h(q)|^2\rangle/\langle|h_{eq}(q)|^2\rangle$. Let us first discuss case I, with tension constraint. For tensionless membranes low-$q$ undulations are suppressed for all values of $\Delta z_{fil}$, while it is slightly enhanced at the intermediate wave-number ($q \approx 0.3\pi$) for large $\Delta z_{fil}$ (Figs. 4a and 5a). A rough measure of the range of $q$ for which the fluctuations are suppressed (i.e., $\langle|h(q)|^2\rangle/\langle|h_{eq}(q)|^2\rangle < 1$) during growth can be obtained from the evaluation of the value of $q$ for which $\langle|h(q)|^2\rangle/\langle|h_{eq}(q)|^2\rangle = 0.6$. This value of $q$ plotted as the solid lines in Fig. 5c, is seen to exhibit a maximum for $\Delta z_{fil} = 0.02$. Hence, the low $q$ range where the modes are suppressed is widest around $\Delta z_{fil} = 0.02$. This value is smaller than that corresponding to the minimum of the gap distance $\Delta z_{fil} = 0.065$. This maximum range of suppressed modes changes only slightly (see Fig. 5c) at higher $\gamma$. However, when $p_{fil}$ is increased, the maximum of the suppression range increases while the corresponding value of $p_{fil}\Delta z_{fil}$ is constant (see Fig. S4). In addition to this suppression effect at small $\Delta z_{fil}$, an enhancement effect is seen for large $\Delta z_{fil}$ for high tensions. For $\gamma = 10$ the undulation increases at $q \ll 1$ with increasing $\Delta z_{fil}$ (see Figs. 4 and 5a), and a steep increase is obtained at $\Delta z_{fil} \sim 1$.

As a summary, for increasing $\Delta z_{fil}$, there are first a range of modes that are suppressed. This range reaches a maximum for some value of $\Delta z_{fil}$ (here $p_{fil}\Delta z_{fil} = 0.01$), and then decreases. Upon a further increase in $\Delta z_{fil}$, the low-$q$ modes are enhanced for non-vanishing tensions.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Dependence of the spectrum shape on $\Delta z_{\text{fil}}$. (a, b) Normalized undulation amplitude $\langle|h(q)|^2\rangle/\langle|h_{\text{eq}}(q)|^2\rangle$ at the lowest wave-number $q = 0.03125\pi$ for (a) case I (tension constraint) and (b) case II (area constraint). (c) The upper four lines represent the wave-number $q$ at $\langle|h(q)|^2\rangle/\langle|h_{\text{eq}}(q)|^2\rangle = 0.6$. The solid lines represent the data at $\gamma = 0, 1$, and $10$ in case I from top to bottom. The upper and lower dashed lines represent the data at $\gamma = 0$ in case II for $\langle|h(q)|^2\rangle/\langle|h_{\text{eq}}(q)|^2\rangle = 0.6$ and the peak of $\langle|h(q)|^2\rangle/\langle|h_{\text{eq}}(q)|^2\rangle$, respectively.</figcaption></figure>

Under the area constraint (case II), a similar trend of suppression of low $q$ modes with a maximum at small $\Delta z_{\text{fil}}$ together with an enhancement of low $q$ modes at large $\Delta z_{\text{fil}}$ is found. However, there is an important superimposed feature: a peak at low $q$ appears in the spectrum ratio for small $\gamma$. The peak can be clearly seen in Fig. 4d for $\gamma = 0$. The position of this peak, shown as the lower dashed line in Fig. 5c, shifts to a higher $q$ with increasing $\Delta z_{\text{fil}}$. For $\gamma = 1$, the undulations are similar to those of $\gamma = 10$ and $\gamma = 0$ in cases I and II, respectively. Globally, the membrane undulation is strikingly changed by interactions with the filaments and strongly depends on the filament growth rate and membrane constraint types.

Moreover, we examined a totally asymmetric membrane (tamb) motion as a limited condition, in which the filaments are in complete contact with the membrane, and all of the downward membrane motion are rejected. The tamb spectra are close to those of the minimum gap distance, as shown in Fig. 4. Thus, the membrane under the minimum gap distance condition is well approximated by the asymmetric membrane motion. In this limit, the undulation due to the bending energy is modified ($\gamma = 0$ and $1$), while that due to the tension is not ($\gamma = 10$). The tamb steady velocity $\langle v_z \rangle$ is $\simeq 50\%$ higher than that of the minimum gap distance condition: $\langle v_z \rangle = 0.0183$ (0.0190), 0.0180 (0.0182), and 0.0162 (0.0164) for $\gamma = 0, 1$, and $10$ in case I (case II), respectively. This difference is reasonable as a smaller minimum distance gives a higher velocity, as shown in Fig. S1.

**Vertical span.** In the studies of surface growth, the time evolution and finite-size scaling of the surface thickness (or width) have been analyzed$^{1–6}$. Here, we call this the surface vertical span to avoid confusion with the thickness of the membrane itself, and defined it as $z_{\text{span}}^2 = \sum_{i}^{N} (z_i - z_G)^2/N$, with $z_G = \sum_{i}^{N} z_i/N$. Due to Parseval’s identity, it corresponds to the sum of the spectrum over all modes.

Figure 6 shows $\langle z_{\text{span}}^2 \rangle$ as a function of $\Delta z_{\text{fil}}$ and the system size $N$. For very small or very large $\Delta z_{\text{fil}}$, the filament vertical span is greater than the membrane span, reflecting the difference of the undulation spectra at high $q$. The $\Delta z_{\text{fil}}$ dependence of $\langle z_{\text{span}}^2 \rangle$ is roughly captured by that of the undulation spectra at a low $q$ (compare Figs. 5a, b and 6a, b). At a small $\Delta z_{\text{fil}}$, $\langle z_{\text{span}}^2 \rangle$ provides slightly large and small values for $p_{\text{fil}} = 0.25$ and $1$, respectively, while they merge well at large $\Delta z_{\text{fil}}$ (see Fig. S1c and f). Under thermal equilibrium at $\gamma = 0$, the membrane vertical span linearly increases as $\langle z_{\text{span}}^2 \rangle \sim N$ following the amplitude of the lowest undulation mode$^{30}$. For the pushed membrane, this increase is reduced for all values of $\Delta z_{\text{fil}}$ (see Fig. 6c). On the other hand, for $\gamma = 10$, a greater increase is obtained at large $\Delta z_{\text{fil}}$ (see Fig. 6d). This observation reflects an amplitude increase in the lowest mode as seen in Fig. 5a.

The excess membrane area increases as $(A)/A_{xy} - 1 = (k_B T/8\pi\kappa) \ln(N) + b$ for the tensionless membrane ($\gamma = 0$) under thermal equilibrium$^{30,38,39}$, where $b$ is a constant. The modification of this size dependence by the filament growth is similar to that in the vertical span, and the excess area increase is largely reduced at $\gamma = 0$ (see Fig. S5). Note that the mean velocity $\langle v_z \rangle$ and gap distance $\langle z_{\text{gap}} \rangle$ exhibit negligibly small size dependence.

**Membrane dynamics.** Finally, we describe the time evolution. After the filaments contact the membrane, both surfaces relax into steady states. The choice of the initial conformations, such as a flat or thermal equilibrium membrane conformation does not lead to any notable difference in the subsequent dynamics. In most of the conditions, the undulations increase monotonously on average despite large fluctuations, as shown in Fig. 7b. However, a characteristic evolution is found at $\Delta z_{\text{fil}} \lesssim 0.01$ for $\gamma = 0$ under the area constraint (see Fig. 7a and c). First, the membrane and filament surface buckle together and form a bump, as shown in the middle snapshot in Fig. 7a. This bump leads to a strong peak in the first mode of the spectrum $|h(q)|^2$ and in $\langle z_{\text{span}}^2 \rangle$. Subsequently,
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Vertical span $\langle z_{\text{span}}^2 \rangle$ of the membrane and filament surfaces. (a, b) Dependence on $\Delta z_{\text{fil}}$ for $\gamma = 0, 1$, and 10 in (a) case I and (b) case II. The solid and dashed lines represent the data for the membrane and filaments, respectively. (c, d) Dependence on system size $N$ for (c) $\gamma = 0$ and (d) $\gamma = 10$ in case I. The thick magenta lines in (c) and (d) represent the data at thermal equilibrium. The dashed black line in (c) indicates $\langle z_{\text{span}}^2 \rangle = N/8\pi^3\kappa$, which overlays the thick magenta line.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Initial time evolution of the membrane and filaments at $\gamma = 0$ and $\Delta z_{\text{fil}} = 0.004$. (a) Sequential snapshots corresponding to the data in (c). (b, c) Time evolution of the amplitudes of the first, second, and third modes along the x-axis and the vertical span $z_{\text{span}}$ in (b) case I and (c) case II. The solid and dashed lines represent the data of the membrane and filaments, respectively. The time unit is MC step. In both cases, a thermal equilibrium membrane conformation is used as an initial state; the initial filament surface is flat and contacts to the minimum position of the membrane.</figcaption></figure>

the higher modes develop, leading to a steady state (the bottom snapshot in Fig. 7a). The initial suppression of the membrane undulation induces this overshoot buckling, because the lowest mode can evolve the fastest. It is known that a similar buckling at the lowest mode is induced by a negative surface tension⁴⁰, so that it may be interpreted that an effective negative tension is yielded by the interaction with the filaments. Such buckling could share also similarities with wrinkles forming when a membrane is confined between two walls⁴¹, which are also be suppressed by tension⁴².
Discussion

We numerically studied the membrane pushed by filament growth. The growth velocity has a maximum at a slightly larger filament-growth step $\Delta z_{\text{fil}}$ than for the minimum gap distance; they are not sensitive to the surface tension and the constraint types. It is found that the membrane undulation spectrum is non-monotonously changed from that of thermal equilibrium. Under the tension constraint, the low-wave-number (long wavelength) undulations are suppressed for the tensionless membrane, and the range of the suppression displays a maximum at a value of $\Delta z_{\text{fil}}$ that is smaller than that corresponding to the maximum velocity and minimum gap. At a high surface tension, this suppression is altered to enhancement with increasing $\Delta z_{\text{fil}}$. Under the membrane-area constraint, we find similar features as in the tension constraint, but with two main differences for low tensions. First, a peak appears at some intermediate wave-number. Second, the membrane dynamics of relaxation to the steady-state is changed. The suppression of the undulation initially induces the buckling of the membrane at the lowest wave-number to maintain the membrane area. Subsequently, other modes are excited. The spectrum of the filament surface is identical to that of the membrane for low wave-numbers but deviates for high wave-numbers. The smaller the gap distance, the wider the identical region. Consequently, the vertical spans of the membrane and filaments deviate at small or large filament-growth steps.

Let us map the present model parameters to the membrane systems in living cells. In the model, $1/\ell_{\text{fil}}^2$ is the lateral density of the filaments under the membrane. The membrane undulation is cut off at a membrane thickness of approximately 5 nm, so that the available range is $\ell_{\text{fil}} \gtrsim 5$ nm. Our simulations clarified that the membrane undulation modes at high wave-numbers are not modified by the filament growth. Hence, we can conclude that the modes higher than $2\pi/\ell_{\text{fil}}$ should also be unchanged, and we can only consider changes in the lower *q* modes.

An actin network pushes the plasma membrane of the leading edge in migrating cells<sup>12,13</sup>, in which the step size for actin growth is $\Delta z_{\text{fil}} = 2.7$ nm<sup>17,18,25</sup>. The density of actin filaments varies among cell types and with the cell state<sup>12,13</sup>. For high and low densities of $\ell_{\text{fil}} = 10$ nm and 100 nm, $\Delta z_{\text{fil}}/\ell_{\text{fil}} = 0.27$ and 0.027 are obtained, respectively. Thus, at a high density, the discreteness of the filament growth may play a significant role in membrane-filament interactions. A microtubule is a hollow cylinder with a diameter of 25 nm, typically consisting of 13 protofilaments, and the growth unit length is $\Delta z_{\text{fil}} = 8$ nm<sup>14,19</sup>. When the microtubules are closely packed, the average distance is $\ell_{\text{fil}} \simeq 7$ nm, such that $\Delta z_{\text{fil}}/\ell_{\text{fil}} \simeq 1$ at the maximum.

Here, we consider the minimal model for filament growth, in which the membrane and filament tips have only a repulsive interaction and no interactions between filaments. In the previous studies, the contact energy between neighboring filaments<sup>20,21</sup>, attractive potentials between the membrane and filament tips<sup>22</sup>, acceleration of filament growth by the membrane contact<sup>27,28</sup>, and filament rigidity<sup>43</sup> have been investigated. These interactions can be easily added to the present model to clarify their effects on the undulations. Actin filaments often bind to membranes via curvature-inducing proteins<sup>12</sup>. Under such conditions, the filament contact is accompanied by the local induction of a membrane spontaneous curvature. Propagation waves can be generated by the coupling of the curvature-inducing proteins with the actin and/or regulatory proteins<sup>44–46</sup>. The effects of such a spontaneous curvature on membrane undulation are also an interesting topic for further studies. The model presented here could therefore be a versatile tool for the investigation of the interactions that may affect the membrane undulations in nonequilibrium conditions.

In previous studies<sup>8–11</sup> on nonequilibrium membrane fluctuations, membrane undulations are always enhanced by active energy inputs. In the present case, the opposite effect (suppression) is also found. Indeed, filament growth can either increase or decrease the undulations depending on the conditions, and in particular, it can induce an excitation at a specific wave-number. Such undulations could stimulate characteristic lengthscales or periodic structures that may give rise to filopodia and microspike coupling with the filament assembly.

Methods

The fluid membrane is modeled by a squared mesh of $N$ vertices with periodic boundary condition, as described in Ref.<sup>33</sup>. The bending energy is given by $U_{\text{bend}} = \int (\kappa/2)(C_1 + C_2)^2 dA$, where $C_1$ and $C_2$ represent the principal curvatures<sup>29,47,48</sup>. The Monge representation ($z = h(x, y)$) is employed, and the curvature is calculated as $C_1 + C_2 = [(1 + h_x^2)h_{yy} + (1 + h_y^2)h_{xx} - 2h_xh_yh_{xy}]/(1 + h_x^2 + h_y^2)^{3/2}$, where the subscripts represent spatial derivatives, such as $h_x = \partial h/\partial x^{29}$. To control the membrane area, the intrinsic tension $\gamma_{\text{it}}$ or the membrane area $A$ is constrained by the addition of $U_{\text{it}}$ or $U_{\text{ar}}$ to the Hamiltonian in case I or II, respectively. To remove the artificial entropy production by the membrane tilt, a correction potential $U_{\text{cor}} = -k_B T \sum_i \ln(\cos \theta_i)$ is also added to the Hamiltonian, where $\theta_i$ is the angle between the normal vector of the $i$-th site and the $z$-axis<sup>33</sup>. Straight filaments are arranged in the squared lattice ($x_i, y_i$) which is shared by the vertices of the membrane (see Fig. 1). Membrane-filament excluded-volume interactions are implemented by forbidding moves leading to inter-penetration ($z_{\text{fil}} < z_{\text{mb}}$).

In the filament growth step, one of the filaments is randomly selected, and its tip moves upward with a probability $p_{\text{fil}}$ for a step of $\Delta z_{\text{fil}}$. If the filament overlaps with the membrane vertex, the trial is rejected. The membrane vertex is moved by a vertical step taken from a uniform random number in $[-\Delta z_{\text{mb}}, \Delta z_{\text{mb}}]$, and the motion is accepted or rejected by the Metropolis MC procedure. In each MC step, $N$ trials are performed for both the filaments and membrane. In this study, $\Delta z_{\text{mb}} = 0.2$ and $K_A/A_0 = 1$ are used. For $\gamma = 0, 1$, and $10$, $\gamma_{\text{it}} = 0.44, 1.44$, and $10.4$ ($A_0/N = 1.036, 1.022$, and $1.01$ for $N = 4096$ in case II) are used, respectively. In our simulations, we have $p_{\text{mb}} \approx 0.5$. More precisely, $p_{\text{mb}} = 0.5226, 0.5128$, and $0.4810$ for $\gamma = 0, 1$, and $10$ in case I and $p_{\text{mb}} = 0.5228, 0.5130$, and $0.4813$ for $\gamma = 0, 1$, and $10$ in case II, respectively. Error bars are calculated from three independent runs.
Data availability

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Received: 8 February 2021; Accepted: 22 March 2021

Published online: 12 April 2021

References

1. Barabasi, A.-L. & Stanley, H. E. *Fractal Concepts in Surface Growth* (Cambridge University Press, 1995).

2. Family, F. & Vicsek, T. Scaling of the active zone in the Eden process on percolation networks and the ballistic deposition model. *J. Phys. A* **18**, L75–L81 (1985).

3. Kardar, M., Parisi, G. & Zhang, Y.-C. Dynamic scaling of growing interfaces. *Phys. Rev. Lett.* **56**, 889–892 (1986).

4. Sasamoto, T. & Spohn, H. One-dimensional Kardar–Parisie–Zhang equation: An exact solution and its universality. *Phys. Rev. Lett.* **104**, 230602 (2010).

5. Corwin, I. The Kardar–Parisie–Zhang equation and universality class. *Random Matrices Theory Appl.* **1**, 1130001 (2012).

6. Halpin-Healy, T. & Takeuchi, K. A. A KPZ Cocktail-shaken, not stirred.... *J. Stat. Phys.* **160**, 794–814 (2015).

7. Cagnetta, F., Evans, M. R. & Marenduzzo, D. Active growth and pattern formation in membrane–protein systems. *Phys. Rev. Lett.* **120**, 258001 (2018).

8. Prost, J. & Bruinsma, R. Shape fluctuations of active membranes. *Europhys. Lett.* **33**, 321–326 (1996).

9. Manneville, J.-B., Bassereau, P., Ramaswamy, S. & Prost, J. Active membrane fluctuations studied by micropipet aspiration. *Phys. Rev. E* **64**, 021908 (2001).

10. Turlier, H. et al. Equilibrium physics breakdown reveals the active nature of red blood cell flickering. *Nat. Phys.* **12**, 513–519 (2016).

11. Almendro-Vedia, V. G. et al. Nonequilibrium fluctuations of lipid membranes by the rotating motor protein F₁F₀-ATP synthase. *Proc. Natl. Acad. Sci. USA* **114**, 11291–11296 (2017).

12. Svitkina, T. The actin cytoskeleton and actin-based motility. *Cold Spring Harb. Perspect. Biol.* **10**, a018267 (2018).

13. Skruber, K. et al. Arp2/3 and Mena/VASP require profilin 1 for actin network assembly at the leading edge. *Curr. Biol.* **30**, 2651–2664 (2020).

14. Wade, R. H. On and around microtubules: An overview. *Mol. Biotechnol.* **43**, 177–191 (2009).

15. Dogterom, M., Kerssemakers, J. W. J., Romet-Lemonne, G. & Janson, M. E. Force generation by dynamic microtubules. *Curr. Opin. Cell Biol.* **17**, 67–74 (2005).

16. Howard, J. *Mechanics of Motor Proteins and the Cytoskeleton* (Sinauer Associates, 2001).

17. Peskin, C. S., Odell, G. M. & Oster, G. F. Cellular motions and thermal fluctuations: The Brownian ratchet. *Biophys. J.* **65**, 316–324 (1993).

18. Mogilner, A. & Oster, G. The physics of lamellipodial protrusion. *Eur. Biophys. J.* **25**, 47–53 (1996).

19. Mogilner, A. & Oster, G. The polymerization ratchet model explains the force-velocity relation for growing microtubules. *Eur. Biophys. J.* **28**, 235–242 (1999).

20. Stukalin, E. B. & Kolomeisky, A. B. Simple growth models of rigid multifilament biopolymers. *J. Chem. Phys.* **121**, 1097–1104 (2004).

21. Krawczyk, J. & Kierfeld, J. Stall force of polymerizing microtubules and filament bundles. *EPL* **93**, 28006 (2011).

22. Motahari, F. & Carlsson, A. E. Thermodynamically consistent treatment of the growth of a biopolymer in the presence of a smooth obstacle interaction potential. *Phys. Rev. E* **100**, 042409 (2019).

23. Whitehouse, J., Blythe, R. A. & Evans, M. R. Width scaling of an interface constrained by a membrane. *Phys. Rev. Lett.* **121**, 058102 (2018).

24. Wood, A. J., Blythe, R. A. & Evans, M. R. Solvable model of a many-filament Brownian ratchet. *Phys. Rev. E* **100**, 042122 (2019).

25. Lan, Y. & Papoian, G. A. The stochastic dynamics of filopodial growth. *Biophys. J.* **94**, 3839–3852 (2008).

26. Narasimhan, S. L. & Baumgaertner, A. Dynamics of a driven surface. *J. Chem. Phys.* **133**, 034702 (2010).

27. Sadhu, R. K. & Chatterjee, S. Actin filaments growing against a barrier with fluctuating shape. *Phys. Rev. E* **93**, 062414 (2016).

28. Sadhu, R. K. & Chatterjee, S. Actin filaments growing against an elastic membrane: Effect of membrane tension. *Phys. Rev. E* **97**, 032408 (2018).

29. Safran, S. A. *Statistical Thermodynamics of Surfaces, Interfaces, and Membranes* (Addison-Wesley, 1994).

30. Helfrich, W. & Servuss, R.-M. Undulations, steric interaction and cohesion of fluid. *Nuovo Cimento D* **3**, 137–151 (1984).

31. Goetz, R., Gompper, G. & Lipowsky, R. Mobility and elasticity of self-assembled membranes. *Phys. Rev. Lett.* **82**, 221–224 (1999).

32. Shiba, H. & Noguchi, H. Estimation of the bending rigidity and spontaneous curvature of fluid membranes in simulations. *Phys. Rev. E* **84**, 031926 (2011).

33. Shiba, H., Noguchi, H. & Fournier, J.-B. Monte Carlo study of the frame, fluctuation and internal tensions of fluctuating membranes with fixed area. *Soft Matter* **12**, 2373–2380 (2016).

34. David, F. & Leibler, S. Vanishing tension of fluctuating membranes. *J. Phys. II* (1), 959–976 (1991).

35. Farago, O. & Pincus, P. The effect of thermal fluctuations on Schulman area elasticity. *Eur. Phys. J. E* **11**, 399–408 (2003).

36. Gueguen, G., Destainville, N. & Manghi, M. Fluctuation tension and shape transition of vesicles: Renormalisation calculations and Monte Carlo simulations. *Soft Matter* **13**, 6100–6117 (2017).

37. Zilman, A. G. & Granek, R. Undulation instability of lamellar phases under shear: A mechanism for onion formation?. *Eur. Phys. J. B* **11**, 593–608 (1999).

38. den Otter, W. K. Area compressibility and buckling of amphiphilic bilayers in molecular dynamics simulations. *J. Chem. Phys.* **123**, 214906 (2005).

39. Noguchi, H. & Gompper, G. Meshless membrane model based on the moving least-squares method. *Phys. Rev. E* **73**, 021903 (2006).

40. Noguchi, H. Anisotropic surface tension of buckled fluid membrane. *Phys. Rev. E* **83**, 061919 (2011).

41. To, T. B. T., Le Goff, T. & Pierre-Louis, O. Adhesion dynamics of confined membranes. *Soft Matter* **14**, 8552–8569 (2018).

42. Le Goff, T., Politi, P. & Pierre-Louis, O. Transition to coarsening for confined one-dimensional interfaces with bending rigidity. *Phys. Rev. E* **92**, 022918 (2015).

43. Weichsel, J. & Geissler, P. L. The more the tubular: Dynamic bundling of actin filaments for membrane tube formation. *PLoS Comput. Biol.* **12**, e1004982 (2016).

44. Peleg, B., Disanza, A., Scita, G. & Gov, N. Propagating cell-membrane waves driven by curved activators of actin polymerization. *PLoS ONE* **6**, e18635 (2011).

45. Wu, Z., Su, N., Tong, C., Wu, M. & Liu, J. Membrane shape-mediated wave propagation of cortical protein dynamics. *Nat. Commun.* **9**, 136 (2018).

46. Tamemoto, N. & Noguchi, H. Pattern formation in reaction-diffusion system on membrane with mechanochemical feedback. *Sci. Rep.* **10**, 19582 (2020).

47. Canham, P. B. The minimum energy of bending as a possible explanation of the biconcave shape of the human red blood cell. *J. Theor. Biol.* **26**, 61–81 (1970).
48. Helfrich, W. Elastic properties of lipid bilayers: Theory and possible experiments. Z. Nat. 28c, 693–703 (1973).

**Acknowledgements**

This work was supported by JSPS KAKENHI Grant Number JP17K05607. HN acknowledges the visiting profes-
sorship program of University of Lyon 1.

**Author contributions**

H.N. and O.P.-L. conceived the research. H.N. performed the simulation and analyzed data. H.N. wrote the
manuscript, which was edited by O.P.-L.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-87073-6.

Correspondence and requests for materials should be addressed to H.N.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021