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.Ics.mit.edu/multicoil.     

Predicting oligomerization states of coiled coils.

An algorithm based on the profile method was developed that faithfully distinguishes between the amino acid sequences of dimeric and trimeric coiled coils. Normalized sequence profiles derived from nonhomologous, two- and three-stranded, coiled-coil sequences with unambiguous registers were used to assign dimer and trimer propensities to test sequences. The difference between the dimer and trimer profile scores accurately reflected the preferred oligomerization state. The method relied on two strategies that may be generally applicable to profile calculations--profile values of solvent-exposed residues and of amino acids that were underrepresented in the data-base were given zero weight. Differences between the dimer and trimer profiles revealed sequence patterns that match and extend experimental studies of oligomer specification.     

The role of position a in determining the stability and oligomerization state of alpha-helical coiled coils: 20 amino acid stability coefficients in the hydrophobic core of proteins

We describe here a systematic investigation into the role of position a in the hydrophobic core of a model coiled-coil protein in determining coiled-coil stability and oligomerization state. We employed a model coiled coil that allowed the formation of an extended three-stranded trimeric oligomerization state for some of the analogs; however, due to the presence of a Cys-Gly-Gly linker, unfolding occurred from the same two-stranded monomeric oligomerization state for all of the analogs. Denaturation from a two-stranded state allowed us to measure the relative contribution of 20 different amino acid side chains to coiled-coil stability from chemical denaturation profiles. In addition, the relative hydrophobicity of the substituted amino acid side chains was assessed by reversed-phase high-performance liquid chromatography and found to correlate very highly (R = 0.95) with coiled-coil stability. We also determined the effect of position a in specifying the oligomerization state using ultracentrifugation as well as high-performance size-exclusion chromatography. We found that nine of the analogs populated one oligomerization state exclusively at peptide concentrations of 50 microM under benign buffer conditions. The Leu-, Tyr-, Gln-, and His-substituted analogs were found to be exclusively three-stranded trimers, while the Asn-, Lys-, Orn-, Arg-, and Trp-substituted analogs formed exclusively two-stranded monomers. Modeling results for the Leu-substituted analog showed that a three-stranded oligomerization state is preferred due to increased side-chain burial, while a two-stranded oligomerization state was observed for the Trp analog due to unfavorable cavity formation in the three-stranded state.     

Sequence divergence of coiled coils--structural rods, myosin filament packing, and the extraordinary conservation of cohesins

The amino acid sequences of the long, anti-parallel coiled coils of the cohesin subunits SMC1 and SMC3 are almost totally conserved in mammals. To understand this exceptional conservation more broadly, we analyzed amino acid sequence variation for several groups of coiled-coil proteins. Some long coiled coils, including giantin, NuMA, and Ndc80p/Nuf2p diverge approximately 20% from humans to rodents, suggesting they function as spacer rods, whose sequence divergence is constrained only by the need to maintain the coiled-coil structure. Other coiled coils such as skeletal muscle myosin, intermediate filaments, and the lamins diverge only 1-3%. We suggest that this sequence divergence is constrained by the extensive packing contacts over the entire surface of the coiled-coil. The coiled coils of SMC5/6 and SMC2/4 (condensin) are slightly more constrained than the presumed spacer rods, diverging 10-15%. Conversely, the coiled coils of SMC1/3 (cohesin) diverge only 0.0-1.0%. This extreme constraint suggests that the entire surface of the coiled coil is intimately involved in the mechanism of sister chromatid cohesion. Direct binding of the coiled coils to chromatin, or perhaps the need to avoid such binding, are two possible mechanisms. Finally, analysis of the heptad repeat shows that the a and d positions are more constrained in spacer rods, and the bcefg positions more constrained in skeletal muscle myosin.     

A coiled-coil oligomerization domain of Bcr is essential for the transforming function of Bcr-Abl oncoproteins.

In Philadelphia chromosome-positive human leukemias, the c-abl proto-oncogene on chromosome 9 becomes fused to the bcr gene on chromosome 22, and chimeric Bcr-Abl proteins are produced. The fused Bcr sequences activate the tyrosine kinase, actin-binding, and transforming functions of Abl. Activation of the Abl transforming function has been shown to require two distinct domains of Bcr: domain 1 (Bcr amino acids 1 to 63) and domain 2 (Bcr amino acids 176 to 242). The amino acid sequence of domain 1 indicates that it may be a coiled-coil oligomerization domain. We show here that domain 1 of Bcr forms a homotetramer. Tetramerization of Bcr-Abl through Bcr domain 1 correlates with activation of the tyrosine kinase and F-actin-binding functions of Abl. Disruption of the coiled coil by insertional mutagenesis inactivates the oligomerization function as well as the ability of Bcr-Abl to transform Rat-1 fibroblasts or to abrogate interleukin-3 dependence in lymphoid cells. These results strongly suggest that Bcr-Abl oligomers are the active entities in transformation.     

Multicoil2: predicting coiled coils and their oligomerization states from sequence in the twilight zone

The alpha-helical coiled coil can adopt a variety of topologies, among the most common of which are parallel and antiparallel dimers and trimers. We present Multicoil2, an algorithm that predicts both the location and oligomerization state (two versus three helices) of coiled coils in protein sequences. Multicoil2 combines the pairwise correlations of the previous Multicoil method with the flexibility of Hidden Markov Models (HMMs) in a Markov Random Field (MRF). The resulting algorithm integrates sequence features, including pairwise interactions, through multinomial logistic regression to devise an optimized scoring function for distinguishing dimer, trimer and non-coiled-coil oligomerization states; this scoring function is used to produce Markov Random Field potentials that incorporate pairwise correlations localized in sequence. Multicoil2 significantly improves both coiled-coil detection and dimer versus trimer state prediction over the original Multicoil algorithm retrained on a newly-constructed database of coiled-coil sequences. The new database, comprised of 2,105 sequences containing 124,088 residues, includes reliable structural annotations based on experimental data in the literature. Notably, the enhanced performance of Multicoil2 is evident when tested in stringent leave-family-out cross-validation on the new database, reflecting expected performance on challenging new prediction targets that have minimal sequence similarity to known coiled-coil families. The Multicoil2 program and training database are available for download from http://multicoil2.csail.mit.edu.     

Removing an interhelical salt bridge abolishes coiled-coil formation in a de novo designed peptide

Alpha-helical coiled coils represent a common protein oligomerization motif that are mainly stabilized by hydrophobic interactions occurring along their coiled-coil interface, the so-called hydrophobic seam. We have recently de novo designed and optimized a series of two-heptad repeat long coiled-coil peptides which are further stabilized by a complex network of inter- and intrahelical salt bridges. Here we have extended the de novo design of such two heptad-repeat long peptides by removing the central and most important g-e' Arg to Glu (g-e'RE) ionic interhelical interaction and replacing these residues by alanine residues. The effect of the missing interhelical ionic interaction on coiled-coil formation and stability has been analyzed by CD spectroscopy, analytical ultracentrifugation, and X-ray crystallography. We show that the peptide, while being highly alpha-helical, is no longer able to form a parallel coiled-coil structure but rather assumes an octameric globular helical assembly devoid of any coiled-coil interactions.     

α-Helical coiled-coil oligomerization domains in extracellular proteins

Subunit oligomerization of many proteins is mediated by α-helical coiled-coil domains. 3,4-Hydrophobic heptad repeat sequences, the characteristic feature of the coiled-coil protein folding motif, have been found in a wide variety of gene products including cytoskeletal, nuclear, muscle, cell surface, extracellular, plasma, bacterial, and viral proteins. Whereas the majority of coiled-coil structures is represented by intracellular α-helical bundles that contain two polypeptide chains, examples of extracellular coiled-coil proteins are fewer in number. Most proteins located in the extracellular space form three-stranded α-helical assemblies. Recently, five-stranded coiled coils have been identified in thrombospondins 3 and 4 in cartilage oligomeric matrix protein, and the formation of a heterotetramer has been observed in in vitro studies with the recombinant asialoglycoprotein receptor oligomerization domain. Coiled-coil domains in laminins and probably also in tenascins and thrombospondins are responsible for the formation of tissue-specific isoforms by selective oligomerization of different polypeptide chains.     

Role of hydrophobic clusters in the stability of alpha-helical coiled coils and their conversion to amyloid-like beta-sheets

We designed a library of short peptides using standard rules for coiled-coil assembly. Depending on the composition of amino acids in the non-interacting region of the coiled coil (positions b, c, and f) these peptides are able to convert from alpha-helical to beta-sheet secondary structure. This type of transition is observed in amyloid-like proteins and is a key feature associated with many types of neurodegenerative diseases. Studies on peptides that are 14 amino acids in length indicated that positioning hydrophobic amino acids at an f position within a heptad repeat accelerated the rate of conformational conversion as compared to that at a c position. We believe that this occurs because of the formation of a hydrophobic pocket that preferentially stabilizes beta-sheets over alpha-helices. This effect was also observed in longer 21 amino acid peptides. Our study shows that the relative rates of structural conversion correlate with the formation of a continuous three-amino-acid hydrophobic patch consisting of amino acids in the d, f, and a positions and not on the secondary structure propensities of the individual amino acids. The sequence-structure relationship observed in this study will be used to help understand the mechanism of amyloid fiber formation and design future coiled-coil and beta-sheet-forming peptide systems.     

Probing designability via a generalized model of helical bundle geometry

Because the space of folded protein structures is highly degenerate, with recurring secondary and tertiary motifs, methods for representing protein structure in terms of collective physically relevant coordinates are of great interest. By collapsing structural diversity to a handful of parameters, such methods can be used to delineate the space of designable structures (i.e., conformations that can be stabilized with a large number of sequences)—a crucial task for de novo protein design. We first demonstrate this on natural α-helical coiled coils using the Crick parameterization. We show that over 95% of known coiled-coil structures are within  1-Å C_α root mean square deviation of a Crick-ideal backbone. Derived parameters show that natural geometric space of coiled coils is highly restricted and can be represented by “allowed” conformations amidst a potential continuum of conformers. Allowed structures have (1) restricted axial offsets between helices, which differ starkly between parallel and anti-parallel structures; (2) preferred superhelical radii, which depend linearly on the oligomerization state; (3) pronounced radius-dependent a- and d-position amino acid propensities; and (4) discrete angles of rotation of helices about their axes, which are surprisingly independent of oligomerization state or orientation. In all, we estimate the space of designable coiled-coil structures to be reduced at least 160-fold relative to the space of geometrically feasible structures. To extend the benefits of structural parameterization to other systems, we developed a general mathematical framework for parameterizing arbitrary helical structures, which reduces to the Crick parameterization as a special case. The method is successfully validated on a set of non-coiled-coil helical bundles, frequent in channels and transporter proteins, which show significant helix bending but not supercoiling. Programs for coiled-coil parameter fitting and structure generation are provided via a web interface at http://www.gevorggrigoryan.com/cccp/, and code for generalized helical parameterization is available upon request.     

A point mutation in the N-terminal coiled-coil domain releases c-Fes tyrosine kinase activity and survival signaling in myeloid leukemia cells

The c-fes locus encodes a 93-kDa non-receptor protein tyrosine kinase (Fes) that regulates the growth and differentiation of hematopoietic and vascular endothelial cells. Unique to Fes is a long N-terminal sequence with two regions of strong homology to coiled-coil oligomerization domains. We introduced leucine-to-proline substitutions into the coiled coils that were predicted to disrupt the coiled-coil structure. The resulting mutant proteins, together with wild-type Fes, were fused to green fluorescent protein and expressed in Rat-2 fibroblasts. We observed that a point mutation in the first coiled-coil domain (L145P) dramatically increased Fes tyrosine kinase and transforming activities in this cell type. In contrast, a similar point mutation in the second coiled-coil motif (L334P) was without effect. However, combining the L334P and L145P mutations reduced transforming and kinase activities by approximately 50% relative to the levels of activity produced with the L145P mutation alone. To study the effects of the coiled-coil mutations in a biologically relevant context, we expressed the mutant proteins in the granulocyte-macrophage colony-stimulating factor (GM-CSF)-dependent myeloid leukemia cell line TF-1. In this cellular context, the L145P mutation induced GM-CSF independence, cell attachment, and spreading. These effects correlated with a marked increase in L145P protein autophosphorylation relative to that of wild-type Fes. In contrast, the double coiled-coil mutant protein showed greatly reduced kinase and biological activities in TF-1 cells. These data are consistent with a role for the first coiled coil in the negative regulation of kinase activity and a requirement for the second coiled coil in either oligomerization or recruitment of signaling partners. Gel filtration experiments showed that the unique N-terminal region interconverts between monomeric and oligomeric forms. Single point mutations favored oligomerization, while the double point mutant protein eluted essentially as the monomer. These data provide new evidence for coiled-coil-mediated regulation of c-Fes tyrosine kinase activity and signaling, a mechanism unique among tyrosine kinases.     

Solution structure of a C-terminal coiled-coil domain from bovine IF(1): the inhibitor protein of F(1) ATPase

Bovine IF(1) is a basic, 84 amino acid residue protein that inhibits the hydrolytic action of the F(1)F(0) ATP synthase in mitochondria under anaerobic conditions. Its oligomerization state is dependent on pH. At a pH value below 6.5 it forms an active dimer. At higher pH values, two dimers associate to form an inactive tetramer. Here, we present the solution structure of a C-terminal fragment of IF(1) (44-84) containing all five of the histidine residues present in the sequence. Most unusually, the molecule forms an anti-parallel coiled-coil in which three of the five histidine residues occupy key positions at the dimer interface.     

Improving coiled-coil stability by optimizing ionic interactions

Alpha-helical coiled coils are a common protein oligomerization motif stabilized mainly by hydrophobic interactions occurring along the coiled-coil interface. We have recently designed and solved the structure of a two-heptad repeat coiled-coil peptide that is stabilized further by a complex network of inter- and intrahelical salt-bridges in addition to the hydrophobic interactions. Here, we extend and improve the de novo design of this two heptad-repeat peptide by four newly designed peptides characterized by different types of ionic interactions. The contribution of these different types of ionic interactions to coiled-coil stability are analyzed by CD spectroscopy and analytical ultracentrifugation. We show that all peptides are highly alpha-helical and two of them are 100% dimeric under physiological conditions. Furthermore, we have solved the X-ray structure of the most stable of these peptides and the rational design principles are verified by comparing this structure to the structure of the parent peptide. We show that by combining the most favorable inter- and intrahelical salt-bridge arrangements it is possible to design coiled-coil oligomerization domains with improved stability properties.     

Crystallization and Preliminary X-Ray Diffraction Analysis of the 190-Å-Long Coiled-Coil Dimerization Domain of the Actin-Bundling Protein Cortexillin I fromDictyostelium discoideum

We have crystallized the ≈190-Å-long parallel two-stranded coiled-coil oligomerization domain of the actin-bundling protein cortexillin I fromDictyostelium discoideum. The orthorhombic crystals belong to the space group C222₁with unit cell dimensions ofa= 71.3 Å,b= 127.8 Å, andc= 91.6 Å. As both native and selenomethionine-substituted protein crystals diffract to 3.0 and 2.85 Å resolution, respectively, using synchrotron radiation, they are suitable for the first high-resolution structural analysis of a two-stranded coiled coil comprising more than six heptad repeats. Moreover, because the polypeptide chain fragment contains a recently identified two-heptad-repeat long sequence that is indispensable for the assembly of the cortexillin I coiled-coil oligomerization domain, its high-resolution structure should enable us to extend our knowledge on the molecular mechanisms underlaying coiled-coil formation and to establish the precise manner in which the two “trigger” sequences interact with one another in the dimer.     

Pitch diversity in α-helical coiled coils

Two complementary methods for measuring local pitch based on heptad position in α-helical coiled coils are described and applied to six crystal structures. The results reveal a diversity of pitch values: two-stranded coiled coils appear to have pitch values near 150 Å the values for three- and four-stranded coiled coils range closer to 200 Å. The methods also provide a rapid and sensitive gauge of local coiled-coil conformation. Polar or charged residues in the apolar interface between coiled-coil helices markedly affect local pitch values, suggesting a connection between pitch uniformity and coiled-coil stability. Moreover, the identification of a skip residue (heptad frame shift) in the hemaglutinin glycoprotein of influenza virus (HA) allows interpretation of local pitch changes. These results on relatively short coiled-coil structures have relevance for the much longer fibrous proteins (many of which have skip residues) whose detailed structures are not yet established. We also show that local pitch values from molecular dynamics predictions of the GCN4 leucine zipper are in striking agreement with the high-resolution crystal structure—a result not readily discerned by direct comparison of atomic coordinates. Taken together, these methods reveal specific aspects of coiled-coil structure which may escape detection by global analyses of pitch. © 1993 Wiley-Liss, Inc.     

Predicting helix orientation for coiled-coil dimers

The alpha-helical coiled coil is a structurally simple protein oligomerization or interaction motif consisting of two or more alpha helices twisted into a supercoiled bundle. Coiled coils can differ in their stoichiometry, helix orientation, and axial alignment. Because of the near degeneracy of many of these variants, coiled coils pose a challenge to fold recognition methods for structure prediction. Whereas distinctions between some protein folds can be discriminated on the basis of hydrophobic/polar patterning or secondary structure propensities, the sequence differences that encode important details of coiled-coil structure can be subtle. This is emblematic of a larger problem in the field of protein structure and interaction prediction: that of establishing specificity between closely similar structures. We tested the behavior of different computational models on the problem of recognizing the correct orientation--parallel vs. antiparallel--of pairs of alpha helices that can form a dimeric coiled coil. For each of 131 examples of known structure, we constructed a large number of both parallel and antiparallel structural models and used these to assess the ability of five energy functions to recognize the correct fold. We also developed and tested three sequence-based approaches that make use of varying degrees of implicit structural information. The best structural methods performed similarly to the best sequence methods, correctly categorizing approximately 81% of dimers. Steric compatibility with the fold was important for some coiled coils we investigated. For many examples, the correct orientation was determined by smaller energy differences between parallel and antiparallel structures distributed over many residues and energy components. Prediction methods that used structure but incorporated varying approximations and assumptions showed quite different behaviors when used to investigate energetic contributions to orientation preference. Sequence based methods were sensitive to the choice of residue-pair interactions scored.     

Orientation and oligomerization specificity of the Bcr coiled-coil oligomerization domain

The Bcr oligomerization domain, from the Bcr-Abl oncoprotein, is an attractive therapeutic target for treating leukemias because it is required for cellular transformation. The domain homodimerizes via an antiparallel coiled coil with an adjacent short, helical swap domain. Inspection of the coiled-coil sequence does not reveal obvious determinants of helix-orientation specificity, raising the possibility that the antiparallel orientation preference and/or the dimeric oligomerization state are due to interactions of the swap domains. To better understand how structural specificity is encoded in Bcr, coiled-coil constructs containing either an N- or C-terminal cysteine were synthesized without the swap domain. When cross-linked to adopt exclusively parallel or antiparallel orientations, these showed similar circular dichroism spectra. Both constructs formed coiled-coil dimers, but the antiparallel construct was approximately 16 degrees C more stable than the parallel to thermal denaturation. Equilibrium disulfide-exchange studies confirmed that the isolated coiled-coil homodimer shows a very strong preference for the antiparallel orientation. We conclude that the orientation and oligomerization preferences of Bcr are not caused by the presence of the swap domains, but rather are directly encoded in the coiled-coil sequence. We further explored possible determinants of structural specificity by mutating residues in the d position of the coiled-coil core. Some of the mutations caused a change in orientation specificity, and all of the mutations led to the formation of higher-order oligomers. In the absence of the swap domain, these residues play an important role in disfavoring alternate states and are especially important for encoding dimeric oligomerization specificity.