Buried lysine, but not arginine, titrates and alters transmembrane helix tilt

Ionization Properties of Histidine Residues in the Lipid Bilayer Membrane Environment

We address the critically important ionization properties of histidine side chains of membrane proteins, when exposed directly to lipid acyl chains within lipid bilayer membranes. The problem is important for addressing general principles that may underlie membrane protein function. To this end, we have employed a favorable host peptide framework provided by GWALP23 (acetyl-GGALW(5)LALALALALALALW(19)LAGA-amide). We inserted His residues into position 12 or 14 of GWALP23 (replacing either Leu(12) or Leu(14)) and incorporated specific [(2)H]Ala labels within the helical core sequence. Solid-state (2)H NMR spectra report the folding and orientation of the core sequence, revealing marked differences in the histidine-containing transmembrane helix behavior between acidic and neutral pH conditions. At neutral pH, the GWALP23-H12 and GWALP23-H14 helices exhibit well defined tilted transmembrane orientations in dioleoylphosphatidylcholine (DOPC)and dilauroylphosphatidylcholine (DLPC) bilayer membranes. Under acidic conditions, when His(12) is protonated and charged, the GWALP23-H12 helix exhibits a major population that moves to the DOPC bilayer surface and a minor population that occupies multiple transmembrane states. The response to protonation of His(14) is an increase in helix tilt, but GWALP23-H14 remains in a transmembrane orientation. The results suggest pKa values of less than 3 for His(12) and about 3-5 for His(14) in DOPC membranes. In the thinner DLPC bilayers, with increased water access, the helices are less responsive to changes in pH. The combined results enable us to compare the ionization properties of lipid-exposed His, Lys, and Arg side chains in lipid bilayer membranes.

Decrypting protein insertion through the translocon with free-energy calculations

Protein insertion into a membrane is a complex process involving numerous players. The most prominent of these players is the Sec translocon complex, a conserved protein-conducting channel present in the cytoplasmic membrane of bacteria and the membrane of the endoplasmic reticulum in eukaryotes. The last decade has seen tremendous leaps forward in our understanding of how insertion is managed by the translocon and its partners, coming from atomic-detailed structures, innovative experiments, and well-designed simulations. In this review, we discuss how experiments and simulations, hand-in-hand, teased out the secrets of the translocon-facilitated membrane insertion process. In particular, we focus on the role of free-energy calculations in elucidating membrane insertion. Amazingly, despite all its apparent complexity, protein insertion into membranes is primarily driven by simple thermodynamic and kinetic principles. This article is part of a Special Issue entitled: Membrane proteins edited by J.C. Gumbart and Sergei Noskov.

The Effect of Membrane Environment on Surfactant Protein C Stability Studied by Constant-pH Molecular Dynamics

Pulmonary surfactant protein C (SP-C) is a small peptide with two covalently linked fatty acyl chains that plays a crucial role in the formation and stabilization of the pulmonary surfactant reservoirs during the compression and expansion steps of the respiratory cycle. Although its function is known to be tightly related to its highly hydrophobic character and key interactions maintained with specific lipid components, much is left to understand about its molecular mechanism of action. Also, although it adopts a mainly helical structure while associated with the membrane, factors as pH variation and deacylation have been shown to affect its stability and function. In this work, the conformational behavior of both the acylated and deacylated SP-C isoforms was studied in a DPPC bilayer under different pH conditions using constant-pH molecular dynamics simulations. Our findings show that both protein isoforms are remarkably stable over the studied pH range, even though the acylated isoform exhibits a labile helix-turn-helix motif rarely observed in the other isoform. We estimate similar tilt angles for the two isoforms over the studied pH range, with a generally higher degree of internalization of the basic N-terminal residues in the deacylated case, and observe and discuss some protonation-conformation coupling effects. Both isoforms establish contacts with the surrounding lipid molecules (preferentially with the sn-2 ester bonds) and have a local effect on the conformational behavior of the surrounding lipid molecules, the latter being more pronounced for acylated SP-C.

Nanotube array method for studying lipid-induced conformational changes of a membrane protein by solid-state NMR

Anodic aluminum oxide substrates with macroscopically aligned homogeneous nanopores of 80 nm in diameter enable two-dimensional, solid-state nuclear magnetic resonance studies of lipid-induced conformational changes of uniformly (15)N-labeled Pf1 coat protein in native-like bilayers. The Pf1 helix tilt angles in bilayers composed of two different lipids are not entirely governed by the membrane thickness but could be rationalized by hydrophobic interactions of lysines at the bilayer interface. The anodic aluminum oxide alignment method is applicable to a broader repertoire of lipids versus bicelle bilayer mimetics currently employed in solid-state nuclear magnetic resonance of oriented samples, thus allowing for elucidation of the role played by lipids in shaping membrane proteins.

Membrane permeation of a peptide: it is better to be positive

A joint experimental and computational study investigates the translocation of a tryptophan molecule through a phospholipid membrane. Time dependent spectroscopy of the tryptophan side chain determines the rate of permeation into 150 nm phospholipid vesicles. Atomically detailed simulations are conducted to calculate the free energy profiles and the permeation coefficient. Different charging conditions of the peptide (positive, negative, or zwitterion) are considered. Both experiment and simulation reproduce the qualitative trend and suggest that the fastest permeation is when the tryptophan is positively charged. The permeation mechanism, which is revealed by molecular dynamics simulations, is of a translocation assisted by a local defect. The influence of long-range electrostatic interactions, such as the membrane dipole potential on the permeation process, is not significant.

Membrane environment modulates the pKa values of transmembrane helices

In this work, we apply the recently developed constant pH molecular dynamics technique to study protonation equilibria of titratable side chains in the context of simple transmembrane (TM) helices and explore the effect of pH on their configurations in membrane bilayers. We observe that, despite a significant shift toward neutral states, considerable population of different side chains stay in the charged state that give rise to pK_a values around 9.6 for Asp and Glu and 4.5 to 6 for His and Lys side chains, respectively. These charged states are highly stabilized by favorable interactions between head groups, water molecules, and the charged side chains that are facilitated by substantial changes in the configuration of the peptides. The pH dependent configurations and the measured pK_a values are in good agreement with relatively recent solid state NMR measurements. Our results presented here demonstrate that all-atom constant pH molecular dynamics can be applied to membrane proteins and peptides to obtain reliable pK_a values and pH dependent behavior for these systems.

Dynamic regulation of lipid-protein interactions

We review the importance of helix motions for the function of several important categories of membrane proteins and for the properties of several model molecular systems. For voltage-gated potassium or sodium channels, sliding, tilting and/or rotational movements of the S4 helix accompanied by a swapping of cognate side-chain ion-pair interactions regulate the channel gating. In the seven-helix G protein-coupled receptors, exemplified by the rhodopsins, collective helix motions serve to activate the functional signaling. Peptides which initially associate with lipid-bilayer membrane surfaces may undergo dynamic transitions from surface-bound to tilted-transmembrane orientations, sometimes accompanied by changes in the molecularity, formation of a pore or, more generally, the activation of biological function. For single-span membrane proteins, such as the tyrosine kinases, an interplay between juxtamembrane and transmembrane domains is likely to be crucial for the regulation of dimer assembly that in turn is associated with the functional responses to external signals. Additionally, we note that experiments with designed single-span transmembrane helices offer fundamental insights into the molecular features that govern protein-lipid interactions. This article is part of a Special Issue entitled: Lipid-protein interactions.

3D hydrophobic moment vectors as a tool to characterize the surface polarity of amphiphilic peptides

The interaction of membranes with peptides and proteins is largely determined by their amphiphilic character. Hydrophobic moments of helical segments are commonly derived from their two-dimensional helical wheel projections, and the same is true for β-sheets. However, to the best of our knowledge, there exists no method to describe structures in three dimensions or molecules with irregular shape. Here, we define the hydrophobic moment of a molecule as a vector in three dimensions by evaluating the surface distribution of all hydrophilic and lipophilic regions over any given shape. The electrostatic potential on the molecular surface is calculated based on the atomic point charges. The resulting hydrophobic moment vector is specific for the instantaneous conformation, and it takes into account all structural characteristics of the molecule, e.g., partial unfolding, bending, and side-chain torsion angles. Extended all-atom molecular dynamics simulations are then used to calculate the equilibrium hydrophobic moments for two antimicrobial peptides, gramicidin S and PGLa, under different conditions. We show that their effective hydrophobic moment vectors reflect the distribution of polar and nonpolar patches on the molecular surface and the calculated electrostatic surface potential. A comparison of simulations in solution and in lipid membranes shows how the peptides undergo internal conformational rearrangement upon binding to the bilayer surface. A good correlation with solid-state NMR data indicates that the hydrophobic moment vector can be used to predict the membrane binding geometry of peptides. This method is available as a web application on http://www.ibg.kit.edu/HM/.

Spontaneous transmembrane helix insertion thermodynamically mimics translocon-guided insertion

The favorable transfer free energy for a transmembrane (TM) α-helix between the aqueous phase and lipid bilayer underlies the stability of membrane proteins. However, the connection between the energetics and process of membrane protein assembly by the Sec61/SecY translocon complex in vivo is not clear. Here, we directly determine the partitioning free energies of a family of designed peptides using three independent approaches: an experimental microsomal Sec61 translocon assay, a biophysical (spectroscopic) characterization of peptide insertion into hydrated planar lipid bilayer arrays, and an unbiased atomic-detail equilibrium folding-partitioning molecular dynamics simulation. Remarkably, the measured free energies of insertion are quantitatively similar for all three approaches. The molecular dynamics simulations show that TM helix insertion involves equilibrium with the membrane interface, suggesting that the interface may play a role in translocon-guided insertion.

Atomistic simulations of pore formation and closure in lipid bilayers

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.

The importance of membrane defects-lessons from simulations

The defects and pores within lipid membranes are scientifically interesting and have a number of biological applications. Although lipid bilayers are extremely thin hydrophobic barriers, just ∼3 nm thick, they include diverse chemistry and have complex structures. Bilayers are soft and dynamic, and as a result, they can bend and deform in response to different stimuli by means of structural changes in their component lipids. Though defects occur within these structures, their transience and small size have made it difficult to characterize them. However, with recent advances in computer power and computational modeling techniques, researchers can now use simulations as a powerful tool to probe the mechanism and energies of defect and pore formation in a number of situations. In this Account, we present results from our detailed molecular dynamics computer simulations of hydrophilic pores and related defects in lipid bilayers at an atomistic level. Electroporation can be used to increase the permeability of cellular membranes, with potential therapeutic applications. Atomistic simulations of electroporation have illustrated the molecular details of this process, including the importance of water dipole interactions at the water-membrane interface. Characterization of the lipid-protein interactions provides an important tool for understanding transmembrane protein structure and thermodynamic stability. Atomistic simulations give a detailed picture of the free energies of model peptides and side chains in lipid membranes; the energetic cost of defect formation strongly influences the energies of interactions between lipids and polar and charged residues. Many antimicrobial peptides form hydrophilic pores in lipid membranes, killing bacteria or cancer cells. On the basis of simulation data, at least some of these peptides form defects and pores near the center of the bilayer, with a common disordered structure where hydrated headgroups form an approximately toroidal shape. The localization and trafficking of lipids supports general membrane structure and a number of important signaling cascades, such as those involving ceramide, diacylglycerol, and cholesterol. Atomistic simulations have determined the rates and free energies of lipid flip-flop. During the flip-flop of most phosphatidylcholine lipids, a hydrophilic pore forms when the headgroup moves near the center of the bilayer. Simulations have provided novel insight into many features of defects and pores in lipid membranes. Simulation data from very different systems and models show how water penetration and defect formation can determine the free energies of many membrane processes. Bilayers can deform and allow transient defects and pores when exposed to a diverse range of stimuli. Future work will explore many aspects of membrane defects with increased resolution and scope, including the study of more complex lipid mixtures, membrane domains, and large-scale membrane remodeling. Such studies will examine processes including vesicle budding and fusion, non-bilayer lipid phases, and interactions between lipid bilayers and other biomolecules. Simulations provide information that complements experimental studies, allowing microscopic insight into experimental observations and suggesting novel hypotheses and experiments. These studies should enable a deeper understanding of the role of lipid bilayers in cellular biology and support the development of future lipid-based biotechnology.

Interaction of Alzheimer's β-amyloid peptides with cholesterol: mechanistic insights into amyloid pore formation

Brain cholesterol plays a critical role in Alzheimer's disease and other neurodegenerative diseases. The molecular mechanisms linking cholesterol to neurotoxicity have remained elusive for a long time, but recent data have allowed the identification of functional cholesterol-binding domains in several amyloidogenic proteins involved in neurodegenerative diseases, including Alzheimer's disease. In this review, we analyze the cholesterol binding properties of β-amyloid (Aβ) peptides and the impact of these interactions on amyloid pore formation. We show that although the cholesterol-binding domains of Aβ peptides and of transmembrane precursor C99 are partially overlapping, they involve distinct amino acid residues, so that cholesterol has a greater affinity for Aβ than for C99. Synthetic 22-35 and 25-35 fragments of Aβ retained the ability of the full-length peptide 1-42 to bind cholesterol and to form zinc-sensitive, calcium-permeable amyloid pores in cultured neural cells. Studies with mutant peptides allowed the identification of key residues involved in cholesterol binding and channel formation. Cholesterol promoted the insertion of Aβ in the plasma membrane, induced α-helical structuration, and forced the peptide to adopt a tilted topology that favored the oligomerization process. Bexarotene, an amphipathic drug currently considered as a potential candidate medication for the treatment of neurodegenerative diseases, competed with cholesterol for binding to Aβ and prevented oligomeric channel formation. These studies indicate that it is possible to prevent the generation of neurotoxic oligomers by targeting the cholesterol-binding domain of Aβ peptides. This original strategy could be used for the treatment of Alzheimer's and other neurodegenerative diseases that involve cholesterol-dependent toxic oligomers.

On the design of supramolecular assemblies made of peptides and lipid bilayers

Peptides confer interesting properties to materials, supramolecular assemblies and to lipid membranes and are used in analytical devices or within delivery vehicles. Their relative ease of production combined with a high degree of versatility make them attractive candidates to design new such products. Here, we review and demonstrate how CD- and solid-state NMR spectroscopic approaches can be used to follow the reconstitution of peptides into membranes and to describe some of their fundamental characteristics. Whereas CD spectroscopy is used to monitor secondary structure in different solvent systems and thereby aggregation properties of the highly hydrophobic domain of p24, a protein involved in vesicle trafficking, solid-state NMR spectroscopy was used to deduce structural information and the membrane topology of a variety of peptide sequences found in nature or designed. (15)N chemical shift solid-state NMR spectroscopy indicates that the hydrophobic domain of p24 as well as a designed sequence of 19 hydrophobic amino acid residues adopt transmembrane alignments in phosphatidylcholine membranes. In contrast, the amphipathic antimicrobial peptide magainin 2 and the designed sequence LK15 align parallel to the bilayer surface. Additional angular information is obtained from deuterium solid-state NMR spectra of peptide sites labelled with (2)H3-alanine, whereas (31)P and (2)H solid-state NMR spectra of the lipids furnish valuable information on the macroscopic order and phase properties of the lipid matrix. Using these approaches, peptides and reconstitution protocols can be elaborated in a rational manner, and the analysis of a great number of peptide sequences is reviewed. Finally, a number of polypeptides with membrane topologies that are sensitive to a variety of environmental conditions such as pH, lipid composition and peptide-to-lipid ratio will be presented.

Comparisons of interfacial Phe, Tyr, and Trp residues as determinants of orientation and dynamics for GWALP transmembrane peptides

Aromatic amino acids often flank the transmembrane alpha helices of integral membrane proteins. By favoring locations within the membrane–water interface of the lipid bilayer, aromatic residues Trp, Tyr, and sometimes Phe may serve as anchors to help stabilize a transmembrane orientation. In this work, we compare the influence of interfacial Trp, Tyr, or Phe residues upon the properties of tilted helical transmembrane peptides. For such comparisons, it has been critical to start with no more than one interfacial aromatic residue near each end of a transmembrane helix, for example, that of GWALP23 (acetyl-GGALW⁵(LA)₆LW¹⁹LAGA-[ethanol]amide). To this end, we have employed ²H-labeled alanines and solid-state NMR spectroscopy to investigate the consequences of moving or replacing W5 or W19 in GWALP23 with selected Tyr, Phe, or Trp residues at the same or proximate locations. We find that GWALP23 peptides having F5, Y5, or W5 exhibit essentially the same average tilt and similar dynamics in bilayer membranes of 1,2-dilauroylphosphatidylcholine (DLPC) or 1,2-dioleoylphosphatidylcholine (DOPC). When double Tyr anchors are present, in Y^4,5GWALP23 the NMR observables are markedly more subject to dynamic averaging and at the same time are less responsive to the bilayer thickness. Decreased dynamics are nevertheless observed when ring hydrogen bonding is removed, such that F^4,5GWALP23 exhibits a similar extent of low dynamic averaging as GWALP23 itself. When F5 is the sole aromatic group in the N-interfacial region, the dynamic averaging is (only) slightly more extensive than with W5, Y5, or Y4 alone or with F4,5, yet it is much less than that observed for Y^4,5GWALP23. Interestingly, moving Y5 to Y4 or W19 to W18, while retaining only one hydrogen-bond-capable aromatic ring at each interface, maintains the low level of dynamic averaging but alters the helix azimuthal rotation. The rotation change is about 40° for Y4 regardless of whether the host lipid bilayer is DLPC or DOPC. The rotational change (Δρ) is more dramatic and more complex when W19 is moved to W18, as Δρ is about +90° in DLPC but about −60° in DOPC. Possible reasons for this curious lipid-dependent helix rotation could include not only the separation distances between flanking aromatic or hydrophobic residues but also the absolute location of the W19 indole ring. For the more usual cases, when the helix azimuthal rotation shows little dependence on the host bilayer identity, excepting W¹⁸GWALP23, the transmembrane helices adapt to different lipids primarily by changing the magnitude of their tilt. We conclude that, in the absence of other functional groups, interfacial aromatic residues determine the preferred orientations and dynamics of membrane-spanning peptides. The results furthermore suggest possibilities for rotational and dynamic control of membrane protein function.

Implicit membrane treatment of buried charged groups: application to peptide translocation across lipid bilayers

The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.

Recovery of a chemically synthesized Japanese encephalitis virus reveals two critical adaptive mutations in NS2B and NS4A

A full-length genome infectious clone is a powerful tool for functional assays in virology. In this study, using a chemical synthesized complete genome of Japanese encephalitis virus (JEV) strain SA14 (GenBank accession no. U14163), we constructed a full-length genomic cDNA clone of JEV. The recovered virus from the cDNA clone replicated poorly in baby hamster kidney (BHK-21) cells and in suckling mice brain. Following serial passage in BHK-21 cells, adaptive mutations within the NS2B and NS4A proteins were recovered in the passaged viruses leading to viruses with a large-plaque phenotype. Mutagenesis analysis, using a genome-length RNA and a replicon of JEV, demonstrated that the adaptive mutations restored replication to different degrees, and the restoration efficiencies were in the order: NS2B-T102M<NS4A-R79K<NS2B-T102M+NS4A-R79K. An in vivo virulence assay in mice showed that the recombinant virus containing double mutations showed similar virulence to the WT SA14 (GenBank accession no. M55506). This study reports the first chemically synthesized JEV. A reverse genetics assay demonstrated that substitutions of NS2B-T102M and NS4A-R79K altered JEV replication.

Effects of Phospholipid Composition on the Transfer of a Small Cationic Peptide Across a Model Biological Membrane

The transfer of a lysine amino acid analogue across phospholipid membrane models was investigated using molecular-dynamics simulations. The evolution of the protonation state of this small peptide as a function of its position inside the membrane was studied by determining the local pKa by means of free-energy calculations. Permeability and mean-first-passage time were evaluated and showed that the transfer occurs on the submillisecond time scale. Comparative studies were conducted to evaluate changes in the pKa arising from differences in the phospholipid chemical structure. We compared, hence, the effect of an ether vs an ester linkage of the lipid headgroup as well as linear vs branched lipid tails. The study reveals that protonated lysine residues can be buried further inside an ether lipid membrane than an ester lipid membrane, while branched lipids are found to stabilize less the charged form compared to their unbranched lipid chain counterparts.

Impact of histidine residues on the transmembrane helices of viroporins

Abstract The role of histidine in channel-forming transmembrane (TM) helices was investigated by comparing the TM helices from Virus protein 'u' (Vpu) and the M2 proton channel. Both proteins are members of the viroporin family of small membrane proteins that exhibit ion channel activity, and have a single TM helix that is capable of forming oligomers. The TM helices from both proteins have a conserved tryptophan towards the C-terminus. Previously, alanine 18 of Vpu was mutated to histidine in order to artificially introduce the same HXXXW motif that is central to the proton channel activity of M2. Interestingly, the mutated Vpu TM resulted in an increase in helix tilt angle of 11° in lipid bilayers compared to the wild-type Vpu TM. Here, we find the reverse, when histidine 37 of the HXXXW motif in M2 was mutated to alanine, it decreased the helix tilt by 10° from that of wild-type M2. The tilt change is independent of both the helix length and the presence of tryptophan. In addition, compared to wild-type M2, the H37A mutant displayed lowered sensitivity to proton concentration. We also found that the solvent accessibility of histidine-containing M2 is greater than without histidine. This suggests that the TM helix may increase the solvent exposure by changing its tilt angle in order to accommodate a polar/charged residue within the hydrophobic membrane region. The comparative results of M2, Vpu and their mutants demonstrated the significance of histidine in a transmembrane helix and the remarkable plasticity of the function and structure of ion channels stemming from changes at a single amino acid site.

Why do arginine and lysine organize lipids differently? Insights from coarse-grained and atomistic simulations

An important puzzle in membrane biophysics is the difference in the behaviors of lysine (Lys) and arginine (Arg) based peptides at the membrane. For example, the translocation of poly-Arg is orders of magnitude faster than that of poly-Lys. Recent experimental work suggests that much of the difference can be inferred from the phase behavior of peptide/lipid mixtures. At similar concentrations, mixtures of phosphatidylethanolamine (PE) and phosphatidylserine (PS) lipids display different phases in the presence of these polypeptides, with a bicontinuous phase observed with poly-Arg peptides and an inverted hexagonal phase observed with poly-Lys peptides. Here we show that simulations with the coarse-grained (CG) BMW-MARTINI model reproduce the experimental results. An analysis using atomistic and CG models reveals that electrostatic and glycerol-peptide interactions play a crucial role in determining the phase behavior of peptide-lipid mixtures, with the difference between Arg and Lys arising from the stronger interactions of the former with lipid glycerols. In other words, the multivalent nature of the guanidinium group allows Arg to simultaneously interact with both phosphate and glycerol groups, while Lys engages solely with phosphate; this feature of amino acid/lipid interactions has not been emphasized in previous studies. The Arg peptides colocalize with PS in regions of high negative Gaussian curvature and stabilize the bicontinuous phase. Decreasing the strength of either the electrostatic interactions or the peptide-glycerol interactions causes the inverted hexagonal phase to become more stable. The results highlight the utility of CG models for the investigation of phase behavior but also emphasize the subtlety of the phenomena, with small changes in specific interactions leading to qualitatively different phases.

The different interactions of lysine and arginine side chains with lipid membranes

The basic amino acids lysine (Lys) and arginine (Arg) play important roles in membrane protein activity, the sensing of membrane voltages, and the actions of antimicrobial, toxin, and cell-penetrating peptides. These roles are thought to stem from the strong interactions and disruptive influences of these amino acids on lipid membranes. In this study, we employ fully atomistic molecular dynamics simulations to observe, quantify, and compare the interactions of Lys and Arg with saturated phosphatidylcholine membranes of different thickness. We make use of both charged (methylammonium and methylguanidinium) and neutral (methylamine and methylguanidine) analogue molecules, as well as Lys and Arg side chains on transmembrane helix models. We find that the free energy barrier experienced by a charged Lys crossing the membrane is strikingly similar to that of a charged Arg (to within 2 kcal/mol), despite the two having different chemistries, H-bonding capability, and hydration free energies that differ by ∼10 kcal/mol. In comparison, the barrier for neutral Arg is higher than that for neutral Lys by around 5 kcal/mol, being more selective than that for the charged species. This can be explained by the different transport mechanisms for charged or neutral amino acid side chains in the membrane, involving membrane deformations or simple dehydration, respectively. As a consequence, we demonstrate that Lys would be deprotonated in the membrane, whereas Arg would maintain its charge. Our simulations also reveal that Arg attracts more phosphate and water in the membrane, and can form extensive H-bonding with its five H-bond donors to stabilize Arg-phosphate clusters. This leads to enhanced interfacial binding and membrane perturbations, including the appearance of a trans-membrane pore in a thinner membrane. These results highlight the special role played by Arg as an amino acid to bind to, disrupt, and permeabilize lipid membranes, as well as to sense voltages for a range of peptide and protein activities in nature and in engineered bionanodevices.

Mechanism of cholesterol-assisted oligomeric channel formation by a short Alzheimer β-amyloid peptide

Alzheimer β-amyloid (Aβ) peptides can self-organize into oligomeric ion channels with high neurotoxicity potential. Cholesterol is believed to play a key role in this process, but the molecular mechanisms linking cholesterol and amyloid channel formation have so far remained elusive. Here, we show that the short Aβ22-35 peptide, which encompasses the cholesterol-binding domain of Aβ, induces a specific increase of Ca(2+) levels in neural cells. This effect is neither observed in calcium-free medium nor in cholesterol-depleted cells, and is inhibited by zinc, a blocker of amyloid channel activity. Double mutations V24G/K28G and N27R/K28R in Aβ22-35 modify cholesterol binding and abrogate channel formation. Molecular dynamic simulations suggest that cholesterol induces a tilted α-helical topology of Aβ22-35. This facilitates the establishment of an inter-peptide hydrogen bond network involving Asn-27 and Lys-28, a key step in the octamerization of Aβ22-35 which proceeds gradually until the formation of a perfect annular channel in a phosphatidylcholine membrane. Overall, these data give mechanistic insights into the role of cholesterol in amyloid channel formation, opening up new therapeutic options for Alzheimer's disease. Aβ22-35 peptide, which encompasses the cholesterol binding domain of Aβ, induces a specific increase of Ca(2+) level in neural cells. Double mutations V24G/K28G and N27R/K28R modify cholesterol binding and abrogate channels formation. Molecular dynamic simulations suggest that cholesterol induces a tilted α-helical peptide topology facilitating the formation of annular octameric channels, as schematically shown in the graphic (with a hydrogen bond shown in green for two vicinal peptides). Overall, the data give insights into the role of cholesterol in amyloid channel formation and open up new therapeutic options for Alzheimer's disease.

Residues of a proposed gate region in type I ATP-binding cassette import systems are crucial for function as revealed by mutational analysis

The type I ATP-binding cassette (ABC) importer for positively charged amino acids of the thermophilic bacterium Geobacillus stearothermophilus consists of the extracellular solute binding protein, ArtJ, and a homodimer each of the transmembrane subunit, ArtM, and the nucleotide-binding and -hydrolyzing subunit, ArtP. We have investigated the functional consequences of mutations affecting conserved residues from two peptide regions in ArtM, recently proposed to form a 'gate' by which access of a substrate to the translocation path is controlled (Hollenstein et al., 2007 [14]). Transporter variants were reconstituted into proteoliposomes and assayed for ArtJ/arginine-stimulated ATPase activity. Replacement of residues from region 1 (Arg-63, Pro-66) caused no or only moderate reduction in ATPase activity. In contrast, mutating residues from gate region 2 (Lys-159, Leu-163) resulted in a substantial increase in ATPase activity which, however, as demonstrated for variants ArtM(K159I) and ArtM(K159E), is not coupled to transport. Replacing homologous residues in the closely related histidine transporter of Salmonella enterica serovar Typhimurium (HisJ-QMP2) caused different phenotypes. Mutation to isoleucine of HisQ(K163) or HisM(H172), both homologous to ArtM(K159), abolished ATPase activity. The mutations most likely caused a structural change as revealed by limited proteolysis. In contrast, substantial, albeit reduced, enzymatic activity was observed with variants of HisQ(L167→G) or HisM(L176→G), both homologous to ArtM(L163). Our study provides the first experimental evidence in favor of a crucial role of residues from the proposed gate region in type I ABC importer function.