A designed protein with packing between left-handed and right-handed helices

A common motif in protein structures is the assembly of alpha-helices. Natural alpha-helical assemblies, such as helical bundles and coiled coils, consist of multiple right-handed alpha-helices. Here we design a protein complex containing both left-handed and right-handed helices, with peptides of D- and L-amino acids, respectively. The two peptides, D-Acid and L-Base, feature hydrophobic heptad repeats and are designed to pack against each other in a "knobs-into-holes" manner. In solution, the peptides form a stable, helical heterotetramer with tight packing in the most solvent-protected core. This motif may be useful for designing protease-resistant, helical D-peptide ligands against biological protein targets.     

Interaction of transmembrane helices by a knobs-into-holes packing characteristic of soluble coiled coils

Membrane-embedded protein domains frequently exist as α-helical bundles, as exemplified by photosynthetic reaction centers, bacteriorhodopsin, and cytochrome C oxidase. The sidechain packing between their transmembrane helices was investigated by a nearest-neighbor analysis which identified sets of interfacial residues for each analyzed helix–helix interface. For the left-handed helix–helix pairs, the interfacial residues almost exclusively occupy positions a, d, e, or g within a heptad motif (abcdefg) which is repeated two to three times for each interacting helical surface. The connectivity between the interfacial residues of adjacent helices conforms to the knobs-into-holes type of sidechain packing known from soluble coiled coils. These results demonstrate on a quantitative basis that the geometry of sidechain packing is similar for left-handed helix–helix pairs embedded in membranes and coiled coils of soluble proteins. The transmembrane helix–helix interfaces studied are somewhat less compact and regular as compared to soluble coiled coils and tolerate all hydrophobic amino acid types to similar degrees. The results are discussed with respect to previous experimental findings which demonstrate that specific interactions between transmembrane helices are important for membrane protein folding and/or oligomerization. Proteins 31:150–159, 1998. © 1998 Wiley-Liss, Inc.     

Comparison of helix interactions in membrane and soluble alpha-bundle proteins

Helix-helix interactions are important for the folding, stability, and function of membrane proteins. Here, two independent and complementary methods are used to investigate the nature and distribution of amino acids that mediate helix-helix interactions in membrane and soluble alpha-bundle proteins. The first method characterizes the packing density of individual amino acids in helical proteins based on the van der Waals surface area occluded by surrounding atoms. We have recently used this method to show that transmembrane helices pack more tightly, on average, than helices in soluble proteins. These studies are extended here to characterize the packing of interfacial and noninterfacial amino acids and the packing of amino acids in the interfaces of helices that have either right- or left-handed crossing angles, and either parallel or antiparallel orientations. We show that the most abundant tightly packed interfacial residues in membrane proteins are Gly, Ala, and Ser, and that helices with left-handed crossing angles are more tightly packed on average than helices with right-handed crossing angles. The second method used to characterize helix-helix interactions involves the use of helix contact plots. We find that helices in membrane proteins exhibit a broader distribution of interhelical contacts than helices in soluble proteins. Both helical membrane and soluble proteins make use of a general motif for helix interactions that relies mainly on four residues (Leu, Ala, Ile, Val) to mediate helix interactions in a fashion characteristic of left-handed helical coiled coils. However, a second motif for mediating helix interactions is revealed by the high occurrence and high average packing values of small and polar residues (Ala, Gly, Ser, Thr) in the helix interfaces of membrane proteins. Finally, we show that there is a strong linear correlation between the occurrence of residues in helix-helix interfaces and their packing values, and discuss these results with respect to membrane protein structure prediction and membrane protein stability.     

Parallel helix bundles and ion channels: molecular modeling via simulated annealing and restrained molecular dynamics.

A parallel bundle of transmembrane (TM) alpha-helices surrounding a central pore is present in several classes of ion channel, including the nicotinic acetylcholine receptor (nAChR). We have modeled bundles of hydrophobic and of amphipathic helices using simulated annealing via restrained molecular dynamics. Bundles of Ala20 helices, with N = 4, 5, or 6 helices/bundle were generated. For all three N values the helices formed left-handed coiled coils, with pitches ranging from 160 A (N = 4) to 240 A (N = 6). Pore radius profiles revealed constrictions at residues 3, 6, 10, 13, and 17. A left-handed coiled coil and a similar pattern of pore constrictions were observed for N = 5 bundles of Leu20. In contrast, N = 5 bundles of Ile20 formed right-handed coiled coils, reflecting loosened packing of helices containing beta-branched side chains. Bundles formed by each of two classes of amphipathic helices were examined: (a) M2a, M2b, and M2c derived from sequences of M2 helices of nAChR; and (b) (LSSLLSL)3, a synthetic channel-forming peptide. Both classes of amphipathic helix formed left-handed coiled coils. For (LSSLLSL)3 the pitch of the coil increased as N increased from 4 to 6. The M2c N = 5 helix bundle is discussed in the context of possible models of the pore domain of nAChR.     

Molecular dynamics (MD) investigations of preformed structures of the transmembrane domain of the oncogenic Neu receptor dimer in a DMPC bilayer

The critical Val/Glu mutation in the membrane spanning domain of the rat Neu receptor confers the ability for ligand-independent signaling and leads to increased dimerization and transforming ability. There is evidence that the two transmembrane interacting helices play a role in receptor activation by imposing orientation constraints to the intracellular tyrosine kinase domains. By using MD simulations we have attempted to discriminate between correct and improper helix-helix packing by examining the structural and energetic properties of preformed left-handed and right-handed structures in a fully hydrated DMPC bilayer. The best energetic balance between the residues at the helix-helix interface and the residues exposed to the lipids is obtained for helices in symmetrical left-handed interactions packed together via Glu side chain/Ala backbone interhelical hydrogen bonds. Analyses demonstrate the importance of the ATVEG motif in helix-helix packing and point to additional contacting residues necessary for association. Our findings, all consistent with experimental data, suggest that a symmetrical left-handed structure of the helices could be the transmembrane domain configuration that promotes receptor activation and transformation. The present study may provide further insight into signal transduction mechanisms of the ErbB/Neu receptors.     

Secondary structure prediction for the spectrin 106-amino acid segment, and a proposed model for tertiary structure

A collective secondary structure prediction for the human erythrocyte spectrin 106-residue repeat segment is developed, based on the sequences of nine segments that have been reported in the literature, utilizing a consensus of several secondary structure prediction methods for locating turn regions. The analysis predicts a five-fold structure, with three alpha-helices and two beta-strand regions, and differs from previous models on the lengths of the helices and the existence of beta-strand structure. We also demonstrate that this structural motif can be folded into tertiary structures that satisfy the experimental spectrin data and several general principles of protein organization.     

Conformational analysis of polynucleotides. I. The favorable left-handed helical model for the poly(8,2''-S-cycloadenylic acid) with high anti conformation.

1) Energy calculations have shown that poly (8,2'-S-cycloadenylic aicd) can form left-handed helices owing to the high anti conformation. 2) Two favorable left-handed helices are characterized by axial translation per residue (Z=4.3 and 3.6A) and by rotations per residue (theta= 40 degrees and -25 degrees). 3) The proposed helical models might be stable in aqueous solution and is well explicable of the optical property of this compound.     

Molecular dynamics simulations of the transmembrane domain of the oncogenic ErbB2 receptor dimer in a DMPC bilayer

Molecular dynamics simulations of an atomic model of the transmembrane domain of the oncogenic ErbB2 receptor dimer embedded in an explicit dimyristoylphosphatidylcholine (DMPC) bilayer were performed for more than 4 ns. The oncogenic Glu mutation in the membrane spanning segment plays a major role in tyrosine kinase activity and receptor dimerization, and is thought to be partly responsible for the structure of the transmembrane domain of the active receptor. MD results show that the interactions between the two transmembrane helices are characteristic of a left-handed packing as previously demonstrated from in vacuo simulations. Moreover, MD results reveal the absence of persistent hydrogen bonds between the Glu side chains in a membrane environment, which raise the question of the ability for Glu alone to stabilize the TM domain of the ErbB2 receptor. Interestingly the formation of the alpha-pi motif in the two ErbB2 transmembrane helices confirms the concept of intrinsic sequence-induced conformational flexibility. From a careful analysis of our MD results, we suggest that the left-handed helix-helix packing could be the key to correctly orient the intracellular domain of the activated receptor dimer. The prediction of such interactions from computer simulations represents a new step towards the understanding of signaling mechanisms.     

Structure of the DNA-binding domain of the response regulator PhoP from Mycobacterium tuberculosis

The PhoP-PhoR two-component signaling system from Mycobacterium tuberculosis is essential for the virulence of the tubercle bacillus. The response regulator, PhoP, regulates expression of over 110 genes. In order to elucidate the regulatory mechanism of PhoP, we determined the crystal structure of its DNA-binding domain (PhoPC). PhoPC exhibits a typical fold of the winged helix-turn-helix subfamily of response regulators. The structure starts with a four-stranded antiparallel beta-sheet, followed by a three-helical bundle of alpha-helices, and then a C-terminal beta-hairpin, which together with a short beta-strand between the first and second helices forms a three-stranded antiparallel beta-sheet. Structural elements are packed through a hydrophobic core, with the first helix providing a scaffold for the rest of the domain to pack. The second and third helices and the long, flexible loop between them form the helix-turn-helix motif, with the third helix being the recognition helix. The C-terminal beta-hairpin turn forms the wing motif. The molecular surfaces around the recognition helix and the wing residues show strong positive electrostatic potential, consistent with their roles in DNA binding and nucleotide sequence recognition. The crystal packing of PhoPC gives a hexamer ring, with neighboring molecules interacting in a head-to-tail fashion. This packing interface suggests that PhoPC could bind DNA in a tandem association. However, this mode of DNA binding is likely to be nonspecific because the recognition helix is partially blocked and would be prevented from inserting into the major groove of DNA. Detailed structural analysis and implications with respect to DNA binding are discussed.     

Structure prediction of the dimeric neu/ErbB-2 transmembrane domain from multi-nanosecond molecular dynamics simulations

Dimerization of the neu/ErbB-2 receptor tyrosine kinase is a necessary but not a sufficient step for signaling. Despite the efforts expended to identify the molecular interactions responsible for receptor-receptor contacts and particularly those involving the transmembrane domain, structural details are still unknown. In this work, molecular dynamics simulations of the helical transmembrane domain (TM) of neu and ErbB-2 receptors are used to predict their dimer structure both in the wild and oncogenic forms. A global conformational search method, applied to define the best orientations of parallel helices, showed an energetically favorable configuration with the specific mutation site within the interface, common for both the nontransforming and the transforming neu/ErbB-2 TM dimers. Starting from this configuration, a total of 10 simulations, about 1.4 ns each, performed in vacuum, without any constraints, show that the two helices preferentially wrap in left-handed interactions with a packing angle at about 20°. The resulting structures are nonsymmetric and the hydrogen bond network analysis shows that helices experience π local distortions that facilitate inter-helix hydrogen bond interactions and may result in a change in the helix packing, leading to a symmetric interface. For the mutated sequences, we show that the Glu side chain interacts directly with its cognate or with carbonyl groups of the facing backbone. We show that the connectivity between interfacial residues conforms to the knobs-into-holes packing mode of transmembrane helices. The dimeric interface described in our models is discussed with respect to mutagenesis studies. Received: 12 March 1999 / Revised version: 23 August 1999 / Accepted: 23 August 1999     

Transmembrane helix packing of ErbB/Neu receptor in membrane environment: a molecular dynamics study

Dimerization or oligomerization of the ErbB/Neu receptors are necessary but not sufficient for initiation of receptor signaling. The two intracellular domains must be properly oriented for the juxtaposition of the kinase domains allowing trans-phosphorylation. This suggests that the transmembrane (TM) domain acts as a guide for defining the proper orientation of the intracellular domains. Two structural models, with the two helices either in left-handed or in right-handed coiling have been proposed as the TM domain structure of the active receptor. Because experimental data do not distinguish clearly helix-helix packing, molecular dynamics (MD) simulations are used to investigate the energetic factors that drive Neu TM-TM interactions of the wild and the oncogenic receptor (Val664/Glu mutation) in DMPC or in POPC environments. MD results indicate that helix-lipid interactions in the bilayer core are extremely similar in the two environments and raise the role of the juxtamembrane residues in helix insertion and helix-helix packing. The TM domain shows a greater propensity to adopt a left-handed structure in DMPC, with helices in optimal position for strong inter-helical Hbonds induced by the Glu mutation. In POPC, the right-handed structure is preferentially formed with the participation of water in inter-helical Hbonds. The two structural arrangements of the Neu(TM) helices both with GG4 residue motif in close contact at the interface are permissible in the membrane environment. According to the hypothesis of a monomer-dimer equilibrium of the proteins it is likely that the bilayer imposes structural constraints that favor dimerization-competent structure responsible of the proper topology necessary for receptor activation.     

Transmembrane Helix Packing of ErbB/Neu Receptor in Membrane Environment: A Molecular Dynamics Study

Abstract

Dimerization or oligomerization of the ErbB/Neu receptors are necessary but not sufficient for initiation of receptor signaling. The two intracellular domains must be properly oriented for the juxtaposition of the kinase domains allowing trans-phosphorylation. This suggests that the transmembrane (TM) domain acts as a guide for defining the proper orientation of the intracellular domains.

Two structural models, with the two helices either in left-handed or in right-handed coiling have been proposed as the TM domain structure of the active receptor. Because experimental data do not distinguish clearly helix-helix packing, molecular dynamics (MD) simulations are used to investigate the energetic factors that drive Neu TM-TM interactions of the wild and the oncogenic receptor (Val664/Glu mutation) in DMPC or in POPC environments. MD results indicate that helix-lipid interactions in the bilayer core are extremely similar in the two environments and raise the role of the juxtamembrane residues in helix insertion and helix-helix packing. The TM domain shows a greater propensity to adopt a left-handed structure in DMPC, with helices in optimal position for strong inter-helical Hbonds induced by the Glu mutation. In POPC, the right-handed structure is preferentially formed with the participation of water in inter-helical Hbonds. The two structural arrangements of the NeuTM helices both with GG4 residue motif in close contact at the interface are permissible in the membrane environment. According to the hypothesis of a monomer-dimer equilibrium of the proteins it is likely that the bilayer imposes structural constraints that favor dimerization- competent structure responsible of the proper topology necessary for receptor activation.     

Association of a model transmembrane peptide containing gly in a heptad sequence motif

A peptide containing glycine at a and d positions of a heptad motif was synthesized to investigate the possibility that membrane-soluble peptides with a Gly-based, left-handed helical packing motif would associate. Based on analytical ultracentrifugation in C14-betaine detergent micelles, the peptide did associate in a monomer-dimer equilibrium, although the association constant was significantly less than that reported for the right-handed dimer of the glycophorin A transmembrane peptide in similar detergents. Fluorescence resonance energy transfer (FRET) experiments conducted on peptides labeled at their N-termini with either tetramethylrhodamine (TMR) or 7-nitrobenz-2-oxa-1,3-diazole (NBD) also indicated association. However, analysis of the FRET data using the usual assumption of complete quenching for NBD-TMR pairs in the dimer could not be quantitatively reconciled with the analytical ultracentrifugation-measured dimerization constant. This led us to develop a general treatment for the association of helices to either parallel or antiparallel structures of any aggregation state. Applying this treatment to the FRET data, constraining the dimerization constant to be within experimental uncertainty of that measured by analytical ultracentrifugation, we found the data could be well described by a monomer-dimer equilibrium with only partial quenching of the dimer, suggesting that the helices are most probably antiparallel. These results also suggest that a left-handed Gly heptad repeat motif can drive membrane helix association, but the affinity is likely to be less strong than the previously reported right-handed motif described for glycophorin A.     

How to untwist an alpha-helix: structural principles of an alpha-helical barrel

Recent crystallographic studies have revealed that 12 alpha-helices can pack in an anti-parallel fashion to form a hollow cylinder of nearly uniform radius. In this architecture, which we refer to as an alpha-barrel, the helices are inclined with respect to the cylindrical axis, and thus they curve and twist. As with conventional coiled-coils, the helices of the barrel associate via "knobs-into-holes" interactions; however, their packing is distinct in several important ways. First, the alpha-barrel helices untwist in comparison with the helices found in two-stranded coiled-coils and, as a consequence of this distortion, their knobs approach closely one end of the complementary holes. This effect defines a requirement for particular size and shape of the protruding residues, and it is associated with a relative axial translation of the paired helices. Second, as each helix packs laterally with two neighbours, the helices have two sequence patterns that are phased to match the two interfaces. The two types of interface are not equivalent and, as one travels around the circumference of the cylinder's interior, they alternate between one type where the knobs approach the holes straight-on, and a second type in which they are inclined. The choice of amino acid depends on the interface type, with small hydrophobic side-chains preferred for the direct contacts and larger aliphatic side-chains for the inclined contacts. Third, small residues are found preferentially on the inside of the tube, in order to make the "wedge" angle between helices compatible with a 12-member tube. Finally, hydrogen-bonding interactions of side-chains within and between helices support the assembly. Using these salient structural features, we present a sequence template that is compatible with some underlying rules for the packing of helices in the barrel, and which may have application to the design of higher-order assemblies from peptides, such as nano-tubes. We discuss the general implications of relative axial translation in coiled-coils and, in particular, the potential role that this movement could play in allosteric mechanisms.     

Aromatic di-alanine repeats (AdAR) are structural motifs characteristic of the soluble N-ethylmaleimide-sensitive factor attachment protein (SNAP) family

The aromatic di-alanine repeat is a novel 12-amino acid-long motif constituting alternate small and large hydrophobic residues that mediate the close packing of alpha-helices. A hidden Markov model profile was constructed from the motifs initially described in Soluble N-ethyl maleimide-sensitive factor attachment proteins (SNAP), a family of soluble proteins involved in intracellular membrane fusion. Scanning different sets of protein sequences showed unambiguously that this profile defines a structural motif independent of the tetratrico peptide repeat, another widespread alpha-helical motif. In addition to SNAP, aromatic di-alanine repeats are found in selective LIM homeodomain binding proteins (SLB) and in proteins from the Pyrococcus and Archaeoglobus prokaryotes.     

General architecture of the alpha-helical globule

A model is presented for the arrangement of alpha-helices in globular proteins. In the model, helices are placed on certain ribs of "quasi-spherical" polyhedra. The polyhedra are chosen so as to allow the close packing of helices around a hydrophobic core and to stress the collective interactions of the individual helices. The model predicts a small set of stable architectures for alpha-helices in globular proteins and describes the geometries of the helix packings. Some of the predicted helix arrangements have already been observed in known protein structures; others are new. An analysis of the three-dimensional structures of all proteins for which co-ordinates are available shows that the model closely approximates the arrangements and packing of helices actually observed. The average deviations of the real helix axes from those in the model polyhedra is +/- 20 degrees in orientation and +/- 2 A in position (1 A = 0.1 nm). We also show that for proteins that are not homologous, but whose helix arrangements are described by the same polyhedron, the root-mean-square difference in the position of the C alpha atoms in the helices is 1.6 to 3.0 A.     

Analysis of alpha-helical coiled coils with the program TWISTER reveals a structural mechanism for stutter compensation

Alpha-helical coiled coils represent a widespread protein structure motif distinguished by a seven-residue periodicity of apolar residues in the primary sequence. A characteristic "knobs-into-holes" packing of these residues into a hydrophobic core results in a superhelical, usually left-handed, rope of two or more alpha-helices. Such a geometry can be parameterized. For this purpose, a new computer program, TWISTER, was developed. With the three-dimensional coordinates as input, TWISTER uses an original algorithm to determine the local coiled-coil parameters as a function of residue number. In addition, heptad positions are assigned based on structural criteria. It is known that frequently encountered discontinuities in the heptad repeat, such as stutters and skips, can be tolerated within a continuous coiled coil but result in a local distortion of its geometry. This was explored in detail with the help of TWISTER for several two- and three-stranded coiled coils. Depending on the particular protein, stutters were found to be compensated locally by an unwinding of the superhelix, alpha-helical unwinding, or both. In the first case, there is often a local switch from a left-handed to a right-handed superhelix. In general, the geometrical distortion is confined to about two alpha-helical turns at either side of the stutter. Furthermore, stutters result in a local increase of the coiled-coil radius.     

Open-and-shut cases in coiled-coil assembly: alpha-sheets and alpha-cylinders

The coiled coil is a ubiquitous protein-folding motif. It generally is accepted that coiled coils are characterized by sequence patterns known as heptad repeats. Such patterns direct the formation and assembly of amphipathic alpha-helices, the hydrophobic faces of which interface in a specific manner first proposed by Crick and termed "knobs-into-holes packing". We developed software, SOCKET, to recognize this packing in protein structures. As expected, in a trawl of the protein data bank, we found examples of canonical coiled coils with a single contiguous heptad repeat. In addition, we identified structures with multiple, overlapping heptad repeats. This observation extends Crick's original postulate: Multiple, offset heptad repeats help explain assemblies with more than two helices. Indeed, we have found that the sequence offset of the multiple heptad repeats is related to the coiled-coil oligomer state. Here we focus on one particular sequence motif in which two heptad repeats are offset by two residues. This offset sets up two hydrophobic faces separated by approximately 150 degrees -160 degrees around the alpha-helix. In turn, two different combinations of these faces are possible. Either similar or opposite faces can interface, which leads to open or closed multihelix assemblies. Accordingly, we refer to these two forms as alpha-sheets and alpha-cylinders. We illustrate these structures with our own predictions and by reference to natural variants on these designs that have recently come to light.     

Discovery of a significant,nontopological preference for antiparallel alignment of helices with parallel regions in sheets

To help elucidate the role of secondary structure packing preferences in protein folding, here we present an analysis of the packing geometry observed between alpha-helices and between alpha-helices and beta-sheets in 1316 diverse, nonredundant protein structures. Finite-length vectors were fit to the alpha-carbon atoms in each of the helices and strands, and the packing angle between the vectors, Omega, was determined at the closest point of approach within each helix-helix or helix-sheet pair. Helix-sheet interactions were found in 391 of the proteins, and the distributions of Omega values were calculated for all the helix-sheet and helix-helix interactions. The packing angle preferences for helix-helix interactions are similar to those previously observed. However, analysis of helix-strand packing preferences uncovered a remarkable tendency for helices to align antiparallel to parallel regions of beta-sheets, independent of the topological constraints or prevalence of beta-alpha-beta motifs in the proteins. This packing angle preference is significantly diminished in helix interactions involving mixed and antiparallel beta-sheets, suggesting a role for helix-sheet dipole alignment in guiding supersecondary structure formation in protein folding. This knowledge of preferred packing angles can be used to guide the engineering of stable protein modules.