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.

Functional analysis of predicted coiled-coil regions in the Escherichia coli K-12 O-antigen polysaccharide chain length determinant Wzz

Wzz is a membrane protein that determines the chain length distribution of the O-antigen lipopolysaccharide by an unknown mechanism. Wzz proteins consist of two transmembrane helices separated by a large periplasmic loop. The periplasmic loop of Escherichia coli K-12 Wzz (244 amino acids from K65 to A308) was purified and found to be a monomer with an extended conformation, as determined by gel filtration chromatography and analytical ultracentrifugation. Circular dichroism showed that the loop has a 60% helical content. The Wzz periplasmic loop also contains three regions with predicted coiled coils. To probe the function of the predicted coiled coils, we constructed amino acid replacement mutants of the E. coli K-12 Wzz protein, which were designed so that the coiled coils could be separate without compromising the helicity of the individual molecules. Mutations in one of the regions, spanning amino acids 108 to 130 (region I), were associated with a partial defect in O-antigen chain length distribution, while mutants with mutations in the region spanning amino acids 209 to 223 (region III) did not have an apparent functional defect. In contrast, mutations in the region spanning amino acids 153 to 173 (region II) eliminated the Wzz function. This phenotype was associated with protein instability, most likely due to conformational changes caused by the amino acid replacements, which was confirmed by limited trypsin proteolysis. Additional mutagenesis based on a three-dimensional model of region I demonstrated that the amino acids implicated in function are all located at the same face of a predicted alpha-helix, suggesting that a coiled coil actually does not exist in this region. Together, our results suggest that the regions predicted to be coiled coils are important for Wzz function because they maintain the native conformation of the protein, although the existence of coiled coils could not be demonstrated experimentally.

The Structure and Functions of Hantavirus Nucleocapsid Protein

Hantaviruses, members of the genus Hantavirus in the Bunyaviridae family, are enveloped single-stranded RNA viruses with tri-segmented genome of negative polarity. In humans, hantaviruses cause two diseases, hemorrhagic fever with renal syndrome (HFRS) and hantavirus pulmonary syndrome (HPS), which vary in severity depending on the causative agent. Each hantavirus is carried by a specific rodent host and is transmitted to humans through excreta of infected rodents. The genome of hantaviruses encodes four structural proteins: the nucleocapsid protein (N), the glycoproteins (Gn and Gc), and the polymerase (L) and also the nonstructural protein (NSs). This thesis deals with the functional characterization of hantavirus N protein with regard to its structure. Structural studies of the N protein have progressed slowly and the crystal structure of the whole protein is still not available, therefore biochemical assays coupled with bioinformatical modeling proved essential for studying N protein structure and functions. Presumably, during RNA encapsidation, the N protein first forms intermediate trimers and then oligomers. First, we investigated the role of N-terminal domain in the N protein oligomerization. The results suggested that the N-terminal region of the N protein forms a coiled-coil, in which two antiparallel alpha helices interact via their hydrophobic seams. Hydrophobic residues L4, I11, L18, L25 and V32 in the first helix and L44, V51, L58 and L65 in the second helix were crucial for stabilizing the structure. The results were consistent with the head-to-head, tail-to-tail model for hantavirus N protein trimerization. We demonstrated that an intact coiled-coil structure of the N terminus is crucial for the oligomerization capacity of the N protein. We also added new details to the head-to-head, tail-to-tail model of trimerization by suggesting that the initial step is based on interaction(s) between intact intra-molecular coiled-coils of the monomers. We further analyzed the importance of charged aa residues located within the coiled-coil for the N protein oligomerization. To predict the interacting surfaces of the monomers we used an upgraded in silico model of the coiled-coil domain that was docked into a trimer. Next the predicted target residues were mutated. The results obtained using the mammalian two-hybrid assay suggested that conserved charged aa residues within the coiled-coil make a substantial contribution to the N protein oligomerization. This contribution probably involves the formation of interacting surfaces of the N monomers and also stabilization of the coiled-coil via intramolecular ionic bridging. We proposed that the tips of the coiled-coils are the first to come into direct contact and thus initiate tight packing of the three monomers into a compact structure. This was in agreement with the previous results showing that an increase in ionic strength abolished the interaction between N protein molecules. We also showed that residues having the strongest effect on the N protein oligomerization are not scattered randomly throughout the coiled-coil 3D model structure, but form clusters. Next we found evidence for the hantaviral N protein interaction with the cytoplasmic tail of the glycoprotein Gn. In order to study this interaction we used the GST pull-down assay in combination with mutagenesis technique. The results demonstrated that intact, properly folded zinc fingers of the Gn protein cytoplasmic tail as well as the middle domain of the N protein (that includes aa residues 80 248 and supposedly carries the RNA-binding domain) are essential for the interaction. Since hantaviruses do not have a matrix protein that mediates the packaging of the viral RNA in other negatve stranded viruses (NSRV), hantaviral RNPs should be involved in a direct interaction with the intraviral domains of the envelope-embedded glycoproteins. By showing the N-Gn interaction we provided the evidence for one of the crucial steps in the virus replication at which RNPs are directed to the site of the virus assembly. Finally we started analysis of the N protein RNA-binding region, which is supposedly located in the middle domain of the N protein molecule. We developed a model for the initial step of RNA-binding by the hantaviral N protein. We hypothesized that the hantaviral N protein possesses two secondary structure elements that initiate the RNA encapsidation. The results suggest that amino acid residues (172-176) presumably act as a hook to catch vRNA and that the positively charged interaction surface (aa residues 144-160) enhances the initial N-RNA interacation. In conclusion, we elucidated new functions of hantavirus N protein. Using in silico modeling we predicted the domain structure of the protein and using experimental techniques showed that each domain is responsible for executing certain function(s). We showed that intact N terminal coiled-coil domain is crucial for oligomerization and charged residues located on its surface form a interaction surface for the N monomers. The middle domain is essential for interaction with the cytoplasmic tail of the Gn protein and RNA binding.

Cell Surface Assembly of HIV gp41 Six-Helix Bundles for Facile, Quantitative Measurements of Hetero-oligomeric Interactions

Helix helix interactions are fundamental to many biological signals and systems and are found in homo- or heteromultimerization of signaling molecules as well as in the process of virus entry into the host. In HIV, virus-host membrane fusion during infection is mediated by the formation of six-helix bundles (6HBs) from homotrimers of gp41, from which a number of synthetic peptides have been derived as antagonists of virus entry. Using a yeast surface two-hybrid (YS2H) system, a platform designed to detect protein-protein interactions occurring through a secretory pathway, we reconstituted 6HB complexes on the yeast surface, quantitatively measured the equilibrium and kinetic constants of soluble 6HB, and delineated the residues influencing homo-oligomeric and hetero-oligomeric coiled-coil interactions. Hence, we present YS2H as a platform for the facile characterization and design of antagonistic peptides for inhibition of HIV and many other enveloped viruses relying on membrane fusion for infection, as well as cellular signaling events triggered by hetero-oligomeric coiled coils.

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.