Saturation with cholesterol increases vertical order and smoothes the surface of the phosphatidylcholine bilayer: a molecular simulation study

Mode of molecular interaction of triterpenoid saponin ginsenoside Rh2 with membrane lipids in liquid-disordered phases

Ginsenoside Rh2 (Rh2) is a ginseng saponin comprising a triterpene core and one unit of glucose and has attracted much attention due to its diverse biological activities. In the present study, we used small-angle X-ray diffraction, solid-state NMR, fluorescence microscopy, and MD simulations to investigate the molecular interaction of Rh2 with membrane lipids in the liquid-disordered (Ld) phase mainly composed of palmitoyloleoylphosphatidylcholine compared with those in liquid-ordered (Lo) phase mainly composed of sphingomyelin and cholesterol. The electron density profiles determined by X-ray diffraction patterns indicated that Rh2 tends to be present in the shallow interior of the bilayer in the Ld phase, while Rh2 accumulation was significantly smaller in the Lo phase. Order parameters at intermediate depths in the bilayer leaflet obtained from 2H NMR spectra and MD simulations indicated that Rh2 reduces the order of the acyl chains of lipids in the Ld phase. The dihydroxy group and glucose moiety at both ends of the hydrophobic triterpene core of Rh2 cause tilting of the molecular axis relative to the membrane normal, which may enhance membrane permeability by loosening the packing of lipid acyl chains. These features of Rh2 are distinct from steroidal saponins such as digitonin and dioscin, which exert strong membrane-disrupting activity.

Interaction of βL- and γ-Crystallin with Phospholipid Membrane Using Atomic Force Microscopy

Highly concentrated lens proteins, mostly β- and γ-crystallin, are responsible for maintaining the structure and refractivity of the eye lens. However, with aging and cataract formation, β- and γ-crystallin are associated with the lens membrane or other lens proteins forming high-molecular-weight proteins, which further associate with the lens membrane, leading to light scattering and cataract development. The mechanism by which β- and γ-crystallin are associated with the lens membrane is unknown. This work aims to study the interaction of β- and γ-crystallin with the phospholipid membrane with and without cholesterol (Chol) with the overall goal of understanding the role of phospholipid and Chol in β- and γ-crystallin association with the membrane. Small unilamellar vesicles made of Chol/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (Chol/POPC) membranes with varying Chol content were prepared using the rapid solvent exchange method followed by probe tip sonication and then dispensed on freshly cleaved mica disk to prepare a supported lipid membrane. The βL- and γ-crystallin from the cortex of the bovine lens was used to investigate the time-dependent association of βL- and γ-crystallin with the membrane by obtaining the topographical images using atomic force microscopy. Our study showed that βL-crystallin formed semi-transmembrane defects, whereas γ-crystallin formed transmembrane defects on the phospholipid membrane. The size of semi-transmembrane defects increases significantly with incubation time when βL-crystallin interacts with the membrane. In contrast, no significant increase in transmembrane defect size was observed in the case of γ-crystallin. Our result shows that Chol inhibits the formation of membrane defects when βL- and γ-crystallin interact with the Chol/POPC membrane, where the degree of inhibition depends upon the amount of Chol content in the membrane. At a Chol/POPC mixing ratio of 0.3, membrane defects were observed when both βL- and γ-crystallin interacted with the membrane. However, at a Chol/POPC mixing ratio of 1, no association of γ-crystallin with the membrane was observed, which resulted in a defect-free membrane, and the severity of the membrane defect was decreased when βL-crystallin interacted with the membrane. The semi-transmembrane or transmembrane defects formed by the interaction of βL- and γ-crystallin on phospholipid membrane might be responsible for light scattering and cataract formation. However, Chol suppressed the formation of such defects in the membrane, likely maintaining lens membrane homeostasis and protecting against cataract formation.

Cholesterol catalyzes unfolding in membrane-inserted motifs of the pore forming protein cytolysin A

Plasma membrane-induced protein folding and conformational transitions play a central role in cellular homeostasis. Several transmembrane proteins are folded in the complex lipid milieu to acquire a specific structure and function. Bacterial pore forming toxins (PFTs) are proteins expressed by a large class of pathogenic bacteria that exploit the plasma membrane environment to efficiently undergo secondary structure changes, oligomerize, and form transmembrane pores. Unregulated pore formation causes ion imbalance, leading to cell death and infection. Determining the free energy landscape of these membrane-driven-driven transitions remains a challenging problem. Although cholesterol recognition is required for lytic activity of several proteins in the PFT family of toxins, the regulatory role of cholesterol for the α-PFT, cytolysin A expressed by Escherichia coli remains unexplained. In a recent free energy computation, we showed that the β tongue, a critical membrane-inserted motif of the ClyA toxin, has an on-pathway partially unfolded intermediate that refolds into the helix-turn-helix motif of the pore state. To understand the molecular role played by cholesterol, we carry out string-method-based computations in membranes devoid of cholesterol, which reveals an increase of ∼30 times in the free energy barrier for the loss of β sheet secondary structure when compared with membranes containing cholesterol. Specifically, the tyrosine-cholesterol interaction was found to be critical to creating the unfolded intermediate. Cholesterol also increases the packing and hydrophobicity of the bilayer, resulting in enhanced interactions of the bound protein before complete membrane insertion. Our study illustrates that cholesterol is critical to catalyzing and stabilizing the membrane-inserted unfolded state of the β tongue motif of ClyA, opening up fresh insights into cholesterol-assisted unfolding of membrane proteins.

Coarse-grained modeling of polystyrene-modified CNTs and their interactions with lipid bilayers

In the present work, we describe Martini3 coarse-grained models of polystyrene and carboxyl-terminated polystyrene functionalized carbon nanotubes (CNTs) and investigate their interactions with lipid bilayers with and without cholesterol (CHOL) using molecular dynamics simulations. By changing the polystyrene chain length and grafting density at the end ring of the CNTs at two different nanotube concentrations, we observe the translocation of nanoparticles as well as changes in the lipid bilayer properties. Our results show that all developed models passively diffuse into the membranes without causing any damage to the membrane integrity, although high concentrations of CNTs induce structural and elastic changes in lipid bilayers. In the presence of CHOL, increasing CNT concentration results in decreased rates of CHOL transmembrane motions. On the other hand, CNTs are prone to lipid and polystyrene blockage, which affects their equilibrated configurations, and tilting behavior within the membranes. Hence, we demonstrate that polystyrene-functionalized CNTs are promising drug-carrier agents. However, polystyrene chain length and grafting density are important factors to consider to enhance the efficiency of drug delivery.

Membrane Models and Experiments Suitable for Studies of the Cholesterol Bilayer Domains

Cholesterol (Chol) is an essential component of animal cell membranes and is most abundant in plasma membranes (PMs) where its concentration typically ranges from 10 to 30 mol%. However, in red blood cells and Schwann cells, PMs Chol content is as high as 50 mol%, and in the PMs of the eye lens fiber cells, it can reach up to 66 mol%. Being amphiphilic, Chol molecules are easily incorporated into the lipid bilayer where they affect the membrane lateral organization and transmembrane physical properties. In the aqueous phase, Chol cannot form free bilayers by itself. However, pure Chol bilayer domains (CBDs) can form in lipid bilayer membranes with the Chol content exceeding 50 mol%. The range of Chol concentrations surpassing 50 mol% is less frequent in biological membranes and is consequently less investigated. Nevertheless, it is significant for the normal functioning of the eye lens and understanding how Chol plaques form in atherosclerosis. The most commonly used membrane models are unilamellar and multilamellar vesicles (MLVs) and supported lipid bilayers (SLBs). CBDs have been observed directly using confocal microscopy, X-ray reflectometry and saturation recovery electron paramagnetic resonance (SR EPR). Indirect evidence of CBDs has also been reported by using atomic force microscopy (AFM) and fluorescence recovery after photobleaching (FRAP) experiments. The overall goal of this review is to demonstrate the advantages and limitations of the various membrane models and experimental techniques suitable for the detection and investigation of the lateral organization, function and physical properties of CBDs.

Membrane binding of pore-forming γ-hemolysin components studied at different lipid compositions

Methicillin-resistant Staphylococcus aureus is among those pathogens currently posing the highest threat to public health. Its host immune evasion strategy is mediated by pore-forming toxins (PFTs), among which the bi-component γ-hemolysin is one of the most common. The complexity of the porogenesis mechanism by γ-hemolysin poses difficulties in the development of antivirulence therapies targeting PFTs from S. aureus, and sparse and apparently contrasting experimental data have been produced. Here, through a large set of molecular dynamics simulations at different levels of resolution, we investigate the first step of pore formation, and in particular the effect of membrane composition on the ability of γ-hemolysin components, LukF and Hlg2, to steadily adhere to the lipid bilayer in the absence of proteinaceous receptors. Our simulations are in agreement with experimental data of γ-hemolysin pore formation on model membranes, which are here explained on the basis of the bilayer properties. Our computational investigation suggests a possible rationale to explain experimental data on phospholipid binding to the LukF component, and to hypothesise a mechanism by which, on purely lipidic bilayers, the stable anchoring of LukF to the cell surface facilitates Hlg2 binding, through the exposure of its N-terminal region. We expect that further insights on the mechanism of transition between soluble and membrane bound-forms and on the role played by the lipid molecules will contribute to the design of antivirulence agents with enhanced efficacy against methicillin-resistant S. aureus infections.

Molecular oxygen as a probe molecule in EPR spin-labeling studies of membrane structure and dynamics

Molecular oxygen (O2) is the perfect probe molecule for membrane studies carried out using the saturation recovery EPR technique. O2 is a small, paramagnetic, hydrophobic enough molecule that easily partitions into a membrane's different phases and domains. In membrane studies, the saturation recovery EPR method requires two paramagnetic probes: a lipid-analog nitroxide spin label and an oxygen molecule. The experimentally derived parameters of this method are the spin-lattice relaxation times (T 1s) of spin labels and rates of bimolecular collisions between O2 and the nitroxide fragment. Thanks to the long T 1 of lipid spin labels (from 1 to 10 μs), the approach is very sensitive to changes of the local (around the nitroxide fragment) O2 diffusion-concentration product. Small variations in the lipid packing affect O2 solubility and O2 diffusion, which can be detected by the shortening of T 1 of spin labels. Using O2 as a probe molecule and a different lipid spin label inserted into specific phases of the membrane and membrane domains allows data about the lateral arrangement of lipid membranes to be obtained. Moreover, using a lipid spin label with the nitroxide fragment attached to its head group or a hydrocarbon chain at different positions also enables data about molecular dynamics and structure at different membrane depths to be obtained. Thus, the method can be used to investigate not only the lateral organization of the membrane (i.e., the presence of membrane domains and phases), but also the depth-dependent membrane structure and dynamics, and, hence, the membrane properties in three dimensions.

Multilamellar Liposomes as a Model for Biological Membranes: Saturation Recovery EPR Spin-Labeling Studies

EPR spin labeling has been used extensively to study lipids in model membranes to understand their structures and dynamics in biological membranes. The lipid multilamellar liposomes, which are the most commonly used biological membrane model, were prepared using film deposition methods and investigated with the continuous wave EPR technique (T2-sensitive spin-labeling methods). These investigations provided knowledge about the orientation of lipids, their rotational and lateral diffusion, and their rate of flip-flop between bilayer leaflets, as well as profiles of membrane hydrophobicity, and are reviewed in many papers and book chapters. In the early 1980s, the saturation recovery EPR technique was introduced to membrane studies. Numerous T1-sensitive spin-label methods were developed to obtain detailed information about the three-dimensional dynamic membrane structure. T1-sensitive methods are advantageous over T2-sensitive methods because the T1 of spin labels (1-10 μs) is 10 to 1000 times longer than the T2, which allows for studies of membrane dynamics in a longer time-space scale. These investigations used multilamellar liposomes also prepared using the rapid solvent exchange method. Here, we review works in which saturation recovery EPR spin-labeling methods were applied to investigate the properties of multilamellar lipid liposomes, and we discuss their relationships to the properties of lipids in biological membranes.

Membrane elasticity modulated by cholesterol in model of porcine eye lens-lipid membrane

Experimental evidence shows that the eye lens loses its elasticity dramatically with age. It has also been reported that the cholesterol (Chol) content in the eye lens fiber cell plasma membrane increases significantly with age. High Chol content leads to the formation of cholesterol bilayer domains (CBDs) in the lens membrane. The role of high Chol associated with lens elasticity is unclear. The purpose of this research is to investigate the membrane elasticity of the model of porcine lens-lipid (MPLL) membrane with increasing Chol content to elucidate the role of high Chol in lens membrane elasticity. In this study, we used atomic force microscopy (AFM) to study the mechanical properties (breakthrough force and area compressibility modulus (K_A)) of the MPLL membrane with increasing Chol content where K_A is the measure of membrane elasticity. We varied Chol concentration in Chol/MPLL membrane from 0 to ∼71 mol%. Supported Chol/MPLL membranes were prepared by fusion of small unilamellar vesicles (SUVs) on top of a flat mica surface. SUVs of the Chol/MPLL lipid mixture were prepared with the rapid solvent exchange method followed by probe-tip sonication. For the Chol/MPLL mixing ratio of 0, AFM image showed the formation of two distinct phases of the membrane, i.e., liquid-disordered phase (l_d) and solid-ordered phase (s_o) membrane. However, with Chol/MPLL mixing ratio of 0.5 and above, only liquid-ordered phase (l_o) membrane was formed. Also, two distinct breakthrough forces corresponding to l_d and s_o were observed for Chol/MPLL mixing ratio of 0, whereas only one breakthrough force was observed for membranes with Chol/MPLL mixing ratio of 0.5 and above. No significant difference in the membrane surface roughness was measured with increasing Chol content for these membranes; however, breakthrough force and K_A for l_o membrane increased when Chol/MPLL mixing ratio was increased from 0.5 to 1. Interestingly above the Chol/MPLL mixing ratio of 1, both breakthrough force and K_A decreased, indicating the formation of CBDs. Furthermore, these results showed that membrane elasticity increases at high Chol content, suggesting that high Chol content in lens membrane might be responsible for maintaining lens membrane elasticity.

Lutein and Zeaxanthin in the Lipid Bilayer-Similarities and Differences Revealed by Computational Studies

Lutein and zeaxanthin are two similar carotenoids of the xanthophyll subgroup. Carotenoids are synthesized almost entirely by plants but are also present in significant amounts in animals. They are essential components of the lipid matrix of biomembranes, and one of their functions is to protect cells from light radiation, free radicals and oxidative stress. Carotenoids, depending on their chemical structure, can locate at various positions and in different orientations in the bilayer. Xanthophylls (XAN) are polar and in the bilayer are positionally restricted. In the case of lutein and zeaxanthin, whose both ionone rings are hydroxy-substituted and as such are anchored in the lipid bilayer interfaces, the position is generally transmembrane. However, both experimental and computer modelling studies indicate that lutein can also locate horizontally below the bilayer interface. This location has never been observed for zeaxanthin. To find a molecular-level explanation for the difference in the orientations of the XAN molecules in the bilayer, a number of phosphatidylcholine-XAN bilayers were constructed and molecular dynamics (MD) simulated for 1.1 µs each. The all-trans XAN molecules were initially placed either parallel or perpendicular to the bilayer surface. With the exception of one lutein, the horizontally placed molecules adopted the transmembrane orientation within 100-600 ns. On the basis of detailed analyses of the XAN orientations and the numbers and lifetimes of their interactions in the bilayer, a plausible explanation is offered as to why a lutein molecule may remain in the horizontal orientation while zeaxanthin does not. Contrary to common believe, lutein horizontal orientation is not related to the ε-ring rotation around the C6'-C7' bond.

Modulation of Phospholipid Bilayer Properties by Simvastatin

Simvastatin (Zocor) is one of the most prescribed drugs for reducing high cholesterol. Although simvastatin is ingested in its inactive lactone form, it is converted to its active dihydroxyheptanoate form by carboxylesterases in the liver. The dihydroxyheptanoate form can also be converted back to its original lactone form. Unfortunately, some of the side effects associated with the intake of simvastatin and other lipophilic statins at higher doses include statin-associated myopathy (SAM) and, in more severe cases, kidney failure. While the cause of SAM is unknown, it is hypothesized that these side effects are dependent on the localization of statins in lipid bilayers and their impact on bilayer properties. In this work, we carry out all-atom molecular dynamics simulations on both the lactone and dihydroxyheptanoate forms of simvastatin (termed "SN" and "SA", respectively) with a pure 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer and a POPC/cholesterol (30 mol %) binary mixture as membrane models. Additional simulations were carried out with multiple simvastatin molecules to mimic in vitro conditions that produced pleiotropic effects. Both SN and SA spontaneously diffused into the lipid bilayer, and a longer simulation time of 4 μs was needed for the complete incorporation of multiple SAs into the bilayer. By constructing potential mean force and electron density profiles, we find that SN localizes deeper within the hydrophobic interior of the bilayer and that SA has a greater tendency to form hydrogen-bonding interactions with neighboring water molecules and lipid headgroups. For the pure POPC bilayer, both SN and SA increase membrane order, while membrane fluidity increases for the POPC/cholesterol bilayer.

On the Structure and Flip-flop of Free Docosahexaenoic Acid in a Model Human Brain Membrane

Among the omega-3 fatty acids, docosahexaenoic acid (DHA, sn22:6) is particularly vital in human brain cell membranes. There is considerable interest in DHA because low-level DHA has been associated with declined cognitive function and poor visual acuity. In this work, atomistic molecular dynamics simulations were used to investigate the effects of free protonated DHA (DHAP) in molar fractions of 0, 17, 30, and 38% in a realistic model of a healthy brain cell membrane comprising 26 lipid types. Numerous flip-flop events of DHAP were observed and categorized as successful or aborted. Novel use of the machine learning technique, density-based spatial clustering of applications with noise (DBSCAN), effectively identified flip-flop events by way of clustering. Our data show that increasing amounts of DHAP in the membrane disorder the bilayer packing, fluidize the membrane, and increase the rates of successful flip-flop from k = 0.2 μs^-1 (17% DHAP) to 0.8 μs^-1 (30% DHAP) and to 1.3 μs^-1 (38% DHAP). In addition, we also provided a detailed understanding of the flip-flop mechanism of DHAP across this complex membrane. Interestingly, we noted the role of hydrogen bonds in two distinct coordinated flip-flop phenomena between two DHAP molecules: double flip-flop and assisted flip-flop. Understanding the effects of various concentrations of DHAP on the dynamics within a lipid membrane and the resulting structural properties of the membrane are important when considering the use of DHAP as a dietary supplement or as a potential therapeutic.

Effects of cholesterol on chlorzoxazone translocation across POPC bilayer

Cholesterol plays a crucial role in modulating the physicochemical properties of membranes, thus influencing the membrane transport of drugs. In this paper, the effects caused by cholesterol on the membrane transport of chlorzoxazone (CZX), a centrally acting muscle relaxant drug, were probed through molecular dynamics simulations. POPC was selected as the model lipid, and three different cholesterol concentrations (0%, 20%, and 50% CHOL) were considered. The outcomes reveal that the area per lipid of POPC decreases and the order parameter increases with enhanced concentration of CHOL. CZX prefers to localize at the interface between the headgroup region and the hydrophobic tail region of POPC, and the main energy barrier occurs in the hydrophobic region. The impact of CHOL on the free energy profile is correlated with concentration: low concentration facilitates CZX permeation, while high concentration hinders CZX permeation. Our findings coincide with experimental results, enhancing the mechanism understanding of how drug molecules are transported through membranes in the presence of CHOL. • The effects caused by cholesterol (CHOL) on the membrane transport of chlorzoxazone (CZX) were studied. • Low CHOL concentration facilitates CZX permeation, while high concentration hinders CZX permeation. • Our findings improve the mechanism understanding of CHOL effects on CZX translocation across membrane.

Functional Group Distributions, Partition Coefficients, and Resistance Factors in Lipid Bilayers Using Site Identification by Ligand Competitive Saturation

Small molecules such as metabolites and drugs must pass through the membrane of the cell, a barrier primarily comprising phospholipid bilayers and embedded proteins. To better understand the process of passive diffusion, knowledge of the ability of various functional groups to partition across bilayers and the associated energetics would be of utility. In the present study, the site identification by ligand competitive saturation (SILCS) methodology has been applied to sample the distributions of a diverse set of chemical solutes representing the functional groups of small molecules across phospholipid bilayers composed of 0.9:0.1 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/cholesterol and a mixture of 0.52:0.18:0.3 1,2-dioleoyl-sn-glycero-3-phospho-l-serine/1,2-dioleoyl-sn-glycero-3-phosphocholine/cholesterol used in parallel artificial membrane permeability assay experiments. A combination of oscillating chemical potential grand canonical Monte Carlo and molecular dynamics in the SILCS simulations was applied to achieve solute sampling through the bilayers and surrounding aqueous environment from which the distribution of solutes and the functional groups they represent were obtained. Results show differential distribution of aliphatic versus aromatic groups with the former having increased sampling in the center of the bilayers versus in the region of the glycerol linker for the latter. Variations in the distribution of different polar groups are evident, with large differences between negative acetate and positive methylammonium with accumulation of the polar-neutral and acetate solutes above the bilayer head groups. Conversion of the distributions to absolute free energies allows for a detailed understanding of energetics of functional groups in different regions of the bilayers and for calculation of absolute free-energy profiles of multifunctional drug-like molecules across the bilayers from which partition coefficients and resistance factors suitable for insertion into the homogenous solubility-diffusion equation for calculation of permeability were obtained. Comparisons of the calculated bilayer/solution partition coefficients with 1-octanol/water experimental data for both drug-like molecules and the solutes show overall good agreement, validating the calculated distributions and associated absolute free-energy profiles.

Mechanical properties of the high cholesterol-containing membrane: An AFM study

Cholesterol (Chol) content in most cellular membranes does not exceed 50 mol%, only in the eye lens's fiber cell plasma membrane, its content surpasses 50 mol%. At this high concentration, Chol induces the formation of pure cholesterol bilayer domains (CBDs), which coexist with the surrounding phospholipid-cholesterol domain (PCD). Here, we applied atomic force microscopy to study the mechanical properties of Chol/phosphatidylcholine membranes where the Chol content was increased from 0 to 75 mol%, relevant to eye lens membranes. The surface roughness of the membrane decreases with an increase of Chol content until it reaches 60 mol%, and roughness increases with a further increment in Chol content. We propose that the increased roughness at higher Chol content results from the formation of CBDs. Force spectroscopy on the membrane with Chol content of 50 mol% or lesser exhibited single breakthrough events, whereas two distinct puncture events were observed for membranes with the Chol content greater than 50 mol%. We propose that the first puncture force corresponds to the membranes containing coexisting PCD and CBDs. In contrast, the second puncture force corresponds to the "CBD water pocket" formed due to coexisting CBDs and PCD. Membrane area compressibility modulus (K_A) increases with an increase in Chol content until it reaches 60 mol%, and with further increment in Chol content, CBDs are formed, and K_A starts to decrease. Our results report the increase in membrane roughness and decrease K_A at very high Chol content (>60 mol%) relevant to the eye lens membrane.

Effects of lipid heterogeneity on model human brain lipid membranes

Cell membranes naturally contain a heterogeneous lipid distribution. However, homogeneous bilayers are commonly preferred and utilised in computer simulations due to their relative simplicity, and the availability of lipid force field parameters. Recently, experimental lipidomics data for the human brain cell membranes under healthy and Alzheimer's disease (AD) conditions were investigated, since disruption to the lipid composition has been implicated in neurodegenerative disorders, including AD [R. B. Chan et al., J. Biol. Chem., 2012, 287, 2678-2688]. In order to observe the effects of lipid complexity on the various bilayer properties, molecular dynamics simulations were used to study four membranes with increasing heterogeneity: a pure POPC membrane, a POPC and cholesterol membrane in a 1 : 1 ratio (POPC-CHOL), and to our knowledge, the first realistic models of a healthy brain membrane and an Alzheimer's diseased brain membrane. Numerous structural, interfacial, and dynamical properties, including the area per lipid, interdigitation, dipole potential, and lateral diffusion of the two simple models, POPC and POPC-CHOL, were analysed and compared to those of the complex brain models consisting of 27 lipid components. As the membranes gain heterogeneity, a number of alterations were found in the structural and dynamical properties, and more significant differences were observed in the lateral diffusion. Additionally, we observed snorkeling behaviour of the lipid tails that may play a role in the permeation of small molecules across biological membranes. In this work, atomistic description of realistic brain membrane models is provided, which can add insight towards the permeability and transport pathways of small molecules across these membrane barriers.

Dynamics and energetics of water transport through aquaporin mutants causing nephrogenic diabetes insipidus (NDI): A molecular dynamics study

Human aquaporin-2 (AQP2) is the principal water channel in the human kidney. Any alteration of its physiological functioning may lead to the water imbalance and consequently diseases in humans, especially nephrogenic diabetes insipidus (NDI). Although many of the mutations associated with NDI are experimentally discovered and examined, a molecular level characterization of the structure and transport mechanism is still missing. In this paper, the gating effects of selectivity filter (SF) as wide/narrow states on the mechanism and dynamics of water permeation within the wild-type AQP2 and two NDI causing mutants as AQP2-V168M and AQP2-G64R are studied for the first time. The analysis of the 200 ns trajectory shows that the SF region in AQP2 is not only a selectivity filter, as previously reported but also it performs as a gating site depending on the side-chain conformation of His172. The assignment of the wide/narrow states of SF is supported by computing the free energy and permeability through the AQP2. Moreover, by exploring the effects of V168M and G64R mutants on the AQP2 structure during 200 ns trajectories, remarkable increases of energy barriers are observed in the middle and cytoplasmic side of the pore, respectively. Interestingly, it is found that due to the variable conformations of the SF region as wide/narrow, the effect of the NDI causing mutants on the average water permeability can be revealed with notably better accuracy by finding the wide states in the wild-type and mutated types of AQP2 and comparing the osmotic permeabilities for this state.Communicated by Ramaswamy H. Sarma.

Highlighting the effect of amyloid beta assemblies on the mechanical properties and conformational stability of cell membrane

Alzheimer disease (AD) is the most common cause of dementia, characterized by a progressive decline in cognitive function due to the abnormal aggregation and deposition of Amyloid beta (Aβ) fibrils in the brain of patients. In this context, the molecular mechanisms of protein misfolding and aggregation that are known to induce significant biophysical alterations in cells, including destabilization of plasma membranes, remain partially unclear. Physical interaction between the Aβ assemblies and the membrane leads to the disruption of the cell membrane in multiple ways including, surface carpeting, generation of transmembrane channels and detergent-like membrane dissolution. Understanding the impact of amyloidogenic protein in different stages of aggregation with the plasma membrane, plays a crucial role to fully elucidate the pathological mechanisms of AD. Within this framework, computer simulations represent a powerful tool able to shed lights on the interactions governing the structural influence of Aβ proteins on biological membrane. In this study, molecular dynamics (MD) simulations have been performed in order to characterize how POPC bilayer conformational and mechanical properties are affected by the interaction with Aβ_11-42 peptide, oligomer and fibril.

A molecular dynamics study proposing the existence of statistical structural heterogeneity due to chain orientation in the POPC-cholesterol bilayer

We performed molecular dynamics simulations of a lipid bilayer consisting of POPC and cholesterol at temperatures from 283 to 308K and cholesterol concentrations from 0 to 50% mol/mol. The purpose of this study was to look for the existence of structural differences in the region delimited by these parameters and, in particular, in a region where coexistence of liquid disordered and liquid ordered phases has been proposed. Our interest in this range of concentration and temperature responds to the fact that polyene ionophore activity varies considerably along it. Two force fields, CHARMM36 and Slipids, were compared in order to determine the most suitable. Both force fields predict non-monotonic behaviors consistent with the existence of phase transitions. We found the presence of lateral structural heterogeneity, statistical in nature, in some of the bilayers occurring in this range of temperatures and sterol concentrations. This heterogeneity was produced by correlated ordering of the POPC tails and not due to cholesterol enrichment, and lasts for tens of nanoseconds. We relate these observations to the action of polyenes in these membranes.

Confocal Microscopy Confirmed that in Phosphatidylcholine Giant Unilamellar Vesicles with very High Cholesterol Content Pure Cholesterol Bilayer Domains Form

The cholesterol (Chol) content in the fiber cell plasma membranes of the eye lens is extremely high, exceeding the solubility threshold in the lenses of old humans. This high Chol content forms pure Chol bilayer domains (CBDs) and Chol crystals in model membranes and membranes formed from the total lipid extracts from human lenses. CBDs have been detected using electron paramagnetic resonance (EPR) spin-labeling approaches. Here, we confirm the presence of CBDs in giant unilamellar vesicles prepared using the electroformation method from Chol/1-palmitoyl-2-oleoylphosphocholine and Chol/distearoylphosphatidylcholine mixtures. Confocal microscopy experiments using phospholipid (PL) analog (1,1'-dioctadecyl-3,3,3',3'-tetramethylindocarbocyanine-5,5'-disulfonic acid) and cholesterol analog fluorescent probes (23-(dipyrrometheneboron difluoride)-24-norcholesterol) were performed, allowing us to make three major conclusions: (1) In all membranes with a Chol/PL mixing ratio (expressed as a molar ratio) >2, pure CBDs were formed within the bulk PL bilayer saturated with Chol. (2) CBDs were present as the pure Chol bilayer and not as separate patches of Chol monolayers in each leaflet of the PL bilayer. (3) CBDs, presented as single large domains, were always located at the top of giant unilamellar vesicles, independent of the change in sample orientation (right-side-up/upside-down). Results obtained with confocal microscopy and fluorescent Chol and PL analogs, combined with those obtained using EPR and spin-labeled Chol and PL analogs, contribute to the understanding of the organization of lipids in the fiber cell plasma membranes of the human eye lens.

Asymmetric Spontaneous Intercalation of Lutein into a Phospholipid Bilayer, a Computational Study

Lutein, a hydroxylated carotenoid, is a pigment synthesised by plants and bacteria. Animals are unable to synthesise lutein, nevertheless, it is present in animal tissues, where its only source is dietary intake. Both in plants and animals, carotenoids are associated mainly with membranes where they carry out important physiological functions. Their trafficking to and insertion into membranes are not well recognised due to experimental difficulties. In this paper, a computational approach is used to elucidate details of the dynamics and energetics of lutein intercalation from the water to the phospholipid bilayer phase. The dynamics is studied using molecular dynamics simulation, and the energetics using umbrella sampling. Lutein spontaneous insertion into the bilayer and translocation across it proceed via formation of hydrogen bonds between its hydroxyl groups and water and/or phospholipid oxygen atoms as well as desolvation of its polyene chain. As lutein molecule is asymmetric, its bilayer intercalation is also asymmetric. The course of events and timescale of the intercalation are different from those of helical peptides. The time of full lutein intercalation ranges from 20 to 100 ns and its final orientation is predominately vertical. Nevertheless, some lutein molecules are in the final horizontal position and some aggregate in the water phase and remain there for the whole simulation time. The highest energy barrier for the intercalation process is ~2.2 kcal/mol and the energy gain is ~18 kcal/mol. The results obtained for lutein can be applied to other xanthophylls and molecules of a similar structure.

•

Lutein as an amphiphilic molecule readily intercalates into a lipid bilayer.

•

MD simulations with dense sampling revealed details of the intercalation process.

•

Time of full lutein membrane intercalation ranges from 20 to 100 ns.

•

Lutein intercalation from the ε ring end is less probable than from the β ring end.

•

Horizontal position of lutein in the bilayer is less probable than vertical.

Transmembrane helix connectivity in Orai1 controls two gates for calcium-dependent transcription

The channel Orai1 requires Ca2+ store depletion in the endoplasmic reticulum and an interaction with the Ca2+ sensor STIM1 to mediate Ca2+ signaling. Alterations in Orai1-mediated Ca2+ influx have been linked to several pathological conditions including immunodeficiency, tubular myopathy, and cancer. We screened large-scale cancer genomics data sets for dysfunctional Orai1 mutants. Five of the identified Orai1 mutations resulted in constitutively active gating and transcriptional activation. Our analysis showed that certain Orai1 mutations were clustered in the transmembrane 2 helix surrounding the pore, which is a trigger site for Orai1 channel gating. Analysis of the constitutively open Orai1 mutant channels revealed two fundamental gates that enabled Ca2+ influx: Arginine side chains were displaced so they no longer blocked the pore, and a chain of water molecules formed in the hydrophobic pore region. Together, these results enabled us to identify a cluster of Orai1 mutations that trigger Ca2+ permeation associated with gene transcription and provide a gating mechanism for Orai1.

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.

Is the cholesterol bilayer domain a barrier to oxygen transport into the eye lens?

In the eye lens, the oxygen partial pressure is very low and the cholesterol (Chol) content in cell membranes is very high. Disturbance of these quantities results in cataract development. In human lens membranes, both bulk phospholipid-Chol domains and the pure Chol bilayer domains (CBDs) were experimentally detected. It is hypothesized that the CBD constitutes a significant barrier to oxygen transport into the lens. Transmembrane profiles of the oxygen diffusion-concentration product, obtained with electron paramagnetic resonance spin-labeling methods, allow evaluation of the oxygen permeability (P_M) of phospholipid membranes but not the CBD. Molecular dynamics simulation can independently provide components of the product across any bilayer domain, thus allowing evaluation of the P_M across the CBD. Therefore, to test the hypothesis, MD simulation was used. Three bilayers containing palmitoyl-oleoyl-phosphorylcholine (POPC) and Chol were built. The pure Chol bilayer modeled the CBD, the 1:1 POPC-Chol bilayer modeled the bulk membrane in which the CBD is embedded, and the POPC bilayer was a reference. To each model, 200 oxygen molecules were added. After equilibration, the oxygen concentration and diffusion profiles were calculated for each model and multiplied by each other. From the respective product profiles, the P_M of each bilayer was calculated. Favorable comparison with experimental data available only for the POPC and POPC-Chol bilayers validated these bilayer models and allowed the conclusion that oxygen permeation across the CBD is ~10 smaller than across the bulk membrane, supporting the hypothesis that the CBD is a barrier to oxygen transport into the eye lens.

Cholesterol Bilayer Domains in the Eye Lens Health: A Review

The most unique biochemical characteristic of the eye lens fiber cell plasma membrane is its extremely high cholesterol content, the need for which is still unclear. It is evident, however, that the disturbance of Chol homeostasis may result in damages associated with cataracts. Electron paramagnetic resonance methods allow discrimination of two types of lipid domains in model membranes overloaded with Chol, namely, phospholipid-cholesterol domains and pure Chol bilayer domains. These domains are also detected in human lens lipid membranes prepared from the total lipids extracted from lens cortices and nuclei of donors from different age groups. Independent of the age-related changes in phospholipid composition, the physical properties of phospholipid-Chol domains remain the same for all age groups and are practically identical for cortical and nuclear membranes. The presence of Chol bilayer domains in these membranes provides a buffering capacity for cholesterol concentration in the surrounding phospholipid-Chol domains, keeping it at a constant saturating level and thus keeping the physical properties of the membrane consistent with and independent of changes in phospholipid composition. It seems that the presence of Chol bilayer domains plays an integral role in the regulation of cholesterol-dependent processes in fiber cell plasm membranes and in the maintenance of fiber cell membrane homeostasis.

High Cholesterol/Low Cholesterol: Effects in Biological Membranes: A Review

Lipid composition determines membrane properties, and cholesterol plays a major role in this determination as it regulates membrane fluidity and permeability, as well as induces the formation of coexisting phases and domains in the membrane. Biological membranes display a very diverse lipid composition, the lateral organization of which plays a crucial role in regulating a variety of membrane functions. We hypothesize that, during biological evolution, membranes with a particular cholesterol content were selected to perform certain functions in the cells of eukaryotic organisms. In this review, we discuss the major membrane properties induced by cholesterol, and their relationship to certain membrane functions.

Glycosylation and Lipids Working in Concert Direct CD2 Ectodomain Orientation and Presentation

Proteins embedded in the plasma membrane mediate interactions with the cell environment and play decisive roles in many signaling events. For cell-cell recognition molecules, it is highly likely that their structures and behavior have been optimized in ways that overcome the limitations of membrane tethering. In particular, the ligand binding regions of these proteins likely need to be maximally exposed. Here we show by means of atomistic simulations of membrane-bound CD2, a small cell adhesion receptor expressed by human T-cells and natural killer cells, that the presentation of its ectodomain is highly dependent on membrane lipids and receptor glycosylation acting in apparent unison. Detailed analysis shows that the underlying mechanism is based on electrostatic interactions complemented by steric interactions between glycans in the protein and the membrane surface. The findings are significant for understanding the factors that render membrane receptors accessible for binding and signaling.

Structural and Functional Effect of an Oscillating Electric Field on the Dopamine-D3 Receptor: A Molecular Dynamics Simulation Study

Dopamine as a neurotransmitter plays a critical role in the functioning of the central nervous system. The structure of D3 receptor as a member of class A G-protein coupled receptors (GPCRs) has been reported. We used MD simulation to investigate the effect of an oscillating electric field, with frequencies in the range 0.6–800 GHz applied along the z-direction, on the dopamine-D3R complex. The simulations showed that at some frequencies, the application of an external oscillating electric field along the z-direction has a considerable effect on the dopamine-D3R. However, there is no enough evidence for prediction of changes in specific frequency, implying that there is no order in changes. Computing the correlation coefficient parameter showed that increasing the field frequency can weaken the interaction between dopamine and D3R and may decrease the Arg128{3.50}-Glu324{6.30} distance. Because of high stability of α helices along the z-direction, applying an oscillating electric field in this direction with an amplitude 10-time higher did not have a considerable effect. However, applying the oscillating field at the frequency of 0.6 GHz along other directions, such as X-Y and Y-Z planes, could change the energy between the dopamine and the D3R, and the number of internal hydrogen bonds of the protein. This can be due to the effect of the direction of the electric field vis-à-vis the ligands orientation and the interaction of the oscillating electric field with the dipole moment of the protein.

Characteristics of a phospholipid DOPC/cholesterol bilayer based on surface free energy and its components

The effect of the cholesterol (Chol) content (x_Chol = 0.25; 0.5; 0.75) on the wettability changes of solid supported bilayers of unsaturated 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) was investigated at 20 °C and 37 °C.

Properties of membranes derived from the total lipids extracted from clear and cataractous lenses of 61-70-year-old human donors

Human lens-lipid membranes prepared from the total lipids extracted from clear and cataractous lens cortexes and nuclei of 61-70-year-old donors by use of a rapid solvent-exchange method were investigated. The measured cholesterol-to-phospholipid (Chol/PL) molar ratio in these membranes was 1.8 and 4.4 for cortex and nucleus of clear lenses, respectively, and 1.14 and 1.45 for cataractous lenses. Properties and organization of the lipid bilayer were investigated by use of electron paramagnetic resonance spin-labeling methods. Formation of Chol crystals was confirmed by use of differential scanning calorimetry. Pure cholesterol bilayer domains (CBDs) were formed in all the membranes investigated. It was shown that in clear lens membranes of the nucleus, Chol exists in three different environments: (1) dispersed in phospholipid bilayers (PCDs), (2) in CBDs, and (3) in Chol crystals. In clear lens membranes of the cortex, and in cortical and nuclear cataractous lens membranes, Chol crystals were not detected, because of the lower Chol content. Profiles of membrane properties (alkyl-chain order, fluidity, oxygen transport, and hydrophobicity) across the PCD were very similar for clear and cataractous membranes. Profiles of the oxygen transport parameter across the CBD were, however, different for cortical clear and cataractous membranes-the amount and size of CBDs was less in cataractous membranes. These results suggest that high Chol content, formation of CBDs, and formation of Chol crystals should not be regarded as major predispositions for the development of age-related cataracts.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

Formation of cholesterol bilayer domains precedes formation of cholesterol crystals in cholesterol/dimyristoylphosphatidylcholine membranes: EPR and DSC studies

Saturation-recovery EPR along with DSC were used to determine the cholesterol content at which pure cholesterol bilayer domains (CBDs) and cholesterol crystals begin to form in dimyristoylphosphatidylcholine (DMPC) membranes. To preserve compositional homogeneity throughout the membrane suspension, lipid multilamellar dispersions were prepared using a rapid solvent exchange method. The cholesterol content increased from 0 to 75 mol %. With spin-labeled cholesterol analogues, it was shown that the CBDs begin to form at ~50 mol % cholesterol. It was confirmed by DSC that the cholesterol solubility threshold for DMPC membranes is detected at ~66 mol % cholesterol. At levels above this cholesterol content, monohydrate cholesterol crystals start to form. The major finding is that the formation of CBDs precedes formation of cholesterol crystals. The region of the phase diagram for cholesterol contents between 50 and 66 mol % is described as a structured one-phase region in which CBDs have to be supported by the surrounding DMPC bilayer saturated with cholesterol. Thus, the phase boundary located at 66 mol % cholesterol separates the structured one-phase region (liquid-ordered phase of DMPC with CBDs) from the two-phase region where the structured liquid-ordered phase of DMPC coexists with cholesterol crystals. It is likely that CBDs are precursors of monohydrate cholesterol crystals.

Comparative computer simulation study of cholesterol in hydrated unary and binary lipid bilayers and in an anhydrous crystal

Models created with molecular dynamics simulations are used to compare the organization and dynamics of cholesterol (Chol) molecules in three different environments: (1) a hydrated pure Chol bilayer that models the Chol bilayer domain, which is a pure Chol domain embedded in the bulk membrane; (2) a 2-palmitoyl-3-oleoyl-d-glycerol-1-phosphorylcholine bilayer saturated with cholesterol (POPC-Chol50) that models the bulk membrane; (3) a Chol crystal. The computer model of the hydrated pure Chol bilayer is stable on the microsecond time scale. Some structural characteristics of Chol molecules in the Chol bilayer are similar to those in the POPC-Chol50 bilayer (e.g., tilt of Chol rings and chains), while others are similar to those in Chol crystals (e.g., surface area per Chol, bilayer thickness). The key result of this study is that the Chol bilayer has, unexpectedly, a dynamic structure, with Chol mobility similar to that in the POPC-Chol50 bilayer though slower. This is the major difference compared to Chol crystals, where Chol molecules are immobile. Also, water accessibility to Chol-OH groups in the Chol bilayer is not limited. On average, each Chol molecule makes 2.3 hydrogen bonds with water in the Chol bilayer, compared with 1.7 hydrogen bonds in the POPC-Col50 bilayer.

Transport of nucleosides in the vcCNT facilitated by sodium gradients from molecular dynamics simulations

Nucleosides are required for DNA and RNA synthesis, and the nucleoside adenosine has a function in a variety of signaling processes. Nucleosides require a specialized class of integral membrane proteins, known as nucleoside transporters (NTs), for specific transport across cell membranes. NTs are also important determinants for the transport of nucleoside-derived drugs across cell membranes. Recently, the crystal structure of the vcCNT (Vibrio cholerae Concentrative Nucleoside Transporter) was reported. Here we perform molecular dynamics (MD) simulations for the vcCNT structure in the presence of various sodium gradients, since CNTs are sodium-coupled transporters. The results highlight the important role of sodium bound to the vcCNT in the transport of uridine. Our MD simulations show that, without NaCl, uridine remains stable in the binding pocket of the vcCNT. In the presence of 20 mM NaCl, uridine moves from the binding pocket and approaches the entrance of the intracellular side. In the presence of 100 mM NaCl, uridine passes through most part of the entrance and approaches the intracellular side. The polar/charged amino acids in the binding pocket are important in the transport process. They first "fix" the ribose and allow the uracil base of uridine to approach the entrance of the intracellular side, and then "release" the ribose to allow uridine to move freely into the intracellular side coupled with the movement of sodium ions and HP1b. Finally, we propose a detailed mechanism of the nucleoside transport from the binding pocket to the intracellular side of the vcCNT.

Properties of fiber cell plasma membranes isolated from the cortex and nucleus of the porcine eye lens

The organization and physical properties of the lipid bilayer portion of intact cortical and nuclear fiber cell plasma membranes isolated from the eye lenses of two-year-old pigs were studied using electron paramagnetic resonance (EPR) spin-labeling. Membrane fluidity, hydrophobicity, and the oxygen transport parameter (OTP) were assessed from the EPR spectra of precisely positioned spin labels. Intact cortical and nuclear membranes, which include membrane proteins, were found to contain three distinct lipid environments. These lipid environments were termed the bulk lipid domain, boundary lipid domain, and trapped lipid domain (lipids in protein aggregates). The amount of boundary and trapped lipids was greater in intact nuclear membranes than in cortical membranes. The properties of intact membranes were compared with the organization and properties of lens lipid membranes made of the total lipid extracts from the lens cortex or nucleus. In cortical lens lipid membranes, only one homogenous environment was detected, which was designated as a bulk lipid domain (phospholipid bilayer saturated with cholesterol). Lens lipid membranes prepared from the lens nucleus possessed two domains, assigned as a bulk lipid domain and a cholesterol bilayer domain (CBD). In intact nuclear membranes, it was difficult to discriminate the CBD, which was clearly detected in nuclear lens lipid membranes, because the OTP measured in the CBD is the same as in the domain formed by trapped lipids. The two domains unique to intact membranes-namely, the domain formed by boundary lipids and the domain formed by trapped lipids-were most likely formed due to the presence of membrane proteins. It is concluded that formation of rigid and practically impermeable domains is enhanced in the lens nucleus, indicating changes in membrane composition that may help to maintain low oxygen concentration in this lens region.

Functions of cholesterol and the cholesterol bilayer domain specific to the fiber-cell plasma membrane of the eye lens

The most unique feature of the eye lens fiber-cell plasma membrane is its extremely high cholesterol content. Cholesterol saturates the bulk phospholipid bilayer and induces formation of immiscible cholesterol bilayer domains (CBDs) within the membrane. Our results (based on EPR spin-labeling experiments with lens-lipid membranes), along with a literature search, have allowed us to identify the significant functions of cholesterol specific to the fiber-cell plasma membrane, which are manifest through cholesterol-membrane interactions. The crucial role is played by the CBD. The presence of the CBD ensures that the surrounding phospholipid bilayer is saturated with cholesterol. The saturating cholesterol content in fiber-cell membranes keeps the bulk physical properties of lens-lipid membranes consistent and independent of changes in phospholipid composition. Thus, the CBD helps to maintain lens-membrane homeostasis when the membrane phospholipid composition changes significantly. The CBD raises the barrier for oxygen transport across the fiber-cell membrane, which should help to maintain a low oxygen concentration in the lens interior. It is hypothesized that the appearance of the CBD in the fiber-cell membrane is controlled by the phospholipid composition of the membrane. Saturation with cholesterol smoothes the phospholipid-bilayer surface, which should decrease light scattering and help to maintain lens transparency. Other functions of cholesterol include formation of hydrophobic and rigidity barriers across the bulk phospholipid-cholesterol domain and formation of hydrophobic channels in the central region of the membrane for transport of small, nonpolar molecules parallel to the membrane surface. In this review, we provide data supporting these hypotheses.