Bilayer deformation by the Kv channel voltage sensor domain revealed by self-assembly simulations

Coarse-grained molecular dynamics simulations are used to explore the interaction with a phospholipid bilayer of the voltage sensor (VS) domain and the S4 helix from the archaebacterial voltage-gated potassium (Kv) channel KvAP. Multiple 2-mus self-assembly simulations reveal that the isolated S4 helix may adopt either interfacial or transmembrane (TM) locations with approximately equal probability. In the TM state, the insertion of the voltage-sensing region of S4 is facilitated via local bilayer deformation that, combined with side chain "snorkeling," enables its Arg side chains to interact with lipid headgroups and water. Multiple 0.2-mus self-assembly simulations of the VS domain are also performed, along with simulations of MscL and KcsA, to permit comparison with more "canonical" integral membrane protein structures. All three stably adopt a TM orientation within a bilayer. For MscL and KcsA, there is no significant bilayer deformation. In contrast, for the VS, there is considerable local deformation, which is again primarily due to the lipid-exposed S4. It is shown that for both the VS and isolated S4 helix, the positively charged side chains of S4 are accommodated within the membrane through a combination of stabilizing interactions with lipid glycerol and headgroup regions, water, and anionic side chains. Our results support the possibility that bilayer deformation around key gating charge residues in Kv channels may result in "focusing" of the electrostatic field, and indicate that, when considering competing models of voltage-sensing, it is essential to consider the dynamics and structure of not only the protein but also of the local lipid environment.     

Side-chain hydrophobicity scale derived from transmembrane protein folding into lipid bilayers

The transfer free energies of the twenty natural amino acid side chains from water to phospholipid bilayers make a major contribution to the assembly and function of membrane proteins. Measurements of those transfer free energies will facilitate the identification of membrane protein sequences and aid in the understanding of how proteins interact with membranes during key biological events. We report the first water-to-bilayer transfer free energy scale (i.e., a "hydrophobicity scale") for the twenty natural amino acid side chains measured in the context of a native transmembrane protein and a phospholipid bilayer. Our measurements reveal parity for apolar side-chain contributions between soluble and membrane proteins and further demonstrate that an arginine side-chain placed near the middle of a lipid bilayer is accommodated with much less energetic cost than predicted by molecular dynamics simulations.     

Structure of the gramicidin A channel: Discrimination between the πL,D and the β helix by electrical measurements with lipid bilayer membranes

Measurements with different chemically modified gramicidins in lipid bilayer membranes were used to discriminate between the dimeric pi(L,D) helix proposed by Urry and the dimeric parallel or antiparallel helices proposed by Veatch and Blout. Evidence for the pi(L,D) helix was obtained on the basis of the different actions of a negatively charged O-pyromellitylgramicidin and a negatively charged N-pyromellityldesformylgramicidin on lipid bilayer membranes. O-Pyromellitylgramicidin forms ionic channels in lipid membranes when it is applied to both sides of the membrane. In contrast to unmodified gramicidin, O-pyromellitylgramicidin is inactive when it is applied only to one side of the membrane. N-Pyromellityldesformylgramicidin does not form ionic channels in lipid bilayer membranes whether it is applied to one or both sides of the membrane. These results support the view that the gramicidin channel is formed by two pi(L,D) helices. Dimer formation by head-to-head association of two pi(L,D) helices needs six intermolecular hydrogen bonds, which are located at the formyl end of the molecule and which occur deep within the lipid membrane. In the head-to-head associated pi(L,D) helix the absence of the formyl group leads to an inactivation of the peptide, whereas in a parallel or antiparallel double-stranded helix the absence of the formyl group should have only minor effects.     

Interface connections of a transmembrane voltage sensor

Voltage-sensitive ion channels open and close in response to changes in transmembrane (TM) potential caused by the motion of the S4 voltage sensors. These sensors are alpha-helices that include four or more positively charged amino acids, most commonly arginine. The so-called paddle model, based on the high-resolution structure of the KvAP K+ channel [Jiang, et al. (2003) Nature 423, 33-41], posits that the S4 sensors move within the membrane bilayer in response to TM voltage changes. Direct exposure of S4 sensors to lipid is contrary to the classical expectation that the dielectric contrast between the membrane hydrocarbon core and water presents an insurmountable energetic penalty to burial of electric charges. Nevertheless, recent experiments have shown that a helix with the sequence of KvAP S4 can be inserted across the endoplasmic reticulum membrane. To reconcile this result with the classical energetics argument, we have carried out a molecular dynamics simulation of an isolated TM S4 helix in a lipid bilayer. The simulation reveals a stabilizing hydrogen-bonded network of water and lipid phosphates around the arginines that reduces the effective thickness of the bilayer hydrocarbon core to approximately 10 A in the vicinity of the helix. It suggests that bilayer phospholipids can adapt locally to strongly perturbing protein elements, causing the phospholipids to become a structural extension of the protein.     

Direct measurement of salt-bridge solvation energies using a peptide model system: implications for protein stability.

The solvation energies of salt bridges formed between the terminal carboxyl of the host pentapeptide AcWL- X-LL and the side chains of Arg or Lys in the guest (X) position have been measured. The energies were derived from octanol-to-buffer transfer free energies determined between pH 1 and pH 9. 13C NMR measurements show that the salt bridges form in the octanol phase, but not in the buffer phase, when the side chains and the terminal carboxyl group are charged. The free energy of salt-bridge formation in octanol is approximately -4 kcal/mol (1 cal = 4.184 J), which is equal to or slightly larger than the sum of the solvation energies of noninteracting pairs of charged side chains. This is about one-half the free energy that would result from replacing a charge pair in octanol with a pair of hydrophobic residues of moderate size. Therefore, salt bridging in octanol can change the favorable aqueous solvation energy of a pair of oppositely charged residues to neutral or slightly unfavorable but cannot provide the same free energy decrease as hydrophobic residues. This is consistent with recent computational and experimental studies of protein stability.     

On the energetics of translocon-assisted insertion of charged transmembrane helices into membranes

The understanding of the mechanism of insertion of transmembrane (TM) helixes through the translocon presents a major open challenge. Although the experimental information about the partition of the inserted helices between the membrane and the solution contains crucial information about this process, it is not clear how to extract this information. In particular, it is not clear how to rationalize the small apparent insertion energy, ΔG_app, of an ionized residue in the center of a TM helix. Here we explore the nature of the insertion energies, asking what should be the value of these parameters if their measurements represent equilibrium conditions. This is done using a coarse-grained model with advanced electrostatic treatment. Estimating the energetics of ionized arginine of a TM helix in the presence of neighboring helixes or the translocon provides a rationale for the observed ΔG_app of ionized residues. It is concluded that the apparent insertion free energy of TM with charged residues reflects probably more than just the free energy of moving the isolate single helix from water into the membrane. The present approach should be effective not only in exploring the mechanism of the operation of the translocon but also for studies of other membrane proteins.     

Lipid interaction of the C terminus and association of the transmembrane segments facilitate atlastin-mediated homotypic endoplasmic reticulum fusion

The homotypic fusion of endoplasmic reticulum (ER) membranes is mediated by atlastin (ATL), which consists of an N-terminal cytosolic domain containing a GTPase module and a three-helix bundle followed by two transmembrane (TM) segments and a C-terminal tail (CT). Fusion depends on a GTP hydrolysis-induced conformational change in the cytosolic domain. Here, we show that the CT and TM segments also are required for efficient fusion and provide insight into their mechanistic roles. The essential feature of the CT is a conserved amphipathic helix. A synthetic peptide corresponding to the helix, but not to unrelated amphipathic helices, can act in trans to restore the fusion activity of tailless ATL. The CT promotes vesicle fusion by interacting directly with and perturbing the lipid bilayer without causing significant lysis. The TM segments do not serve as mere membrane anchors for the cytosolic domain but rather mediate the formation of ATL oligomers. Point mutations in either the C-terminal helix or the TMs impair ATL's ability to generate and maintain ER morphology in vivo. Our results suggest that protein-lipid and protein-protein interactions within the membrane cooperate with the conformational change of the cytosolic domain to achieve homotypic ER membrane fusion.     

Transmembrane orientation and possible role of the fusogenic peptide from parainfluenza virus 5 (PIV5) in promoting fusion

Membrane fusion is required for diverse biological functions ranging from viral infection to neurotransmitter release. Fusogenic proteins increase the intrinsically slow rate of fusion by coupling energetically downhill conformational changes of the protein to kinetically unfavorable fusion of the membrane-phospholipid bilayers. Class I viral fusogenic proteins have an N-terminal hydrophobic fusion peptide (FP) domain, important for interaction with the target membrane, plus a C-terminal transmembrane (C-term-TM) helical membrane anchor. The role of the water-soluble regions of fusogenic proteins has been extensively studied, but the contributions of the membrane-interacting FP and C-term-TM peptides are less well characterized. Typically, FPs are thought to bind to membranes at an angle that allows helix penetration but not traversal of the lipid bilayer. Here, we show that the FP from the paramyxovirus parainfluenza virus 5 fusogenic protein, F, forms an N-terminal TM helix, which self-associates into a hexameric bundle. This FP also interacts strongly with the C-term-TM helix. Thus, the fusogenic F protein resembles SNARE proteins involved in vesicle fusion by having water-soluble coiled coils that zipper during fusion and TM helices in both membranes. By analogy to mechanosensitive channels, the force associated with zippering of the water-soluble coiled-coil domain is expected to lead to tilting of the FP helices, promoting interaction with the C-term-TM helices. The energetically unfavorable dehydration of lipid headgroups of opposing bilayers is compensated by thermodynamically favorable interactions between the FP and C-term-TM helices as the coiled coils zipper into the membrane phase, leading to a pore lined by both lipid and protein.     

Arginine residues at internal positions in a protein are always charged

Many functionally essential ionizable groups are buried in the hydrophobic interior of proteins. A systematic study of Lys, Asp, and Glu residues at 25 internal positions in staphylococcal nuclease showed that their pK_a values can be highly anomalous, some shifted by as many as 5.7 pH units relative to normal pK_a values in water. Here we show that, in contrast, Arg residues at the same internal positions exhibit no detectable shifts in pK_a; they are all charged at pH ≤ 10. Twenty-three of these 25 variants with Arg are folded at both pH 7 and 10. The mean decrease in thermodynamic stability from substitution with Arg was 6.2 kcal/mol at this pH, comparable to that for substitution with Lys, Asp, or Glu at pH 7. The physical basis behind the remarkable ability of Arg residues to remain protonated in environments otherwise incompatible with charges is suggested by crystal structures of three variants showing how the guanidinium moiety of the Arg side chain is effectively neutralized through multiple hydrogen bonds to protein polar atoms and to site-bound water molecules. The length of the Arg side chain, and slight deformations of the protein, facilitate placement of the guanidinium moieties near polar groups or bulk water. This unique capacity of Arg side chains to retain their charge in dehydrated environments likely contributes toward the important functional roles of internal Arg residues in situations where a charge is needed in the interior of a protein, in a lipid bilayer, or in similarly hydrophobic environments.     

Anionic lipid headgroups as a proton-conducting pathway along the surface of membranes: a hypothesis.

Evidence has been gathering from several laboratories that protons in proton-pumping membranes move along or within the bilayer rather than exchange with the bulk phase. These experiments are typically conducted on the natural membrane in vivo or in vitro or on fragments of natural membrane. Anionic lipids are present in all proton-pumping membranes. Model studies on the protonation state of the fatty acids of liposomes containing entrapped water show that the bilayers always contain mixtures of protonated and deprotonated carboxylates. Protonated fatty acids form stable acid-anion pairs with deprotonated fatty acids through unusually strong hydrogen bonds. Such acid-anion dimers have a single negative charge, which is shared by the four negative oxygens of both headgroups. The two pK values of the resulting dimer will be significantly different from the pK of the monomeric species, so that the dimer will be stable over a wide pH range. It is proposed that anionic lipid headgroups in biological membranes share protons as acid-anion dimers and that anionic lipids thus trap and conduct protons along the headgroup domain of bilayers that contain such anionic lipids. Protons pumped from the other side of the membrane may enter and move within the headgroup sheet because the protonation rate of negatively charged proton acceptors is 5 orders of magnitude faster than that of water. Protons trapped in the acidic headgroup sheet need not leave this region in order to be utilized by a responsive proton-translocating pore (a transport protein using the proton gradient). Experiments suggest the proton concentration in the headgroup domain may vary widely and the anionic lipid headgroup sheet may therefore function as a proton buffer. Due to the Gouy-Chapman-Stern layer at polyanionic surfaces, anionic lipids will also sequester protons from the bulk solution at low and moderate ionic strengths. At high ionic strength metal cations may replace protons sequestered near the headgroups, but these cations cannot substitute for protons in the "proton-conducting pathway," which is based on hydrogen bonding.     

Molecular dynamics simulations suggest a mechanism for translocation of the HIV-1 TAT peptide across lipid membranes

The recombinant HIV-1 Tat protein contains a small region corresponding to residues ⁴⁷YGRKKRRQRR⁵⁷R, which is capable of translocating cargoes of different molecular sizes, such as proteins, DNA, RNA, or drugs, across the cell membrane in an apparently energy-independent manner. The pathway that these peptides follow for entry into the cell has been the subject of strong controversy for the last decade. This peptide is highly basic and hydrophilic. Therefore, a central question that any candidate mechanism has to answer is how this highly hydrophilic peptide is able to cross the hydrophobic barrier imposed by the cell membrane. We propose a mechanism for the spontaneous translocation of the Tat peptides across a lipid membrane. This mechanism involves strong interactions between the Tat peptides and the phosphate groups on both sides of the lipid bilayer, the insertion of charged side chains that nucleate the formation of a transient pore, followed by the translocation of the Tat peptides by diffusing on the pore surface. This mechanism explains how key ingredients, such as the cooperativity among the peptides, the large positive charge, and specifically the arginine amino acids, contribute to the uptake. The proposed mechanism also illustrates the importance of membrane fluctuations. Indeed, mechanisms that involve large fluctuations of the membrane structure, such as transient pores and the insertion of charged amino acid side chains, may be common and perhaps central to the functions of many membrane protein functions.     

Membrane physical properties influence transmembrane helix formation

The pHLIP peptide has three states: (I) soluble in aqueous buffer, (II) bound to the bilayer surface at neutral pH, and (III) inserted as a transmembrane (TM) helix at acidic pH. The membrane insertion of pHLIP at low pH can be used to target the acidic tissues characteristic of different diseases, such as cancer. We find that the α-helix content of state II depends on lipid acyl chain length but not cholesterol, suggesting the helicity of the bound state may be controlled by the bilayer elastic bending modulus. Experiments with the P20G variant show the proline residue in pHLIP reduces the α-helix content of both states II and III. We also observe that the membrane insertion pKa is influenced by membrane physical properties, following a biphasic pattern similar to the membrane thickness optima observed for the function of eukaryotic membrane proteins. Because tumor cells exhibit altered membrane fluidity, we suggest this might influence pHLIP tumor targeting. We used a cell insertion assay to determine the pKa in live cells, observing that the properties in liposomes held in the more complex plasma membrane. Our results show that the formation of a TM helix is modulated by both the conformational propensities of the peptide and the physical properties of the bilayer. These results suggest a physical role for helix-membrane interactions in optimizing the function of more complex TM proteins.     

Structural fluctuations of a helical polypeptide traversing a lipid bilayer.

Time-resolved fluorescence anisotropy (FA) measurements are reported for five helical bilayer-spanning henicosapeptides, each containing one tryptophan at sequence position 1, 6, 11, 16, or 21. The FA decay reflects two molecular processes in all cases: local internal fluctuations of the tryptophan side chain with a relaxation time of 200-500 ps, and motions of the whole polypeptide molecule with a relaxation time of 9-10 ns. The amplitudes of the fast fluctuation are largest at the helix ends and decrease toward the center of the helix. A similar observation was made for the helical polypeptides dissolved in organic solvents showing that the mobility gradient along the polypeptide sequence is an inherent property of the polypeptide helix and not due to the lipid bilayer. However, the amplitudes of the fast fluctuations can be modulated by the physical state of the lipid bilayer. With decreasing temperature, the internal mobility of the tryptophan in all positions decreases with an abrupt change at the lipid-phase transition. Concomitant molecular dynamics calculations indicate a correlation between the fast FA decay and conformational fluctuations within the helix. According to the simulation, the conformation of the polypeptide is on average predominantly helical, but actually the molecule can fluctuate between a variety of different substructures. The conformational fluctuations are largest at the helix ends and provide the free space required for rotation of the indole ring around the tryptophan side chain bonds.     

Membrane protein thermodynamic stability may serve as the energy sink for sorting in the periplasm

Thermodynamic stabilities are pivotal for understanding structure–function relationships of proteins, and yet such determinations are rare for membrane proteins. Moreover, the few measurements that are available have been conducted under very different experimental conditions, which compromises a straightforward extraction of physical principles underlying stability differences. Here, we have overcome this obstacle and provided structure–stability comparisons for multiple membrane proteins. This was enabled by measurements of the free energies of folding and the m values for the transmembrane proteins PhoP/PhoQ-activated gene product (PagP) and outer membrane protein W (OmpW) from Escherichia coli. Our data were collected in the same lipid bilayer and buffer system we previously used to determine those parameters for E. coli outer membrane phospholipase A (OmpLA). Biophysically, our results suggest that the stabilities of these proteins are strongly correlated to the water-to-bilayer transfer free energy of the lipid-facing residues in their transmembrane regions. We further discovered that the sensitivities of these membrane proteins to chemical denaturation, as judged by their m values, was consistent with that previously observed for water-soluble proteins having comparable differences in solvent exposure between their folded and unfolded states. From a biological perspective, our findings suggest that the folding free energies for these membrane proteins may be the thermodynamic sink that establishes an energy gradient across the periplasm, thus driving their sorting by chaperones to the outer membranes in living bacteria. Binding free energies of these outer membrane proteins with periplasmic chaperones support this energy sink hypothesis.     

Monoclinic uncomplexed double-stranded, antiparallel, left-handed beta 5.6-helix (increases decreases beta 5.6) structure of gramicidin A: alternate patterns of helical association and deformation.

A comparison of the monoclinic and orthorhombic crystal structures of the uncomplexed double-stranded, antiparallel, left-handed beta-helix (5.6 amino acid residues per turn) (increases decreases beta 5.6) conformers of gramicidin A reveals marked differences in the tryptophan side-chain orientations and the degree of helical uniformity of the dimer and in the manner in which these helical dimers associate with one another in the crystal. The helix of the orthorhombic dimer exhibits a regular pattern of bulges and constrictions that appears to be induced by crystal packing forces affecting tryptophan side chains that are aligned parallel to the helix axis. The monoclinic dimer is more uniform than the orthorhombic dimer as a consequence of pi stacking interactions between dimers in which orientation of tryptophan side chains is normal to the helix axis to relieve the lateral crystal packing forces that may locally twist and deform the helix. It may be inferred from these observations that lipid interactions may be expected to destabilize the increases decreases beta 5.6 helix when it is inserted into a membrane bilayer.     

Multiple topogenic sequences in bovine opsin.

Bovine opsin, a polytopic integral membrane protein, contains seven transmembrane segments connecting eight hydrophilic domains alternating on each side of the membrane. To localize topogenic sequences that might specify the distinct topology of opsin in the membrane, we constructed various opsin mutants, each containing only one transmembrane segment. Messenger RNAs transcribed from these mutants were translated in a cell-free system supplemented with microsomal membranes. Among six of the seven transmembrane segments of opsin that were analyzed, five were able to function as signal sequences and also expressed stop-transfer sequences of variable strength. By the criteria of extractability at pH 11 and protease sensitivity, the presence of a signal sequence in combination with a "strong" stop-transfer sequence yielded integration into the lipid bilayer of the majority of chains. However, in combination with a "weak" stop-transfer sequence, we observed integration into the lipid bilayer of only some chains, with the others either completely translocated across the membrane or retained in a water-accessible space in the membrane.     

Protein area occupancy at the center of the red blood cell membrane

In the Fluid Mosaic Model for biological membrane structure, proposed by Singer and Nicolson in 1972, the lipid bilayer is represented as a neutral two-dimensional solvent in which the proteins of the membrane are dispersed and distributed randomly. The model portrays the membrane as dominated by a membrane lipid bilayer, directly exposed to the aqueous environment, and only occasionally interrupted by transmembrane proteins. This view is reproduced in virtually every textbook in biochemistry and cell biology, yet some critical features have yet to be closely examined, including the key parameter of the relative occupancy of protein and lipid at the center of a natural membrane. Here we show that the area occupied by protein and lipid at the center of the human red blood cell (RBC) plasma membrane is at least approximately 23% protein and less than approximately 77% lipid. This measurement is in close agreement with previous estimates for the RBC plasma membrane and the recently published measurements for the synaptic vesicle. Given that transmembrane proteins are surrounded by phospholipids that are perturbed by their presence, the occupancy by protein of more than approximately 20% of the RBC plasma membrane and the synaptic vesicle plasma membrane implies that natural membrane bilayers may be more rigid and less fluid than has been thought for the past several decades, and that studies of pure lipid bilayers do not fully reveal the properties of lipids in membranes. Thus, it appears to be the case that membranes may be more mosaic than fluid, with little unperturbed phospholipid bilayer.     

A method for the modulation of membrane fluidity: homogeneous catalytic hydrogenation of phospholipids and phospholipids and phospholipid-water model biomembranes.

The fatty acids associated with phospholipids of cell membranes, and particularly their degree of unsaturation, contribute to the fluidity of their structure and hance determine many of their biological properties. We describe a technique for modulating membrane fluidity which consists of hydrogenating the unsaturated double bonds of membrane phospholipids. This has been accomplished using a homogeneous catalyst. The process has been applied to phospholipids in organic solvents, to phospholipids dispersed as multibilayers in aqueous systems, and also to sonicated preparations of phospholipids arranged as single bilayer vesicles. Preliminary experiments have also been performed with biological membranes. These results indicate that the process of homogeneous catalysis for the modulation of lipid fluidity of biological cell membranes may have considerable future biological and biochemical application.     

Intramembrane cavitation as a unifying mechanism for ultrasound-induced bioeffects

The purpose of this study was to develop a unified model capable of explaining the mechanisms of interaction of ultrasound and biological tissue at both the diagnostic nonthermal, noncavitational (<100 mW · cm(-2)) and therapeutic, potentially cavitational (>100 mW · cm(-2)) spatial peak temporal average intensity levels. The cellular-level model (termed "bilayer sonophore") combines the physics of bubble dynamics with cell biomechanics to determine the dynamic behavior of the two lipid bilayer membrane leaflets. The existence of such a unified model could potentially pave the way to a number of controlled ultrasound-assisted applications, including CNS modulation and blood-brain barrier permeabilization. The model predicts that the cellular membrane is intrinsically capable of absorbing mechanical energy from the ultrasound field and transforming it into expansions and contractions of the intramembrane space. It further predicts that the maximum area strain is proportional to the acoustic pressure amplitude and inversely proportional to the square root of the frequency (ε A,max ∝ P(A)(0.8f - 0.5) and is intensified by proximity to free surfaces, the presence of nearby microbubbles in free medium, and the flexibility of the surrounding tissue. Model predictions were experimentally supported using transmission electron microscopy (TEM) of multilayered live-cell goldfish epidermis exposed in vivo to continuous wave (CW) ultrasound at cavitational (1 MHz) and noncavitational (3 MHz) conditions. Our results support the hypothesis that ultrasonically induced bilayer membrane motion, which does not require preexistence of air voids in the tissue, may account for a variety of bioeffects and could elucidate mechanisms of ultrasound interaction with biological tissue that are currently not fully understood.     

Hydrophobic association of alpha-helices, steric dewetting, and enthalpic barriers to protein folding

Efficient protein folding implies a microscopic funnel-like multidimensional free-energy landscape. Macroscopically, conformational entropy reduction can manifest itself as part of an empirical barrier in the traditional view of folding, but experiments show that such barriers can also entail significant unfavorable enthalpy changes. This observation raises the puzzling possibility, irrespective of conformational entropy, that individual microscopic folding trajectories may encounter large uphill moves and thus the multidimensional free-energy landscape may not be funnel-like. Here, we investigate how nanoscale hydrophobic interactions might underpin this salient enthalpic effect in biomolecular assembly by computer simulations of the association of two preformed polyalanine or polyleucine helices in water. We observe a high, positive enthalpic signature at room temperature when the helix separation is less than a single layer of water molecules. Remarkably, this unfavorable enthalpy change, with a parallel increase in void volume, is largely compensated for by a concomitant increase in solvent entropy, netting only a small or nonexistent microscopic free-energy barrier. Thus, our findings suggest that high enthalpic folding barriers can be consistent with a funnel picture of folding and are mainly a desolvation phenomenon indicative of a cooperative mechanism of simultaneous formation of multiple side-chain contacts at the rate-limiting step.