Cholesterol Promotes Protein Binding by Affecting Membrane Electrostatics and Solvation Properties

HIV-1 Gag Polyprotein Affinity to the Lipid Membrane Is Independent of Its Surface Charge

The binding of the HIV-1 Gag polyprotein to the plasma membrane is a critical step in viral replication. The association with membranes depends on the lipid composition, but its mechanisms remain unclear. Here, we report the binding of non-myristoylated Gag to lipid membranes of different lipid compositions to dissect the influence of each component. We tested the contribution of phosphatidylserine, PI(4,5)P2, and cholesterol to membrane charge density and Gag affinity to membranes. Taking into account the influence of the membrane surface potential, we quantitatively characterized the adsorption of the protein onto model lipid membranes. The obtained Gag binding constants appeared to be the same regardless of the membrane charge. Furthermore, Gag adsorbed on uncharged membranes, suggesting a contribution of hydrophobic forces to the protein-lipid interaction. Charge-charge interactions resulted in an increase in protein concentration near the membrane surface. Lipid-specific interactions were observed in the presence of cholesterol, resulting in a two-fold increase in binding constants. The combination of cholesterol with PI(4,5)P2 showed cooperative effects on protein adsorption. Thus, we suggest that the affinity of Gag to lipid membranes results from a combination of electrostatic attraction to acidic lipids, providing different protein concentrations near the membrane surface, and specific hydrophobic interactions.

Small-angle X-ray and neutron scattering applied to lipid-based nanoparticles: Recent advancements across different length scales

Lipid-based nanoparticles (LNPs), ranging from nanovesicles to non-lamellar assemblies, have gained significant attention in recent years, as versatile carriers for delivering drugs, vaccines, and nutrients. Small-angle scattering methods, employing X-rays (SAXS) or neutrons (SANS), represent unique tools to unveil structure, dynamics, and interactions of such particles on different length scales, spanning from the nano to the molecular scale. This review explores the state-of-the-art on scattering methods applied to unveil the structure of lipid-based nanoparticles and their interactions with drugs and bioactive molecules, to inform their rational design and formulation for medical applications. We will focus on complementary information accessible with X-rays or neutrons, ranging from insights on the structure and colloidal processes at a nanoscale level (SAXS) to details on the lipid organization and molecular interactions of LNPs (SANS). In addition, we will review new opportunities offered by Time-resolved (TR)-SAXS and -SANS for the investigation of dynamic processes involving LNPs. These span from real-time monitoring of LNPs structural evolution in response to endogenous or external stimuli (TR-SANS), to the investigation of the kinetics of lipid diffusion and exchange upon interaction with biomolecules (TR-SANS). Finally, we will spotlight novel combinations of SAXS and SANS with complementary on-line techniques, recently enabled at Large Scale Facilities for X-rays and neutrons. This emerging technology enables synchronized multi-method investigation, offering exciting opportunities for the simultaneous characterization of the structure and chemical or mechanical properties of LNPs.

Impact of Lipid Composition on Vesicle Protein Adsorption: A BSA Case Study

Investigating the interaction between liposomes and proteins is of paramount importance in the development of liposomal formulations with real potential for bench-to-bedside transfer. Upon entering the body, proteins are immediately adsorbed on the liposomal surface, changing the nanovehicles' biological identity, which has a significant impact on their biodistribution and pharmacokinetics and ultimately on their therapeutic effect. Albumin is the most abundant plasma protein and thus usually adsorbs immediately on the liposomal surface. We herein report a comprehensive investigation on the adsorption of model protein bovine serum albumin (BSA) onto liposomal vesicles containing the zwitterionic lipid 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), in combination with either cholesterol (CHOL) or the cationic lipid 1,2-dioleoyl-3-trimethylammoniumpropane (DOTAP). While many studies regarding protein adsorption on the surface of liposomes with different compositions have been performed, to the best of our knowledge, the differential responses of CHOL and DOTAP upon albumin adsorption on vesicles have not yet been investigated. UV-vis spectroscopy and sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) revealed a strong influence of the phospholipid membrane composition on protein adsorption. Hence, it was found that DOTAP-containing vesicles adsorb proteins more robustly but also aggregate in the presence of BSA, as confirmed by DLS and TEM. Separation of liposome-protein complexes from unadsorbed proteins performed by means of centrifugation and size exclusion chromatography (SEC) was also investigated. Our results show that neither method can be regarded as a golden experimental setup to study the protein corona of liposomes. Yet, SEC proved to be more successful in the separation of unbound proteins, although the amount of lipid loss upon liposome elution was higher than expected. In addition, coarse-grained molecular dynamics simulations were employed to ascertain key membrane parameters, such as the membrane thickness and area per lipid. Overall, this study highlights the importance of surface charge and membrane fluidity in influencing the extent of protein adsorption. We hope that our investigation will be a valuable contribution to better understanding protein-vesicle interactions for improved nanocarrier design.

Efficient calculation of orientation-dependent lipid dynamics from membrane simulations

Molecular dynamics simulations of lipid membranes have become increasingly impactful in biophysics because they offer atomistic resolution of structural fluctuations in relation to their functional outputs. Yet quantitative characterization of multiscale processes is a formidable challenge due to the distribution of motions that evade analysis of discrete simulation data. Here we investigate the efficient calculation of CH bond relaxation rates from membrane simulations. Widely used computational approaches offer numerical simplicity but fall short of capturing crucial aspects of the orientation dependence of the dynamics. To circumvent this problem, we introduced a robust framework based on liquid crystal theory which considers explicitly the CH bond motions with respect to the director axis (bilayer normal). Analysis of the orientation dependence of the dynamics shows excellent agreement with experiment, illustrating how the ordering potential affects the calculated relaxation rates. Furthermore, a fit-based resampling of the autocorrelation function of the bond fluctuations validates the new approach for low-temporal resolution data. The recovered relaxation rates indicate that at short timescales, both with and without cholesterol, the local motions of CH bonds describe the bilayer microviscosity and resemble liquid hydrocarbons. Our results establish the critical role of the orientational anisotropy in analysis of membrane simulations, explain fundamental aspects of lipid dynamics, and provide guidelines for extracting information that can be compared to experimental data.

A myristoyl switch at the plasma membrane triggers cleavage and oligomerization of Mason-Pfizer monkey virus matrix protein

For most retroviruses, including HIV, association with the plasma membrane (PM) promotes the assembly of immature particles, which occurs simultaneously with budding and maturation. In these viruses, maturation is initiated by oligomerization of polyprotein precursors. In contrast, several retroviruses, such as Mason-Pfizer monkey virus (M-PMV), assemble in the cytoplasm into immature particles that are transported across the PM. Therefore, protease activation and specific cleavage must not occur until the pre-assembled particle interacts with the PM. This interaction is triggered by a bipartite signal consisting of a cluster of basic residues in the matrix (MA) domain of Gag polyprotein and a myristoyl moiety N-terminally attached to MA. Here, we provide evidence that myristoyl exposure from the MA core and its insertion into the PM occurs in M-PMV. By a combination of experimental methods, we show that this results in a structural change at the C-terminus of MA allowing efficient cleavage of MA from the downstream region of Gag. This suggests that, in addition to the known effect of the myristoyl switch of HIV-1 MA on the multimerization state of Gag and particle assembly, the myristoyl switch may have a regulatory role in initiating sequential cleavage of M-PMV Gag in immature particles.

Antiviral efficacy of nanomaterial-treated textiles in real-life like exposure conditions

Due to the growing interest towards reducing the number of potentially infectious agents on critical high-touch surfaces, the popularity of antimicrobially and antivirally active surfaces, including textiles, has increased. The goal of this study was to create antiviral textiles by spray-depositing three different nanomaterials, two types of CeO2 nanoparticles and quaternary ammonium surfactant CTAB loaded SiO2 nanocontainers, onto the surface of a knitted polyester textile and assess their antiviral activity against two coronaviruses, porcine transmissible gastroenteritis virus (TGEV) and severe acute respiratory syndrome virus (SARS CoV-2). Antiviral testing was carried out in small droplets in semi-dry conditions and in the presence of organic soiling, to mimic aerosol deposition of viruses onto the textiles. In such conditions, SARS CoV-2 stayed infectious at least for 24 h and TGEV infected cells even after 72h of semi-dry deposition suggesting that textiles exhibiting sufficient antiviral activity before or at 24 h, can be considered promising. The antiviral efficacy of nanomaterial-deposited textiles was compared with the activity of the same nanomaterials in colloidal form and with positive control textiles loaded with copper nitrate and CTAB. Our results indicated that after deposition onto the textile, CeO2 nanoparticles lost most of their antiviral activity, but antiviral efficacy of CTAB-loaded SiO2 nanocontainers was retained also after deposition. Copper nitrate deposited textile that was used as a positive control, showed relatively high antiviral activity as expected. However, as copper was effectively washed away from the textile already during 1 h, the use of copper for creating antiviral textiles would be impractical. In summary, our results indicated that antiviral activity of textiles cannot be predicted from antiviral efficacy of the deposited compounds in colloid and attention should be paid on prolonged efficacy of antivirally coated textiles.

Graph-Based Analyses of Dynamic Water-Mediated Hydrogen-Bond Networks in Phosphatidylserine: Cholesterol Membranes

Phosphatidylserine lipids are anionic molecules present in eukaryotic plasma membranes, where they have essential physiological roles. The altered distribution of phosphatidylserine in cells such as apoptotic cancer cells, which, unlike healthy cells, expose phosphatidylserine, is of direct interest for the development of biomarkers. We present here applications of a recently implemented Depth-First-Search graph algorithm to dissect the dynamics of transient water-mediated lipid clusters at the interface of a model bilayer composed of 1-palmytoyl-2-oleoyl-sn-glycero-2-phosphatidylserine (POPS) and cholesterol. Relative to a reference POPS bilayer without cholesterol, in the POPS:cholesterol bilayer there is a somewhat less frequent sampling of relatively complex and extended water-mediated hydrogen-bond networks of POPS headgroups. The analysis protocol used here is more generally applicable to other lipid:cholesterol bilayers.

Building Asymmetric Lipid Bilayers for Molecular Dynamics Simulations: What Methods Exist and How to Choose One?

The compositional asymmetry of biological membranes has attracted significant attention over the last decade. Harboring more differences from symmetric membranes than previously appreciated, asymmetric bilayers have proven quite challenging to study with familiar concepts and techniques, leaving many unanswered questions about the reach of the asymmetry effects. One particular area of active research is the computational investigation of composition- and number-asymmetric lipid bilayers with molecular dynamics (MD) simulations. Offering a high level of detail into the organization and properties of the simulated systems, MD has emerged as an indispensable tool in the study of membrane asymmetry. However, the realization that results depend heavily on the protocol used for constructing the asymmetric bilayer models has sparked an ongoing debate about how to choose the most appropriate approach. Here we discuss the underlying source of the discrepant results and review the existing methods for creating asymmetric bilayers for MD simulations. Considering the available data, we argue that each method is well suited for specific applications and hence there is no single best approach. Instead, the choice of a construction protocol-and consequently, its perceived accuracy-must be based primarily on the scientific question that the simulations are designed to address.

Effect of cholesterol on the ion-membrane interaction: Zeta potential and dynamic light scattering study

Cholesterol in a bio-membrane plays a significant role in many cellular event and is known to regulate the functional activity of protein and ion channel. In this study we report a significant effect of cholesterol on the ion-membrane interaction. We prepare large unilamellar vesicles, composed of zwitterionic lipid DOPC and anionic lipid DOPG with different cholesterol concentration. Electrostatics of anionic membranes containing cholesterol in the presence of NaCl has systematically been explored using dynamic light scattering and zeta potential. Negative zeta potential of the membrane decreases its negative value with increasing ion concentration for all cholesterol concentrations. However, zeta potential itself decreases with increasing cholesterol content even in the absence of monovalent ions. Electrostatic behaviour of the membrane is determined from well-known Gouy Chapmann model. Negative surface charge density of the membrane decreases with increasing cholesterol content. Binding constant, estimated from the electrostatic double layer theory, is found to increase significantly in the presence of cholesterol. Comparison of electrostatic parameters of the membrane in the presence and absence of cholesterol suggests that cholesterol significantly alter the electrostatic behaviour of the membrane.

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent κC modulus for bilayer bending. Our results show a strong correlation with κC values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

Cholesterol promotes clustering of PI(4,5)P2 driving unconventional secretion of FGF2

FGF2 is a cell survival factor involved in tumor-induced angiogenesis that is secreted through an unconventional secretory pathway based upon direct protein translocation across the plasma membrane. Here, we demonstrate that both PI(4,5)P2-dependent FGF2 recruitment at the inner plasma membrane leaflet and FGF2 membrane translocation into the extracellular space are positively modulated by cholesterol in living cells. We further revealed cholesterol to enhance FGF2 binding to PI(4,5)P2-containing lipid bilayers. Based on extensive atomistic molecular dynamics (MD) simulations and membrane tension experiments, we proposed cholesterol to modulate FGF2 binding to PI(4,5)P2 by (i) increasing head group visibility of PI(4,5)P2 on the membrane surface, (ii) increasing avidity by cholesterol-induced clustering of PI(4,5)P2 molecules triggering FGF2 oligomerization, and (iii) increasing membrane tension facilitating the formation of lipidic membrane pores. Our findings have general implications for phosphoinositide-dependent protein recruitment to membranes and explain the highly selective targeting of FGF2 toward the plasma membrane, the subcellular site of FGF2 membrane translocation during unconventional secretion of FGF2.

Small-Angle Neutron Scattering for Studying Lipid Bilayer Membranes

Small-angle neutron scattering (SANS) is a powerful tool for studying biological membranes and model lipid bilayer membranes. The length scales probed by SANS, being from 1 nm to over 100 nm, are well-matched to the relevant length scales of the bilayer, particularly when it is in the form of a vesicle. However, it is the ability of SANS to differentiate between isotopes of hydrogen as well as the availability of deuterium labeled lipids that truly enable SANS to reveal details of membranes that are not accessible with the use of other techniques, such as small-angle X-ray scattering. In this work, an overview of the use of SANS for studying unilamellar lipid bilayer vesicles is presented. The technique is briefly presented, and the power of selective deuteration and contrast variation methods is discussed. Approaches to modeling SANS data from unilamellar lipid bilayer vesicles are presented. Finally, recent examples are discussed. While the emphasis is on studies of unilamellar vesicles, examples of the use of SANS to study intact cells are also presented.

Effects of cholesterol on the anionic magnetite nanoparticle-induced deformation and poration of giant lipid vesicles

We have investigated the effects of cholesterol on the deformation and poration of giant unilamellar vesicles (GUVs) induced by anionic magnetite nanoparticles (NPs). Negatively charged lipid, neutral lipid, and cholesterol were used to prepare the charged GUVs (surface charge density of membranes - 0.16 C m-2), while only neutral lipid and cholesterol were used to prepare the neutral GUVs. Cholesterol content varied from 0 to 40 mole% for preparing the biologically relevant membranes. The degree of deformation has been characterized by compactness, the value of which remains at 1.0 for spherical GUVs. The value of compactness increases with time for both membranes, but this increase depends on cholesterol content. The average compactness decreases with cholesterol content, and at 60 min, the values are 1.280 ± 0.002 and 1.131 ± 0.010 for 0 and 40 mole% cholesterol containing charged GUVs. The average compactness is relatively lower for neutral GUVs for the corresponding cholesterol. Membrane poration has been investigated by the leakage of calcein, which indicates a two-state transition model. The fraction of deformation is higher for charged GUVs than for neutral ones, while the fraction of poration shows the opposite result. Both the fractions decrease with cholesterol content.

Measles and Nipah virus assembly: Specific lipid binding drives matrix polymerization

Measles virus, Nipah virus, and multiple other paramyxoviruses cause disease outbreaks in humans and animals worldwide. The paramyxovirus matrix (M) protein mediates virion assembly and budding from host cell membranes. M is thus a key target for antivirals, but few high-resolution structures of paramyxovirus M are available, and we lack the clear understanding of how viral M proteins interact with membrane lipids to mediate viral assembly and egress that is needed to guide antiviral design. Here, we reveal that M proteins associate with phosphatidylserine and phosphatidylinositol 4,5-bisphosphate [PI(4,5)P2] at the plasma membrane. Using x-ray crystallography, electron microscopy, and molecular dynamics, we demonstrate that PI(4,5)P2 binding induces conformational and electrostatic changes in the M protein surface that trigger membrane deformation, matrix layer polymerization, and virion assembly.

Effects of cholesterol on the size distribution and bending modulus of lipid vesicles

The influence of cholesterol fraction in the membranes of giant unilamellar vesicles (GUVs) on their size distributions and bending moduli has been investigated. The membranes of GUVs were synthesized by a mixture of two elements: electrically neutral lipid 1, 2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and cholesterol and also a mixture of three elements: electrically charged lipid 1,2-dioleoyl-sn-glycero-3-phospho-(1′-rac-glycerol) (DOPG), DOPC and cholesterol. The size distributions of GUVs have been presented by a set of histograms. The classical lognormal distribution is well fitted to the histograms, from where the average size of vesicle is obtained. The increase of cholesterol content in the membranes of GUVs increases the average size of vesicles in the population. Using the framework of Helmholtz free energy of the system, the theory developed by us is extended to explain the experimental results. The theory determines the influence of cholesterol on the bending modulus of membranes from the fitting of the proper histograms. The increase of cholesterol in GUVs increases both the average size of vesicles in population and the bending modulus of membranes.

Interaction of Spike protein and lipid membrane of SARS-CoV-2 with Ursodeoxycholic acid, an in-silico analysis

Numerous repositioned drugs have been sought to decrease the severity of SARS-CoV-2 infection. It is known that among its physicochemical properties, Ursodeoxycholic Acid (UDCA) has a reduction in surface tension and cholesterol solubilization, it has also been used to treat cholesterol gallstones and viral hepatitis. In this study, molecular docking was performed with the SARS-CoV-2 Spike protein and UDCA. In order to confirm this interaction, we used Molecular Dynamics (MD) in “SARS-CoV-2 Spike protein-UDCA”. Using another system, we also simulated MD with six UDCA residues around the Spike protein at random, naming this “SARS-CoV-2 Spike protein-6UDCA”. Finally, we evaluated the possible interaction between UDCA and different types of membranes, considering the possible membrane conformation of SARS-CoV-2, this was named “SARS-CoV-2 membrane-UDCA”. In the “SARS-CoV-2 Spike protein-UDCA”, we found that UDCA exhibits affinity towards the central region of the Spike protein structure of − 386.35 kcal/mol, in a region with 3 alpha helices, which comprises residues from K986 to C1032 of each monomer. MD confirmed that UDCA remains attached and occasionally forms hydrogen bonds with residues R995 and T998. In the presence of UDCA, we observed that the distances between residues atoms OG1 and CG2 of T998 in the monomers A, B, and C in the prefusion state do not change and remain at 5.93 ± 0.62 and 7.78 ± 0.51 Å, respectively, compared to the post-fusion state. Next, in “SARS-CoV-2 Spike protein-6UDCA”, the three UDCA showed affinity towards different regions of the Spike protein, but only one of them remained bound to the region between the region's heptad repeat 1 and heptad repeat 2 (HR1 and HR2) for 375 ps of the trajectory. The RMSD of monomer C was the smallest of the three monomers with a value of 2.89 ± 0.32, likewise, the smallest RMSF was also of the monomer C (2.25 ± 056). In addition, in the simulation of “SARS-CoV-2 membrane-UDCA”, UDCA had a higher affinity toward the virion-like membrane; where three of the four residues remained attached once they were close (5 Å, to the centre of mass) to the membrane by 30 ns. However, only one of them remained attached to the plasma-like membrane and this was in a cluster of cholesterol molecules. We have shown that UDCA interacts in two distinct regions of Spike protein sequences. In addition, UDCA tends to stay bound to the membrane, which could potentially reduce the internalization of SARS-CoV-2 in the host cell.

An infectious Rous Sarcoma Virus Gag mutant that is defective in nuclear cycling

During retroviral replication, unspliced viral genomic RNA (gRNA) must escape the nucleus for translation into viral proteins and packaging into virions. "Complex" retroviruses such as Human Immunodeficiency Virus (HIV) use cis-acting elements on the unspliced gRNA in conjunction with trans-acting viral proteins to facilitate this escape. "Simple" retroviruses such as Mason-Pfizer Monkey Virus (MPMV) and Murine Leukemia Virus (MLV) exclusively use cis-acting elements on the gRNA in conjunction with host nuclear export proteins for nuclear escape. Uniquely, the simple retrovirus Rous Sarcoma Virus (RSV) has a Gag structural protein that cycles through the nucleus prior to plasma membrane binding. This trafficking has been implicated in facilitating gRNA nuclear export and is thought to be a required mechanism. Previously described mutants that abolish nuclear cycling displayed enhanced plasma membrane binding, enhanced virion release, and a significant loss in genome incorporation resulting in loss of infectivity. Here, we describe a nuclear cycling deficient RSV Gag mutant that has similar plasma membrane binding and genome incorporation to WT virus and surprisingly, is replication competent albeit with a slower rate of spread compared to WT. This mutant suggests that RSV Gag nuclear cycling is not strictly required for RSV replication. Importance While mechanisms for retroviral Gag assembly at the plasma membrane are beginning to be characterized, characterization of intermediate trafficking locales remain elusive. This is in part due to the difficulty of tracking individual proteins from translation to plasma membrane binding. RSV Gag nuclear cycling is a unique phenotype that may provide comparative insight to viral trafficking evolution and may present a model intermediate to cis- and trans-acting mechanisms for gRNA export.

Structural Insights into the Mechanism of Human T-cell Leukemia Virus Type 1 Gag Targeting to the Plasma Membrane for Assembly

Retroviral Gag targeting to the plasma membrane (PM) for assembly is mediated by the N-terminal matrix (MA) domain. For many retroviruses, Gag-PM interaction is dependent on phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂). However, it has been shown that for human T-cell leukemia virus type 1 (HTLV-1), Gag binding to membranes is less dependent on PI(4,5)P₂ than HIV-1, suggesting that other factors may modulate Gag assembly. To elucidate the mechanism by which HTLV-1 Gag binds to the PM, we employed NMR techniques to determine the structure of unmyristoylated MA (myr(-)MA) and to characterize its interactions with lipids and liposomes. The MA structure consists of four α-helices and unstructured N- and C-termini. We show that myr(-)MA binds to PI(4,5)P₂ via the polar head and that binding to inositol phosphates (IPs) is significantly enhanced by increasing the number of phosphate groups on the inositol ring, indicating that the MA-IP binding is governed by charge-charge interactions. The IP binding site was mapped to a well-defined basic patch formed by lysine and arginine residues. Using an NMR-based liposome binding assay, we show that PI(4,5)P₂and phosphatidylserine enhance myr(-)MA binding in a synergistic fashion. Confocal microscopy data revealed formation of puncta on the PM of Gag expressing cells. However, G2A-Gag mutant, lacking myristoylation, is diffuse and cytoplasmic. These results suggest that although myr(-)MA binds to membranes, myristoylation appears to be key for formation of HTLV-1 Gag puncta on the PM. Altogether, these findings advance our understanding of a key mechanism in retroviral assembly.

Binding of Ca2+-independent C2 domains to lipid membranes: A multi-scale molecular dynamics study

C2 domains facilitate protein interactions with lipid bilayers in either a Ca²⁺-dependent or -independent manner. We used molecular dynamics (MD) simulations to explore six Ca²⁺-independent C2 domains, from KIBRA, PI3KC2α, RIM2, PTEN, SHIP2, and Smurf2. In coarse-grained MD simulations these C2 domains formed transient interactions with zwitterionic bilayers, compared with longer-lived interactions with anionic bilayers containing phosphatidylinositol bisphosphate (PIP₂). Type I C2 domains bound non-canonically via the front, back, or side of the β sandwich, whereas type II C2 domains bound canonically, via the top loops. C2 domains interacted strongly with membranes containing PIP₂, causing bound anionic lipids to cluster around the protein. Binding modes were refined via atomistic simulations. For PTEN and SHIP2, CG simulations of their phosphatase plus C2 domains with PIP₂-containing bilayers were also performed, and the roles of the two domains in membrane localization compared. These studies establish a simulation protocol for membrane-recognition proteins.

•

Binding of Ca²⁺-independent C2 domains to membranes was explored by MD simulation

•

C2 domains from KIBRA, PI3KC2α, RIM2, PTEN, SHIP2, and Smurf2 were compared

•

C2 domains formed longer-lived interactions with lipid bilayers containing PIP₂

•

For PTEN and SHIP2, simulations of their phosphatase plus C2 domains were performed

Unusual mode of dimerization of retinitis pigmentosa-associated F220C rhodopsin

Mutations in the G protein-coupled receptor (GPCR) rhodopsin are a common cause of autosomal dominant retinitis pigmentosa, a blinding disease. Rhodopsin self-associates in the membrane, and the purified monomeric apo-protein opsin dimerizes in vitro as it transitions from detergent micelles to reconstitute into a lipid bilayer. We previously reported that the retinitis pigmentosa-linked F220C opsin mutant fails to dimerize in vitro, reconstituting as a monomer. Using fluorescence-based assays and molecular dynamics simulations we now report that whereas wild-type and F220C opsin display distinct dimerization propensities in vitro as previously shown, they both dimerize in the plasma membrane of HEK293 cells. Unexpectedly, molecular dynamics simulations show that F220C opsin forms an energetically favored dimer in the membrane when compared with the wild-type protein. The conformation of the F220C dimer is unique, with transmembrane helices 5 and 6 splayed apart, promoting widening of the intracellular vestibule of each protomer and influx of water into the protein interior. FRET experiments with SNAP-tagged wild-type and F220C opsin expressed in HEK293 cells are consistent with this conformational difference. We speculate that the unusual mode of dimerization of F220C opsin in the membrane may have physiological consequences.

Role of the Lipid Membrane and Membrane Proteins in Tau Pathology

Abnormal accumulation of misfolded tau aggregates is a pathological hallmark of various tauopathies including Alzheimer’s disease (AD). Although tau is a cytosolic microtubule-associated protein enriched in neurons, it is also found in extracellular milieu, such as interstitial fluid, cerebrospinal fluid, and blood. Accumulating evidence showed that pathological tau spreads along anatomically connected areas in the brain through intercellular transmission and templated misfolding, thereby inducing neurodegeneration and cognitive dysfunction. In line with this, the spatiotemporal spreading of tau pathology is closely correlated with cognitive decline in AD patients. Although the secretion and uptake of tau involve multiple different pathways depending on tau species and cell types, a growing body of evidence suggested that tau is largely secreted in a vesicle-free forms. In this regard, the interaction of vesicle-free tau with membrane is gaining growing attention due to its importance for both of tau secretion and uptake as well as aggregation. Here, we review the recent literature on the mechanisms of the tau-membrane interaction and highlights the roles of lipids and proteins at the membrane in the tau-membrane interaction as well as tau aggregation.

Structural and functional consequences of reversible lipid asymmetry in living membranes

Maintenance of lipid asymmetry across the two leaflets of the plasma membrane (PM) bilayer is a ubiquitous feature of eukaryotic cells. Loss of this asymmetry has been widely associated with cell death. However, increasing evidence points to the physiological importance of non-apoptotic, transient changes in PM asymmetry. Such transient scrambling events are associated with a range of biological functions, including intercellular communication and intracellular signaling. Thus, regulation of interleaflet lipid distribution in the PM is a broadly important but underappreciated cellular process with key physiological and structural consequences. Here, we compile the mounting evidence revealing multifaceted, functional roles of PM asymmetry and transient loss thereof. We discuss the consequences of reversible asymmetry on PM structure, biophysical properties and interleaflet coupling. We argue that despite widespread recognition of broad aspects of membrane asymmetry, its importance in cell biology demands more in-depth investigation of its features, regulation, and physiological and pathological implications.

The influence of myristoylation, liposome surface charge and nucleic acid interaction in the partition properties of HIV-1 Gag-N-terminal peptides to membranes

The group-specific antigen (GAG) polyprotein of HIV-1 is the main coordinator of the virus assembly process at the plasma membrane (PM) and is directed by its N-terminal matrix domain (MA). MA is myristoylated and possess a highly basic region (HBR) responsible for the interaction with the negative lipids of the PM, especially with PIP₂. In addition, MA binds RNA molecules proposed as a regulatory step of the assembly process. Here we study the interaction of a synthetic peptide (N-terminal 21 amino acids of MA) and liposomes of different compositions using a variety of biophysical techniques. Particularly, we use the fluorescence properties of the single tryptophan of the peptide to analyze its partition to membranes, where we harness for first time the analytical ability of spectral phasors method to study this interaction. We found that electrostatic interactions play an important role for peptide partition to membranes and myristoylation reduces the free energy of the process. Interestingly, we observe that while the presence of PIP₂ does not cause measurable changes on the peptide-membrane interaction, the interaction is favored by cholesterol. Additionally, we found that the partition process goes through a transition state involving peptide disaggregation and changes in the peptide secondary structure. On the other hand, we found that the presence of oligonucleotides competes with the interaction with lipids by increasing peptide solubility. In summary, we think that our results, in context of the current knowledge of the role of HIV-1 MA, contribute to a better molecular understanding of the membrane association process.

Molecular Structure of Sphingomyelin in Fluid Phase Bilayers Determined by the Joint Analysis of Small-Angle Neutron and X-ray Scattering Data

We have determined the fluid bilayer structure of palmitoyl sphingomyelin (PSM) and stearoyl sphingomyelin (SSM) by simultaneously analyzing small-angle neutron and X-ray scattering data. Using a newly developed scattering density profile (SDP) model for sphingomyelin lipids, we report structural parameters including the area per lipid, total bilayer thickness, and hydrocarbon thickness, in addition to lipid volumes determined by densitometry. Unconstrained all-atom simulations of PSM bilayers at 55 °C using the C36 CHARMM force field produced a lipid area of 56 Ų, a value that is 10% lower than the one determined experimentally by SDP analysis (61.9 Ų). Furthermore, scattering form factors calculated from the unconstrained simulations were in poor agreement with experimental form factors, even though segmental order parameter (S_CD) profiles calculated from the simulations were in relatively good agreement with S_CD profiles obtained from NMR experiments. Conversely, constrained area simulations at 61.9 Ų resulted in good agreement between the simulation and experimental scattering form factors, but not with S_CD profiles from NMR. We discuss possible reasons for the discrepancies between these two types of data that are frequently used as validation metrics for molecular dynamics force fields.

Phospholipid packing defects and oxysterols in atherosclerosis: Dietary prevention and the French paradox

The research literature on atherosclerosis includes findings investigating the atherosclerotic effect of oxysterols, which are the oxidation products of cholesterol; and the literature on oxysterols refers to mechanisms by which oxysterols cause phospholipid packing defects in cell membranes. This review synthesizes these two bodies of research findings to describe how oxysterols cause phospholipid packing defects within the membranes of vascular endothelial cells, potentially increasing cell permeability of low-density lipoprotein cholesterol which may lead to atheroma formation. Exogenous sources of oxysterols are provided by dietary intake of animal-based foods that contain cholesterol oxidation products. This review proposes an explanation for the anti-atherosclerotic effect of plant-based dietary patterns, which is attributed to restriction or avoidance of dietary oxysterol intake from animal-based foods. Furthermore, raw-milk cheeses play an important role in the traditional French diet-low oxysterol content in these unheated foods may contribute to the French paradox, in which reduced coronary heart disease is associated with a diet high in saturated fat and cholesterol.

HIV-1 Gag specifically restricts PI(4,5)P2 and cholesterol mobility in living cells creating a nanodomain platform for virus assembly

HIV-1 Gag protein assembles at the plasma membrane of infected cells for viral particle formation. Gag targets lipids, mainly PI(4,5)P2, at the inner leaflet of this membrane. Here, we address the question whether Gag is able to trap specifically PI(4,5)P2 or other lipids during HIV-1 assembly in the host CD4+ T lymphocytes. Lipid dynamics within and away from HIV-1 assembly sites were determined using super-resolution microscopy coupled with scanning fluorescence correlation spectroscopy in living cells. Analysis of HIV-1-infected cells revealed that, upon assembly, HIV-1 is able to specifically trap PI(4,5)P2 and cholesterol, but not phosphatidylethanolamine or sphingomyelin. Furthermore, our data showed that Gag is the main driving force to restrict the mobility of PI(4,5)P2 and cholesterol at the cell plasma membrane. This is the first direct evidence highlighting that HIV-1 creates its own specific lipid environment by selectively recruiting PI(4,5)P2 and cholesterol as a membrane nanoplatform for virus assembly.

Concurrent Compression of Phospholipid Membranes by Calcium and Cholesterol

Regulation of cell metabolism, membrane fusion, association of proteins with cellular membranes, and cellular signaling altogether would not be possible without Ca²⁺ ions. The distribution of calcium within the cell is uneven with the negatively charged inner leaflet of the plasma membrane being one of the primary targets of its accumulation. Therefore, we decided to map the influence of Ca²⁺ on the properties of lipid bilayers closely resembling natural lipid membranes. We combined fluorescence spectroscopy (analysis of time-resolved emission spectra of Laurdan probe and derived parameters: integrated relaxation time related to local lipid mobility, and total emission shift reflecting membrane polarity and hydration) with molecular dynamics simulations to determine the effect of the increasing CaCl₂ concentration on model lipid membranes containing POPC, POPS, and cholesterol. On top of that, the impact of calcium on the plasma membranes isolated from HEK293 cells was investigated using the steady-state fluorescence of Laurdan. We found that calcium increases rigidity of all the model lipid membranes used, elevates their thickness, increases lipid packing and ordering, and impedes the local lipid mobility. All these effects were to a great extent similar to those elicited by cholesterol. However, the changes of the membrane properties induced by calcium and cholesterol seem largely independent from each other. At sufficiently high concentrations of calcium or cholesterol, the steric effects hindered a further alteration of membrane organization, i.e., the compressibility limit of membrane structures was reached. We found no indication for mutual interaction between Ca²⁺ and cholesterol, nor competition of Ca²⁺ ions and hydroxyl groups of cholesterol for binding to phospholipids. Fluorescence measurements indicated that Ca²⁺ adsorption decreases mobility within the carbonyl region of model bilayers more efficiently than monovalent ions do (Ca²⁺ ≫ Li⁺ > Na⁺ > K⁺ > Cs⁺). The effects of calcium ions were to a great extent mitigated in the plasma membranes isolated from HEK293 cells when compared to the model lipid membranes. Noticeably, the plasma membranes showed remarkably higher resistance toward rigidification induced by calcium ions even when compared with the model membranes containing cholesterol.

Preparation of asymmetric phospholipid vesicles for use as cell membrane models

Freely suspended liposomes are widely used as model membranes for studying lipid-lipid and protein-lipid interactions. Liposomes prepared by conventional methods have chemically identical bilayer leaflets. By contrast, living cells actively maintain different lipid compositions in the two leaflets of the plasma membrane, resulting in asymmetric membrane properties that are critical for normal cell function. Here, we present a protocol for the preparation of unilamellar asymmetric phospholipid vesicles that better mimic biological membranes. Asymmetry is generated by methyl-β-cyclodextrin-catalyzed exchange of the outer leaflet lipids between vesicle pools of differing lipid composition. Lipid destined for the outer leaflet of the asymmetric vesicles is provided by heavy-donor multilamellar vesicles containing a dense sucrose core. Donor lipid is exchanged into extruded unilamellar acceptor vesicles that lack the sucrose core, facilitating the post-exchange separation of the donor and acceptor pools by centrifugation because of differences in vesicle size and density. We present two complementary assays allowing quantification of each leaflet's lipid composition: the overall lipid composition is determined by gas chromatography-mass spectrometry, whereas the lipid distribution between the two leaflets is determined by NMR, using the lanthanide shift reagent Pr³⁺. The preparation protocol and the chromatographic assay can be applied to any type of phospholipid bilayer, whereas the NMR assay is specific to lipids with choline-containing headgroups, such as phosphatidylcholine and sphingomyelin. In ~12 h, the protocol can produce a large yield of asymmetric vesicles (up to 20 mg) suitable for a wide range of biophysical studies.

Mesoscale Computational Modeling of Protein-Membrane Interactions Based on Continuum Mean-Field Theory

Quantitative computational modeling of protein-membrane interactions is of great importance as it aids in the interpretation of experimental results and enables design and exploration of new experimental systems. This review describes one such computational approach conceived specifically to treat electrostatically driven interactions between a lipid membrane and a protein (or protein domains) adsorbing onto the membrane. The methodology is based on self-consistent minimization of the governing free energy functional which is expressed in the mean-field approximation and has contributions from electrostatic interactions as well as from mixing entropy of lipids in the membrane and ions in the solution. The method enables calculation of the free energy of the binding process and quantification of the steady-state lipid distribution around the adsorbing protein. The extension of the method to include membrane deformation degrees of freedom further allows calculation of the equilibrium bilayer shape upon the protein binding.

The matrix domain of the Gag protein from avian sarcoma virus contains a PI(4,5)P2-binding site that targets Gag to the cell periphery

The Gag protein of avian sarcoma virus (ASV) lacks an N-myristoyl (myr) group, but contains structural domains similar to those of HIV-1 Gag. Similarly to HIV-1, ASV Gag accumulates on the plasma membrane (PM) before egress; however, it is unclear whether the phospholipid PI(4,5)P₂ binds directly to the matrix (MA) domain of ASV Gag, as is the case for HIV-1 Gag. Moreover, the role of PI(4,5)P₂ in ASV Gag localization and budding has been controversial. Here, we report that substitution of residues that define the PI(4,5)P₂-binding site in the ASV MA domain (reported in an accompanying paper) interfere with Gag localization to the cell periphery and inhibit the production of virus-like particles (VLPs). We show that co-expression of Sprouty2 (Spry2) or the pleckstrin homology domain of phospholipase Cδ (PH-PLC), two proteins that bind PI(4,5)P₂, affects ASV Gag trafficking to the PM and budding. Replacement of the N-terminal 32 residues of HIV-1 MA, which encode its N-terminal myr signal and its PI(4,5)P₂-binding site, with the structurally equivalent N-terminal 24 residues of ASV MA created a chimera that localized at the PM and produced VLPs. In contrast, the homologous PI(4,5)P₂-binding signal in ASV MA could target HIV-1 Gag to the PM when substituted, but did not support budding. Collectively, these findings reveal a basic patch in both ASV and HIV-1 Gag capable of mediating PM binding and budding for ASV but not for HIV-1 Gag. We conclude that PI(4,5)P₂ is a strong determinant of ASV Gag targeting to the PM and budding.

Structural basis for targeting avian sarcoma virus Gag polyprotein to the plasma membrane for virus assembly

For most retroviruses, including HIV-1, binding of the Gag polyprotein to the plasma membrane (PM) is mediated by interactions between Gag's N-terminal myristoylated matrix (MA) domain and phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂) in the PM. The Gag protein of avian sarcoma virus (ASV) lacks the N-myristoylation signal but contains structural domains having functions similar to those of HIV-1 Gag. The molecular mechanism by which ASV Gag binds to the PM is incompletely understood. Here, we employed NMR techniques to elucidate the molecular determinants of the membrane-binding domain of ASV MA (MA₈₇) to lipids and liposomes. We report that MA₈₇ binds to the polar head of phosphoinositides such as PI(4,5)P₂ We found that MA₈₇ binding to inositol phosphates (IPs) is significantly enhanced by increasing the number of phosphate groups, indicating that the MA₈₇-IP binding is governed by charge-charge interactions. Using a sensitive NMR-based liposome-binding assay, we show that binding of MA₈₇ to liposomes is enhanced by incorporation of PI(4,5)P₂ and phosphatidylserine. We also show that membrane binding is mediated by a basic surface formed by Lys-6, Lys-13, Lys-23, and Lys-24. Substitution of these residues to glutamate abolished binding of MA₈₇ to both IPs and liposomes. In an accompanying paper, we further report that mutation of these lysine residues diminishes Gag assembly on the PM and inhibits ASV particle release. These findings provide a molecular basis for ASV Gag binding to the inner leaflet of the PM and advance our understanding of the basic mechanisms of retroviral assembly.

Interaction of the β amyloid - Aβ(25-35) - peptide with zwitterionic and negatively charged vesicles with and without cholesterol

The interactions of the Alzheimer's β-amyloid peptide, Aβ(25-35), with 18:1 (Δ9-Cis) PC 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), L-α-phosphatidylcholine (EPC), 1,2-dioleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (sodium salt) (DOPG), and L-α-phosphatidylglycerol (EPG) phospholipid vesicles with and without cholesterol (Ch) are studied by the nitroxide spin probe electron paramagnetic resonance (EPR) method. Two nitroxide spin probes, 2,2,6,6-tetramethyl-piperidin-1-oxyl-4-yl hexadecanoate (TP, TEMPO-Palmitate) and 2-Ethyl-2-(15-methoxy-15-oxopentadecyl)-4,4-dimethyl-3-oxazolidinyloxy (16-DSE), are utilized in the study. TEMPO-Palmitate has the reporting EPR moiety located at the top of this spin probe, while 16-DSE has the reporting EPR moiety located at the tail of the spin probe. These two probes enable us to sample the surface and the middle of the phospholipid bilayer, respectively. All EPR measurements are done above the melting points of all four phospholipids when the bilayer is in the liquid crystal phase, the physiologically relevant phase. Due to non-linear spectral line fitting, the EPR spectral parameters are extracted with high precision. The results show that there are two populations of Aβ(25-35) and that one of them is located in the hydrophobic phospholipid layer below the hydrophilic headgroup region. The second population appears to be weakly coupled to the surface of the bilayer. Both hydrophobic and electrostatic interactions affect the insertion of Aβ(25-35) in the bilayer. Also, there is strong evidence for an interaction between cholesterol and Aβ(25-35), which affects the dielectric and dynamic properties of the bilayer.