Structural properties of docosahexaenoyl phospholipid bilayers investigated by solid-state 2H NMR spectroscopy

Glycerophospholipid polyunsaturation modulates resveratrol action on biomimetic membranes

The phytoalexin resveratrol has received increasing attention for its potential to prevent oxidative damages in human organism. To shed further light on molecular mechanisms of its interaction with lipid membranes we study resveratrol influence on the organisation and mechanical properties of biomimetic lipid systems composed of synthetic phosphatidylcholines with mixed aliphatic chains and different degree of unsaturation at sn-2 position (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine, POPC, and 1-palmitoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine, PDPC). High-sensitivity isothermal titration calorimetric measurements reveal stronger spontaneous resveratrol association to polyunsaturated phosphatidylcholine bilayers compared to the monounsaturated ones resulting from hydrophobic interactions, conformational changes of the interacting species and desolvation of molecular surfaces. The latter is supported by the results from Laurdan spectroscopy of large unilamellar vesicles providing data on hydration at the glycerol backbones of glycerophospholipides. Higher degree of lipid order is reported for POPC membranes compared to PDPC. While resveratrol mostly enhances the hydration of PDPC membranes, increasing POPC dehydration is reported upon treatment with the polyphenol. Dehydration of the polyunsaturated lipid bilayers is measured only at the highest phytoalexin content studied (resveratrol/lipid 0.5 mol/mol) and is less pronounced than the effect reported for POPC membranes. The polyphenol effect on membrane mechanics is probed by thermal shape fluctuation analysis of quasispherical giant unilamellar vesicles. Markedly different trend of the bending elasticity with increasing resveratrol concentration is reported for the two types of phospholipid bilayers studied. POPC membranes become more rigid in the presence of resveratrol, whereas PDPC-containing bilayers exhibit softening at lower concentrations of the polyphenol followed by a slight growth without bilayer stiffening even at the highest resveratrol content explored. The new data on the structural organization and membrane properties of resveratrol-treated phosphatidylcholine membranes may underpin the development of future liposomal applications of the polyphenol in medicinal chemistry.

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

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

Cholesterol Stiffening of Lipid Membranes

Biomembrane order, dynamics, and other essential physicochemical parameters are controlled by cholesterol, a major component of mammalian cell membranes. Although cholesterol is well known to exhibit a condensing effect on fluid lipid membranes, the extent of stiffening that occurs with different degrees of lipid acyl chain unsaturation remains an enigma. In this review, we show that cholesterol locally increases the bending rigidity of both unsaturated and saturated lipid membranes, suggesting there may be a length-scale dependence of the bending modulus. We review our published data that address the origin of the mechanical effects of cholesterol on unsaturated and polyunsaturated lipid membranes and their role in biomembrane functions. Through a combination of solid-state deuterium NMR spectroscopy and neutron spin-echo spectroscopy, we show that changes in molecular packing cause the universal effects of cholesterol on the membrane bending rigidity. Our findings have broad implications for the role of cholesterol in lipid–protein interactions as well as raft-like mixtures, drug delivery applications, and the effects of antimicrobial peptides on lipid membranes.

The structures of polyunsaturated lipid bilayers by joint refinement of neutron and X-ray scattering data

We present the detailed structural analysis of polyunsaturated fatty acid-containing phospholipids namely, 1-palmitoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (PDPC) and 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (SDPC). A newly developed molecular dynamics (MD) simulation parsing scheme for lipids containing fatty acids with multiple double bonds was implemented into the scattering density profile (SDP) model to simultaneously refine differently contrasted neutron and X-ray scattering data. SDP analyses of scattering data at 30 °C yielded lipid areas of 71.1 Ų and 70.4 Ų for PDPC and SDPC bilayers, respectively, and a model free analysis of PDPC at 30 °C resulted in a lipid area of 72 Ų. In addition to bilayer structural parameters, using area-constrained MD simulations we determined the area compressibility modulus, K_A, to be 246.4 mN/m, a value similar to other neutral phospholipids.

Flexible lipid nanomaterials studied by NMR spectroscopy

Advances in solid-state nuclear magnetic resonance spectroscopy inform the emergence of material properties from atomistic-level interactions in membrane lipid nanostructures.

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

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

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

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

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

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

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

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

α-Tocopherol's Location in Membranes Is Not Affected by Their Composition

To this day, α-tocopherol's (aToc) role in humans is not well known. In previous studies, we have tried to connect aToc's biological function with its location in a lipid bilayer. In the present study, we have determined, by means of small-angle neutron diffraction, that not only is aToc's hydroxyl group located high in the membrane but its tail also resides far from the center of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. In addition, we located aToc's hydroxyl group above the lipid backbone in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine (POPS), and sphingomyelin bilayers, suggesting that aToc's location near the lipid-water interface may be a universal property of vitamin E. In light of these data, how aToc efficiently terminates lipid hydroperoxy radicals at the membrane center remains an open question.

The glycan role in the glycopeptide immunogenicity revealed by atomistic simulations and spectroscopic experiments on the multiple sclerosis biomarker CSF114(Glc)

Glycoproteins are often recognized as not-self molecules by antibodies triggering the onset of severe autoimmune diseases such as Multiple Sclerosis (MS). Thus, the development of antigen-mimicking biomarkers represents an attractive strategy for an early diagnosis of the disease. An example is the synthetic glycopeptide CSF114(Glc), which was designed and tested as MS biomarker and whose clinical application was limited by its reduced ability to detect autoantibodies in MS patients. In the attempt to improve the efficacy of CSF114(Glc), we have characterized all the events leading to the final binding of the biomarker to the autoantibody using atomistic simulations, ESR and NMR experiments. The glycosydic moiety plays a primary role in the whole process. In particular, in an environment mimicking that used in the clinical tests the glycopeptide assumes a α-helix structure that is functional for the interaction with the antibody. In this conformation CSF114(Glc) binds the monoclonal antibody mAb8-18C5 similarly to the myelin oligodendrocyte glycoprotein MOG, which is a known MS auto-antigen, thus explaining its diagnostic activity. Our study offers new molecular bases to design more effective biomarkers and provides a most valid protocol to investigate other systems where the environment effect is determinant for the biological activity.

Effects of lipid interactions on model vesicle engulfment by alveolar macrophages

The engulfment function of macrophages relies on complex molecular interactions involving both lipids and proteins. In particular, the clearance of apoptotic bodies (efferocytosis) is enabled by externalization on the cell target of phosphatidylserine lipids, which activate receptors on macrophages, suggesting that (local) specific lipid-protein interactions are required at least for the initiation of efferocytosis. However, in addition to apoptotic cells, macrophages can engulf foreign bodies that vary substantially in size from a few nanometers to microns, suggesting that nonspecific interactions over a wide range of length scales could be relevant. Here, we use model lipid membranes (made of phosphatidylcholine, phosphatidylserine, and ceramide) and rat alveolar macrophages to show how lipid bilayer properties probed by small-angle x-ray scattering and solid-state (2)H NMR correlate with engulfment rates measured by flow cytometry. We find that engulfment of protein-free model lipid vesicles is promoted by the presence of phosphatidylserine lipids but inhibited by ceramide, in accord with a previous study of apoptotic cells. We conclude that the roles of phosphatidylserine and ceramide in phagocytosis is based, at least in part, on lipid-mediated modification of membrane physical properties, including interactions at large length scales as well as local lipid ordering and possible domain formation.

Elastic deformation and area per lipid of membranes: atomistic view from solid-state deuterium NMR spectroscopy

This article reviews the application of solid-state ²H nuclear magnetic resonance (NMR) spectroscopy for investigating the deformation of lipid bilayers at the atomistic level. For liquid-crystalline membranes, the average structure is manifested by the segmental order parameters (SCD) of the lipids. Solid-state ²H NMR yields observables directly related to the stress field of the lipid bilayer. The extent to which lipid bilayers are deformed by osmotic pressure is integral to how lipid-protein interactions affect membrane functions. Calculations of the average area per lipid and related structural properties are pertinent to bilayer remodeling and molecular dynamics (MD) simulations of membranes. To establish structural quantities, such as area per lipid and volumetric bilayer thickness, a mean-torque analysis of ²H NMR order parameters is applied. Osmotic stress is introduced by adding polymer solutions or by gravimetric dehydration, which are thermodynamically equivalent. Solid-state NMR studies of lipids under osmotic stress probe membrane interactions involving collective bilayer undulations, order-director fluctuations, and lipid molecular protrusions. Removal of water yields a reduction of the mean area per lipid, with a corresponding increase in volumetric bilayer thickness, by up to 20% in the liquid-crystalline state. Hydrophobic mismatch can shift protein states involving mechanosensation, transport, and molecular recognition by G-protein-coupled receptors. Measurements of the order parameters versus osmotic pressure yield the elastic area compressibility modulus and the corresponding bilayer thickness at an atomistic level. Solid-state ²H NMR thus reveals how membrane deformation can affect protein conformational changes within the stress field of the lipid bilayer.

A Hybrid Approach for Highly Coarse-grained Lipid Bilayer Models

We present a systematic methodology to develop highly coarse-grained (CG) lipid models for large scale bio-membrane simulations, in which we derive CG interactions using a powerful combination of the multiscale coarse-graining (MS-CG) method, and an analytical form of the CG potential to model interactions at short range. The resulting hybrid coarse-graining (HCG) methodology is used to develop a three-site solvent-free model for 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), and a 1:1 mixture of 1,2-dioleoyl-sn-glycero-3-phospho-L-serine (DOPS) and DOPC. In addition, we developed a four-site model of DOPC, demonstrating the capability of the HCG methodology in designing model lipid systems of a desired resolution. We carried out microsecond-scale molecular dynamics (MD) simulations of large vesicles, highlighting the ability of the model to study systems at mesoscopic length and time scales. The models of DLPC, DOPC and DOPC-DOPS have elastic properties consistent with experiment and structural properties such as the radial distribution functions (RDF), bond and angle distributions, and the z-density distributions that compare well with reference all-atom systems.

Molecular dynamics simulations reveal specific interactions of post-translational palmitoyl modifications with rhodopsin in membranes

We present a detailed analysis of the behavior of the highly flexible post-translational lipid modifications of rhodopsin from multiple-microsecond all-atom molecular dynamics simulations. Rhodopsin was studied in a realistic membrane environment that includes cholesterol, as well as saturated and polyunsaturated lipids with phosphocholine and phosphoethanolamine headgroups. The simulation reveals striking differences between the palmitoylations at Cys322 and Cys323 as well as between the palmitoyl chains and the neighboring lipids. Notably the palmitoyl group at Cys322 shows considerably greater contact with helix H1 of rhodopsin, yielding frequent chain upturns with longer reorientational correlation times, and relatively low order parameters. While the palmitoylation at Cys323 makes fewer protein contacts and has increased order compared to Cys322, it nevertheless exhibits greater flexibility with smaller order parameters than the stearoyl chains of the surrounding lipids. The dynamical structure of the palmitoylations-as well as their extensive fluctuations-suggests a complex function for the post-translational modifications in rhodopsin and potentially other G protein-coupled receptors, going beyond their role as membrane anchoring elements. Rather, we propose that the palmitoylation at Cys323 has a potential role as a lipid anchor, whereas the palmitoyl-protein interaction observed for Cys322 suggests a more specific interaction that affects the stability of the dark state of rhodopsin.

Molecular simulations and solid-state NMR investigate dynamical structure in rhodopsin activation

Rhodopsin has served as the primary model for studying G protein-coupled receptors (GPCRs)-the largest group in the human genome, and consequently a primary target for pharmaceutical development. Understanding the functions and activation mechanisms of GPCRs has proven to be extraordinarily difficult, as they are part of a complex signaling cascade and reside within the cell membrane. Although X-ray crystallography has recently solved several GPCR structures that may resemble the activated conformation, the dynamics and mechanism of rhodopsin activation continue to remain elusive. Notably solid-state ((2))H NMR spectroscopy provides key information pertinent to how local dynamics of the retinal ligand change during rhodopsin activation. When combined with molecular mechanics simulations of proteolipid membranes, a new paradigm for the rhodopsin activation process emerges. Experiment and simulation both suggest that retinal isomerization initiates the rhodopsin photocascade to yield not a single activated structure, but rather an ensemble of activated conformational states. This article is part of a Special Issue entitled: Membrane protein structure and function.

Solid-state ²H NMR shows equivalence of dehydration and osmotic pressures in lipid membrane deformation

Lipid bilayers represent a fascinating class of biomaterials whose properties are altered by changes in pressure or temperature. Functions of cellular membranes can be affected by nonspecific lipid-protein interactions that depend on bilayer material properties. Here we address the changes in lipid bilayer structure induced by external pressure. Solid-state ²H NMR spectroscopy of phospholipid bilayers under osmotic stress allows structural fluctuations and deformation of membranes to be investigated. We highlight the results from NMR experiments utilizing pressure-based force techniques that control membrane structure and tension. Our ²H NMR results using both dehydration pressure (low water activity) and osmotic pressure (poly(ethylene glycol) as osmolyte) show that the segmental order parameters (S(CD)) of DMPC approach very large values of ≈ 0.35 in the liquid-crystalline state. The two stresses are thermodynamically equivalent, because the change in chemical potential when transferring water from the interlamellar space to the bulk water phase corresponds to the induced pressure. This theoretical equivalence is experimentally revealed by considering the solid-state ²H NMR spectrometer as a virtual osmometer. Moreover, we extend this approach to include the correspondence between osmotic pressure and hydrostatic pressure. Our results establish the magnitude of the pressures that lead to significant bilayer deformation including changes in area per lipid and volumetric bilayer thickness. We find that appreciable bilayer structural changes occur with osmotic pressures in the range of 10-100 atm or lower. This research demonstrates the applicability of solid-state ²H NMR spectroscopy together with bilayer stress techniques for investigating the mechanism of pressure sensitivity of membrane proteins.

An NMR database for simulations of membrane dynamics

Computational methods are powerful in capturing the results of experimental studies in terms of force fields that both explain and predict biological structures. Validation of molecular simulations requires comparison with experimental data to test and confirm computational predictions. Here we report a comprehensive database of NMR results for membrane phospholipids with interpretations intended to be accessible by non-NMR specialists. Experimental ¹³C-¹H and ²H NMR segmental order parameters (S(CH) or S(CD)) and spin-lattice (Zeeman) relaxation times (T(1Z)) are summarized in convenient tabular form for various saturated, unsaturated, and biological membrane phospholipids. Segmental order parameters give direct information about bilayer structural properties, including the area per lipid and volumetric hydrocarbon thickness. In addition, relaxation rates provide complementary information about molecular dynamics. Particular attention is paid to the magnetic field dependence (frequency dispersion) of the NMR relaxation rates in terms of various simplified power laws. Model-free reduction of the T(1Z) studies in terms of a power-law formalism shows that the relaxation rates for saturated phosphatidylcholines follow a single frequency-dispersive trend within the MHz regime. We show how analytical models can guide the continued development of atomistic and coarse-grained force fields. Our interpretation suggests that lipid diffusion and collective order fluctuations are implicitly governed by the viscoelastic nature of the liquid-crystalline ensemble. Collective bilayer excitations are emergent over mesoscopic length scales that fall between the molecular and bilayer dimensions, and are important for lipid organization and lipid-protein interactions. Future conceptual advances and theoretical reductions will foster understanding of biomembrane structural dynamics through a synergy of NMR measurements and molecular simulations.

Liposome-incorporated DHA increases neuronal survival by enhancing non-amyloidogenic APP processing

The fluidity of neuronal membranes plays a pivotal role in brain aging and neurodegeneration. In this study, we investigated the role of the omega-3 fatty acid docosahexaenoic acid (DHA) in modulation of membrane fluidity, APP processing, and protection from cytotoxic stress. To this end, we applied unilamellar transfer liposomes, which provided protection from oxidation and effective incorporation of DHA into cell membranes. Liposomes transferring docosanoic acid (DA), the completely saturated form of DHA, to the cell cultures served as controls. In HEK-APP cells, DHA significantly increased membrane fluidity and non-amyloidogenic processing of APP, leading to enhanced secretion of sAPPα. This enhanced secretion of sAPPα was associated with substantial protection against apoptosis induced by ER Ca(2+) store depletion. sAPPα-containing supernatants obtained from HEK-APP cells exerted similar protective effects as DHA in neuronal PC12 cells and HEK293 control cells. Correlating to further increased sAPPα levels, supernatants obtained from DHA-treated HEK-APP cells enhanced protection, whereas supernatants obtained from DHA-treated HEK293 control cells did not inhibit apoptosis, likely due to the low expression of endogenous APP and negligible sAPPα secretion in these cells. Further experiments with the small molecule inhibitors LY294002 and SP600125 indicated that sAPPα-induced cytoprotection relied on activation of the anti-apoptotic PI3K/Akt pathway and inhibition of the stress-triggered JNK signaling pathway in PC12 cells. Our data suggest that liposomal DHA is able to restore or maintain physiological membrane properties, which are required for neuroprotective sAPPα secretion and autocrine modulation of neuronal survival.

A systematically coarse-grained solvent-free model for quantitative phospholipid bilayer simulations

We present an implicit solvent coarse-grained (CG) model for quantitative simulations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. The absence of explicit solvent enables membrane simulations on large length and time scales at moderate computational expense. Despite improved computational efficiency, the model preserves chemical specificity and quantitative accuracy. The bonded and nonbonded interactions together with the effective cohesion mimicking the hydrophobic effect were systematically tuned by matching structural and mechanical properties from experiments and all-atom bilayer simulations, such as saturated area per lipid, radial distribution functions, density and pressure profiles across the bilayer, P(2) order, etc. The CG lipid model is shown to self-assemble into a bilayer starting from a random dispersion. Its line tension and elastic properties, such as bending and stretching modulus, are semiquantitatively consistent with experiments. The effects of (i) reduced molecular friction and (ii) more efficient integration combine to an overall speed-up of 3-4 orders of magnitude compared to all-atom bilayer simulations. Our CG lipid model is especially useful for studies of large-scale phenomena in membranes that nevertheless require a fair description of chemical specificity, e.g., membrane patches interacting with movable and transformable membrane proteins and peptides.

Systematic implicit solvent coarse-graining of bilayer membranes: lipid and phase transferability of the force field

We study the lipid and phase transferability of our recently developed systematically coarse-grained solvent-free membrane model. The force field was explicitly parameterized to describe a fluid 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayer at 310 K with correct structure and area per lipid, while gaining at least three orders of magnitude in computational efficiency (see Wang and Deserno 2010 J. Phys. Chem. B 114 11207-20). Here, we show that exchanging CG tails, without any subsequent re-parameterization, creates reliable models of 1,2-dioleoylphosphatidylcholine (DOPC) and 1,2-dipalmitoylphosphatidylcholine (DPPC) lipids in terms of structure and area per lipid. Furthermore, all CG lipids undergo a liquid-gel transition upon cooling, with characteristics like those observed in experiments and all-atom simulations during phase transformation. These studies suggest a promising transferability of our force field parameters to different lipid species and thermodynamic state points, properties that are a prerequisite for even more complex systems, such as mixtures.

Insights from biophysical studies on the role of polyunsaturated fatty acids for function of G-protein coupled membrane receptors

The composition of the lipid matrix is critical for function of membrane proteins. Perhaps one of the best studied examples is the function of the G-protein-coupled membrane receptor (GPCR) rhodopsin which is located in membranes with high content of phospholipids with polyunsaturated docosahexaenoic acid chains (DHA, 22:6n-3). Technological advances enabled a more detailed study of structure and dynamics of DHA chains and their interaction with rhodopsin. It was established that polyunsaturated DHA differs from saturated and monounsaturated hydrocarbon chains by far more rapid structural conversions. Furthermore, DHA chains tend to have higher density near the lipid/water interface while density of saturated chains is higher in the bilayer center. The interface of rhodopsin has a small number of sites for tighter interaction with DHA. Polyunsaturated phosphatidylethanolamines accumulate preferentially near the protein. Surprisingly, the high conformational freedom of most DHA chains is not measurably reduced upon interaction with rhodopsin. While some observations point at an involvement of continuum elastic properties of membranes in modulation of rhodopsin function, there is growing evidence for a role of weakly specific DHA-rhodopsin interactions.

Flexibility of ras lipid modifications studied by 2H solid-state NMR and molecular dynamics simulations

Human posttranslationally modified N-ras oncogenes are known to be implicated in numerous human cancers. Here, we applied a combination of experimental and computational techniques to determine structural and dynamical details of the lipid chain modifications of an N-ras heptapeptide in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) membranes. Experimentally, 2H NMR spectroscopy was used to study oriented membranes that incorporated ras heptapeptides with two covalently attached perdeuterated hexadecyl chains. Atomistic molecular dynamics simulations of the same system were carried out over 100 ns including 60 DMPC and 4 ras molecules. Several structural and dynamical experimental parameters could be directly compared to the simulation. Experimental and simulated 2H NMR order parameters for the methylene groups of the ras lipid chains exhibited a systematic difference attributable to the absence of collective motions in the simulation and to geometrical effects. In contrast, experimental 2H NMR spin-lattice relaxation rates for Zeeman order were well reproduced in the simulation. The lack of slower collective motions in the simulation did not appreciably influence the relaxation rates at a Larmor frequency of 115.1 MHz. The experimental angular dependence of the 2H NMR relaxation rates with respect to the external magnetic field was also relatively well simulated. These relaxation rates showed a weak angular dependence, suggesting that the lipid modifications of ras are very flexible and highly mobile in agreement with the low order parameters. To quantify these results, the angular dependence of the 2H relaxation rates was calculated by an analytical model considering both molecular and collective motions. Peptide dynamics in the membrane could be modeled by an anisotropic diffusion tensor with principal values of Dparallel=2.1x10(9) s(-1) and Dperpendicular=4.5x10(5) s(-1). A viscoelastic fitting parameter describing the membrane elasticity, viscosity, and temperature was found to be relatively similar for the ras peptide and the DMPC host matrix. Large motional amplitudes and relatively short correlation times facilitate mixing and dispersal with the lipid bilayer matrix, with implications for the role of the full-length ras protein in signal transduction and oncogenesis.

Polyunsaturation in lipid membranes: dynamic properties and lateral pressure profiles

We elucidate the influence of unsaturation on single-component membrane properties, focusing on their dynamical aspects and lateral pressure profiles across the membrane. To this end, we employ atomistic molecular dynamics simulations to study five different membrane systems with varying degrees of unsaturation, starting from saturated membranes and systematically increasing the level of unsaturation, ending up with a bilayer of phospholipids containing the docosahexaenoic acid. For an increasing level of unsaturation, we find considerable effects on dynamical properties, such as accelerated dynamics of the phosphocholine head groups and glycerol backbones and speeded up rotational dynamics of the lipid molecules. The lateral pressure profile is found to be altered by the degree of unsaturation. For an increasing number of double bonds, the peak in the middle of the bilayer decreases. This is compensated for by changes in the membrane-water interface region in terms of increasing peak heights of the lateral pressure profile. Implications of the findings are briefly discussed.

Effect of chain unsaturation on bilayer response to pressure

Using wide-line deuterium NMR, the effects of pressure on saturated-chain orientational order and gel-to-liquid-crystal phase transition temperature were observed in bilayers of 16:0-18:1 PC-d31 (POPC-d31) and 16:0-18:2 PC-d31 (PLPC-d31). Spectra were recorded for a range of pressures at selected temperatures and for a range of temperatures at selected pressures up to 193 MPa. For 16:0-18:1 PC-d31, the main transition temperature increased by approximately 0.18 K/MPa, a rate that is similar to what is found for bilayers of disaturated PC's. For 16:0-18:2 PC-d31, the increase in transition temperature with pressure was slightly smaller at approximately 0.13 K/MPa. To investigate the isothermal response of chain orientational order to pressure, spectra for each lipid were obtained for three pressures (ambient, 55 MPa, and 110 MPa) near room temperature (approximately 25 degrees C) and for three pressures (ambient, 110 MPa, and 193 MPa) at higher temperature (approximately 40 degrees C). These temperatures were chosen such that the difference between the higher observation temperature and the main transition of 16:0-18:1 PC-d31 would be similar to the difference between the lower observation temperature and the main transition of 16:0-18:2 PC-d31. Application of a given pressure was found to raise the orientational order for all methylene groups on the saturated chain of a particular lipid by roughly similar amounts. For comparable pressure differences, the pressure-induced ordering of the 16:0-18:1 PC-d31 saturated chain at approximately 40 degrees C was greater than that of the corresponding chain in 16:0-18:2 PC- d31 at approximately 25 degrees C. These observations suggest that increasing levels of chain unsaturation may reduce the sensitivity of bilayer order to variations in pressure at corresponding temperatures relative to their ambient pressure transitions.

The structure of polyunsaturated lipid bilayers important for rhodopsin function: a neutron diffraction study

The structure of oriented 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine bilayers with perdeuterated stearoyl- or docosahexaenoyl hydrocarbon chains was investigated by neutron diffraction. Experiments were conducted at two different relative humidities, 66 and 86%. At both humidities we observed that the polyunsaturated docosahexaenoyl chain has a preference to reside near the lipid water interface. That leaves voids in the bilayer center that are occupied by saturated stearoyl chain segments. This uneven distribution of saturated- and polyunsaturated chain densities is likely to result in membrane elastic stress that modulates function of integral receptor proteins like rhodopsin.

Lipid Modifications of a Ras Peptide Exhibit Altered Packing and Mobility versus Host Membrane as Detected by 2H Solid-State NMR

The human N-ras protein binds to cellular membranes by insertion of two covalently bound posttranslational lipid modifications, which is crucial for its function in signal transduction and cell proliferation. Mutations in ras may lead to unregulated cell growth and eventually cancer, making it an important therapeutic target. Here we have investigated the molecular details of the membrane binding mechanism. A heptapeptide derived from the C-terminus of the human N-ras protein was synthesized including two hexadecyl modifications. Solid-state 2H NMR was used to determine the packing and molecular dynamics of the ras lipid chains as well as the phospholipid matrix. Separately labeling the chains of the peptide and the phospholipids with 2H enabled us to obtain atomically resolved parameters relevant to their structural dynamics. While the presence of ras only marginally affected the packing of DMPC membranes, dramatically lower order parameters (S(CD)) were observed for the ras acyl chains indicating modified packing properties. Essentially identical projected lengths of the 16:0 ras chains and the 14:0 DMPC chains were found, implying that the polypeptide backbone is located at the lipid-water interface. Dynamical properties of both the ras and phospholipid chains were determined from spin-lattice 2H relaxation (R1Z) measurements. Plots of R1Z rates versus the corresponding squared segmental order parameters revealed striking differences. We propose the ras peptide is confined to microdomains containing DMPC chains which are in exchange with the bulk bilayer on the 2H NMR time scale (approximately 10(-5) s). Compared to the host DMPC matrix, the ras lipid modifications are extremely flexible and undergo relatively large amplitude motions. It is hypothesized that this flexibility is a requirement for the optimal anchoring of lipid-modified proteins to cellular membranes.

Specific volumes of unsaturated phosphatidylcholines in the liquid crystalline lamellar phase

The specific volumes of seven 1,2-diacyl-sn-glycero-3-phosphocholines with symmetric, unbranched acyl chains containing one, four, or six cis double bonds per chain, or with a saturated sn-1 chain and one, four, or six cis double bonds in the sn-2 chain were determined by the neutral buoyancy method. Experiments were conducted in the liquid crystalline lamellar phase over the temperature range from 5 to 35 degrees C. It is demonstrated that the molecular volume of phosphatidylcholines can be well approximated as the sum of a constant volume of the polar lipid head region and the temperature-dependent volumes of hydrocarbon chain CH2, CH, and terminal CH3 groups. A linear dependence of chain segment volumes on temperature was observed. A self-consistent set of partially temperature-dependent volumes is obtained that allows prediction of phosphatidylcholine molecular volumes within very tight error margins.

Oleic and docosahexaenoic acid differentially phase separate from lipid raft molecules: a comparative NMR, DSC, AFM, and detergent extraction study

We have previously suggested that the omega-3 polyunsaturated fatty acid, docosahexaenoic acid (DHA) may in part function by enhancing membrane lipid phase separation into lipid rafts. Here we further tested for differences in the molecular interactions of an oleic (OA) versus DHA-containing phospholipid with sphingomyelin (SM) and cholesterol (CHOL) utilizing (2)H NMR spectroscopy, differential scanning calorimetry, atomic force microscopy, and detergent extractions in model bilayer membranes. (2)H NMR and DSC (differential scanning calorimetry) established the phase behavior of the OA-containing 1-[(2)H(31)]palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (16:0-18:1PE-d(31))/SM (1:1) and the DHA-containing 1-[(2)H(31)]palmitoyl-2-docosahexaenoyl-sn-glycero-3-phosphoethanolamine (16:0-22:6PE-d(31))/SM (1:1) in the absence and presence of equimolar CHOL. CHOL was observed to affect the OA-containing phosphatidylethanolamine (PE) more than the DHA-containing PE, as exemplified by >2 x greater increase in order measured for the perdeuterated palmitic chain in 16:0-18:1PE-d(31)/SM (1:1) compared to 16:0-22:6PE-d(31)/SM (1:1) bilayers in the liquid crystalline phase. Atomic force microscopy (AFM) experiments showed less lateral phase separation between 16:0-18:1PE-rich and SM/CHOL-rich raft domains in 16:0-18:1PE/SM/CHOL (1:1:1) bilayers than was observed when 16:0-22:6PE replaced 16:0-18:1PE. Differences in the molecular interaction of 16:0-18:1PE and 16:0-22:6PE with SM/CHOL were also found using biochemical detergent extractions. In the presence of equimolar SM/CHOL, 16:0-18:1PE showed decreased solubilization in comparison to 16:0-22:6PE, indicating greater phase separation with the DHA-PE. Detergent experiments were also conducted with cardiomyocytes fed radiolabeled OA or DHA. Although both OA and DHA were found to be largely detergent solubilized, the amount of OA that was found to be associated with raft-rich detergent-resistant membranes exceeded DHA by almost a factor of 2. We conclude that the OA-PE phase separates from rafts far less than DHA-PE, which may have implications for cellular signaling.

Packing and Viscoelasticity of Polyunsaturated ω-3 and ω-6 Lipid Bilayers as Seen by2H NMR and X-ray Diffraction

Polyunsaturated phospholipids of the omega-3 and omega-6 classes play key roles in cellular functions, yet their mechanisms of biological action are still a matter of debate. Using deuterium ((2)H) NMR spectroscopy and small-angle X-ray diffraction, we show how membrane properties are modified by docosahexaenoic (DHA; 22:6) and arachidonic (AA; 20:4) acyl chains of the omega-3 and the omega-6 families, respectively. Structural and dynamical differences due to polyunsaturation are evident in both the ordered and disordered phases of mixed-chain (16:0)(22:6)PC and (16:0)(20:4)PC bilayers. Due to the lower chain melting temperature, the omega-6 AA bilayer is more disordered in the fluid (L(alpha)) state than the omega-3 DHA bilayer; it is thinner with a larger area per lipid. The thermal hysteresis observed for the DHA bilayer may represent the influences of angle-iron conformers in the gel state and back-bended, hairpinlike conformers in the fluid state, consistent with molecular dynamics studies. Interpretation of the (2)H NMR order profiles of (16:0-d(31))(22:6)PC and (16:0-d(31))(20:4)PC together with X-ray electron density profiles reveals an uneven distribution of mass; i.e., the sn-1 saturated chain is displaced toward the membrane center, whereas the sn-2 polyunsaturated chain is shifted toward the bilayer aqueous interface. Moreover, the (2)H NMR relaxation rates are increased by the presence of omega-6 AA chains compared to omega-3 DHA chains. When evaluated at the same amplitude of motion, relaxation parameters give a naturally calibrated scale for comparison of fluid lipid bilayers. Within this framework, polyunsaturated bilayers are relatively soft to bending and area fluctuations on the mesoscale approaching molecular dimensions. Significant differences are evident in the viscoelastic properties of the omega-3 and omega-6 bilayers, a possibly biologically relevant feature that distinguishes between the two phospholipid classes.

Membrane model for the G-protein-coupled receptor rhodopsin: hydrophobic interface and dynamical structure

Rhodopsin is the only member of the pharmacologically important superfamily of G-protein-coupled receptors with a known structure at atomic resolution. A molecular dynamics model of rhodopsin in a POPC phospholipid bilayer was simulated for 15 ns, revealing a conformation significantly different from the recent crystal structures. The structure of the bilayer compared with a protein-free POPC control indicated hydrophobic matching with the nonpolar interface of the receptor, in agreement with deuterium NMR experiments. A new generalized molecular surface method, based on a three-dimensional Voronoi cell construction for atoms with different radii, was developed to quantify cross-sectional area profiles for the protein, lipid acyl chains and headgroups, and water. Thus, it was possible to investigate the bilayer deformation due to curvature of the individual lipid monolayers. Moreover, the generalized molecular surface derived hydrophobic interface allowed benchmarking of the hydropathy sequence analysis, an important structural genomics tool. Five water molecules diffused into internal hydration sites during the simulation, yielding a total of 12 internal waters. The cytoplasmic loops and the C-terminal tail, containing the G-protein recognition and protein sorting sequences, exhibited a high mobility, in marked contrast to the extracellular and transmembrane domains. The proposed functional coupling of the highly conserved ERY motif to the lipid-water interface via the cytoplasmic loops provides insight into lipid effects on G-protein-coupled receptor activation in terms of a flexible surface model, involving the spontaneous monolayer curvature.

Polyunsaturated docosahexaenoic vs docosapentaenoic acid-differences in lipid matrix properties from the loss of one double bond

Insufficient supply to the developing brain of docosahexaenoic acid (22:6n3, DHA), or its omega-3 fatty acid precursors, results in replacement of DHA with docosapentaenoic acid (22:5n6, DPA), an omega-6 fatty acid that is lacking a double bond near the chain's methyl end. We investigated membranes of 1-stearoyl(d(35))-2-docosahexaenoyl-sn-glycero-3-phosphocholine and 1-stearoyl(d(35))-2-docosapentaenoyl-sn-glycero-3-phosphocholine by solid-state NMR, X-ray diffraction, and molecular dynamics simulations to determine if the loss of this double bond alters membrane physical properties. The low order parameters of polyunsaturated chains and the NMR relaxation data indicate that both DHA and DPA undergo rapid conformational transitions with correlation times of the order of nanoseconds at carbon atom C(2) and of picoseconds near the terminal methyl group. However, there are important differences between DHA- and DPA-containing lipids: the DHA chain with one additional double bond is more flexible at the methyl end and isomerizes with shorter correlation times. Furthermore, the stearic acid paired with the DHA in mixed-chain lipids has lower order, in particular in the middle of the chain near carbons C(10)(-)(12), indicating differences in the packing of hydrocarbon chains. Such differences are also reflected in the electron density profiles of the bilayers and in the simulation results. The DHA chain has a higher density near the lipid-water interface, whereas the density of the stearic acid chain is higher in the bilayer center. The loss of a single double bond from DHA to DPA results in a more even distribution of chain densities along the bilayer normal. We propose that the function of integral membrane proteins such as rhodopsin is sensitive to such a redistribution.

Lipid-protein interactions in biological membranes: a structural perspective

Lipid molecules bound to membrane proteins are resolved in some high-resolution structures of membrane proteins. An analysis of these structures provides a framework within which to analyse the nature of lipid-protein interactions within membranes. Membrane proteins are surrounded by a shell or annulus of lipid molecules, equivalent to the solvent layer surrounding a water-soluble protein. The lipid bilayer extends right up to the membrane protein, with a uniform thickness around the protein. The surface of a membrane protein contains many shallow grooves and protrusions to which the fatty acyl chains of the surrounding lipids conform to provide tight packing into the membrane. An individual lipid molecule will remain in the annular shell around a protein for only a short period of time. Binding to the annular shell shows relatively little structural specificity. As well as the annular lipid, there is evidence for other lipid molecules bound between the transmembrane alpha-helices of the protein; these lipids are referred to as non-annular lipids. The average thickness of the hydrophobic domain of a membrane protein is about 29 A, with a few proteins having significantly smaller or greater thicknesses than the average. Hydrophobic mismatch between a membrane protein and the surrounding lipid bilayer generally leads to only small changes in membrane thickness. Possible adaptations in the protein to minimise mismatch include tilting of the helices and rotation of side chains at the ends of the helices. Packing of transmembrane alpha-helices is dependent on the chain length of the surrounding phospholipids. The function of membrane proteins is dependent on the thickness of the surrounding lipid bilayer, sometimes on the presence of specific, usually anionic, phospholipids, and sometimes on the phase of the phospholipid.

Elastic deformation of membrane bilayers probed by deuterium NMR relaxation

In deuterium ((2)H) NMR spectroscopy of fluid lipid bilayers, the average structure is manifested in the segmental order parameters (S(CD)) of the flexible molecules. The corresponding spin-lattice relaxation rates (R(1Z) depend on both the amplitudes and the rates of the segmental fluctuations, and indicate the types of lipid motions. By combining (2)H NMR order parameter measurements with relaxation studies, we have obtained a more comprehensive picture of lipids in the liquid-crystalline (L(alpha)) state than formerly possible. Our data suggest that a lipid bilayer constitutes an ordered fluid, in which the phospholipids are grafted to the aqueous interface via their polar headgroups, whereas the fatty acyl chains are in effect liquid hydrocarbon. Studies of (2)H-labeled saturated lipids indicate their R(1Z) rates and S(CD) order parameters are correlated by a model-free, square-law functional dependence, signifying the presence of relatively slow bilayer fluctuations. A new composite membrane deformation model explains simultaneously the frequency (magnetic field) dependence and the angular anisotropy of the relaxation. The results imply the R(1Z) rates are due to a broad spectrum of 3-D collective bilayer excitations, together with effective axial rotations of the lipids. For the first time, NMR relaxation studies show that the viscoelastic properties of membrane lipids at megahertz frequencies are modulated by the lipid acyl length (bilayer thickness), polar headgroups (bilayer interfacial area), inclusion of a nonionic detergent (C(12)E(8)), and the presence of cholesterol, leading to a range of bilayer softness. Our findings imply the concept of elastic deformation is relevant on lengths approaching the bilayer thickness and less (the mesoscopic scale), and suggest that application of combined R(1Z) and S(CD) studies of phospholipids can be used as a simple membrane elastometer. Heuristic estimates of the bilayer bending rigidity kappa and the area elastic modulus K(a) enable comparison to other biophysical studies, involving macroscopic deformation of thin membrane lipid films. Finally, the bilayer softness may be correlated with the lipid diversity of biomembranes, for example, with regard to membrane curvature, repulsive interactions between bilayers, and lipid-protein interactions.

Electrostatic properties of membrane lipids coupled to metarhodopsin II formation in visual transduction

Changes in lipid composition have recently been shown to exert appreciable influences on the activities of membrane-bound proteins and peptides. We tested the hypothesis that the conformational states of rhodopsin linked to visual signal transduction are related to biophysical properties of the membrane lipid bilayer. For bovine rhodopsin, the meta I-meta II conformational transition was studied in egg phosphatidylcholine (PC) recombinants versus the native rod outer segment (ROS) membranes by means of flash photolysis. Formation of metarhodopsin II was observed by the change in absorbance at 478 nm after a single actinic flash was delivered to the sample. The meta I/meta II ratio was investigated as a function of both temperature and pH. The data clearly demonstrated thermodynamic reversibility of the transition for both the egg PC recombinants and the native ROS membranes. A significant shift of the apparent pK(a) for the acid-base equilibrium to lower values was evident in the egg PC recombinant, with little meta II produced under physiological conditions. Calculations of the membrane surface pH using a Poisson-Boltzmann model suggested the free energies of the meta I and meta II states were significantly affected by electrostatic properties of the bilayer lipids. In the ROS membranes, phosphatidylserine (PS) is needed for full formation of meta II, in combination with phosphatidylethanolamine (PE) and polyunsaturated docosahexaenoic acid (DHA; 22:6omega3) chains. We propose that the PS surface potential leads to an accumulation of hydronium ions, H(3)O(+), in the electrical double layer, which drive the reaction together with the large negative spontaneous curvature (H(0)) conferred by PE plus DHA chains. The elastic stress/strain of the bilayer arises from an interplay of the approximately zero H(0) from PS and the negative H(0) due to the PE headgroups and polyunsaturated chains. The lipid influences are further explained in terms of matching of the bilayer spontaneous curvature to the curvature at the lipid/rhodopsin interface, as formulated by the Helfrich bending energy. These new findings guide current ideas as to how bilayer properties govern the conformational energetics of integral membrane proteins. Moreover, they yield knowledge of how membrane lipid-protein interactions involving acidic phospholipids such as PS and neutral polyunsaturated DHA chains are implicated in key biological functions such as vision.