Molecular modeling of proteinlike inclusions in lipid bilayers: lipid-mediated interactions

Self-consistent field modeling of mesomorphic phase changes of monoolein and phospholipids in response to additives

Mapping the topological phase behaviour of lipids in aqueous solution is time consuming and finding the ideal lipid system for a desired application is often a matter of trial and error. Modelling techniques that can accurately predict the mesomorphic phase behaviour of lipid systems are therefore of paramount importance. Here, the self-consistent field theory of Scheutjens and Fleer (SF-SCF) in which a lattice refinement has been implemented, is used to scrutinize how various additives modify the self-assembled phase behaviour of monoolein (MO) and 1,2-dioleoyl-phosphatidylcholine (DOPC) lipids in water. The mesomorphic behaviour is inferred from trends in the mechanical properties of equilibrium lipid bilayers with increasing additive content. More specifically, we focus on the Helfrich parameters, that is, the mean and Gaussian bending rigidities (κ and [small kappa, Greek, macron], respectively) supplemented with the spontaneous curvature of the monolayer (Jm0). We use previously established interaction parameters that position the unperturbed DOPC system in the lamellar Lα phase ([small kappa, Greek, macron] < 0, κ > 0 and Jm0 ≈ 0). Similar interaction parameters position the MO system firmly in a bicontinuous cubic phase ([small kappa, Greek, macron] > 0). In line with experimental data, a mixture of MO and DOPC tends to be in one of these two phases, depending on the mixing ratio. Moreover we find good correlations between predicted trends and experimental data concerning the phase changes of MO in response to a wide range of additives. These correlations give credibility to the use of SF-SCF modelling as a valuable tool to quickly explore the mesomorphic phase space of (phospho)lipid bilayer systems including additives.

Structural and mechanical parameters of lipid bilayer membranes using a lattice refined self-consistent field theory

The self-consistent field theory of Scheutjens and Fleer is implemented on a grid with (lattice) sites that are smaller than the segment size. In this quasi lattice-free implementation we consider united atom-like molecular models and study bilayer self-assembly of phospholipids in a selective solvent (water). We find structural as well as mechanical parameters for these bilayers. The mean (κ) and Gaussian ([small kappa, Greek, macron]) bending moduli, as well as the spontaneous curvature of the monolayer (Jm0), are computed for the first time following a grand canonical ensemble route. Results are in line with previous estimates for mechanical parameters that at the time could not be made following this correct route. This proves that the mean bending modulus is only a very weak function of the membrane tension. We performed a systematic study on the effects of model parameter variations. The mean bending modulus generally grows with increasing bilayer thickness. As expected Jm0 and [small kappa, Greek, macron] behave oppositely with respect to each other and for classical phospholipids assumes values near zero. As an example, an increase in the lipophilic to hydrophilic ratio in the lipids, may cause the Gaussian bending rigidity to switch sign from negative to positive, while - not necessarily at the same point - the spontaneous curvature of the monolayer may switch sign from positive to negative. Together with other investigated trends, these results point to mechanisms of how topological phase transitions of the lipid bilayer membranes may be regulated in the biological context, which correlates with known lipid phase behaviour.

The physics of microemulsions extracted from modeling balanced tensionless surfactant-loaded liquid-liquid interfaces

Microemulsions are explored using the self-consistent field approach. We consider a balanced model that features two solvents of similar size and a symmetric surfactant. Interaction parameter χ and surfactant concentration φ_s ^b complement the model definition. The phase diagram in χ-φ_s ^b coordinates is known to feature two lines of critical points, the Scott and Leibler lines. Only upon imposing a finite distance between the interfaces, we observe that the Scott line meets the Leibler line. We refer to this as a Lifshitz point (LP) for real systems. We add regions that are relevant for microemulsions to this phase diagram by considering the saturation line, which connects (χ, φ_s ^b)-points for which the interface becomes tensionless. Crossing this line implies a first-order phase transition as internal interfaces develop, characteristic for one-phase microemulsions. The saturation line ends at the so-called microemulsion point (MP). The MP is shown to connect with the LP by a line of MP-like critical points, found by searching for a "MP" while the distance between interfaces is fixed. A pair of binodal lines that envelop the three-phase (Winsor III) microemulsion region is shown to connect to the MP. The cohesiveness of the middle phase in Winsor III is related to non-monotonic, inverse DLVO-type interaction curves between the surfactant-loaded tensionless interfaces. The mean and Gaussian bending modulus, relevant for the shape fluctuations and the topology of interfaces, respectively, are evaluated along the saturation line. Near the MP, both rigidities are positive and vanish in a power-law fashion with coefficient unity at the MP. Overseeing these results proves that the MP has a pivoting role in the physics of microemulsions.

Elastic properties of symmetric liquid-liquid interfaces

The mean (κ) and Gaussian (κ[over ¯]) bending rigidities of liquid-liquid interfaces, of importance for shape fluctuations and topology of interfaces, respectively, are not yet established: Even their signs are debated. Using the Scheutjens Fleer variant of the self-consistent field theory, we implemented a model for a symmetric L-L interface and obtained high-precision (mean-field) results in the grand-canonical (μ,V,T) ensemble. We report positive values for both moduli when the system is close to critical where the rigidities show the same scaling behavior as the interfacial tension γ. At strong segregation, when the interfacial width becomes of the order of the segment size, κ[over ¯] turns negative. The length scale λ≡sqrt[κ/γ] remains of the order of segment size for all strengths of interaction; yet the 1/sqrt[N] chain length correction reduces λ significantly when the chain length N is small.

Sign Switch of Gaussian Bending Modulus for Microemulsions: A Self-Consistent Field Analysis Exploring Scale Invariant Curvature Energies

Bending rigidities of tensionless balanced liquid-liquid interfaces as occurring in microemulsions are predicted using self-consistent field theory for molecularly inhomogeneous systems. Considering geometries with scale invariant curvature energies gives unambiguous bending rigidities for systems with fixed chemical potentials: the minimal surface Im3m cubic phase is used to find the Gaussian bending rigidity κ[over ¯], and a torus with Willmore energy W=2π^{2} allows for direct evaluation of the mean bending modulus κ. Consistent with this, the spherical droplet gives access to 2κ+κ[over ¯]. We observe that κ[over ¯] tends to be negative for strong segregation and positive for weak segregation, a finding which is instrumental for understanding phase transitions from a lamellar to a spongelike microemulsion. Invariably, κ remains positive and increases with increasing strength of segregation.

On the edge energy of lipid membranes and the thermodynamic stability of pores

To perform its barrier function, the lipid bilayer membrane requires a robust resistance against pore formation. Using a self-consistent field (SCF) theory and a molecularly detailed model for membranes composed of charged or zwitterionic lipids, it is possible to predict structural, mechanical, and thermodynamical parameters for relevant lipid bilayer membranes. We argue that the edge energy in membranes is a function of the spontaneous lipid monolayer curvature, the mean bending modulus, and the membrane thickness. An analytical Helfrich-like model suggests that most bilayers should have a positive edge energy. This means that there is a natural resistance against pore formation. Edge energies evaluated explicitly in a two-gradient SCF model are consistent with this. Remarkably, the edge energy can become negative for phosphatidylglycerol (e.g., dioleoylphosphoglycerol) bilayers at a sufficiently low ionic strength. Such bilayers become unstable against the formation of pores or the formation of lipid disks. In the weakly curved limit, we study the curvature dependence of the edge energy and evaluate the preferred edge curvature and the edge bending modulus. The latter is always positive, and the former increases with increasing ionic strength. These results point to a small window of ionic strengths for which stable pores can form as too low ionic strengths give rise to lipid disks. Higher order curvature terms are necessary to accurately predict relevant pore sizes in bilayers. The electric double layer overlap across a small pore widens the window of ionic strengths for which pores are stable.

Coverage and disruption of phospholipid membranes by oxide nanoparticles

We studied the interactions of silica and titanium dioxide nanoparticles with phospholipid membranes and show how electrostatics plays an important role. For this, we systematically varied the charge density of both the membranes by changing their lipid composition and the oxide particles by changing the pH. For the silica nanoparticles, results from our recently presented fluorescence vesicle leakage assay are combined with data on particle adsorption onto supported lipid bilayers obtained by optical reflectometry. Because of the strong tendency of the TiO2 nanoparticles to aggregate, the interaction of these particles with the bilayer was studied only in the leakage assay. Self-consistent field (SCF) modeling was applied to interpret the results on a molecular level. At low charge densities of either the silica nanoparticles or the lipid bilayers, no electrostatic barrier to adsorption exists. However, the adsorption rate and adsorbed amounts drop with increasing (negative) charge densities on particles and membranes because of electric double-layer repulsion, which is confirmed by the effect of the ionic strength. SCF calculations show that charged particles change the structure of lipid bilayers by a reorientation of a fraction of the zwitterionic phosphatidylcholine (PC) headgroups. This explains the affinity of the silica particles for pure PC lipid layers, even at relatively high particle charge densities. Particle adsorption does not always lead to the disruption of the membrane integrity, as is clear from a comparison of the leakage and adsorption data for the silica particles. The attraction should be strong enough, and in line with this, we found that for positively charged TiO2 particles vesicle disruption increases with increasing negative charge density on the membranes. Our results may be extrapolated to a broader range of oxide nanoparticles and ultimately may be used for establishing more accurate nanoparticle toxicity assessments and drug-delivery systems.