Dioxygen transmembrane distributions and partitioning thermodynamics in lipid bilayers and micelles

Cellular respiration, mediated by the passive diffusion of oxygen across lipid membranes, is key to many basic cellular processes. In this work, we report the detailed distribution of oxygen across lipid bilayers and examine the thermodynamics of oxygen partitioning via NMR studies of lipids in a small unilamellar vesicle (SUV) morphology. Dissolved oxygen gives rise to paramagnetic chemical shift perturbations and relaxation rate enhancements, both of which report on local oxygen concentration. From SUVs containing the phospholipid sn-2-perdeuterio-1-myristelaidoyl, 2-myristoyl-sn-glycero-3-phosphocholine (MLMPC), an analogue of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), we deduced the complete trans-bilayer oxygen distribution by measuring (13)C paramagnetic chemical shifts perturbations for 18 different sites on MLMPC arising from oxygen at a partial pressure of 30 bar. The overall oxygen solubility at 45 °C spans a factor of 7 between the bulk water (23.7 mM) and the bilayer center (170 mM) and is lowest in the vicinity of the phosphocholine headgroup, suggesting that oxygen diffusion across the glycerol backbone should be the rate-limiting step in diffusion-mediated passive transport of oxygen across the lipid bilayer. Lowering of the temperature from 45 to 25 °C gave rise to a slight decrease of the oxygen solubility within the hydrocarbon interior of the membrane. An analysis of the temperature dependence of the oxygen solubility profile, as measured by (1)H paramagnetic relaxation rate enhancements, reveals that oxygen partitioning into the bilayer is entropically favored (ΔS° = 54 ± 3 J K(-1) mol(-1)) and must overcome an enthalpic barrier (ΔH° = 12.0 ± 0.9 kJ mol(-1)).     

Structure and dynamics of helix-0 of the N-BAR domain in lipid micelles and bilayers

Bin/Amphiphysin/Rvs-homology (BAR) domains generate and sense membrane curvature by binding the negatively charged membrane to their positively charged concave surfaces. N-BAR domains contain an N-terminal extension (helix-0) predicted to form an amphipathic helix upon membrane binding. We determined the NMR structure and nano-to-picosecond dynamics of helix-0 of the human Bin1/Amphiphysin II BAR domain in sodium dodecyl sulfate and dodecylphosphocholine micelles. Molecular dynamics simulations of this 34-amino acid peptide revealed electrostatic and hydrophobic interactions with the detergent molecules that induce helical structure formation from residues 8-10 toward the C-terminus. The orientation in the micelles was experimentally confirmed by backbone amide proton exchange. The simulation and the experiment indicated that the N-terminal region is disordered, and the peptide curves to adopted the micelle shape. Deletion of helix-0 reduced tubulation of liposomes by the BAR domain, whereas the helix-0 peptide itself was fusogenic. These findings support models for membrane curving by BAR domains in which helix-0 increases the binding affinity to the membrane and enhances curvature generation.     

The transmembrane homotrimer of ADAM 1 in model lipid bilayers

Fertilin is a transmembrane protein heterodimer formed by the two subunits fertilin alpha and fertilin beta that plays an important role in sperm-egg fusion. Fertilin alpha and beta are members of the ADAM family, and contain each one transmembrane alpha-helix, and are termed ADAM 1 and ADAM 2, respectively. ADAM 1 is the subunit that contains a putative fusion peptide, and we have explored the possibility that the transmembrane alpha-helical domain of ADAM 1 forms homotrimers, in common with other viral fusion proteins. Although this peptide was found to form various homooligomers in SDS, the infrared dichroic data obtained with the isotopically labeled peptide at specific positions is consistent with the presence of only one species in DMPC or POPC lipid bilayers. Comparison of the experimental orientational data with molecular dynamics simulations performed with sequence homologues of ADAM 1 show that the species present in lipid bilayers is only consistent with an evolutionarily conserved homotrimeric model for which we provide a backbone structure. These results support a model where ADAM 1 forms homotrimers as a step to create a fusion active intermediate.     

Effect of pathogenic cysteine mutations on FGFR3 transmembrane domain dimerization in detergents and lipid bilayers

Mutations in fibroblast growth factor receptors are known as the genetic basis of skeletal growth disorders. The mechanism of pathogenesis, as determined by mutation-induced changes in receptor structure, interactions, and function, is elusive. Here we study three pathogenic Cys mutations, associated with either thanatophoric dysplasia or achondroplasia, in the TM domain of fibroblast growth factor receptors 3 (FGFR3). We characterize the dimerization propensities of the mutant TM domains in detergents and in lipid bilayers, in the presence and absence of reducing agents, and compare them to previous measurements of wild-type. We find that the Cys mutations increase the propensity for dimerization in detergent, with the Cys370 mutant exhibiting the highest propensity for disulfide bond formation, the Cys371 mutant having an intermediate propensity, and Cys375 the lowest. Thus, disulfide bonds readily form in detergents, with efficiency that correlates with the severity of the phenotype. In lipid bilayers, however, the Cys370 mutant, which dimerizes strongly in detergent, behaves as the wild-type, suggesting that Cys370-mediated disulfide bonds do not form between the isolated TM domains in bilayers. Thus, the nature of the hydrophobic environment plays an important role in defining the structure and flexibility of transmembrane dimers. These results and previous findings from cellular studies lead us to propose a conformational flexibility mechanism of receptor stabilization as a basis for disregulated FGFR3 signaling in thanatophoric dysplasia and achondroplasia.     

Structure and plasticity of the human immunodeficiency virus gp41 fusion domain in lipid micelles and bilayers

A thorough understanding of the structure of fusion domains of enveloped viruses in changing lipid environments helps us to formulate mechanistic models on how they might function in mediating viral entry by membrane fusion. We have expressed the N-terminal fusion domain of HIV-1 gp41 as a construct that is water-soluble in the absence of membranes, but that also binds with high affinity to lipid micelles and bilayers in their presence. We have solved the structure and studied the dynamics of this domain bound to dodecylphosphocholine micelles by homo- and heteronuclear NMR spectroscopy. The fusion peptide forms a stable hydrophobic helix from Ile(4) to Ala(14), but is increasingly more disordered and dynamic in a segment of intermediate polarity that stretches from Ala(15) to Ser(23). When bound to lipid bilayers at low concentration, the HIV fusion domain is also largely alpha-helical, as determined by CD and FTIR spectroscopy. However, at higher protein/lipid ratios, the domain is partially converted to form beta-structures in lipid bilayers. Controlled lipid mixing occurs at concentrations that support the alpha-helical, but not the beta-strand conformation.     

Hydrophobic mismatch and lipid sorting near OmpA in mixed bilayers: atomistic and coarse-grained simulations

To understand the effects of lipid composition on membrane protein function in a mixture as complex as a biomembrane, one must know whether the lipid composition local to the protein differs from the mean lipid composition. In this study, we simulated the transmembrane domain of a β-barrel protein, OmpA, in mixtures of lipids of different tail lengths under conditions of negative hydrophobic mismatch, i.e., local bilayer thinning. We modeled the influence of OmpA on the local lipid composition both at a coarse-grained (CG) resolution using conventional molecular dynamics, and at an atomistic resolution within the semi-grand canonical ensemble using mutation moves to rapidly approach an equilibrium lateral distribution of lipids. Moderate enrichment of the shorter tail component (either DDPC in DDPC/DMPC mixtures or DMPC in DMPC/DSPC mixtures) extending 2-3 nm away from the protein surface was observed with both the atomistic and CG models. The similarity in trends suggests that the more computationally economical CG models capture the essential features of lipid sorting induced by hydrophobic mismatch.     

Quantification of helix-helix binding affinities in micelles and lipid bilayers

A theoretical approach for estimating association free energies of alpha-helices in nonpolar media has been developed. The parameters of energy functions have been derived from DeltaDeltaG values of mutants in water-soluble proteins and partitioning of organic solutes between water and nonpolar solvents. The proposed approach was verified successfully against three sets of published data: (1) dissociation constants of alpha-helical oligomers formed by 27 hydrophobic peptides; (2) stabilities of 22 bacteriorhodopsin mutants, and (3) protein-ligand binding affinities in aqueous solution. It has been found that coalescence of helices is driven exclusively by van der Waals interactions and H-bonds, whereas the principal destabilizing contributions are represented by side-chain conformational entropy and transfer energy of atoms from a detergent or lipid to the protein interior. Electrostatic interactions of alpha-helices were relatively weak but important for reproducing the experimental data. Immobilization free energy, which originates from restricting rotational and translational rigid-body movements of molecules during their association, was found to be less than 1 kcal/mole. The energetics of amino acid substitutions in bacteriorhodopsin was complicated by specific binding of lipid and water molecules to cavities created in certain mutants.     

A transmembrane peptide from the human EGF receptor: behaviour of the cytoplasmic juxtamembrane domain in lipid bilayers

Solid state (2)H NMR spectroscopy was employed to study peptides related to the transmembrane domain of the human epidermal growth factor receptor, for insight into the interaction of its cytoplasmic juxtamembrane domain with the membrane surface. Since such receptors have clusters of (+)charged amino acids in this region, the effect of (-)charged phosphatidylserine at the concentration found naturally in the cytoplasmic leaflet (15 mol%) was considered. Each peptide contained 34 amino acids, which included the hydrophobic 23 amino acid stretch thought to span the membrane and a ten amino acid segment beyond the 'cytoplasmic' surface. Non-perturbing deuterium probe nuclei were located within alanine side chains in intramembranous and extramembranous portions. (2)H NMR spectra were recorded at 35 degrees C and 65 degrees C in fluid lipid bilayers consisting of (zwitterionic) 1-palmitoyl-2-oleoylphosphatidylcholine, with and without 15 mol% (anionic) phosphatidylserine. The cationic extramembranous portion of the receptor backbone was found to be highly rotationally mobile on a time scale of 10(-4)-10(-5) s in both types of membrane - as was the alpha-helical intramembranous portion. Deuterium nuclei in alanine side chains (-CD(3)) detected modest changes in peptide backbone orientation and/or dynamics related to the presence of 1-stearoyl-2-oleoylphosphatidylserine: in the case of the extramembranous portion of the peptide these seemed related to lipid charge. Temperature effects on the peptide backbone external to the membrane were qualitatively different from effects on the helical transmembrane domain - likely reflecting the different physical constraints on these peptide regions and the greater flexibility of the extramembranous domain. Effects related to lipid charge could be detected in the spectrum of CD(3) groups on the internally mobile side chain of Val(650), six residues beyond the membrane surface.     

Myelin basic protein ability to organize lipid bilayers: structural transition in bilayers of lysophosphatidylcholine micelles

Myelin basic protein isolated by a single step with the cationic detergent cethyltrimethylammonium bromide in a lipid-bound form is able to induce structural transition of lysophosphatydilcholine micelles into multi-laminar vesicles. This finding, observed through electron microscopy, is discussed in the light of the assumed ability of the basic protein to organize myelin lipids.     

Fluorescence quenching in lipid micelles and bilayers: Evaluation of bimolecular rate constants

Dynamic quenching of fluorophores and quenchers in lipid micelles and bilayers can yield information about the bimolecular rate constant for the quenching reaction, and hence information about the microviscosity of the fluorophore-quencher environment. When the fluorophore and quencher have relatively fixed transverse positions in the bilayer, the analysis of Sikaris et al. (Chem. Phys. Lipids. 29 (1981) 23) can be used to separate the microviscosity and proximity contributions to quenching. We now extend this method to show explicitly the effect of static quenching on the analysis. We show by simulation and experiment that a correction factor must be included when static quenching contributes to the observed quenching efficiency.     

Organization of model helical peptides in lipid bilayers: insight into the behavior of single-span protein transmembrane domains

Selectively deuterated transmembrane peptides comprising alternating leucine-alanine subunits were examined in fluid bilayer membranes by solid-state nuclear magnetic resonance (NMR) spectroscopy in an effort to gain insight into the behavior of membrane proteins. Two groups of peptides were studied: 21-mers having a 17-amino-acid hydrophobic domain calculated to be close in length to the hydrophobic thickness of 1-palmitoyl-2-oleoyl phosphatidylcholine and 26-mers having a 22-amino-acid hydrophobic domain calculated to exceed the membrane hydrophobic thickness. (2)H NMR spectral features similar to ones observed for transmembrane peptides from single-span receptors of higher animal cells were identified which apparently correspond to effectively monomeric peptide. Spectral observations suggested significant distortion of the transmembrane alpha-helix, and/or potential for restriction of rotation about the tilted helix long axis for even simple peptides. Quadrupole splittings arising from the 26-mer were consistent with greater peptide "tilt" than were those of the analogous 21-mer. Quadrupole splittings associated with monomeric peptide were relatively insensitive to concentration and temperature over the range studied, indicating stable average conformations, and a well-ordered rotation axis. At high peptide concentration (6 mol% relative to phospholipid) it appeared that the peptide predicted to be longer than the membrane thickness had a particular tendency toward reversible peptide-peptide interactions occurring on a timescale comparable with or faster than approximately 10(-5) s. This interaction may be direct or lipid-mediated and was manifest as line broadening. Peptide rotational diffusion rates within the membrane, calculated from quadrupolar relaxation times, T(2e), were consistent with such interactions. In the case of the peptide predicted to be equal to the membrane thickness, at low peptide concentration spectral lineshape indicated the additional presence of a population of peptide having rotational motion that was restricted on a timescale of 10(-5) s.     

The transmembrane domain of Neu in a lipid bilayer: molecular dynamics simulations

The results of full-atom molecular dynamics simulations of the transmembrane domains (TMDs) of both native, and Glu⁶⁶⁴-mutant (either protonated or unprotonated) Neu in an explicit fully hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayer are presented. For the native TMD peptide, a 10.05 ns trajectory was collected, while for the mutant TMD peptides 5.05 ns trajectories were collected for each. The peptides in all three simulations display stable predominantly -helical hydrogen bonding throughout the trajectories. The only significant exception occurs near the C-terminal end of the native and unprotonated mutant TMDs just outside the level of the lipid headgroups, where -helical hydrogen bonding develops, introducing a kink in the backbone structure. However, there is no indication of the formation of a bulge within the hydrophobic region of either native or mutant peptides. Over the course of the simulation of the mutant peptide, it is found that a significant number of water molecules penetrate the hydrophobic region of the surrounding lipid molecules, effectively hydrating Glu⁶⁶⁴. If the energy cost of such water penetration is significant enough, this may be a factor in the enhanced dimerization affinity of Glu⁶⁶⁴-mutant Neu.     

Structure and dynamics of the gammaM4 transmembrane domain of the acetylcholine receptor in lipid bilayers: insights into receptor assembly and function

A 28-mer peptide (gammaM4) corresponding to the fourth transmembrane segment of the nicotinic acetylcholine receptor (AChR) gamma-subunit, with a single tryptophan residue (Trp6), was reconstituted into lipid bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), loaded with either high or low amounts of cholesterol, i.e., in the conjugated liquid-ordered and liquid-disordered phases, respectively, at room temperature. By making use of the Trp intrinsic fluorescence, both steady-state and time-resolved fluorescence techniques were employed, namely, red-edge excitation shift effect, decay-associated spectra (DAS), and time-resolved anisotropy. The results obtained here, together with previous studies on the same reconstituted peptide, indicate that: (i) Trp6 is strongly anchored in the bilayer with a defined transverse location; (ii) the modifications in the measured DAS are related to the complex result of a self-quenching process on the decay parameters; (iii) the wobbling movement of the indole moiety of Trp6 is fast but severely restricted in amplitude; and, (iv) in the liquid-ordered phase, the bilayer properties and the tilt angle of the peptide enhance peptide-peptide interactions, with the formation of peptide rich patches and possibly some anti-parallel helix-helix aggregates, showing different dynamics from that of the peptide in the liquid-disordered phase where the peptide is randomly distributed.     

The molecular reorganization of lipid bilayers by osmium tetroxide. A spin-label study of orientation and restricted y-axis anisotropic motion in model membrane systems

Phospholipid dispersions spontaneously form oriented lamellar multilayers when dried onto glass slides. These oriented multilayers form useful model systems for studying the molecular dynamics of lipid bilayers. In order to examine the effects of osmium tetroxide on the orientation and motion of hydrocarbon chains in lipid bilayers, lecithin multilayers containing the spin label 3-doxyl-5α-cholestane (the 4′,4′-dimethyloxazolidine-N-oxyl derivative of 5α-cholestan-3-one) were prepared and examined by electron spin resonance spectroscopy. In egg lecithin multilayers at room temperature and 81% relative humidity the osmium tetroxide causes nearly complete loss of orientation and severe reduction of molecular motion. In contrast, the high degree of order in l-α-dipalmitoyl lecithin multilayers is not affected by exposure to osmium tetroxide vapors. Experiments are also reported on macroscopically disordered lecithin preparations, and the data support the conclusions drawn from the ordered lecithin multilayers that rotational mobility of the probe is severely reduced by fixation of the lipid chains.A simple mathematical model has been developed to account for the amplitude of the high-frequency (τ < 10^?8 sec) restricted y-axis anisotropic motion occurring in the bilayer plane. Since the y-axis is roughly parallel to the molecular axis of the rigid steroid spin label, this model enables quantitative comparisons of various degrees of restricted motion about the molecular axis.     

Spatial structure and dimer--monomer equilibrium of the ErbB3 transmembrane domain in DPC micelles

In present work the interaction of two TM α-helices of the ErbB3 receptor tyrosine kinase from the ErbB or HER family (residues 639-670) was studied by means of NMR spectroscopy in a membrane-mimicking environment provided by the DPC micelles. The ErbB3 TM segment appeared to form a parallel symmetric dimer in a left-handed orientation. The interaction between TM spans is accomplished via the non-standard motif and is supported by apolar contacts of bulky side chains and by stacking of aromatic rings together with π-cation interactions of Phe and Arg side chains. The investigation of the dimer--monomer equilibrium revealed thermodynamic properties of the assembly and the presence of two distinct regimes of the dimerization at low and at high peptide/detergent ratio. It was found that the detergent in case of ErbB3 behaves not as an ideal solvent, thus affecting the dimer--monomer equilibrium. Such behavior may account for the problems occurring with the refolding and stability of multispan helical membrane proteins in detergent solutions. The example of ErbB3 allows us to conclude that the thermodynamic parameters of dimerization, measured in micelles for two different helical pairs, cannot be compared without the investigation of their dependence on detergent concentration.     

Topological stability and self-association of a completely hydrophobic model transmembrane helix in lipid bilayers

Investigation of interactions between hydrophobic model peptides and lipid bilayers is perhaps the only way to elucidate the principles of the folding and stability of membrane proteins (White, S. H., and Wimley, W. C. (1998) Biochim. Biophys. Acta 1367, 339-352). We designed the completely hydrophobic "inert" peptide modeling a transmembrane (TM) helix without any of the specific side-chain interactions expected, X-(LALAAAA)(3)-NH(2) [X = Ac (I), 7-nitro-2-1,3-benzoxadiazol-4-yl (II), or 5(6)-carboxytetramethylrhodamine (III)]. Fourier transform infrared-polarized attenuated total reflection measurements revealed that I as well as II assume a TM helix in hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayers. Dithionite quenching experiments detected no topological change (flip-flop) in the helix II for at least 24 h. Thus, the TM helix itself is a highly stable structure, even in the absence of flanking hydrophilic or aromatic amino acids which are suggested to play important roles in stable TM positioning. Helix self-association in lipid bilayers was detected by fluorescence resonance energy transfer between II and III. The peptide was in a monomer-antiparallel dimer equilibrium with an association free energy of approximately -13 kJ/mol. Electron spin resonance spectra of 1-palmitoyl-2-stearoyl-(14-doxyl)-sn-glycero-3-phosphocholine demonstrated the presence of a motionally restricted component at lower temperatures.     

Expression and characterization of the FXYD ion transport regulators for NMR structural studies in lipid micelles and lipid bilayers

The proteins PLM (phospholemman), CHIF (channel inducing factor), and Mat8 (mammary tumor protein 8 kDa) are members of the FXYD family of ion transport regulatory membrane proteins. Here we describe their cloning and expression in Escherichia coli, and their purification for NMR structural studies in lipid micelles and lipid bilayers. The molecular masses of the purified recombinant FXYD proteins, determined from SDS-PAGE and from MALDI TOF mass spectrometry, reflect monomeric species. The solution NMR and CD spectra in SDS micelles show that they adopt helical conformations. The solid-state NMR spectra in lipid bilayers give the first view of their transmembrane architecture.     

Micropolarities of lipid bilayers and micelles 5. Localization of pyrene in small unilamellar phosphatidylcholine vesicles

The excimer/monomer ratio of emission intensities (I_E/I_M) and the enhancement of the 0-0 vibronic transition in the fluorescence spectra of pyrene (PY) and 16-(1-pyrenyl)hexadecanoic acid (C₁₆PY) were used to investigate the localization of PY in the bilayers of small unilamellar vesicles constituted of phosphatidylcholine (SUV-PC). First, from comparison of the fluorescence characteristics of PY in water with those of PY incorporated into the SUV-PC membranes, we concluded that the probe is incorporated preferentially in the lipid phase of the vesicles and not in the bulk aqueous phase. In addition, we found that, contrary to what happens with the pyrenyl moiety of C₁₆PY the location of PY varies with its relative concentration in the membrane space. The critical concentration was observed to be around 1.0 mol% of incorporated PY. At concentrations below this value, PY is located in the hydrocarbon core of the lipid bilayers. Above 1.0 mol%, the PY molecules reside preferentially in the neighbourhood of the glyceryl moiety region of the PC vesicles.     

The transmembrane protein KpOmpA anchoring the outer membrane of Klebsiella pneumoniae unfolds and refolds in response to tensile load

In Klebsiella pneumoniae the transmembrane β-barrel forming outer membrane protein KpOmpA mediates adhesion to a wide range of immune effector cells, thereby promoting respiratory tract and urinary infections. As major transmembrane protein OmpA stabilizes Gram-negative bacteria by anchoring their outer membrane to the peptidoglycan layer. Adhesion, osmotic pressure, hydrodynamic flow, and structural deformation apply mechanical stress to the bacterium. This stress can generate tensile load to the peptidoglycan-binding domain (PGBD) of KpOmpA. To investigate how KpOmpA reacts to mechanical stress, we applied a tensile load to the PGBD and observed a detailed unfolding pathway of the transmembrane β-barrel. Each step of the unfolding pathway extended the polypeptide connecting the bacterial outer membrane to the peptidoglycan layer and absorbed mechanical energy. After relieving the tensile load, KpOmpA reversibly refolded back into the membrane. These results suggest that bacteria may reversibly unfold transmembrane proteins in response to mechanical stress.