MultiCoil: a program for predicting two- and three-stranded coiled coils

A new multidimensional scoring approach for identifying and distinguishing trimeric and dimeric coiled coils is implemented in the MultiCoil program. The program extends the two-stranded coiled-coil prediction program PairCoil to the identification of three-stranded coiled coils. The computations are based upon data gathered from a three-stranded coiled-coil database comprising 6,319 amino acid residues, as well as from the previously constructed two-stranded coiled-coil database. In addition to identifying coiled coils not predicted by the two-stranded database programs, MultiCoil accurately classifies the oligomerization states of known dimeric and trimeric coiled coils. Analysis of the MultiCoil scores provides insight into structural features of coiled coils, and yields estimates that 0.9% of all protein residues form three-stranded coiled coils and that 1.5% form two-stranded coiled coils. The MultiCoil program is available at http:/(/)theory.lcs.mit.edu/multicoil.     

The coiled-coil helix in the neck of kinesin

Kinesin is a microtubule-dependent motor protein. We have recently determined the X-ray structure of monomeric and dimeric kinesin from rat brain. The dimer consists of two motor domains, held together by their alpha-helical neck domains forming a coiled coil. Here we analyze the nature of the interactions in the neck domain (residues 339-370). Overall, the neck helix shows a heptad repeat (abcdefg)n typical of coiled coils, with mostly nonpolar residues in positions a and d. However, the first segment (339-355) contains several nonclassical residues in the a and d positions which tend to weaken the hydrophobic interaction along the common interface. Instead, stabilization is achieved by a hydrophobic "coat" formed by the a and d residues and the long aliphatic moieties of lysines and glutamates, extending away from the coiled-coil core. By contrast, the second segment of the kinesin neck (356-370) shows a classical leucine zipper pattern in which most of the hydrophobic residues are buried at the highly symmetrical dimer interface. The end of the neck reveals the structure of a potential coiled-coil "trigger" sequence.     

Self-association of herpes simplex virus type 1 ICP35 is via coiled-coil interactions and promotes stable interaction with the major capsid protein

The ordered copolymerization of viral proteins to form the herpes simplex virus (HSV) capsid occurs within the nucleus of the infected cell and is a complex process involving the products of at least six viral genes. In common with capsid assembly in double-stranded DNA bacteriophages, HSV capsid assembly proceeds via the assembly of an outer capsid shell around an interior scaffold. This capsid intermediate matures through loss of the scaffold and packaging of the viral genomic DNA. The interior of the HSV capsid intermediate contains the viral protease and assembly protein which compose the scaffold. Proteolytic processing of these proteins is essential for and accompanies capsid maturation. The assembly protein (ICP35) is the primary component of the scaffold, and previous studies have demonstrated it to be capable of intermolecular association with itself and with the major capsid protein, VP5. We have defined structural elements within ICP35 which are responsible for intermolecular self-association and for interaction with VP5. Yeast (Saccharomyces cerevisiae) two-hybrid assays and far-Western studies with purified recombinant ICP35 mapped a core self-association domain between Ser165 and His219. Site-directed mutations in this domain implicate a putative coiled coil in ICP35 self-association. This coiled-coil motif is highly conserved within the assembly proteins of other alpha herpesviruses. In the two-hybrid assay the core self-association domain was sufficient to mediate stable self-association only in the presence of additional structural elements in either N- or C-terminal flanking regions. These regions also contain conserved sequences which exhibit a high propensity for alpha helicity and may contribute to self-association by forming additional short coiled coils. Our data supports a model in which ICP35 molecules have an extended conformation and associate in parallel orientation through homomeric coiled-coil interactions. In additional two-hybrid experiments we evaluated ICP35 mutants for association with VP5. We discovered that in addition to the C-terminal 25 amino acids of ICP35, previously shown to be required for VP5 binding, an additional upstream region was required. This region is between Ser165 and His234 and contains the core self-association domain. Site-directed mutations and construction of chimeric molecules in which the self-association domain of ICP35 was replaced by the GCN4 leucine zipper indicated that this region contributes to VP5 binding through mediating self-association of ICP35 and not through direct binding interactions. Our results suggest that self-association of ICP35 strongly promotes stable association with VP5 in vivo and are consistent with capsid formation proceeding via formation of stable subassemblies of ICP35 and VP5 which subsequently assemble into capsid intermediates in the nucleus.     

Structural stability of short subsequences of the tropomyosin chain

The native tropomyosin molecule is a parallel, registered, alpha-helical coiled coil made from two 284-residue chains. Long excised subsequences (> or = 95 residues) form the same structure with comparable thermal stability. Here, we investigate local stability using shorter subsequences (20-50 residues) that are chemically synthesized or excised from various regions along the protein chain. Thermal unfolding studies of such shorter peptides by CD in the same solvent medium used in extant studies of the parent protein indicate very low helix content, almost no coiled-coil formation, and high thermal lability of such secondary structure as does form. This behavior is in stark contrast to extant data on leucine-zipper peptides and short "designed" synthetic peptides, many of which have high alpha-helix content and form highly stable coiled coils. The existence of short coiled coils calls into question the older idea that short subsequences of a protein have little structure. The present study supports the older view, at least in its application to tropomyosin. The intrinsic local alpha-helical propensity and helix-helix interaction in this prototypical alpha-helical protein is sufficiently weak as to require not only dimerization, but macro-molecular amplification in order to attain its native conformation in common benign media near neutral pH.     

Subdomain folding of the coiled coil leucine zipper from the bZIP transcriptional activator GCN4

One popular model for protein folding, the framework model, postulates initial formation of secondary structure elements, which then assemble into the native conformation. However, short peptides that correspond to secondary structure elements in proteins are often only marginally stable in isolation. A 33-residue peptide (GCN4-p1) corresponding to the GCN4 leucine zipper folds as a parallel, two-stranded coiled coil [O'Shea, E.K., Klemm, J.D., Kim, P.S., & Alber, T.A. (1991) Science 254, 539-544]. Deletion of the first residue (Arg 1) results in local, N-terminal unfolding of the coiled coil, suggesting that a stable subdomain of GCN4-p1 can form. N- and C-terminal deletion studies result in a 23-residue peptide, corresponding to residues 8-30 of GCN4-p1, that folds as a parallel, two-stranded coil with substantial stability (the melting temperature of a 1 mM solution is 43 degrees C at pH 7). In contrast, a closely related 23-residue peptide (residues 11-33 of GCN4-p1) is predominantly unfolded, even at 0 degrees C, as observed previously for many isolated peptides of similar length. Thus, specific tertiary packing interactions between two short units of secondary structure can be energetically more important in stabilizing folded structure than secondary structure propensities. These results provide strong support for the notion that stable, cooperatively folded subdomains are the important determinants of protein folding.     

Why fibrous proteins are romantic

Here I give a personal account of the great history of fibrous protein structure. I describe how Astbury first recognized the essential simplicity of fibrous proteins and their paradigmatic role in protein structure. The poor diffraction patterns yielded by these proteins were then deciphered by Pauling, Crick, Ramachandran and others (in part by model building) to reveal alpha-helical coiled coils, beta-sheets, and the collagen triple helical coiled coil-all characterized by different local sequence periodicities. Longer-range sequence periodicities (or "magic numbers") present in diverse fibrous proteins, such as collagen, tropomyosin, paramyosin, myosin, and were then shown to account for the characteristic axial repeats observed in filaments of these proteins. More recently, analysis of fibrous protein structure has been extended in many cases to atomic resolution, and some systems, such as "leucine zippers," are providing a deeper understanding of protein design than similar studies of globular proteins. In the last sections, I provide some dramatic examples of fibrous protein dynamics. One example is the so-called "spring-loaded" mechanism for viral fusion by the hemagglutinin protein of influenza. Another is the possible conformational changes in prion proteins, implicated in "mad cow disease," which may be related to similar transitions in a variety of globular and fibrous proteins.     

Geometrical criteria for formation of coiled-coil structures of polypeptide chains

Crick's general formulas describing a coiled coil are expressed in a different form to combine the parameters of a coiled coil with the backbone dihedral angles of a polypeptide chain, assuming that the bond lengths and bond angles of the chain are fixed. While the existence of a low-energy coiled-coil conformation depends on energetic considerations, these formulas, which pertain to single-stranded structures and, by application of symmetry operations, to multistranded structures, provide the geometrical criteria for the existence of coiled coils. The concept of "the averaged structure of the minor helix", introduced here, makes it possible to relate the shape of the major helix to that of the minor helix. It is shown, in the analysis of a simple model of a single-stranded coiled-coil beta structure, that strong geometrical restrictions exist for the formation of coiled-coil structures from a given minor helix conformation of a polypeptide chain; these restrictions are expressed in a general form that is applicable to any coiled-coil of any number of residues in a repeat unit. As an application, the possible existence of a two-stranded coiled-coil antiparallel beta structure is considered, both geometrically and energetically, and discussed in relation to the observed twisted beta structures in globular proteins. The proposed coiled-coil models of alpha-helical proteins are also examined briefly.     

Surface salt bridges stabilize the GCN4 leucine zipper

We present a study of the role of salt bridges in stabilizing a simplified tertiary structural motif, the coiled-coil. Changes in GCN4 sequence have been engineered that introduce trial patterns of single and multiple salt bridges at solvent exposed sites. At the same sites, a set of alanine mutants was generated to provide a reference for thermodynamic analysis of the salt bridges. Introduction of three alanines stabilizes the dimer by 1.1 kcal/mol relative to the wild-type. An arrangement corresponding to a complex type of salt bridge involving three groups stabilizes the dimer by 1.7 kcal/ mol, an apparent elevation of the melting temperature relative to wild type of about 22 degrees C. While identifying local from nonlocal contributions to protein stability is difficult, stabilizing interactions can be identified by use of cycles. Introduction of alanines for side chains of lower helix propensity and complex salt bridges both stabilize the coiled-coil, so that combining the two should yield melting temperatures substantially higher than the starting species, approaching those of thermophilic sequences.     

Computational learning reveals coiled coil-like motifs in histidine kinase linker domains

The recent rapid growth of protein sequence databases is outpacing the capacity of researchers to biochemically and structurally characterize new proteins. Accordingly, new methods for recognition of motifs and homologies in protein primary sequences may be useful in determining how these proteins might function. We have applied such a method, an iterative learning algorithm, to analyze possible coiled coil domains in histidine kinase receptors. The potential coiled coils have not yet been structurally characterized in any histidine kinase, and they appear outside previously noted kinase homology regions. The learning algorithm uses a combination of established sequence patterns in known coiled coil proteins and histidine kinase sequence data to learn to recognize efficiently this coiled coil-like motif in the histidine kinases. The common appearance of the structural motif in a functionally important part of the receptors suggests hypotheses for kinase regulation and signal transduction.     

C repeats of the streptococcal M1 protein achieve the human serum albumin binding ability by flanking regions which stabilize the coiled-coil conformation

The M and M-like proteins of Streptococcus pyogenes are fibrous cell surface proteins. They have multiple binding sites for several human proteins and are composed of the C-terminal anchor domain, the alpha-helical coiled-coil domain, and the N-terminal non-coiled-coil domain. The coiled-coil domain of the M1 protein consists of repeat units called B, C, and D and a spacer unit S between B and C. Recombinant fragments A-B-S-C-D, A-B-S, B-S-C, S-C, S-C-D, C-D, and C of the coiled-coil domain were studied by analyzing their secondary structures and binding affinities to human serum albumin (HSA). As shown by circular dichroism, all fragments are in an alpha-helical conformation. C-D and S-C-D form coiled coils at room temperature and bind below 37 degrees C with high affinity to HSA. C-D and S-C-D unfold in two steps with Tm values of approximately 31 and approximately 65 degrees C; complex formation with HSA increases the unfolding temperatures. B-S-C has a lower alpha-helical content, a less pronounced coiled-coil conformation, and a reduced thermal stability, binds HSA weaker, and is only slightly stabilized by HSA binding in comparison to C-D and S-C-D. C and S-C are less stable than the other fragments and are not organized as coiled coils showing some features of alpha-helical single strands only below 20 degrees C, and binding of HSA was not observed. The results indicate that the formation of coiled-coil structures, supported by flanking D regions and, to a lesser extent also B regions, is essential for the binding of C repeat units to HSA.     

Structural adaptations in the specialized bacteriophage T4 co-chaperonin Gp31 expand the size of the Anfinsen cage

The Gp31 protein from bacteriophage T4 functionally substitutes for the bacterial co-chaperonin GroES in assisted protein folding reactions both in vitro and in vivo. But Gp31 is required for the folding and/or assembly of the T4 major capsid protein Gp23, and this requirement cannot be satisfied by GroES. The 2.3 A crystal structure of Gp31 shows that its tertiary and quaternary structures are similar to those of GroES despite the existence of only 14% sequence identity between the two proteins. However, Gp31 shows a series of structural adaptations which will increase the size and the hydrophilicity of the "Anfinsen cage," the enclosed cavity within the GroEL/GroES complex that is the location of the chaperonin-assisted protein folding reaction.     

Crystal structures of a single coiled-coil peptide in two oligomeric states reveal the basis for structural polymorphism

Each protein sequence generally adopts a single native fold, but the sequence features that confer structural uniqueness are not well understood. To define the basis for structural heterogeneity, we determined the high resolution X-ray crystal structures of a single GCN4 leucine-zipper mutant (Asn 16 to aminobutyric acid) in both dimeric and trimeric coiled-coil conformations. The mutant sequence is accommodated in two distinct structures by forming similarly-shaped packing surfaces with different sets of atoms. The trimer structure, in comparison to a previously-characterized trimeric mutant with substitutions in eight core residues, shows that the twist of individual helices and the helix-helix crossing angles can vary significantly to produce the most favoured packing arrangement.     

Thermodynamic characterization of the coupled folding and association of heterodimeric coiled coils (leucine zippers)

Folding thermodynamics of nine heterodimeric, parallel coiled coils were studied by isothermal titration calorimetry (ITC) and thermal unfolding circular dichroism measurements. The heterodimers were composed of an acidic and a basic 30-residue peptide, which when in isolation were monomeric and essentially unstructured. The reaction followed a two-state mechanism indicating that folding and association were coupled. delta Hfold, delta Sfold and delta Cp normalized per mol of residue were of the same magnitude as for monomeric globular proteins, hence the energetics of folding and association of the heterodimeric coiled coils was balanced similarly to the folding of a single polypeptide chain. Cavity creating Leu/Ala substitutions revealed strong and position-dependent energetic coupling between leucine residues in the hydrophobic core of the coiled coil. delta Gunfold (equivalent to -delta Gfold in the two-state reaction) was determined from thermal unfolding. Global stability curves were calculated according to the Gibbs-Helmholtz equation and using the combined free energy data from ITC and thermal unfolding. Maximum stabilities were between 15 and 37 degrees C and cold denaturation could be demonstrated by direct calorimetry. The stability curves were based on free energies of folding measured between 10 and 85 degrees C and under identical solvent conditions. This represents a novel experimental approach which circumvents the use of varying solvent conditions as is typically required to measure protein stability curves. Discrepancies were noticed between van't Hoff enthalpies deduced from thermal unfolding and measured by direct calorimetry. The discrepancies are thought to be due to residual ordered structure in the denatured single chains around room temperature but not near the transition midpoint temperature Tm. This demonstrates that over an extended temperature range the assumption of a common denatured state implicit in the van't Hoff analysis may not always be valid.     

Paramyxovirus phosphoproteins form homotrimers as determined by an epitope dilution assay, via predicted coiled coils

When HA epitope-tagged and untagged Sendai virus (SeV) P proteins are coexpressed and the products reacted with anti-HA, the untagged P protein is also selected because this protein is found as an oligomer. The oligomer was determined to be a homotrimer by coselection studies in which increasing amounts of untagged versus tagged protein were coexpressed, and these findings were extended to mumps virus, a member of the rubulavirus genus. The region of the SeV protein responsible for the oligomerization was localized to residues 344-411. Computer analysis of the 13 Paramyxovirus P proteins in the database revealed that all but one are predicted to form coiled coils in this region, the first of only two regions that can be aligned throughout the entire virus subfamily. The predicted coiled-coil region of the measles virus P protein, when grafted onto the C-terminus of the normally monomeric La protein, led to the efficient oligomerization of this reporter protein. The predicted coiled-coil region of these P proteins thus appears to be sufficient for oligomerization.     

Activity patterns in Glossina longipennis: a field study using different sampling methods

Protein folding in the cell is controlled at the levels of translation and post-translational modification, depends on a number of conserved proteins known as chaperones, and is catalyzed by specific enzymes, such as protein disulfide isomerase and peptidyl prolyl cis-trans isomerase. The chaperones stabilize folding intermediates and participate in assembly and disaggregation of supramolecular structures. Bacteriophage T4 is an especially convenient system for studying of protein folding mechanisms, since its genome encodes several virus-specific chaperones. In this review, the chaperones of phage T4 that take part in capsid formation (gp31 and gp40) and in folding and assembly of virion tail fibers (gp38, gp57A) have been considered. Protein encoded by gene 31 completely substitutes co-chaperonin GroES of the host cell in folding of the major capsid protein, gp23, aided by chaperonin GroEL. The product of gene 40, which is homologous to analogs of eukaryotic GroEL and peptidyl prolyl cis-trans isomerase, participates in assembly of gp20 while the formation of procapsid connector. The chaperone encoded by gene 57A is essential for folding and oligomerization of both long and short phage tail fibers. gp38, together with gp57A, participates in the formation of the distal part of the long fibers. This protein seems to represent a principally new group of chaperones that change steric structure of folded polypeptide. One phage chaperone, fibritin, encoded by gene wac (whiskers antigen control) and taking part in assembly the subunits of the long tail fibers is a constituent of the virion. Fibritin is a convenient model for studying mechanisms of folding and oligomerization of fibrous proteins due to its labile triple-stranded alpha-helical coiled-coil structure.     

Ribosome-binding domain of eukaryotic initiation factor-2 kinase GCN2 facilitates translation control

A family of protein kinases regulate translation initiation in response to cellular stresses by phosphorylation of eukaryotic initiation factor-2 (eIF-2). One family member from yeast, GCN2, contains a region homologous to histidyl-tRNA synthetases juxtaposed to the kinase catalytic domain. It is thought that uncharged tRNA accumulating during amino acid starvation binds to the synthetase-related sequences and stimulates phosphorylation of the alpha subunit of eIF-2. In this report, we define another domain in GCN2 that functions to target the kinase to ribosomes. A truncated version of GCN2 containing only amino acid residues 1467 to 1590 can independently associate with the translational machinery. Interestingly, this region of GCN2 shares sequence similarities with the core of the double-stranded RNA-binding domain (DRBD). Substitutions of the lysine residues conserved among DRBD sequences block association of GCN2 with ribosomes and impaired the ability of the kinase to stimulate translational control in response to amino acid limitation. Additionally, as found for other DRBD sequences, recombinant protein containing GCN2 residues 1467-1590 can bind double-stranded RNA in vitro, suggesting that interaction with rRNA mediates ribosome targeting. These results indicate that appropriate ribosome localization of the kinase is an obligate step in the mechanism leading to translational control by GCN2.     

Protein folding: optimized sequences obtained by simulated breeding in a minimalist model

For a minimalist model of protein folding, which we introduced recently, we investigate various methods to obtain folding sequences. A detailed study of random sequences shows that, for this model, such sequences usually do not fold to their ground states during simulations. Straight-forward techniques for the construction of folding sequences, based solely on the target structure, fail. We describe in detail an optimization algorithm, based on genetic algorithms, for the "simulated breeding" of folding sequences in this model. We find that, for any target structure studied, there is not only a single folding sequence but a patch of sequences in sequence space that fold to this structure. In addition, we show that, much as in real proteins, nonhomologous sequences may fold to the same target structure.     

Identification and mutagenesis of a highly conserved domain in troponin T responsible for troponin I binding: potential role for coiled coil interaction

Troponin T (TnT), a thin filament myofibrillar protein, is essential for the Ca2+ regulation of striated muscle contraction in vertebrates, both in vivo and in vitro. To understand the role of TnT in this process, its interaction with two other troponin components, troponin I (TnI) and troponin C (TnC) was examined by using the yeast two hybrid system, which is a genetic approach to detect protein-protein interactions. Computer assisted analysis of phylogenetically distant TnT amino acid sequences unveiled a highly conserved protein domain that is characterized by a heptad repeat (HR) motif with a potential for alpha-helical coiled coil formation. A similar, potentially coiled coil forming domain is also conserved in all known TnI sequences. These protein motifs appeared to be the regions where TnI-TnT interaction may take place. Deletions and point mutations in TnT, which disrupted its HR motif, severely reduced or abolished TnI binding, but binding to TnC was not affected, indicating that the TnT-TnI and TnT-TnC binary interactions can be uncoupled. Remarkably, the truncated fragments of TnT and TnI in which the HR motifs were retained showed binary interaction in the yeast two hybrid system. It was also observed that the formation of the TnT-TnI heterodimers is favored over the homodimers TnT-TnT and TnI-TnI. These results indicate that the evolutionarily conserved HR motifs may play a role in TnT-TnI dimerization, presumably through the formation of alpha-helical coiled coils.     

Extremely fast folding of a very stable leucine zipper with a strengthened hydrophobic core and lacking electrostatic interactions between helices

The dimer interface of a leucine zipper involves hydrophobic as well as electrostatic interactions between the component helices. Here we ask how hydrophobic effects and electrostatic repulsion balance the rate of folding and thermodynamic stability of a designed dimeric leucine zipper formed by the acidic peptide A that contains four repeating sequence units, (abcdefg)4. The aliphatic a and d residues of peptide A were the same as in the GCN4 leucine zipper but the e and g positions were occupied by Glu, which prevented folding above pH 6 because of electrostatic repulsion. Leucine zipper A2 was formed by protonation of the e and g side chains with a sharp transition midpoint at pH 5.2. Folding could be described by a two-state transition from two unfolded random coil monomers to a coiled coil dimer. There was a linear relationship between the logarithm of the rate constants and the number of repulsive charges on the folded leucine zipper dimer. The same linear relationship applied to the free energy of unfolding and the number of repulsive charges at thermodynamic equilibrium. Fully protonated peptide A folded at a near diffusion-limited rate (kon = 3 x 10(8) M-1 s-1), and the free energy of folding was -55 kJ mol-1 at 25 degrees C. The present work shows that protonation of Glu in positions e and g increases both the folding rate and the stability of the leucine zipper in the absence of any interhelical electrostatic interactions. Protonated Glu is proposed to act like a nonpolar residue and to strengthen the hydrophobic core by folding back toward the core residues in the a and d positions. This effect adds more to the free energy of unfolding and to the rate of folding than maximizing the number of salt bridges across the helix interface in an electrostatically stabilized heterodimeric leucine zipper [Wendt, H., Leder, L., H?rm?, H., Jelesarov, I., Baici, A., and Bosshard, H. R. (1997) Biochemistry 36, 204-213].