Acyl chain order parameter profiles in phospholipid bilayers: computation from molecular dynamics simulations and comparison with 2H NMR experiments

Localization and Ordering of Lipids Around Aquaporin-0: Protein and Lipid Mobility Effects

Hydrophobic matching, lipid sorting, and protein oligomerization are key principles by which lipids and proteins organize in biological membranes. The Aquaporin-0 channel (AQP0), solved by electron crystallography (EC) at cryogenic temperatures, is one of the few protein-lipid complexes of which the structure is available in atomic detail. EC and room-temperature molecular dynamics (MD) of dimyristoylglycerophosphocholine (DMPC) annular lipids around AQP0 show similarities, however, crystal-packing and temperature might affect the protein surface or the lipids distribution. To understand the role of temperature, lipid phase, and protein mobility in the localization and ordering of AQP0-lipids, we used MD simulations of an AQP0-DMPC bilayer system. Simulations were performed at physiological and at DMPC gel-phase temperatures. To decouple the protein and lipid mobility effects, we induced gel-phase in the lipids or restrained the protein. We monitored the lipid ordering effects around the protein. Reducing the system temperature or inducing lipid gel-phase had a marginal effect on the annular lipid localization. However, restraining the protein mobility increased the annular lipid localization around the whole AQP0 surface, resembling EC. The distribution of the inter-phosphate and hydrophobic thicknesses showed that stretching of the DMPC annular layer around AQP0 surface is the mechanism that compensates the hydrophobic mismatch in this system. The distribution of the local area-per-lipid and the acyl-chain order parameters showed particular fluid- and gel-like areas that involved several lipid layers. These areas were in contact with the surfaces of higher and lower protein mobility, respectively. We conclude that the AQP0 surfaces induce specific fluid- and gel-phase prone areas. The presence of these areas might guide the AQP0 lipid sorting interactions with other membrane components, and is compatible with the squared array oligomerization of AQP0 tetramers separated by a layer of annular lipids.

Biomimetic Phospholipid Membrane Organization on Graphene and Graphene Oxide Surfaces: A Molecular Dynamics Simulation Study

Supported phospholipid membrane patches stabilized on graphene surfaces have shown potential in sensor device functionalization, including biosensors and biocatalysis. Lipid dip-pen nanolithography (L-DPN) is a method useful in generating supported membrane structures that maintain lipid functionality, such as exhibiting specific interactions with protein molecules. Here, we have integrated L-DPN, atomic force microscopy, and coarse-grained molecular dynamics simulation methods to characterize the molecular properties of supported lipid membranes (SLMs) on graphene and graphene oxide supports. We observed substantial differences in the topologies of the stabilized lipid structures depending on the nature of the surface (polar graphene oxide vs nonpolar graphene). Furthermore, the addition of water to SLM systems resulted in large-scale reorganization of the lipid structures, with measurable effects on lipid lateral mobility within the supported membranes. We also observed reduced lipid ordering within the supported structures relative to free-standing lipid bilayers, attributed to the strong hydrophobic interactions between the lipids and support. Together, our results provide insight into the molecular effects of graphene and graphene oxide surfaces on lipid bilayer membranes. This will be important in the design of these surfaces for applications such as biosensor devices.

Detailed Structural Analysis of a Self-Assembled Vesicular Amphiphilic NCN-Pincer Palladium Complex by Using Wide-Angle X-Ray Scattering and Molecular Dynamics Calculations

Wide-angle X-ray scattering experiments and all-atomistic molecular dynamics calculations were performed to elucidate the detailed structure of bilayer vesicles constructed by self-assembly of an amphiphilic palladium NCN-pincer complex. We found an excellent agreement between the experimental and calculated X-ray spectra, and between the membrane thickness determined from a TEM image and that calculated from an electron-density profile, which indicated that the calculated structure was highly reliable. The analysis of the simulated bilayer structure showed that in general the membrane was softer than other phospholipid bilayer membranes. In this bilayer assemblage, the degree of alignment of complex molecules in the bilayer membrane was quite low. An analysis of the electron-density profile shows that the bilayer assemblage contains a space through which organic molecules can exit. Furthermore, the catalytically active center is near this space and is easily accessible by organic molecules, which permits the bilayer membrane to act as a nanoreactor. The free energy of permeation of water through the bilayer membrane of the amphiphilic complex was 12 kJ mol^-1 , which is much lower than that for phospholipid bilayer membranes in general. Organic molecules are expected to pass though the bilayer membrane. The self-assembled vesicles were shown to be catalytically active in a Miyaura-Michael reaction in water.

Molecular Dynamics of POPC Phospholipid Bilayers through the Gel to Fluid Phase Transition: An Incoherent Quasi-Elastic Neutron Scattering Study

The microscopic dynamics for the gel and liquid-crystalline phase of highly aligned D2O-hydrated bilayers of 1-palmitoyl-oleoyl-sn-glycero-phosphocholine (POPC) were investigated in the temperature range from 248 to 273 K by using incoherent quasi-elastic neutrons scattering (QENS). We develop a model for describing the molecular motions of the liquid phase occurring in the 0.3 to 350 ps time range. Accordingly, the complex dynamics of hydrogen are described in terms of simple dynamical processes involving different parts of the phospholipid chain. The analysis of the data evidences the existence of three different motions: the fast motion of hydrogen vibrating around the carbon atoms, the intermediate motion of carbon atoms in the acyl chains, and the slower translational motion of the entire phospholipid molecule. The influence of the temperature on these dynamical processes is investigated. In particular, by going from gel to liquid-crystalline phase, we reveal an increase of the segmental motion mainly affecting the terminal part of the acyl chains and a change of the diffusional dynamics from a localized rattling-like motion to a confined diffusion.

Inclusion of lipopeptides into the DMPC lipid bilayers prevents Aβ peptide insertion

Lipopeptides prevent penetration of Alzheimer's Aβ peptides into lipid bilayers.

Highly Hydrophilic Segments Attached to Hydrophobic Peptides Translocate Rapidly across Membranes

Hydrophilic segments attached to transmembrane helices often cross membranes. In an increasing number of cases, it has become apparent that this occurs in a biologically relevant post-translational event. In this study, we investigate whether juxta-membrane (JM) hydrophilic sequences attached to hydrophobic helices are able to rapidly cross lipid bilayers via their ability or inability to block hydrophobic helix interconversion between a transmembrane (TM) and non-TM membrane-associated state. Interconversion was triggered by changing the protonation state of an Asp residue in the hydrophobic core of the peptides, and peptide configuration was monitored by the fluorescence of a Trp residue at the center of the hydrophobic sequence. In POPC vesicles, conversion of the TM to non-TM state at high pH and the non-TM to TM state at low pH was rapid (seconds or less) for KK, KKNN, and the KKNNNNNN flanking sequences on both N- and C-termini and the KLFAGHQ sequence that flanks the spontaneously TM-inserting 3A protein of polio virus. In vesicles composed of 6:4 (mol/mol) POPC/cholesterol, interconversion was still rapid, with the exception of the peptide flanked by KKNNNNNN sequences, for which the half time of interconversion slowed to minutes. This behavior suggests that, at least in membranes with low levels of cholesterol, movement of hydrophilic JM segments (and analogous hydrophobic loops in multipass TM proteins) across membranes may be more facile than previously thought. This may have important biological implications.

Analytical calculation of the lipid bilayer bending modulus

Bending and Gaussian moduli of a homogenious single-component lipid bilayer are calculated analytically using microscopic model of the lipid hydrocarbon chains. The approach allows for thermodynamic averaging over different chains conformations. Each chain is modeled as a flexible string with finite bending rigidity and an incompressible cross-section area. The interchain steric repulsion is accounted for self-consistently determined single-chain confining parabolic potential. The model provides a simple analytical expression for the membrane bending modulus, which falls within a range of experimental values. An observed dependence of the modulus on hydrocarbon chain length is also reproduced. Correspondence between our microscopic model and the membrane theory of elasticity is established.

Algorithm for Designing Nanoscale Supramolecular Therapeutics with Increased Anticancer Efficacy

In the chemical world, evolution is mirrored in the origin of nanoscale supramolecular structures from molecular subunits. The complexity of function acquired in a supramolecular system over a molecular subunit can be harnessed in the treatment of cancer. However, the design of supramolecular nanostructures is hindered by a limited atomistic level understanding of interactions between building blocks. Here, we report the development of a computational algorithm, which we term Volvox after the first multicellular organism, that sequentially integrates quantum mechanical energy-state- and force-field-based models with large-scale all-atomistic explicit water molecular dynamics simulations to design stable nanoscale lipidic supramolecular structures. In one example, we demonstrate that Volvox enables the design of a nanoscale taxane supramolecular therapeutic. In another example, we demonstrate that Volvox can be extended to optimizing the ratio of excipients to form a stable nanoscale supramolecular therapeutic. The nanoscale taxane supramolecular therapeutic exerts greater antitumor efficacy than a clinically used taxane in vivo. Volvox can emerge as a powerful tool in the design of nanoscale supramolecular therapeutics for effective treatment of cancer.

Naratriptan aggregation in lipid bilayers: perspectives from molecular dynamics simulations

In order to understand the interaction between naratriptan and a fully hydrated bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidyl-choline (POPC), we carried out molecular dynamics simulations. The simulations were performed considering neutral and protonated ionization states, starting from different initial conditions. At physiological pH, the protonated state of naratriptan is predominant. It is expected that neutral compounds could have larger membrane partition than charged compounds. However, for the specific case of triptans, it is difficult to study neutral species in membranes experimentally, making computer simulations an interesting tool. When the naratriptan molecules were originally placed in water, they partitioned between the bilayer/water interface and water phase, as has been described for similar compounds. From this condition, the drugs displayed low access to the hydrophobic environment, with no significant effects on bilayer organization. The molecules anchored in the interface, due mainly to the barrier function of the polar and oriented lipid heads. On the other hand, when placed inside the bilayer, both neutral and protonated naratriptan showed self-aggregation in the lipid tail environment. In particular, the protonated species exhibited a pore-like structure, dragging water through this environment. Graphical Abstract Different behaviour of Naratriptan and Sumatriptan, when the drugs were originally placed in the lipid core.

How the spontaneous insertion of amphiphilic imidazolium-based cations changes biological membranes: a molecular simulation study

The insertion of 1-octyl-3-methylimidazolium cations (OMIM(+)) from a diluted aqueous ionic liquid (IL) solution into a model of a bacterial cell membrane is investigated. Subsequently, the mutual interactions of cations inside the membrane and their combined effect on membrane properties are derived. The ionic liquid solution and the membrane model are simulated using molecular dynamics in combination with empirical force fields. A high propensity of OMIM(+) for membrane insertion is observed, with a cation concentration at equilibrium inside the membrane 47 times larger than in the solvent. Once inserted, cations exhibit a weak effective attraction inside the membrane at a distance of 1.3 nm. At this free energy minimum, negatively charged phosphates of the phospholipids are sandwiched between two OMIM(+) to form energetically favorable OMIM(+)-phosphate-OMIM(+) types of coordination. The cation-cation association free energy is 5.9 kJ mol(-1), whereas the activation barrier for dissociation is 10.1 kJ mol(-1). Subsequently, OMIM(+) are inserted into the leaflet of the membrane bilayer that represents the extracellular side. The cations are evenly distributed with mutual cation distances according to the found optimum distance of 1.3 nm. Because of the short length of the cation alkyl chains compared to lipid fatty acids, voids are generated in the hydrophobic core of the membrane. These voids disorder the fatty acids, because they enable fatty acids to curl into these empty spaces and also cause a thinning of the membrane by 0.6 nm. Additionally, the membrane density increases at its center. The presence of OMIM(+) in the membrane facilitates the permeation of small molecules such as ammonia through the membrane, which is chosen as a model case for small polar solutes. The permeability coefficient of the membrane with respect to ammonia increases substantially by a factor of seven. This increase is caused by a reduction of the involved free energy barriers, which is effected by the cations through the thinning of the membrane and favorable interactions of the delocalized OMIM(+) charge with ammonia inside the membrane. Overall, the results indicate the antimicrobial effect of amphiphilic imidazolium-based cations that are found in various common ILs. This effect is caused by an alteration of the permeability of the bacterial membrane and other property changes.

Rationally Designed 2-in-1 Nanoparticles Can Overcome Adaptive Resistance in Cancer

The development of resistance is the major cause of mortality in cancer. Combination chemotherapy is used clinically to reduce the probability of evolution of resistance. A similar trend toward the use of combinations of drugs is also emerging in the application of cancer nanomedicine. However, should a combination of two drugs be delivered from a single nanoparticle or should they be delivered in two different nanoparticles for maximal efficacy? We explored these questions in the context of adaptive resistance, which emerges as a phenotypic response of cancer cells to chemotherapy. We studied the phenotypic dynamics of breast cancer cells under cytotoxic chemotherapeutic stress and analyzed the data using a phenomenological mathematical model. We demonstrate that cancer cells can develop adaptive resistance by entering into a predetermined transitional trajectory that leads to phenocopies of inherently chemoresistant cancer cells. Disrupting this deterministic program requires a unique combination of inhibitors and cytotoxic agents. Using two such combinations, we demonstrate that a 2-in-1 nanomedicine can induce greater antitumor efficacy by ensuring that the origins of adaptive resistance are terminated by deterministic spatially constrained delivery of both drugs to the target cells. In contrast, a combination of free-form drugs or two nanoparticles, each carrying a single payload, is less effective, arising from a stochastic distribution to cells. These findings suggest that 2-in-1 nanomedicines could emerge as an important strategy for targeting adaptive resistance, resulting in increased antitumor efficacy.

Pair interactions among ternary DPPC/POPC/cholesterol mixtures in liquid-ordered and liquid-disordered phases

Saturated phospholipids, unsaturated phospholipids, and cholesterol are essential components of cell membranes, making the understanding of their mutual interactions of great significance. We have performed microsecond molecular dynamics simulations on the ternary mixtures of DPPC/POPC/cholesterol to systematically examine lipid-lipid and cholesterol-lipid interactions in the liquid-ordered and the liquid-disordered phases. The results show that there exists a competition between the tighter packing of cholesterol-lipid and the looser packing of lipid-lipid as the membrane changes from the liquid-disordered phase to the liquid-ordered phase. Depending on the lipid saturation, the favor of lipid-lipid interactions is in the order of saturated-saturated > monounsaturated-monounsaturated > saturated-monounsaturated. Cholesterol-saturated lipid interactions are more favorable than cholesterol-monounsaturated lipid ones. The results are consistent with the push-pull forces derived from experiments and give general insights into the interactions among membrane components.

Effect of lipid composition and amino acid sequence upon transmembrane peptide-accelerated lipid transleaflet diffusion (flip-flop)

We examined how hydrophobic peptide-accelerated transleaflet lipid movement (flip-flop) was affected by peptide sequence and vesicle composition and properties. A peptide with a completely hydrophobic sequence had little if any effect upon flip-flop. While peptides with a somewhat less hydrophobic sequence accelerated flip-flop, the half-time remained slow (hours) with substantial (0.5mol%) peptide in the membranes. It appears that peptide-accelerated lipid flip-flop involves a rare event that may reflect a rare state of the peptide or lipid bilayer. There was no simple relationship between peptide overall hydrophobicity and flip-flop. In addition, flip-flop was not closely linked to whether the peptides were in a transmembrane or non-transmembrane (interfacial) inserted state. Flip-flop was also not associated with peptide-induced pore formation. We found that peptide-accelerated flip-flop is initially faster in small (highly curved) unilamellar vesicles relative to that in large unilamellar vesicles. Peptide-accelerated flip-flop was also affected by lipid composition, being slowed in vesicles with thick bilayers or those containing 30% cholesterol. Interestingly, these factors also slow spontaneous lipid flip-flop in the absence of peptide. Combined with previous studies, the results are most consistent with acceleration of lipid flip-flop by peptide-induced thinning of bilayer width.

Photopolymerization of Dienoyl Lipids Creates Planar Supported Poly(lipid) Membranes with Retained Fluidity

Polymerization of substrate-supported bilayers composed of dienoylphosphatidylcholine (PC) lipids is known to greatly enhance their chemical and mechanical stability; however, the effects of polymerization on membrane fluidity have not been investigated. Here planar supported lipid bilayers (PSLBs) composed of dienoyl PCs on glass substrates were examined to assess the degree to which UV-initiated polymerization affects lateral lipid mobility. Fluorescence recovery after photobleaching (FRAP) was used to measure the diffusion coefficients (D) and mobile fractions of rhodamine-DOPE in unpolymerized and polymerized PSLBs composed of bis-sorbyl phosphatidylcholine (bis-SorbPC), mono-sorbyl-phosphatidylcholine (mono-SorbPC), bis-dienoyl-phosphatidylcholine (bis-DenPC), and mono-dienoyl phosphatidylcholine (mono-DenPC). Polymerization was performed in both the Lα and Lβ phase for each lipid. In all cases, polymerization reduced membrane fluidity; however, measurable lateral diffusion was retained which is attributed to a low degree of polymerization. The D values for sorbyl lipids were less than those of the denoyl lipids; this may be a consequence of the distal location of polymerizable group in the sorbyl lipids which may facilitate interleaflet bonding. The D values measured after polymerization were 0.1-0.8 of those measured before polymerization, a range that corresponds to fluidity intermediate between that of a Lα phase and a Lβ phase. This D range is comparable to ratios of D values reported for liquid-disordered (Ld) and liquid-ordered (Lo) lipid phases and indicates that the effect of UV polymerization on lateral diffusion in a dienoyl PSLB is similar to the transition from a Ld phase to a Lo phase. The partial retention of fluidity in UV-polymerized PSLBs, their enhanced stability, and the activity of incorporated transmembrane proteins and peptides is discussed.

Force fields for simulating the interaction of surfaces with biological molecules

The interaction of biomolecules with solid interfaces is of fundamental importance to several emerging biotechnologies such as medical implants, anti-fouling coatings and novel diagnostic devices. Many of these technologies rely on the binding of peptides to a solid surface, but a full understanding of the mechanism of binding, as well as the effect on the conformation of adsorbed peptides, is beyond the resolution of current experimental techniques. Nanoscale simulations using molecular mechanics offer potential insights into these processes. However, most models at this scale have been developed for aqueous peptide and protein simulation, and there are no proven models for describing biointerfaces. In this review, we detail the current research towards developing a non-polarizable molecular model for peptide-surface interactions, with a particular focus on fitting the model parameters as well as validation by choice of appropriate experimental data.

Living on the edge: Simulations of bacterial outer-membrane proteins

Gram-negative bacteria are distinguished in part by a second, outer membrane surrounding them. This membrane is distinct from others, possessing an outer leaflet composed not of typical phospholipids but rather large, highly charged molecules known as lipopolysaccharides. Therefore, modeling the structure and dynamics of proteins embedded in the outer membrane requires careful consideration of their native environment. In this review, we examine how simulations of such outer-membrane proteins have evolved over the last two decades, culminating most recently in detailed, highly accurate atomistic models of the outer membrane. We also draw attention to how the simulations have coupled with experiments to produce novel insights unattainable through a single approach. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Atomistic resolution structure and dynamics of lipid bilayers in simulations and experiments

Accurate details on the sampled atomistic resolution structures of lipid bilayers can be experimentally obtained by measuring C-H bond order parameters, spin relaxation rates and scattering form factors. These parameters can be also directly calculated from the classical atomistic resolution molecular dynamics simulations (MD) and compared to the experimentally achieved results. This comparison measures the simulation model quality with respect to 'reality'. If agreement is sufficient, the simulation model gives an atomistic structural interpretation of the acquired experimental data. Significant advance of MD models is made by jointly interpreting different experiments using the same structural model. Here we focus on phosphatidylcholine lipid bilayers, which out of all model membranes have been studied mostly by experiments and simulations, leading to the largest available dataset. From the applied comparisons we conclude that the acyl chain region structure and rotational dynamics are generally well described in simulation models. Also changes with temperature, dehydration and cholesterol concentration are qualitatively correctly reproduced. However, the quality of the underlying atomistic resolution structural changes is uncertain. Even worse, when focusing on the lipid bilayer properties at the interfacial region, e.g. glycerol backbone and choline structures, and cation binding, many simulation models produce an inaccurate description of experimental data. Thus extreme care must be applied when simulations are applied to understand phenomena where the interfacial region plays a significant role. This work is done by the NMRlipids Open Collaboration project running at https://nmrlipids.blogspot.fi and https://github.com/NMRLipids. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Force Field Development for Lipid Membrane Simulations

With the rapid development of computer power and wide availability of modelling software computer simulations of realistic models of lipid membranes, including their interactions with various molecular species, polypeptides and membrane proteins have become feasible for many research groups. The crucial issue of the reliability of such simulations is the quality of the force field, and many efforts, especially in the latest several years, have been devoted to parametrization and optimization of the force fields for biomembrane modelling. In this review, we give account of the recent development in this area, covering different classes of force fields, principles of the force field parametrization, comparison of the force fields, and their experimental validation. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Interactions of the N-terminal domain of human islet amyloid polypeptide with lipid membranes: the effect of cholesterol

Cholesterol facilitates the insertion and aggregation of hIAPP_1–19 in membrane and the CARC motif mediates the peptide–cholesterol interaction.

Effect of shape and bending modulus on the properties of nematic lyotropic liquid crystals

Variation in the structure of the molecular aggregate associated with the increase of the TTAC concentration in the liquid crystal.

Binding Orientations and Lipid Interactions of Human Amylin at Zwitterionic and Anionic Lipid Bilayers

Increasing evidence suggests that the interaction of human islet amyloid polypeptide (hIAPP) with lipids may facilitate hIAPP aggregation and cause the death of pancreatic islet β-cells. However, the detailed hIAPP-membrane interactions and the influences of lipid compositions are unclear. In this study, as a first step to understand the mechanism of membrane-mediated hIAPP aggregation, we investigate the binding behaviors of hIAPP monomer at zwitterionic palmitoyloleoyl-phosphatidylcholine (POPC) bilayer by performing atomistic molecular dynamics simulations. The results are compared with those of hIAPP at anionic palmitoyloleoyl-phosphatidylglycerol (POPG) bilayers. We find that the adsorption of hIAPP to POPC bilayer is mainly initiated from the C-terminal region and the peptide adopts a helical structure with multiple binding orientations, while the adsorption to POPG bilayer is mostly initiated from the N-terminal region and hIAPP displays one preferential binding orientation, with its hydrophobic residues exposed to water. hIAPP monomer inserts into POPC lipid bilayers more readily than into POPG bilayers. Peptide-lipid interaction analyses show that the different binding features of hIAPP at POPC and POPG bilayers are attributed to different magnitudes of electrostatic and hydrogen-bonding interactions with lipids. This study provides mechanistic insights into the different interaction behaviors of hIAPP with zwitterionic and anionic lipid bilayers.

Alteration of lipid membrane structure and dynamics by diacylglycerols with unsaturated chains

Diacylglycerols (DAGs) with unsaturated acyl chains play many important roles in biomembranes, such as a second messenger and activator for protein kinase C. In this study, three DAGs of distinctly different chain unsaturations (i.e. di16:0DAG (DPG), 16:0-18:1DAG (POG), and di18:1DAG (DOG)) are studied using atomistic MD simulation to compare their roles in the structure and dynamics of 16:0-18:1phosphatidylcholine (POPC) membranes. All three DAGs are able to produce the so-called 'condensing effect' in POPC membranes: decreasing area-per-lipid, and increasing acyl chain order and bilayer thickness. Our visual and quantitative analyses clearly show that DAG with unsaturated chains induce larger spacing between POPC headgroups, compared with DAG with saturated chains; this particular effect has long been hypothesized to be crucial for activating enzymes and receptors in cell membranes. DAGs with unsaturated chains are also located closer to the bilayer/aqueous interface than DPG and are more effective in slowing down lateral diffusion of molecules. We show that DAG molecules seek the "umbrella coverage" from neighboring phospholipid headgroups - similar to cholesterol. Unlike cholesterol, DAGs also hide their chains from water by laterally inserting their chains into the surrounding. Thus, acyl chains of DAG are more spread and disordered than those of PC due to the insertion. By calculating the potential of mean force (PMF) for POPC in POPC/DAG bilayers, we found that all three DAGs can significantly increase the free energy barrier for POPC to flip-flop, but only DAGs with unsaturated chains can additionally increase the free energy of POPC desorption.

Effects of ion interactions with a cholesterol-rich bilayer

Previous molecular dynamics (MD) simulations of ion-lipid interactions have focused on pure phospholipid bilayers. Many functional microdomains in membranes have a complex composition of cholesterol and phospholipids. Here, we reveal the distinctiveness of the interactions and the effects of the ions on a cholesterol-rich bilayer by performing MD simulations of a cholesterol-rich bilayer with a Na(+)/K(+) mixture or a Na(+)/K(+)/Ca(2+)/Mg(2+) mixture. The simulations reveal that Ca(2+) maintains its dominant role in the interaction with the cholesterol-rich bilayer, but the binding affinity of Mg(2+) to the cholesterol-rich bilayer is even weaker than the affinities of Na(+) and K(+), whereas its interaction with pure phospholipid bilayers is strong and is only slightly weaker than that of Ca(2+). Additionally, it was found that the presence of additional divalent cations induces the headgroups of phospholipids to be more perpendicular to the membrane surface, reducing the lateral movement of lipids and slightly altering the ordering and packing of the cholesterol-rich bilayer, different from divalent cations, which strongly influence that ordering and packing of pure phospholipid bilayers. Therefore, this study indicates that cholesterol in the membrane could affect the interactions between membrane and cations. The findings could be helpful in understanding the biological processes relevant to regulation of cations in cholesterol-rich regions.

α-Tocopherol Is Well Designed to Protect Polyunsaturated Phospholipids: MD Simulations

The presumptive function for alpha-tocopherol (αtoc) in membranes is to protect polyunsaturated lipids against oxidation. Although the chemistry of the process is well established, the role played by molecular structure that we address here with atomistic molecular-dynamics simulations remains controversial. The simulations were run in the constant particle NPT ensemble on hydrated lipid bilayers composed of SDPC (1-stearoyl-2-docosahexaenoylphosphatidylcholine, 18:0-22:6PC) and SOPC (1-stearoyl-2-oleoylphosphatidylcholine, 18:0-18:1PC) in the presence of 20 mol % αtoc at 37°C. SDPC with SA (stearic acid) for the sn-1 chain and DHA (docosahexaenoic acid) for the sn-2 chain is representative of polyunsaturated phospholipids, while SOPC with OA (oleic acid) substituted for the sn-2 chain serves as a monounsaturated control. Solid-state (2)H nuclear magnetic resonance and neutron diffraction experiments provide validation. The simulations demonstrate that high disorder enhances the probability that DHA chains at the sn-2 position in SDPC rise up to the bilayer surface, whereby they encounter the chromanol group on αtoc molecules. This behavior is reflected in the van der Waals energy of interaction between αtoc and acyl chains, and illustrated by density maps of distribution for acyl chains around αtoc molecules that were constructed. An ability to more easily penetrate deep into the bilayer is another attribute conferred upon the chromanol group in αtoc by the high disorder possessed by DHA. By examining the trajectory of single molecules, we found that αtoc flip-flops across the SDPC bilayer on a submicrosecond timescale that is an order-of-magnitude greater than in SOPC. Our results reveal mechanisms by which the sacrificial hydroxyl group on the chromanol group can trap lipid peroxyl radicals within the interior and near the surface of a polyunsaturated membrane. At the same time, water-soluble reducing agents that regenerate αtoc can access the chromanol group when it locates at the surface.

Cyclic and Linear Monoterpenes in Phospholipid Membranes: Phase Behavior, Bilayer Structure, and Molecular Dynamics

Monoterpenes are abundant in essential oils extracted from plants. These relatively small and hydrophobic molecules have shown important biological functions, including antimicrobial activity and membrane penetration enhancement. The interaction between the monoterpenes and lipid bilayers is considered important to the understanding of the biological functions of monoterpenes. In this study, we investigated the effect of cyclic and linear monoterpenes on the structure and dynamics of lipids in model membranes. We have studied the ternary system 1,2-dimyristoyl-sn-glycero-3-phosphocholine-monoterpene-water as a model with a focus on dehydrated conditions. By combining complementary techniques, including differential scanning calorimetry, solid-state nuclear magnetic resonance, and small- and wide-angle X-ray scattering, bilayer structure, phase transitions, and lipid molecular dynamics were investigated at different water contents. Monoterpenes cause pronounced melting point depression and phase segregation in lipid bilayers, and the extent of these effects depends on the hydration conditions. The addition of a small amount of thymol to the fluid bilayer (volume fraction of 0.03 in the bilayer) leads to an increased order in the acyl chain close to the bilayer interface. The findings are discussed in relation to biological systems and lipid formulations.

A Parameterization of Cholesterol for Mixed Lipid Bilayer Simulation within the Amber Lipid14 Force Field

The Amber Lipid14 force field is expanded to include cholesterol parameters for all-atom cholesterol and lipid bilayer molecular dynamics simulations. The General Amber and Lipid14 force fields are used as a basis for assigning atom types and basic parameters. A new RESP charge derivation for cholesterol is presented, and tail parameters are adapted from Lipid14 alkane tails. 1,2-Dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers are simulated at a range of cholesterol contents. Experimental bilayer structural properties are compared with bilayer simulations and are found to be in good agreement. With this parameterization, another component of complex membranes is available for molecular dynamics with the Amber Lipid14 force field.

Structure and Dynamics of Phospholipid Nanodiscs from All-Atom and Coarse-Grained Simulations

We investigated structural and dynamical properties of nanodiscs comprising dimyristoylphosphatidylcholine (DMPC) lipids and major scaffold protein MSP1Δ(1-22) from human apolipoprotein A-1 using combined all-atom and coarse-grained (CG) molecular dynamics (MD) simulations. The computational efficiency of the Martini-CG force field enables the spontaneous self-assembly of lipids and scaffold proteins into stable nanodisc structures on time scales up to tens of microseconds. Subsequent all-atom and CG-MD simulations reveal that the lipids in the nanodisc have lower configurational entropy and higher acyl tail order than in a lamellar bilayer phase. These altered average properties arise from rather differential behavior of lipids, depending on their location in the nanodisc. Since the scaffold proteins exert constrictive forces from the outer rim of the disc toward its center, lipids at the center of the nanodisc are highly ordered, whereas annular lipids that are in contact with the MSP proteins are remarkably disordered due to perturbed packing. Although specific differences between all-atom and CG simulations are also evident, the results obtained at both levels of resolution are in overall good agreement with each other and provide atomic level interpretations of recent experiments. Thus, the present study highlights the applicability of multiscale simulation approaches for nanodisc systems and opens the way for future applications, including the study of nanodisc-embedded membrane proteins.

Membrainy: a 'smart', unified membrane analysis tool

BACKGROUND

The study of biological membranes using Molecular Dynamics has become an increasingly popular means by which to investigate the interactions of proteins, peptides and potentials with lipid bilayers. These interactions often result in changes to the properties of the lipids which can modify the behaviour of the membrane. Membrainy is a unified membrane analysis tool that contains a broad spectrum of analytical techniques to enable: measurement of acyl chain order parameters; presentation of 2D surface and thickness maps; determination of lateral and axial headgroup orientations; measurement of bilayer and leaflet thickness; analysis of the annular shell surrounding membrane-embedded objects; quantification of gel percentage; time evolution of the transmembrane voltage; area per lipid calculations; and quantification of lipid mixing/demixing entropy.

RESULTS

Each analytical component within Membrainy has been tested on a variety of lipid bilayer systems and was found to be either comparable to or an improvement upon existing software. For the analytical techniques that have no direct comparable software, our results were confirmed with experimental data.

CONCLUSIONS

Membrainy is a user-friendly, intelligent membrane analysis tool that automatically interprets a variety of input formats and force fields, is compatible with both single and double bilayers, and capable of handling asymmetric bilayers and lipid flip-flopping. Membrainy has been designed for ease of use, requiring no installation or configuration and minimal user-input to operate.

Distinct helix propensities and membrane interactions of human and rat IAPP(1-19) monomers in anionic lipid bilayers

Islet amyloid polypeptide, IAPP or amylin, is a 37-residue peptide hormone coexpressed with insulin by pancreatic β-cells. The aggregation of human IAPP (hIAPP) into amyloid deposits is associated with type II diabetes. Substantial evidence suggests that the interaction of anionic membranes with hIAPP may facilitate peptide aggregation and the N-terminal 1-19 fragment (IAPP(1-19)) plays an important role in peptide-membrane interaction. As a first step to understand how structural differences between human and rat IAPP peptides in membranes may influence the later oligomerization process, we have investigated the structures and orientations of hIAPP(1-19) and the less toxic rIAPP(1-19) (i.e., the H18R mutant of hIAPP(1-19)) monomers in anionic POPG bilayers by performing replica exchange molecular dynamics (REMD) simulations. On the basis of ∼20 μs REMD simulations started from a random coil conformation of the peptide placed in water, we find that unfolded h(r)IAPP(1-19) can insert into the bilayers and the membrane-bound peptide stays mainly at the lipid head-tail interface. hIAPP(1-19) displays higher helix propensity than rIAPP(1-19), especially in the L12-L16 region. The helix is oriented parallel to the bilayer surface and buried in the membrane 0.3-0.8 nm below the phosphorus atoms, consistent with previous electron paramagnetic resonance data. The helical conformation is an amphiphilic helix with its hydrophilic and hydrophobic faces pointing, respectively, to the lipid head and tail regions. The H18R substitution enhances the electrostatic interactions of IAPP(1-19) with the membrane, while it weakens the intrapeptide interactions crucial for helix formation, thus leading to lower helix propensity of rIAPP(1-19). Implications of our simulation results on the membrane-mediated IAPP(1-19) oligomerization are discussed.

Solid state NMR of lipid model membranes

In this chapter we describe the use of solid state nuclear magnetic spectroscopy to study the structure of lyotropic phases and lipid model membranes and show its ability to probe, site specifically, at a sub-Ångstrom resolution. Here, we demonstrate the immense versatility of the technique and its ability to provide information on the different liquid crystalline phases present. A multinuclear for example (31)P, (1)H, and (13)C approach is able to elucidate both the structure and dynamics over a wide variety of timescales. This coupled with a non-perturbing label (2)H is able to provide information such as the order parameters for a wide variety of different liquid phases.

Behaviour of NBD-head group labelled phosphatidylethanolamines in POPC bilayers: a molecular dynamics study

An NBD-diC₁₆PE/POPC bilayer with typical fluorophore inverted-snorkelling configurations, and mass density profiles across the membrane. The wide distribution of the NBD fluorophore lies at the origin of its complex emission kinetics.

Molecular dynamics of dibenz[a,h]anthracene and its metabolite interacting with lung surfactant phospholipid bilayers

Dibenz[a,h]anthracene and its metabolite may form aggregates, which have implications in the clearance process of the lung surfactant phospholipid bilayers.

Dynamics of lipids, cholesterol, and transmembrane α-helices from microsecond molecular dynamics simulations

Extensive all-atom molecular dynamics (∼24 μs total) allowed exploration of configurational space and calculation of lateral diffusion coefficients of the components of a protein-embedded, cholesterol-containing model bilayer. The three model membranes are composed of an ∼50/50 (by mole) dipalmitoylphosphatidylcholine (DPPC)/cholesterol bilayer and contained an α-helical transmembrane protein (HIV-1 gp41 TM). Despite the high concentration of cholesterol, normal Brownian motion was observed and the calculated diffusion coefficients (on the order of 10(-9) cm(2)/s) are consistent with experiments. Diffusion is sensitive to a variety of parameters, and a temperature difference of ∼4 K from thermostat artifacts resulted in 2-10-fold differences in diffusion coefficients and significant differences in lipid order, membrane thickness, and unit cell area. Also, the specific peptide sequence likely underlies the consistently observed faster diffusion in one leaflet. Although the simulations here present molecular dynamics (MD) an order of magnitude longer than those from previous studies, the three systems did not approach ergodicity. The distributions of cholesterol and DPPC around the peptides changed on the microsecond time scale, but not significantly enough to thoroughly explore configurational space. These simulations support conclusions of other recent microsecond MD in that even longer time scales are needed for equilibration of model membranes and simulations of more realistic cellular or viral bilayers.

Long-timescale motions in glycerol-monopalmitate lipid bilayers investigated using molecular dynamics simulation

The occurrence of long-timescale motions in glycerol-1-monopalmitate (GMP) lipid bilayers is investigated based on previously reported 600 ns molecular dynamics simulations of a 2×8×8 GMP bilayer patch in the temperature range 302-338 K, performed at three different hydration levels, or in the presence of the cosolutes methanol or trehalose at three different concentrations. The types of long-timescale motions considered are: (i) the possible phase transitions; (ii) the precession of the relative collective tilt-angle of the two leaflets in the gel phase; (iii) the trans-gauche isomerization of the dihedral angles within the lipid aliphatic tails; and (iv) the flipping of single lipids across the two leaflets. The results provide a picture of GMP bilayers involving a rich spectrum of events occurring on a wide range of timescales, from the 100-ps range isomerization of single dihedral angles, via the 100-ns range of tilt precession motions, to the multi-μs range of phase transitions and lipid-flipping events.