Coarse-grained simulations of the membrane-active antimicrobial Peptide maculatin 1.1

Maculatin 1.1 (M1.1) is a membrane-active antimicrobial peptide (AMP) from an Australian tree frog that forms a kinked amphipathic α-helix in the presence of a lipid bilayer or bilayer-mimetic environment. To help elucidate its mechanism of membrane-lytic activity, we performed a total of ∼8 μs of coarse-grained molecular dynamics (CG-MD) simulations of M1.1 in the presence of zwitterionic phospholipid membranes. Several systems were simulated in which the peptide/lipid ratio was varied. At a low peptide/lipid ratio, M1.1 adopted a kinked, membrane-interfacial location, consistent with experiment. At higher peptide/lipid ratios, we observed spontaneous, cooperative membrane insertion of M1.1 peptide aggregates. The minimum size for formation of a transmembrane (TM) aggregate was just four peptides. The absence of a simple and well-defined central channel, along with the exclusion of lipid headgroups from the aggregates, suggests that a pore-like model is an unlikely explanation for the mechanism of membrane lysis by M1.1. We also performed an extended 1.25 μs simulation of the permeabilization of a complete liposome by multiple peptides. Consistent with the simpler bilayer simulations, formation of monomeric interfacial peptides and TM peptide clusters was observed. In contrast, major structural changes were observed in the vesicle membrane, implicating induced membrane curvature in the mechanism of active antimicrobial peptide lysis. This contrasted with the behavior of the nonpore-forming model peptide WALP23, which inserted into the vesicle to form extended clusters of TM α-helices with relatively little perturbation of bilayer properties.     

Hydrophobic organization of α-helix membrane bundle in bacteriorhodopsin

The hydrophobic organization of the intramembraneα-helical bundle in bacteriorhodopsin (BRh) was assessed based on a new approach to characterization of spatial hydrophobic properties of transmembrane (TM)α-helical peptides. The method employs two independent techniques: Monte Carlo simulations of nonpolar solvent around TM peptides and analysis of molecular hydrophobicity potential on their surfaces. The results obtained by the two methods agree with each other and permit precise hydrophobicity mapping of TM peptides. Superimposition of such data on the experimentally derived spatial model of the membrane moiety together with 2D maps of hydrophobic hydrophilic contacts provide considerable insight into the hydrophobic organization of BRh. The helix bundle is stabilized to a large extent by hydrophobic interactions between helices—neighbors in the sequence of BRh, by long-range interactions in helix pairs C-E, C-F, and C-G, and by nonpolar contracts between retinal and helices C, D, E, F. Unlike globular proteins, no polar contacts between residues distantly separated in the sequence of BRh were found in the bundle. One of the most striking results of this study is the finding that the hydrophobic organization of BRh is significantly different from those in bacterial photoreaction centers. Thus, TMα-helices in BRh expose their most nonpolar sides to the bilayer as well as to the neighboring helices and to the interior of the bundle. Some of them contact lipids with their relatively hydrophilic surfaces. No correlation was found between disposition of the most hydrophobic and the most variable sides of the TM helices.     

Transmembrane helix dynamics of bacterial chemoreceptors supports a piston model of signalling

Transmembrane α-helices play a key role in many receptors, transmitting a signal from one side to the other of the lipid bilayer membrane. Bacterial chemoreceptors are one of the best studied such systems, with a wealth of biophysical and mutational data indicating a key role for the TM2 helix in signalling. In particular, aromatic (Trp and Tyr) and basic (Arg) residues help to lock α-helices into a membrane. Mutants in TM2 of E. coli Tar and related chemoreceptors involving these residues implicate changes in helix location and/or orientation in signalling. We have investigated the detailed structural basis of this via high throughput coarse-grained molecular dynamics (CG-MD) of Tar TM2 and its mutants in lipid bilayers. We focus on the position (shift) and orientation (tilt, rotation) of TM2 relative to the bilayer and how these are perturbed in mutants relative to the wildtype. The simulations reveal a clear correlation between small (ca. 1.5 Å) shift in position of TM2 along the bilayer normal and downstream changes in signalling activity. Weaker correlations are seen with helix tilt, and little/none between signalling and helix twist. This analysis of relatively subtle changes was only possible because the high throughput simulation method allowed us to run large (n = 100) ensembles for substantial numbers of different helix sequences, amounting to ca. 2000 simulations in total. Overall, this analysis supports a swinging-piston model of transmembrane signalling by Tar and related chemoreceptors.     

Helix-helix interactions in membrane proteins: coarse-grained simulations of glycophorin a helix dimerization

Oligomerization of transmembrane (TM) helices is a key stage in the folding of membrane proteins. Glycophorin A (GpA) is a well-documented test system for this process. Coarse-grained molecular dynamics (CG-MD) allows us to simulate the self-assembly of TM helices into dimers, for both wild-type (WT) and mutant GpA sequences. For the WT sequences, dimers formed rapidly and remained stable in all simulations. The resultant dimers exhibited right-handed crossing and the same interhelix contacts as in NMR structures. Simulations of disruptive mutants revealed the dimers were less stable, with values of DeltaDelta G dimerization consistent with experimental data. The dimers of disruptive mutants were distorted relative to the WT and showed left-handed crossing of their helices. CG-MD can therefore be used to explore the interactions of TM helices, an important stage in the folding of membrane proteins. In particular, CG-MD has been shown to be sensitive enough to detect disruptions introduced by mutation. Future refinement of such models via atomistic simulations will enable a multiscale approach to predict the folding of membrane proteins.     

Differential repositioning of the second transmembrane helices from E. coli Tar and EnvZ upon moving the flanking aromatic residues

Aromatic tuning, i.e. repositioning aromatic residues found at the cytoplasmic end of transmembrane (TM) domains within bacterial receptors, has been previously shown to modulate signal output from the aspartate chemoreceptor (Tar) and the major osmosensor EnvZ of Escherichia coli. In the case of Tar, changes in signal output consistent with the vertical position of the native Trp-Tyr aromatic tandem within TM2 were observed. In contrast, within EnvZ, where a Trp-Leu-Phe aromatic triplet was repositioned, the surface that the triplet resided upon was the major determinant governing signal output. However, these studies failed to determine whether moving the aromatic residues was sufficient to physically reposition the TM helix within a membrane. Recent coarse-grained molecular dynamics (CG-MD) simulations predicted displacement of Tar TM2 upon moving the aromatic residues at the cytoplasmic end of the helix. Here, we demonstrate that repositioning the Trp-Tyr tandem within Tar TM2 displaces the C-terminal boundary of the helix relative to the membrane. In a similar analysis of EnvZ, an abrupt initial displacement of TM2 was observed but no subsequent movement was seen, suggesting that the vertical position of TM2 is not governed by the location of the Trp-Leu-Phe triplet. Our results also provide another set of experimental data, i.e. the resistance of EnvZ TM2 to being displaced upon aromatic tuning, which could be useful for subsequent refinement of the initial CG-MD simulations. Finally, we discuss the limitations of these methodologies, how moving flanking aromatic residues might impact steady-state signal output and the potential to employ aromatic tuning in other bacterial membrane-spanning receptors.     

Characterization of Membrane Protein Interactions by Isothermal Titration Calorimetry

Understanding the structure, folding, and interaction of membrane proteins requires experimental tools to quantify the association of transmembrane (TM) helices. Here, we introduce isothermal titration calorimetry (ITC) to measure integrin αIIbβ3 TM complex affinity, to study the consequences of helix–helix preorientation in lipid bilayers, and to examine protein-induced lipid reorganization. Phospholipid bicelles served as membrane mimics. The association of αIIbβ3 proceeded with a free energy change of − 4.61 ± 0.04 kcal/mol at bicelle conditions where the sampling of random helix–helix orientations leads to complex formation. At bicelle conditions that approach a true bilayer structure in effect, an entropy saving of > 1 kcal/mol was obtained from helix–helix preorientation. The magnitudes of enthalpy and entropy changes increased distinctly with bicelle dimensions, indicating long-range changes in bicelle lipid properties upon αIIbβ3 TM association. NMR spectroscopy confirmed ITC affinity measurements and revealed αIIbβ3 association and dissociation rates of 4500 ± 100 s^− 1 and 2.1 ± 0.1 s^− 1, respectively. Thus, ITC is able to provide comprehensive insight into the interaction of membrane proteins.     

Thermodynamic and kinetic stabilities of transmembrane helix bundles as revealed by single-pair FRET analysis: Effects of the number of membrane-spanning segments and cholesterol

The tertiary structures and conformational dynamics of transmembrane (TM) helical proteins are maintained by the interhelical interaction network in membranes, although it is complicated to analyze the underlying driving forces because the amino acid sequences can involve multiple and various types of interactions. To obtain insights into basal and common effects of the number of membrane-spanning segments and membrane cholesterol, we measured stabilities of helix bundles composed of simple TM helices (AALALAA)₃ (1TM) and (AALALAA)₃-G₅-(AALALAA)₃ (2TM). Association–dissociation dynamics for 1TM–1TM, 1TM–2TM, and 2TM–2TM pairs were monitored to compare stabilities of 2-, 3-, and 4-helical bundles, respectively, with single-pair fluorescence resonance energy transfer (sp-FRET) in liposome membranes. Both thermodynamic and kinetic stabilities of the helix bundles increased with a greater number of membrane-spanning segments in POPC. The presence of 30 mol% cholesterol strongly enhanced the formation of 1TM–1TM and 1TM–2TM bundles (~ ? 9 kJ mol^?1), whereas it only weakly stabilized the 2TM–2TM bundle (~ ? 3 kJ mol^?1). Fourier transform infrared-polarized attenuated total reflection (ATR-FTIR) spectroscopy revealed an ~30° tilt of the helix axis relative to bilayer normal for the 1TM–2TM pair in the presence of cholesterol, suggesting the formation of a tilted helix bundle to release high lateral pressure at the center of cholesterol-containing membranes. These results demonstrate that the number of membrane-spanning segments affects the stability and structure of the helix bundle, and their cholesterol-dependences. Such information is useful to understand the basics of folding and assembly of multispanning TM proteins.     

Revisiting Hydrophobic Mismatch with Free Energy Simulation Studies of Transmembrane Helix Tilt and Rotation

Protein-lipid interaction and bilayer regulation of membrane protein functions are largely controlled by the hydrophobic match between the transmembrane (TM) domain of membrane proteins and the surrounding lipid bilayer. To systematically characterize responses of a TM helix and lipid adaptations to a hydrophobic mismatch, we have performed a total of 5.8-μs umbrella sampling simulations and calculated the potentials of mean force (PMFs) as a function of TM helix tilt angle under various mismatch conditions. Single-pass TM peptides called WALPn (n = 16, 19, 23, and 27) were used in two lipid bilayers with different hydrophobic thicknesses to consider hydrophobic mismatch caused by either the TM length or the bilayer thickness. In addition, different flanking residues, such as alanine, lysine, and arginine, instead of tryptophan in WALP23 were used to examine their influence. The PMFs, their decomposition, and trajectory analysis demonstrate that 1), tilting of a single-pass TM helix is the major response to a hydrophobic mismatch; 2), TM helix tilting up to ∼10° is inherent due to the intrinsic entropic contribution arising from helix precession around the membrane normal even under a negative mismatch; 3), the favorable helix-lipid interaction provides additional driving forces for TM helix tilting under a positive mismatch; 4), the minimum-PMF tilt angle is generally located where there is the hydrophobic match and little lipid perturbation; 5), TM helix rotation is dependent on the specific helix-lipid interaction; and 6), anchoring residues at the hydrophilic/hydrophobic interface can be an important determinant of TM helix orientation.     

Effect of sequence hydrophobicity and bilayer width upon the minimum length required for the formation of transmembrane helices in membranes

The minimum hydrophobic length necessary to form a transmembrane (TM) helix in membranes was investigated using model membrane-inserted hydrophobic helices. The fluorescence of a Trp at the center of the sequence and its sensitivity to quenching were used to ascertain helix position within the membrane. Peptides with hydrophobic cores composed of poly(Leu) were compared to sequences containing a poly 1:1 Leu:Ala core (which have a hydrophobicity typical of natural TM helices). Studies varying bilayer width revealed that the poly(Leu) core peptides predominately formed a TM state when the bilayer width exceeded hydrophobic sequence length by (i.e. when negative mismatch was) up to ∼ 11-12 Å (e.g. the case of a 11-12 residue hydrophobic sequence in bilayers with a biologically relevant width, i.e. dioleoylphosphatidylcholine (DOPC) bilayers), while poly(LeuAla) core peptides formed predominantly TM state with negative mismatch of up to 9 Å (a 13 residue hydrophobic sequence in DOPC bilayers). This indicates that minimum length necessary to form a predominating amount of a TM state (minimum TM length) is only modestly hydrophobicity-dependent for the sequences studied here, and a formula that defines the minimum TM length as a function of hydrophobicity for moderately-to-highly hydrophobic sequences was derived. The minimum length able to form a stable TM helix for alternating LeuAla sequences, and that for sequences with a Leu block followed by an Ala block, was similar, suggesting that a hydrophobicity gradient along the sequence may not be an important factor in TM stability. TM stability was also similar for sequences flanked by different charged ionizable residues (Lys, His, Asp). However, ionizable flanking residues destabilized the TM configuration much more when charged than when uncharged. The ability of short hydrophobic sequences to form TM helices in membranes in the presence of substantial negative mismatch implies that lipid bilayers have a considerable ability to adjust to negative mismatch, and that short TM helices may be more common than generally believed. Factors that modulate the ability of bilayers to adjust to mismatch may strongly affect the configuration of short hydrophobic helices.     

Coarse-grained molecular dynamics simulations of membrane proteins and peptides

Molecular dynamics (MD) simulations provide a valuable approach to the dynamics, structure, and stability of membrane-protein systems. Coarse-grained (CG) models, in which small groups of atoms are treated as single particles, enable extended (>100 ns) timescales to be addressed. In this study, we explore how CG-MD methods that have been developed for detergents and lipids may be extended to membrane proteins. In particular, CG-MD simulations of a number of membrane peptides and proteins are used to characterize their interactions with lipid bilayers. CG-MD is used to simulate the insertion of synthetic model membrane peptides (WALPs and LS3) into a lipid (PC) bilayer. WALP peptides insert in a transmembrane orientation, whilst the LS3 peptide adopts an interfacial location, both in agreement with experimental biophysical data. This approach is extended to a transmembrane fragment of the Vpu protein from HIV-1, and to the coat protein from fd phage. Again, simulated protein/membrane interactions are in good agreement with solid state NMR data for these proteins. CG-MD has also been applied to an M3-M4 fragment from the CFTR protein. Simulations of CFTR M3-M4 in a detergent micelle reveal formation of an alpha-helical hairpin, consistent with a variety of biophysical data. In an I231D mutant, the M3-M4 hairpin is additionally stabilized via an inter-helix Q207/D231 interaction. Finally, CG-MD simulations are extended to a more complex membrane protein, the bacterial sugar transporter LacY. Comparison of a 200 ns CG-MD simulation of LacY in a DPPC bilayer with a 50 ns atomistic simulation of the same protein in a DMPC bilayer shows that the two methods yield comparable predictions of lipid-protein interactions. Taken together, these results demonstrate the utility of CG-MD simulations for studies of membrane/protein interactions.     

Sequence determinants of a transmembrane proton channel: an inverse relationship between stability and function

The driving forces behind the folding processes of integral membrane proteins after insertion into the bilayer, is currently under debate. The M2 protein from the influenza A virus is an ideal system to study lateral association of transmembrane helices. Its proton selective channel is essential for virus functioning and a target of the drug amantadine. A 25 residue transmembrane fragment of M2, M2TM, forms a four-helix bundle in vivo and in various detergents and phospholipid bilayers. Presented here are the energetic consequences for mutations made to the helix/helix interfaces of the M2TM tetramer. Analytical ultracentrifugation has been used to determine the effect of ten single-site mutations, to either alanine or phenylalanine, on the oligomeric state and the free energy of M2TM in the absence and the presence of amantadine. It was expected that many of these mutations would perturb the M2TM stability and tetrameric integrity. Interestingly, none of the mutations destabilize tetramerization. This finding suggests that M2 sacrifices stability to preserve its functions, which require rapid and specific interchange between distinct conformations involved in gating and proton conduction. Mutations might therefore restrict the full range of conformations by stabilizing a given native or non-native conformational state. In order to assess one specific conformation of the tetramer, we measured the binding of amantadine to the resting state of the channel, and examined the overall free energy of assembly of the amantadine bound tetramer. All of the mutations destabilized amantadine binding or were isoenergetic. We also find that large to small residue changes destabilize the amantadine bound tetramer whereas mutations to side-chains of similar volume stabilize this conformation. A structural model of the amantadine bound state of M2TM was generated using a novel protocol that optimizes a structure for an ensemble of neutral and disruptive mutations. The model structure is consistent with the mutational data.     

An alamethicin channel in a lipid bilayer: molecular dynamics simulations

We present the results of 2-ns molecular dynamics (MD) simulations of a hexameric bundle of Alm helices in a 1-palmitoyl-2-oleoylphosphatidylcholine bilayer. These simulations explore the dynamic properties of a model of a helix bundle channel in a complete phospholipid bilayer in an aqueous environment. We explore the stability and conformational dynamics of the bundle in a phospholipid bilayer. We also investigate the effect on bundle stability of the ionization state of the ring of Glu18 side chains. If all of the Glu18 side chains are ionised, the bundle is unstable; if none of the Glu18 side chains are ionized, the bundle is stable. pKA calculations suggest that either zero or one ionized Glu18 is present at neutral pH, correlating with the stable form of the helix bundle. The structural and dynamic properties of water in this model channel were examined. As in earlier in vacuo simulations (Breed et al., 1996 .Biophys. J. 70:1643-1661), the dipole moments of water molecules within the pore were aligned antiparallel to the helix dipoles. This contributes to the stability of the helix bundle.     

Conformational Changes in Talin on Binding to Anionic Phospholipid Membranes Facilitate Signaling by Integrin Transmembrane Helices

Integrins are heterodimeric (αβ) cell surface receptors that are activated to a high affinity state by the formation of a complex involving the α/β integrin transmembrane helix dimer, the head domain of talin (a cytoplasmic protein that links integrins to actin), and the membrane. The talin head domain contains four sub-domains (F0, F1, F2 and F3) with a long cationic loop inserted in the F1 domain. Here, we model the binding and interactions of the complete talin head domain with a phospholipid bilayer, using multiscale molecular dynamics simulations. The role of the inserted F1 loop, which is missing from the crystal structure of the talin head, PDB:3IVF, is explored. The results show that the talin head domain binds to the membrane predominantly via cationic regions on the F2 and F3 subdomains and the F1 loop. Upon binding, the intact talin head adopts a novel V-shaped conformation which optimizes its interactions with the membrane. Simulations of the complex of talin with the integrin α/β TM helix dimer in a membrane, show how this complex promotes a rearrangement, and eventual dissociation of, the integrin α and β transmembrane helices. A model for the talin-mediated integrin activation is proposed which describes how the mutual interplay of interactions between transmembrane helices, the cytoplasmic talin protein, and the lipid bilayer promotes integrin inside-out activation.     

Effect of lipid composition on the topography of membrane-associated hydrophobic helices: stabilization of transmembrane topography by anionic lipids

To investigate the effect of lipid structure upon the membrane topography of hydrophobic helices, the behavior of hydrophobic peptides was studied in model membrane vesicles. To define topography, fluorescence and fluorescence quenching methods were used to determine the location of a Trp at the center of the hydrophobic sequence. For peptides with cationic residues flanking the hydrophobic sequence, the stability of the transmembrane (TM) configuration (relative to a membrane-bound non-TM state) increased as a function of lipid composition on the order: 1:1 (mol:mol) 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC):1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine ∼ 6:4 POPC:cholesterol < POPC ∼ dioleoylphosphatidylcholine (DOPC) < 1,2-dioleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] sodium salt (DOPG) ≤ 1,2-dioleoyl-sn-glycero-3-[phospho-l-serine] sodium salt (DOPS), indicating that the anionic lipids DOPG and DOPS most strongly stabilized the TM configuration. TM stabilization was near maximal at 20-30 mol% anionic lipid, which are physiologically relevant values. TM stabilization by anionic lipid was observed for hydrophobic sequences with a diverse set of sequences (including polyAla), diverse lengths (from 12 to 22 residues), and various cationic flanking residues (H, R, or K), but not when the flanking residues were uncharged. TM stabilization by anionic lipid was also dependent on the number of cationic residues flanking the hydrophobic sequence, but was still significant with only one cationic residue flanking each end of the peptide. These observations are consistent with TM-stabilizing effects being electrostatic in origin. However, Trp located more deeply in DOPS vesicles relative to DOPG vesicles, and peptides in DOPS vesicles showed increased helix formation relative to DOPG and all other lipid compositions. These observations fit a model in which DOPS anchors flanking residues near the membrane surface more strongly than does DOPG and/or increases the stability of the TM state to a greater degree than DOPG. We conclude that anionic lipids can have significant and headgroup structure-specific effects upon membrane protein topography.     

Defining the transmembrane helix of M2 protein from influenza A by molecular dynamics simulations in a lipid bilayer

Integral membrane proteins containing at least one transmembrane (TM) alpha-helix are believed to account for between 20% and 30% of most genomes. There are several algorithms that accurately predict the number and position of TM helices within a membrane protein sequence. However, these methods tend to disagree over the beginning and end residues of TM helices, posing problems for subsequent modeling and simulation studies. Molecular dynamics (MD) simulations in an explicit lipid and water environment are used to help define the TM helix of the M2 protein from influenza A virus. Based on a comparison of the results of five different secondary structure prediction algorithms, three different helix lengths (an 18mer, a 26mer, and a 34mer) were simulated. Each simulation system contained 127 POPC molecules plus approximately 3500-4700 waters, giving a total of approximately 18,000-21,000 atoms. Two simulations, each of 2 ns duration, were run for the 18mer and 26mer, and five separate simulations were run for the 34mer, using different starting models generated by restrained in vacuo MD simulations. The total simulation time amounted to 11 ns. Analysis of the time-dependent secondary structure of the TM segments was used to define the regions that adopted a stable alpha-helical conformation throughout the simulation. This analysis indicates a core TM region of approximately 20 residues (from residue 22 to residue 43) that remained in an alpha-helical conformation. Analysis of atomic density profiles suggested that the 18mer helix revealed a local perturbation of the lipid bilayer. Polar side chains on either side of this region form relatively long-lived H-bonds to lipid headgroups and water molecules.     

Membrane insertion of a voltage sensor helix

Most membrane proteins contain a transmembrane (TM) domain made up of a bundle of lipid-bilayer-spanning α-helices. TM α-helices are generally composed of a core of largely hydrophobic amino acids, with basic and aromatic amino acids at each end of the helix forming interactions with the lipid headgroups and water. In contrast, the S4 helix of ion channel voltage sensor (VS) domains contains four or five basic (largely arginine) side chains along its length and yet adopts a TM orientation as part of an independently stable VS domain. Multiscale molecular dynamics simulations are used to explore how a charged TM S4 α-helix may be stabilized in a lipid bilayer, which is of relevance in the context of mechanisms of translocon-mediated insertion of S4. Free-energy profiles for insertion of the S4 helix into a phospholipid bilayer suggest that it is thermodynamically favorable for S4 to insert from water to the center of the membrane, where the helix adopts a TM orientation. This is consistent with crystal structures of Kv channels, biophysical studies of isolated VS domains in lipid bilayers, and studies of translocon-mediated S4 helix insertion. Decomposition of the free-energy profiles reveals the underlying physical basis for TM stability, whereby the preference of the hydrophobic residues of S4 to enter the bilayer dominates over the free-energy penalty for inserting charged residues, accompanied by local distortion of the bilayer and penetration of waters. We show that the unique combination of charged and hydrophobic residues in S4 allows it to insert stably into the membrane.     

Membrane tension, lipid adaptation, conformational changes, and energetics in MscL gating

This study aims to explore gating mechanisms of mechanosensitive channels in terms of membrane tension, membrane adaptation, protein conformation, and energetics. The large conductance mechanosensitive channel from Mycobacterium tuberculosis (Tb-MscL) is used as a model system; Tb-MscL acts as a safety valve by releasing small osmolytes through the channel opening under extreme hypoosmotic conditions. Based on the assumption that the channel gating involves tilting of the transmembrane (TM) helices, we have performed free energy simulations of Tb-MscL as a function of TM helix tilt angle in a dimyristoylphosphatidylcholine bilayer. Based on the change in system dimensions, TM helix tilting is shown to be essentially equivalent to applying an excess surface tension to the membrane, causing channel expansion, lipid adaptation, and membrane thinning. Such equivalence is further corroborated by the observation that the free energy cost of Tb-MscL channel expansion is comparable to the work done by the excess surface tension. Tb-MscL TM helix tilting results in an expanded water-conducting channel of an outer dimension similar to the proposed fully open MscL structure. The free energy decomposition indicates a possible expansion mechanism in which tilting and expanding of TM2 facilitates the iris-like motion of TM1, producing an expanded Tb-MscL.     

Evidence that the translocon may function as a hydropathy partitioning filter

Developing a greater understanding of the function of the translocon-and the source of its selectivity for transmembrane helix insertion-are important steps toward deciphering the role of disease-causing mutations in membrane regions. To address these phenomena, we have prepared a library of helix-loop-helix (“hairpin”) constructs derived from helices 3 and 4 of the first membrane domain of CFTR, in which position 232 was mutated individually to each of the 20 commonly-occurring amino acids. Using retention times on a reverse phase-HPLC C18 column to mimic the process of hairpin partitioning, we have quantitatively determined a hydropathy scale in the context of a bona fide membrane protein fragment that correlates to an in vivo hydropathy scale with r = −0.78—a value that rises to r = −0.92 when Asp and Glu are excluded due to protonation effects. Our results provide evidence that the translocon may act as a facilitator in the insertion selection process, effectively allowing the bilayer to “decide” through favorable non-polar solvation whether or not to allow a translocating helix to enter the membrane.