Electrostatic interactions across a charged lipid bilayer

Alexander J. Wagner,Sylvio May

DOI: https://doi.org/10.48550/arXiv.cond-mat/0607733

2006-07-28

Abstract:We present theoretical work in which the degree of electrostatic coupling across a charged lipid bilayer in aqueous solution is analyzed on the basis of nonlinear Poisson-Boltzmann theory. In particular, we consider the electrostatic interaction of a single, large macroion with the two apposed leaflets of an oppositely charged lipid bilayer where the macroion is allowed to optimize its distance to the membrane. Three regimes are identified: a weak and a high macroion charge regime, separated by a regime of close macroion-membrane contact for intermediate charge densities. The corresponding free energies are used to estimate the degree of electrostatic coupling in a lamellar cationic lipid-DNA complex. That is, we calculate to what extent the one-dimensional DNA arrays in a sandwich-like lipoplex interact across the cationic membranes. We find that, in spite of the low dielectric constant inside a lipid membranes, there can be a significant electrostatic contribution to the experimentally observed cross-bilayer orientational ordering of the DNA arrays. Our approximate analytical model is complemented and supported by numerical calculations of the electrostatic potentials and free energies of the lamellar lipoplex geometry. To this end, we solve the nonlinear Poisson-Boltzmann equation within a unit cell of the lamellar lipoplex using a new lattice Boltzmann method.

Soft Condensed Matter

What problem does this paper attempt to address?

The core problem that this paper attempts to solve is to analyze the electrostatic interactions in charged lipid bilayer membranes in aqueous solutions, especially to study the electrostatic coupling between macromolecular ions and lipid bilayer membranes within the framework of the nonlinear Poisson - Boltzmann theory. Specifically, the article focuses on:

1. **Electrostatic interactions between macromolecular ions and lipid bilayer membranes**: The electrostatic interactions between a single macromolecular ion (macroion) and two adjacent monolayers of a lipid bilayer membrane with opposite charges are studied. In particular, the situation where the macromolecular ion can optimize its distance from the membrane is considered in the paper.

2. **Three electrostatic interaction mechanisms**:
- **Weak - charge macromolecular ion mechanism**: When the charge of the macromolecular ion is low, the distance between it and the membrane is large.
- **High - charge macromolecular ion mechanism**: When the charge of the macromolecular ion is high, it will be in close contact with the membrane.
- **Intermediate mechanism**: For macromolecular ions with medium charge density, there is an intermediate mechanism where the macromolecular ion is in close contact with the membrane.

3. **Free energy calculation**: By calculating the free energy under different mechanisms, the degree of electrostatic coupling is evaluated and applied to the electrostatic interactions between DNA arrays in lamellar cationic lipid - DNA complexes (lipoplex structures).

4. **Explanation of experimental phenomena**: The paper points out that although the dielectric constant inside the lipid membrane is low, there may still be a significant electrostatic contribution, which helps to explain the orientation order of transmembrane DNA arrays observed in experiments.

5. **Numerical verification**: To support the theoretical model, the author also uses a new lattice Boltzmann method to numerically solve the two - dimensional Poisson - Boltzmann equation and verifies the accuracy of the theoretical prediction.

### Formula summary

- **Bjerrum length** \( l_B=\frac{e^2}{4\pi\epsilon_0\epsilon_W k_B T} \)
- **Debye shielding length** \( l_D = \frac{1}{\kappa}, \quad \kappa^2=\frac{2 e^2 n_0}{\epsilon_0\epsilon_W k_B T} \)
- **Coupling parameter** \( H=\frac{\epsilon_L l_D}{\epsilon_W d} \)
- **Surface potential** \( \Psi_E=-2 \text{arsinh}(p_0 \bar{\phi}_E), \quad \Psi_I=-2 \text{arsinh}(p_0 \bar{\phi}_I) \)
- **Effective composition** \( \bar{\phi}_E=\phi_E+\frac{H \Delta \Psi}{2 p_0}, \quad \bar{\phi}_I=\phi_I-\frac{H \Delta \Psi}{2 p_0} \)
- **Transcendental equation of potential difference** \( \frac{\Delta \Psi}{2}=\text{arsinh}\left(\frac{p_0 \phi_I-\frac{H}{2} \Delta \Psi}{\sqrt{1+(p_0 \bar{\phi}_I)^2}}\right)-\text{arsinh}\left(\frac{p_0 \phi_E+\frac{H}{2} \Delta \Psi}{\sqrt{1+(p_0 \bar{\phi}_E)^2}}\right) \)

These formulas and theoretical analyses help to understand and explain the complex electrostatic interaction mechanisms in charged lipid bilayer membranes.