Surface instability of confined elastic bilayers: Theory and simulations

The surface of a soft elastic film becomes unstable and forms a self-organized undulating pattern because of adhesive interactions when it comes in contact proximity with a rigid surface. For a single film, the pattern length scale lambda, which is governed by the minimization of the elastic stored energy, gives lambda similar to 3h, where h is the film thickness. Based on a linear stability analysis and simulations of adhesion and debonding, we consider the contact instability of an elastic bilayer, which provides greater flexibility in the morphological control of interfacial instability. Unlike the case of a single film, the morphology of the contact instability patterns, debonding distance, and debonding force in a bilayer can be controlled in a nonlinear way by varying the thicknesses and shear moduli of the films. Interestingly, the pattern wavelength in a bilayer can be greatly increased or decreased compared to a single film when the adhesive contact is formed by the stiffer or the softer of the two films, respectively. In particular, lambda as small as 0.5h can be obtained. This indicates a new strategy for pattern miniaturization in elastic contact lithography.

Observation of reduced activation energy and the possible existence of decoupled pancake vortices in superconductor/ferromagnet bilayers

Investigations of different superconducting (S)/ferromagnetic (F) heterostructures grown by pulsed laser deposition reveal that the activation energy (U) for the vortex motion in a high T-c superconductor is reduced remarkably by the presence of F layers. The U exhibits a logarithmic dependence on the applied magnetic field in the S/F bilayers suggesting the existence of decoupled two-dimensional (2D) pancake vortices. This result is discussed in terms of the reduction in the effective S layer thickness and the weakening of the S coherence length due to the presence of F layers. In addition, the U and the superconducting T-c in YBa2Cu3O7-delta/La0.5Sr0.5CoO3 bilayers are observed to be much lower than in the YBa2Cu3O7-delta/La0.7Sr0.3MnO3 ones. This in turn suggests that the degree of spin polarization of the F layer might not play a crucial role for the suppression of superconductivity due to a spin polarized induced pair-breaking effect in S/F bilayers.

Investigations on the melting and bending modulus of polymer grafted bilayers using dissipative particle dynamics

Understanding the influence of polymer grafted bilayers on the physicomechanical properties of lipid membranes is important while developing liposomal based drug delivery systems. The melting characteristics and bending moduli of polymer grafted bilayers are investigated using dissipative particle dynamics simulations as a function of the amount of grafted polymer and lipid tail length. Simulations are carried out using a modified Andersen barostat, whereby the membrane is maintained in a tensionless state. For lipids made up of four to six tail beads, the transition from the low temperature L-beta phase to the L-alpha phase is lowered only above a grafting fraction of G(f)=0.12 for polymers made up of 20 beads. Below G(f)=0.12 small changes are observed only for the HT4 bilayer. The bending modulus of the bilayers is obtained as a function of G(f) from a Fourier analysis of the height fluctuations. Using the theory developed by Marsh Biochim. Biophys. Acta 1615, 33 (2003)] for polymer grafted membranes, the contributions to the bending modulus due to changes arising from the grafted polymer and bilayer thinning are partitioned. The contributions to the changes in kappa from bilayer thinning were found to lie within 11% for the lipids with four to six tail beads, increasing to 15% for the lipids containing nine tail beads. The changes in the area stretch modulus were also assessed and were found to have a small influence on the overall contribution from membrane thinning. The increase in the area per head group of the lipids was found to be consistent with the scalings predicted by self-consistent mean field results. (C) 2010 American Institute of Physics.

Microscopic origin of the chirality driven morphologies of the amphiphilic monolayers and bilayers

It is widely known that the compressed monolayers and bilayers of chiral lipids or fatty acids form helical morphologies, while the corresponding racemic modification gives only flat platelets without twist. No molecular explanation of this phenomenon is yet available, although subtle interactions at the chiral centers have often been proposed as the driving force behind the morphology of the aggregate to form a particular shape. In the present study, the morphologies of the chiral amphiphilic assemblies have been predicted on the basis of an effective pair potential between the molecules, which depends on the relative sizes of the groups attached to the chiral centers, the orientation of the amphiphilic molecules and also on the distance between them. It is shown that fur a pair of same kind of enantiomers, the minimum energy conformation favours a twist angle between them. This twist between the neighbouring molecules gives rise to the helicity of the aggregate. The present theory also shows from the molecular considerations that for a pair of mirror-image isomers (i.e. the racemic modification) the minimum energy conformation corresponds to the zero angle between the molecules, thus giving rise to flat platelets as observed in experiments. Another fascinating aspect of such chirality driven helical structures is that the sense (or the handedness) of the helix is highly specific about the chirality of the monomer concerned. The molecular theory shows, for the first time, that the sense of the helical structures in many cases is determined by the sizes of the groups attached to the chiral centers and the effective potential between them. The predicted senses of the helical structures are in complete agreement with the experimental results.

Prediction of the senses of helical amphiphilic assemblies from effective intermolecular pair potential: Studies on chiral monolayers and bilayers

It is well-known that the senses (or the handedness) of the helical assemblies formed from compressed monolayers and bilayers of chiral amphiphiles are highly specific about the chirality of the monomers concerned. We present here a molecular approach that can successfully predict the senses of such helical morphologies. The present approach is based on a reduced tractable description in terms of an effective pair potential (EPP) which depends on the distance of separation and the relative orientations of the two amphiphiles. This approach explicitly considers the pairwise intermolecular interactions between the groups attached to the chiral centers of the two neighboring amphiphiles. It is found that for a pair of the same kind of enantiomers the minimum energy configuration favors a twist angle between molecules and that this twist from neighbor to neighbor gives rise to the helicity of the aggregate. From the known twist angles at the minimum energy configuration the successive arrangement of an array of molecules can be predicted. Therefore, the sense of the helicity can be predicted from the molecular interactions. The predicted senses of the helical structures are in complete agreement with all known experimental results.

Phase transitions in the anchored organic bilayers of long-chain alkylammonium lead iodides (CnH2n+1NH3)(2)PbI4; n=12, 16, 18

The single perovskite slab alkylammonium lead iodides (CnH2n+1NH3)(2)PbI4, n = 12, 16, 18, display two phase transitions, just above room temperature, associated with changes in the alkylammonium chains. We have followed these two phase transitions using scanning calorimetry, X-ray powder diffraction, and IR and Raman spectroscopies. We find the first phase transition to be associated with symmetry changes arising from a dynamic rotational disordering of the ammonium headgroup of the chain whereas the second transition, the melting of the chains in two dimensions, is characterized by an increased conformational disorder of the methylene units of the alkyl chains. We examine these phase transitions in light of the interesting optical properties of these materials, as well as the relevance of these systems as models for phase transitions in lipid bilayers.

Estimation of activation energy for electroporation and pore growth rate in liquid crystalline and gel phases of lipid bilayers using molecular dynamics simulations

Molecular dynamics simulations of electroporation in POPC and DPPC lipid bilayers have been carried out at different temperatures ranging from 230 K to 350 K for varying electric fields. The dynamics of pore formation, including threshold field, pore initiation time, pore growth rate, and pore closure rate after the field is switched off, was studied in both the gel and liquid crystalline (L-alpha) phases of the bilayers. Using an Arrhenius model of pore initiation kinetics, the activation energy for pore opening was estimated to be 25.6 kJ mol(-1) and 32.6 kJ mol(-1) in the L-alpha phase of POPC and DPPC lipids respectively at a field strength of 0.32 V nm(-1). The activation energy decreases to 24.2 kJ mol(-1) and 23.7 kJ mol(-1) respectively at a higher field strength of 1.1 V nm(-1). At temperatures below the melting point, the activation energy in the gel phase of POPC and DPPC increases to 28.8 kJ mol(-1) and 34.4 kJ mol(-1) respectively at the same field of 1.1 V nm(-1). The pore closing time was found to be higher in the gel than in the L-alpha phase. The pore growth rate increases linearly with temperature and quadratically with field, consistent with viscosity limited growth models.

Differential dynamics of the serotonin(1A) receptor in membrane bilayers of varying cholesterol content revealed by all atom molecular dynamics simulation

The serotonin(1A) receptor belongs to the superfamily of G protein-coupled receptors (GPCRs) and is a potential drug target in neuropsychiatric disorders. The receptor has been shown to require membrane cholesterol for its organization, dynamics and function. Although recent work suggests a close interaction of cholesterol with the receptor, the structural integrity of the serotonin(1A) receptor in the presence of cholesterol has not been explored. In this work, we have carried out all atom molecular dynamics simulations, totaling to 3s, to analyze the effect of cholesterol on the structure and dynamics of the serotonin(1A) receptor. Our results show that the presence of physiologically relevant concentration of membrane cholesterol alters conformational dynamics of the serotonin(1A) receptor and, on an average lowers conformational fluctuations. Our results show that, in general, transmembrane helix VII is most affected by the absence of membrane cholesterol. These results are in overall agreement with experimental data showing enhancement of GPCR stability in the presence of membrane cholesterol. Our results constitute a molecular level understanding of GPCR-cholesterol interaction, and represent an important step in our overall understanding of GPCR function in health and disease.

Molecular Origin of the Intrinsic Bending Force for Helical Morphology Observed in Chiral Amphiphilic Assemblies: Concentration and Size Dependence

The formation of the helical morphology in monolayers and bilayers of chiral amphiphilic assemblies is believed to be driven at least partly by the interactions at the chiral centers of the amphiphiles. However, a detailed microscopic understanding of these interactions and their relation with the helix formation is still not clear. In this article a study of the molecular origin of the chirality-driven helix formation is presented by calculating, for the first time, the effective pair potential between a pair of chiral molecules. This effective potential depends on the relative sizes of the groups attached to the two chiral centers, on the orientation of the amphiphile molecules, and also on the distance between them. We find that for the mirror-image isomers (in the racemic modification) the minimum energy conformation is a nearly parallel alignment of the molecules. On the other hand, the same for a pair of molecules of one kind of enantiomer favors a tilt angle between them, thus leading to the formation of a helical morphology of the aggregate. The tilt angle is determined by the size of the groups attached to the chiral centers of the pair of molecules considered and in many cases predicted it to be close to 45 degrees. The present study, therefore, provides a molecular origin of the intrinsic bending force, suggested by Helfrich (J. Chem. Phys. 1986, 85, 1085-1087), to be responsible for the formation of helical structure. This effective potential may explain many of the existing experimental results, such as the size and the concentration dependence of the formation of helical morphology. It is further found that the elastic forces can significantly modify the pitch predicted by the chiral interactions alone and that the modified real pitch is close to the experimentally observed value. The present study is expected to provide a starting point for future microscopic studies.

Chiral molecular self-assembly of phospholipid tubules: A circular dichroism study

We report on spectroscopic studies of the chiral structure in phospholipid tubules formed in mixtures of alcohol and water. Synthetic phospholipids containing diacetylenic moieties in the acyl chains self-assemble into hollow, cylindrical tubules in appropriate conditions. Circular dichroism provides a direct measure of chirality of the molecular structure. We find that the CD spectra of tubules formed in mixtures of alcohol and water depends strongly on the alcohol used and the lipid concentration. The relative spectral intensity of different circular dichroism bands correlates with the number of bilayers observed using microscopy. The results provide experimental evidence that tubule formation is based on chiral packing of the lipid molecules and that interbilayer interactions are important to the tubule structure

X-Ray Studies On The Bilayer Structure Of Trypsin-Treated Rat-Brain Myelin

Trypsin-treated rat brain myelin was subjected to biochemical and X-ray studies. Untreated myelin gave rise to a pattern of three rings with a fundamental repeat period of 155 Angstrom consisting of two bilayers per repeat period, whereas myelin treated with trypsin showed a fundamental repeat period of 75 Angstrom with one bilayer per repeat period. The integrated raw intensity of the h=4 reflection with respect to the h=2 reflection is 0.38 for untreated myelin. The corresponding value reduced to 0.23, 0.18, 0.17 for myelin treated with 5, 10, 40 units of trypsin per mg of myelin, respectively, for 30 min at 30 degrees C. The decrease in relative raw intensity of the higher-order reflection relative to the lower-order reflection is suggestive of a disordering of the phosphate groups upon trypsin treatment or an increased mosaicity of the membrane or a combination of both these effects, However, trypsin treatment does not lead to a complete breakdown of the membrane, The integrated intensity of the h=1 reflection, though weak, is above the measurable threshold for untreated myelin, whereas the corresponding intensity is below the measurable threshold for trypsin-treated myelin, indicating a possible asymmetric to symmetric transition of the myelin bilayer structure about its centre after trypsin treatment.

Synthesis and vesicle formation from hybrid bolaphile/amphiphile ion-pairs. Evidence of membrane property modulation by molecular design

Four new hybrid (bolaphile/amphiphile) ion-pairs were synthesized. Electron microscopy indicated that each of these forms bilayer membranes upon dispersion in aqueous media. Membrane properties have also been examined by differential scanning calorimetry, microcalorimetry, temperature-dependent fluorescence anisotropy measurements, and UV-vis spectroscopy. The T-m values for the vesicular 1, 2, 3, 4, and 5 were 38, 12, 85, 31.3, and 41.6 degrees C, respectively. Interestingly the T-m values for 1 and 3 were found to depend on their concentration. The entrapment of small solute and the release capability have also been examined to demonstrate that these bilayers form enclosed vesicles. X-ray diffraction of the cast films has been performed to understand the nature and the thickness of these membrane organizations. The membrane widths ranged from 33 to 47 Angstrom. Finally, the above observations have been analyzed in light of the results obtained from molecular modeling studies. Thus we have demonstrated that membrane properties can be modulated by simple structural changes at the amphiphile level. It was shown that by judicious incorporation of central, isomeric, disubstituted aromatic units as structural anchors into different bolaphiles, one can modulate the properties of the resulting vesicles.

Synthesis of macrocyclic diacyl/dialkyl glycerols containing disulfide tether and studies of their effects upon incorporation in DPPC membranes. Implications in the design of phospholipase A(2) modulators

A general method for the preparation of novel disulfide-tethered macrocyclic diacylglycerols (DAGs) has been described. Overall synthesis involved stepwise protection, acylation, and deprotection to yield the bis(omega-bromoacyl) glycerols. In the crucial macrocyclization step, a unique reagent, benzyltriethylammonium tetrathiomolybdate (BTAT), has been used to convert individual bis(omega-bromoacyl) glycerols to their respective macrocyclic disulfides. DAG 6, which had ether linkages between hydrocarbon chains and the glycerol backbone, was also synthesized from an appropriate precursor using a similar protocol. One of the DAGs (DAG 5) had a carbon-carbon tether instead of a disulfide one and was synthesized using modified Glaser coupling. Preparation of alpha-disulfide-tethered DAG (DAG 4) required an alternative method, as treatment of the bisbromo precursor with BTAT gave a mixture of several compounds from which separation of the target molecule was cumbersome. To avoid this problem, the bisbromide was converted to its corresponding dithiocyanate, which on further treatment with BTAT yielded the desired DAG (DAG 4) in good yield. Upon treatment with the reducing agent dithiothreitol (DTT), the DAGs that contain a disulfide tether could be quantitatively converted to their "open-chain" thiol analogues. These macrocyclic DAGs and their reduced "open-chain" analogues have been incorporated in DPPC vesicles to study their effect on model membranes. Upon incorporation of DAG 1 in DPPC vesicles, formation of new isotropic phases was observed by P-31 NMR, These isotropic phases disappeared completely on opening the macrocyclic ring by a reducing agent. The thermotropic properties of DPPC bilayers having DAGs (1-6) incorporated at various concentrations were studied by differential scanning calorimetry. Incorporation of DAGs in general reduced the cooperativity unit (CU) of the vesicles. Similar experiments with reduced "open-chain" DAGs incorporated in a DPPC bilayer indicated a recovery of CU with respect to their macrocyclic "disulfide" counterparts. The effect of inclusion of these DAGs on the activity of phospholipase A(2) (PLA(2)) was studied in vitro. Incorporation of DAC 1 in DPPC membranes potentiated both bee venom and cobra venom PLA(2) activities.

Vortex dynamics at rf frequencies in Bi2Sr2CaCu2O8 single crystals

We have studied the low magnetic field high temperature region of the H-T phase diagram of Bi2Sr2CaCu2O8 single crystals using the technique of non-resonant rf response at a frequency of 20 MHz. With H(rf)parallel to a, H parallel to c, the isothermal magnetic field scans below T-c show that the frequency f(H) of the tank circuit decreases continuously with increase in H before saturating at H similar to H-D(T). Such a decrease in f(H) reflects increasing rf penetration into the weakly screened region between CuO bilayers. The saturation of f(H) at its lowest value for H similar to H-D(T) indicates complete rf penetration land hence the disappearance of field dependence) due to the vanishing of the screening rf currents I-rf(c) in those regions or equivalently when the phase coherence between adjacent superconducting layers vanishes. Therefore H,(T) represents the decoupling of the adjacent superconducting bilayers, and hence also a 3D to 2D decoupling transition of the vortex structure. Simultaneous monitoring of the field dependent rf power dissipation P(H) shows a maximum in dP/dH at H-D(T). The observed H-D(T) line in many crystals is in excellent agreement with the (l/t-1) behavior proposed for decoupling.