Correlations in simulated model bilayers

Molecular dynamics simulations of a model bilayer in a solvent are used to determine the structure factors S(q) and other correlations at vanishing lateral tension of the system and at positive and negative tensions as well. The role played by curvature in the absence of surface tension, metastability, and phase transitions, and the dependence on the specific area, are discussed. Possible height-height correlation functions S(q) are examined and tested.     

Balance of enthalpy and entropy in depletion forces

Solutes added to solutions often dramatically impact molecular processes ranging from the suspension or precipitation of colloids to biomolecular associations and protein folding. Here we revisit the origins of the effective attractive interactions that emerge between and within macromolecules immersed in solutions containing cosolutes that are preferentially excluded from the macromolecular interfaces. Until recently, these depletion forces were considered to be entropic in nature, resulting primarily from the tendency to increase the space available to the cosolute. However, recent experimental evidence indicates the existence of additional, energetically-dominated mechanisms. In this review we follow the emerging characteristics of these different mechanisms. By compiling a set of available thermodynamic data for processes ranging from protein folding to protein–protein interactions, we show that excluded cosolutes can act through two distinct mechanisms that correlate to a large extent with their molecular properties. For many polymers at low to moderate concentrations the steric interactions and molecular crowding effects dominate, and the mechanism is entropic. To contrast, for many small excluded solutes, such as naturally occurring osmolytes, the mechanism is dominated by favorable enthalpy, whereas the entropic contribution is typically unfavorable. We review the available models for these thermodynamic mechanisms, and comment on the need for new models that would be able to explain the full range of observed depletion forces.     

Hydration forces between bilayers in the presence of dissolved or surface-linked sugars

We analyse the experimental evidence of the hydration force near phospholipid bilayers when the “solvent” is a solution of carbohydrates. Two cases must be clearly distinguished: when sugar is dissolved, depletion causes a supplementary attractive force, while in the case of sugar linked to the lipid the contact pressure increases by orders of magnitude. Attractive interaction inferred between bilayers is sometimes derived from indirect evidence, i.e. scattering, attraction between layers adsorbed, shape of phase boundary limits, and without the simultaneous determination of the osmotic compressibility. Generally, water molecules in the first hydration shell of sugar compete with water molecules bound (by more than one kT in free energy) to lipid head-groups. A general result is that the decay length of any repulsive effect remains close to 0.2 nm, even in concentrated sugar solutions. A tentative general explanation of this experimental fact is given together with consequences, such as the possibility of several types of critical points appearing in bilayer stacks. Decay length as well as effective contact pressure is considered with respect to carbohydrate activity.     

Hydration forces between surfactant bilayers: An equilibrium binding description

Water vapour pressure measurements at several temperatures are reported for non-ionic surfactant/water liquid crystals. The results are fitted to a simple model involving an equilibrium between head-group-bound and free water. This model, which uses parameters previously derived from NMR data, allows “hydration” forces to be estimated without invoking layers of “structured” water at an interface.     

Effects of long-range electrostatic forces on simulated protein folding kinetics

Molecular dynamics simulations are a useful tool for characterizing protein folding pathways. There are several methods of treating electrostatic forces in these simulations with varying degrees of physical fidelity and computational efficiency. In this article, we compare the reaction field (RF) algorithm, particle-mesh Ewald (PME), and tapered cutoffs with increasing cutoff radii to address the impact of the electrostatics method employed on the folding kinetics. We quantitatively compare different methods by a correlation of quantitative measures of protein folding kinetics. The results of these comparisons show that for protein folding kinetics, the RF algorithm can quantitatively reproduce the kinetics of the more costly PME algorithm. These results not only assist the selection of appropriate algorithms for future simulations, but also give insight on the role that long-range electrostatic forces have in protein folding.     

Interaction of oleic acid with dipalmitoylphosphatidylcholine (DPPC) bilayers simulated by molecular dynamics

The fatty acid oleic acid (OA) is known to modulate the structure of membranes, which forms the basis for a number of its important applications including its use as a therapeutic supplement to reduce the risk of cardiovascular disease, in molecule delivery systems such as liposomes, and as a skin permeability enhancer. While a number of studies have investigated the effect of OA on lipid membranes, our understanding of its mechanisms of action at the molecular level remains rudimentary. We have carried out molecular dynamics simulations using coarse-grained models to investigate the interactions of OA at a range of concentrations with a dipalmitoylphosphatidylcholine (DPPC) bilayer in the liquid-crystalline phase. We have also investigated the relative permeability of the bilayers to model hydrophilic and hydrophobic penetrants by means of chemical potential calculations. The results indicate that OA is able to disperse homogeneously into the bilayer at all concentrations without much perturbation. OA appears to slightly weaken the lateral forces between lipid headgroups, and as the concentration of OA increases this manifests itself as a slight decrease in the area compressibility modulus and a minor increase in the diffusion rate of the OA molecules. While the chemical potential profiles showed little or no variation as a function of OA concentration, the frequency of water permeation events was found to double, indicating some OA-induced permeability enhancement. The study suggests that physiological effects of OA are probably more subtle rather than via gross perturbation of the structure, or that its significant effects are restricted to more condensed membrane structures such as the gel phase.     

Balance of hydrophobic and electrostatic forces in the pH response of weak polyelectrolyte capsules

A detailed study of the role of solution pH and ionic strength on the swelling behavior of capsules composed of the weak polyelectrolytes poly(4-vinylpyridine) (P4VP) and poly(methacrylic acid) (PMA) with different numbers of layers was carried out. The polyelectrolyte layers were assembled onto silicon oxide particles and multilayer formation was followed by zeta-potential measurements. Hollow capsules were investigated by scanning electron microscopy and atomic force microscopy. The pH-dependent behavior of P4VP/PMA capsules was probed in aqueous media using confocal laser scanning microscopy. All systems exhibited a pronounced swelling at the edges of stability, at pHs of 2 and 8.1. The swelling degree increased when more polymer material was adsorbed. The swollen state can be attributed to uncompensated positive and negative charges within the multilayers, and it is stabilized by counteracting hydrophobic interactions. The swelling was related to the electrostatic interactions by infrared spectroscopy and zeta-potential measurements. The stability of the capsules as well as the swelling degree at a given pH could be tuned, when the ionic strength of the medium was altered.     

Polymerisation in supported bilayers

The intercalation of the twin-tailed amphiphilic dioctadecyldimethylammonium (DODA) ions in smectite clay minerals provides a well-defined supported bilayer system. Polymerisation of styrene in these robust model bilayers allows one to focus on the effect of the constrained medium on the polymerisation process and the polymer. Small-angle X-ray scattering analysis and differential scanning calorimetry data indicated that the polymerisation of styrene in DODA–montmorillonite leads to phase separation between polystyrene and the surfactant matrix. Received: 14 April 2000/Accepted: 19 June 2000     

Flip-flop in adsorbed bilayers

There is very rapid exchange between molecules between an adsorbed (supported) bilayer and solution. Complete exchange of a bilayer occurs within seconds; the actual exchange process may be faster but is limited by the resolution of our technique. The exchange process was monitored by infrared spectroscopy, which can independently monitor deuterium-labeled and unlabeled surfactants, and surface selectivity was achieved by using an internal reflection configuration. Total exchange requires exchange of both the inner and outer layer. We demonstrate that the desorption of a monolayer into water is extremely slow, yet exchange of the inner layer into surfactant solution is extremely fast when a bilayer is present. This shows that flip-flop for an adsorbed bilayer is very fast. Rapid flip-flop at interfaces is in sharp contrast to the extremely slow rate of flip-flop in vesicles. We hypothesize that the rapid rate of flip-flop in adsorbed layers is a consequence of surface-induced defects; in free vesicles, the defects are closed by the film pressure.     

Hydrophobic forces, electrostatic steering, and acid-base bridging between atomically smooth self-assembled monolayers and end-functionalized PEGolated lipid bilayers

A molecular level understanding of interaction forces and dynamics between asymmetric apposing surfaces (including end-functionalized polymers) in water plays a key role in the utilization of molecular structures for smart and functional surfaces in biological, medical, and materials applications. To quantify interaction forces and binding dynamics between asymmetric apposing surfaces in terms of their chemical structure and molecular design we developed a novel surface forces apparatus experiment, using self-assembled monolayers (SAMs) on atomically smooth gold substrates. Varying the SAM head group functionality allowed us to quantitatively identify, rationalize, and therefore control which interaction forces dominated between the SAM surfaces and a surface coated with short-chain, amine end-functionalized polyethylene glycol (PEG) polymers extending from a lipid bilayer. Three different SAM-terminations were chosen for this study: (a) carboxylic acid, (b) alcohol, and (c) methyl head group terminations. These three functionalities allowed for the quantification of (a) specific acid-base bindings, (b) steric effects of PEG chains, and (c) adhesion of hydrophobic segments of the polymer backbone, all as a function of the solution pH. The pH-dependent acid-base binding appears to be a specific and charge mediated hydrogen bonding interaction between oppositely charged carboxylic acid and amine functionalities, at pH values above the acid pK(A) and below the amine pK(A). The long-range electrostatic "steering" of acid and base pairs leads to remarkably rapid binding formation and high binding probability of this specific binding even at distances close to full extension of the PEG tethers, a result which has potentially important implications for protein folding processes and enzymatic catalysis.     

Dimetallic functionalities in liposome bilayers

A dicopper(II) complex can be covalently linked to palmitate/palmitoyl-oleoyl-phosphatidylcholine (PA/POPC) liposomes using the following one-pot strategy: preformed [Cu₂(bpbp)(PA)](ClO₄)₂ (bpbp^?  = 2,6-bis((N,N′-bis(2-picolyl)amino)methyl)-4-tertbutylphenolato) was incorporated into POPC liposomes with a loading of up to 10 mol%. Despite its shape and charge, the decoration of PA/POPC liposomes with {Cu₂(bpbp)}³⁺ did not disrupt the liposome structure; however, the mean liposome diameter increased from about 130 nm (0 mol% dicopper complex) to about 150 nm (10 mol% dicopper complex). Single-crystal X-ray structures furnish ‘snapshots’ of the pH-dependent solution state derivatives of {Cu₂(bpbp)}³⁺, and model the structure of the [Cu₂(bpbp)(PA)]²⁺ head group at the surface of the liposomes. An impressive plasticity in the intramolecular non-bonded Cu….Cu distance for these ions, ranging from 3 to 4 Å, in [Cu₂(bpbp)OH]²⁺, [Cu₂(bpbp)(OAc)(H₂O)]²⁺ and [Cu₂(bpbp)(H₂O)₂]³⁺ allows for their utility as labile reagents in water. Remarkably, the flexible dicopper site is selective for a single carboxylate ligand, so that [Cu₂(bpbp)(PA)]²⁺ is favoured even in the presence of other chemically similar oxoanions, such as CO₃^2 ? , HCO₃^? , NO₃^? , ClO₄^? , ReO₄^?  and CF₃SO₃^? .     

Lateral diffusion in sphingomyelin bilayers

Sphingomyelin (SM) is an important lipid of eukaryotic cellular membranes and neuronal tissues. We studied lateral diffusion in macroscopically oriented bilayers of synthetic palmitoylsphingomyelin (PSM) and natural sphingomyelins of egg yolk (eSM), bovine brain (bSM) and bovine milk (mSM) by pulsed field gradient NMR (PFG NMR) in the temperature range 45-60 °C. We found that the mean values of lateral diffusion coefficients (LDCs) of SMs are 1.9-fold lower compared with those of dipalmitoylphosphatidylcholine (DPPC), which is similar in molecular structure. This discrepancy could be explained by the characteristics of intermolecular SM interactions. The LDCs of different SMs differ: egg SM is most similar to PSM; both of them have a 10% higher LDC value compared with the other two natural SMs. Besides, all natural SMs show a complicated form of the spin-echo diffusion decay (DD), which is an indicator of a distribution of LDC values in bilayers. This peculiarity is explained by the broad distributions of hydrocarbon chain lengths of the natural SMs studied here, especially mSM and bSM. We confirmed the relationship between chain length and LDC in the bilayers by computer analysis of a set of (1)H NMR spectra obtained by scanning the value of the pulsed field gradient. There is a correlation between lower LDC values and SM molecules with longer acyl chains. The most probable mechanisms by which long-chain SM molecules decrease their lateral diffusion relative to the average value are protrusion into the other side of the bilayer or lateral separation into areas that diverge with their LDCs.     

Supported lipopolysaccharide bilayers

In this report, the formation of supported lipopolysaccharide bilayers (LPS-SLBs) is studied with extracted native and glycoengineered LPS from Escherichia coli ( E. coli ) and Salmonella enterica sv typhimurium ( S. typhimurium ) to assemble a platform that allows measurement of LPS membrane structure and the detection of membrane tethered saccharide-protein interactions. We present quartz crystal microbalance with dissipation monitoring (QCM-D) and fluorescence recovery after photobleaching (FRAP) characterization of LPS-SLBs with different LPS species, having, for example, different molecular weights, that show successful formation of SLBs through vesicle fusion on SiO(2) surfaces with LPS fractions up to 50 wt %. The thickness of the LPS bilayers were investigated with AFM force-distance measurements which showed only a slight thickness increase compared to pure POPC SLBs. The E. coli LPS were chosen to study the saccharide-protein interaction between the Htype II glycan epitope and the Ralstonia solanacearum lectin (RSL). RSL specifically recognizes fucose sugars, which are present in the used Htype II glycan epitope and absent in the epitopes LPS1 and EY2. We show via fluorescence microscopy that the specific, but weak and multivalent interaction can be detected and discriminated on the LPS-SLB platform.     

Flexoelectricity of lipid bilayers

Abstract

The direct flexoeffect in single lipid bilayers in the form of black lipid membranes has been investigated experimentally by the oscillating pressure technique in the regime of voltage measurement. Black lipid membranes of various composition have been studied in order to check the effect of lipid surface charge on the curvature-electric response and its frequency dependence; these include egg yolk lecithin (low negative charge); egg yolk lecithin plus phosphatidyl serine (high negative charge); egg yolk lecithin with surface adsorbed ions of uranyl acetate (high positive charge). An increase of the response has been found by increasing the surface charge and a reversal of the sign of the flexoelectric coefficient from positive to negative has been obtained by changing the sign of the surface charge from negative to positive. These results underline the leading role of the contribution of the surface charge to the flexoelectricity of lyotropics. Their theoretical interpretation provides further insight into the molecular mechanism of this phenomenon.     

Self-assembling bilayers of palladiumthiolates in organic media

Alkylthiolates of palladium forming a homologous series (butyl to octadecyl) have been prepared and characterized using X-ray diffraction and STM. The thiolates adopt an unusual bilayered lamellar structure, whose thickness is governed by the length of the alkyl chain. These mesophases melt in the temperature range, 60° to 100°C, with the melting point increasing linearly with the thiol chain length. There is evidence to suggest that the alkyl chains are orientationally disordered especially prior to melting.     

Selective sterol-phospholipid associations in fluid bilayers

The nearest-neighbor preferences of three exchangeable lipid monomers (two phospholipids that differ in length, A and B, and a derivative of cholesterol, C) have been quantified in the fluid bilayer state by use of the nearest-neighbor recognition method (Davidson, S. K. M.; Regen, S. L. Chem. Rev. 1997, 97, 1269). Thus, an analysis of the equilibrium dimer distributions has shown that (i) the sterol favors both phospholipids as nearest neighbors relative to other sterol molecules, (ii) that this recognition is selective (i.e., the sterol favors the longer phospholipid as a nearest neighbor over the shorter one, especially when the sterol concentrations in the bilayer is high (e.g., 40 mol %), and (iii) the phospholipids, themselves, are unable to recognize each other. Taken together, these findings indicate that the probable mechanism by which cholesterol induces homoassociation of A and B in analogous bilayers is one in which the sterol "pulls" two or more of the longer phospholipid monomers (B) out of a "sea" of randomly mixed A and B. These findings also lend support for the notion of cholesterol-phospholipid complexation in fluid bilayers. The biological implications of these findings are briefly discussed.     

Subgel studies of dimyristoylphosphatidylcholine bilayers

It is known that when bilayers of some saturated phosphatidylcholines are stored for 3 or more days at approximately 0 degrees C, a lamellar subgel (Lc) phase is detected at temperatures below the pretransition by differential scanning calorimetry (DSC). However, the subgel (Lc) phase and the corresponding subtransition (Lc--> Lbeta') for dimyristoylphosphatidylcholine (DMPC) has not been clearly characterized. In this study, using the temperature jump protocol first developed by Tristram-Nagle et al. for the dipalmitoylphosphatidylcholine (DPPC) system, new and accurate data characterizing the subgel formation and subtransition of DMPC were obtained through DSC and fluorescence spectroscopy with 1,6-diphenyl-1,3,5-hexatriene (DPH). It was discovered that the formation of the DMPC subgel phase requires incubation at temperatures of -5 degrees C or lower for 2 h or more. Kinetics of the subgel formation indicate that it is a very complex process and demonstrates that the planar gel phase is merely metastable below the subtransition, and not the thermodynamically stable phase. The subgel growth of DMPC is proven to be the dehydration of the headgroup region, and the subtransition is a process in which poorly hydrated DMPC becomes hydrated.