This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Stockner, T. Articles by Tieleman, D. P. Search for Related Content PubMed PubMed Citation Articles by Stockner, T. Articles by Tieleman, D. P.

Direct Simulation of Transmembrane Helix Association: Role of Asparagines

Department of Biological Sciences, University of Calgary, Calgary, Alberta, Canada

Correspondence: Address reprint requests to D. Peter Tieleman, Dept. of Biological Sciences, University of Calgary, 2500 University Dr., NW Calgary, Alberta T2N 4N1 Canada. Tel.: 403-220-2996; Fax: 403-298-9311; E-mail: tieleman{at}ucalgary.ca.

	ABSTRACT

TOP ABSTRACT INTRODUCTION METHODS RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
The forces contributing to the association of transmembrane helices in folded membrane proteins have received considerable attention recently. In this study we investigate the importance of hydrogen bonding by studying the effect of a single Asn residue in the center of an otherwise hydrophobic transmembrane peptide using computer simulations. We use the model peptide MS1 which has been derived from the leucine zipper coiled-coil dimer of the transcription factor peptide GCN4-P1. We follow the trajectory of 36 initially monomeric MS1 transmembrane helical peptides in a membrane-mimicking octane layer as they associate into larger structures. These peptides predominately form dimers. The interaction between the polar asparagine residues, capable of simultaneously being a hydrogen-bond donor and acceptor, contributes strongly to the stability of associated helices. Only dimers with interhelical hydrogen bonds form stable structures, whereas aggregates without any hydrogen-bonding interactions form very transient structures. We examine the hydrogen-bonding patterns and find that there are two forms of dimer, one with symmetric hydrogen bonds and one with asymmetric hydrogen bonds. Based on the structures in our simulation we propose a model with a monomer symmetric dimer asymmetric dimer trimer equilibrium.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION METHODS RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
Membrane proteins probably represent more than a quarter of all proteins encoded in genomes (Liu and Rost, 2001). They play an important role in cell biology and physiology, carrying out many essential cell functions. The forces stabilizing membrane protein structures have received considerable attention, although detailed structural information is difficult to obtain. A two-state model for membrane protein folding proposes that folding occurs in two separate steps (Popot and Engelman, 1990). The first step involves insertion of the transmembrane helices; the second involves the formation of specific interactions between these helices, forming tightly packed native structures. The association of TM helices can occur as either an inter- or intramolecular process. For example, bacteriorhodopsin (Kataoka et al., 1992) has been shown to fold into a functional protein by refolding from peptides representing its transmembrane helices, and many membrane proteins consist of several subunits which must associate to form the fully functional protein. Known motifs of closely packed helices in membranes are e.g. the GxxxG or the heptad repeat motif. In the GxxxG packing motif (Rees et al., 1989; Fleming and Engelman, 2001) the close approach of the backbones at the glycine positions facilitates weak hydrogen bonds between -carbon hydrogens and carbonyl groups stabilizing the dimer structure. In contrast, the soluble GCN4-P1 peptide has been found to form a left-handed coiled-coil structure (Crick, 1953; O'Shea et al., 1991) associating with a heptad repeat motif (Bowie, 1997; Walther et al., 1998) as shown in Fig. 1. In the heptad repeat motif the number of residues per turn is reduced to 3.5 compared to the normal 3.6 residues per turn, thereby improving the packing of the coiled-coil structures and repeating the packing motif each seven residues. The side chains of each residue of the dimer interface protrude into cavities formed by four side chains of the neighboring helices in a regular manner termed knobs into holes or leucine zipper packing (Langosch and Heringa, 1998; Zhou et al., 2001).

View larger version (24K): [in this window] [in a new window]   FIGURE 1  Heptad repeat of the coiled-coil structures. The dimer interface consists of residues in the a and d positions. The polar Asn residue of MS1 occupies an a position in the dimer interface. The sequence of the MS1 peptide is shown together with its heptad repeat positions.

Polar residues such as Asn and Gln are not equally distributed over the surface of known TM helices, but exhibit a strong tendency to be buried in the interface in the middle of the bilayer (Lear et al., 2003). If not buried within the protein structure these residues could potentially lead to nonspecific interactions between membrane proteins leading to misfolding and misfunction (Bowie, 2000; Zhou et al., 2000). There are several examples of nonpolar to hydrogen bonding residues leading to disease including the single mutation (V664E) in the neu/erb-2 protooncogene (Bargmann et al., 1986; Smith et al., 1996) leading to constitutive activation of the encoded tyrosine kinase receptor, or the V232D mutation in the TM4 helix of cystic fibrosis transmembrane conductance regulator (Partridge et al., 2002) altering the structure and function of the mature protein. Two Asn mutants (M701N and G708N) of the transmembrane helix ß3 of integrin IIbß3 have been identified to drive the association by homooligomerization (Li et al., 2003) activating the mutant integrin to constitutively bind fibrinogen.

In this study we are investigating the association behavior of single-folded transmembrane helices into higher aggregates and the importance of polar residues in the center of the membrane for specificity and stability using the model peptide MS1, designed by Choma et al. (2000). MS1 is based on a hydrophobic version of the leucine-zipper GCN4-P1 peptide, a 32-residue long homodimeric coiled-coil peptide from the yeast transcription factor GCN4 (O'Shea et al., 1991). The soluble parent peptide GCN4-P1 forms coiled-coil dimers, whereas mutants have been shown to exist in different association states (monomers, dimers, and trimers) (O'Shea et al., 1991; Harbury et al., 1993; Gonzalez et al., 1996a,b; Lino et al., 1996). The membrane soluble MS1 peptide conserves the seven residue heptad repeat (Langosch and Heringa, 1998; Gurezka et al., 1999) of GCN4-P1, with Leu residues occupying the "d" positions, and Val residues at all but one "a" position. The hydrophobicity of the interface is interrupted at the central "a" position by an Asn residue that is capable of forming hydrogen bonds across the dimer interface. MS1 associates in micelles in a reversible monomer-dimer-trimer equilibrium. This equilibrium depends on the type of micelle (Gratkowski et al., 2002) with a larger degree of trimers in C14-betaine and SDS micelles than in DPC micelles.

We are using MD simulations to investigate the association behavior of MS1 by following the trajectory of 36 peptides in a model membrane system over time. The dominant association state present in our simulation is the dimeric form, stabilized by interchain hydrogen bonds between the carboxamido groups of the Asn side chains. These dimer structures are on average left-handed coiled-coil structures indicating that the membrane soluble MS1 peptide assumes indeed a coiled-coil structure comparable to GCN4-P1 in water.

	METHODS

TOP ABSTRACT INTRODUCTION METHODS RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
Simulation setup
The starting structure was produced by inserting 36 identical, ideal -helices built from the sequence of MS1 (BQLLI AVLLL IAVNL ILLIA VARLR YLVG, B=ß-Ala) on a regular grid into a 3.0 nm thick layer of octane. The entire system was then solvated, and sodium and chloride ions were added to give a charge neutral system with a salt concentration of 250 mM. The final simulation box measured 17.0 x 17.0 x 7.5 nm and contained 3024 octane molecules, 37,764 water molecules, 180 Na⁺, and 288 Cl^– ions.

Molecular dynamics simulations were carried out using the GROMACS 3.0 MD package (Berendsen et al., 1995; Lindahl et al., 2001) applying periodic boundary conditions. The peptide and octane were represented using the GROMOS96 43a2 force field (Scott et al., 1999), whereas water was represented using the SPC model (Berendsen et al., 1981). The simulation was carried out for 45 ns with the temperature maintained at 300 K using a Berendsen (_T = 0.1 ps) thermostat (Berendsen et al., 1984), coupling the protein, the octane, and the water/ions separately. The area of the octane/water interface was held fixed and the pressure in the direction normal to the interface was maintained at 1 bar using the weak coupling algorithm (Berendsen et al., 1984) with a coupling constant of 1.0 ps and a compressibility of 4.6 x 10^–5 bar^–1. The electrostatic interactions were evaluated using the smooth particle mesh Ewald method (Darden et al., 1993; Essmann et al., 1995), with a real-space cutoff of 0.9 nm. The long-range electrostatic interactions were calculated with fourth-order B-spline interpolation and a Fourier spacing of 0.16 nm. The Lennard-Jones interactions were evaluated using a twin-range cutoff (0.9 and 1.4 nm) with the neighbor list updated every five steps. All bonds in the peptides and octane were constrained using LINCS (Hess et al., 1997). The bonds and angles of water were constrained using the SETTLE algorithm (Miyamoto and Kollman, 1992). Additionally, the hydrogen atoms of the peptides were treated as dummy atoms (Feenstra et al., 1999) allowing a time step of 5 fs.

Model membrane system
The octane layer used in our simulation models the aliphatic region of the lipids in the bilayer as it has approximately the same thickness, whereas the headgroup region is not included in our model. The use of octane as a membrane-mimicking environment enables us to investigate the association process and observe the formation of dimers, trimers, and higher-order aggregates. The lateral diffusion properties of octane allow single peptides to come into contact within a few nanoseconds, whereas a full representation of the lipid environment would make the association almost impossible to simulate because of a much lower lateral diffusion coefficient, slowing down the diffusion and therefore shifting the timescale of the association process far beyond the time window accessible by simulations.

The starting configuration of our simulations is shown in Fig. 2 A, with the 36 peptides equally distributed over the simulation box. The peptide/octane ratio is 1:84, which is roughly equivalent to a peptide/DPC ratio of 1:56 (on the basis of the lipid hydrophobic tail). This concentration is only 3–6 times larger than the concentration used by Gratkowski et al. in their experiments, and allows us to follow association events in a reasonable time window obtaining results that are close to the experimental data in DPC (Gratkowski et al., 2002).

View larger version (19K): [in this window] [in a new window]   FIGURE 2  Snapshots of the trajectory (core region of peptides only) showing the starting structure (A) and the final structure (B) from our simulation. The MS1 helices are shown in blue, the Asn residues are shown in orange, and the rectangle represents the simulation box. The periodic images are added for clarity. The final structure is a snapshot of a very dynamic and transient system with nine stable dimers whereas all other close helices are rapidly diffusing.

	RESULTS

TOP ABSTRACT INTRODUCTION METHODS RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
Structure of associated peptides
The initial aggregation process of MS1 (Fig. 3) is characterized by fast formation of transient aggregates. Contacts between peptides can be identified after the first nanosecond with their numbers increasing over the entire 45 ns trajectory. The majority of these contact events do not lead to the formation of stable structures, but are transient, disassociating on a picosecond timescale, and involve only a few atoms on each helix. The average number of dimers with at least two atoms in van der Waal's contact increases with simulation time as shown in Fig. 3, reaching a final molecular fraction near fifty percent with the dimeric state being the dominant association form (Fig. 3 A middle panel). The inclusion of all atoms of each peptide chain shows a higher aggregation (Fig. 3 A, upper panel) rate that is not representative for the real fraction of associated peptides, and likely related to the limitations of our model membrane system that does not include the lipid headgroups. The fraction of close contacts of the Asn residue correlates well with the fraction of close contacts between peptides indicating that a large percentage of the associates are stabilized by the formation of interchain hydrogen bonds (see Fig. 3 A, lower panel). The final structure of our 45 ns trajectory is shown in Fig. 2 B. This figure represents a snapshot of a highly dynamic system with transiently associated aggregates and it contains a number of very stable hydrogen-bonded dimers.

View larger version (37K): [in this window] [in a new window]   FIGURE 3  Association state distribution of MS1 peptides (A) listed as monomeric (black), dimeric (red), trimeric (green), and tetrameric and larger (blue) focusing on the complete length of the peptides (upper panel), the hydrophobic core region (middle panel), and the Asn residue inside the membrane (lower panel). Dimer formation (B) showing the number of hydrogen bonds (black) between the Asn residues, and the center of mass distance (red) of the two helices of a dimer with the actual dimer structures displayed in the inserts. The Asn residues are highlighted in orange.

Only one hydrogen-bonded and tightly packed trimer formed during the 45 ns time window. This symmetric trimer developed from a dimer by the inclusion of a third peptide. The stability of this trimeric complex seems to be lower than that of dimers, as after several ns this complex changed to a more open structure and converted into two dimers with the involvement of a fourth peptide. A longer simulation would be required to show the formation of more of these trimeric structures allowing a more thorough analysis.

Dimer conformations
In our simulations we see a clear tendency toward the formation of left-handed coiled-coil structures over time. At the moment of dimer formation we find a slightly higher number of right-handed structures. Over time the amount of left-handed dimers increases at the expense of the right-handed conformation. The left-handed coiled-coil dimers are highly stable, and once formed these complexes generally remain associated for the remaining simulation time. The right-handed coiled-coil structures are less stable, as two of these initially right-handed coiled-coils change their conformation—making a transition toward a left-handed coiled-coil structure without disrupting the hydrogen bonds between the Asn residues as shown in Fig. 4. No transition from a left-handed toward a right-handed coiled-coil could be identified. At the endpoint of our simulation 5 left-handed coiled-coils and only three right-handed can be observed.

View larger version (55K): [in this window] [in a new window]   FIGURE 4  Snapshots showing a transition from the right- to a left-handed conformation. (A) Snapshot at 20 ns showing a right-handed crossing angle. (B) Snapshot at 40 ns showing an approximate parallel orientation. (C) Snapshot at 45 ns showing a left-handed coiled-coil structure. The Asn residues are highlighted in orange.

Types of interhelical hydrogen bonds
The Asn-mediated dimers are not restricted to one single conformation. Several conformations with hydrogen bonds across the dimer interface can be identified, with all keeping close packing and many hydrophobic interactions over the entire dimer interface (Fig. 4). The interhelical hydrogen bonding pattern differs between single dimers: i), the main group of dimeric structures is characterized by hydrogen bonds across the dimer interface involving only Asn side chains, discussed below; and ii), a second group of dimers is characterized by interhelical hydrogen bonds between one Asn side chain and the backbone of the second helix. The stability of these structures is lower compared to the previously described structures and these dimers are less frequent. The hydrogen bonds are much weaker, as the N-O distances are generally large, the hydrogen bonding angle is not ideal and because this type of interchain hydrogen bond disturbs the helical backbone hydrogen bonds. Additionally the inherent asymmetry leaves the Asn of the second peptide chain without stabilizing interactions.

Interhelical Asn mediated hydrogen bonds
We can identify two distinct populations of interhelical hydrogen-bond patterns of dimers as displayed in Fig. 5. Both species are very stable in the time window of our simulation showing hardly any interconversion. The first group consists of structures that are characterized by only one hydrogen bond between the two Asn side chains in line with the asymmetric hydrogen-bond described in the crystal structure of GCN4-P1 (O'Shea et al., 1991). The other potential hydrogen bond donor (amide group) or acceptor (carbonyl oxygen) is occasionally involved in a hydrogen bond to its own backbone or in an interchain hydrogen bond to the backbone of the second chain. The second dimeric form is a more symmetric structure with two interhelical hydrogen bonds across the dimer interface between the side chains of the Asn residues. The formation of these two hydrogen bonds permits both amide and carbonyl groups to be involved in one hydrogen bond each.

View larger version (12K): [in this window] [in a new window]   FIGURE 5  Model describing the monomer, dimer, and trimer structures and the interhelical hydrogen-bonding pattern present in our simulations.

	DISCUSSION

TOP ABSTRACT INTRODUCTION METHODS RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

The most common feature of the associated MS1 peptides is the Asn-mediated dimerization with many van der Waal's contacts along the remainder of the helices and some interaction between N-terminal amphipathic regions outside the octane slab. The polar Asn residue, situated in the dimer interface in the middle of the membrane spanning stretch, was identified by experiments to be crucial for the stability of MS1 oligomers. Indeed mutating this polar residue to Val, removing the only polar residue, largely abolished the association properties of MS1. In our simulations the same polar residue has been found to be important for stability. A clear correlation can be identified between the fraction of close contacts of the Asn residues and the presence of dimers and trimers (Fig. 3 A, middle and lower panel), outlining the important role this polar residue is playing in the association. These interhelical hydrogen bonds are a strong structural feature, determining stability with an additional contribution from van der Waal's contacts. The property of the Asn residue to be at the same time hydrogen-bonding donor and acceptor seems to be determinant for the structure-stabilizing effect, with the strength of a hydrogen bond expected to be large inside the hydrophobic environment of the bilayer.

The conformations of the side chains of this Asn residue have been investigated by NMR experiments by Zhou et al. (2000) using a GCN4-based leucine zipper peptide with identical residues in the dimer interface ("a" and "d" position) and Leu residues at all other positions. In micelles at least two distinguishable populations of the side chain amide protons of the Asn residue can be identified that slightly differ in their chemical shift as a result of differences in their local chemical environment on the NMR timescale. The Asn side chain shows fast rotations in our simulation, suggesting an averaging effect. Therefore smaller differences will probably not be resolved in the NMR spectra. The effect of the hydrogen bonding on the chemical shift should nevertheless be visible in an NMR spectrum in the micellar environment and could give rise to the two populations observed by Zhou et al. (2000), because hydrogen bonds have a large effect on the chemical shift and the rotation of the Asn side chains does not affect the stability of these hydrogen bonds that remain intact. Based on our simulations, we postulate a model of the monomer, dimer, and trimer structures and the hydrogen-bonding interactions of MS1 (shown in Fig. 5) that is similar to the model proposed by Zhou et al. (2000), the difference being that we suggest the existence of two different stable species of dimers. We propose that the monomeric MS1 peptides associate into two types of dimers, both relatively stable. One is a symmetric structure with two hydrogen bonds across the dimer interface, whereas the second conformation is characterized by only one hydrogen bond. The latter asymmetric structure is in line with the asymmetric hydrogen-bonding pattern described in the crystal structure of GCN4-P1. The relatively limited sampling of trimeric structures leaves some uncertainty on possible trimeric states. The trimer structure identified in our simulation is a very symmetric structure with the Asn residues in the middle of the helix bundle forming hydrogen bonds to the other peptides.

The formation of stable dimers without the involvement of hydrogen bonds is unexpected, as ultracentrifugation and SDS-PAGE experiments could not identify clearly detectable amounts of dimers in the N14V mutant (14), although showing traces of dimers in the case of the N14L mutation (Gratkowski et al., 2001). Our finding of these dimers are supported by results obtained by Gurezka et al. (1999; Gurezka and Langosch, 2001) who have shown that almost any combination of the hydrophobic residues LMIVF in the "a" and "d" positions of the heptad motif can form stable dimers in bacterial Escherichia coli membranes. Leu, Ile, and Val were particularly common in sequences that promoted dimer formation. Both these experimental data and our result suggest the presence of a low fraction of dimers stabilized by hydrophobic interactions without contributions from hydrogen bonds.

	CONCLUSIONS

TOP ABSTRACT INTRODUCTION METHODS RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
The conformation of MS1 has been expected to assume a left-handed coiled-coil structure like the parent water soluble GCN4-P1 structure. Although our data are not sufficient to determine the equilibrium distribution of structures, our simulations give a first indication of the structural preferences of MS1 as we can identify a clear tendency toward the formation of left-handed coiled-coils, confirming this prediction.

Hydrogen bonding between buried polar residues capable of being simultaneously hydrogen-bond donor and acceptor can clearly be identified to contribute dominantly to the stability of associated helices in membranes. In our simulations, only dimers that are held together by hydrogen bonds between the carboxamido group of the Asn side chains form very stable structures, whereas structures excluding such interactions are less stable, less frequent, and have short lifetimes. Hydrophobic interactions additionally contribute to the stability of tightly packed dimers, but do not seem to be determinant as seen from the time window of our simulations.

Two types of dimers with different interhelical hydrogen-bonding patterns, one symmetric, the other asymmetric, can be identified and both seem equally likely, suggesting a complex situation with two contemporaneously populated configurations with different binding modes.

The results of this study make us confident of the validity and usefulness of this large-scale approach to investigate the association behavior of membrane peptides. Further studies will investigate the equilibrium properties of MS1 and address the effect of single amino acid mutation.

	ACKNOWLEDGEMENTS

TOP ABSTRACT INTRODUCTION METHODS RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
This work was supported by the Canadian Institutes of Health Research. D.P.T. is an Alberta Heritage Foundation for Medical Research Scholar. J.L.M. is supported by studentships from Alberta Ingenuity and the Natural Sciences and Engineering Research Council of Canada. W.L.A. is supported by a studentship from the Alberta Heritage Foundation for Medical Research. We acknowledge computational support from WestGrid.

	FOOTNOTES

 
Abbreviations used: GCN4-P1, peptide from yeast transcription factor GCN4; MS1, model transmembrane peptide based on GCN4-P1; TM, transmembrane; DPC, dodecyl phosphocholine.

Submitted on May 1, 2004; accepted for publication June 23, 2004.

	REFERENCES

TOP ABSTRACT INTRODUCTION METHODS RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
Arkin, I. T., and A. T. Brunger. 1998. Statistical analysis of predicted transmembrane alpha-helices. Biochim. Biophys. Acta. 1429:113–128.[CrossRef][Medline]

Bargmann, C. I., M. C. Hung, and R. A. Weinberg. 1986. Multiple independent activations of the neu oncogene by a point mutation altering the transmembrane domain of p185. Cell. 45:649–657.[CrossRef][Medline]

Berendsen, H. J. C., J. P. M. Postma, W. F. van Gunsteren, A. DiNola, and J. R. Haak. 1984. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 81:3684–3689.[CrossRef]

Berendsen, H. J. C., J. P. M. Postma, W. F. van Gunsteren, and J. Hermans. 1981. In: Intermolecular Forces. Pullman B, editor. Reidel, Dordrecht, The Netherlands.

Berendsen, H. J. C., D. Van der Spoel, and R. Van Drunen. 1995. Gromacs—a message-passing parallel molecular-dynamics implementation. Comput. Phys. Comm. 91:43–56.[CrossRef]

Bowie, J. U. 1997. Helix packing angle preferences. Nat. Struct. Biol. 4:915–917.[CrossRef][Medline]

Bowie, J. U. 2000. Understanding membrane protein structure by design. Nat. Struct. Biol. 7:91–94.[CrossRef][Medline]

Choma, C., H. Gratkowski, J. D. Lear, and W. F. DeGrado. 2000. Asparagine-mediated self-association of a model transmembrane helix. Nat. Struct. Biol. 7:161–166.[CrossRef][Medline]

Crick, F. H. C. 1953. The packing of alpha-helices: simple coiled coils. Acta Crystallogr. 6:689–697.[CrossRef]

Darden, T., D. York, and L. Pedersen. 1993. Particle mesh Ewald—an N-log(N) method for Ewald sums in large systems. J. Chem. Phys. 98:10089–10092.[CrossRef]

Daura, X., K. Gademann, B. Jaun, D. Seebach, W. F. van Gunsteren, and A. E. Mark. 1999. Peptide folding: when simulation meets experiment. Angew. Chem. Int. Ed. 38:236–240.[CrossRef]

Essmann, U., L. Perera, M. L. Berkowitz, T. Darden, H. Lee, and L. G. Pedersen. 1995. A smooth particle mesh Ewald method. J. Chem. Phys. 103:8577–8593.[CrossRef]

Feenstra, K. A., B. Hess, and H. J. C. Berendsen. 1999. Improving efficiency of large time-scale molecular dynamics simulations of hydrogen-rich systems. J. Comput. Chem. 20:786–798.[CrossRef]

Fleming, K. G., and D. M. Engelman. 2001. Specificity in transmembrane helix-helix interactions can define a hierarchy of stability for sequence variants. Proc. Natl. Acad. Sci. USA. 98:14340–14344.[Abstract/Free Full Text]

Gonzalez, L., R. A. Brown, D. Richardson, and T. Alber. 1996a. Crystal structures of a single coiled-coil peptide in two oligomeric states reveal the basis for structural polymorphism. Nat. Struct. Biol. 3:1002–1010.[CrossRef][Medline]

Gonzalez, L., J. J. Plecs, and T. Alber. 1996b. An engineered allosteric switch in leucine-zipper oligomerization. Nat. Struct. Biol. 3:510–515.[CrossRef][Medline]

Gratkowski, H., Q. H. Dai, A. J. Wand, W. F. DeGrado, and J. D. Lear. 2002. Cooperativity and specificity of association of a designed transmembrane peptide. Biophys. J. 83:1613–1619.[Abstract/Free Full Text]

Gratkowski, H., J. D. Lear, and W. F. DeGrado. 2001. Polar side chains drive the association of model transmembrane peptides. Proc. Natl. Acad. Sci. USA. 98:880–885.[Abstract/Free Full Text]

Gurezka, R., R. Laage, B. Brosig, and D. Langosch. 1999. A heptad motif of leucine residues found in membrane proteins can drive self-assembly of artificial transmembrane segments. J. Biol. Chem. 274:9265–9270.[Abstract/Free Full Text]

Gurezka, R., and D. Langosch. 2001. In vitro selection of membrane-spanning leucine zipper protein-protein interaction motifs using POSSYCCAT. J. Biol. Chem. 276:45580–45587.[Abstract/Free Full Text]

Harbury, P. B., T. Zhang, P. S. Kim, and T. Alber. 1993. A switch between 2-stranded, 3-stranded and 4-stranded coiled coils in Gcn4 leucine-zipper mutants. Science. 262:1401–1407.[Abstract/Free Full Text]

Hess, B., H. Bekker, H. J. C. Berendsen, and J. G. E. M. Fraaije. 1997. LINCS: a linear constraint solver for molecular simulations. J. Comput. Chem. 18:1463–1472.[CrossRef]

Jones, D. T., W. R. Taylor, and J. M. Thornton. 1994. A mutation data matrix for transmembrane proteins. FEBS Lett. 339:269–275.[CrossRef][Medline]

Kataoka, M., T. W. Kahn, Y. Tsujiuchi, D. M. Engelman, and F. Tokunaga. 1992. Bacteriorhodopsin reconstituted from 2 individual helices and the complementary 5-helix fragment is photoactive. Photochem. Photobiol. 56:895–901.[Medline]

Langosch, D., and J. Heringa. 1998. Interaction of transmembrane helices by a knobs-into-holes packing characteristic of soluble coiled coils. Proteins. 31:150–159.[CrossRef][Medline]

Lear, J. D., H. Gratkowski, L. Adamian, J. Liang, and W. F. DeGrado. 2003. Position-dependence of stabilizing polar interactions of asparagine in transmembrane helical bundles. Biochemistry. 42:6400–6407.[CrossRef][Medline]

Li, R., N. Mitra, H. Gratkowski, G. Vilaire, R. Litvinov, C. Nagasami, J. W. Weisel, J. D. Lear, W. F. DeGrado, and J. S. Bennett. 2003. Activation of integrin alphaIIbbeta3 by modulation of transmembrane helix associations. Science. 300:795–798.[Abstract/Free Full Text]

Lindahl, E., B. Hess, and D. van der Spoel. 2001. GROMACS 3.0: a package for molecular simulation and trajectory analysis. J. Mol. Mod. 7:306–317.

Lino, G., D. N. Woolfson, and T. Alber. 1996. Buried polar residues and structural specificity in the GCN4 leucine zipper. Nat. Struct. Biol. 3:1011–1018.[CrossRef][Medline]

Liu, J. F., and B. Rost. 2001. Comparing function and structure between entire proteomes. Protein Sci. 10:1970–1979.[Abstract/Free Full Text]

Miyamoto, S., and P. A. Kollman. 1992. SETTLE—an analytical version of the SHAKE and RATTLE algorithm for rigid water models. J. Comput. Chem. 13:952–962.[CrossRef]

O'Shea, E. K., J. D. Klemm, P. S. Kim, and T. Alber. 1991. X-ray structure of the Gcn4 leucine zipper, a 2-stranded, parallel coiled coil. Science. 254:539–544.[Abstract/Free Full Text]

Partridge, A. W., R. A. Melnyk, and C. M. Deber. 2002. Polar residues in membrane domains of proteins: Molecular basis for helix-helix association in a mutant CFTR transmembrane segment. Biochemistry. 41:3647–3653.[CrossRef][Medline]

Popot, J. L., and D. M. Engelman. 1990. Membrane-protein folding and oligomerization—the 2-stage model. Biochemistry. 29:4031–4037.[CrossRef][Medline]

Popot, J. L., and D. M. Engelman. 2000. Helical membrane protein folding, stability, and evolution. Annu. Rev. Biochem. 69:881–922.[CrossRef][Medline]

Rees, D. C., H. Komiya, T. O. Yeates, J. P. Allen, and G. Feher. 1989. The bacterial photosynthetic reaction center as a model for membrane proteins. Annu. Rev. Biochem. 58:607–633.[Medline]

Scott, W. R. P., P. H. Hunenberger, I. G. Tironi, A. E. Mark, S. R. Billeter, J. Fennen, A. E. Torda, T. Huber, P. Kruger, and W. F. van Gunsteren. 1999. The GROMOS biomolecular simulation program package. J. Phys. Chem. A. 103:3596–3607.[CrossRef]

Smith, S. O., C. S. Smith, and B. J. Bormann. 1996. Strong hydrogen bonding interactions involving a buried glutamic acid in the transmembrane sequence of the neu/erbB-2 receptor. Nat. Struct. Biol. 3:252–258.[CrossRef][Medline]

Strelkov, S. V., and P. Burkhard. 2002. Analysis of alpha-helical coiled coils with the program TWISTER reveals a structural mechanism for stutter compensation. J. Struct. Biol. 137:54–64.[CrossRef][Medline]

Walther, D., C. Springer, and F. E. Cohen. 1998. Helix-helix packing angle preferences for finite helix axes. Proteins. 33:457–459.[CrossRef][Medline]

White, S. H., and W. C. Wimley. 1999. Membrane protein folding and stability: physical principles. Annu. Rev. Biophys. Biomol. Struct. 28:319–365.[CrossRef][Medline]

Zhou, F. X., M. J. Cocco, W. P. Russ, A. T. Brunger, and D. M. Engelman. 2000. Interhelical hydrogen bonding drives strong interactions in membrane proteins. Nat. Struct. Biol. 7:154–160.[CrossRef][Medline]

Zhou, F. X., H. J. Merianos, A. T. Brunger, and D. M. Engelman. 2001. Polar residues drive association of polyleucine transmembrane helices. Proc. Natl. Acad. Sci. USA. 98:2250–2255.[Abstract/Free Full Text]

This article has been cited by other articles:

				T. Carpenter, P. J. Bond, S. Khalid, and M. S. P. Sansom Self-Assembly of a Simple Membrane Protein: Coarse-Grained Molecular Dynamics Simulations of the Influenza M2 Channel Biophys. J., October 15, 2008; 95(8): 3790 - 3801. [Abstract] [Full Text] [PDF]

				J. L. MacCallum, W.F. D. Bennett, and D. P. Tieleman Partitioning of Amino Acid Side Chains into Lipid Bilayers: Results from Computer Simulations and Comparison to Experiment J. Gen. Physiol., April 30, 2007; 129(5): 371 - 377. [Full Text] [PDF]

				S. Esteban-Martin and J. Salgado Self-Assembling of Peptide/Membrane Complexes by Atomistic Molecular Dynamics Simulations Biophys. J., February 1, 2007; 92(3): 903 - 912. [Abstract] [Full Text] [PDF]

				A. C. V. Johansson and E. Lindahl Amino-Acid Solvation Structure in Transmembrane Helices from Molecular Dynamics Simulations Biophys. J., December 15, 2006; 91(12): 4450 - 4463. [Abstract] [Full Text] [PDF]

				P. A. Lee, G. L. Orriss, G. Buchanan, N. P. Greene, P. J. Bond, C. Punginelli, R. L. Jack, M. S. P. Sansom, B. C. Berks, and T. Palmer Cysteine-scanning Mutagenesis and Disulfide Mapping Studies of the Conserved Domain of the Twin-arginine Translocase TatB Component J. Biol. Chem., November 10, 2006; 281(45): 34072 - 34085. [Abstract] [Full Text] [PDF]

				J. M. Johnston, G. A. Cook, J. M. Tomich, and M. S. P. Sansom Conformation and Environment of Channel-Forming Peptides: A Simulation Study Biophys. J., March 15, 2006; 90(6): 1855 - 1864. [Abstract] [Full Text] [PDF]

This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Stockner, T. Articles by Tieleman, D. P. Search for Related Content PubMed PubMed Citation Articles by Stockner, T. Articles by Tieleman, D. P.