Hydrophobic mismatch of mobile transmembrane helices: Merging theory and experiments

Hemoglobin Binding to the Red Blood Cell (RBC) Membrane Is Associated with Decreased Cell Deformability

The deformability of red blood cells (RBCs), expressing their ability to change their shape as a function of flow-induced shear stress, allows them to optimize oxygen delivery to the tissues and minimize their resistance to flow, especially in microcirculation. During physiological aging and blood storage, or under external stimulations, RBCs undergo metabolic and structural alterations, one of which is hemoglobin (Hb) redistribution between the cytosol and the membrane. Consequently, part of the Hb may attach to the cell membrane, and although this process is reversible, the increase in membrane-bound Hb (MBHb) can affect the cell's mechanical properties and deformability in particular. In the present study, we examined the correlation between the MBHb levels, determined by mass spectroscopy, and the cell deformability, determined by image analysis. Six hemoglobin subunits were found attached to the RBC membranes. The cell deformability was negatively correlated with the level of four subunits, with a highly significant inter-correlation between them. These data suggest that the decrease in RBC deformability results from Hb redistribution between the cytosol and the cell membrane and the respective Hb interaction with the cell membrane.

A Coronin 1-Derived Peptide Inhibits Membrane Fusion by Modulating Membrane Organization and Dynamics

Membrane fusion is considered the first step in the entry of enveloped viruses into the host cell. Several targeted strategies have been implemented to block viral entry by limiting the fusion protein to form a six-helix bundle, which is a prerequisite for fusion. Nonetheless, the development of broad-spectrum fusion inhibitors is essential to combat emerging and re-emerging viral infections. TG-23, a coronin 1, a tryptophan-aspartate-rich phagosomal protein-derived peptide, demonstrated inhibition of fusion between small unilamellar vesicles (SUVs) by modulating the membrane's physical properties. However, its inhibitory efficacy reduces with an increasing concentration of membrane cholesterol. The present work aims to develop a fusion inhibitor whose efficacy would be unaltered in the presence of membrane cholesterol. A stretch of the tryptophan-aspartic acid-containing peptide with a similar secondary structure and hydrophobicity profile of TG-23 from coronin 1 was synthesized, and its ability to inhibit SUV-SUV fusion with varying concentrations of membrane cholesterol was evaluated. Our results demonstrate that the GG-21 peptide inhibits fusion irrespective of the cholesterol content of the membrane. We have further evaluated the peptide-induced change in the membrane organization and dynamics utilizing arrays of steady-state and time-resolved fluorescence measurements and correlated these results with their effect on fusion. Interestingly, GG-21 displays inhibitory efficacy in a wide variety of lipid compositions despite having a secondary structure and physical properties similar to those of TG-23. Overall, our results advocate that the secondary structure and physical properties of the peptide may not be sufficient to predict its inhibitory efficacy.

Sensing membrane voltage by reorientation of dipolar transmembrane peptides

Membrane voltage plays a vital role in the behavior and functions of the lipid bilayer membrane. For instance, it regulates the exchange of molecules across the membrane through transmembrane proteins such as ion channels. In this paper, we study the membrane voltage-sensing mechanism, which entails the reorientation of α-helices with a change in the membrane voltage. We consider a helix having a large electrical macrodipole embedded in a lipid bilayer as a model system. We performed extensive molecular dynamics simulations to study the effect of variation of membrane voltage on the tilt angle of peptides and ascertain the optimal parameters for designing such a voltage-sensing peptide. A theoretical model for the system is also developed to investigate the interplay of competing effects of hydrophobic mismatch and dipole-electric field coupling on the tilt of the peptide and further explore the parameter space. This work opens the possibility for the design and fabrication of artificial dipolar membrane voltage-sensing elements for biomedical applications.

Double-Transmembrane Domain of SNAREs Decelerates the Fusion by Increasing the Protein-Lipid Mismatch

SNARE is the essential mediator of membrane fusion that highly relies on the molecular structure of SNAREs. For instance, the protein syntaxin-1 involved in neuronal SNAREs, has a single transmembrane domain (sTMD) leading to fast fusion, while the syntaxin 17 has a V-shape double TMDs (dTMDs), taking part in the autophagosome maturation. However, it is not clear how the TMD structure influences the fusion process. Here, we demonstrate that the dTMDs significantly reduce fusion rate compared with the sTMD by using an in vitro reconstitution system. Through theoretical analysis, we reveal that the V-shape dTMDs can significantly increase protein-lipid mismatch, thereby raising the energy barrier of the fusion, and that increasing the number of SNAREs can reduce the energy barrier or protein-lipid mismatch. This study provides a physicochemical mechanistic understanding of SNARE-regulated membrane fusion.

Length matters: Functional flip of the short TatA transmembrane helix

The twin arginine translocase (Tat) exports folded proteins across bacterial membranes. The putative pore-forming or membrane-weakening component (TatA_d in B. subtilis) is anchored to the lipid bilayer via an unusually short transmembrane α-helix (TMH), with less than 16 residues. Its tilt angle in different membranes was analyzed under hydrophobic mismatch conditions, using synchrotron radiation circular dichroism and solid-state NMR. Positive mismatch (introduced either by reconstitution in short-chain lipids or by extending the hydrophobic TMH length) increased the helix tilt of the TMH as expected. Negative mismatch (introduced either by reconstitution in long-chain lipids or by shortening the TMH), on the other hand, led to protein aggregation. These data suggest that the TMH of TatA is just about long enough for stable membrane insertion. At the same time, its short length is a crucial factor for successful translocation, as demonstrated here in native membrane vesicles using an in vitro translocation assay. Furthermore, when reconstituted in model membranes with negative spontaneous curvature, the TMH was found to be aligned parallel to the membrane surface. This intrinsic ability of TatA to flip out of the membrane core thus seems to play a key role in its membrane-destabilizing effect during Tat-dependent translocation.

Free-energy landscapes and insertion pathways for peptides in membrane environment

Free-energy landscapes for short peptides-specifically for variants of the pH low insertion peptide (pHLIP)-in the heterogeneous environment of a lipid bilayer or cell membrane are constructed, taking into account a set of dominant interactions and the conformational preferences of the peptide backbone. Our methodology interprets broken internal H-bonds along the backbone of a polypeptide as statistically interacting quasiparticles, activated from the helix reference state. The favored conformation depends on the local environment (ranging from polar to nonpolar), specifically on the availability of external H-bonds (with H_{2}O molecules or lipid headgroups) to replace internal H-bonds. The dominant side-chain contribution is accounted for by residue-specific transfer free energies between polar and nonpolar environments. The free-energy landscape is sensitive to the level of pH in the aqueous environment surrounding the membrane. For high pH, we identify pathways of descending free energy that suggest a coexistence of membrane-adsorbed peptides with peptides in solution. A drop in pH raises the degree of protonation of negatively charged residues and thus increases the hydrophobicity of peptide segments near the C terminus. For low pH, we identify insertion pathways between the membrane-adsorbed state and a stable trans-membrane state with the C terminus having crossed the membrane.

Structure and dynamics of the SARS-CoV-2 envelope protein monomer

Coronaviruses, especially severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), present an ongoing threat to human wellbeing. Consequently, elucidation of molecular determinants of their function and interaction with the host is an important task. Whereas some of the coronaviral proteins are extensively characterized, others remain understudied. Here, we use molecular dynamics simulations to analyze the structure and dynamics of the SARS-CoV-2 envelope (E) protein (a viroporin) in the monomeric form. The protein consists of the hydrophobic α-helical transmembrane domain (TMD) and amphiphilic α-helices H2 and H3, connected by flexible linkers. We show that TMD has a preferable orientation in the membrane, while H2 and H3 reside at the membrane surface. Orientation of H2 is strongly influenced by palmitoylation of cysteines Cys40, Cys43, and Cys44. Glycosylation of Asn66 affects the orientation of H3. We also observe that the monomeric E protein both generates and senses the membrane curvature, preferably localizing with the C-terminus at the convex regions of the membrane; the protein in the pentameric form displays these properties as well. Localization to curved regions may be favorable for assembly of the E protein oligomers, whereas induction of curvature may facilitate the budding of the viral particles. The presented results may be helpful for a better understanding of the function of the coronaviral E protein and viroporins in general, and for overcoming the ongoing SARS-CoV-2 pandemic.

Spiers Memorial Lecture: Analysis and de novo design of membrane-interactive peptides

Membrane-peptide interactions play critical roles in many cellular and organismic functions, including protection from infection, remodeling of membranes, signaling, and ion transport. Peptides interact with membranes in a variety of ways: some associate with membrane surfaces in either intrinsically disordered conformations or well-defined secondary structures. Peptides with sufficient hydrophobicity can also insert vertically as transmembrane monomers, and many associate further into membrane-spanning helical bundles. Indeed, some peptides progress through each of these stages in the process of forming oligomeric bundles. In each case, the structure of the peptide and the membrane represent a delicate balance between peptide-membrane and peptide-peptide interactions. We will review this literature from the perspective of several biologically important systems, including antimicrobial peptides and their mimics, α-synuclein, receptor tyrosine kinases, and ion channels. We also discuss the use of de novo design to construct models to test our understanding of the underlying principles and to provide useful leads for pharmaceutical intervention of diseases.

Illuminating Disorder Induced by Glu in a Stable Arg-Anchored Transmembrane Helix

Membrane proteins are vital for biological function and are complex to study. Even in model peptide-lipid systems, the combined influence or interaction of pairs of chemical groups still is not well understood. Disordered proteins, whether in solution or near lipid membranes, are an emerging paradigm for the initiation and control of biological function. The disorder can involve molecular orientation as well as molecular folding. This paper reports an astonishing induction of disorder when one Glu residue is introduced into a highly stable 23-residue transmembrane helix. The parent helix is anchored by a single Arg residue, tilted at a well-defined angle with respect to the DOPC bilayer normal and undergoes rapid cone precession. When Glu is introduced two residues away from Arg, with 200° (or 160°) radial separation, the helix properties change radically to exhibit a multiplicity of three or more disordered states. The helix characteristics have been monitored by deuterium (²H) NMR spectroscopy as functions of the pH and lipid bilayer composition. The disordered multistate behavior of the (Glu, Arg)-containing helix varies with the lipid bilayer thickness and pH. The results highlight a fundamental induction of protein multistate properties by a single Glu residue in a lipid membrane environment.

A photoswitchable helical peptide with light-controllable interface/transmembrane topology in lipidic membranes

The spontaneous insertion of helical transmembrane (TM) polypeptides into lipid bilayers is driven by three sequential equilibria: solution-to-membrane interface (MI) partition, unstructured-to-helical folding, and MI-to-TM helix insertion. A bottleneck for understanding these three steps is the lack of experimental approaches to perturb membrane-bound hydrophobic polypeptides out of equilibrium rapidly and reversibly. Here, we report on a 24-residues-long hydrophobic α-helical polypeptide, covalently coupled to an azobenzene photoswitch (KCALP-azo), which displays a light-controllable TM/MI equilibrium in hydrated lipid bilayers. FTIR spectroscopy reveals that trans KCALP-azo folds as a TM α-helix (TM topology). After trans-to-cis photoisomerization of the azobenzene moiety with UV light (reversed with blue light), the helical structure of KCALP-azo is maintained, but its helix tilt increased from 32 ± 5° to 79 ± 8°, indication of a reversible TM-to-MI transition. Further analysis indicates that this transition is incomplete, with cis KCALP-azo existing in a ∼90% TM and ∼10% MI mixture.

•

We present an α-helical transmembrane peptide modified with a molecular photoswitch

•

The peptide exhibits reversible photocontrol of its membrane topology

•

A fraction moves to the membrane interface with UV and inserts back with blue light

•

This system will be useful to address the molecular mechanism for membrane insertion

Lipid-Dependent Titration of Glutamic Acid at a Bilayer Membrane Interface

The ionization properties of protein side chains in lipid-bilayer membranes will differ from the canonical values of side chains exposed to an aqueous solution. While the propensities of positively charged side chains of His, Lys, and Arg to release a proton in lipid membranes have been rather well characterized, the propensity for a negatively charged Glu side chain to receive a proton and achieve the neutral state in a bilayer membrane has been less well characterized. Indeed, the ionization of the glutamic acid side chain has been predicted to depend on its depth of burial in a lipid membrane but has been difficult to verify experimentally. To address the issue, we incorporated an interfacial Glu residue at position 4 of a distinct 23-residue transmembrane helix and used ²H NMR to examine the helix properties as a function of pH. We observe that the helix tilt and azimuthal rotation vary little with pH, but the extent of helix unraveling near residues 3 and 4 changes as the Glu residue E4 titrates. Remarkably, the ²H quadrupolar splitting for the side chain of alanine A3 responds to pH with an apparent pK _a of 4.8 in 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) and 6.3 in 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC), but is unchanged up to pH 8.0 in 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) in the presence of residue E4. With bilayers composed of alkali-stable ether-linked lipids, the side chain of A3 responds to pH with an apparent pK _a of 11.0 in the ether analogue of DOPC. These results suggest that the depth dependence of Glu ionization in lipid-bilayer membranes may be steeper than previously predicted or envisioned.

Examination of pH dependency and orientation differences of membrane spanning alpha helices carrying a single or pair of buried histidine residues

We have employed the peptide framework of GWALP23 (acetyl-GGALWLALALALALALALWLAGA-amide) to examine the orientation, dynamics and pH dependence of peptides having buried single or pairs of histidine residues. When residue L8 is substituted to yield GWALP23-H⁸, acetyl-GGALWLAH⁸ALALALALALWLAGA-amide, the deuterium NMR spectra of ²H-labeled core alanine residues reveal a helix that occupies a single transmembrane orientation in DLPC, or in DMPC at low pH, yet shows multiple states at higher pH or in bilayers of DOPC. Moreover, a single histidine at position 8 or 16 in the GWALP23 framework is sensitive to pH. Titration points are observed near pH 3.5 for the deprotonation of H8 in lipid bilayers of DLPC or DMPC, and for H16 in DOPC. When residues L8 and L16 both are substituted to yield GWALP23-H^8,16, the ²H NMR spectra show, interestingly, no titration dependence from pH 2-8, yet bilayer thickness-dependent orientation differences. The helix with H8 and H16 is found to adopt a transmembrane orientation in thin bilayers of DLPC, a combination of transmembrane and surface orientations in DMPC, and then a complete transition to a surface bound orientation in the thicker DPoPC and DOPC lipid bilayers. In the surface orientations, alanine A7 no longer fits within the core helix. These results along with previous studies with different locations of histidine residues suggest that lipid hydrophobic thickness is a first determinant and pH a second determinant for the helical orientation, along with possible side-chain snorkeling, when the His residues are incorporated into the hydrophobic region of a lipid membrane-associated helix.

Molecular design of antimicrobial conjugated oligoelectrolytes with enhanced selectivity toward bacterial cells

A series of cationic conjugated oligoelectrolytes (COEs) was designed to understand how variations in molecular dimensions impact the relative activity against bacteria and mammalian cells. These COEs kept a consistent distyrylbenzene framework but differed in the length of linker between the core and the cationic site and the length of substitute on the quaternary ammonium functioned group. Their antimicrobial efficacy, mammalian cell cytotoxicity, hemolytic activity, and cell association were determined. We find that hydrophobicity is a factor that controls the degree of COE association to cells, but in vitro efficacy and cytotoxicity depend on more subtle structural features. COE2-3C-C4butyl was found to be the optimal structure with a minimum inhibitory concentration (MIC) of 4 μg mL⁻¹ against E. coli K12, low cytotoxicity against HepG2 cells and negligible hemolysis of red blood cells, even at 1024 μg mL⁻¹. A time-kill kinetics study of COE2-3C-C4butyl against E. coli K12 demonstrates bactericidal activity. These findings provide the first systematic investigation of how COEs may be modulated to achieve low mammalian cell cytotoxicity with the long-range perspective of finding candidates suitable for developing a broad-spectrum antimicrobial agent.

A series of cationic conjugated oligoelectrolytes (COEs) was designed to understand how variations in molecular dimensions impact the relative activity against bacteria and mammalian cells.

Computer simulations of protein-membrane systems

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

Comparing Interfacial Trp, Interfacial His and pH Dependence for the Anchoring of Tilted Transmembrane Helical Peptides

Charged and aromatic amino acid residues, being enriched toward the terminals of membrane-spanning helices in membrane proteins, help to stabilize particular transmembrane orientations. Among them, histidine is aromatic and can be positively charge at low pH. To enable investigations of the underlying protein-lipid interactions, we have examined the effects of single or pairs of interfacial histidine residues using the constructive low-dynamic GWALP23 (acetyl-GG²ALW⁵LALALALALALALW¹⁹LAG²²A-amide) peptide framework by incorporating individual or paired histidines at locations 2, 5, 19 or 22. Analysis of helix orientation by means of solid-state ²H NMR spectra of labeled alanine residues reveals marked differences with H^2,22 compared to W^2,22. Nevertheless, the properties of membrane-spanning H^2,22WALP23 helices show little pH dependence and are similar to those having Gly, Arg or Lys at positions 2 and 22. The presence of H5 or H19 influences the helix rotational preference but not the tilt magnitude. H5 affects the helical integrity, as residue 7 unwinds from the core helix; yet once again the helix orientation and dynamic properties show little sensitivity to pH. The overall results reveal that the detailed properties of transmembrane helices depend upon the precise locations of interfacial histidine residues.

Influence of interfacial tryptophan residues on an arginine-flanked transmembrane helix

The transmembrane helices of membrane proteins often are flanked by interfacial charged or aromatic residues that potentially help to anchor the membrane-spanning protein. For isolated single-span helices, the interfacial residues may be especially important for stabilizing particular tilted transmembrane orientations. The peptide RWALP23 (acetyl-GR²ALW(LA)₆LWLAR²²A-amide) has been employed to investigate the interplay between interfacial arginines and tryptophans. Here we replace the tryptophans of RWALP23 with A5 and A19, to investigate arginines alone with respect to helix fraying and orientation in varying lipid bilayers. Deuterated alanines incorporated into the central sequence allow the orientation and stability of the core helix to be assessed by means of solid -state ²H NMR in bilayers of DOPC, DMPC and DLPC. The helix tilt from the bilayer normal is found to increase slightly when R2 and R22 are present, and increases still further when the tryptophans W5 and W19 are replaced by alanines. The extent of helix dynamic averaging remains low in all cases. The preferred helix azimuthal rotation is essentially constant for all of the helices in each of the lipid membranes considered here. The alanines located outside of the core region of the peptide are sensitive to helical integrity. The new alanines, A5 and A19, therefore, provide new information about the length of the core helix and the onset of unraveling of the terminals. Residue A19 remains essentially on the central helix in each lipid membrane, while residues A3, A5 and A21 deviate from the core helix to an extent that depends on the membrane thickness. Differential unraveling of the two ends to expose peptide backbone groups for hydrogen bonding therefore acts together with specific interfacial side chains to stabilize a transmembrane helix.

Bilayer thickness determines the alignment of model polyproline helices in lipid membranes

Our understanding of protein folds relies fundamentally on the set of secondary structures found in the proteomes. Yet, there also exist intriguing structures and motifs that are underrepresented in natural biopolymeric systems. One example is the polyproline II helix, which is usually considered to have a polar character and therefore does not form membrane spanning sections of membrane proteins. In our work, we have introduced specially designed polyproline II helices into the hydrophobic membrane milieu and used 19F NMR to monitor the helix alignment in oriented lipid bilayers. Our results show that these artificial hydrophobic peptides can adopt several different alignment states. If the helix is shorter than the thickness of the hydrophobic core of the membrane, it is submerged into the bilayer with its long axis parallel to the membrane plane. The polyproline helix adopts a transmembrane alignment when its length exceeds the bilayer thickness. If the peptide length roughly matches the lipid thickness, a coexistence of both states is observed. We thus show that the lipid thickness plays a determining role in the occurrence of a transmembrane polyproline II helix. We also found that the adaptation of polyproline II helices to hydrophobic mismatch is in some notable aspects different from α-helices. Finally, our results prove that the polyproline II helix is a competent structure for the construction of transmembrane peptide segments, despite the fact that no such motif has ever been reported in natural systems.

Influence of Lipid Saturation, Hydrophobic Length and Cholesterol on Double-Arginine-Containing Helical Peptides in Bilayer Membranes

Membrane proteins are essential for many cell processes yet are more difficult to investigate than soluble proteins. Charged residues often contribute significantly to membrane protein function. Model peptides such as GWALP23 (acetyl-GGALW⁵ LAL⁸ LALALAL¹⁶ ALW¹⁹ LAGA-amide) can be used to characterize the influence of specific residues on transmembrane protein domains. We have substituted R8 and R16 in GWALP23 in place of L8 and L16, equidistant from the peptide center, and incorporated specific ² H-labeled alanine residues within the central sequence for detection by solid-state ² H NMR spectroscopy. The resulting pattern of [² H]Ala quadrupolar splitting (Δν_q ) magnitudes indicates the core helix for R^8,16 GWALP23 is significantly tilted to give a similar transmembrane orientation in thinner bilayers with either saturated C12:0 or C14:0 acyl chains (1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) or 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC)) or unsaturated C16:1 Δ9 cis acyl chains. In bilayers of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC; C18:1 Δ9 cis) multiple orientations are indicated, whereas in longer, unsaturated 1,2-dieicosenoyl-sn-glycero-3-phosphocholine (DEiPC; C20:1 Δ11 cis) bilayers, the R^8,16 GWALP23 helix adopts primarily a surface orientation. The inclusion of 10-20 mol % cholesterol in DOPC bilayers drives more of the R^8,16 GWALP23 helix population to the membrane surface, thereby allowing both charged arginines access to the interfacial lipid head groups. The results suggest that hydrophobic thickness and cholesterol content are more important than lipid saturation for the arginine peptide dynamics and helix orientation in lipid membranes.

Breaking the Backbone: Central Arginine Residues Induce Membrane Exit and Helix Distortions within a Dynamic Membrane Peptide

Transmembrane domains of membrane proteins sometimes contain conserved charged or ionizable residues which may be essential for protein function and regulation. This work examines the molecular interactions of single Arg residues within a highly dynamic transmembrane peptide helix. To this end, we have modified the GW^4,20ALP23 (acetyl-GGAW⁴(AL)₇AW²⁰AGA-amide) model peptide framework to incorporate Arg residues near the center of the peptide. Peptide helix formation, orientation and dynamics were analyzed by means of solid-state NMR spectroscopy to monitor specific ²H- or ¹⁵N-labeled residues. GW^4,20ALP23 itself adopts a tilted orientation within lipid bilayer membranes. Nevertheless, the GW^4,20ALP23 helix exhibits moderate to high dynamic averaging of NMR observables, such as ²H quadrupolar splittings or ¹⁵N-¹H dipolar couplings, due to competition between the interfacial Trp residues on opposing helix faces. Here we examine how the helix dynamics are impacted by the introduction of a single Arg residue at position 12 or 14. Residue R14 restricts the helix to low dynamic averaging and a well-defined tilt that varies inversely with the lipid bilayer thickness. To compensate for the dominance of R14, the competing Trp residues cause partial unwinding of the helix at the C-terminal. By contrast, R¹²GW^4,20ALP23 exits the DOPC bilayer to an interfacial surface-bound location. Interestingly, multiple orientations are exhibited by a single residue, Ala-9. Quadrupolar splittings generated by ²H-labeled residues A3, A5, A7, and A9 do not fit to the α-helical quadrupolar wave plot defined by residues A11, A13, A15, A17, A19, and A21. The discontinuity at residue A9 implicates a helical swivel distortion and an apparent 3₁₀-helix involving the N-terminal residues preceding A11. These molecular features suggest that, while arginine residues are prominent factors controlling transmembrane helix dynamics, the influence of interfacial tryptophan residues cannot be ignored.

Modulating Hinge Flexibility in the APP Transmembrane Domain Alters γ-Secretase Cleavage

Intramembrane cleavage of the β-amyloid precursor protein C99 substrate by γ-secretase is implicated in Alzheimer's disease pathogenesis. Biophysical data have suggested that the N-terminal part of the C99 transmembrane domain (TMD) is separated from the C-terminal cleavage domain by a di-glycine hinge. Because the flexibility of this hinge might be critical for γ-secretase cleavage, we mutated one of the glycine residues, G38, to a helix-stabilizing leucine and to a helix-distorting proline. Both mutants impaired γ-secretase cleavage and also altered its cleavage specificity. Circular dichroism, NMR, and backbone amide hydrogen/deuterium exchange measurements as well as molecular dynamics simulations showed that the mutations distinctly altered the intrinsic structural and dynamical properties of the substrate TMD. Although helix destabilization and/or unfolding was not observed at the initial ε-cleavage sites of C99, subtle changes in hinge flexibility were identified that substantially affected helix bending and twisting motions in the entire TMD. These resulted in altered orientation of the distal cleavage domain relative to the N-terminal TMD part. Our data suggest that both enhancing and reducing local helix flexibility of the di-glycine hinge may decrease the occurrence of enzyme-substrate complex conformations required for normal catalysis and that hinge mobility can thus be conducive for productive substrate-enzyme interactions.

Transmembrane Helix Integrity versus Fraying To Expose Hydrogen Bonds at a Membrane-Water Interface

Transmembrane helices dominate the landscape for many membrane proteins. Often flanked by interfacial aromatic residues, these transmembrane helices also contain loops and interhelix segments, which could help in stabilizing a transmembrane orientation. Using ²H nuclear magnetic resonance spectroscopy to monitor bilayer-incorporated model GWALP23 family peptides, we address systematically the issue of helix fraying in relation to the dynamics and orientation of highly similar individual transmembrane helices. We inserted aromatic (Phe, Trp, Tyr, and His) or non-aromatic residues (Ala and Gly) into positions 4 and 5 adjacent to a core transmembrane helix to examine the side-chain dependency of the transmembrane orientation, dynamics, and helix integrity (extent and location of unraveling). Incorporation of [²H]alanine labels enables one to assess the helicity of the core sequence and the peptide termini. For most of the helices, we observed substantial unwinding involving at least three residues at both ends. For the unique case of histidine at positions 4 and 5, an extended N-terminal unwinding was observed up to residue 7. For further investigation of the onset of fraying, we employed A^4,5GWALP23 with ²H labels at residues 4 and 5 and found that the number of terminal residues involved in the unwinding depends on bilayer thicknesses and helps to govern the helix dynamics. The combined results enable us to compare and contrast the extent of fraying for each related helix, as reflected by the deviation of experimental ²H quadrupolar splitting magnitudes of juxta-terminal alanines A3 and A21 from those represented by an ideal helix geometry.

Control of Transmembrane Helix Dynamics by Interfacial Tryptophan Residues

Transmembrane protein domains often contain interfacial aromatic residues, which may play a role in the insertion and stability of membrane helices. Residues such as Trp or Tyr, therefore, are often found situated at the lipid-water interface. We have examined the extent to which the precise radial locations of interfacial Trp residues may influence peptide helix orientation and dynamics. To address these questions, we have modified the GW5,19ALP23 (acetyl-GGALW5(LA)6LW19LAGA-[ethanol]amide) model peptide framework to relocate the Trp residues. Peptide orientation and dynamics were analyzed by means of solid-state nuclear magnetic resonance (NMR) spectroscopy to monitor specific 2H- and 15N-labeled residues. GW5,19ALP23 adopts a defined, tilted orientation within lipid bilayer membranes with minimal evidence of motional averaging of NMR observables, such as 2H quadrupolar or 15N-1H dipolar splittings. Here, we examine how peptide dynamics are impacted by relocating the interfacial Trp (W) residues on both ends and opposing faces of the helix, for example by a 100° rotation on the helical wheel for positions 4 and 20. In contrast to GW5,19ALP23, the modified GW4,20ALP23 helix experiences more extensive motional averaging of the NMR observables in several lipid bilayers of different thickness. Individual and combined Gaussian analyses of the 2H and 15N NMR signals confirm that the extent of dynamic averaging, particularly rotational "slippage" about the helix axis, is strongly coupled to the radial distribution of the interfacial Trp residues as well as the bilayer thickness. Additional 2H labels on alanines A3 and A21 reveal partial fraying of the helix ends. Even within the context of partial unwinding, the locations of particular Trp residues around the helix axis are prominent factors for determining transmembrane helix orientation and dynamics within the lipid membrane environment.

Informed Molecular Design of Conjugated Oligoelectrolytes To Increase Cell Affinity and Antimicrobial Activity

Membrane-intercalating conjugated oligoelectrolytes (COEs) are emerging as potential alternatives to conventional, yet increasingly ineffective, antibiotics. Three readily accessible COEs, belonging to an unreported series containing a stilbene core, namely D4, D6, and D8, were designed and synthesized so that the hydrophobicity increases with increasing side-chain length. Decreased aqueous solubility correlates with increased uptake by E. coli. The minimum inhibitory concentration (MIC) of D8 is 4 μg mL^-1 against both E. coli and E. faecalis, with an effective uptake of 72 %. In contrast, the MIC value of the shortest COE, D4, is 128 μg mL^-1 owing to the low cellular uptake of 3 %. These findings demonstrate the application of rational design to generate efficacious antimicrobial COEs that have potential as low-cost antimicrobial agents.

Transmembrane Polyproline Helix

The third most abundant polypeptide conformation in nature, the polyproline-II helix, is a polar, extended secondary structure with a local organization stabilized by intercarbonyl interactions within the peptide chain. Here we design a hydrophobic polyproline-II helical peptide based on an oligomeric octahydroindole-2-carboxylic acid scaffold and demonstrate its transmembrane alignment in model lipid bilayers by means of solid-state ¹⁹F NMR. As result, we provide a first example of a purely artificial transmembrane peptide with a structural organization that is not based on hydrogen-bonding.

New applications of solid-state NMR in structural biology

Various recent developments in solid-state nuclear magnetic resonance (ssNMR) spectroscopy have enabled an array of new insights regarding the structure, dynamics, and interactions of biomolecules. In the ever more integrated world of structural biology, ssNMR studies provide structural and dynamic information that is complementary to the data accessible by other means. ssNMR enables the study of samples lacking a crystalline lattice, featuring static as well as dynamic disorder, and does so independent of higher-order symmetry. The present study surveys recent applications of biomolecular ssNMR and examines how this technique is increasingly integrated with other structural biology techniques, such as (cryo) electron microscopy, solution-state NMR, and X-ray crystallography. Traditional ssNMR targets include lipid bilayer membranes and membrane proteins in a lipid bilayer environment. Another classic application has been in the area of protein misfolding and aggregation disorders, where ssNMR has provided essential structural data on oligomers and amyloid fibril aggregates. More recently, the application of ssNMR has expanded to a growing array of biological assemblies, ranging from non-amyloid protein aggregates, protein–protein complexes, viral capsids, and many others. Across these areas, multidimensional magic angle spinning (MAS) ssNMR has, in the last decade, revealed three-dimensional structures, including many that had been inaccessible by other structural biology techniques. Equally important insights in structural and molecular biology derive from the ability of MAS ssNMR to probe information beyond comprehensive protein structures, such as dynamics, solvent exposure, protein–protein interfaces, and substrate–enzyme interactions.

Helix formation and stability in membranes

In this article we review current understanding of basic principles for the folding of membrane proteins, focusing on the more abundant alpha-helical class. Membrane proteins, vital to many biological functions and implicated in numerous diseases, fold into their active conformations in the complex environment of the cell bilayer membrane. While many membrane proteins rely on the translocon and chaperone proteins to fold correctly, others can achieve their functional form in the absence of any translation apparatus or other aides. Nevertheless, the spontaneous folding process is not well understood at the molecular level. Recent findings suggest that helix fraying and loop formation may be important for overall structure, dynamics and regulation of function. Several types of membrane helices with ionizable amino acids change their topology with pH. Additionally we note that some peptides, including many that are rich in arginine, and a particular analogue of gramicidin, are able passively to translocate across cell membranes. The findings indicate that a final protein structure in a lipid-bilayer membrane is sequence-based, with lipids contributing to stability and regulation. While much progress has been made toward understanding the folding process for alpha-helical membrane proteins, it remains a work in progress. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

The role of hydrophobic matching on transmembrane helix packing in cells

Folding and packing of membrane proteins are highly influenced by the lipidic component of the membrane. Here, we explore how the hydrophobic mismatch (the difference between the hydrophobic span of a transmembrane protein region and the hydrophobic thickness of the lipid membrane around the protein) influences transmembrane helix packing in a cellular environment. Using a ToxRED assay in Escherichia coli and a Bimolecular Fluorescent Complementation approach in human-derived cells complemented by atomistic molecular dynamics simulations we analyzed the dimerization of Glycophorin A derived transmembrane segments. We concluded that, biological membranes can accommodate transmembrane homo-dimers with a wide range of hydrophobic lengths. Hydrophobic mismatch and its effects on dimerization are found to be considerably weaker than those previously observed in model membranes, or under in vitro conditions, indicating that biological membranes (particularly eukaryotic membranes) can adapt to structural deformations through compensatory mechanisms that emerge from their complex structure and composition to alleviate membrane stress. Results based on atomistic simulations support this view, as they revealed that Glycophorin A dimers remain stable, despite of poor hydrophobic match, using mechanisms based on dimer tilting or local membrane thickness perturbations. Furthermore, hetero-dimers with large length disparity between their monomers are also tolerated in cells, and the conclusions that one can draw are essentially similar to those found with homo-dimers. However, large differences between transmembrane helices length hinder the monomer/dimer equilibrium, confirming that, the hydrophobic mismatch has, nonetheless, biologically relevant effects on helix packing in vivo.

Critical Comparison of Biomembrane Force Fields: Protein-Lipid Interactions at the Membrane Interface

Molecular dynamics (MD) simulations offer the possibility to study biological processes at high spatial and temporal resolution often not reachable by experiments. Corresponding biomolecular force field parameters have been developed for a wide variety of molecules ranging from inorganic ligands and small organic molecules over proteins and lipids to nucleic acids. Force fields have typically been parametrized and validated on thermodynamic observables and structural characteristics of individual compounds, e.g. of soluble proteins or lipid bilayers. Less strictly, due to the added complexity and missing experimental data to compare to, force fields have hardly been tested on the properties of mixed systems, e.g. on protein-lipid systems. Their selection and combination for mixed systems is further complicated by the partially differing parametrization strategies. Additionally, the presence of other compounds in the system may shift the subtle balance of force field parameters. Here, we assessed the protein-lipid interactions as described in the four atomistic force fields GROMOS54a7, CHARMM36 and the two force field combinations Amber14sb/Slipids and Amber14sb/Lipid14. Four observables were compared, focusing on the membrane-water interface: the conservation of the secondary structure of transmembrane proteins, the positioning of transmembrane peptides relative to the lipid bilayer, the insertion depth of side chains of unfolded peptides absorbed at the membrane interface, and the ability to reproduce experimental insertion energies of Wimley-White peptides at the membrane interface. Significant differences between the force fields were observed that affect e.g. membrane insertion depths and tilting of transmembrane peptides.

Influence of glutamic acid residues and pH on the properties of transmembrane helices

Negatively charged side chains are important for the function of particular ion channels and certain other membrane proteins. To investigate the influence of single glutamic acid side chains on helices that span lipid-bilayer membranes, we have employed GWALP23 (acetyl-GGALW⁵LALALALALALALW¹⁹LAGA-amide) as a favorable host peptide framework. We substituted individual Leu residues with Glu residues (L12E or L14E or L16E) and incorporated specific ²H-labeled alanine residues within the core helical region or near the ends of the sequence. Solid-state ²H NMR spectra reveal little change for the core labels in GWALP23-E12, -E14 and -E16 over a pH range of 4 to 12.5, with the spectra being broader for samples in DOPC compared to DLPC bilayers. The spectra for samples with deuterium labels near the helix ends on alanines 3 and 21 show modest pH-dependent changes in the extent of unwinding of the helix terminals in DLPC and DOPC bilayers. The combined results indicate minor overall responses of these transmembrane helices to changes in pH, with the most buried residue E12 showing no pH dependence. While the Glu residues E14 and E16 may have high pK_a values in the lipid bilayer environment, it is also possible that a paucity of helix response is masking the pK_a values. Interestingly, when E16 is present, spectral changes at high pH report significant local unwinding of the core helix. Our results are consistent with the expectation that buried carboxyl groups aggressively hold their protons and/or waters of hydration.

Lipid bilayer thickness determines cholesterol's location in model membranes

Cholesterol is an essential biomolecule of animal cell membranes, and an important precursor for the biosynthesis of certain hormones and vitamins. It is also thought to play a key role in cell signaling processes associated with functional plasma membrane microdomains (domains enriched in cholesterol), commonly referred to as rafts. In all of these diverse biological phenomena, the transverse location of cholesterol in the membrane is almost certainly an important structural feature. Using a combination of neutron scattering and solid-state ²H NMR, we have determined the location and orientation of cholesterol in phosphatidylcholine (PC) model membranes having fatty acids of different lengths and degrees of unsaturation. The data establish that cholesterol reorients rapidly about the bilayer normal in all the membranes studied, but is tilted and forced to span the bilayer midplane in the very thin bilayers. The possibility that cholesterol lies flat in the middle of bilayers, including those made from PC lipids containing polyunsaturated fatty acids (PUFAs), is ruled out. These results support the notion that hydrophobic thickness is the primary determinant of cholesterol's location in membranes.

Revisiting the curvature-mediated interactions between proteins in biological membranes

Proteins embedded in soft biological membranes experience a long-range force mediated by elastic curvature deformations. The classical linearized Helfrich-Canham Hamiltonian based derivations reveal the nature of the force between a pair of proteins to be repulsive in the zero-temperature limit and the interaction potential is inversely proportional to the fourth power of the distance separating the inclusions. Such a result is the starting point to understand many-body interactions between proteins in biological membranes and the study of their clustering or, more broadly, self-organization. A key observation regarding this widely quoted result is that any two (mechanically rigid) proteins will experience an identical force. In other words, there is no specificity in the currently employed continuum models that purport to explain protein interactions. In this work we argue that each protein has a unique mechanical signature based on its interaction with the surrounding lipid bilayer membrane and cannot be treated as a non-specific rigid object. We modify the classical Helfrich-Canham theory of curvature elasticity to incorporate protein-membrane specificity, discuss the estimation of the new model parameters via atomistic simulations and re-evaluate the curvature-mediated force between proteins. We find that the incorporation of protein-specificity can reduce the interaction force by several orders of magnitude. Our result may provide at least one plausible reason behind why in some computational and experimental studies, a net attractive force between proteins is in evidence.

Influence of High pH and Cholesterol on Single Arginine-Containing Transmembrane Peptide Helices

An essential component of mammalian cells, cholesterol exerts significant influence on the physical properties of the cell membrane and in turn its constituents, including membrane proteins. Although sparse, polar amino acid residues are highly conserved in membrane proteins and play pivotal roles in determining specific structural and functional properties. To improve our understanding of particular polar residues in the membrane environment, we have examined two specific “guest” Arg residues within a well-defined and deuterium-labeled “host” framework provided by the transmembrane helical peptide GWALP23 (acetyl-GGALWLALALALALALALWLAGA-amide). Solid-state ²H nuclear magnetic resonance (NMR) spectra from aligned bilayer membrane samples effectively report changes in the host helix properties because of the incorporation of the guest residues. The focus of this work is two-pronged. First, GWALP23-R14 was examined over a pH range of 2–13 in 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) ester- or ether-linked bilayer membranes. Our findings indicate that the Arg guanidinium side chain remains charged over this entire range, in agreement with numerous molecular dynamics simulations. Second, GWALP23-R12 and GWALP23-R14 peptides were characterized in DOPC bilayers with varying cholesterol content. Our findings suggest that 10 or 20% cholesterol content has minimal impact on the orientation of the R14 peptide. Although the NMR signals are broader and weaker in the presence of 20% cholesterol, the deuterium quadrupolar splittings for [²H]Ala residues in GWALP23-R14 change very little. Conversely, cholesterol appears to modulate the multistate behavior of GWALP23-R12 and to favor a major interfacial state for the helix, bound at the bilayer surface. These results indicate a conditional sensitivity of a complex multistate transmembrane Arg-containing peptide helix to the presence of cholesterol.

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.

Helix perturbations in membrane proteins assist in inter-helical interactions and optimal helix positioning in the bilayer

Transmembrane (TM) helices in integral membrane proteins are primarily α-helical in structure. Here we analyze 1134 TM helices in 90 high resolution membrane proteins and find that apart from the widely prevalent α-helices, TM regions also contain stretches of 3₁₀ (3 to 8 residues) and π-helices (5 to 19 residues) with distinct sequence signatures. The various helix perturbations in TM regions comprise of helices with kinked geometry, as well as those with an interspersed 3₁₀/π-helical fragment and show high occurrence in a few membrane proteins. Proline is frequently present at sites of these perturbations, but it is neither a necessary nor a sufficient requirement. Helix perturbations are also conserved within a family of membrane proteins despite low sequence identity in the perturbed region. Furthermore, a perturbation influences the geometry of the TM helix, mediates inter-helical interactions within and across protein chains and avoids hydrophobic mismatch of the helix termini with the bilayer. An analysis of π-helices in the TM regions of the heme copper oxidase superfamily shows that interspersed π-helices can vary in length from 6 to 19 amino acids or be entirely absent, depending upon the protein function. The results presented here would be helpful for prediction of 3₁₀ and π-helices in TM regions and can assist the computational design of membrane proteins.

How Tolerant are Membrane Simulations with Mismatch in Area per Lipid between Leaflets?

Difficulties in estimating the correct number of lipids in each leaflet of complex bilayer membrane simulation systems make it inevitable to introduce a mismatch in lipid packing (i.e., area per lipid) and thus alter the lateral pressure of each leaflet. To investigate potential impacts of such mismatch on simulation results, we performed molecular dynamics simulations of saturated and monounsaturated lipid bilayers with and without gramicidin A or WALP23 at various mismatches by adjusting the number of lipids in the lower leaflet from no mismatch to a 25% reduction compared to that in the upper leaflet. All simulations were stable under the constant pressure barostat, but the mismatch induces asymmetric lipid packing between the leaflets, so that the upper leaflet becomes more ordered, and the lower leaflet becomes less ordered. The mismatch impacts on various bilayer properties are mild up to 5–10% mismatch, and bilayers with fully saturated chains appear to be more prone to these impacts than those with unsaturated tails. The nonvanishing leaflet surface tensions and the free energy derivatives with respect to the bilayer curvature indicate that the bilayer would be energetically unstable in the presence of mismatch. We propose a quantitative criterion for allowable mismatch based on the energetics derived from a continuum elastic model, which grows as a square root of the number of the lipids in the system. On the basis of this criterion, we infer that the area per lipid mismatch up to 5% would be tolerable in various membrane simulations of reasonable all-atom system sizes (40–160 lipids per leaflet).

Membrane Thinning and Thickening Induced by Membrane-Active Amphipathic Peptides

Membrane thinning has been discussed as a fundamental mechanism by which antimicrobial peptides can perturb cellular membranes. To understand which factors play a role in this process, we compared several amphipathic peptides with different structures, sizes and functions in their influence on the lipid bilayer thickness. PGLa and magainin 2 from X. laevis were studied as typical representatives of antimicrobial cationic amphipathic α-helices. A 1:1 mixture of these peptides, which is known to possess synergistically enhanced activity, allowed us to evaluate whether and how this synergistic interaction correlates with changes in membrane thickness. Other systems investigated here include the α-helical stress-response peptide TisB from E. coli (which forms membrane-spanning dimers), as well as gramicidin S from A. migulanus (a natural antibiotic), and BP100 (designer-made antimicrobial and cell penetrating peptide). The latter two are very short, with a circular β-pleated and a compact α-helical structure, respectively. Solid-state (2)H-NMR and grazing incidence small angle X-ray scattering (GISAXS) on oriented phospholipid bilayers were used as complementary techniques to access the hydrophobic thickness as well as the bilayer-bilayer repeat distance including the water layer in between. This way, we found that magainin 2, gramicidin S, and BP100 induced membrane thinning, as expected for amphiphilic peptides residing in the polar/apolar interface of the bilayer. PGLa, on the other hand, decreased the hydrophobic thickness only at very high peptide:lipid ratios, and did not change the bilayer-bilayer repeat distance. TisB even caused an increase in the hydrophobic thickness and repeat distance. When reconstituted as a mixture, PGLa and magainin 2 showed a moderate thinning effect which was less than that of magainin 2 alone, hence their synergistically enhanced activity does not seem to correlate with a modulation of membrane thickness. Overall, the absence of a typical thinning response in the case of PGLa, and the increase in the repeat distance and membrane thickening observed for TisB, demonstrate that the concept of peptide-induced membrane thinning cannot be generalized. Instead, these results suggest that different factors contribute to the resulting changes in membrane thickness, such as the peptide orientation in the bilayer, and/or bilayer adaptation to hydrophobic mismatch.

Hydrophobic Mismatch Drives the Interaction of E5 with the Transmembrane Segment of PDGF Receptor

The oncogenic E5 protein from bovine papillomavirus is a short (44 amino acids long) integral membrane protein that forms homodimers. It activates platelet-derived growth factor receptor (PDGFR) β in a ligand-independent manner by transmembrane helix-helix interactions. The nature of this recognition event remains elusive, as numerous mutations are tolerated in the E5 transmembrane segment, with the exception of one hydrogen-bonding residue. Here, we examined the conformation, stability, and alignment of the E5 protein in fluid lipid membranes of substantially varying bilayer thickness, in both the absence and presence of the PDGFR transmembrane segment. Quantitative synchrotron radiation circular dichroism analysis revealed a very long transmembrane helix for E5 of ∼26 amino acids. Oriented circular dichroism and solid-state ¹⁵N-NMR showed that the alignment and stability of this unusually long segment depend critically on the membrane thickness. When reconstituted alone in exceptionally thick DNPC lipid bilayers, the E5 helix was found to be inserted almost upright. In moderately thick bilayers (DErPC and DEiPC), it started to tilt and became slightly deformed, and finally it became aggregated in conventional DOPC, POPC, and DMPC membranes due to hydrophobic mismatch. On the other hand, when E5 was co-reconstituted with the transmembrane segment of PDGFR, it was able to tolerate even the most pronounced mismatch and was stabilized by binding to the receptor, which has the same hydrophobic length. As E5 is known to activate PDGFR within the thin membranes of the Golgi compartment, we suggest that the intrinsic hydrophobic mismatch of these two interaction partners drives them together. They seem to recognize each other by forming a closely packed bundle of mutually aligned transmembrane helices, which is further stabilized by a specific pair of hydrogen-bonding residues.

Effects of Transmembrane α-Helix Length and Concentration on Phase Behavior in Four-Component Lipid Mixtures: A Molecular Dynamics Study

We used coarse-grained molecular dynamics simulations to examine the effects of transmembrane α-helical WALP peptides on the behavior of four-component lipid mixtures. These mixtures contain a high-melting temperature (high-Tm) lipid, a nanodomain-inducing low-Tm lipid, a macrodomain-inducing low-Tm lipid and cholesterol to model the outer leaflet of cell plasma membranes. In a series of simulations, we incrementally replace the nanodomain-inducing low-Tm lipid by the macrodomain-inducing low-Tm lipid and measure how lipid and phase properties are altered by the addition of WALPs of different length. Regardless of the ratio of the two low-Tm lipids, shorter WALPs increase domain size and all WALPs increase domain alignment between the two leaflets. These effects are smallest for the longest WALP tested, and increase with increasing WALP concentration. Thus, our simulations explain the experimental observation that WALPs induce macroscopic domains in otherwise nanodomain-forming lipid-only mixtures (unpublished). Since the cell plasma membrane contains a large fraction of transmembrane proteins, these findings link the behavior of lipid-only model membranes in vitro to phase behavior in vivo.

Extending the Hydrophobic Mismatch Concept to Amphiphilic Membranolytic Peptides

A series of nine amphiphilic, pore-forming α-helical KIA peptides (KIAGKIA repeats) with lengths between 14 and 28 residues were studied by solid-state (15)N NMR to determine their alignment in oriented lipid bilayers. In a 2:1 mixture of 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) with its corresponding 1-myristoyl-2-hydroxy-sn-glycero-3-phosphocholine (lyso-MPC), which has a highly positive spontaneous curvature, the helix tilt angle was found to vary steadily with peptide length. The shortest peptide was aligned transmembrane and upright, while the longer ones successively became tilted away from the membrane normal. This behavior is in agreement with the hydrophobic matching concept, conceived so far only for hydrophobic helices. In 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine, with a negative spontaneous curvature, all KIA peptides remained flat on the bilayer surface, while the cylindrical DMPC lipids permitted a slight tilt. Peptide insertion thus depends critically on the intrinsic lipid curvature, and helix orientation is then fine-tuned by membrane thickness. A refined toroidal pore model is proposed.

Juxta-terminal Helix Unwinding as a Stabilizing Factor to Modulate the Dynamics of Transmembrane Helices

Transmembrane helices of integral membrane proteins often are flanked by interfacial aromatic residues that can serve as anchors to aid the stabilization of a tilted transmembrane orientation. Yet, physical factors that govern the orientation or dynamic averaging of individual transmembrane helices are not well understood and have not been adequately explained. Using solid-state (2) H NMR spectroscopy to examine lipid bilayer-incorporated model peptides of the GWALP23 (acetyl-GGALW(LA)6 LWLAGA-amide) family, we observed substantial unwinding at the terminals of several tilted helices spanning the membranes of DLPC, DMPC, or DOPC lipid bilayers. The fraying of helix ends might be vital for defining the dynamics and orientations of transmembrane helices in lipid bilayer membranes.

Alanine scan and (2)H NMR analysis of the membrane-active peptide BP100 point to a distinct carpet mechanism of action

The short membrane-active peptide BP100 [KKLFKKILKYL-NH2] is known as an effective antimicrobial and cell penetrating agent. For a functional alanine scan each of the 11 amino acids was replaced with deuterated Ala-d3, one at a time. MIC assays showed that a substitution of Lys did not affect the antimicrobial activity, but it decreased when a hydrophobic residue was replaced. In most cases, a reduction in hydrophobicity led to a decrease in hemolysis, and some peptide analogues had an improved therapeutic index. Circular dichroism showed that BP100 folds as an amphiphilic α-helix in a bilayer. Its alignment was determined from (2)H NMR in oriented membranes of different composition. The azimuthal rotation angle was the same under all conditions, but the average helix tilt angle and the dynamical behavior of the peptide varied in a systematic manner. In POPC/POPG bilayers, with a negative spontaneous curvature, the peptide was found to lie flat on the bilayer surface, and with little wobble. In DMPC/DMPG, with a positive spontaneous curvature, BP100 at higher concentrations became tilted obliquely into the membrane, with the uncharged C-terminus inserted more deeply into the lipid bilayer, experiencing significant fluctuations in tilt angle. In DMPC/DMPG/lyso-MPC, with a pronounced positive spontaneous curvature, the helix tilted even further and became even more mobile. The 11-mer BP100 is obviously too short to form transmembrane pores. We conclude that BP100 operates via a carpet mechanism, whereby the C-terminus gets inserted into the hydrophobic core of the bilayer, which leads to membrane perturbation and induces transient permeability.

Solid-State NMR-Restrained Ensemble Dynamics of a Membrane Protein in Explicit Membranes

Solid-state NMR has been used to determine the structures of membrane proteins in native-like lipid bilayer environments. Most structure calculations based on solid-state NMR observables are performed using simulated annealing with restrained molecular dynamics and an energy function, where all nonbonded interactions are represented by a single, purely repulsive term with no contributions from van der Waals attractive, electrostatic, or solvation energy. To our knowledge, this is the first application of an ensemble dynamics technique performed in explicit membranes that uses experimental solid-state NMR observables to obtain the refined structure of a membrane protein together with information about its dynamics and its interactions with lipids. Using the membrane-bound form of the fd coat protein as a model membrane protein and its experimental solid-state NMR data, we performed restrained ensemble dynamics simulations with different ensemble sizes in explicit membranes. For comparison, a molecular dynamics simulation of fd coat protein was also performed without any restraints. The average orientation of each protein helix is similar to a structure determined by traditional single-conformer approaches. However, their variations are limited in the resulting ensemble of structures with one or two replicas, as they are under the strong influence of solid-state NMR restraints. Although highly consistent with all solid-state NMR observables, the ensembles of more than two replicas show larger orientational variations similar to those observed in the molecular dynamics simulation without restraints. In particular, in these explicit membrane simulations, Lys(40), residing at the C-terminal side of the transmembrane helix, is observed to cause local membrane curvature. Therefore, compared to traditional single-conformer approaches in implicit environments, solid-state NMR restrained ensemble simulations in explicit membranes readily characterize not only protein dynamics but also protein-lipid interactions in detail.

Increased Microbial Butanol Tolerance by Exogenous Membrane Insertion Molecules

Butanol is an ideal biofuel, although poor titers lead to high recovery costs by distillation. Fluidization of microbial membranes by butanol is one of the major factors limiting titers in butanol-producing bioprocesses. Starting with the hypothesis that certain membrane insertion molecules would stabilize the lipid bilayer in the presence of butanol, we applied a combination of in vivo and in vitro techniques within an in silico framework to describe a new approach to achieve solvent tolerance in bacteria. Single-molecule tracking of a model supported bilayer showed that COE1-5C, a five-ringed oligo-polyphenylenevinylene conjugated oligoelectrolyte (COE), reduced the diffusion rate of phospholipids in a microbially derived lipid bilayer to a greater extent than three-ringed and four-ringed COEs. Furthermore, COE1-5C treatment increased the specific growth rate of E. coli K12 relative to a control at inhibitory butanol concentrations. Consequently, to confer butanol tolerance to microbes by exogenous means is complementary to genetic modification of strains in industrial bioprocesses, extends the physiological range of microbes to match favorable bioprocess conditions, and is amenable with complex and undefined microbial consortia for biobutanol production. Molecular dynamics simulations indicated that the π-conjugated aromatic backbone of COE1-5C likely acts as a hydrophobic tether for glycerophospholipid acyl chains by enhancing bilayer integrity in the presence of high butanol concentrations, which thereby counters membrane fluidization. COE1-5C-mitigated E. coli K12 membrane depolarization by butanol is consistent with the hypothesis that improved growth rates in the presence of butanol are a consequence of improved bilayer stability.

A simple and transferable all-atom/coarse-grained hybrid model to study membrane processes

We present an efficient all-atom/coarse-grained hybrid model and apply it to membrane processes. This model is an extension of the all-atom/ELBA model applied previously to processes in water. Here, we improve the efficiency of the model by implementing a multiple-time step integrator that allows the atoms and the coarse-grained beads to be propagated at different timesteps. Furthermore, we fine-tune the interaction between the atoms and the coarse-grained beads by computing the potential of mean force of amino acid side chain analogs along the membrane normal and comparing to atomistic simulations. The model was independently validated on the calculation of small-molecule partition coefficients. Finally, we apply the model to membrane peptides. We studied the tilt angle of the Walp23 and Kalp23 helices in two different model membranes and the stability of the glycophorin A dimer. The model is efficient, accurate, and straightforward to use, as it does not require any extra interaction particles, layers of atomistic solvent molecules or tabulated potentials, thus offering a novel, simple approach to study membrane processes.