Interleaflet interaction and asymmetry in phase separated lipid bilayers: molecular dynamics simulations

In order to investigate experimentally inaccessible, molecular-level detail regarding interleaflet interaction in membranes, we have run an extensive series of coarse-grained molecular dynamics simulations of phase separated lipid bilayers. The simulations are motivated by differences in lipid and cholesterol composition in the inner and outer leaflets of biological membranes. Over the past several years, this phenomenon has inspired a series of experiments in model membrane systems which have explored the effects of lipid compositional asymmetry in the two leaflets. The simulations are directed at understanding one potential consequence of compositional asymmetry, that being regions of bilayers where liquid-ordered (L(o)) domains in one leaflet are opposite liquid-disordered (L(d)) domains in the other leaflet (phase asymmetry). The simulated bilayers are of two sorts: 1) Compositionally symmetric leaflets where each of the two leaflets contains an identical, phase separated (L(o)/L(d)) mixture of cholesterol, saturated and unsaturated phospholipid; and 2) Compositionally asymmetric leaflets, where one leaflet contains a phase separated (L(o)/L(d)) mixture while the other contains only unsaturated lipid, which on its own would be in the L(d) phase. In addition, we have run simulations where the lengths of the saturated lipid chains as well as the mole ratios of the three lipid components are varied. Collectively, we report on three types of interleaflet coupling within a bilayer. First, we show the effects of compositional asymmetry on acyl chain tilt and order, lipid rotational dynamics, and lateral diffusion in regions of leaflets that are opposite L(o) domains. Second, we show substantial effects of compositional asymmetry on local bilayer curvature, with the conclusion that phase separated leaflets resist curvature, while inducing large degrees of curvature in an opposing L(d) leaflet. Finally, in compositionally symmetric, phase separated bilayers, we find phase asymmetry (domain antiregistration) between the two leaflets occurs as a consequence of mismatched acyl chain-lengths in the saturated and unsaturated lipids.     

Dynamic scaling in phase separation kinetics for quasi-two-dimensional membranes

We consider the dynamics of phase separation in lipid bilayer membranes, modeled as flat two-dimensional liquid sheets within a bulk fluid, both in the creeping flow approximation. We present scaling arguments that suggest asymptotic coarsening in these systems is characterized by a length scale R(t) ~ t(1/2) for critical (bicontinuous) phase separation and R(t) ~t(1/3) for off-critical concentrations (droplet morphology). In this limit, the bulk fluid is the primary source of dissipation. We also address these questions with continuum stochastic hydrodynamic simulations. We see evidence of scaling violation in critical phase separation, where isolated circular domains coarsen slower than elongated ones. However, we also find a region of apparent scaling where R(t) ~ t(1/2) is observed. This appears to be due to the competition of thermal and hydrodynamic effects. We argue that the diversity of scaling exponents measured in experiment and prior simulations can in part be attributed to certain measurements lying outside the asymptotic long-length-scale regime, and provide a framework to help understand these results. We also discuss a few simple generalizations to confined membranes and membranes in which inertia is relevant.     

Domain Size and Fluctuations at Domain Interfaces in Lipid Mixtures

Biomembranes consist of a complex mixture of a large number of lipids and proteins. In such mixtures, microscopic domains and macroscopically separated phases may exist. Here, we discuss phase behavior and domains formation of binary lipid mixtures. We show that the domain formation is accompanied by large fluctuations at the domain boundaries, resulting in altered physical properties at the boundaries, for instance in a pronounced increase of the elastic constants. Therefore, we argue that the physics of the membrane depends on the overall length scale of its domains interfaces. We present here confocal microscopy images, calorimetric melting profiles and Monte-Carlo simulations to understand the factors that determine domain formation, their sizes and the role of the domain interfaces.     

Unforced translocation of a polymer chain through a nanopore: the solvent effect

The authors have performed the Langevin dynamics simulation to investigate the unforced polymer translocation through a narrow nanopore in an impermeable membrane. The effects of solvent quality controlled by the attraction strength lambda of the Lennard-Jones cosine potential between polymer beads and beads on two sides of the membrane on the translocation processes are extensively examined. For polymer translocation under the same solvent quality on both sides of the membrane, the two-dimensional and three-dimensional simulations confirm the scaling law of tautrans approximately N1+2upsilon for the translocation in the good solvent, where tautrans is the translocation time, N is the chain length, and upsilon is the Flory exponent. For the three-dimensional polymer translocation under different solvent qualities on two sides of the membrane, the translocation efficiency may be notably improved. The scaling law between tautrans and N varies from tautrans approximately N1+2upsilon to tautrans approximately N with the increase of the difference of solvent qualities, and the crossover occurs at the theta temperature point, where a scaling law of tautrans approximately N1.27 is found. The simulation results here also show that the translocation time changes from a wide and asymmetric distribution with a long tail to a narrow and symmetric distribution with the increase of the difference of the solvent qualities.     

Giant phospholipid/block copolymer hybrid vesicles: mixing behavior and domain formation

Lipids and block copolymers can be individually assembled into unsupported, spherical membranes (liposomes or polymersomes), each having their own particular benefits and limitations. Here we demonstrate the preparation of microscale, hybrid "lipopolymersomes" composed of the common lipid POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine) and the commercially available copolymer PBd-b-PEO (polybutadiene-b-poly(ethylene oxide)) with the goal of incorporating the advantageous qualities of the unitary systems into mixed-membrane capsules. We investigate the lipopolymersomes using confocal fluorescence microscopy and demonstrate that these hybrid membranes are well mixed on nanoscopic length scales within the permittable compositional windows for hybrid vesicle formation. We measure the intramembrane dynamics and mechanical properties of these hybrid membranes by fluorescence recovery after photobleaching (FRAP) and micropipet aspiration, respectively. For the first time, we demonstrate the demixing of lipid-rich and polymer-rich membrane domains within the same vesicle membrane. This is achieved by the biotinylation of one of the constituent species and cross linking with the protein NeutrAvidin. The resultant domain patterning is dependent upon which component carries the biotin functionality: cross linking of the copolymer species results in domains that ripen into a single, large, copolymer-rich island, and cross linking of the lipids yields many small, "spot-like", lipid-rich domains within a copolymer-rich matrix. We discuss these morphological differences in terms of the fluidity and mechanical properties of the membrane phases and the possible resultant interdomain interactions within the membrane. These heterogeneous hybrid lipopolymersomes could find applications in fields such as targeted delivery, controlled release, and environmental detection assays where these capsules possess the characteristics of biocompatible lipid membranes combined with enhanced mechanical strength and stability from the copolymer matrix.     

Molecular dynamics of nicotinic acetylcholine receptor correlating biological functions

The nicotinic acetylcholine receptor (nAChR) that mediates fast intercellular communication in response to neurotransmitters is a paradigm of ligand-gated ion channels. Molecular dynamics (MD) simulations are valuable in understanding membrane protein function at atomic level, providing useful clues for further experimental/theoretical studies. In this brief review, recent progress in MD simulations of the nAChR has been illustrated, mainly focusing on the latest simulation of the whole transmembrane domain of the receptor. On the basis of MD simulations, asymmetrical and asynchronous motions of five subunits were observed both in the ligand binding and transmembrane domains; a closed-to-open conformational shift of the gate was captured in different simulation systems; the contributions from the lipid molecules and other transmembrane segments rather than M2 to the gate switch as well as the conformational change of the whole channel were assessed; the dynamic behavior and related physical/chemical properties of the water molecules and cations within the ion channel were examined; and an experimentally comparable single-channel conductance and ion selectivity were obtained.     

Membrane protein topology probed by (1)H spin diffusion from lipids using solid-state NMR spectroscopy

We describe a two-dimensional solid-state NMR technique to investigate membrane protein topology under magic-angle spinning conditions. The experiment detects the rate of (1)H spin diffusion from the mobile lipids to the rigid protein. While spin diffusion within the rigid protein is fast, magnetization transfer in the mobile lipids is an inefficient and slow process. Qualitative analysis of (1)H spin-diffusion build-up curves from the lipid chain-end methyl groups to the protein allows the identification of membrane-embedded domains in the protein. Numerical simulations of spin-diffusion build-up curves yield the approximate insertion depth of protein segments in the membrane. The experiment is demonstrated on the selectively (13)C labeled colicin Ia channel domain, known to have a membrane-embedded domain, and on DNA/cationic lipid complexes where the DNA rods are bound to the membrane surface. The experiment is designed for X-nucleus detection, which could be (13)C or (15)N in the protein and (31)P for the DNA. Finally, we show that a qualitative distinction between membrane proteins with and without a membrane-embedded domain can be made even by using an unlabeled protein, by detection of lipid signals. This spin-diffusion experiment is simple to perform and requires no oriented bilayer preparations and only standard NMR hardware.     

DPD Simulations of PMMA-Oleic Acid Mixture Behaviour in Organic Capped Nanoparticle Based Polymer Nanocomposite

Dissipative Particle Dynamics has been used to investigate the different morphology of polymer nanocomposites. Such a study was addressed to the definition of a suitable tool for understanding the distribution of oleic acid (OA) capped nanoparticles embedded into poly-methylmethacrylate (PMMA) matrix for the formation of nanocomposite materials. In particular, simulations of PMMA/OA mixtures at different composition have exhibited the self-assembly of amphiphiles to form separated nanosized domains with different morphologies going from spheres, to tubules up to the formation of continuous planar sheets as the OA composition increases. On the other hand, simulations carried out on nanocomposite systems have shown that NPs do not perturb the observed phase behaviour of PMMA/OA mixtures. In fact, at low OA compositions nanoparticles are confined in the spherical lipid domains to form NP clusters, while at high OA composition NPs appear homogeneously distributed in the continuous lipid domain.     

Flow-induced translocation of polymers through a fluidic channel: a dissipative particle dynamics simulation study

The dynamics of flow-induced translocation of polymers through a fluidic channel has been studied by dissipative particle dynamics (DPD) approach. Unlike implicit solvent models, the many-body energetic and hydrodynamic interactions are preserved naturally by incorporating explicit solvent particles in this approach. The no-slip wall boundary and the adaptive boundary conditions have been implemented in the modified DPD approach to model the hydrodynamic flow within a specific wall structure of fluidic channel and control the particles' density fluctuations. The results show that the average translocation time versus polymer chain length satisfies a power-law scaling of τ ～N(1.152). The conformational changes and translocation dynamics of polymers through the fluidic channel have also been investigated in our simulations, and two different translocation processes, i.e., the single-file and double-folded translocation events, have been observed in detail. These findings may be helpful in understanding the conformational and dynamic behaviors of such polymer and/or DNA molecules during the translocation processes.     

Lateral reorganization of myelin lipid domains by myelin basic protein studied at the air-water interface

It has been speculated that adsorption of myelin basic protein (MBP) to the myelin lipid membrane leads to lateral reorganization of the lipid molecules within the myelin membrane. This hypothesis was tested in this study by surface pressure measurement and fluorescent imaging of a monolayer composed of a myelin lipid mixture. The properties of the lipid monolayer before and after addition of MBP into the subphase were monitored. Upon addition of MBP to the monolayer subphase, the surface pressure rose and significant rearrangement of the lipid domains was observed. These results suggest that binding and partial insertion of MBP into the lipid monolayer led to dramatic rearrangement and morphological changes of the lipid domains. A model of adsorption of MBP to the lipid domains and subsequent domain fusion promoted by minimization of electrostatic repulsion between the domains was proposed to account for the experimental observations. The significance of these results in light of the role of MBP in maintaining the myelin structural integrity is discussed.     

A generalized Born formalism for heterogeneous dielectric environments: application to the implicit modeling of biological membranes

Reliable computer simulations of complex biological environments such as integral membrane proteins with explicit water and lipid molecules remain a challenging task. We propose a modification of the standard generalized Born theory of homogeneous solvent for modeling the heterogeneous dielectric environments such as lipid/water interfaces. Our model allows the representation of biological membranes in the form of multiple layered dielectric regions with dielectric constants that are different from the solute cavity. The proposed new formalism is shown to predict the electrostatic component of solvation free energy with a relative error of 0.17% compared to exact finite-difference solutions of the Poisson equation for a transmembrane helix test system. Molecular dynamics simulations of melittin and bacteriorhodopsin are carried out and performed over 10 ns and 7 ns of simulation time, respectively. The center of melittin along the membrane normal in these stable simulations is in excellent agreement with the relevant experimental data. Simulations of bacteriorhodopsin started from the experimental structure remained stable and in close agreement with experiment. We also examined the free energy profiles of water and amino acid side chain analogs upon membrane insertion. The results with our implicit membrane model agree well with the experimental transfer free energy data from cyclohexane to water as well as explicit solvent simulations of water and selected side chain analogs.     

Heterogeneous molecular distribution in supported multicomponent lipid bilayers

Membrane domains contribute important structural and functional attributes to biological membranes. We describe the heterogeneous nanoscale distribution of lipid molecules within microscale membrane domains in multicomponent lipid bilayers composed of dipalmitoylphosphatidylcholine (DPPC), dilauroylphosphatidylcholine (DLPC), and cholesterol (chol). The lipids were labeled with the fluorescent lipid analogues Bodipy-PC and DiI-C20:0 to identify the distribution of individual membrane components. We used a near-field scanning optical microscope (NSOM) at room temperature to identify the nanoscale structures in the membrane. Simultaneous multicolor NSOM imaging at the emission maxima of the fluorescent analogues revealed a patchy distribution of Bodipy-PC and DiI-C20:0 indicative of phase separations in the bilayer. In a cholesterol-free system (DPPC/DLPC = 1:1), NSOM images proved that the two phosphatidylcholine molecules can coexist in domains at the micrometer level but form nanoscopic patches within the domains; DPPC occurs at the edge of the domains, whereas DLPC is present throughout the domains. In the presence of cholesterol (DPPC/DLPC = 7:3, chol = 18.9%), the two lipid molecules were more miscible but incomplete phase separations also occurred. The average domain sizes were 140-200 nm, well below the resolution capabilities of diffraction-limited light microscopy techniques; the domains were unresolvable by confocal microscopy. Our high-resolution NSOM studies of membrane domain behavior provide a better understanding of complex membrane phase phenomena in multicomponent biological membranes.     

Partitioning of dual-lipidated peptides into membrane microdomains: lipid sorting vs peptide aggregation

The lateral membrane organization and phase behavior of the lipid mixture DMPC(di-C(14))/DSPC(di-C(18))/cholesterol (0-33 mol %) with and without an incorporated fluorescence-labeled palmitoyl/farnesyl dual-lipidated peptide, BODIPY-Gly-Cys(Pal)-Met-Gly-Leu-Pro-Cys(Far)-OMe, which represents a membrane recognition model system for Ras proteins, was studied by two-photon excitation fluorescence microscopy. Measurements were performed on giant unilamellar vesicles (GUVs) over a large temperature range, ranging from 30 to 80 degrees C to cover different lipid phase states (all-gel, fluid/gel, liquid-ordered, all-fluid). At temperatures where the fluid-gel coexistence region of the pure binary phospholipid system occurs, large-scale concentration fluctuations appear. Incorporation of cholesterol levels up to 33 mol % leads to a significant increase of conformational order in the membrane system and a reduction of large domain structures. Adding the peptide leads to dramatic changes in the lateral organization of the membrane. With cholesterol present, a phase separation is induced by a lipid sorting mechanism owing to the high affinity of the lipidated peptide to a fluid, DMPC-rich environment. This phase separation leads to the formation of peptide-containing domains with high fluorescence intensity that become progressively smaller with decreasing temperature. As a result, the local concentration of the peptide increases steadily within the confines of the shrinking domains. At the lowest temperatures, where the acyl-chain order parameter of the membrane has already drastically increased and the membrane achieves a liquid-ordered character, an efficient lipid sorting mechanism is no longer supported and aggregation of the peptide into small clusters prevails. We can conclude that palmitoyl/farnesyl dual-lipidated peptides do not associate with liquid-ordered or gel-like domains in phase-separated bilayer membranes. In particular, the study shows the interesting ability of the peptide to induce formation of fluid microdomains at physiologically relevant cholesterol concentrations, and this effect very much depends on the concentration of fluid vs ordered lipid molecules.     

A second generation mesoscopic lipid bilayer model: connections to field-theory descriptions of membranes and nonlocal hydrodynamics

A new mesoscopic membrane model is developed in order to examine long-wavelength structural and dynamical membrane phenomena. Two different explicit mesoscopic solvent models are employed. The first mesoscopic solvent is denoted the big liquid oscillating blob system, which is parametrized to model water at a coarse-grained level and is motivated by a Langevin-like approach; the resulting membrane dynamics predict a solvent viscosity dependence consistent with the known viscosity of water. The second mesoscopic solvent is a Weeks-Chandler-Anderson model. Here, it is found that the correct mesoscopic hydrodynamic scaling of the membrane undulation dynamics is still preserved, although accelerated. When the behavior of the two membranes in close proximity to one another is examined, very little correlated motion is observed. However, the theoretically predicted scaling of the entropic undulation energy is confirmed, demonstrating that the entropic interaction between two membranes becomes increasingly repulsive with decreasing separation.     

Floret-shaped solid domains on giant fluid lipid vesicles induced by pH

Lateral lipid phase separation of titratable PS or PA lipids and their assembly in domains induced by changes in pH are significant in liposome-based drug delivery: environmentally responsive lipid heterogeneities can be tuned to alter collective membrane properties such as permeability (altering drug release) and surface topography (altering drug carrier reactivity) impacting, therefore, the therapeutic outcomes. At the micrometer scale fluorescence microscopy on giant unilamellar fluid vesicles (GUVs) shows that lowering pH (from 7.0 to 5.0) promotes condensation of titratable PS or PA lipids into beautiful floret-shaped domains in which lipids are tightly packed via hydrogen-bonding and van der Waals interactions. The order of lipid packing within domains increases radially toward the domain center. Lowering pH enhances the lipid packing order, and at pH 5.0 domains appear to be entirely in the solid (gel) phase. Domains phenomenologically comprise a circular "core" cap beyond which interfacial instabilities emerge resembling leaf-like stripes. At pH 5.0 stripes are of almost vanishing Gaussian curvature independent of GUVs' preparation path and in agreement with a general condensation mechanism. Increasing incompressibility of domains is strongly correlated with a larger number of thinner stripes per domain and increasing relative rigidity of domains with smaller core cap areas. Line tension drives domain ripening; however, the final domain shape is a result of enhanced incompressibility and rigidity maximized by domain coupling across the bilayer. Introduction of a transmembrane osmotic gradient (hyperosmotic on the outer lipid leaflet) allows the domain condensation process to reach its maximum extent which, however, is limited by the minimal expansivity of the continuous fluid membrane.     

Local concentration of gel phase domains in supported lipid bilayers under light irradiation in binary mixture of phospholipids doped with dyes for photoinduced activation

Recently, lipid bilayers supported on solid substrates are considered to offer potential as biological devices utilizing biological membranes and membrane proteins. In particular, artificially patterned supported bilayers hold great promise for the development of biological devices. In this study, we show control of the formation and location of phase-separated domain structures by light irradiation for gel phase and liquid-crystalline phase separation structures in a DMPC-DOPC binary lipid bilayer tagged with dye molecules on SiO2/Si substrates. Upon light irradiation, the gel phase domain structures disappeared from the phase-separated bilayers. This disappearance indicates that the light irradiation causes a local increase in the temperature of the lipid bilayer. In this disappearance phenomenon, the photoinduced activation of dye lipids, e.g. fluorescent lipids, is considered to play an important role, since the same phenomenon does not occur in lipid bilayers that have a low concentration of dye lipids. Thus, the local increase in temperature is propagated by light absorption of the dye lipid and subsequent photoinduced activation of nonradiative molecular vibrations. Subsequent interruption of the photoinduced activation for molecular motion allowed the gel phase domain structures to precipitate and grow again. Moreover, the domain area fraction remaining after the photoinduced activation was higher than that before the photoinduced activation. This result indicates that the local increase in temperature propagated by dye-excitation enhances formation of the gel phase domains. By utilizing this phenomenon, we could preferentially induce formation of domain structures within the light-irradiated regions. This technique could be the basis for a new patterning technique based on domain structures. Moreover, these domain structure patterns can be eliminated by increasing the temperature, allowing rewritable patterning.     

Diffusion of nanoparticles in dense fluids

The diffusion process of a single spherical nanoparticle immersed in a fluid solvent is studied by molecular dynamics simulations. When the nanoparticle mass stays constant, it is shown that, at short times, the decay of the nanoparticle velocity autocorrelation function is strongly modified when the particle diameter increases. It is also shown that, at large times, the characteristic algebraic decay induced by the hydrodynamic correlations between the solvated particle and the solvent presents a scaling behavior depending on the particle diameter.     

Complementary matching in domain formation within lipid bilayers

Domain structure and formation in lipid bilayers are investigated by molecular dynamics simulations using a coarse-grained lipid model. The lipid bilayers consist of two lipid types that are identical except for tail length. At a temperature intermediate to the two melting temperatures of the constituent lipid types, gel domains spontaneously form from an initial random structure. The simulations reveal that the gel domains consist of both lipid types in a complementary match. If a long lipid is in the top monolayer, then a short lipid is underneath and vice versa. The gel domains have a larger thickness than the surrounding liquid phase. The thickness of the gel domains is close to that of the pure long lipid gel phase bilayers. However, since in the mixed gel domains the lipids are not tilted and in the pure gel phase the lipids are tilted, the two thicknesses are similar, and the underlying structure is therefore not distinguishable solely by thickness measurements.     

Structural characterization of the voltage-sensor domain and voltage-gated K+-channel proteins vectorially oriented within a single bilayer membrane at the solid/vapor and solid/liquid interfaces via neutron interferometry

The voltage-sensor domain (VSD) is a modular four-helix bundle component that confers voltage sensitivity to voltage-gated cation channels in biological membranes. Despite extensive biophysical studies and the recent availability of X-ray crystal structures for a few voltage-gated potassium (Kv) channels and a voltage-gate sodium (Nav) channel, a complete understanding of the cooperative mechanism of electromechanical coupling, interconverting the closed-to-open states (i.e., nonconducting to cation conducting) remains undetermined. Moreover, the function of these domains is highly dependent on the physical-chemical properties of the surrounding lipid membrane environment. The basis for this work was provided by a recent structural study of the VSD from a prokaryotic Kv-channel vectorially oriented within a single phospholipid (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)) membrane investigated by X-ray interferometry at the solid/moist He (or solid/vapor) and solid/liquid interfaces, thus achieving partial to full hydration, respectively (Gupta et al. Phys. Rev. E2011, 84, 031911-1-15). Here, we utilize neutron interferometry to characterize this system in substantially greater structural detail at the submolecular level, due to its inherent advantages arising from solvent contrast variation coupled with the deuteration of selected submolecular membrane components, especially important for the membrane at the solid/liquid interface. We demonstrate the unique vectorial orientation of the VSD and the retention of its molecular conformation manifest in the asymmetric profile structure of the protein within the profile structure of this single bilayer membrane system. We definitively characterize the asymmetric phospholipid bilayer solvating the lateral surfaces of the VSD protein within the membrane. The profile structures of both the VSD protein and phospholipid bilayer depend upon the hydration state of the membrane. We also determine the distribution of water and exchangeable hydrogen throughout the profile structure of both the VSD itself and the VSD:POPC membrane. These two experimentally determined water and exchangeable hydrogen distribution profiles are in good agreement with molecular dynamics simulations of the VSD protein vectorially oriented within a fully hydrated POPC bilayer membrane, supporting the existence of the VSD's water pore. This approach was extended to the full-length Kv-channel (KvAP) at a solid/liquid interface, providing the separate profile structures of the KvAP protein and the POPC bilayer within the reconstituted KvAP:POPC membrane.