Abstract:We have constructed a model for the kinetics of rupture of membranes under tension, applying physical principles relevant to lipid bilayers held together by hydrophobic interactions. The membrane is characterized by the bulk compressibility (for expansion), the thickness of the hydrophobic part of the bilayer, the hydrophobicity and a parameter characterizing the tail rigidity of the lipids. The model is a lattice model which incorporates strain relaxation, and considers the nucleation of pores at constant area, constant temperature, and constant particle number. The particle number is conserved by allowing multiple occupancy of the sites. An equilibrium ``phase diagram'' is constructed as a function of temperature and strain with the total pore surface and distribution as the order parameters. A first order rupture line is found with increasing tension, and a continuous increase in proto-pore concentration with rising temperature till instability. The model explains current results on saturated and unsaturated PC lipid bilayers and thicker artificial bilayers made of diblock copolymers. Pore size distributions are presented for various values of area expansion and temperature, and the fractal dimension of the pore edge is evaluated.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
This paper aims to explore the rupture dynamics of lipid bilayer membranes under tension. Specifically, the author constructs a lattice model to explain the rupture process of the membrane under stretching conditions, especially by considering the following aspects:
1. **Physical properties of the membrane**: The volume compressibility \(K\) of the membrane, the thickness \(2h\) of the hydrophobic part, the hydrophobicity \(\sigma\), and the lipid tail stiffness \(\gamma\).
2. **Hole nucleation mechanism**: In the case of constant area, temperature, and number of particles, study the process of hole nucleation and consider strain relaxation.
3. **Thermodynamic phase diagram**: Construct an equilibrium phase diagram based on temperature and strain to determine the pore - free, sub - critical pore (or native pore), and rupture regions.
4. **Explanation of rupture phenomena**: Explain the experimental results of saturated and unsaturated phosphatidylcholine (PC) lipid bilayer membranes and thicker artificial bilayer membranes composed of diblock copolymers.
### Specific problem description
This research mainly focuses on the following problems:
- **Mechanical mechanism of membrane rupture**: When the membrane is stretched, how to release stress by forming holes.
- **Energy barrier for hole nucleation**: Calculate the energy barrier for hole formation and explore its relationship with the physical properties of the membrane (such as compressibility, hydrophobicity, and tail stiffness).
- **Effect of temperature on membrane stability**: Study how temperature affects the stability of the membrane, especially the hole formation and distribution under different stretching conditions.
- **Different types of lipid bilayer membranes**: Compare the differences in rupture behavior between saturated and unsaturated lipid bilayer membranes and bilayer membranes composed of diblock copolymers.
### Summary of mathematical formulas
The main formulas involved in the paper include:
- Membrane stretching energy:
\[
E_m=\frac{1}{2}Ka_m\left(\frac{\Delta a}{a_0}\right)^2
\]
where \(K\) is the compressibility, \(a_m\) is the total relaxation area, \(a_0\) is the molecular area, and \(\Delta a = a - a_0\) is the extra area per molecule.
- Membrane stretching energy in the presence of holes:
\[
E_m(a_p)=2\lambda\sqrt{\pi a_m}\left(\frac{a_p}{a_m}\right)^{1/2}+\frac{1}{2}Ka_m\left(\frac{\Delta a}{a_0}\right)^2
\]
where \(\lambda\) is the edge tension (effective edge energy per unit length), and \(a_p\) is the hole area.
- Edge tension expression:
\[
\lambda = 2h_t\sigma
\]
where \(h_t\) is the hydrophobic thickness of the lipid tail, and \(\sigma\) is the hydrophobic energy per unit area.
- Surface energy:
\[
E_s=\frac{1}{2}Ka_m\left(\frac{\Delta a}{a_0}\right)^2
\]
- Edge energy:
\[
E_{\text{edge}}=\left(\frac{N}{N - n_h}\right)^{1/2}E_{\text{model, edge}}
\]
These formulas together describe the dynamic process of membrane rupture under stretching conditions and its energy changes.