Chapter 8 Influence of Nonlamellar-Forming Lipids on Rhodopsin

Biophysics of Membrane Stiffening by Cholesterol and Phosphatidylinositol 4,5-bisphosphate (PIP2)

Cell membranes regulate a wide range of phenomena that are implicated in key cellular functions. Cholesterol, a critical component of eukaryotic cell membranes, is responsible for cellular organization, membrane elasticity, and other critical physicochemical parameters. Besides cholesterol, other lipid components such as phosphatidylinositol 4,5-bisphosphate (PIP2) are found in minor concentrations in cell membranes yet can also play a major regulatory role in various cell functions. In this chapter, we describe how solid-state deuterium nuclear magnetic resonance (²H NMR) spectroscopy together with neutron spin-echo (NSE) spectroscopy can inform synergetic changes to lipid molecular packing due to cholesterol and PIP2 that modulate the bending rigidity of lipid membranes. Fundamental structure-property relations of molecular self-assembly are illuminated and point toward a length and time-scale dependence of cell membrane mechanics, with significant implications for biological activity and membrane lipid-protein interactions.

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.

Altered Membrane Mechanics Provides a Receptor-Independent Pathway for Serotonin Action

Serotonin, an important signaling molecule in humans, has an unexpectedly high lipid membrane affinity. The significance of this finding has evoked considerable speculation. Here we show that membrane binding by serotonin can directly modulate membrane properties and cellular function, providing an activity pathway completely independent of serotonin receptors. Atomic force microscopy shows that serotonin makes artificial lipid bilayers softer, and induces nucleation of liquid disordered domains inside the raft‐like liquid‐ordered domains. Solid‐state NMR spectroscopy corroborates this data at the atomic level, revealing a homogeneous decrease in the order parameter of the lipid chains in the presence of serotonin. In the RN46A immortalized serotonergic neuronal cell line, extracellular serotonin enhances transferrin receptor endocytosis, even in the presence of broad‐spectrum serotonin receptor and transporter inhibitors. Similarly, it increases the membrane binding and internalization of oligomeric peptides. Our results uncover a mode of serotonin–membrane interaction that can potentiate key cellular processes in a receptor‐independent fashion.

Curvature Energetics Determined by Alchemical Simulation on Four Topologically Distinct Lipid Phases

The relative curvature energetics of two lipids are tested using thermodynamic integration (TI) on four topologically distinct lipid phases. Simulations use TI to switch between choline headgroup lipids (POPC; that prefers to be flat) and ethanolamine headgroup lipids (POPE; that prefer, for example, the inner monolayer of vesicles). The thermodynamical moving of the lipids between planar, inverse hexagonal (H_II), cubic (Q_II; Pn3m space group), and vesicle topologies reveals differences in material parameters that were previously challenging to access. The methodology allows for predictions of two important lipid material properties: the difference in POPC/POPE monolayer intrinsic curvature (ΔJ₀) and the difference in POPC/POPE monolayer Gaussian curvature modulus (Δκ̅_m), both of which are connected to the energetics of topological variation. Analysis of the TI data indicates that, consistent with previous experiment and simulation, the J₀ of POPE is more negative than POPC (ΔJ₀ = -0.018 ± 0.001 Å^-1). The theoretical framework extracts significant differences in κ̅_m of which POPE is less negative than POPC by 2.0 to 4.0 kcal/mol. The range of these values is determined by considering subsets of the simulations, and disagreement between these subsets suggests separate mechanical parameters at very high curvature. Finally, the fit of the TI data to the model indicates that the position of the pivotal plane of curvature is not constant across topologies at high curvature. Overall, the results offer insights into lipid material properties, the limits of a single HC model, and how to test them using simulation.

Membrane Curvature Revisited-the Archetype of Rhodopsin Studied by Time-Resolved Electronic Spectroscopy

G-protein-coupled receptors (GPCRs) comprise the largest and most pharmacologically targeted membrane protein family. Here, we used the visual receptor rhodopsin as an archetype for understanding membrane lipid influences on conformational changes involved in GPCR activation. Visual rhodopsin was recombined with lipids varying in their degree of acyl chain unsaturation and polar headgroup size using 1-palmitoyl-2-oleoyl-sn-glycero- and 1,2-dioleoyl-sn-glycerophospholipids with phosphocholine (PC) or phosphoethanolamine (PE) substituents. The receptor activation profile after light excitation was measured using time-resolved ultraviolet-visible spectroscopy. We discovered that more saturated POPC lipids back shifted the equilibrium to the inactive state, whereas the small-headgroup, highly unsaturated DOPE lipids favored the active state. Increasing unsaturation and decreasing headgroup size have similar effects that combine to yield control of rhodopsin activation, and necessitate factors beyond proteolipid solvation energy and bilayer surface electrostatics. Hence, we consider a balance of curvature free energy with hydrophobic matching and demonstrate how our data support a flexible surface model (FSM) for the coupling between proteins and lipids. The FSM is based on the Helfrich formulation of membrane bending energy as we previously first applied to lipid-protein interactions. Membrane elasticity and curvature strain are induced by lateral pressure imbalances between the constituent lipids and drive key physiological processes at the membrane level. Spontaneous negative monolayer curvature toward water is mediated by unsaturated, small-headgroup lipids and couples directly to GPCR activation upon light absorption by rhodopsin. For the first time to our knowledge, we demonstrate this modulation in both the equilibrium and pre-equilibrium evolving states using a time-resolved approach.

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.

Perturbations of Native Membrane Protein Structure in Alkyl Phosphocholine Detergents: A Critical Assessment of NMR and Biophysical Studies

Membrane proteins perform a host of vital cellular functions. Deciphering the molecular mechanisms whereby they fulfill these functions requires detailed biophysical and structural investigations. Detergents have proven pivotal to extract the protein from its native surroundings. Yet, they provide a milieu that departs significantly from that of the biological membrane, to the extent that the structure, the dynamics, and the interactions of membrane proteins in detergents may considerably vary, as compared to the native environment. Understanding the impact of detergents on membrane proteins is, therefore, crucial to assess the biological relevance of results obtained in detergents. Here, we review the strengths and weaknesses of alkyl phosphocholines (or foscholines), the most widely used detergent in solution-NMR studies of membrane proteins. While this class of detergents is often successful for membrane protein solubilization, a growing list of examples points to destabilizing and denaturing properties, in particular for α-helical membrane proteins. Our comprehensive analysis stresses the importance of stringent controls when working with this class of detergents and when analyzing the structure and dynamics of membrane proteins in alkyl phosphocholine detergents.

Differential Effect of Bilayer Thickness on Sticholysin Activity

In this study, we examined the influence of bilayer thickness on the activity of the actinoporin toxins sticholysin I and II (StnI and StnII) at 25 °C. Bilayer thickness was varied using dimonounsaturated phosphatidylcholine (PC) analogues (with 14:1, 16:1, 18:1, 20:1, and 22:1 acyl chains). In addition, N-14:0-sphingomyelin (SM) was always included because StnI and StnII are SM specific. Cholesterol was also incorporated as indicated. In cholesterol-free large unilamellar vesicles (LUVs) the PC:SM molar ratio was 4:1, and when cholesterol was included, the complete molar ratio was 4:1:0.5 (PC:SM:cholesterol, respectively). Stn toxins promote bilayer leakage through pores formed by oligomerized toxin monomers. Initial calcein leakage was moderately dependent on bilayer PC acyl chain length (and thus bilayer thickness), with higher rates observed with di-16:1 and di-18:1 PC bilayers. In the presence of cholesterol, the maximum rates of calcein leakage were observed in di-14:1 and di-16:1 PC bilayers. Using isothermal titration calorimetry to study the Stn-LUV interaction, we observed that the bilayer affinity constant (K_a) peaked with LUVs containing di-18:1 PC, and was lower in shorter and longer PC acyl chain bilayers. The presence of cholesterol increased the binding affinity approximately 30-fold at the optimal bilayer thickness (di-18:1-PC). We conclude that bilayer thickness affects both functional and conformational aspects of Stn membrane binding and pore formation. Moreover, the length of the actinoporins' N-terminal α-helix, which penetrates the membrane to form a functional pore, appears to be optimal for the membrane thickness represented by di-18:1 PC.

Gramicidin A Channel Formation Induces Local Lipid Redistribution II: A 3D Continuum Elastic Model

To change conformation, a protein must deform the surrounding bilayer. In this work, a three-dimensional continuum elastic model for gramicidin A in a lipid bilayer is shown to describe the sensitivity to thickness, curvature stress, and the mechanical properties of the lipid bilayer. A method is demonstrated to extract the gramicidin-lipid boundary condition from all-atom simulations that can be used in the three-dimensional continuum model. The boundary condition affects the deformation dramatically, potentially much more than typical variations in the material stiffness do as lipid composition is changed. Moreover, it directly controls the sensitivity to curvature stress. The curvature stress and hydrophobic surfaces of the all-atom and continuum models are found to be in excellent agreement. The continuum model is applied to estimate the enrichment of hydrophobically matched lipids near the channel in a mixture, and the results agree with single-channel experiments and extended molecular dynamics simulations from the companion article by Beaven et al. in this issue of Biophysical Journal.

Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are central to many fundamental cellular signaling pathways. They transduce signals from the outside to the inside of cells in physiological processes ranging from vision to immune response. It is extremely challenging to look at them individually using conventional experimental techniques. Recently, a pseudo atomistic molecular model has emerged as a valuable tool to access information on GPCRs, more specifically on their interactions with their environment in their native cell membrane and the consequences on their supramolecular organization. This approach uses the Martini coarse grain (CG) model to describe the receptors, lipids, and solvent in molecular dynamics (MD) simulations and in enough detail to allow conserving the chemical specificity of the different molecules. The elimination of unnecessary degrees of freedom has opened up large-scale simulations of the lipid-mediated supramolecular organization of GPCRs. Here, after introducing the Martini CGMD method, we review these studies carried out on various members of the GPCR family, including rhodopsin (visual receptor), opioid receptors, adrenergic receptors, adenosine receptors, dopamine receptor, and sphingosine 1-phosphate receptor. These studies have brought to light an interesting set of novel biophysical principles. The insights range from revealing localized and heterogeneous deformations of the membrane bilayer at the surface of the protein, specific interactions of lipid molecules with individual GPCRs, to the effect of the membrane matrix on global GPCR self-assembly. The review ends with an overview of the lessons learned from the use of the CGMD method, the biophysical-chemical findings on lipid-protein interplay.

Can Specific Protein-Lipid Interactions Stabilize an Active State of the Beta 2 Adrenergic Receptor?

G-protein-coupled receptors are eukaryotic membrane proteins with broad biological and pharmacological relevance. Like all membrane-embedded proteins, their location and orientation are influenced by lipids, which can also impact protein function via specific interactions. Extensive simulations totaling 0.25 ms reveal a process in which phospholipids from the membrane's cytosolic leaflet enter the empty G-protein binding site of an activated β2 adrenergic receptor and form salt-bridge interactions that inhibit ionic lock formation and prolong active-state residency. Simulations of the receptor embedded in an anionic membrane show increased lipid binding, providing a molecular mechanism for the experimental observation that anionic lipids can enhance receptor activity. Conservation of the arginine component of the ionic lock among Rhodopsin-like G-protein-coupled receptors suggests that intracellular lipid ingression between receptor helices H6 and H7 may be a general mechanism for active-state stabilization.

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.

General and specific lipid-protein interactions in Na,K-ATPase

The molecular activity of Na,K-ATPase and other P2 ATPases like Ca(2+)-ATPase is influenced by the lipid environment via both general (physical) and specific (chemical) interactions. Whereas the general effects of bilayer structure on membrane protein function are fairly well described and understood, the importance of the specific interactions has only been realized within the last decade due particularly to the growing field of membrane protein crystallization, which has shed new light on the molecular details of specific lipid-protein interactions. It is a remarkable observation that specific lipid-protein interactions seem to be evolutionarily conserved, and conformations of specifically bound lipids at the lipid-protein surface within the membrane are similar in crystal structures determined with different techniques and sources of the protein, despite the rather weak lipid-protein interaction energy. Studies of purified detergent-soluble recombinant αβ or αβFXYD Na,K-ATPase complexes reveal three separate functional effects of phospholipids and cholesterol with characteristic structural selectivity. The observations suggest that these three effects are exerted at separate binding sites for phophatidylserine/cholesterol (stabilizing), polyunsaturated phosphatidylethanolamine (stimulatory), and saturated PC or sphingomyelin/cholesterol (inhibitory), which may be located within three lipid-binding pockets identified in recent crystal structures of Na,K-ATPase. The findings point to a central role of direct and specific interactions of different phospholipids and cholesterol in determining both stability and molecular activity of Na,K-ATPase and possible implications for physiological regulation by membrane lipid composition. This article is part of a special issue titled "Lipid-Protein Interactions."

Investigation of rhodopsin dynamics in its signaling state by solid-state deuterium NMR spectroscopy

Site-directed deuterium NMR spectroscopy is a valuable tool to study the structural dynamics of biomolecules in cases where solution NMR is inapplicable. Solid-state (2)H NMR spectral studies of aligned membrane samples of rhodopsin with selectively labeled retinal provide information on structural changes of the chromophore in different protein states. Moreover (2)H NMR relaxation time measurements allow one to study the dynamics of the ligand during the transition from the inactive to the active state. Here we describe the methodological aspects of solid-state (2)H NMR spectroscopy for functional studies of rhodopsin, with an emphasis on the dynamics of the retinal cofactor. We provide complete protocols for the preparation of NMR samples of rhodopsin with 11-cis-retinal selectively deuterated at the methyl groups in aligned membranes. In addition we review optimized conditions for trapping the rhodopsin photointermediates; and we address the challenging problem of trapping the signaling state of rhodopsin in aligned membrane films.

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.

Coarse-grained molecular dynamics provides insight into the interactions of lipids and cholesterol with rhodopsin

Protein function is a complicated interplay between structure and dynamics, which can be heavily influenced by environmental factors and conditions. This is particularly true in the case of membrane proteins, such as the visual receptor rhodopsin. It has been well documented that lipid headgroups, polyunsaturated tails, and the concentration of cholesterol in membranes all play a role in the function of rhodopsin. Recently, we used all-atom simulations to demonstrate that different lipid species have preferential interactions and possible binding sites on rhodopsin's surface, consistent with experiment. However, the limited timescales of the simulations meant that the statistical uncertainty of these results was substantial. Accordingly, we present here 32 independent 1.6 μs coarse-grained simulations exploring lipids and cholesterols surrounding rhodopsin and opsin, in lipid bilayers mimicking those found naturally. Our results agree with those found experimentally and in previous simulations, but with far better statistical certainty. The results demonstrate the value of combining all-atom and coarse-grained models with experiment to provide a well-rounded view of lipid-protein interactions.

Retinal conformation governs pKa of protonated Schiff base in rhodopsin activation

We have explored the relationship between conformational energetics and the protonation state of the Schiff base in retinal, the covalently bound ligand responsible for activating the G protein-coupled receptor rhodopsin, using quantum chemical calculations. Guided by experimental structural determinations and large-scale molecular simulations on this system, we examined rotation about each bond in the retinal polyene chain, for both the protonated and deprotonated states that represent the dark and photoactivated states, respectively. Particular attention was paid to the torsional degrees of freedom that determine the shape of the molecule, and hence its interactions with the protein binding pocket. While most torsional degrees of freedom in retinal are characterized by large energetic barriers that minimize structural fluctuations under physiological temperatures, the C6-C7 dihedral defining the relative orientation of the β-ionone ring to the polyene chain has both modest barrier heights and a torsional energy surface that changes dramatically with protonation of the Schiff base. This surprising coupling between conformational degrees of freedom and protonation state is further quantified by calculations of the pKa as a function of the C6-C7 dihedral angle. Notably, pKa shifts of greater than two units arise from torsional fluctuations observed in molecular dynamics simulations of the full ligand-protein-membrane system. It follows that fluctuations in the protonation state of the Schiff base occur prior to forming the activated MII state. These new results shed light on important mechanistic aspects of retinal conformational changes that are involved in the activation of rhodopsin in the visual process.

UV-visible and infrared methods for investigating lipid-rhodopsin membrane interactions

We describe experimental UV-visible and Fourier transform infrared (FTIR) spectroscopic methods for characterizing lipid-protein interactions for rhodopsin in a membrane bilayer environment. The combination of FTIR and UV-visible difference spectroscopy is used to monitor the structural and functional changes during rhodopsin activation. Investigations of how membrane lipids stabilize various rhodopsin photoproducts are analogous to mutating the protein in terms of gain or loss of function. Interpretation of the results entails a flexible surface model for explaining membrane lipid-protein interactions through material properties relevant to biological activity.

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 2H NMR relaxation illuminates functional dynamics of retinal cofactor in membrane activation of rhodopsin

Rhodopsin is a canonical member of the family of G protein-coupled receptors, which transmit signals across cellular membranes and are linked to many drug interventions in humans. Here we show that solid-state (2)H NMR relaxation allows investigation of light-induced changes in local ps-ns time scale motions of retinal bound to rhodopsin. Site-specific (2)H labels were introduced into methyl groups of the retinal ligand that are essential to the activation process. We conducted solid-state (2)H NMR relaxation (spin-lattice, T(1Z), and quadrupolar-order, T(1Q)) experiments in the dark, Meta I, and Meta II states of the photoreceptor. Surprisingly, we find the retinylidene methyl groups exhibit site-specific differences in dynamics that change upon light excitation--even more striking, the C9-methyl group is a dynamical hotspot that corresponds to a crucial functional hotspot of rhodopsin. Following 11-cis to trans isomerization, the (2)H NMR data suggest the β-ionone ring remains in its hydrophobic binding pocket in all three states of the protein. We propose a multiscale activation mechanism with a complex energy landscape, whereby the photonic energy is directed against the E2 loop by the C13-methyl group, and toward helices H3 and H5 by the C5-methyl of the β-ionone ring. Changes in retinal structure and dynamics initiate activating fluctuations of transmembrane helices H5 and H6 in the Meta I-Meta II equilibrium of rhodopsin. Our proposals challenge the Standard Model whereby a single light-activated receptor conformation yields the visual response--rather an ensemble of substates is present, due to the entropy gain produced by photolysis of the inhibitory retinal lock.

Reconstitution of rhodopsin into polymerizable planar supported lipid bilayers: influence of dienoyl monomer structure on photoactivation

G-protein-coupled receptors (GPCRs) play key roles in cellular signal transduction and many are pharmacologically important targets for drug discovery. GPCRs can be reconstituted in planar supported lipid bilayers (PSLBs) with retention of activity, which has led to development of GPCR-based biosensors and biochips. However, PSLBs composed of natural lipids lack the high stability desired for many technological applications. One strategy is to use synthetic lipid monomers that can be polymerized to form robust bilayers. A key question is how lipid polymerization affects GPCR structure and activity. Here we have investigated the photochemical activity of bovine rhodopsin (Rho), a model GPCR, reconstituted into PSLBs composed of lipids having one or two polymerizable dienoyl moieties located in different regions of the acyl chains. Plasmon waveguide resonance spectroscopy was used to compare the degree of Rho photoactivation in fluid and poly(lipid) PSLBs. The position of the dienoyl moiety was found to have a significant effect: polymerization near the glycerol backbone significantly attenuates Rho activity whereas polymerization near the acyl chain termini does not. Differences in cross-link density near the acyl chain termini also do not affect Rho activity. In unpolymerized PSLBs, an equimolar mixture of phosphatidylethanolamine and phosphatidylcholine (PC) lipids enhances activity relative to pure PC; however after polymerization, the enhancement is eliminated which is attributed to stabilization of the membrane lamellar phase. These results should provide guidance for the design of robust lipid bilayers functionalized with transmembrane proteins for use in membrane-based biochips and biosensors.

A glance at the structural and functional diversity of membrane lipids

In the postgenomic era, spatially and temporally regulated molecular interactions as signals are beginning to take center stage in the understanding of fundamental biological events. For years, reductionism derived from the "fluid mosaic" model of the cell membrane has portrayed membrane lipids as rather passive molecules that, whereas separating biologically relevant aqueous phases, provided an environment so that membrane proteins could fulfill the specificity and selectivity required for proper cell signaling. Whereas these roles for membrane lipids still stand, the structural diversity of lipids and their complex arrangement in supramolecular assemblies have expanded such limited, although fundamental roles. Growing developments in the field of membrane lipids help to understand biological phenomena at the nanoscale domain, and reveal this heterogeneous group of organic compounds as a long underestimated group of key regulatory molecules. In this introductory chapter, brief overviews of the structural diversity of membrane lipids, the impact of different lipids on membrane properties, the vertical organization of lipids into rafts and caveolae, and the functional role of lipids as mediators of inter- and intracellular signals are provided. Any comprehensive review on membrane lipids, whether emphasizing structural or functional aspects, will require several volumes. The purpose of this chapter is to provide both introduction and rationale for the selection of topics that lie ahead in this book. For this reason, the list of references primarily includes reviews on particular issues dealing with membrane lipids wherein the reader can find further references.

Solid-state 2H NMR spectroscopy of retinal proteins in aligned membranes

Solid-state 2H NMR spectroscopy gives a powerful avenue to investigating the structures of ligands and cofactors bound to integral membrane proteins. For bacteriorhodopsin (bR) and rhodopsin, retinal was site-specifically labeled by deuteration of the methyl groups followed by regeneration of the apoprotein. 2H NMR studies of aligned membrane samples were conducted under conditions where rotational and translational diffusion of the protein were absent on the NMR time scale. The theoretical lineshape treatment involved a static axial distribution of rotating C-C2H3 groups about the local membrane frame, together with the static axial distribution of the local normal relative to the average normal. Simulation of solid-state 2H NMR lineshapes gave both the methyl group orientations and the alignment disorder (mosaic spread) of the membrane stack. The methyl bond orientations provided the angular restraints for structural analysis. In the case of bR the retinal chromophore is nearly planar in the dark- and all-trans light-adapted states, as well upon isomerization to 13-cis in the M state. The C13-methyl group at the "business end" of the chromophore changes its orientation to the membrane upon photon absorption, moving towards W182 and thus driving the proton pump in energy conservation. Moreover, rhodopsin was studied as a prototype for G protein-coupled receptors (GPCRs) implicated in many biological responses in humans. In contrast to bR, the retinal chromophore of rhodopsin has an 11-cis conformation and is highly twisted in the dark state. Three sites of interaction affect the torsional deformation of retinal, viz. the protonated Schiff base with its carboxylate counterion; the C9-methyl group of the polyene; and the beta-ionone ring within its hydrophobic pocket. For rhodopsin, the strain energy and dynamics of retinal as established by 2H NMR are implicated in substituent control of activation. Retinal is locked in a conformation that is twisted in the direction of the photoisomerization, which explains the dark stability of rhodopsin and allows for ultra-fast isomerization upon absorption of a photon. Torsional strain is relaxed in the meta I state that precedes subsequent receptor activation. Comparison of the two retinal proteins using solid-state 2H NMR is thus illuminating in terms of their different biological functions.

Structural analysis and dynamics of retinal chromophore in dark and meta I states of rhodopsin from 2H NMR of aligned membranes

Rhodopsin is a prototype for G protein-coupled receptors (GPCRs) that are implicated in many biological responses in humans. A site-directed (2)H NMR approach was used for structural analysis of retinal within its binding cavity in the dark and pre-activated meta I states. Retinal was labeled with (2)H at the C5, C9, or C13 methyl groups by total synthesis, and was used to regenerate the opsin apoprotein. Solid-state (2)H NMR spectra were acquired for aligned membranes in the low-temperature lipid gel phase versus the tilt angle to the magnetic field. Data reduction assumed a static uniaxial distribution, and gave the retinylidene methyl bond orientations plus the alignment disorder (mosaic spread). The dark-state (2)H NMR structure of 11-cis-retinal shows torsional twisting of the polyene chain and the beta-ionone ring. The ligand undergoes restricted motion, as evinced by order parameters of approximately 0.9 for the spinning C-C(2)H(3) groups, with off-axial fluctuations of approximately 15 degrees . Retinal is accommodated within the rhodopsin binding pocket with a negative pre-twist about the C11=C12 double bond that explains its rapid photochemistry and the trajectory of 11-cis to trans isomerization. In the cryo-trapped meta I state, the (2)H NMR structure shows a reduction of the polyene strain, while torsional twisting of the beta-ionone ring is maintained. Distortion of the retinal conformation is interpreted through substituent control of receptor activation. Steric hindrance between trans retinal and Trp265 can trigger formation of the subsequent activated meta II state. Our results are pertinent to quantum and molecular mechanics simulations of ligands bound to GPCRs, and illustrate how (2)H NMR can be applied to study their biological mechanisms of action.

Effects of lipid composition on the membrane activity and lipid phase behaviour of Vibrio sp. DSM14379 cells grown at various NaCl concentrations

The membrane lipid composition of living cells generally adjusts to the prevailing environmental and physiological conditions. In this study, membrane activity and lipid composition of the Gram-negative bacterium Vibrio sp. DSM14379, grown aerobically in a peptone-yeast extract medium supplemented with 0.5, 1.76, 3, 5 or 10% (w/v) NaCl, was determined. The ability of the membrane to reduce a spin label was studied by EPR spectroscopy under different salt concentrations in cell suspensions labeled with TEMPON. For lipid composition studies, cells were harvested in a late exponential phase and lipids were extracted with chloroform-methanol-water, 1:2:0.8 (v/v). The lipid polar head group and acyl chain compositions were determined by thin-layer and gas-liquid chromatographies. (31)P-NMR spectroscopy was used to study the phase behaviour of the cell lipid extracts with 20 wt.% water contents in a temperature range from -10 to 50 degrees C. The results indicate that the ability of the membrane to reduce the spin label was highest at optimal salt concentrations. The composition of both polar head groups and acyl chains changed markedly with increasing salinity. The fractions of 16:0, 16:1 and 18:0 acyl chains increased while the fraction of 18:1 acyl chains decreased with increasing salinity. The phosphatidylethanolamine fraction correlated inversely with the lysophosphatidylethanolamine fraction, with phosphatidylethanolamine exhibiting a minimum, and lysophosphatidylethanolamine a maximum, at the optimum growth rate. The fraction of lysophosphatidylethanolamine was surprisingly high in the lipid extracts. This lipid can form normal micellar and hexagonal phases and it was found that all lipid extracts form a mixture of lamellar and normal isotropic liquid crystalline phases. This is an anomalous behaviour since the nonlamellar phases formed by total lipid extracts are generally of the reversed type.

How lipids affect the activities of integral membrane proteins

The activities of integral membrane proteins are often affected by the structures of the lipid molecules that surround them in the membrane. One important parameter is the hydrophobic thickness of the lipid bilayer, defined by the lengths of the lipid fatty acyl chains. Membrane proteins are not rigid entities, and deform to ensure good hydrophobic matching to the surrounding lipid bilayer. The structure of the lipid headgroup region is likely to be important in defining the structures of those parts of a membrane protein that are located in the lipid headgroup region. A number of examples are given where the conformation of the headgroup-embedded region of a membrane protein changes during the reaction cycle of the protein; activities of such proteins might be expected to be particularly sensitive to lipid headgroup structure. Differences in hydrogen bonding potential and hydration between the headgroups of phosphatidycholines and phosphatidylethanolamines could be important factors in determining the effects of these lipids on protein activities, as well as any effects related to the tendency of the phosphatidylethanolamines to form a curved, hexagonal H(II) phase. Effects of lipid structure on protein aggregation and helix-helix interactions are also discussed, as well as the effects of charged lipids on ion concentrations close to the surface of the bilayer. Interpretations of lipid effects in terms of changes in protein volume, lipid free volume, and curvature frustration are also described. Finally, the role of non-annular, or 'co-factor' lipids, tightly bound to membrane proteins, is described.

Lipids do influence protein function-the hydrophobic matching hypothesis revisited

A topical review of the current state of lipid-protein interactions is given with focus on the physical interactions between lipids and integral proteins in lipid-bilayer membranes. The concepts of hydrophobic matching and curvature stress are revisited in light of recent data obtained from experimental and theoretical studies which demonstrate that not only do integral proteins perturb the lipids, but the physical state of the lipids does also actively influence protein function. The case of the trans-membrane water-channel protein aquaporin GlpF from E. coli imbedded in lipid-bilayer membranes is discussed in some detail. Numerical data obtained from Molecular Dynamics simulations show on the one side that the lipid bilayer adapts to the channel by a hydrophobic matching condition which reflects the propensity of the lipid molecules for forming curved structures. On the other side, it is demonstrated that the transport function of the channel is modulated by the matching condition and/or the curvature stress in a lipid-specific manner.

Regulation of sodium channel function by bilayer elasticity: the importance of hydrophobic coupling. Effects of Micelle-forming amphiphiles and cholesterol

Membrane proteins are regulated by the lipid bilayer composition. Specific lipid–protein interactions rarely are involved, which suggests that the regulation is due to changes in some general bilayer property (or properties). The hydrophobic coupling between a membrane-spanning protein and the surrounding bilayer means that protein conformational changes may be associated with a reversible, local bilayer deformation. Lipid bilayers are elastic bodies, and the energetic cost of the bilayer deformation contributes to the total energetic cost of the protein conformational change. The energetics and kinetics of the protein conformational changes therefore will be regulated by the bilayer elasticity, which is determined by the lipid composition. This hydrophobic coupling mechanism has been studied extensively in gramicidin channels, where the channel–bilayer hydrophobic interactions link a “conformational” change (the monomer↔dimer transition) to an elastic bilayer deformation. Gramicidin channels thus are regulated by the lipid bilayer elastic properties (thickness, monolayer equilibrium curvature, and compression and bending moduli). To investigate whether this hydrophobic coupling mechanism could be a general mechanism regulating membrane protein function, we examined whether voltage-dependent skeletal-muscle sodium channels, expressed in HEK293 cells, are regulated by bilayer elasticity, as monitored using gramicidin A (gA) channels. Nonphysiological amphiphiles (β-octyl-glucoside, Genapol X-100, Triton X-100, and reduced Triton X-100) that make lipid bilayers less “stiff”, as measured using gA channels, shift the voltage dependence of sodium channel inactivation toward more hyperpolarized potentials. At low amphiphile concentration, the magnitude of the shift is linearly correlated to the change in gA channel lifetime. Cholesterol-depletion, which also reduces bilayer stiffness, causes a similar shift in sodium channel inactivation. These results provide strong support for the notion that bilayer–protein hydrophobic coupling allows the bilayer elastic properties to regulate membrane protein function.

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.

Molecular Dynamics Simulation of Dark-adapted Rhodopsin in an Explicit Membrane Bilayer: Coupling between Local Retinal and Larger Scale Conformational Change

The light-driven photocycle of rhodopsin begins the photoreceptor cascade that underlies visual response. In a sequence of events, the retinal covalently attached to the rhodopsin protein undergoes a conformational change that communicates local changes to a global conformational change throughout the whole protein. In turn, the large-scale protein change then activates G-proteins and signal amplification throughout the cell. The nature of this change, involving a coupling between a local process and larger changes throughout the protein, may be important for many membrane proteins. In addition, functional work has shown that this coupling occurs with different efficiency in different lipid settings. To begin to understand the nature of the efficiency of this coupling in different lipid settings, we present a molecular dynamics study of rhodopsin in an explicit dioleoyl-phosphatidylcholine bilayer. Our system was simulated for 40 ns and provides insights into the very early events of the visual cascade, before the full transition and activation have occurred. In particular, we see an event near 10 ns that begins with a change in hydrogen bonding near the retinal and that leads through a series of coupled changes to a shift in helical tilt. This type of event, though rare on the molecular dynamics time-scale, could be an important clue to the types of coupling that occur between local and large-scale conformational change in many membrane proteins.

Determination of L(alpha)-H(II) phase transition temperature for 1,2-dioleoyl-sn-glycero-3-phosphatidylethanolamine

The thermodynamic properties of fully-hydrated lipids provide important information about the stability of membranes and the energetic interactions of lipid bilayers with membrane proteins (Nagle and Scott, Physics Today, 2:39, 1978). The lamellar/inverse hexagonal (L(alpha)-H(II)) phase transition of 1,2-dioleoyl-sn-glycero-3-phosphatidylethanolamine (DOPE) water mixtures is a first-order transition and, therefore, at constant pressure, must have a thermodynamically well-defined equilibrium transition temperature. The observed transition temperature is known to be dependent upon the rate at which the temperature is changed, which accounts for the many different values in the literature. X-ray diffraction was used to study the phase transition of fully-hydrated DOPE to determine the rate-independent transition temperature, T(LH). Samples were heated or cooled for a range of rates, 0.212 < r < 225 degrees C/hr, and the rate-dependent apparent phase transition temperatures, T(A)(r) were determined from the x-ray data. By use of a model-free extrapolation method, the transition temperature was found to be T(LH) = 3.33 +/- 0.16 degrees C. The hysteresis, /T(A)(r) - T(LH)/, was identical for heating and cooling rates, +/-r, and varied as /r/beta for beta approximately 1/4. This unexpected power-law relationship is consistent with a previous study (Tate et al., Biochemistry, 31:1081-1092, 1992) but differs markedly from the exponential behavior of activation barrier kinetics. The methods used in this study are general and provide a simple way to determine the true mesomorphic phase transition temperatures of other lipid and lyotropic systems.

Breaking the Meyer-Overton rule: predicted effects of varying stiffness and interfacial activity on the intrinsic potency of anesthetics

Exceptions to the Meyer-Overton rule are commonly cited as evidence against indirect, membrane-mediated mechanisms of general anesthesia. However, another interpretation is possible within the context of an indirect mechanism in which solubilization of an anesthetic in the membrane causes a redistribution of lateral pressures in the membrane, which in turn shifts the conformational equilibrium of membrane proteins such as ligand-gated ion channels. It is suggested that compounds of different stiffness and interfacial activity have different intrinsic potencies, i.e., they cause widely different redistributions of the pressure profile (and thus different effects on protein conformational equilibria) per unit concentration of the compound in the membrane. Calculations incorporating the greater stiffness of perfluoromethylenic chains and the large interfacial attraction of hydroxyl groups predict the higher intrinsic potency of short alkanols than alkanes, the cutoffs in potency of alkanes and alkanols and the much shorter cutoffs for their perfluorinated analogues. Both effects, increased stiffness and interfacial activity, are present in unsaturated hydrocarbon solutes, and the intrinsic potencies are predicted to depend on the magnitude of both effects and on the number and locations of multiple bonds within the molecule. Most importantly, the intrinsic potencies of polymeric alkanols with regularly spaced hydroxyl groups are predicted to rise with increasing chain length, without cutoff; such molecules should serve to distinguish unambiguously between indirect mechanisms and direct binding mechanisms of anesthesia.

Physical properties of liposomes and proteoliposomes prepared from Escherichia coli polar lipids

Reconstituted proteoliposomes serve as experimental systems for the study of membrane enzymes. Osmotic shifts and other changes in the solution environment may influence the structures and membrane properties of phospholipid vesicles (including liposomes, proteoliposomes and biological membrane vesicles) and hence the activities of membrane-associated proteins. Polar lipid extracts from Escherichia coli are commonly used in membrane protein reconstitution. The solution environment influenced the phase transition temperature and the diameter of liposomes and proteoliposomes prepared from E. coli polar lipid by extrusion. Liposomes prepared from E. coli polar lipids differed from dioleoylphosphatidylglycerol liposomes in Young's elastic modulus, yield point for solute leakage and structural response to osmotic shifts, the latter indicated by static light scattering spectroscopy. At high concentrations, NaCl caused aggregation of E. coli lipid liposomes that precluded detailed interpretation of light scattering data. Proteoliposomes and liposomes prepared from E. coli polar lipids were similar in size, yield point for solute leakage and structural response to osmotic shifts imposed with sucrose as osmolyte. These results will facilitate studies of bacterial enzymes implicated in osmosensing and of other enzymes that are reconstituted in E. coli lipid vesicles.

Curvature properties of novel forms of phosphatidylcholine with branched acyl chains

We studied the properties of a series of phosphatidylcholine molecules with branched acyl chains. These lipids have previously been shown to have marked stimulatory effects on the side-chain cleavage activity of cytochrome P450SCC (CYP11A1), an enzyme of the inner mitochondrial membrane. The synthetic lipids used were diacyl phosphatidylcholines with the decanoyl, dodecanoyl or tetradecanoyl chain having a hexyl, octyl or decyl straight chain aliphatic branch at the 2-position. All three lipids lowered the bilayer to hexagonal phase transition temperature of dielaidoyl phosphatidylethanolamine, the lipids with longer acyl chains being more effective in this regard. As pure lipids all of the forms were found by X-ray diffraction to be predominantly in the hexagonal phase (HII) over the entire temperature range of 7-75 degrees C. The properties of the HII phase were unusual with regard to the small size of the lattice spacings and the small temperature dependence of the spacings. We used tetradecane to relieve hydrocarbon packing constraints to determine the intrinsic radius of curvature of the lipid monolayer. The elastic bending modulus was measured in the presence of tetradecane by introducing an osmotic gradient across the hexagonal phase cylinders with aqueous solutions of poly(ethylene glycol). The elastic bending modulus was found to be higher than that observed with other lipids and to increase with temperature. Both the small intrinsic radius of curvature and the high elastic bending modulus indicate that the presence of these lipids in bilayer membranes will impose a high degree of negative curvature strain.

Solute modulation of conformational equilibria in intrinsic membrane proteins: apparent "cooperativity" without binding

The activity of many membrane proteins depends on a conformational transition that is often strongly influenced by small membrane-soluble solutes. This allosteric modulation may be direct, involving binding to the protein at localized sites of varying specificity, or may be indirect, resulting from altered membrane properties. In the present paper, a general expression for solute-protein titration curves is predicted, using an indirect mechanism that couples solute-induced changes in the lateral pressure profile of the bilayer to a shift in protein conformational equilibrium. When the common practice of fitting dose-response data to the Hill equation is applied to these curves, the fits are found to be reasonably good, with large Hill coefficients. Because this would commonly be interpreted as evidence of the existence of multiple sites with strong positive cooperativity, it is argued that caution must therefore be exercised in the interpretation of titration data in the absence of direct evidence of the existence of binding sites. The form of the titration curve predicted from this lateral pressure mechanism is shown to be quite general for indirect mechanisms. It is also shown that this form is the same as would be obtained from classical models of binding cooperativity, such as that of Monod, Wyman, and Changeux, in the limit of an infinite number of sites with vanishingly small site affinity.

The influence of membrane lateral pressures on simple geometric models of protein conformational equilibria

The function of many intrinsic membrane proteins requires a conformational transition that is often strongly influenced by the molecular composition of the bilayer in which the protein is embedded. Recently, a mechanism for this shift in conformational equilibrium was suggested, in which it is argued that a shift in distribution of lateral pressures of the bilayer resulting from a change in lipid composition alters the amount of mechanical work of the protein conformational transition, if the change in the cross-sectional area profile of the protein varies with depth within the bilayer. As there is little information on the change in shape of the transmembrane region of any protein, various simple geometric models are considered. For both a generic model, and more specific models that approximate likely cooperative rearrangements of alpha-helices in bundles, it is found that the conformational equilibrium depends on the first and second integral moments of the lateral pressure distribution. In addition to revealing the possible physical underpinnings of the well-known correlation between protein activity and the 'nonlamellar' tendency of bilayer lipids, this dependence on moments of the pressure profile allows for prediction of the relative effects of different lipid compositional changes even in the absence of information on specific protein shape changes. Effects of variation in acyl chain length, degree and position of cis-unsaturation, and addition of cholesterol and small interfacially-active solutes (n-alkanols) are compared.

Plasmon resonance spectroscopy: probing molecular interactions within membranes

Surface plasmon resonance (SPR) has become a popular method for investigating biomolecular interactions. A new variant of this technique, coupled plasmon-waveguide resonance (CPWR) spectroscopy, allows the characterization of anisotropic biological membranes. Plasmon resonance can therefore be used to study the molecular events involved in a wide variety of membrane processes, including energy conversion and signal transduction.

Lipid composition and the lateral pressure profile in bilayers

The mechanisms by which variations in the lipid composition of cell membranes influence the function of membrane proteins are not yet well understood. In recent work, a nonlocal thermodynamic mechanism was suggested in which changes in lipid composition cause a redistribution of lateral pressures that in turn modulates protein conformational (or aggregation) equilibria. In the present study, results of statistical thermodynamic calculations of the equilibrium pressure profile and bilayer thickness are reported for a range of lipids and lipid mixtures. Large redistributions of lateral pressure are predicted to accompany variation in chain length, degree and position of chain unsaturation, head group repulsion, and incorporation of cholesterol and interfacially active solutes. Combinations of compositional changes are found that compensate with respect to bilayer thickness, thus eliminating effects of hydrophobic mismatch, while still effecting significant shifts of the pressure profile. It is also predicted that the effect on the pressure profile of addition of short alkanols can be reproduced with certain unnatural lipids. These results suggest possible roles of cholesterol, highly unsaturated fatty acids and small solutes in modulating membrane protein function and suggest unambiguous experimental tests of the pressure profile hypothesis. As a test of the methodology, calculated molecular areas and area elastic moduli are compared with experimental and simulation results.

Lipid polymorphism and protein-lipid interactions

Non-lamellar-forming lipids play an important role in determining the physical properties of membranes. They affect the activity of membrane proteins and peptides. In addition, peptides which lyse membranes as well as those which promote membrane fusion facilitate the formation of non-lamellar phases, either micelles, cubic or hexagonal phases. The relationship of these diverse effects on membrane curvature is discussed in relation to the function of certain peptides and proteins. Specific examples of ionophoric peptides, cytotoxic peptides and viral fusion peptides are given. In addition, we compare the modulation of the rate of photoisomerisation of an integral membrane protein, rhodopsin, by non-lamellar-forming lipids with the effects of these lipids on an amphitropic protein, protein kinase C. Among these diverse systems it is frequently observed that the modulation of biological activity can be described in terms of the effect of the peptide or protein on the relative stability of lamellar and non-lamellar structures.