Structure of a novel 13 nm dodecahedral nanocage assembled from a redesigned bacterial microcompartment shell protein.

We report the crystal structure of a novel 60-subunit dodecahedral cage that results from self-assembly of a re-engineered version of a natural protein (PduA) from the Pdu microcompartment shell. Biophysical data illustrate the dependence of assembly on solution conditions, opening up new applications in microcompartment studies and nanotechnology.

Self-assembly in the carboxysome: a viral capsid-like protein shell in bacterial cells.

Many proteins self-assemble to form large supramolecular complexes. Numerous examples of these structures have been characterized, ranging from spherical viruses to tubular protein assemblies. Some new kinds of supramolecular structures are just coming to light, while it is likely there are others that have not yet been discovered. The carboxysome is a subcellular structure that has been known for more than 40 years, but whose structural and functional details are just now emerging. This giant polyhedral body is constructed as a closed shell assembled from several thousand protein subunits. Within this protein shell, the carboxysome encapsulates the CO(2)-fixing enzymes, Rubisco (ribulose-1,5-bisphosphate carboxylase/oxygenase) and carbonic anhydrase; this arrangement enhances the efficiency of cellular CO(2) fixation. The carboxysome is present in many photosynthetic and chemoautotrophic bacteria, and so plays an important role in the global carbon cycle. It also serves as the prototypical member of what appears to be a large class of primitive protein-based organelles in bacteria. A series of crystal structures is beginning to reveal the secrets of how the carboxysome is assembled and how it enhances the efficiency of CO(2) fixation. Some of the assembly principles revealed in the carboxysome are reminiscent of those seen in icosahedral viral capsids. In addition, the shell appears to be perforated by pores for metabolite transport into and out of the carboxysome, suggesting comparisons to the pores through oligomeric transmembrane proteins, which serve to transport small molecules across the membrane bilayers of cells and eukaryotic organelles.

Crystal structure of a protein repair methyltransferase from Pyrococcus furiosus with its L-isoaspartyl peptide substrate.

Protein L-isoaspartyl (D-aspartyl) methyltransferases (EC 2.1.1.77) are found in almost all organisms. These enzymes catalyze the S-adenosylmethionine (AdoMet)-dependent methylation of isomerized and racemized aspartyl residues in age-damaged proteins as part of an essential protein repair process. Here, we report crystal structures of the repair methyltransferase at resolutions up to 1.2 A from the hyperthermophilic archaeon Pyrococcus furiosus. Refined structures include binary complexes with the active cofactor AdoMet, its reaction product S-adenosylhomocysteine (AdoHcy), and adenosine. The enzyme places the methyl-donating cofactor in a deep, electrostatically negative pocket that is shielded from solvent. Across the multiple crystal structures visualized, the presence or absence of the methyl group on the cofactor correlates with a significant conformational change in the enzyme in a loop bordering the active site, suggesting a role for motion in catalysis or cofactor exchange. We also report the structure of a ternary complex of the enzyme with adenosine and the methyl-accepting polypeptide substrate VYP(L-isoAsp)HA at 2.1 A. The substrate binds in a narrow active site cleft with three of its residues in an extended conformation, suggesting that damaged proteins may be locally denatured during the repair process in cells. Manual and computer-based docking studies on different isomers help explain how the enzyme uses steric effects to make the critical distinction between normal L-aspartyl and age-damaged L-isoaspartyl and D-aspartyl residues.

Identification of a subunit interface in transthyretin amyloid fibrils: evidence for self-assembly from oligomeric building blocks.

Amyloid and prion diseases appear to stem from the conversion of normally folded proteins into insoluble, fiber-like assemblies. Despite numerous structural studies, a detailed molecular characterization of amyloid fibrils remains elusive. In particular, models of amyloid fibrils proposed thus far have not adequately defined the constituent protein subunit interactions. To further our understanding of amyloid structure, we employed thiol-specific cross-linking and site-directed spin labeling to identify specific protein-protein associations in transthyretin (TTR) amyloid fibrils. We find that certain cysteine mutants of TTR, when dimerized by chemical cross-linkers, still form fibers under typical in vitro fibrillogenic conditions. In addition, site-directed spin labeling of many residues at the natural dimer interface reveals that their spatial proximity is preserved in the fibrillar state even in the absence of cross-linking constraints. Here, we present the first view of a subunit interface in TTR fibers and show that it is very similar to one of the natural dimeric interchain associations evident in the structure of soluble TTR. The results clarify varied models of amyloidogenesis by demonstrating that transthyretin amyloid fibrils may assemble from oligomeric protein building blocks rather than structurally rearranged monomers.

Structures of cytochrome c-549 and cytochrome c6 from the cyanobacterium Arthrospira maxima.

Cytochrome c(6) and cytochrome c-549 are small (89 and 130 amino acids, respectively) monoheme cytochromes that function in photosynthesis. They appear to have descended relatively recently from the same ancestral gene but have diverged to carry out very different functional roles, underscored by the large difference between their midpoint potentials of nearly 600 mV. We have determined the X-ray crystal structures of both proteins isolated from the cyanobacterium Arthrospira maxima. The two structures are remarkably similar, superimposing on backbone atoms with an rmsd of 0.7 A. Comparison of the two structures suggests that differences in solvent exposure of the heme and the electrostatic environment of the heme propionates, as well as in heme iron ligation, are the main determinants of midpoint potential in the two proteins. In addition, the crystal packing of both A. maxima cytochrome c-549 and cytochrome c(6) suggests that the proteins oligomerize. Finally, the cytochrome c-549 dimer we observe can be readily fit into the recently described model of cyanobacterial photosystem II.

Nanohedra: using symmetry to design self assembling protein cages, layers, crystals, and filaments.

A general strategy is described for designing proteins that self assemble into large symmetrical nanomaterials, including molecular cages, filaments, layers, and porous materials. In this strategy, one molecule of protein A, which naturally forms a self-assembling oligomer, A(n), is fused rigidly to one molecule of protein B, which forms another self-assembling oligomer, B(m). The result is a fusion protein, A-B, which self assembles with other identical copies of itself into a designed nanohedral particle or material, (A-B)(p). The strategy is demonstrated through the design, production, and characterization of two fusion proteins: a 49-kDa protein designed to assemble into a cage approximately 15 nm across, and a 44-kDa protein designed to assemble into long filaments approximately 4 nm wide. The strategy opens a way to create a wide variety of potentially useful protein-based materials, some of which share similar features with natural biological assemblies.

Evaluation of phase accuracy via topological and geometrical analysis of electron-density maps.

An empirical function is developed to measure the protein-like character of electron-density maps. The function is based upon a systematic analysis of numerous local and global map properties or descriptors. Local descriptors measure the occurrence throughout the unit cell of unique patterns on various defined templates, while global descriptors enumerate topological characteristics that define the connectivity and complexity of electron-density isosurfaces. We examine how these quantitative descriptors vary as error is introduced into the phase sets used to generate maps. Informative descriptors are combined in an optimal fashion to arrive at a predictive function. When the topological and geometrical analysis is applied to protein maps generated from phase sets with varying amounts of error, the function is able to estimate changes in average phase error with an accuracy of better than 10 degrees. Additionally, when used to monitor maps generated with experimental phases from different heavy-atom models, the analysis clearly distinguishes between the correct heavy-atom substructure solution and incorrect heavy-atom solutions. The function is also evaluated as a tool to monitor changes in map quality and phase error before and after density-modification procedures.

The crystal structure of adenylosuccinate lyase from Pyrobaculum aerophilum reveals an intracellular protein with three disulfide bonds.

Adenylosuccinate lyase catalyzes two separate reactions in the de novo purine biosynthetic pathway. Through its dual action in this pathway, adenylosuccinate lyase plays an integral part in cellular replication and metabolism. Mutations in the human enzyme can result in severe neurological disorders, including mental retardation with autistic features. The crystal structure of adenylosuccinate lyase from the hyperthermophilic archaebacterium Pyrobaculum aerophilum has been determined to 2.1 A resolution. Although both the fold of the monomer and the architecture of the tetrameric assembly are similar to adenylosuccinate lyase from the thermophilic eubacterium Thermotoga maritima, the archaebacterial lyase contains unique features. Surprisingly, the structure of adenylosuccinate lyase from P. aerophilum reveals that this intracellular protein contains three disulfide bonds that contribute significantly to its stability against thermal and chemical denaturation. The observation of multiple disulfide bonds in the recombinant form of the enzyme suggests the need for further investigations into whether the intracellular environment of P. aerophilum, and possibly other hyperthermophiles, may be compatible with protein disulfide bond formation. In addition, the protein is shorter in P. aerophilum than it is in other organisms. This abbreviation results from an internal excision of a cluster of helices that may be involved in protein-protein interactions in other organisms and may relate to the observed clinical effects of human mutations in that region.

Protein function in the post-genomic era.

Faced with the avalanche of genomic sequences and data on messenger RNA expression, biological scientists are confronting a frightening prospect: piles of information but only flakes of knowledge. How can the thousands of sequences being determined and deposited, and the thousands of expression profiles being generated by the new array methods, be synthesized into useful knowledge? What form will this knowledge take? These are questions being addressed by scientists in the field known as 'functional genomics'.

The structure of adenylosuccinate lyase, an enzyme with dual activity in the de novo purine biosynthetic pathway.

BACKGROUND: Adenylosuccinate lyase is an enzyme that plays a critical role in both cellular replication and metabolism via its action in the de novo purine biosynthetic pathway. Adenylosuccinate lyase is the only enzyme in this pathway to catalyze two separate reactions, enabling it to participate in the addition of a nitrogen at two different positions in adenosine monophosphate. Both reactions catalyzed by adenylosuccinate lyase involve the beta-elimination of fumarate. Enzymes that catalyze this type of reaction belong to a superfamily, the members of which are homotetramers. Because adenylosuccinate lyase plays an integral part in maintaining proper cellular metabolism, mutations in the human enzyme can have severe clinical consequences, including mental retardation with autistic features. RESULTS: The 1.8 A crystal structure of adenylosuccinate lyase from Thermotoga maritima has been determined by multiwavelength anomalous dispersion using the selenomethionine-substituted enzyme. The fold of the monomer is reminiscent of other members of the beta-elimination superfamily. However, its active tetrameric form exhibits striking differences in active-site architecture and cleft size. CONCLUSIONS: This first structure of an adenylosuccinate lyase reveals that, along with the catalytic base (His141) and the catalytic acid (His68), Gln212 and Asn270 might play a vital role in catalysis by properly orienting the succinyl moiety of the substrates. We propose a model for the dual activity of adenylosuccinate lyase: a single 180 degrees bond rotation must occur in the substrate between the first and second enzymatic reactions. Modeling of the pathogenic human S413P mutation indicates that the mutation destabilizes the enzyme by disrupting the C-terminal extension.

Searching for frameshift evolutionary relationships between protein sequence families.

The protein sequence database was analyzed for evidence that some distinct sequence families might be distantly related in evolution by changes in frame of translation. Sequences were compared using special amino acid substitution matrices for the alternate frames of translation. The statistical significance of alignment scores were computed in the true database and shuffled versions of the database that preserve any potential codon bias. The comparison of results from these two databases provides a very sensitive method for detecting remote relationships. We find a weak but measurable relatedness within the database as a whole, supporting the notion that some proteins may have evolved from others through changes in frame of translation. We also quantify residual homology in the ordinary sense within a database of generally unrelated sequences.

A combined algorithm for genome-wide prediction of protein function.

The availability of over 20 fully sequenced genomes has driven the development of new methods to find protein function and interactions. Here we group proteins by correlated evolution, correlated messenger RNA expression patterns and patterns of domain fusion to determine functional relationships among the 6,217 proteins of the yeast Saccharomyces cerevisiae. Using these methods, we discover over 93,000 pairwise links between functionally related yeast proteins. Links between characterized and uncharacterized proteins allow a general function to be assigned to more than half of the 2,557 previously uncharacterized yeast proteins. Examples of functional links are given for a protein family of previously unknown function, a protein whose human homologues are implicated in colon cancer and the yeast prion Sup35.

A census of protein repeats.

In this study, we analyzed all known protein sequences for repeating amino acid segments. Although duplicated sequence segments occur in 14 % of all proteins, eukaryotic proteins are three times more likely to have internal repeats than prokaryotic proteins. After clustering the repetitive sequence segments into families, we find repeats from eukaryotic proteins have little similarity with prokaryotic repeats, suggesting most repeats arose after the prokaryotic and eukaryotic lineages diverged. Consequently, protein classes with the highest incidence of repetitive sequences perform functions unique to eukaryotes. The frequency distribution of the repeating units shows only weak length dependence, implicating recombination rather than duplex melting or DNA hairpin formation as the limiting mechanism underlying repeat formation. The mechanism favors additional repeats once an initial duplication has been incorporated. Finally, we show that repetitive sequences are favored that contain small and relatively water-soluble residues. We propose that error-prone repeat expansion allows repetitive proteins to evolve more quickly than non-repeat-containing proteins.

Detecting protein function and protein-protein interactions from genome sequences.

A computational method is proposed for inferring protein interactions from genome sequences on the basis of the observation that some pairs of interacting proteins have homologs in another organism fused into a single protein chain. Searching sequences from many genomes revealed 6809 such putative protein-protein interactions in Escherichia coli and 45,502 in yeast. Many members of these pairs were confirmed as functionally related; computational filtering further enriches for interactions. Some proteins have links to several other proteins; these coupled links appear to represent functional interactions such as complexes or pathways. Experimentally confirmed interacting pairs are documented in a Database of Interacting Proteins.

A fast algorithm for genome-wide analysis of proteins with repeated sequences.

We present a fast algorithm to search for repeating fragments within protein sequences. The technique is based on an extension of the Smith-Waterman algorithm that allows the calculation of sub-optimal alignments of a sequence against itself. We are able to estimate the statistical significance of all sub-optimal alignment scores. We also rapidly determine the length of the repeating fragment and the number of times it is found in a sequence. The technique is applied to sequences in the Swissprot database, and to 16 complete genomes. We find that eukaryotic proteins contain more internal repeats than those of prokaryotic and archael organisms. The finding that 18% of yeast sequences and 28% of the known human sequences contain detectable repeats emphasizes the importance of internal duplication in protein evolution.

Protein crystals and their evil twins.

Different types of crystal twinning are reviewed with an emphasis on how to detect the phenomenon from protein diffraction data. The recent literature and a database survey both serve as reminders to perform routine checks whenever twinning is a possibility.

Assigning protein functions by comparative genome analysis: protein phylogenetic profiles.

Determining protein functions from genomic sequences is a central goal of bioinformatics. We present a method based on the assumption that proteins that function together in a pathway or structural complex are likely to evolve in a correlated fashion. During evolution, all such functionally linked proteins tend to be either preserved or eliminated in a new species. We describe this property of correlated evolution by characterizing each protein by its phylogenetic profile, a string that encodes the presence or absence of a protein in every known genome. We show that proteins having matching or similar profiles strongly tend to be functionally linked. This method of phylogenetic profiling allows us to predict the function of uncharacterized proteins.

Crystal structure and possible dimerization of the high-potential iron-sulfur protein from Chromatium purpuratum.

The crystal structure of the high-potential iron-sulfur protein (HiPIP) isolated from Chromatium purpuratum is reported at 2.7 A resolution. The three HiPIP molecules in the asymmetric unit of the crystals form one and one-half dimers. Two molecules are related by a noncrystallographic symmetry rotation of approximately 175 degrees with negligible translation along the dyad axis. The third molecule in the asymmetric unit also forms a dimer with a second HiPIP molecule across the crystallographic 2-fold symmetry axis. The Fe4S4 clusters in both the crystallographic and noncrystallographic dimers are separated by approximately 13.0 A. Solution studies give mixed results regarding the oligomeric state of the C. purpuratum HiPIP. A comparison with crystal structures of HiPIPs from other species shows that HiPIP tends to associate rather nonspecifically about a conserved, relatively hydrophobic surface patch to form dimers.

The 1.8 A crystal structure of the ycaC gene product from Escherichia coli reveals an octameric hydrolase of unknown specificity.

BACKGROUND: The ycaC gene comprises a 621 base pair open reading frame in Escherichia coli. The ycaC gene product (ycaCgp) is uncharacterized and has no assigned function. The closest sequence homologs with an assigned function belong to a family of bacterial hydrolases that catalyze isochorismatase-like reactions, but these have only low sequence similarity to ycaCgp (approximately 20% amino acid identity). The ycaCgp was obtained and identified during crystallization trials of an unrelated E. coli protein with which it co-purified. RESULTS: The 1.8 A crystal structure of ycaCgp reveals an octameric complex comprised of two tetrameric rings. A large three-layer (alphabetaalpha) sandwich domain and a small helical domain form the folded structure of the monomeric unit. Comparisons with sequence and structure databases suggest that ycaCgp belongs to a diverse family of bacterial hydrolases. The most closely related three-dimensional structure is that of the D2 tetrameric N-carbamoylsarcosine amidohydrolase (CSHase) from an Arthrobacter species. A conspicuous cleft between two ycaCgp subunits contains several conserved residues including Cys118, which we propose to be catalytic. In the active site, a nonprolyl cis peptide bond precedes Val114 and coincides with a cis peptide bond in CSHase in a region of dissimilar sequence. The crystal structure reveals a probable error or mutation relative to the reported genomic sequence. CONCLUSIONS: Although the specific function of ycaCgp is not yet known, structural studies solidify the relationship of this protein to other hydrolases and illuminate its active site and key elements of the catalytic mechanism.

Highly constrained multiple-copy refinement of protein crystal structures.

In the course of refining atomic protein structures, one often encounters difficulty with molecules that are unusually flexible or otherwise disordered. We approach the problem by combining two relatively recent developments: simultaneous refinement of multiple protein conformations and highly constrained refinement. A constrained Langevin dynamics refinement is tested on two proteins: neurotrophin-3 and glutamine synthetase. The method produces closer agreement between the calculated and observed scattering amplitudes than standard, single-copy, Gaussian atomic displacement parameter refinement. This is accomplished without significantly increasing the number of fitting parameters in the model. These results suggest that loop motion in proteins within a crystal lattice can be extensive and that it is poorly modeled by isotropic Gaussian distributions for each atom.