Lipid conformation in model membranes and biological membranes

Calcium ions do not influence the Aβ(25-35) triggered morphological changes of lipid membranes

We have studied the effect of calcium ions (Ca2+) at various concentrations on the structure of lipid vesicles in the presence of amyloid-beta peptide Aβ(25-35). In particular, we have investigated the influence of calcium ions on the formation of recently documented bicelle-like structures (BLSs) emerged as a result of Aβ(25-35) triggered membrane disintegration. First, we have shown by using small-angle X-ray and neutron scattering that peptide molecules rigidify the lipid bilayer of gel phase DPPC unilamellar vesicles (ULVs), while addition of the calcium ions to the system hinders this effect of Aβ(25-35). Secondly, the Aβ(25-35) demonstrates a critical peptide concentration at which the BLSs reorganize from ULVs due to heating and cooling the samples through the lipid main phase transition temperature (Tm). However, addition of calcium ions does not affect noticeably the Aβ-induced formation of BLSs and their structural parameters, though the changes in peptide's secondary structure, e.g. the increased α-helix fraction, has been registered by circular dichroism spectroscopy. Finally, according to 31P nuclear magnetic resonance (NMR) measurements, calcium ions do not affect the lipid-peptide arrangement in BLSs and their ability to align in the magnetic field of NMR spectrometer. The influences of various concentrations of calcium ions on the lipid-peptide interactions may prove biologically important because their local concentrations vary widely in in-vivo conditions. In the present work, calcium ions were investigated as a possible tool aimed at regulating the lipid-peptide interactions that demonstrated the disruptive effect of Aβ(25-35) on lipid membranes.

Sum-frequency vibrational spectroscopy, a tutorial: Applications for the study of lipid membrane structure and dynamics

Planar supported lipid bilayers (PSLBs) are an ideal model for the study of lipid membrane structures and dynamics when using sum-frequency vibrational spectroscopy (SFVS). In this paper, we describe the construction of asymmetric PSLBs and the basic SFVS theory needed to understand and make measurements on these membranes. Several examples are presented, including the determination of phospholipid orientation and measuring phospholipid transmembrane translocation (flip-flop).

Orientational Nanoconjugation with Gold Endows Marked Antimicrobial Potential and Drugability of Ultrashort Dipeptides

Antibiotic resistance is a global threat. Antimicrobial peptides (AMPs) are highly desirable to treat multidrug-resistant pathogen infection. However, few AMPs are clinically available, due to high cost, instability, and poor selectivity. Here, ultrashort AMPs (2-3 residues with an N-terminal cysteine) are designed and assembled as gold nanoparticles. Au-S conjugation and ultrashort size restrict nonspecific reactions and peptide orientation, thus concentrating positively charged residues on the surface. The nanostructured assemblies enormously enhance antimicrobial abilities by 1000-6000-fold and stability. One representative (Au-Cys-Arg-NH2, Au_CR) shows selective antibacterial activity against Staphylococcus aureus with 10 nM minimal inhibitory concentration. Au_CR has comparable or better in vivo antimicrobial potency than vancomycin and methicillin, with low propensity to induce resistance, little side effects, and high stability (17.5 h plasma half-life). Au_CR acts by inducing collapse of membrane potential and rupture of the bacterial membrane. The report provides insights for developing AMP-metal nanohybrids, particularly tethering nonspecific reactions and AMP orientation on the metal surface.

Passive and reversible area regulation of supported lipid bilayers in response to fluid flow

Biological and model membranes are frequently subjected to fluid shear stress. However, membrane mechanical responses to flow remain incompletely described. This is particularly true of membranes supported on a solid substrate, and the influences of membrane composition and substrate roughness on membrane flow responses remain poorly understood. Here, we combine microfluidics, fluorescence microscopy, and neutron reflectivity to explore how supported lipid bilayer patches respond to controlled shear stress. We demonstrate that lipid membranes undergo a significant, passive, and partially reversible increase in membrane area due to flow. We show that these fluctuations in membrane area can be constrained, but not prevented, by increasing substrate roughness. Similar flow-induced changes to membrane structure may contribute to the ability of living cells to sense and respond to flow.

Phospholipid headgroups govern area per lipid and emergent elastic properties of bilayers

Phospholipid bilayers are liquid-crystalline materials whose intermolecular interactions at mesoscopic length scales have key roles in the emergence of membrane physical properties. Here we investigated the combined effects of phospholipid polar headgroups and acyl chains on biophysical functions of membranes with solid-state 2H NMR spectroscopy. We compared the structural and dynamic properties of phosphatidylethanolamine and phosphatidylcholine with perdeuterated acyl chains in the solid-ordered (so) and liquid-disordered (ld) phases. Our analysis of spectral lineshapes of 1,2-diperdeuteriopalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE-d62) and 1,2-diperdeuteriopalmitoyl-sn-glycero-3-phosphocholine (DPPC-d62) in the so (gel) phase indicated an all-trans rotating chain structure for both lipids. Greater segmental order parameters (SCD) were observed in the ld (liquid-crystalline) phase for DPPE-d62 than for DPPC-d62 membranes, while their mixtures had intermediate values irrespective of the deuterated lipid type. Our results suggest the SCD profiles of the acyl chains are governed by methylation of the headgroups and are averaged over the entire system. Variations in the acyl chain molecular dynamics were further investigated by spin-lattice (R1Z) and quadrupolar-order relaxation (R1Q) measurements. The two acyl-perdeuterated lipids showed distinct differences in relaxation behavior as a function of the order parameter. The R1Z rates had a square-law dependence on SCD, implying collective mesoscopic dynamics, with a higher bending rigidity for DPPE-d62 than for DPPC-d62 lipids. Remodeling of lipid average and dynamic properties by methylation of the headgroups thus provides a mechanism to control the actions of peptides and proteins in biomembranes.

Improved chemical and isotopic labeling of biomembranes in Bacillus subtilis by leveraging CRISPRi inhibition of beta-ketoacyl-ACP synthase (fabF)

Assessing the structure of living microbial cell membranes is a challenging analytical goal. The cell membrane is defined by its transverse structure, an approximately 5 nm-thick selectively permeable bilayer that serves many important cellular functions. Compositionally complex, dynamic, and organized in both the transverse and lateral dimensions, understanding the cell membrane structure-and the role that structure plays in cellular function, communication, and environmental sensing is an active scientific effort. Previously, we have devised a novel isotopic labeling approach for membrane lipids to enable direct in vivo structural studies of the cell membrane in the Gram-positive bacterium, Bacillus subtilis, using small-angle neutron scattering. This was accomplished through a genetic inhibition of fatty acid (FA) degradation (ΔfadN) and a chemical inhibition of FA biosynthesis using cerulenin, an irreversible inhibitor of type II fatty acid synthases. Here, we improve upon the previous system by introducing a dCas9/sgRNA-fabF complex that blocks transcription of the essential fabF gene when under xylose induction. This leads to greater sensitivity to cerulenin in the mutant strain (JEBS102) and more robust cell growth when supplementary FAs are introduced to the culture medium. A subtle change in FA uptake is noted when compared to the prior labeling strategy. This is seen in the gas chromatography/mass spectrometry (GC/MS) data as a higher ratio of n16:0 to a15:0, and manifests in an apparent increase in the membrane thickness determined via neutron scattering. This represents an improved method of isotopic labeling for the cell membrane of Bacillus subtilis; enabling improved investigations of cellular uptake and utilization of FAs, cell membrane structure and organization as a phenotypic response to metabolic and environmental changes.

Location of the TEMPO moiety of TEMPO-PC in phosphatidylcholine bilayers is membrane phase-dependent

The (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO) moiety tethered to the headgroup of phosphatidylcholine (PC) lipid is employed in spin labeling electron paramagnetic resonance (EPR) spectroscopy to probe the water dynamics near lipid bilayer interfaces. Due to its amphiphilic character, however, TEMPO spin label could partition between aqueous and lipid phases, and may even be stabilized in the lipid phase. Accurate assessment of the TEMPO-PC configuration in bilayer membranes is essential for correctly interpreting the data from measurements. Here, we carry out all-atom molecular dynamics (MD) simulations of TEMPO-PC probe in single-component lipid bilayers at varying temperatures, using two standard MD force fields. We find that for DPPC membrane whose gel-to-fluid lipid phase transition occurs at 314 K, while the TEMPO spin label is stabilized above the bilayer interface in the gel phase, there is a preferential location of TEMPO below the membrane interface in the fluid phase. For bilayers made of unsaturated lipids, DOPC and POPC, which adopt the fluid phase at ambient temperature, TEMPO is unequivocally stabilized inside the bilayers. Our finding of membrane phase-dependent positioning of TEMPO moiety highlights the importance of assessing the packing order and fluidity of lipids under a given measurement condition.

Bending Modulus of Lipid Membranes from Density Correlation Functions

The bending modulus κ quantifies the elasticity of biological membranes in terms of the free energy cost of increasing the membrane corrugation. Molecular dynamics (MD) simulations provide a powerful approach to quantify κ by analyzing the thermal fluctuations of the lipid bilayer. However, existing methods require the identification and filtering of non-mesoscopic fluctuation modes. State of the art methods rely on identifying a smooth surface to describe the membrane shape. These methods introduce uncertainties in calculating κ since they rely on different criteria to select the relevant fluctuation modes. Here, we present a method to compute κ using molecular simulations. Our approach circumvents the need to define a mesoscopic surface or an orientation field for the lipid tails explicitly. The bending and tilt moduli can be extracted from the analysis of the density correlation function (DCF). The method introduced here builds on the Bedeaux and Weeks (BW) theory for the DCF of fluctuating interfaces and on the coupled undulatory (CU) mode introduced by us in previous work. We test the BW-DCF method by computing the elastic properties of lipid membranes with different system sizes (from 500 to 6000 lipid molecules) and using coarse-grained (for POPC and DPPC lipids) and fully atomistic models (for DPPC). Further, we quantify the impact of cholesterol on the bending modulus of DPPC bilayers. We compare our results with bending moduli obtained with X-ray diffraction data and different computer simulation methods.

Temperature regulation of membrane composition in the Firmicute, Enterococcus faecalis, parallels that of Escherichia coli

Both Enterococcus faecalis and Escherichia coli can undergo abrupt temperature transitions in nature. E. coli changes the composition of its phospholipid acyl chains in response to shifts growth temperature. This is mediated by a naturally temperature sensitive enzyme, FabF (3-ketoacyl-acyl carrier protein synthase II), that elongates the 16 carbon unsaturated acyl chain palmitoleate to the 18 carbon unsaturated acyl chain, cis-vaccenate. FabF is more active at low temperatures resulting in increased incorporation of cis-vaccenoyl acyl chains into the membrane phospholipids. This response to temperature is an intrinsic property of FabF and does not require increased synthesis of the enzyme. We report that the FabF of the very divergent bacterium, E. faecalis, has properties very similar to E. coli FabF and is responsible for changing E. faecalis membrane phospholipid acyl chain composition in response to temperature. Moreover, expression E. faecalis FabF in an E. coli ∆fabF strain restores temperature regulation to the E. coli strain.

Cholesterol Effects on the Physical Properties of Lipid Membranes Viewed by Solid-state NMR Spectroscopy

In this chapter, we review the physical properties of lipid/cholesterol mixtures involving studies of model membranes using solid-state NMR spectroscopy. The approach allows one to quantify the average membrane structure, fluctuations, and elastic deformation upon cholesterol interaction. Emphasis is placed on understanding the membrane structural deformation and emergent fluctuations at an atomistic level. Lineshape measurements using solid-state NMR spectroscopy give equilibrium structural properties, while relaxation time measurements study the molecular dynamics over a wide timescale range. The equilibrium properties of glycerophospholipids, sphingolipids, and their binary and tertiary mixtures with cholesterol are accessible. Nonideal mixing of cholesterol with other lipids explains the occurrence of liquid-ordered domains. The entropic loss upon addition of cholesterol to sphingolipids is less than for glycerophospholipids, and may drive formation of lipid rafts. The functional dependence of ²H NMR spin-lattice relaxation (R _1Z) rates on segmental order parameters (S _CD) for lipid membranes is indicative of emergent viscoelastic properties. Addition of cholesterol shows stiffening of the bilayer relative to the pure lipids and this effect is diminished for lanosterol. Opposite influences of cholesterol and detergents on collective dynamics and elasticity at an atomistic scale can potentially affect lipid raft formation in cellular membranes.

Lipid bilayer position and orientation of novel carprofens, modulators of γ-secretase in Alzheimer's disease

γ-Secretase is an integral membrane protein complex and is involved in the cleavage of the amyloid precursor protein APP to produce amyloid-β peptides. Amyloid-β peptides are considered causative agents for Alzheimer's disease and drugs targeted at γ-secretase are investigated as therapeutic treatments. We synthesized new carprofen derivatives, which showed γ-secretase modulating activity and determined their precise position, orientation, and dynamics in lipid membranes by combining neutron diffraction, solid-state NMR spectroscopy, and molecular dynamics simulations. Our data indicate that the carprofen derivatives are inserted into the membrane interface, where the exact position and orientation depends on the lipid phase. This knowledge will help to understand the docking of carprofen derivatives to γ-secretase and in the design of new potent drugs. The approach presented here promises to serve as a general guideline how drug/target interactions in membranes can be analyzed in a comprehensive manner.

Aryl-Alkyl-Lysines Interact with Anionic Lipid Components of Bacterial Cell Envelope Eliciting Anti-Inflammatory and Antibiofilm Properties

The emergence of bacterial resistance and hesitance in approving new drugs has bolstered research on membrane-active agents such as antimicrobial peptides and their synthetic derivatives as therapeutic alternatives against bacterial infections. Herein, we document the action of aryl-alkyl-lysines on liposomes mimicking bacterial membranes using solid-state nuclear magnetic resonance spectroscopy. A significant perturbation of the lipid thickness and order parameter of the lipid membrane was observed upon treatment with this class of compounds. Encouraged by these results, the ability of the most active compound (NCK-10) to interact with aggregates of lipopolysaccharides (LPSs) was studied. In vitro experiments showed that NCK-10 was able to prevent the LPS-induced stimulation of proinflammatory cytokines such as tumor necrosis factor-α and interleukin-6. The compound could also disrupt the biofilms of Pseudomonas aeruginosa in vitro and bring down the bacterial burden by more than 99% in a mice model of burn infections caused by the biofilms of P. aeruginosa.

Spectroscopic Studies as a Toolbox for Biophysical and Chemical Characterization of Lipid-Based Nanotherapeutics

The goal of this study is to provide tools to minimize trial-and-error in the development of novel lipid-based nanotherapeutics, in favor of a rational design process. For this purpose, we present case-study examples of biophysical assays that help addressing issues of lipid-based nanotherapeutics' profiling and assist in the design of lipid nanocarriers for therapeutic usage. The assays presented are rooted in spectroscopic methods (steady-state and time-resolved fluorescence; UV-Vis derivative spectroscopy; fluorescence anisotropy and fluorescence lifetime image microscopy) and allow accessing physical-chemical interactions between drugs and lipid nanocarriers, as well as studying interactions between lipid-based nanotherapeutics and membranes and/or proteins, as this is a key factor in predicting their therapeutic and off target effects. Derivative spectroscopy revealed Naproxen's high distribution (LogD ≈ 3) in different lipid-based nanocarriers (micelles and unilamellar or multilamellar vesicles) confirming the adequacy of such systems for encapsulating this anti-inflammatory drug. Fluorescence quenching studies revealed that the anti-inflammatory drugs Acemetacin and Indomethacin can reach an inner location at the lipid nanocarrier while being anchored with its carboxylic moiety at the polar headgroup. The least observed quenching effect suggested that Tolmetin is probably located at the polar headgroup region of the lipid nanocarriers and this superficial location may translate in a fast drug release from the nanocarriers. Fluorescent anisotropy measurements indicated that the drugs deeply buried within the lipid nanocarrier where the ones that had a greater fluidizing effect which can also translate in a faster drug release. The drug binding strength to serum albumin was also compared for a free drug (Clonixin) or for the same drug after encapsulation in a lipid nanocarrier DSPC:DODAP (2:1). Under both conditions there is a strong binding to serum albumin, at one binding site, suggesting the need to produce a stealth nanosystem. Finally the cellular uptake of lipid nanocarriers loaded with Daunorubicin was investigated in cancer cells using fluorescence lifetime imaging microscopy. From the images obtained it was possible to conclude that even at short incubation times (15 min) there was a distribution of the drug in the cytoplasm, whereas for longer incubation periods (4 h) the drug has reached the nucleus.

Effect of ergosterol on the fungal membrane properties. All-atom and coarse-grained molecular dynamics study

Cell membranes are complex multicomponent systems consisting of thousands of different lipids with numerous embedded membrane proteins and many types of sterols. We used all-atom and coarse-grained molecular dynamics simulations to study the structural and dynamical properties of phospholipid bilayers containing four types of phospholipids and different amount of ergosterol, main sterol component in the fungal membranes. To characterize the influence of ergosterol on the membrane properties we analyzed the surface area per lipid, bilayer thickness, area compressibility modulus, mass density profiles, deuterium order parameters, and lateral diffusion coefficients. The presence of ergosterol induces the ordering of lipids leading to their denser packing, to reducing the lateral diffusion of lipids and lipid surface area, to increasing the thickness of bilayer and compressibility modulus. In addition, we evaluated each calculated property between the two simulation methods.

On the Calculation of Acyl Chain Order Parameters from Lipid Simulations

For molecular dynamics simulations of biological membrane systems to live up to the potential of providing accurate atomic level detail into membrane properties and functions, it is essential that the force fields used to model such systems are as accurate as possible. One membrane property that is often used to assess force field accuracy is the carbon-hydrogen (or carbon-deuterium) order parameters of the lipid tails, which can be accurately measured using experimental NMR techniques. There are a variety of analysis tools available to calculate these order parameters from simulations and it is essential that these computational tools work correctly to ensure the accurate assessment of the simulation force fields. In this work we compare many of these computational tools for calculating the order parameters of POPC membranes. While tools that work on all-atom systems and tools that work on saturated lipid tails in general work extremely well, we demonstrate that the majority of the tested tools that calculate the order parameters for unsaturated united-atom lipid tails do so incorrectly. We identify tools that do perform accurate calculations and include one such program with this work, enabling rapid and accurate calculation of united-atom lipid order parameters. Furthermore, we discuss cases in which it is nontrivial to appropriately predict the unsaturated carbon order parameters in united-atom systems. Finally, we examine order parameter splitting for carbon 2 in sn-2 lipid chains, demonstrating substantial deviations from experimental values in several all-atom and united-atom lipid force fields.

Concepts and Methods of Solid-State NMR Spectroscopy Applied to Biomembranes

Concepts of solid-state NMR spectroscopy and applications to fluid membranes are reviewed in this paper. Membrane lipids with ²H-labeled acyl chains or polar head groups are studied using ²H NMR to yield knowledge of their atomistic structures in relation to equilibrium properties. This review demonstrates the principles and applications of solid-state NMR by unifying dipolar and quadrupolar interactions and highlights the unique features offered by solid-state ²H NMR with experimental illustrations. For randomly oriented multilamellar lipids or aligned membranes, solid-state ²H NMR enables direct measurement of residual quadrupolar couplings (RQCs) due to individual C-²H-labeled segments. The distribution of RQC values gives nearly complete profiles of the segmental order parameters S_CD⁽ⁱ⁾ as a function of acyl segment position (i). Alternatively, one can measure residual dipolar couplings (RDCs) for natural abundance lipid samples to obtain segmental S_CH order parameters. A theoretical mean-torque model provides acyl-packing profiles representing the cumulative chain extension along the normal to the aqueous interface. Equilibrium structural properties of fluid bilayers and various thermodynamic quantities can then be calculated, which describe the interactions with cholesterol, detergents, peptides, and integral membrane proteins and formation of lipid rafts. One can also obtain direct information for membrane-bound peptides or proteins by measuring RDCs using magic-angle spinning (MAS) in combination with dipolar recoupling methods. Solid-state NMR methods have been extensively applied to characterize model membranes and membrane-bound peptides and proteins, giving unique information on their conformations, orientations, and interactions in the natural liquid-crystalline state.

pH-Dependent Membrane Interactions of the Histidine-Rich Cell-Penetrating Peptide LAH4-L1

The histidine-rich designer peptide LAH4-L1 exhibits antimicrobial and potent cell-penetrating activities for a wide variety of cargo including nucleic acids, polypeptides, adeno-associated viruses, and nanodots. The non-covalent complexes formed between the peptide and cargo enter the cell via an endosomal pathway where the pH changes from neutral to acidic. Here, we investigated the membrane interactions of the peptide with phospholipid bilayers and its membrane topology using static solid-state NMR spectroscopy. Oriented ¹⁵N solid-state NMR indicates that in membranes composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine (POPS) 3:1 mol/mole and at neutral pH, the peptide adopts transmembrane topologies. Furthermore, ³¹P and ²H solid-state NMR spectra show that liquid crystalline 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and POPC/POPS 3:1 liposomes retain a bilayer macroscopic phase even at the highest peptide concentrations investigated, with an oblate orientational distribution of the phospholipids at a peptide/lipid ratio of 1:5. At pH 5, as it occurs in the endosome, the alignment of LAH4-L1 at a peptide/lipid ratio of 1:25 is predominantly parallel to POPC/POPS 3:1 bilayers (prolate deformation) when at the same time it induces a considerable decrease of the deuterium order parameter of POPC/²H₃₁-POPS 3:1. In addition, when studied in mechanically supported lipid membranes, a pronounced disordering of the phospholipid alignment is observed. In the presence of even higher peptide concentrations, lipid spectra are observed that suggest the formation of magnetically oriented or isotropic bicelles. This membrane-disruptive effect is enhanced for gel phase DMPC membranes. By protonation of the four histidines in acidic environments, the overall charge and hydrophobic moment of LAH4-L1 considerably change, and much of the peptide is released from the cargo. Thus, the amphipathic peptide sequences become available to disrupt the endosomal membrane and to assure highly efficient release from this organelle.

The in vivo structure of biological membranes and evidence for lipid domains

Examining the fundamental structure and processes of living cells at the nanoscale poses a unique analytical challenge, as cells are dynamic, chemically diverse, and fragile. A case in point is the cell membrane, which is too small to be seen directly with optical microscopy and provides little observational contrast for other methods. As a consequence, nanoscale characterization of the membrane has been performed ex vivo or in the presence of exogenous labels used to enhance contrast and impart specificity. Here, we introduce an isotopic labeling strategy in the gram-positive bacterium Bacillus subtilis to investigate the nanoscale structure and organization of its plasma membrane in vivo. Through genetic and chemical manipulation of the organism, we labeled the cell and its membrane independently with specific amounts of hydrogen (H) and deuterium (D). These isotopes have different neutron scattering properties without altering the chemical composition of the cells. From neutron scattering spectra, we confirmed that the B. subtilis cell membrane is lamellar and determined that its average hydrophobic thickness is 24.3 ± 0.9 Ångstroms (Å). Furthermore, by creating neutron contrast within the plane of the membrane using a mixture of H- and D-fatty acids, we detected lateral features smaller than 40 nm that are consistent with the notion of lipid rafts. These experiments-performed under biologically relevant conditions-answer long-standing questions in membrane biology and illustrate a fundamentally new approach for systematic in vivo investigations of cell membrane structure.

Biophysical study of resin acid effects on phospholipid membrane structure and properties

Hydrophobic resin acids (RAs) are synthesized by conifer trees as part of their defense mechanisms. One of the functions of RAs in plant defense is suggested to be the perturbation of the cellular membrane. However, there is a vast diversity of chemical structures within this class of molecules, and there are no clear correlations to the molecular mechanisms behind the RA's toxicity. In this study we unravel the molecular interactions of the three closely related RAs dehydroabietic acid, neoabietic acid, and the synthetic analogue dichlorodehydroabietic acid with dipalmitoylphosphatidylcholine (DPPC) model membranes and the polar lipid extract of soybeans. The complementarity of the biophysical techniques used (NMR, DLS, NR, DSC, Cryo-TEM) allowed correlating changes at the vesicle level with changes at the molecular level and the co-localization of RAs within DPPC monolayer. Effects on DPPC membranes are correlated with the physical chemical properties of the RA and their toxicity.

Supercomputing Studies of Biomembranes

Although computer simulation of biological molecules has seen widespread growth and is widely accepted as an important biochemical tool, it is hampered by lim ited computing resources. Biomolecular systems, by necessity, contain a large number of interaction sites. In many cases, these sites interact over quite large dis tances. Further, the time scales of biological interest are long, which requires that simulations of dynamical properties at the atomic level must be lengthy to ade quately probe these motions. We address these issues through discussions of atomic-level molecular dynam ics simulations of biological lipid bilayer membranes, which are key constructs in biochemistry. These simu lations reproduce many experimental observables and provide a degree of resolution currently unavailable experimentally. The lengths of these simulations, the longest of which was 2 nanoseconds, were sufficient to effectively sample many of the motions governing the behavior of biomembranes. Examples are given showing the importance of long-range interactions. The number of interaction sites required by these sim ulations is discussed, particularly the need for explicit representation of solvent molecules.

Tyrosine Residues from the S4-S5 Linker of Kv11.1 Channels Are Critical for Slow Deactivation

Slow deactivation of Kv11.1 channels is critical for its function in the heart. The S4-S5 linker, which joins the voltage sensor and pore domains, plays a critical role in this slow deactivation gating. Here, we use NMR spectroscopy to identify the membrane-bound surface of the S4S5 linker, and we show that two highly conserved tyrosine residues within the KCNH subfamily of channels are membrane-associated. Site-directed mutagenesis and electrophysiological analysis indicates that Tyr-542 interacts with both the pore domain and voltage sensor residues to stabilize activated conformations of the channel, whereas Tyr-545 contributes to the slow kinetics of deactivation by primarily stabilizing the transition state between the activated and closed states. Thus, the two tyrosine residues in the Kv11.1 S4S5 linker play critical but distinct roles in the slow deactivation phenotype, which is a hallmark of Kv11.1 channels.

Effect of Sterol Structure on the Physical Properties of 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine Membranes Determined Using (2)H Nuclear Magnetic Resonance

The effect of a series of phytosterols on lipid chain ordering in 1-palmitoyl((2)H31)-2-oleoyl-sn-glycero-3-phosphocholine (POPC-d31) multibilayer vesicles was examined by (2)H NMR spectroscopy at 25 °C. These results, along with existing data for other sterols, indicate that the ordering power of sterols in POPC-d31 depends on subtle aspects of sterol structure. Cholesterol, 7-dehydrocholesterol (7-DHC), campesterol, β-sitosterol, ergosterol, brassicasterol, and stigmasterol all increase the lipid chain order as sterol concentration is increased. However, saturation of the ordering occurs at different sterol concentrations for ergosterol (as previously reported), brassicasterol, β-sitosterol, and stigmasterol. Here our interest lies in finding which part of the sterol structure is responsible for the observed saturation of the palmitoyl chain order as a function of sterol concentration. In particular, we propose that the saturation of the ordering of POPC-d31/brassicasterol and POPC-d31/stigmasterol membranes at quite low sterol concentrations is due to the presence of a double bond at C22. We also discuss how the structural differences between the sterols affect their ability to intercalate between the POPC acyl chains. Furthermore, the effective solubility of sterols in POPC is discussed in relation to the dependence of maximum POPC-d31 chain order vs sterol concentration.

Cholesterol-induced suppression of membrane elastic fluctuations at the atomistic level

Applications of solid-state NMR spectroscopy for investigating the influences of lipid-cholesterol interactions on membrane fluctuations are reviewed in this paper. Emphasis is placed on understanding the energy landscapes and fluctuations at an emergent atomistic level. Solid-state (2)H NMR spectroscopy directly measures residual quadrupolar couplings (RQCs) due to individual C-(2)H labeled segments of the lipid molecules. Moreover, residual dipolar couplings (RDCs) of (13)C-(1)H bonds are obtained in separated local-field NMR spectroscopy. The distributions of RQC or RDC values give nearly complete profiles of the order parameters as a function of acyl segment position. Measured equilibrium properties of glycerophospholipids and sphingolipids including their binary and tertiary mixtures with cholesterol show unequal mixing associated with liquid-ordered domains. The entropic loss upon addition of cholesterol to sphingolipids is less than for glycerophospholipids and may drive the formation of lipid rafts. In addition relaxation time measurements enable one to study the molecular dynamics over a wide time-scale range. For (2)H NMR the experimental spin-lattice (R1Z) relaxation rates follow a theoretical square-law dependence on segmental order parameters (SCD) due to collective slow dynamics over mesoscopic length scales. The functional dependence for the liquid-crystalline lipid membranes is indicative of viscoelastic properties as they emerge from atomistic-level interactions. A striking decrease in square-law slope upon addition of cholesterol denotes stiffening relative to the pure lipid bilayers that is diminished in the case of lanosterol. Measured equilibrium properties and relaxation rates infer opposite influences of cholesterol and detergents on collective dynamics and elasticity at an atomistic scale that potentially affects lipid raft formation in cellular membranes.

Computational studies of plasma lipoprotein lipids

Plasma lipoproteins are macromolecular assemblies of proteins and lipids found in the blood. The lipid components of lipoproteins are amphipathic lipids such as phospholipids (PLs), and unesterified cholesterols (UCs) and hydrophobic lipids such as cholesteryl esters (CEs) and triglycerides (TGs). Since lipoproteins are soft matter supramolecular assemblies easily deformable by thermal fluctuations and they also exist in varying densities and protein/lipid components, a detailed understanding of their structure/function is experimentally difficult. Molecular dynamics (MD) simulation has emerged as a particularly promising way to explore the structure and dynamics of lipoproteins. The purpose of this review is to survey the current status of computational studies of the lipid components of the lipoproteins. Computational studies aim to explore three levels of complexity for the 3-dimensional structural dynamics of lipoproteins at various metabolic stages: (i) lipoprotein particles consist of protein with minimal lipid; (ii) lipoprotein particles consist of PL-rich discoidal bilayer-like lipid particles; (iii) mature circulating lipoprotein particles consist of CE-rich or TG-rich spheroidal lipid-droplet-like particles. Due to energy barriers involved in conversion between these species, other biomolecules also participate in lipoprotein biological assembly. For example: (i) lipid-poor apolipoprotein A-I (apoA-I) interacts with ATP-binding cassette transporter A1 (ABCA1) to produce nascent discoidal high density lipoprotein (dHDL) particles; (ii) lecithin-cholesterol acyltransferase (LCAT) mediates the conversion of UC to CE in dHDL, driving spheroidal HDL (sHDL) formation; (iii) transfer proteins, cholesterol ester transfer protein (CETP) and phospholipid transfer protein (PLTP), transfer both CE and TG and PL, respectively, between lipoprotein particles. Computational studies have the potential to explore different lipoprotein particles at each metabolic stage in atomistic detail. This review discusses the current status of computational methods including all-atom MD (AAMD), coarse-grain MD (CGMD), and MD-simulated annealing (MDSA) and their applications in lipoprotein structural dynamics and biological assemblies. Results from MD simulations are discussed and compared across studies in order to identify key findings, controversies, issues and future directions. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

C60 fullerene localization and membrane interactions in RAW 264.7 immortalized mouse macrophages

There continues to be a significant increase in the number and complexity of hydrophobic nanomaterials that are engineered for a variety of commercial purposes making human exposure a significant health concern. This study uses a combination of biophysical, biochemical and computational methods to probe potential mechanisms for uptake of C60 nanoparticles into various compartments of living immune cells. Cultures of RAW 264.7 immortalized murine macrophage were used as a canonical model of immune-competent cells that are likely to provide the first line of defense following inhalation. Modes of entry studied were endocytosis/pinocytosis and passive permeation of cellular membranes. The evidence suggests marginal uptake of C60 clusters is achieved through endocytosis/pinocytosis, and that passive diffusion into membranes provides a significant source of biologically-available nanomaterial. Computational modeling of both a single molecule and a small cluster of fullerenes predicts that low concentrations of fullerenes enter the membrane individually and produce limited perturbation; however, at higher concentrations the clusters in the membrane causes deformation of the membrane. These findings are bolstered by nuclear magnetic resonance (NMR) of model membranes that reveal deformation of the cell membrane upon exposure to high concentrations of fullerenes. The atomistic and NMR models fail to explain escape of the particle out of biological membranes, but are limited to idealized systems that do not completely recapitulate the complexity of cell membranes. The surprising contribution of passive modes of cellular entry provides new avenues for toxicological research that go beyond the pharmacological inhibition of bulk transport systems such as pinocytosis.

A computer model of a polyunsaturated monogalactolipid bilayer

1,2-di-O-acyl-3-O-β-D-galactopyranosyl-sn-glycerol (MGDG) is the main lipid component of thylakoid membranes of higher plants and algae. This monogalactolipid is thought of as a non-bilayer lipid but actually it can form both lamellar and nonlamellar phases. In this study, molecular dynamics (MD) simulations of the fully hydrated di-18:3 MGDG bilayer in the lamellar phase were carried out at 310 and 295 K for 200 and 450 ns, respectively, using the GROMACS 4 software package and OPLS-AA force field. At both temperatures, the lamellar phase of the systems was stable. The pure di-18:3 MGDG bilayer is the first step towards creating a computer model of the lipid matrix of the thylakoid membrane and the main aim of this study was to validate the computer model of di-18:3 MGDG in the bilayer and also to assess the properties of the bilayer. However, only a few of the predicted properties could be compared with those derived experimentally and in other MD simulations because of insufficient amount of such data. Thus, direct validation of the MGDG bilayer proved difficult. Therefore, in the validation process also an indirect approach was used, in which a computer model of the 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) bilayer simulated at the same temperatures using the same force field as the MGDG bilayer was assessed. Successful validation of the DOPC bilayer parameterized in the OPLS-AA force field and similar properties of the MGDG molecules in the pure 18:3 MGDG and in binary 18:3 MGDG-PC bilayers indicate that the computer model of the MDGD molecule is faithful and the MGDG bilayer is representative on the time scales covered in these MD simulations.

Oriented samples: a tool for determining the membrane topology and the mechanism of action of cationic antimicrobial peptides by solid-state NMR

Overuse and misuse of antibiotics have led bacteria to acquire several mechanisms of resistance. In response to this, researchers have identified natural antimicrobial peptides as promising candidates to fight against multiresistant bacteria. However, their mode of action is still unclear, and a better understanding of the mode of action of these peptides is of primary importance to develop new peptides displaying high antibacterial activity and low hemolytic activity. One of the main features that defines the mechanism of action is the membrane topology of the peptide. Among the spectroscopic techniques, solid-state NMR is the technique of choice for determining the location of the peptide within the membrane. It can be achieved by performing experiments with oriented samples. In the literature, the two most common types of oriented samples are bicelles and phospholipids mechanically oriented between glass plates. The mode of perturbation of the membrane-active peptide can be studied by phosphorus-31 and deuterium NMR. On the other hand, several experiments such as nitrogen-15 and fluorine solid-state NMR, that require labeled peptides, can give valuable information on the membrane topology of the peptide. The combination of the latter techniques allows the determination of a precise topology, thus a better knowledge of the molecular determinants involved in the membrane interactions of antimicrobial peptides.

Area per lipid and cholesterol interactions in membranes from separated local-field (13)C NMR spectroscopy

Investigations of lipid membranes using NMR spectroscopy generally require isotopic labeling, often precluding structural studies of complex lipid systems. Solid-state (13)C magic-angle spinning NMR spectroscopy at natural isotopic abundance gives site-specific structural information that can aid in the characterization of complex biomembranes. Using the separated local-field experiment DROSS, we resolved (13)C-(1)H residual dipolar couplings that were interpreted with a statistical mean-torque model. Liquid-disordered and liquid-ordered phases were characterized according to membrane thickness and average cross-sectional area per lipid. Knowledge of such structural parameters is vital for molecular dynamics simulations, and provides information about the balance of forces in membrane lipid bilayers. Experiments were conducted with both phosphatidylcholine (dimyristoylphosphatidylcholine (DMPC) and palmitoyloleoylphosphatidylcholine (POPC)) and egg-yolk sphingomyelin (EYSM) lipids, and allowed us to extract segmental order parameters from the (13)C-(1)H residual dipolar couplings. Order parameters were used to calculate membrane structural quantities, including the area per lipid and bilayer thickness. Relative to POPC, EYSM is more ordered in the ld phase and experiences less structural perturbation upon adding 50% cholesterol to form the lo phase. The loss of configurational entropy is smaller for EYSM than for POPC, thus favoring its interaction with cholesterol in raftlike lipid systems. Our studies show that solid-state (13)C NMR spectroscopy is applicable to investigations of complex lipids and makes it possible to obtain structural parameters for biomembrane systems where isotope labeling may be prohibitive.

Mimicking and Understanding the Agglutination Effect of the Antimicrobial Peptide Thanatin Using Model Phospholipid Vesicles

Thanatin is a cationic 21-residue antimicrobial and antifongical peptide found in the spined soldier bug Podisus maculiventris. It is believed that it does not permeabilize membranes but rather induces the agglutination of bacteria and inhibits cellular respiration. To clarify its mode of action, lipid vesicle organization and aggregation propensity as well as peptide secondary structure have been studied using different membrane models. Dynamic light scattering and turbidimetry results show that specific mixtures of negatively charged and zwitterionic phospholipid vesicles are able to mimic the agglutination effect of thanatin observed on Gram-negative and Gram-positive bacterial cells, while monoconstituent ("conventional") models cannot reproduce this phenomenon. The model of eukaryotic cell reveals no particular interaction with thanatin, which is consistent with the literature. Infrared spectroscopy shows that under the conditions under which vesicle agglutination occurs, thanatin exhibits a particular spectral pattern in the amide I' region and in the region associated with Arg side chains. The data suggest that thanatin mainly retains its hairpin structure, Arg residues being involved in strong interactions with anionic groups of phospholipids. In the absence of vesicle agglutination, the peptide conformation and Arg side-chain environment are similar to those observed in solution. The data show that a negatively charged membrane is required for thanatin to be active, but this condition is insufficient. The activity of thanatin seems to be modulated by the charge surface density of membranes and thanatin concentration.

Interaction between the NS4B amphipathic helix, AH2, and charged lipid headgroups alters membrane morphology and AH2 oligomeric state--Implications for the Hepatitis C virus life cycle

The non-structural protein 4B (NS4B) from Hepatitis C virus (HCV) plays a pivotal role in the remodelling of the host cell's membranes, required for the formation of the viral replication complex where genome synthesis occurs. NS4B is an integral membrane protein that possesses a number of domains vital for viral replication. Structural and biophysical studies have revealed that one of these, the second amphipathic N-terminal helix (AH2), plays a key role in these remodelling events. However, there is still limited understanding of the mechanism through which AH2 promotes these changes. Here we report on solid-state NMR and molecular dynamics studies that demonstrate that AH2 promotes the clustering of negatively charged lipids within the bilayer, a process that reduces the strain within the bilayer facilitating the remodelling of the lipid bilayer. Furthermore, the presence of negatively charged lipids within the bilayer appears to promote the disassociation of AH2 oligomers, highlighting a potential role for lipid recruitment in regulating NS protein interactions.

•

Changes in membrane morphology studied by ²H and ³¹P Solid-state NMR.

•

Bilayer charge influences the oligomeric state of the amphipathic helix AH2 from NS4B.

•

Interaction of AH2 with charged lipid membranes reduces strain within bilayer.

•

AH2 from NS4B is involved in membrane remodelling and membranous web formation.

•

Lipid bilayer/NS4B interactions may regulate Hepatitis C virus lifecycle.

Membrane interactions in small fast-tumbling bicelles as studied by 31P NMR

Small fast-tumbling bicelles are ideal for studies of membrane interactions at molecular level; they allow analysis of lipid properties using solution-state NMR. In the present study we used 31P NMR relaxation to obtain detailed information on lipid head-group dynamics. We explored the effect of two topologically different membrane-interacting peptides on bicelles containing either dimyristoylphosphocholine (DMPC), or a mixture of DMPC and dimyristoylphosphoglycerol (DMPG), and dihexanoylphosphocholine (DHPC). KALP21 is a model transmembrane peptide, designed to span a DMPC bilayer and dynorphin B is a membrane surface active neuropeptide. KALP21 causes significant increase in bicelle size, as evidenced by both dynamic light scattering and 31P T2 relaxation measurements. The effect of dynorphin B on bicelle size is more modest, although significant effects on T2 relaxation are observed at higher temperatures. A comparison of 31P T1 values for the lipids with and without the peptides showed that dynorphin B has a greater effect on lipid head-group dynamics than KALP21, especially at elevated temperatures. From the field-dependence of T1 relaxation data, a correlation time describing the overall lipid motion was derived. Results indicate that the positively charged dynorphin B decreases the mobility of the lipid molecules--in particular for the negatively charged DMPG--while KALP21 has a more modest influence. Our results demonstrate that while a transmembrane peptide has severe effects on overall bilayer properties, the surface bound peptide has a more dramatic effect in reducing lipid head-group mobility. These observations may be of general importance for understanding peptide-membrane interactions.

Molecular simulation of ibuprofen passing across POPC membrane

Permeability assessment is an important procedure in the drug development process, and drug partitioning in membrane bilayer is related to permeability. To investigate the pH dependence on drug partitioning, the process of different ionization state of ibuprofen passing across POPC bilayer was studied using molecular dynamics simulation. The results show that both atomic charge scheme and ionization state of the drug affect the value and shape of energy profile when passing across the POPC bilayer. The neutral ibuprofen (ibuprofen under acidic condition) has a much lower energy barrier as compared with the anionic ibuprofen (ibuprofen under basic condition). Meantime, hydrogen bond analysis also certifies that it is easy for neutral ibuprofen to pass from bulk water to bilayer center. Our calculation suggests that the ionization state of ibuprofen may be changed between neutral and anionic state when passing across membrane: it may be ionized outside the membrane and neutralized inside the membrane.

Interaction Study of an Amorphous Solid Dispersion of Cyclosporin A in Poly-Alpha-Cyclodextrin with Model Membranes by (1)H-, (2)H-, (31)P-NMR and Electron Spin Resonance

The properties of an amorphous solid dispersion of cyclosporine A (ASD) prepared with the copolymer alpha cyclodextrin (POLYA) and cyclosporine A (CYSP) were investigated by ¹H-NMR in solution and its membrane interactions were studied by ¹H-NMR in small unilamellar vesicles and by ³¹P ²H NMR in phospholipidic dispersions of DMPC (dimyristoylphosphatidylcholine) in comparison with those of POLYA and CYSP alone. ¹H-NMR chemical shift variations showed that CYSP really interacts with POLYA, with possible adduct formation, dispersion in the solid matrix of the POLYA, and also complex formation. A coarse approach to the latter mechanism was tested using the continuous variations method, indicating an apparent 1 : 1 stoichiometry. Calculations gave an apparent association constant of log Ka = 4.5. A study of the interactions with phospholipidic dispersions of DMPC showed that only limited interactions occurred at the polar head group level (³¹P). Conversely, by comparison with the expected chain rigidification induced by CYSP, POLYA induced an increase in the fluidity of the layer while ASD formation led to these effects almost being overcome at 298 K. At higher temperature, while the effect of CYSP seems to vanish, a resulting global increase in chain fluidity was found in the presence of ASD.