Determination of elastic parameters of lipid membranes from simulation under varied external pressure

Many cellular processes such as endocytosis, exocytosis, and vesicle trafficking involve membrane deformations, which can be analyzed in the framework of the elastic theories of lipid membranes. These models operate with phenomenological elastic parameters. A connection between these parameters and the internal structure of lipid membranes can be provided by three-dimensional (3D) elastic theories. Considering a membrane as a 3D layer, Campelo et al. [F. Campelo et al., Adv. Colloid Interface Sci. 208, 25 (2014)10.1016/j.cis.2014.01.018] developed a theoretical basis for the calculation of elastic parameters. In this work we generalize and improve this approach by considering a more general condition of global incompressibility instead of local incompressibility. Crucially, we find an important correction to the theory of Campelo et al., which if not taken into account leads to a significant miscalculation of elastic parameters. With the total volume conservation taken into account, we derive an expression for the local Poisson's ratio, which determines how the local volume changes upon stretching and permits a more precise determination of elastic parameters. Also, we substantially simplify the procedure by calculating the derivatives of the moments of the local tension with respect to stretching instead of calculating the local stretching modulus. We obtain a relation between the Gaussian curvature modulus as a function of stretching and the bending modulus, showing that these two elastic parameters are not independent, as was previously assumed. The proposed algorithm is applied to membranes composed of pure dipalmitoylphosphatidylcholine (DPPC), dioleoylphosphatidylcholine (DOPC), and their mixture. The following elastic parameters of these systems are obtained: the monolayer bending and stretching moduli, spontaneous curvature, neutral surface position, and local Poisson's ratio. It is shown that the bending modulus of the DPPC/DOPC mixture follows a more complex trend than predicted by the classical Reuss averaging, which is often employed in theoretical frameworks.

Determination of Elastic Parameters of Lipid Membranes with Molecular Dynamics: A Review of Approaches and Theoretical Aspects

Lipid membranes are abundant in living organisms, where they constitute a surrounding shell for cells and their organelles. There are many circumstances in which the deformations of lipid membranes are involved in living cells: fusion and fission, membrane-mediated interaction between membrane inclusions, lipid-protein interaction, formation of pores, etc. In all of these cases, elastic parameters of lipid membranes are important for the description of membrane deformations, as these parameters determine energy barriers and characteristic times of membrane-involved phenomena. Since the development of molecular dynamics (MD), a variety of in silico methods have been proposed for the determination of elastic parameters of simulated lipid membranes. These MD methods allow for the consideration of details unattainable in experimental techniques and represent a distinct scientific field, which is rapidly developing. This work provides a review of these MD approaches with a focus on theoretical aspects. Two main challenges are identified: (i) the ambiguity in the transition from the continuum description of elastic theories to the discrete representation of MD simulations, and (ii) the determination of intrinsic elastic parameters of lipid mixtures, which is complicated due to the composition-curvature coupling effect.

Free volume and dynamics in a lipid bilayer

The lateral diffusion of lipids and of small molecules inside a membrane is strictly related to the arrangement of acyl chains and to their mobility. In this study, we use FTIR and time resolved 2D-IR spectroscopic techniques to characterize the structure and dynamics of the hydrophobic region of palmitoyl-oleylphosphatidylcholine/cholesterol vesicles dispersed in water/dimethylsulfoxide solutions. By means of a non-polar probe, hexacarbonyl tungsten, we monitor the distribution of free volumes inside the bilayer and the conformational dynamics of hydrophobic tails in relation to the different compositions of the membrane or the different compositions of the solvent. Despite the important structural changes induced by the presence of DMSO in the solvating medium, the picosecond dynamics of the membrane is preserved under the different conditions.

Chloroplast-selective gene delivery and expression in planta using chitosan-complexed single-walled carbon nanotube carriers

Plant genetic engineering is an important tool used in current efforts in crop improvement, pharmaceutical product biosynthesis and sustainable agriculture. However, conventional genetic engineering techniques target the nuclear genome, prompting concerns about the proliferation of foreign genes to weedy relatives. Chloroplast transformation does not have this limitation, since the plastid genome is maternally inherited in most plants, motivating the need for organelle-specific and selective nanocarriers. Here, we rationally designed chitosan-complexed single-walled carbon nanotubes, utilizing the lipid exchange envelope penetration mechanism. The single-walled carbon nanotubes selectively deliver plasmid DNA to chloroplasts of different plant species without external biolistic or chemical aid. We demonstrate chloroplast-targeted transgene delivery and transient expression in mature Eruca sativa, Nasturtium officinale, Nicotiana tabacum and Spinacia oleracea plants and in isolated Arabidopsis thaliana mesophyll protoplasts. This nanoparticle-mediated chloroplast transgene delivery tool provides practical advantages over current delivery techniques as a potential transformation method for mature plants to benefit plant bioengineering and biological studies.

Application of liposome membrane as the reaction field: A case study using the Horner-Wadsworth-Emmons reaction

The properties of the liposome membrane as a reaction field were investigated by focusing on the Horner-Wadsworth-Emmons reaction as a case study. Use of the liposomes existing in the gel phase resulted in the enhanced activity of the substrates and furnished the products with same E/Z stereoselectivity as in the liposome-free system. The membrane environment in the gel phase most likely assisted the formation of adducts that induced selective generation of the E-isomer. The possible role of liposomes is to assist the proton removal from the reactant, rather than providing the basic interfacial environment.

Rational Design Principles for the Transport and Subcellular Distribution of Nanomaterials into Plant Protoplasts

The ability to control the subcellular localization of nanoparticles within living plants offers unique advantages for targeted biomolecule delivery and enables important applications in plant bioengineering. However, the mechanism of nanoparticle transport past plant biological membranes is poorly understood. Here, a mechanistic study of nanoparticle cellular uptake into plant protoplasts is presented. An experimentally validated mathematical model of lipid exchange envelope penetration mechanism for protoplasts, which predicts that the subcellular distribution of nanoparticles in plant cells is dictated by the particle size and the magnitude of the zeta potential, is advanced. The mechanism is completely generic, describing nanoparticles ranging from quantum dots, gold and silica nanoparticles, nanoceria, and single-walled carbon nanotubes (SWNTs). In addition, the use of imaging flow cytometry to investigate the influence of protoplasts' morphological characteristics on nanoparticle uptake efficiency is demonstrated. Using DNA-wrapped SWNTs as model nanoparticles, it is found that glycerolipids, the predominant lipids in chloroplast membranes, exhibit stronger lipid-nanoparticle interaction than phospholipids, the major constituent in protoplast membrane. This work can guide the rational design of nanoparticles for targeted delivery into specific compartments within plant cells without the use of chemical or mechanical aid, potentially enabling various plant engineering applications.

Controlled release of nitric oxide from liposomes

We report the ability to readily tune NO release from N-diazeniumdiolate-encapsulated liposomal structures by altering the NO donor molecule structure and/or phospholipid composition (independently or in combination). While encapsulating more stable NO donors expectedly enhanced the NO release (up to 48 h) from the liposomes, the phospholipid headgroup surface area proved equally useful in controlling NO-release kinetics by influencing the water uptake and concomitant N-diazeniumdiolate NO donor breakdown (to NO). The potential therapeutic utility of the NO-releasing liposomes was further assessed in biological/proteinaceous fluids. The NO-release kinetics were similar in buffer and serum.

Solubility and diffusion of oxygen in phospholipid membranes

The transport of oxygen and other nonelectrolytes across lipid membranes is known to depend on both diffusion and solubility in the bilayer, and to be affected by changes in the physical state and by the lipid composition, especially the content of cholesterol and unsaturated fatty acids. However, it is not known how these factors affect diffusion and solubility separately. Herein we measured the partition coefficient of oxygen in liposome membranes of dilauroyl-, dimiristoyl- and dipalmitoylphosphatidylcholine in buffer at different temperatures using the equilibrium-shift method with electrochemical detection. The apparent diffusion coefficient was measured following the fluorescence quenching of 1-pyrenedodecanoate inserted in the liposome bilayers under the same conditions. The partition coefficient varied with the temperature and the physical state of the membrane, from below 1 in the gel state to above 2.8 in the liquid-crystalline state in DMPC and DPPC membranes. The partition coefficient was directly proportional to the partial molar volume and was then associated to the increase in free-volume in the membrane as a function of temperature. The apparent diffusion coefficients were corrected by the partition coefficients and found to be nearly the same, with a null dependence on viscosity and physical state of the membrane, probably because the pyrene is disturbing the surrounding lipids and thus becoming insensitive to changes in membrane viscosity. Combining our results with those of others, it is apparent that both solubility and diffusion increase when increasing the temperature or when comparing a membrane in the gel to one in the fluid state.

How Tolerant are Membrane Simulations with Mismatch in Area per Lipid between Leaflets?

Difficulties in estimating the correct number of lipids in each leaflet of complex bilayer membrane simulation systems make it inevitable to introduce a mismatch in lipid packing (i.e., area per lipid) and thus alter the lateral pressure of each leaflet. To investigate potential impacts of such mismatch on simulation results, we performed molecular dynamics simulations of saturated and monounsaturated lipid bilayers with and without gramicidin A or WALP23 at various mismatches by adjusting the number of lipids in the lower leaflet from no mismatch to a 25% reduction compared to that in the upper leaflet. All simulations were stable under the constant pressure barostat, but the mismatch induces asymmetric lipid packing between the leaflets, so that the upper leaflet becomes more ordered, and the lower leaflet becomes less ordered. The mismatch impacts on various bilayer properties are mild up to 5–10% mismatch, and bilayers with fully saturated chains appear to be more prone to these impacts than those with unsaturated tails. The nonvanishing leaflet surface tensions and the free energy derivatives with respect to the bilayer curvature indicate that the bilayer would be energetically unstable in the presence of mismatch. We propose a quantitative criterion for allowable mismatch based on the energetics derived from a continuum elastic model, which grows as a square root of the number of the lipids in the system. On the basis of this criterion, we infer that the area per lipid mismatch up to 5% would be tolerable in various membrane simulations of reasonable all-atom system sizes (40–160 lipids per leaflet).

Lipid Exchange Envelope Penetration (LEEP) of Nanoparticles for Plant Engineering: A Universal Localization Mechanism

Nanoparticles offer clear advantages for both passive and active penetration into biologically important membranes. However, the uptake and localization mechanism of nanoparticles within living plants, plant cells, and organelles has yet to be elucidated.1 Here, we examine the subcellular uptake and kinetic trapping of a wide range of nanoparticles for the first time, using the plant chloroplast as a model system, but validated in vivo in living plants. Confocal visible and near-infrared fluorescent microscopy and single particle tracking of gold-cysteine-AF405 (GNP-Cys-AF405), streptavidin-quantum dot (SA-QD), dextran and poly(acrylic acid) nanoceria, and various polymer-wrapped single-walled carbon nanotubes (SWCNTs), including lipid-PEG-SWCNT, chitosan-SWCNT and 30-base (dAdT) sequence of ssDNA (AT)15 wrapped SWCNTs (hereafter referred to as ss(AT)15-SWCNT), are used to demonstrate that particle size and the magnitude, but not the sign, of the zeta potential are key in determining whether a particle is spontaneously and kinetically trapped within the organelle, despite the negative zeta potential of the envelope. We develop a mathematical model of this lipid exchange envelope and penetration (LEEP) mechanism, which agrees well with observations of this size and zeta potential dependence. The theory predicts a critical particle size below which the mechanism fails at all zeta potentials, explaining why nanoparticles are critical for this process. LEEP constitutes a powerful particulate transport and localization mechanism for nanoparticles within the plant system.

Insights into chitosan multiple functional properties: the role of chitosan conformation in the behavior of liposomal membrane

Interactions of chitosan with liposomes correlate with multiple functionalities. Chitosan chains can self-aggregate above a critical aggregation concentration. The physical properties of liposomes are affected by chitosan conformation. Chitosan displays “polymeric surfactant property” in the form of coils.

Ceramide increases free volume voids in DPPC membranes

We have measured by positron annihilation lifetime spectroscopy (PALS) that ceramide increases the size of the free volume holes in DPPC lipid membranes.

Agomelatine strongly interacts with zwitterionic DPPC and charged DPPG membranes

Depression is one of the most common psychiatric diseases in the population. Agomelatine is a novel antidepressant drug with melatonin receptor agonistic and serotonin 5-HT2C antagonistic properties. Furthermore, being a melatonergic drug, agomelatine has the potential of being used in therapeutic applications like melatonin as an antioxidant, anti-inflammatory and antiapoptotic drug. The action mechanism of agomelatine on the membrane structure has not been clarified yet. In the present study, we aimed to investigate the interaction of agomelatine with model membranes of dipalmitoylphosphatidylcholine (DPPC) and dipalmitoylphosphatidylgylcerol (DPPG) by Fourier transform infrared (FTIR) spectroscopy and differential scanning calorimetry (DSC). We found that agomelatine interacts with the head group in such a manner that it destabilizes the membrane architecture to a large extent. Thus, agomelatine causes alterations in the order, packing and dynamics of the DPPC and DPPG model membranes. Our results suggest that agomelatine strongly interacts with zwitterionic and charged membrane phospholipids. Because lipid structure and dynamics may have influence on the structure of membrane bound proteins and affect the signal transduction systems of membranes, these effects of agomelatine may be important in its action mechanism.

Poly(ethylene glycol) lipid-shelled microbubbles: abundance, stability, and mechanical properties

Poly(ethylene glycol) (PEG) is widely used on the outside of biomedical delivery vehicles to impart stealth properties. Encapsulated gas microbubbles (MBs) are being increasingly considered as effective carriers for therapeutic intervention to deliver drug payloads or genetic vectors. MBs have the advantage that they can be imaged and manipulated by ultrasound fields with great potential for targeted therapy and diagnostic purposes. Lipid-shelled MBs are biocompatible and can be functionalized on the outer surface for tissue targeting and new therapeutic methods. As MBs become a key route for drug delivery, exploring the effect of PEG-ylation on the MB properties is important. Here, we systematically investigate the effect of PEG-lipid solution concentration ranging between 0 and 35 mol % on the formation of MBs in a microfluidic flow-focusing device. The abundance of the MBs is correlated with the MB lifetime and the whole MB mechanical response, as measured by AFM compression using a tipless cantilever. The maximal MB concentration and stability (lifetime) occurs at a low concentration of PEG-lipid (∼5 mol %). For higher PEG-lipid concentrations, the lifetime and MB concentration decrease, and are accompanied by a correlation between the predicted surface PEG configuration and the whole MB stiffness, as measured at higher compression loads. These results inform the rationale design and fabrication of lipid-based MBs for therapeutic applications and suggest that only relatively small amounts of PEG incorporation are required for optimizing MB abundance and stability while retaining similar mechanical response at low loads.

Nanoscale packing differences in sphingomyelin and phosphatidylcholine revealed by BODIPY fluorescence in monolayers: physiological implications

Phosphatidycholines (PC) with two saturated acyl chains (e.g., dipalmitoyl) mimic natural sphingomyelin (SM) by promoting raft formation in model membranes. However, sphingoid-based lipids, such as SM, rather than saturated-chain PCs have been implicated as key components of lipid rafts in biomembranes. These observations raise questions about the physical packing properties of the phase states that can be formed by these two major plasma membrane lipids with identical phosphocholine headgroups. To investigate, we developed a monolayer platform capable of monitoring changes in surface fluorescence by acquiring multiple spectra during measurement of a lipid force-area isotherm. We relied on the concentration-dependent emission changes of 4,4-difluoro-4-bora-3a,4a-diaza-s-indacene (BODIPY)-labeled PC to detect nanoscale alterations in lipid packing and phase state induced by monolayer lateral compression. The BODIPY-PC probe contained an indacene ring with four symmetrically located methyl (Me) substituents to enhance localization to the lipid hydrocarbon region. Surface fluorescence spectra indicated changes in miscibility even when force-area isotherms showed no deviation from ideal mixing behavior in the surface pressure versus cross-sectional molecular area response. We detected slightly better mixing of Me4-BODIPY-8-PC with the fluid-like, liquid expanded phase of 1-palmitoyl-2-oleoyl-PC compared to N-oleoyl-SM. Remarkably, in the gel-like, liquid condensed phase, Me4-BODIPY-8-PC mixed better with N-palmitoyl-SM than dipalmitoyl-PC, suggesting naturally abundant SMs with saturated acyl chains form gel-like lipid phase(s) with enhanced ability to accommodate deeply embedded components compared to dipalmitoyl-PC gel phase. The findings reveal a fundamental difference in the lateral packing properties of SM and PC that occurs even when their acyl chains match.

Extracting curvature preferences of lipids assembled in flat bilayers shows possible kinetic windows for genesis of bilayer asymmetry and domain formation in biological membranes

Studies on the assembly of pure lipid components allow mechanistic insights toward understanding the structural and functional aspects of biological membranes. Molecular dynamic (MD) simulations on membrane systems provide molecular details on membrane dynamics that are difficult to obtain experimentally. A large number of MD studies have remained somewhat disconnected from a key concept of amphipathic assembly resulting in membrane structures--shape parameters of lipid molecules in those structures in aqueous environments. This is because most of the MD studies have been done on flat lipid membranes. With the above in view, we analyzed MD simulations of 26 pure lipid patches as a function of (1) lipid type(s) and (2) time of MD simulations along with 35-40 ns trajectories of five pure lipids. We report, for the first time, extraction of curvature preferences of lipids from MD simulations done on flat bilayers. Our results may lead to mechanistic insights into the possible origins of bilayer asymmetries and domain formation in biological membranes.

Pressure modulation of Ras-membrane interactions and intervesicle transfer

Proteins attached to the plasma membrane frequently encounter mechanical stresses, including high hydrostatic pressure (HHP) stress. Signaling pathways involving membrane-associated small GTPases (e.g., Ras) have been identified as critical loci for pressure perturbation. However, the impact of mechanical stimuli on biological outputs is still largely terra incognita. The present study explores the effect of HHP on the membrane association, dissociation, and intervesicle transfer process of N-Ras by using a FRET-based assay to obtain the kinetic parameters and volumetric properties along the reaction path of these processes. Notably, membrane association is fostered upon pressurization. Conversely, depending on the nature and lateral organization of the lipid membrane, acceleration or retardation is observed for the dissociation step. In addition, HHP can be inferred as a positive regulator of N-Ras clustering, in particular in heterogeneous membranes. The susceptibility of membrane interaction to pressure raises the idea of a role of lipidated signaling molecules as mechanosensors, transducing mechanical stimuli to chemical signals by regulating their membrane binding and dissociation. Finally, our results provide first insights into the influence of pressure on membrane-associated Ras-controlled signaling events in organisms living under extreme environmental conditions such as those that are encountered in the deep sea and sub-seafloor environments, where pressures reach the kilobar (100 MPa) range.

Early stages of interactions of cell-penetrating peptide penetratin with a DPPC bilayer

Cell-penetrating-peptides (CPPs) can deliver themselves together with a macromolecular cargo into cells and, hence, have promising applications in drug delivery. The detailed physical mechanisms that underlie and determine their cellular uptake remain unknown. We used molecular dynamics (MD) simulations to study the interaction of a well-known CPP, namely penetratin, with a zwitterionic di-palmitoyl-phosphatidyl-choline (DPPC) bilayer. Our study shows that the arginine and lysine residues play a crucial role in peptide-membrane binding through charge-pair and hydrogen bond interactions. We also characterize peptide conformation and show that it remains helical near the N-terminus but can fold into a variety of other conformations in the residues close to the C-terminus. The response of membrane to the peptide is also investigated.

Binding of serotonin to lipid membranes

Serotonin (5-hydroxytryptamine, 5-HT) is a prevalent neurotransmitter throughout the animal kingdom. It exerts its effect through the specific binding to the serotonin receptor, but recent research has suggested that neural transmission may also be affected by its nonspecific interactions with the lipid matrix of the synaptic membrane. However, membrane-5-HT interactions remain controversial and superficially investigated. Fundamental knowledge of this interaction appears vital in discussions of putative roles of 5-HT, and we have addressed this by thermodynamic measurements and molecular dynamics (MD) simulations. 5-HT was found to interact strongly with lipid bilayers (partitioning coefficient ~1200 in mole fraction units), and this is highly unusual for a hydrophilic solute like 5-HT which has a bulk, oil-water partitioning coefficient well below unity. It follows that membrane affinity must rely on specific interactions, and the MD simulations identified the salt-bridge between the primary amine of 5-HT and the lipid phosphate group as the most important interaction. This interaction anchored cationic 5-HT in the membrane interface with the aromatic ring system pointing inward and a prevailing residence between the phosphate and the carbonyl groups of the lipid. The unprotonated form of 5-HT shows the opposite orientation, with the primary amine pointing toward the membrane core. Partitioning of 5-HT was found to decrease lipid chain order. These distinctive interactions of 5-HT and model membranes could be related to nonspecific effects of this neurotransmitter.

Molecular dynamic studies of transportan interacting with a DPPC lipid bilayer

Translocation of peptides through cellular membranes is a fundamental problem in developing antimicrobial peptides and in drug delivery. There is a class of peptides, known as cell-penetrating peptides, that are able to penetrate membranes without disrupting them. They can carry pharmacological compounds, thus a promising strategy for drug delivery. The physical mechanisms that facilitate translocation are not known. We have used large-scale molecular dynamics simulations to study the penetration of transportan across a zwitterionic dipalmitoyl-phosphatidyl-choline (DPPC) bilayer. We obtained the free energy profile for one peptide inside the bilayer and discuss the response of the bilayer to the presence of transportan. We also discuss the importance of lysine residues and speculate on the possible penetration mechanism of the peptide and propose a graded-like penetration process.

Fluorescence study on aggregated lysozyme and lipid bilayer interactions

Fluorescent probes 1,6-diphenyl-1,3,5-hexatriene (DPH), pyrene, 4-dimethylaminochalcone (DMC) and 4-p-(dimethylaminostyryl)-1-dodecylpyridinium (DSP-12) have been utilized to monitor the impact of lysozyme (Lz) oligomers on physicochemical properties of phosphatidylcholine/cardiolipin (PC/CL) membranes. Analysis of spectral responses of the employed probes revealed the reduction of membrane free volume and dehydration of lipid bilayer surface upon incorporation of Lz self-assemblies. Hydrophobic interactions were found to control the binding of Lz oligomers to the lipid bilayer. Comparison of the effects of Lz monomers, oligomers and fibrils showed that soluble oligomeric intermediates exert the most destructive influence on membrane properties.

Interaction of hexadecylbetainate chloride with biological relevant lipids

The present work investigates the interaction of hexadecylbetainate chloride (C(16)BC), a glycine betaine-based ester with palmitoyl-oleoyl-phosphatidylcholine (POPC), sphingomyelin (SM), and cholesterol (CHOL), three biological relevant lipids present in the outer leaflet of the mammalian plasma membrane. The binding affinity and the mixing behavior between the lipids and C(16)BC are discussed based on experimental (isothermal titration calorimetry (ITC) and Langmuir film balance) and molecular modeling studies. The results show that the interaction between C(16)BC and each lipid is thermodynamically favorable and does not affect the integrity of the lipid vesicles. The primary adsorption of C(16)BC into the lipid film is mainly governed by a hydrophobic effect. Once C(16)BC is inserted in the lipid film, the polar component of the interaction energy between C(16)BC and the lipid becomes predominant. Presence of CHOL increases the affinity of C(16)BC for membrane. This result can be explained by the optimal matching between C(16)BC and CHOL within the film rather by a change of membrane fluidity due to the presence of CHOL. The interaction between C(16)BC and SM is also favorable and gives rise to highly stable monolayers probably due to hydrogen bonds between their hydrophilic groups. The interaction of C(16)BC with POPC is less favorable but does not destabilize the mixed monolayer from a thermodynamic point of view. Interestingly, for all the monolayers investigated, the exclusion surface pressures are above the presumed lateral pressure of the plasma membranes suggesting that C(16)BC would be able to penetrate into mammalian plasma membranes in vivo. These results may serve as a useful basis in understanding the interaction of C(16)BC with real membranes.

Monomeric fullerenes in lipid membranes: effects of molecular shape and polarity

We report a combined theoretical and experimental study on the single-molecule interaction of fullerenes with phospholipid membranes. We studied pristine C(60) (1) and two N-substituted fulleropyrrolidines (2 and 3), one of which (3) bore a paramagnetic nitroxide group. Theoretical predictions of fullerene distribution and permeability across lipid bilayers were combined with electron paramagnetic resonance (EPR) experiments in aligned DMPC/DHPC bicelles containing the paramagnetic fulleropyrrolidine 3 or either one of the diamagnetic fullerenes together with spin-labeled lipids. We found that, at low concentrations, fullerenes are present in the bilayer as single molecules. Their preferred location in the membrane is only slightly influenced by the derivatization: all derivatives were confined just below the hydrophilic/hydrophobic interface, because of the key role played by dispersion interactions between the highly polarizable fullerene cage and the hydrocarbon chains, which are especially tight within this region. However, the deviation from spherical shape is sufficient to induce a preferential orientation of 2 and 3 in the membrane. We predict that monomeric fullerenes spontaneously penetrate the bilayer, in agreement with the results of molecular dynamics simulations, but we point out the limits of the currently used permeability model when applied to hydrophobic solutes.

Polarity-sensitive fluorescent probes in lipid bilayers: bridging spectroscopic behavior and microenvironment properties

We have studied the emission features of the fluorescent polarity-sensitive probes known as Prodan and Laurdan in a liquid-crystalline DPPC bilayer. To this purpose, we have combined high-level quantum mechanical electronic structure calculations with a molecular field theory for the positional-orientational-conformational distribution of the probes, in their ground and excited states, inside of the lipid bilayer, taking into account at both levels the nonuniformity and anisotropy of the environment. Thus, we can interpret the features of the fluorescence spectra of Prodan and Laurdan in relation to the position and orientation of their chromophore in the bilayer. We have found that the environment polarity is not sufficient to explain the large red shifts experimentally observed and that specific effects due to hydrogen bonding must be considered. We show that the orientation of the probe is important in determining the accessibility to water of the H-bond-acceptor group; in the case of Laurdan interesting conformational effects are highlighted.

Free volume theory applied to lateral diffusion in Langmuir monolayers: atomistic simulations for a protein-free model of lung surfactant

We hereby present a study on lateral diffusion of lipids in Langmuir monolayers. We apply atomistic molecular dynamics simulations to a model system whose composition is consistent with protein-free lung surfactant. Our main focus is on the assessment of the validity of the free volume theory for lateral diffusion and on the interpretation of the cross-sectional area and activation energy parameters appearing in the theory. We find that the diffusion results can be fitted to the description given by the free volume theory, but the interpretation of its parameters is not straightforward. While the cross-sectional area appears to be related to the hard-core cross-sectional area of a lipid, its role in the lateral diffusion process is unclear. Also, the activation energy derived using the free volume theory is different from the activation energy found through Arrhenius analysis, and its physical interpretation remains elusive. Finally, we find that lipid diffusion does not occur via rapid single-particle "jumps". Instead, lipids move in a concerted manner as loosely defined transient clusters, as observed earlier for lipid bilayers.

Self-organization of a stable pore structure in a phospholipid bilayer

We demonstrate the self-organization process of a stable pore structure in a phospholipid bilayer by unsteady and nonequilibrium molecular dynamics simulations. The simulation is started from an initial state including some amount of water molecules in its hydrophobic region, which is a model of a cell membrane stimulated by ultrasound radiation for the membrane permeabilization (sonoporation). We show that, in several nanoseconds, the bilayer-water system can spontaneously develop into a water-filled pore structure without any mechanical and electrical forcing from outside, when the initial number of water molecules in the hydrophobic region exceeds a critical value. The increase in the initial number of water molecules enhances the probability of pore formation, and sometimes induces the formation of transient micellelike structures of phospholipid molecules.

Temperature-induced structural transition in-situ in porcine lens--changes observed in void size distribution

The function of mammalian ocular lens is to provide a sharp image to the retina. Accordingly, the lens needs to be transparent and minimize light scattering. To do so the lens fiber cells first loose intracellular organelles, organize the cytoplasm and arrange the fiber cell membranes. Because the fiber cells are metabolically inactive, the plasma membrane becomes the only cellular organelle and consequently, the phase behavior of these membranes determines the physiological state of the lens. Previous studies have shown that lipids extracted from the nuclear and cortical region of human lens show a temperature-induced phase transition close to the body temperature. Yet, the physiological function of this phase transition is not known, and even the presence of the phase transition in intact lenses is unknown. Positron annihilation lifetime spectroscopy (PALS) was used to characterize the sub-nanometer-sized local structure of intact porcine lens and these studies were complemented with differential scanning calorimeter and mass spectrometric analysis in extracted porcine lens lipids. Using PALS, we present evidence for the presence of a temperature-dependent structural transition centered at 35.5 degrees C in-situ in clear extracted porcine lenses. Further studies employing extracted lens lipids and purified egg-yolk sphingomyelin and cholesterol mixtures suggest that the nano-scale transition emerges from the phase behavior of lens lipids. Based on our results, PALS seems to be a viable method for gaining additional information on biological tissues, especially since it enables non-destructive studies on intact tissues.

Atomistic simulations of phosphatidylcholines and cholesteryl esters in high-density lipoprotein-sized lipid droplet and trilayer: clues to cholesteryl ester transport and storage

Cholesteryl esters (CEs) are the water-insoluble transport and storage form of cholesterol. For both transport and storage, phospholipids and proteins embrace the CEs to form an amphipathic monolayer that surrounds the CEs. CEs are transported extracellularly in lipoproteins and are stored intracellularly as cytoplasmic lipid droplets. To clarify the molecular phenomena related to the above structures, we conducted atomistic molecular-dynamics simulations for a spherical, approximately high density lipoprotein sized lipid droplet comprised of palmitoyl-oleoyl-phosphatidylcholine (POPC) and cholesteryl oleate (CO) molecules. An additional simulation was conducted for a lamellar lipid trilayer consisting of the same lipid constituents. The density profiles showed that COs were located in the core of the spherical droplet. In trilayer simulations, CO molecules were also in the core and formed two denser strata. This is remarkable because the intra- and intermolecular behaviors of the COs were similar to previous findings from bulk COs in the fluid phase. In accordance with previous experimental studies, the solubility of COs in the POPC monolayers was found to be low. The orientation distribution of the sterol moiety with respect to the normal of the system was found to be broad, with mainly isotropic or slightly parallel orientations observed deep in the core of the lipid droplet or the trilayer, respectively. In both systems, the orientation of the sterol moiety changed to perpendicular with respect to the normal close to the phopsholipid monolayers. Of interest, within the POPC monolayers, the intramolecular conformation of the COs varied from the previously proposed horseshoe-like conformation to a more extended one. From a metabolic point of view, the observed solubilization of CEs into the phospholipid monolayers, and the conformation of CEs in the phospholipid monolayers are likely to be important regulatory factors of CE transport and hydrolysis.

Atom-scale molecular interactions in lipid raft mixtures

We review the relationship between molecular interactions and the properties of lipid environments. A specific focus is given on bilayers which contain sphingomyelin (SM) and sterols due to their essential role for the formation of lipid rafts. The discussion is based on recent atom-scale molecular dynamics simulations, complemented by extensive comparison to experimental data. The discussion is divided into four sections. The first part investigates the properties of one-component SM bilayers and compares them to bilayers with phosphatidylcholine (PC), the focus being on a detailed analysis of the hydrogen bonding network in the two bilayers. The second part deals with binary mixtures of sterols with either SM or PC. The results show how the membrane properties may vary substantially depending on the sterol and SM type available, the membrane order and interdigitation being just two of the many examples of this issue. The third part concentrates on the specificity of intermolecular interactions in three-component mixtures of SM, PC and cholesterol (CHOL) under conditions where the concentrations of SM and CHOL are dilute with respect to that of PC. The results show how SM and CHOL favor one another, thus acting as nucleation sites for the formation of highly ordered nanosized domains. Finally, the fourth part discusses the large-scale properties of raft-like membrane environments and compares them to the properties of non-raft membranes. The differences turn out to be substantial. As a particularly intriguing example of this, the lateral pressure profiles of raft-like and non-raft systems indicate that the lipid composition of membrane domains may have a major impact on membrane protein activation.

Liquid ordered and gel phases of lipid bilayers: fluorescent probes reveal close fluidity but different hydration

Hydration and fluidity of lipid bilayers in different phase states were studied using fluorescent probes selectively located at the interface. The probe of hydration was a recently developed 3-hydroxyflavone derivative, which is highly sensitive to the environment, whereas the probe of fluidity was the diphenylhexatriene derivative, 1-[4-(trimethylamino)phenyl]-6-phenylhexa-1,3,5-triene. By variation of the cholesterol content and temperature in large unilamellar vesicles composed of sphingomyelin or dipalmitoylphosphatidlycholine, we generated different phases: gel, liquid ordered (raft), liquid crystalline, and liquid disordered (considered as liquid crystalline phase with cholesterol). For these four phases, the hydration increases in the following order: liquid ordered << gel approximately liquid disordered < liquid crystalline. The membrane fluidity shows a somewhat different trend, namely liquid ordered approximately gel < liquid disordered < liquid crystalline. Thus, gel and liquid ordered phases exhibit similar fluidity, whereas the last phase is significantly less hydrated. We expect that cholesterol due to its specific H-bonding interactions with lipids and its ability to fill the voids in lipid bilayers expels efficiently water molecules from the highly ordered gel phase to form the liquid ordered phase. In this study, the liquid ordered (raft) and gel phases are for the first time clearly distinguished by their strong difference in hydration.