Electrons, because of their strong interaction with matter, produce high-resolution diffraction patterns from tiny 3D crystals only a few hundred nanometers thick in a frozen-hydrated state. This discovery offers the prospect of facile structure determination of complex biological macromolecules, which cannot be coaxed to form crystals large enough for conventional crystallography or cannot easily be produced in sufficient quantities. Two potential obstacles stand in the way. The first is a phenomenon known as dynamical scattering, in which multiple scattering events scramble the recorded electron diffraction intensities so that they are no longer informative of the crystallized molecule. The second obstacle is the lack of a proven means of de novo phase determination, as is required if the molecule crystallized is insufficiently similar to one that has been previously determined. We show with four structures of the amyloid core of the Sup35 prion protein that, if the diffraction resolution is high enough, sufficiently accurate phases can be obtained by direct methods with the cryo-EM method microelectron diffraction (MicroED), just as in X-ray diffraction. The success of these four experiments dispels the concern that dynamical scattering is an obstacle to ab initio phasing by MicroED and suggests that structures of novel macromolecules can also be determined by direct methods.

Structural studies of amyloidogenic segments by X-ray crystallography have revealed a novel packing motif, consisting of out-of-register β sheets, which may constitute one of the toxic species in aggregation related diseases. Here we sought to determine the presence of such a motif in islet amyloid polypeptide (IAPP), whose amyloidogenic properties are associated with type 2 diabetes. We determined four new crystal structures of segments within IAPP, all forming steric zippers. Most interestingly, one of the segments in the fibril core of IAPP forms an out-of-register steric zipper. Analysis of this structure reveals several commonalities with previously solved out-of-register fibrils. Our results provide additional evidence of out-of-register β sheets as a common structural motif in amyloid aggregates.

Amyloid fibers, once exclusively associated with disease, are acquiring utility as a class of biological nanomaterials. Here we introduce a method that utilizes the atomic structures of amyloid peptides, to design materials with versatile applications. As a model application, we designed amyloid fibers capable of capturing carbon dioxide from flue gas, to address the global problem of excess anthropogenic carbon dioxide. By measuring dynamic separation of carbon dioxide from nitrogen, we show that fibers with designed amino acid sequences double the carbon dioxide binding capacity of the previously reported fiber formed by VQIVYK from Tau protein. In a second application, we designed fibers that facilitate retroviral gene transfer. By measuring lentiviral transduction, we show that designed fibers exceed the efficiency of polybrene, a commonly used enhancer of transduction. The same procedures can be adapted to the design of countless other amyloid materials with a variety of properties and uses.

It has long been known that toxins produced by Bacillus thuringiensis (Bt) are stored in the bacterial cells in crystalline form. Here we describe the structure determination of the Cry3A toxin found naturally crystallized within Bt cells. When whole Bt cells were streamed into an X-ray free-electron laser beam we found that scattering from other cell components did not obscure diffraction from the crystals. The resolution limits of the best diffraction images collected from cells were the same as from isolated crystals. The integrity of the cells at the moment of diffraction is unclear; however, given the short time (∼5 µs) between exiting the injector to intersecting with the X-ray beam, our result is a 2.9-Å-resolution structure of a crystalline protein as it exists in a living cell. The study suggests that authentic in vivo diffraction studies can produce atomic-level structural information.

ALS is a terminal disease of motor neurons that is characterized by accumulation of proteinaceous deposits in affected cells. Pathological deposition of mutated Cu/Zn superoxide dismutase (SOD1) accounts for ∼20% of the familial ALS (fALS) cases. However, understanding the molecular link between mutation and disease has been difficult, given that more than 140 different SOD1 mutants have been observed in fALS patients. In addition, the molecular origin of sporadic ALS (sALS) is unclear. By dissecting the amino acid sequence of SOD1, we identified four short segments with a high propensity for amyloid fibril formation. We find that fALS mutations in these segments do not reduce their propensity to form fibrils. The atomic structures of two fibril-forming segments from the C terminus, ¹⁰¹DSVISLS¹⁰⁷ and ¹⁴⁷GVIGIAQ¹⁵³, reveal tightly packed β-sheets with steric zipper interfaces characteristic of the amyloid state. Based on these structures, we conclude that both C-terminal segments are likely to form aggregates if available for interaction. Proline substitutions in ¹⁰¹DSVISLS¹⁰⁷ and ¹⁴⁷GVIGIAQ¹⁵³ impaired nucleation and fibril growth of full-length protein, confirming that these segments participate in aggregate formation. Our hypothesis is that improper protein maturation and incompletely folded states that render these aggregation-prone segments available for interaction offer a common molecular pathway for sALS and fALS.

New materials capable of binding carbon dioxide are essential for addressing climate change. Here, we demonstrate that amyloids, self-assembling protein fibers, are effective for selective carbon dioxide capture. Solid-state NMR proves that amyloid fibers containing alkylamine groups reversibly bind carbon dioxide via carbamate formation. Thermodynamic and kinetic capture-and-release tests show the carbamate formation rate is fast enough to capture carbon dioxide by dynamic separation, undiminished by the presence of water, in both a natural amyloid and designed amyloids having increased carbon dioxide capacity. Heating to 100 °C regenerates the material. These results demonstrate the potential of amyloid fibers for environmental carbon dioxide capture.

Although amyloid fibers are found in neurodegenerative diseases, evidence points to soluble oligomers of amyloid-forming proteins as the cytotoxic species. Here, we establish that our preparation of toxic amyloid-β1–42 (Abeta42) fibrillar oligomers (TABFOs) shares with mature amyloid fibrils the cross-β structure, in which adjacent β-sheets adhere by interpenetration of protein side chains. We study the structure and properties of TABFOs by powder X-ray diffraction, EM, circular dichroism, FTIR spectroscopy, chromatography, conformational antibodies, and celluar toxicity. In TABFOs, Abeta42 molecules stack into short protofilaments consisting of pairs of helical β-sheets that wrap around each other to form a superhelix. Wrapping results in a hole along the superhelix axis, providing insight into how Abeta may form pathogenic amyloid pores. Our model is consistent with numerous properties of Abeta42 fibrillar oligomers, including heterogenous size, ability to seed new populations of fibrillar oligomers, and fiber-like morphology.

Although aberrant protein aggregation has been conclusively linked to dozens of devastating amyloid diseases, scientists remain puzzled about the molecular features that render amyloid fibrils or small oligomers toxic. Here, we report a previously unobserved type of amyloid fibril that tests as cytotoxic: one in which the strands of the contributing β-sheets are out of register. In all amyloid fibrils previously characterized at the molecular level, only in-register β-sheets have been observed, in which each strand makes its full complement of hydrogen bonds with the strands above and below it in the fibril. In out-of-register sheets, strands are sheared relative to one another, leaving dangling hydrogen bonds. Based on this finding, we designed out-of-register β-sheet amyloid mimics, which form both cylindrin-like oligomers and fibrils, and these mimics are cytotoxic. Structural and energetic considerations suggest that out-of-register fibrils can readily convert to toxic cylindrins. We propose that out-of-register β-sheets and their related cylindrins are part of a toxic amyloid pathway, which is distinct from the more energetically favored in-register amyloid pathway.

Protein amyloid oligomers have been strongly linked to amyloid diseases and can be intermediates to amyloid fibers. β-Sheets have been identified in amyloid oligomers. However, because of their transient and highly polymorphic properties, the details of their self-association remain elusive. Here we explore oligomer structure using a model system: macrocyclic peptides. Key amyloidogenic sequences from Aβ and tau were incorporated into macrocycles, thereby restraining them to β-strands, but limiting the growth of the oligomers so they may crystallize and cannot fibrillate. We determined the atomic structures for four such oligomers, and all four reveal tetrameric interfaces in which β-sheet dimers pair together by highly complementary, dry interfaces, analogous to steric zippers found in fibers, suggesting a common structure for amyloid oligomers and fibers. In amyloid fibers, the axes of the paired sheets are either parallel or antiparallel, whereas the oligomeric interfaces display a variety of sheet-to-sheet pairing angles, offering a structural explanation for the heterogeneity of amyloid oligomers.

Prion represents a unique class of pathogens devoid of nucleic acid. The deadly diseases transmitted by it between members of one species and, in certain instances, to members of other species present a public health concern. Transmissibility and the barriers to transmission between species have been suggested to arise from the degree to which a pathological protein conformation from an individual of one species can seed a pathological conformation in another species. However, this hypothesis has never been illustrated at an atomic level. Here we present three X-ray atomic structures of the same segment from human, mouse, and hamster PrP, which is critical for forming amyloid and confers species specificity in PrP seeding experiments. The structures reveal that different sequences encode different steric zippers and suggest that the degree of dissimilarity of these zipper structures gives rise to transmission barriers in prion disease, such as those that protect humans from acquiring bovine spongiform encephalopathy (BSE) and chronic wasting disease (CWD).

Amyloid-beta (Aβ) aggregates are the main constituent of senile plaques, the histological hallmark of Alzheimer’s disease. Aβ molecules form β-sheet containing structures that assemble into a variety of polymorphic oligomers, protofibers, and fibers that exhibit a range of lifetimes and cellular toxicities. This polymorphic nature of Aβ has frustrated its biophysical characterization, its structural determination, and our understanding of its pathological mechanism. To elucidate Aβ polymorphism in atomic detail, we determined eight new microcrystal structures of fiber-forming segments of Aβ. These structures, all of short, self-complementing pairs of β-sheets termed steric zippers, reveal a variety of modes of self-association of Aβ. Combining these atomic structures with previous NMR studies allows us to propose several fiber models, offering molecular models for some of the repertoire of polydisperse structures accessible to Aβ. These structures and molecular models contribute fundamental information for understanding Aβ polymorphic nature and pathogenesis.

The amylome is the universe of proteins that are capable of forming amyloid-like fibrils. Here we investigate the factors that enable a protein to belong to the amylome. A major factor is the presence in the protein of a segment that can form a tightly complementary interface with an identical segment, which permits the formation of a steric zipper—two self-complementary beta sheets that form the spine of an amyloid fibril. Another factor is sufficient conformational freedom of the self-complementary segment to interact with other molecules. Using RNase A as a model system, we validate our fibrillogenic predictions by the 3D profile method based on the crystal structure of NNQQNY and demonstrate that a specific residue order is required for fiber formation. Our genome-wide analysis revealed that self-complementary segments are found in almost all proteins, yet not all proteins form amyloids. The implication is that chaperoning effects have evolved to constrain self-complementary segments from interaction with each other.

In their recent Perspective, Cusick et al.1 state “...curation can be error-prone and possibly of lower quality than commonly assumed.” Although we welcome rigorous scrutiny of curation efforts, Cusick et al.1 had arrived at their conclusions by misunderstanding the difference between the reliability of experimental data supporting protein interactions and the correctness of the curation process itself.

In the rare medical condition termed injection amyloidosis, extracellular fibrils of insulin are observed. We found that the segment of the insulin B-chain with sequence LVEALYL is the smallest segment that both nucleates and inhibits the fibrillation of full-length insulin in a molar ratio–dependent manner, suggesting that this segment is central to the cross-β spine of the insulin fibril. In isolation from the rest of the protein, LVEALYL forms microcrystalline aggregates with fibrillar morphology, the structure of which we determined to 1 Å resolution. The LVEALYL segments are stacked into pairs of tightly interdigitated β-sheets, each pair displaying the dry steric zipper interface typical of amyloid-like fibrils. This structure leads to a model for fibrils of human insulin consistent with electron microscopic, x-ray fiber diffraction, and biochemical studies.

As genome sequencing outstrips the rate of high-quality, low-throughput biochemical and genetic experimentation, accurate annotation of protein function becomes a bottleneck in the progress of the biomolecular sciences. Most gene products are now annotated by homology, in which an experimentally determined function is applied to a similar sequence. This procedure becomes error-prone between more divergent sequences and can contaminate biomolecular databases. Here, we propose a computational method of assignment of function, termed Generalized Functional Linkages (GFL), that combines nonhomology-based methods with other types of data. Functional linkages describe pairwise relationships between proteins that work together to perform a biological task. GFL provides a Bayesian framework that improves annotation by arbitrating a competition among biological process annotations to best describe the target protein. GFL addresses the unequal strengths of functional linkages among proteins, the quality of existing annotations, and the similarity among them while incorporating available knowledge about the cellular location or individual molecular function of the target protein. We demonstrate GFL with functional linkages defined by an algorithm known as zorch that quantifies connectivity in protein–protein interaction networks. Even when using proteins linked only by indirect or high-throughput interactions, GFL predicts the biological processes of many proteins in Saccharomyces cerevisiae, improving the accuracy of annotation by 20% over majority voting.

Amyloid fibrils formed from different proteins, each associated with a particular disease, contain a common cross-β spine. The atomic architecture of a spine, from the fibril-forming segment GNNQQNY of the yeast prion protein Sup35, was recently revealed by X-ray microcrystallography. It is a pair of β-sheets, with the facing side chains of the two sheets interdigitated in a dry ‘steric zipper’. Here we report some 30 other segments from fibril-forming proteins that form amyloid-like fibrils, microcrystals, or usually both. These include segments from the Alzheimer’s amyloid-β and tau proteins, the PrP prion protein, insulin, islet amyloid polypeptide (IAPP), lysozyme, myoglobin, α-synuclein and β2-microglobulin, suggesting that common structural features are shared by amyloid diseases at the molecular level. Structures of 13 of these microcrystals all reveal steric zippers, but with variations that expand the range of atomic architectures for amyloid-like fibrils and offer an atomic-level hypothesis for the basis of prion strains.

Identifying sequence determinants of fibril-forming proteins is crucial for understanding the processes causing >20 proteins to form pathological amyloid depositions. Our approach to identifying which sequences form amyloid-like fibrils is to screen the amyloid-forming proteins human insulin and β2-microglobulin for segments that form fibrils. Our screen is of 60 sequentially overlapping peptides, 59 being six residues in length and 1 being five residues, covering every noncysteine-containing segment in these two proteins. Each peptide was characterized as amyloid-like or nonfibril-forming. Amyloid-like peptides formed fibrils visible in electron micrographs or needle-like microcrystals showing a cross-β diffraction pattern. Eight of the 60 peptides (three from insulin and five from β2-microglobulin) were identified as amyloid-like. The results of the screen were used to assess the computational method, and good agreement between prediction and experiments was found. This agreement suggests that the pair-of-sheets, zipper spine model on which the computational method is based is at least approximately correct for the structure of the fibrils and suggests the nature of the sequence signal for formation of amyloid-like fibrils.

Based on the crystal structure of the cross-β spine formed by the peptide NNQQNY, we have developed a computational approach for identifying those segments of amyloidogenic proteins that themselves can form amyloid-like fibrils. The approach builds on experiments showing that hexapeptides are sufficient for forming amyloid-like fibrils. Each six-residue peptide of a protein of interest is mapped onto an ensemble of templates, or 3D profile, generated from the crystal structure of the peptide NNQQNY by small displacements of one of the two intermeshed β-sheets relative to the other. The energy of each mapping of a sequence to the profile is evaluated by using rosettadesign, and the lowest energy match for a given peptide to the template library is taken as the putative prediction. If the energy of the putative prediction is lower than a threshold value, a prediction of fibril formation is made. This method can reach an accuracy of ≈80% with a P value of ≈10−12 when a conservative energy threshold is used to separate peptides that form fibrils from those that do not. We see enrichment for positive predictions in a set of fibril-forming segments of amyloid proteins, and we illustrate the method with applications to proteins of interest in amyloid research.

Numerous soluble proteins convert to insoluble amyloid-like fibrils that have common properties. Amyloid fibrils are associated with fatal diseases such as Alzheimer's, and amyloid-like fibrils can be formed in vitro. For the yeast protein Sup35, conversion to amyloid-like fibrils is associated with a transmissible infection akin to that caused by mammalian prions. A seven-residue peptide segment from Sup35 forms amyloid-like fibrils and closely related microcrystals, from which we have determined the atomic structure of the cross-β spine. It is a double β-sheet, with each sheet formed from parallel segments stacked in register. Side chains protruding from the two sheets form a dry, tightly self-complementing steric zipper, bonding the sheets. Within each sheet, every segment is bound to its two neighbouring segments through stacks of both backbone and side-chain hydrogen bonds. The structure illuminates the stability of amyloid fibrils, their self-seeding characteristic and their tendency to form polymorphic structures.

In humans suffering from dialysis-related amyloidosis, the protein β2-microglobulin (β2M) is deposited as an amyloid; however, an amyloid of β2M is unknown in mice. β2M sequences from human and mouse are 70% identical, but there is a seven-residue peptide in which six residues differ. This peptide from human β2M forms amyloid in vitro, whereas the mouse peptide does not. Substitution of the human peptide for its counterpart in the mouse sequence results in the formation of amyloid in vitro. These results show that a seven-residue segment of human β2M is sufficient to convert β2M to the amyloid state, and that specific residue interactions are crucial to the conversion. These observations are consistent with a proposed Zipper-spine model for β2M amyloid, in which the spine of the fibril consists of an anhydrous β-sheet.

The infectious form of prion protein, PrPSc, self-propagates by its conversion of the normal, cellular prion protein molecule PrPC to another PrPSc molecule. It has not yet been demonstrated that recombinant prion protein can convert prion protein molecules from PrPC to PrPSc. Here we show that recombinant hamster prion protein is converted to a second form, PrPRDX, by a redox process in vitro and that this PrPRDX form seeds the conversion of other PrPC molecules to the PrPRDX form. The converted form shows properties of oligomerization and seeded conversion that are characteristic of PrPSc. We also find that the oligomerization can be reversed in vitro. X-ray fiber diffraction suggests an amyloid-like structure for the oligomerized prion protein. A domain-swapping model involving intermolecular disulfide bonds can account for the stability and coexistence of two molecular forms of prion protein and the capacity of the second form for self-propagation.

To better understand the roles of Sm proteins in forming the cores of many RNA-processing ribonucleoproteins, we determined the crystal structure of an atypical Sm-like archaeal protein (SmAP3) in which the conserved Sm domain is augmented by a previously uncharacterized, mixed α/β C-terminal domain. The structure reveals an unexpected SmAP3 14-mer that is perforated by a cylindrical pore and is bound to 14 cadmium (Cd2+) ions. Individual heptamers adopt either “apical” or “equatorial” conformations that chelate Cd2+ differently. SmAP3 forms supraheptameric oligomers (SmAP3)n = 7,14,28 in solution, and assembly of the asymmetric 14-mer is modulated by differential divalent cation-binding in apical and equatorial subunits. Phylogenetic and sequence analyses substantiate SmAP3s as a unique subset of SmAPs. These results distinguish SmAP3s from other Sm proteins and provide a model for the structure and properties of Sm proteins >100 residues in length, e.g., several human Sm proteins.

The Directional Atomic Solvation EnergY (DASEY) is an atom-based description of the environment of an amino acid position within a known 3D protein structure. The DASEY has been developed to align and score a probe amino acid sequence to a library of template protein structures for fold assignment. DASEY is computed by summing the atomic solvation parameters of atoms falling within a tetrahedral sector, or petal, extending 16 Å along each of the four bond axes of each α-carbon atom of the protein. The DASEY discriminates between pairs of structurally equivalent positions and random pairs in protein structures sharing a fold but belonging to different superfamilies, unlike some previous descriptors of protein environments, such as buried area. Furthermore, the DASEY values have characteristic patterns of residue replacement, an essential feature of a successful fold assignment method. Benchmarking fold assignment with DASEY achieves coverage of 56% of sequences with 90% accuracy when probe sequences are matched to protein structural templates belonging to the same fold but to a different superfamily, an improvement of greater than 200% over a previous method.

Sm proteins form the core of small nuclear ribonucleoprotein particles (snRNPs), making them key components of several mRNA-processing assemblies, including the spliceosome. We report the 1.75-Å crystal structure of SmAP, an Sm-like archaeal protein that forms a heptameric ring perforated by a cationic pore. In addition to providing direct evidence for such an assembly in eukaryotic snRNPs, this structure (i) shows that SmAP homodimers are structurally similar to human Sm heterodimers, (ii) supports a gene duplication model of Sm protein evolution, and (iii) offers a model of SmAP bound to single-stranded RNA (ssRNA) that explains Sm binding-site specificity. The pronounced electrostatic asymmetry of the SmAP surface imparts directionality to putative SmAP–RNA interactions.

Three-dimensional (3D) domain-swapped proteins are intermolecularly folded analogs of monomeric proteins; both are stabilized by the identical interactions, but the individual domains interact intramolecularly in monomeric proteins, whereas they form intermolecular interactions in 3D domain-swapped structures. The structures and conditions of formation of several domain-swapped dimers and trimers are known, but the formation of higher order 3D domain-swapped oligomers has been less thoroughly studied. Here we contrast the structural consequences of domain swapping from two designed three-helix bundles: one with an up-down-up topology, and the other with an up-down-down topology. The up-down-up topology gives rise to a domain-swapped dimer whose structure has been determined to 1.5 Å resolution by x-ray crystallography. In contrast, the domain-swapped protein with an up-down-down topology forms fibrils as shown by electron microscopy and dynamic light scattering. This demonstrates that design principles can predict the oligomeric state of 3D domain-swapped molecules, which should aid in the design of domain-swapped proteins and biomaterials.

X-ray diffraction and other biophysical tools reveal features of the atomic structure of an amyloid-like crystal. Sup35, a prion-like protein in yeast, forms fibrillar amyloid assemblies intrinsic to its prion function. We have identified a polar peptide from the N-terminal prion-determining domain of Sup35 that exhibits the amyloid properties of full-length Sup35, including cooperative kinetics of aggregation, fibril formation, binding of the dye Congo red, and the characteristic cross-β x-ray diffraction pattern. Microcrystals of this peptide also share the principal properties of the fibrillar amyloid, including a highly stable, β-sheet-rich structure and the binding of Congo red. The x-ray powder pattern of the microcrystals, extending to 0.9-Å resolution, yields the unit cell dimensions of the well-ordered structure. These dimensions restrict possible atomic models of this amyloid-like structure and demonstrate that it forms packed, parallel-stranded β-sheets. The unusually high density of the crystals shows that the packed β-sheets are dehydrated, despite the polar character of the side chains. These results suggest that amyloid is a highly intermolecularly bonded, dehydrated array of densely packed β-sheets. This dry β-sheet could form as Sup35 partially unfolds to expose the peptide, permitting it to hydrogen-bond to the same peptide of other Sup35 molecules. The implication is that amyloid-forming units may be short segments of proteins, exposed for interactions by partial unfolding.

d-Lactate dehydrogenase (d-LDH) of Escherichia coli is a peripheral membrane respiratory enzyme involved in electron transfer, located on the cytoplasmic side of the inner membrane. d-LDH catalyzes the oxidation of d-lactate to pyruvate, which is coupled to transmembrane transport of amino acids and sugars. Here we describe the crystal structure at 1.9 Å resolution of the three domains of d-LDH: the flavin adenine dinucleotide (FAD)-binding domain, the cap domain, and the membrane-binding domain. The FAD-binding domain contains the site of d-lactate reduction by a noncovalently bound FAD cofactor and has an overall fold similar to other members of a recently discovered FAD-containing family of proteins. This structural similarity extends to the cap domain as well. The most prominent difference between d-LDH and the other members of the FAD-containing family is the membrane-binding domain, which is either absent in some of these proteins or differs significantly. The d-LDH membrane-binding domain presents an electropositive surface with six Arg and five Lys residues, which presumably interacts with the negatively charged phospholipid head groups of the membrane. Thus, d-LDH appears to bind the membrane through electrostatic rather than hydrophobic forces.

Three-dimensional protein folds were assigned to all ORFs of the recently sequenced genome of the hyperthermophilic archaeon Pyrobaculum aerophilum. Binary hypothesis testing was used to estimate a confidence level for each assignment. A separate test was conducted to assign a probability for whether each sequence has a novel fold—i.e., one that is not yet represented in the experimental database of known structures. Of the 2,130 predicted nontransmembrane proteins in this organism, 916 matched a fold at a cumulative 90% confidence level, and 245 could be assigned at a 99% confidence level. Likewise, 286 proteins were predicted to have a previously unobserved fold with a 90% confidence level, and 14 at a 99% confidence level. These statistically based tools are combined with homology searches against the Online Mendelian Inheritance in Man (OMIM) human genetics database and other protein databases for the selection of attractive targets for crystallographic or NMR structure determination. Results of these studies have been collated and placed at http://www.doe-mbi.ucla.edu/people/parag/PA_HOME/, the University of California, Los Angeles–Department of Energy Pyrobaculum aerophilum web site.

We introduce a computational method for identifying subcellular locations of proteins from the phylogenetic distribution of the homologs of organellar proteins. This method is based on the observation that proteins localized to a given organelle by experiments tend to share a characteristic phylogenetic distribution of their homologs—a phylogenetic profile. Therefore any other protein can be localized by its phylogenetic profile. Application of this method to mitochondrial proteins reveals that nucleus-encoded proteins previously known to be destined for mitochondria fall into three groups: prokaryote-derived, eukaryote-derived, and organism-specific (i.e., found only in the organism under study). Prokaryote-derived mitochondrial proteins can be identified effectively by their phylogenetic profiles. In the yeast Saccharomyces cerevisiae, 361 nucleus-encoded mitochondrial proteins can be identified at 50% accuracy with 58% coverage. From these values and the proportion of conserved mitochondrial genes, it can be inferred that ≈630 genes, or 10% of the nuclear genome, is devoted to mitochondrial function. In the worm Caenorhabditis elegans, we estimate that there are ≈660 nucleus-encoded mitochondrial genes, or 4% of its genome, with ≈400 of these genes contributed from the prokaryotic mitochondrial ancestor. The large fraction of organism-specific and eukaryote-derived genes suggests that mitochondria perform specialized roles absent from prokaryotic mitochondrial ancestors. We observe measurably distinct phylogenetic profiles among proteins from different subcellular compartments, allowing the general use of prokaryotic genomes in learning features of eukaryotic proteins.

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 evolution1, correlated messenger RNA expression patterns2 and patterns of domain fusion3 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.

Fatty acid biosynthesis is essential for the survival of Mycobacterium tuberculosis and acetyl-coenzyme A (acetyl-CoA) is an essential precursor in this pathway. We have determined the 3-D crystal structure of M. tuberculosis citrate lyase β-subunit (CitE), which as annotated should cleave protein bound citryl-CoA to oxaloacetate and a protein-bound CoA derivative. The CitE structure has the (β/α)8 TIM barrel fold with an additional α-helix, and is trimeric. We have determined the ternary complex bound with oxaloacetate and magnesium, revealing some of the conserved residues involved in catalysis. While the bacterial citrate lyase is a complex with three subunits, the M. tuberculosis genome does not contain the α and γ subunits of this complex, implying that M. tuberculosis CitE acts differently from other bacterial CitE proteins. The analysis of gene clusters containing the CitE protein from 168 fully sequenced organisms has led us to identify a grouping of functionally related genes preserved in M. tuberculosis, Rattus norvegicus, Homo sapiens, and Mus musculus. We propose a novel enzymatic function for M. tuberculosis CitE in fatty acid biosynthesis that is analogous to bacterial citrate lyase but producing acetyl-CoA rather than a protein-bound CoA derivative.

In the study of protein complexes, is there a computational method for inferring which combinations of proteins in an organism are likely to form a crystallizable complex? Here we attempt to answer this question, using the Protein Data Bank (PDB) to assess the usefulness of inferred functional protein linkages from the Prolinks database. We find that of the 242 nonredundant prokaryotic protein complexes shared between the current PDB and Prolinks, 44% (107/242) contain proteins linked at high confidence by one or more methods of computed functional linkages. Similarly, high-confidence linkages detect 47% of known Escherichia coli protein complexes, with 45% accuracy. Together these findings suggest that functional linkages will be useful in defining protein complexes for structural studies, including for structural genomics. We offer a database of inferred linkages corresponding to likely protein complexes for some 629,952 pairs of proteins in 154 prokaryotes and archaea.