Dynamic simulations of multicomponent lipid membranes over long length and time scales

Brian A. Camley,Frank L.H. Brown

DOI: https://doi.org/10.1103/PhysRevLett.105.148102

2011-05-25

Abstract:We present a stochastic phase-field model for multicomponent lipid bilayers that explicitly accounts for the quasi-two-dimensional hydrodynamic environment unique to a thin fluid membrane immersed in aqueous solution. Dynamics over a wide range of length scales (from nanometers to microns) for durations up to seconds and longer are readily accessed and provide a direct comparison to fluorescence microscopy measurements in ternary lipid/cholesterol mixtures. Simulations of phase separation kinetics agree with experiment and elucidate the importance of hydrodynamics in the coarsening process.

Soft Condensed Matter

What problem does this paper attempt to address?

This paper attempts to solve the following problems:

1. **Dynamic simulation of multi - component lipid bilayer membranes on long - time and large - space scales**
The author proposes a stochastic phase - field model to describe the dynamic behavior of multi - component lipid bilayer membranes. This model can explicitly consider the unique quasi - two - dimensional (quasi - 2D) hydrodynamic environment of thin fluid membranes immersed in aqueous solutions.

2. **Cross - scale dynamic research**
The model can cover spatial scales from nanometers to micrometers and time scales from seconds to longer times, thus providing a direct basis for comparison with experimental measurements. Specifically, the paper focuses on the results of fluorescence microscopy measurements in ternary lipid/cholesterol mixtures.

3. **Phase separation kinetics and experimental consistency**
The simulation results are consistent with the experiments and reveal the importance of hydrodynamics in the coarsening process. By combining quasi - two - dimensional hydrodynamics and thermal fluctuations, the paper develops a new simulation method that can be explored on experimental - relevant length and time scales.

4. **Verification of theory and experiment**
The paper verifies the consistency between theoretical predictions and experimental observations, such as the relationship of the lipid domain translational diffusion coefficient $D_a=\frac{k_BT}{4\pi\eta_m}F(a/L_{sd})$, and shows the agreement between the relaxation times $\tau_n$ in different modes and theoretical predictions.

5. **Research on complex non - linear dynamics**
Through numerical simulation, the paper shows how to directly study the complex non - linear bilayer membrane dynamics on micrometer - scale and second - level time scales. In particular, for the dynamics of quasi - two - dimensional membrane phase separation, the paper finds that there are significant differences from the pure two - dimensional case and proposes a mechanism for natural explanation.

6. **Extension to more complex biomembrane systems**
The method proposed by the author is not only applicable to simple model membrane systems but also can be extended to more complex biomembrane systems, such as the coupling with the cytoskeleton, the dynamics of membrane proteins, etc.

### Summary of core formulas
- Lipid domain diffusion coefficient:
$$
D_a = \frac{k_BT}{4\pi\eta_m}F\left(\frac{a}{L_{sd}}\right)
$$
where $L_{sd}=\frac{\eta_m}{2\eta_f}$ is the Saffman - Delbrück length scale.

- Relaxation time:
- Stone - McConnell theory:
$$
\tau^{\text{Stone}}_n=\frac{2\pi R_o^2\eta_f}{\sigma(n^2 - 1/4)n^2(n^2 - 1)}
$$
- Mann theory:
$$
\tau^{\text{Mann}}_n=\frac{4\eta_mR_o}{n\sigma(4)}
$$
- Generalized theory:
$$
\tau^{\text{general}}_n=\frac{\eta_mR_o}{\sigma}\left[\int_0^\infty dx\frac{J_n^2(x)}{x^2(x + R_o/L_{sd})}\right]^{-1}
$$

Through these studies, the paper aims to provide a powerful tool for understanding biomembrane dynamics and lay the foundation for future research on complex systems.