An autonomous folding unit mediates the assembly of two-stranded coiled coils

Subunit oligomerization of many proteins is mediated by coiled-coil domains. Although the basic features contributing to the thermodynamic stability of coiled coils are well understood, the mechanistic details of their assembly have not yet been dissected. Here we report a 13-residue sequence pattern that occurs with limited sequence variations in many two-stranded coiled coils and that is absolutely required for the assembly of the Dictyostelium discoideum actin-bundling protein cortexillin I and the yeast transcriptional activator GCN4. The functional relationship between coiled-coil “trigger” sequences was manifested by replacing the intrinsic trigger motif of GCN4 with the related sequence from cortexillin I. We demonstrate that these trigger sequences represent autonomous helical folding units that, in contrast to arbitrarily chosen heptad repeats, can mediate coiled-coil formation. Aside from being of general interest for protein folding, trigger motifs should be of particular importance in the protein de novo design.

A computationally directed screen identifying interacting coiled coils from Saccharomyces cerevisiae

Computational methods can frequently identify protein-interaction motifs in otherwise uncharacterized open reading frames. However, the identification of candidate ligands for these motifs (e.g., so that partnering can be determined experimentally in a directed manner) is often beyond the scope of current computational capabilities. One exception is provided by the coiled-coil interaction motif, which consists of two or more α helices that wrap around each other: the ligands for coiled-coil sequences are generally other coiled-coil sequences, thereby greatly simplifying the motif/ligand recognition problem. Here, we describe a two-step approach to identifying protein–protein interactions mediated by two-stranded coiled coils that occur in Saccharomyces cerevisiae. Coiled coils from the yeast genome are first predicted computationally, by using the multicoil program, and associations between coiled coils are then determined experimentally by using the yeast two-hybrid assay. We report 213 unique interactions between 162 putative coiled-coil sequences. We evaluate the resulting interactions, focusing on associations identified between components of the spindle pole body (the yeast centrosome).

Predicting coiled coils by use of pairwise residue correlations.

A method is presented that predicts coiled-coil domains in protein sequences by using pairwise residue correlations obtained from a (two-stranded) coiled-coil database of 58,217 amino acid residues. A program called PAIRCOIL implements this method and is significantly better than existing methods at distinguishing coiled coils from alpha-helices that are not coiled coils. The database of pairwise residue correlations suggests structural features that stabilize or destabilize coiled coils.

Repacking protein cores with backbone freedom: structure prediction for coiled coils.

Progress in homology modeling and protein design has generated considerable interest in methods for predicting side-chain packing in the hydrophobic cores of proteins. Present techniques are not practically useful, however, because they are unable to model protein main-chain flexibility. Parameterization of backbone motions may represent a general and efficient method to incorporate backbone relaxation into such fixed main-chain models. To test this notion, we introduce a method for treating explicitly the backbone motions of alpha-helical bundles based on an algebraic parameterization proposed by Francis Crick in 1953 [Crick, F. H. C. (1953) Acta Crystallogr. 6, 685-689]. Given only the core amino acid sequence, a simple calculation can rapidly reproduce the crystallographic main-chain and core side-chain structures of three coiled coils (one dimer, one trimer, and one tetramer) to within 0.6-A root-mean-square deviations. The speed of the predictive method [approximately 3 min per rotamer choice on a Silicon Graphics (Mountain View, CA) 4D/35 computer] permits it to be used as a design tool.

Sequence conservation in Plasmodium falciparum alpha-helical coiled coil domains proposed for vaccine development.

BACKGROUND: The availability of the P. falciparum genome has led to novel ways to identify potential vaccine candidates. A new approach for antigen discovery based on the bioinformatic selection of heptad repeat motifs corresponding to alpha-helical coiled coil structures yielded promising results. To elucidate the question about the relationship between the coiled coil motifs and their sequence conservation, we have assessed the extent of polymorphism in putative alpha-helical coiled coil domains in culture strains, in natural populations and in the single nucleotide polymorphism data available at PlasmoDB. METHODOLOGY/PRINCIPAL FINDINGS: 14 alpha-helical coiled coil domains were selected based on preclinical experimental evaluation. They were tested by PCR amplification and sequencing of different P. falciparum culture strains and field isolates. We found that only 3 out of 14 alpha-helical coiled coils showed point mutations and/or length polymorphisms. Based on promising immunological results 5 of these peptides were selected for further analysis. Direct sequencing of field samples from Papua New Guinea and Tanzania showed that 3 out of these 5 peptides were completely conserved. An in silico analysis of polymorphism was performed for all 166 putative alpha-helical coiled coil domains originally identified in the P. falciparum genome. We found that 82% (137/166) of these peptides were conserved, and for one peptide only the detected SNPs decreased substantially the probability score for alpha-helical coiled coil formation. More SNPs were found in arrays of almost perfect tandem repeats. In summary, the coiled coil structure prediction was rarely modified by SNPs. The analysis revealed a number of peptides with strictly conserved alpha-helical coiled coil motifs. CONCLUSION/SIGNIFICANCE: We conclude that the selection of alpha-helical coiled coil structural motifs is a valuable approach to identify potential vaccine targets showing a high degree of conservation.

Controlling leucine zipper specificity with interfacial hydrophobic residues

Dimerisation of leucine zippers results from the parallel association of alpha-helices to form a coiled coil. Coiled coils comprise a heptad repeat, denoted as (abcdefg)(n), where residues at positions a and d are hydrophobic and constitute the core of the dimer interface. Charged amino acids at the e and g positions of the coiled coil are thought to be the major influence on dimerisation specificity through the formation of attractive and repulsive interhelical electrostatic interactions. However, the variability of a-position residues in leucine zipper transcription factors prompted us to investigate their influence on dimerisation specificity. We demonstrate that mutation of a single interfacial a-position Ala residue to either Val, Ile or Leu significantly alters the homo- and heterodimerisation specificities of the leucine zipper domain from the c-Jun transcription factor. These results illustrate the importance of a-position residues in controlling leucine zipper dimerisation specificity in addition to providing substantial contributions to dimer stability.

Deciphering the role of the electrostatic interactions in the alpha-tropomyosin head-to-tail complex

Skeletal alpha-tropomyosin (Tm) is a dimeric coiled-coil protein that forms linear assemblies under low ionic strength conditions in vitro through head-to-tail interactions. A previously published NMR structure of the Tin head-to-tail complex revealed that it is formed by the insertion of the N-terminal coiled-coil of one molecule into a cleft formed by the separation of the helices at the C-terminus of a second molecule. To evaluate the contribution of charged residues to complex stability, we employed single and double-mutant Tm fragments in which specific charged residues were changed to alanine in head-to-tail binding assays, and the effects of the mutations were analyzed by thermodynamic double-mutant cycles and protein-protein docking. The results show that residues K5, K7, and D280 are essential to the stability of the complex. Though D2, K6, D275, and H276 are exposed to the solvent and do not participate in intermolecular contacts in the NMR structure, they may contribute to head-to-tail complex stability by modulating the stability of the helices at the Tm termini.

Septin C-Terminal Domain Interactions: Implications for Filament Stability and Assembly

Septins form a conserved family of filament forming GTP binding proteins found in a wide range of eukaryotic cells. They share a common structural architecture consisting of an N-terminal domain, a central GTP binding domain and a C-terminal domain, which is often predicted to adopt a coiled-coil conformation, at least in part. The crystal structure of the human SEPT2/SEPT6/SEPT7 heterocomplex has revealed the importance of the GTP binding domain in filament formation, but surprisingly no electron density was observed for the C-terminal domains and their function remains obscure. The dearth of structural information concerning the C-terminal region has motivated the present study in which the putative C-terminal domains of human SEPT2, SEPT6 and SEPT7 were expressed in E. coli and purified to homogeneity. The thermal stability and secondary structure content of the domains were studied by circular dichroism spectroscopy, and homo- and hetero-interactions were investigated by size exclusion chromatography, chemical cross-linking, analytical ultracentrifugation and surface plasmon resonance. Our results show that SEPT6-C and SEPT7-C are able to form both homo- and heterodimers with a high alpha-helical content in solution. The heterodimer is elongated and considerably more stable than the homodimers, with a K (D) of 15.8 nM. On the other hand, the homodimer SEPT2-C has a much lower affinity, with a K (D) of 4 mu M, and a moderate alpha-helical content. Our findings present the first direct experimental evidence toward better understanding the biophysical properties and coiled-coil pairings of such domains and their potential role in filament assembly and stability.

Computational methods for the analysis of protein structure and function

The vast majority of known proteins have not yet been experimentally characterized and little is known about their function. The design and implementation of computational tools can provide insight into the function of proteins based on their sequence, their structure, their evolutionary history and their association with other proteins. Knowledge of the three-dimensional (3D) structure of a protein can lead to a deep understanding of its mode of action and interaction, but currently the structures of <1% of sequences have been experimentally solved. For this reason, it became urgent to develop new methods that are able to computationally extract relevant information from protein sequence and structure. The starting point of my work has been the study of the properties of contacts between protein residues, since they constrain protein folding and characterize different protein structures. Prediction of residue contacts in proteins is an interesting problem whose solution may be useful in protein folding recognition and de novo design. The prediction of these contacts requires the study of the protein inter-residue distances related to the specific type of amino acid pair that are encoded in the so-called contact map. An interesting new way of analyzing those structures came out when network studies were introduced, with pivotal papers demonstrating that protein contact networks also exhibit small-world behavior. In order to highlight constraints for the prediction of protein contact maps and for applications in the field of protein structure prediction and/or reconstruction from experimentally determined contact maps, I studied to which extent the characteristic path length and clustering coefficient of the protein contacts network are values that reveal characteristic features of protein contact maps. Provided that residue contacts are known for a protein sequence, the major features of its 3D structure could be deduced by combining this knowledge with correctly predicted motifs of secondary structure. In the second part of my work I focused on a particular protein structural motif, the coiled-coil, known to mediate a variety of fundamental biological interactions. Coiled-coils are found in a variety of structural forms and in a wide range of proteins including, for example, small units such as leucine zippers that drive the dimerization of many transcription factors or more complex structures such as the family of viral proteins responsible for virus-host membrane fusion. The coiled-coil structural motif is estimated to account for 5-10% of the protein sequences in the various genomes. Given their biological importance, in my work I introduced a Hidden Markov Model (HMM) that exploits the evolutionary information derived from multiple sequence alignments, to predict coiled-coil regions and to discriminate coiled-coil sequences. The results indicate that the new HMM outperforms all the existing programs and can be adopted for the coiled-coil prediction and for large-scale genome annotation. Genome annotation is a key issue in modern computational biology, being the starting point towards the understanding of the complex processes involved in biological networks. The rapid growth in the number of protein sequences and structures available poses new fundamental problems that still deserve an interpretation. Nevertheless, these data are at the basis of the design of new strategies for tackling problems such as the prediction of protein structure and function. Experimental determination of the functions of all these proteins would be a hugely time-consuming and costly task and, in most instances, has not been carried out. As an example, currently, approximately only 20% of annotated proteins in the Homo sapiens genome have been experimentally characterized. A commonly adopted procedure for annotating protein sequences relies on the "inheritance through homology" based on the notion that similar sequences share similar functions and structures. This procedure consists in the assignment of sequences to a specific group of functionally related sequences which had been grouped through clustering techniques. The clustering procedure is based on suitable similarity rules, since predicting protein structure and function from sequence largely depends on the value of sequence identity. However, additional levels of complexity are due to multi-domain proteins, to proteins that share common domains but that do not necessarily share the same function, to the finding that different combinations of shared domains can lead to different biological roles. In the last part of this study I developed and validate a system that contributes to sequence annotation by taking advantage of a validated transfer through inheritance procedure of the molecular functions and of the structural templates. After a cross-genome comparison with the BLAST program, clusters were built on the basis of two stringent constraints on sequence identity and coverage of the alignment. The adopted measure explicity answers to the problem of multi-domain proteins annotation and allows a fine grain division of the whole set of proteomes used, that ensures cluster homogeneity in terms of sequence length. A high level of coverage of structure templates on the length of protein sequences within clusters ensures that multi-domain proteins when present can be templates for sequences of similar length. This annotation procedure includes the possibility of reliably transferring statistically validated functions and structures to sequences considering information available in the present data bases of molecular functions and structures.

Serodiagnosis of Echinococcus spp. infection: explorative selection of diagnostic antigens by peptide microarray

BACKGROUND: Production of native antigens for serodiagnosis of helminthic infections is laborious and hampered by batch-to-batch variation. For serodiagnosis of echinococcosis, especially cystic disease, most screening tests rely on crude or purified Echinococcus granulosus hydatid cyst fluid. To resolve limitations associated with native antigens in serological tests, the use of standardized and highly pure antigens produced by chemical synthesis offers considerable advantages, provided appropriate diagnostic sensitivity and specificity is achieved. METHODOLOGY/PRINCIPAL FINDINGS: Making use of the growing collection of genomic and proteomic data, we applied a set of bioinformatic selection criteria to a collection of protein sequences including conceptually translated nucleotide sequence data of two related tapeworms, Echinococcus multilocularis and Echinococcus granulosus. Our approach targeted alpha-helical coiled-coils and intrinsically unstructured regions of parasite proteins potentially exposed to the host immune system. From 6 proteins of E. multilocularis and 5 proteins of E. granulosus, 45 peptides between 24 and 30 amino acids in length were designed. These peptides were chemically synthesized, spotted on microarrays and screened for reactivity with sera from infected humans. Peptides reacting above the cut-off were validated in enzyme-linked immunosorbent assays (ELISA). Peptides identified failed to differentiate between E. multilocularis and E. granulosus infection. The peptide performing best reached 57% sensitivity and 94% specificity. This candidate derived from Echinococcus multilocularis antigen B8/1 and showed strong reactivity to sera from patients infected either with E. multilocularis or E. granulosus. CONCLUSIONS/SIGNIFICANCE: This study provides proof of principle for the discovery of diagnostically relevant peptides by bioinformatic selection complemented with screening on a high-throughput microarray platform. Our data showed that a single peptide cannot provide sufficient diagnostic sensitivity whereas pooling several peptide antigens improved sensitivity; thus combinations of several peptides may lead the way to new diagnostic tests that replace, or at least complement conventional immunodiagnosis of echinococcosis. Our strategy could prove useful for diagnostic developments in other pathogens.

SIRE-1, a copia/Ty1-like retroelement from soybean, encodes a retroviral envelope-like protein

The soybean genome hosts a family of several hundred, relatively homogeneous copies of a large, copia/Ty1-like retroelement designated SIRE-1. A copy of this element has been recovered from a Glycine max genomic library. DNA sequence analysis of two SIRE-1 subclones revealed that SIRE-1 contains a long, uninterrupted, ORF between the 3′ end of the pol ORF and the 3′ long terminal repeat (LTR), a region that harbors the env gene in retroviral genomes. Conceptual translation of this second ORF produces a 70-kDa protein. Computer analyses of the amino acid sequence predicted patterns of transmembrane domains, α-helices, and coiled coils strikingly similar to those found in mammalian retroviral envelope proteins. In addition, a 65-residue, proline-rich domain is characterized by a strong amino acid compositional bias virtually identical to that of the 60-amino acid, proline-rich neutralization domain of the feline leukemia virus surface protein. The assignment of SIRE-1 to the copia/Ty1 family was confirmed by comparison of the conceptual translation of its reverse transcriptase-like domain with those of other retroelements. This finding suggests the presence of a proretrovirus in a plant genome and is the strongest evidence to date for the existence of a retrovirus-like genome closely related to copia/Ty1 retrotransposons.

Identification of a Second Myosin-II in Schizosaccharomyces pombe:: Myp2p Is Conditionally Required for Cytokinesis

As in many eukaryotic cells, fission yeast cytokinesis depends on the assembly of an actin ring. We cloned myp2+, a myosin-II in Schizosaccharomyces pombe, conditionally required for cytokinesis. myp2+, the second myosin-II identified in S. pombe, does not completely overlap in function with myo2+. The catalytic domain of Myp2p is highly homologous to known myosin-IIs, and phylogenetic analysis places Myp2p in the myosin-II family. The Myp2p sequence contains well-conserved ATP- and actin-binding motifs, as well as two IQ motifs. However, the tail sequence is unusual, since it is predicted to form two long coiled-coils separated by a stretch of sequence containing 19 prolines. Disruption of myp2+ is not lethal but under nutrient limiting conditions cells lacking myp2+ function are multiseptated, elongated, and branched, indicative of a defect in cytokinesis. The presence of salt enhances these morphological defects. Additionally, Δmyp2 cells are cold sensitive in high salt, failing to form colonies at 17°C. Thus, myp2+ is required under conditions of stress, possibly linking extracellular growth conditions to efficient cytokinesis and cell growth. GFP-Myp2p localizes to a ring in the middle of late mitotic cells, consistent with a role in cytokinesis. Additionally, we constructed double mutants of Δmyp2 with temperature-sensitive mutant strains defective in cytokinesis. We observed synthetic lethal interactions between Δmyp2 and three alleles of cdc11ts, as well as more modest synthetic interactions with cdc14ts and cdc16ts, implicating myp2+ function for efficient cytokinesis under normal conditions.

An integral membrane component of coatomer-coated transport vesicles defines a family of proteins involved in budding.

We have isolated a major integral membrane protein from Golgi-derived coatomer-coated vesicles. This 24-kDa protein, p24, defines a family of integral membrane proteins with homologs present in yeast and humans. In addition to sequence similarity, all p24 family members contain a motif with the characteristic heptad repeats found in coiled coils. When the yeast p24 isoform, yp24A, is knocked out in a strain defective for vesicle fusion, a dramatic reduction in the accumulation of transport vesicles is observed. Together, these results indicate a role for this protein family in the budding of coatamer-coated and other species of coated vesicles.

Mechanical properties of interfacial films formed by lysozyme self-assembly at the air-water interface

We present the first characterization of the mechanical properties of lysozyme films formed by self-assembly at the air-water interface using the Cambridge interfacial tensiometer (CIT), an apparatus capable of subjecting protein films to a much higher level of extensional strain than traditional dilatational techniques. CIT analysis, which is insensitive to surface pressure, provides a direct measure of the extensional stress-strain behavior of an interfacial film without the need to assume a mechanical model (e.g., viscoelastic), and without requiring difficult-to-test assumptions regarding low-strain material linearity. This testing method has revealed that the bulk solution pH from which assembly of an interfacial lysozyme film occurs influences the mechanical properties of the film more significantly than is suggested by the observed differences in elastic moduli or surface pressure. We have also identified a previously undescribed pH dependency in the effect of solution ionic strength on the mechanical strength of the lysozyme films formed at the air-water interface. Increasing solution ionic strength was found to increase lysozyme film strength when assembly occurred at pH 7, but it caused a decrease in film strength at pH 11, close to the pI of lysozyme. This result is discussed in terms of the significant contribution made to protein film strength by both electrostatic interactions and the hydrophobic effect. Washout experiments to remove protein from the bulk phase have shown that a small percentage of the interfacially adsorbed lysozyme molecules are reversibly adsorbed. Finally, the washout tests have probed the role played by additional adsorption to the fresh interface formed by the application of a large strain to the lysozyme film and have suggested the movement of reversibly bound lysozyme molecules from a subinterfacial layer to the interface.