PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosbiolPLoS BiologySubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
 
PLoS Biol. Dec 2004; 2(12): e380.
Published online Nov 2, 2004. doi:  10.1371/journal.pbio.0020380
PMCID: PMC524472
Components of Coated Vesicles and Nuclear Pore Complexes Share a Common Molecular Architecture
Damien Devos,1 Svetlana Dokudovskaya,2 Frank Alber,1 Rosemary Williams,2 Brian T Chait,3 Andrej Sali,corresponding author1 and Michael P Routcorresponding author2
1Departments of Biopharmaceutical Sciences and Pharmaceutical Chemistry and California Institute for Quantitative Biomedical Research, University of California, San Francisco, California, United States of America
2Laboratory of Cellular and Structural Biology, Rockefeller University, New York, New York, United States of America
3Laboratory of Mass Spectrometry and Gaseous Ion Chemistry, Rockefeller University, New York, New York, United States of America
corresponding authorCorresponding author.
Andrej Sali: sali/at/salilab.org; Michael P Rout: rout/at/rockefeller.edu
Received July 13, 2004; Accepted August 7, 2004.
Numerous features distinguish prokaryotes from eukaryotes, chief among which are the distinctive internal membrane systems of eukaryotic cells. These membrane systems form elaborate compartments and vesicular trafficking pathways, and sequester the chromatin within the nuclear envelope. The nuclear pore complex is the portal that specifically mediates macromolecular trafficking across the nuclear envelope. Although it is generally understood that these internal membrane systems evolved from specialized invaginations of the prokaryotic plasma membrane, it is not clear how the nuclear pore complex could have evolved from organisms with no analogous transport system. Here we use computational and biochemical methods to perform a structural analysis of the seven proteins comprising the yNup84/vNup107–160 subcomplex, a core building block of the nuclear pore complex. Our analysis indicates that all seven proteins contain either a β-propeller fold, an α-solenoid fold, or a distinctive arrangement of both, revealing close similarities between the structures comprising the yNup84/vNup107–160 subcomplex and those comprising the major types of vesicle coating complexes that maintain vesicular trafficking pathways. These similarities suggest a common evolutionary origin for nuclear pore complexes and coated vesicles in an early membrane-curving module that led to the formation of the internal membrane systems in modern eukaryotes.
The ability to sharply curve membranes was a defining event in the evolution of early eukaryotes, allowing the formation of endomembrane systems (Blobel 1980). In modern eukaryotes, these systems have become elaborate internal membranes, such as the Golgi apparatus, the endoplasmic reticulum (ER), and the nuclear envelope (NE). To date three major kinds of transport vesicles, distinguished by the compositions of their protein coat complexes, have been shown to traffic between these internal membranes and the plasma membrane: First, the clathrin/adaptin complexes are responsible for endocytosis and vesicular trafficking between the Golgi, lysosomes, and endosomes; second, the COPI complex mediates intra-Golgi and Golgi-to-ER trafficking; and lastly, the COPII complex supports vesicle movement from the ER to the Golgi (reviewed in Kirchhausen 2000a, 2000b; Boehm and Bonifacino 2001; Bonifacino and Lippincott-Schwartz 2003; Lippincott-Schwartz and Liu 2003).
The NE is contiguous with the ER and delineates the nucleus. It is made of an inner and outer membrane that together form a barrier between the nucleoplasm and the cytoplasm. This barrier is perforated by nuclear pore complexes (NPCs), which form pores between the inner and outer NE membranes by stabilizing a sharply curved section of connecting pore membrane. NPCs are approximately 50-MDa octagonally symmetric cylinders that function as the only known mediators of nucleocytoplasmic exchange; while permitting the free flow of small molecules, they restrict macromolecular trafficking to selected cargoes that are recognized by cognate transport factors. NPCs are found in all eukaryotic cells and are composed of a broadly conserved set of proteins, termed nups (reviewed in Rout and Aitchison 2001; Bednenko et al. 2003; Rout et al. 2003; Suntharalingam and Wente 2003; Fahrenkrog et al. 2004). Although the nups have been fully cataloged for both yeast (Saccharomyces) (Rout et al. 2000) and vertebrates (Cronshaw et al. 2002), there is currently little information concerning their origin and evolution. To this end, protein structures are helpful because it is easier to recognize similarities in structure than in sequence, especially for distantly related proteins. Thus, we have characterized the structures of seven proteins forming a core building block of the NPC, termed the yNup84 subcomplex in Saccharomyces and the vNup107–160 subcomplex in vertebrates. These structures reveal how the nuclear pore complex could have evolved from organisms with no analogous transport system.
The yNup84/vNup107–160 subcomplex has a molecular weight of approximately 600 kDa and has been shown in yeast to be flexible (Siniossoglou et al. 1996; Siniossoglou et al. 2000; Lutzmann et al. 2002), presenting a considerable challenge to conventional experimental methods for structure determination; thus, we used a computational approach that relies on a database of experimentally determined structures (Marti-Renom et al. 2000). We first focused on the component nups of the yNup84 subcomplex: ySeh1, ySec13, yNup84, yNup85, yNup120, yNup133, and yNup145C, whose corresponding vertebrate homologs are, respectively, vSec13 l, vSec13R, vNup107, vNup75, vNup160, vNup133, and vNup96 (Siniossoglou et al. 1996; Fontoura et al. 1999; Siniossoglou et al. 2000; Cronshaw et al. 2002; Lutzmann et al. 2002; Boehmer et al. 2003; Harel et al. 2003; Walther et al. 2003; Loiodice et al. 2004). For putative domains in each of these nups, we first applied two threading programs to assign structure folds based on similarity to known protein structures (templates) (Marti-Renom et al. 2000) (see Materials and Methods). The corresponding sequence-structure alignments were refined and used to generate three-dimensional models of the nup domains, followed by evaluation of the models. Our analyses predicted that every nup in the yNup84/vNup107–160 subcomplex consists of a β-propeller domain, an α-solenoid domain, or both (Figure 1; Table 1). β-propellers contain several blades arranged radially around a central axis, each blade consisting of a four-stranded antiparallel β-sheet; α-solenoid domains are composed of numerous pairs of antiparallel α-helices stacked to form a solenoid (Figure 1) (Neer et al. 1994; Andrade et al. 2001a; Andrade et al. 2001b). While we have not defined the precise details of each domain, such as the exact shapes and numbers of propeller blades and solenoid repeats, the overall fold assignments for each nup are clear. These predictions indicate that yNup84, yNup85, and yNup145C all mainly consist of an α-solenoid domain, whereas yNup120 and yNup133 contain both an amino-terminal β-propeller and a large carboxyl-terminal α-solenoid region. Both ySec13 and ySeh1 are predicted to be almost entirely single-domain β-propellers of six and seven blades, respectively. These latter two proteins fall into the well-conserved class of tryptophan/aspartic acid (WD) repeat-containing β-propeller proteins. For both proteins, homology with the WD-repeat β-propellers has been reported previously (Saxena et al. 1996; Siniossoglou et al. 1996; Yu et al. 2000) and is confirmed here.
Figure 1
Figure 1
Ribbon Representation of Nup Models
Table 1
Table 1
Nup84 Subcomplex Proteins are Composed of Two Fold Types
We support our fold assignments using four considerations (Figure 2; Tables 1 and S1–S7). First, both fold assignment programs returned their predictions with highly significant scores (Tables S1–S7), and they predominantly assigned only the two predicted folds out of the approximately 1,000 different known fold types (Tables S1–S7) (Orengo et al. 1997). Moreover, while there are numerous variations corresponding to different proteins within each predicted fold type, the two different methods used for fold recognition often selected the same template proteins (Tables S1–S7). Second, the evaluation of the atomic model for each nup was statistically significant when compared against the best models generated for random sequences of identical amino acid composition and length; all the nup models were at least six standard deviations away from the mean score of the random models (Figure S1; Tables 1 and S1–S7) (Melo et al. 2002). Third, secondary structure predictions from amino acid sequences alone indicate that all seven nups consist mainly of repetitive structures that largely match the secondary structures observed in their corresponding three-dimensional models (Figure 3 and Figure S2). The agreement ranges from 58% to 87% of the residues for a three-state assignment (helix, strand, other). This agreement is the maximum possible level of consistency, given the approximately 75% accuracy of the secondary structure prediction methods (Koh et al. 2003).
Figure 2
Figure 2
Proteolytic Domain Map of the Yeast Nup84 Subcomplex Proteins
Figure 3
Figure 3
Predicted Secondary Structure Maps of the Nup84 Subcomplex Proteins
Finally, we provide direct biochemical evidence in support of our fold assignments, using proteolytic mapping of domain boundaries and loop locations in the seven nups (see Figure 2). Tagged nups were purified from yeast extracts and incubated with the endoproteinases Asp-N (which hydrolyzes peptide bonds at the amino side of aspartic acid) or Lys-C (which hydrolyzes peptide bonds at the carboxylic side of lysines) while still attached to the magnetic beads via their proteolytically resistant tags. After digestion, proteolytic fragments that remained attached to the beads were separated by SDS-PAGE, and cleavage sites were determined either by molecular weight estimation of the fragments or by amino-terminal Edman sequencing (Table 2). The regions predicted to form β-propellers were, as expected, extremely resistant to proteolysis (see Figure 2) (Kirchhausen and Harrison 1984; Saxena et al. 1996). On the whole, the predicted α-solenoid regions were also resistant to proteolysis, although less so than the β-propellers. However, the major cleavages were found toward the end of the predicted α-solenoid domains, even in the most susceptible nup (yNup133). Strikingly, the strongest cleavages generally occurred in the border regions between the predicted domains, as is particularly evident for yNup133 and yNup120 (Figure 3). Hence, in every case, the regions that we predicted to form compact folded structures were proteolytically resistant, and the predicted linkers between these domains were proteolytically sensitive. This correlation provides support for all seven of our structural models. In addition, circular dichroism and Fourier transform infrared spectra reported for Nup85 are in agreement with our predictions, indicating a composition characteristic of α-solenoids (approximately 50% α-helical, 23% loops, 5% turns, and 10% β-sheet) (Hirano et al. 1990; Denning et al. 2003). We expect our findings will spur efforts to determine the detailed atomic structures of nups; the rapid proteolytic domain mapping and molecular modeling techniques we have utilized here should aid these efforts.
Table 2
Table 2
Proteolytically Sensitive Sites of yNup84 Subcomplex Proteins
Having established the domain folds for the yNup84 subcomplex, we also assigned domain folds in its vertebrate (i.e., human) and plant (i.e., Arabidopsis) homologs. All seven nups from both human and Arabidopsis yielded identical domain fold assignments to their yeast counterparts (Table S7), despite low primary sequence conservation among them (Suntharalingam and Wente 2003). These findings indicate that the overall architecture of the yNup84/vNup107–160 subcomplex has been preserved throughout the eukaryotes. Hence, the yNup84/vNup107–160 subcomplex, which contributes nearly one-quarter of the mass of the NPC, is composed in the main of repetitive β-propellers and α-solenoids; taken together with other repetitive domain nups (such as the FG repeat nups), this suggests that a significant percentage of the NPC's bulk is composed of protein repeats (Rout and Aitchison 2001; Suntharalingam and Wente 2003).
To gain insight into the function and origin of the yNup84/vNup107–160 subcomplex, we asked whether there are other known subcomplexes that share similar compositions and fold arrangements. A search of the entire SwissProt/TrEMBL database for entries that contain an amino-terminal β-propeller followed by an α-solenoid revealed that this specific architectural combination is absent from both bacteria and archaebacteria, and is found only in eukaryotic proteins, whose role (where known) is as components either of coated vesicles or of the yNup84/vNup107–160 subcomplex. Thus, the clathrin heavy chain, a major component of clathrin-coated vesicles, appears remarkably similar in domain architecture (ter Haar et al. 1998; Kirchhausen 2000b) to both yNup120/vNup160 and yNup133/vNup133. All three proteins are composed of an amino-terminal β-propeller followed by an extended α-solenoid. Proteolysis of assembled clathrin cages leads to the release of an amino-terminal fragment of 52–59 kDa (Kirchhausen and Harrison 1984). This result is similar to our domain mapping results, where the proteolysis of yNup120 and yNup133 resulted in amino-terminal fragments of 45 kDa and 60 kDa, respectively. Strikingly, one component of the yNup84/vNup107–160 subcomplex, ySec13/vSec13R, is also a known vesicle-coating protein. Similarly, ySeh1/vSec13L, a close homolog of ySec13/vSec13R, is also associated with both the yNup84/vNup107–160 subcomplex and the vesicle-coating proteins (Siniossoglou et al. 1996; Kirchhausen 2000b; Cronshaw et al. 2002; Gavin et al. 2002; Harel et al. 2003). Together, these results point to an intimate connection between vesicle-coating complexes and the yNup84/vNup107–160 subcomplex.
In clathrin-coated vesicles, clathrin is attached via its amino-terminal domain to an adaptin complex. There are four types of adaptin complexes, all made of two large subunits that wrap around two small subunits. The bulk of each large subunit is made of an α-solenoid trunk (Figure 4) (Collins et al. 2002; Evans and Owen 2002). Similarly, the bulk of yNup84/vNup107, yNup85/vNup75, and yNup145C/vNup96 are also composed of α-solenoid trunks. Hence, the yNup84/vNup107–160 subcomplex resembles the clathrin/adaptin complex, in that the clathrin-like yNup120/vNup160 and yNup133/vNup133 are attached to the adaptin-like proteins yNup84/vNup107, yNup85/vNup75, and yNup145C/vNup96. This resemblance is further strengthened by our observation that the preferred templates for modeling the α-solenoid domains in the yNup84/vNup107–160 subcomplex were derived from proteins in vesicle coating complexes (Figure S1; Tables S1–S7).
Figure 4
Figure 4
The Nup84 Complex and Coated Vesicles Share a Common Architecture
Our analyses showed that the yNup84/vNup107–160 subcomplex and all three major classes of vesicle coating complexes can be linked together through their common architecture. As summarized in Figure 4, these similarities include both previously reported relationships (e.g., between the clathrin/adaptin complexes and the COPI complexes) (Schledzewski et al. 1999), and previously unsuspected relationships (e.g., between the COPII component Sec31 [Salama et al. 1997; Shugrue et al. 1999; Belden and Barlowe 2001; Boehm and Bonifacino 2001; Lederkremer et al. 2001] and clathrin).
The common architecture of the yNup84/vNup107–160 subcomplex and all three major classes of vesicle-coating complexes suggests that all of these complexes have common function in curving membranes. There is, in fact, circumstantial evidence for a role of the yNup84/vNup107–160 subcomplex in the establishment and maintenance of pore membrane curvature. Members of this complex, when disrupted in yeast, cause the uniformly distributed NPCs to cluster into patches in the plane of the NE (Siniossoglou et al. 1996; Siniossoglou et al. 2000; Ryan and Wente 2002; Teixeira et al. 2002), suggesting that impairment of yNup84 subcomplex function results in a suboptimal interaction of the NPC with its surrounding nuclear membranes.
As shown here, protein structure modeling is particularly useful in uncovering potential evolutionary and functional relationships that are refractory to classical approaches based on comparison of protein sequences alone. Our results show that clathrin/adaptin complexes, COPI complexes, COPII complexes, and the yNup84/vNup107–160 subcomplex all share a common molecular architecture. This commonality could have arisen by either convergent or divergent evolutionary pathways.
In a convergent pathway, β-propeller and α-solenoid folds could have been independently utilized by both NPCs and vesicle-coating complexes at different stages of eukaryotic evolution. This possibility is supported by the high abundance of both fold types in eukaryotic genomes (which could potentially make their fusion in proteins or complexes relatively frequent) (Yanai et al. 2002) and the low sequence similarities between proteins of the NPC and vesicle coating complexes (which may suggest that they are not related).
In a divergent pathway, NPCs and vesicle-coating complexes share these folds because both complex types could have originated from a common ancestor. In this scenario, a single “protocoatomer” would have been the progenitor for numerous vesicle coating complexes, as well as the yNup84/vNup107–160 subcomplex. Several lines of evidence support this latter hypothesis. First, the most confident models of the yNup84/vNup107–160 subcomplex proteins are based on structures of coated vesicle proteins (Figure S1; Tables S1–S7). Second, the particular arrangement of an amino-terminal β-propeller followed by an α-solenoid appears to be unique to components of either vesicle coating complexes or of the yNup84/vNup107–160 subcomplex (Protocol S1). Third, the overall composition of both complex types is similar, being mainly composed of proteins containing comparable distributions of β-propellers and α-solenoids (Figure 4). Fourth, both vesicle coating complexes and NPCs apparently share a common function: the bending and stabilizing of curved membranes. Fifth, the yNup84/vNup107–160 subcomplex actually contains bona fide vesicle coat components, Sec13 and Seh1. In light of these considerations, we favor the “protocoatomer” hypothesis, in which the NPCs and vesicle-coating complexes arose by a process of divergent evolution.
The lack of detectable sequence similarity between the proteins in the yNup84/vNup107–160 subcomplex and the coated vesicles is not surprising. Sequence comparisons of α-solenoid- and β-propeller-containing proteins suggest that these folds arose just before or around the time of the origin of eukaryotes, then rapidly duplicated and diversified (Cingolani et al. 1999; Smith et al. 1999; Andrade et al. 2001b). Both folds consist of repetitive structures, so the functional constraints on an individual repeat are weak, compared with the whole fold domain. It has been proposed that the robustness of these folds with respect to changes in their sequences permits their component repeats to individually lose their sequence similarity, eventually allowing the proteins they comprise to drift into new functions (Malik et al. 1997; Smith et al. 1999; Andrade et al. 2001a; Andrade et al. 2001b). Moreover, the lack of detectable sequence similarity for members of the same fold family is not necessarily an indicator of convergent evolution; obvious sequence similarities are often lost during long periods of evolution (e.g., FtsZ and tubulin or MreB and actin [Amos et al. 2004]). The divergent pathway is also consistent with the conservation among members of the syntaxin family (key components of the vesicular transport machinery), which points to a similar early origin and rapid diversification of the eukaryotic endomembrane system (Dacks and Doolittle 2002; Dacks and Field 2004). Based on these observations, we propose a single evolutionary origin for the structures maintaining both the endomembrane systems and the nucleus (Figure 5) over models suggesting separate or even endosymbiotic origins for these structures.
Figure 5
Figure 5
Proposed Model for the Evolution of Coated Vesicles and Nuclear Pore Complexes
The current protocoatomer hypothesis posits that a simple coating module containing minimal copies of the two conserved folds evolved in protoeukaryotes as a mechanism to bend membranes into sharply curved sheets and invaginated tubules (Figure 5). The ability to so manipulate cell membranes represented a major evolutionary innovation that allowed, among other possibilities, the elaboration of internal membranes, phagotrophy, and endosymbiosis (Maynard Smith and Szathmâary 1997); the importance of this ability is underscored by the presence of numerous types of membrane-curving devices in modern eukaryotes. As with clathrin, the flexibility of the α-solenoid in this simple module enabled the formation of curved membranes of various sizes. In addition, the α-solenoid repeat structure, together with the repeats in the β-propeller fold, provided the coating module with a large binding area. These features allowed the membrane-curving module to polymerize and form a coat, as well as to interact with other membrane-associated proteins. The endomembranes and their membrane-coating modules subsequently evolved to become more elaborate and specialized, with the partitioning of different functions into separate, interconnected compartments such as the ER, the Golgi, and the nucleus (Figure 5), each with their own specialized set of coating modules.
In conclusion, we suggest that the progenitor of the NPC arose from a membrane-coating module that wrapped extensions of an early ER around the cell's chromatin. In this primitive NE, the coating modules would have originally formed the sharply curved membrane, creating large and freely permeable pores (Figure 5). These pores then closed to form the selectively permeable NPCs of modern eukaryotes (Rout et al. 2003). In doing so, they retained at their core a coating module as a relic of their evolutionary origins. This module, the yNup84/vNup107–160 subcomplex, may still serve to curve and stabilize the nuclear pore membrane in modern eukaryotes; as such, it would function as a key scaffold to form the NPCs, the portals of the nucleus. Our findings could thus provide an explanation for the origin of the nuclear pore complex (which until now has been a mystery) and may fill a significant gap in our understanding of the evolution of eukaryotes.
Only two domains in the seven nups have their folds assigned by sequence comparison to proteins of known structure (Saxena et al. 1996; Siniossoglou et al. 1996). Therefore, to assign folds for as many target domains comprising the yNup84/vNup107–160 subcomplex as possible, we applied a structure-based approach consisting of iterative detection of potential template structures, their alignment to the target sequence, model building, and model assessment (Marti-Renom et al. 2000). Secondary structure was predicted from sequence by the PredictProtein (Rost 1996) and PSI-Pred (McGuffin et al. 2000) servers.
Detection of potential template structures
For each of the seven yeast nups and representative homologs, potentially related known structures were detected by the mGenThreader (McGuffin and Jones 2003) and FUGUE (Shi et al. 2001) web servers (Tables S1–S7). Several other servers gave similar results (unpublished data). To find out whether or not mGenThreader frequently identifies the β-propeller and α-solenoid folds as false positives, we randomly selected 20 sequences of known structure from each one of the structural classes and submitted them to mGenThreader. Using the same parameters as in our analysis of the nups, only two of these 140 sequences were incorrectly predicted to contain β-propeller or α-solenoid folds (unpublished data). Thus, we estimate the false positives rate for the nup fold assignments based on mGenThreader alone to be approximately 1%–2%.
Alignment of the matched target-template pairs
The matches obtained in the previous step provided an operational definition of a domain. They were either accepted or refined by manual and automated alignment. Manual realignment relied on sequence conservation and secondary structure predictions by PROF (Rost 1996) and PSI-PRED (McGuffin et al. 2000). The automatic realignments were obtained by SALIGN (Marti-Renom et al. 2004) and T-Coffee (Notredame et al. 2000). In the last iteration, the alignments and the models were refined by MOULDER, a genetic algorithm method for iterative alignment, model building, and model assessment (John and Sali 2003).
Model building
For each alignment, an all-atom model was built by comparative modeling based on satisfaction of spatial restraints as implemented in MODELLER (Sali and Blundell 1993).
Model assessment.
The fold assignment, alignment, and model building were repeated by varying the domain boundaries, target sequences for modeling, template structures, and their alignments. The aim was to improve model assessment by statistical potentials of ProsaII (Sippl 1993) and DFIRE (Zhou and Zhou 2002), and by a composite model evaluation criterion (Melo et al. 2002; John and Sali 2003). The only importance of explicit model building in this analysis was to provide another semi-independent way to validate the fold assignments: If a model was assessed to have the correct fold, the initial fold assignment must have been correct. Beyond that, the models were not used.
Domain combination search.
To search for proteins that resemble the domain architecture of clathrin, we queried MODBASE (Pieper et al. 2004), our relational database of annotated comparative protein structure models, and Superfamily (Gough et al. 2001), a database of HMM-based structural assignments. Both databases assign folds to all available protein sequences that match at least one known protein structure. We first searched for any protein sequences that were matched to both β-propeller and α-solenoid structures. We used the broadest definitions of the β-propeller folds (b.66, b.67, b.68, b.69, b.70, for 4-, 5-, 6-, 7- and 8-bladed β-propellers, respectively) and α-solenoid folds (a.118) from the SCOP database (v1.65) (Lo Conte et al. 2002). In MODBASE, we found 95 proteins predicted to contain both β-propeller and α-solenoid domains (Protocol S1). Of these 95 proteins, 37 passed the following filters, ensuring clathrin-like characteristics: they had to be 800–2,000 residues long, the amino-terminal β-propeller domain had to be followed by a carboxyl-terminal α-solenoid domain, the β-propeller and α-solenoid domains each had to span at least 35% of the total length, and no other domain could be more than 25% of the total length. All of the 37 proteins were from eukaryotes. Their functions were assigned either as clathrin or unknown in the Swiss-Prot/TrEMBL database (O'Donovan et al. 2002). Similar results were obtained by querying the Superfamily database (Gough et al. 2001).
Proteolytic domain laddering.
Magnetic beads (2.8 μm Dynabeads M-270 Epoxy [#143.02; Dynal, Oslo, Norway]) were conjugated to rabbit IgG (#55944; ICN Biochemicals, Costa Mesa, California, United States) according to the manufacturer's instructions. Yeast cells carrying PrA-tagged versions of nups were grown and harvested as described previously (Rout et al. 2000). Cell pellets were frozen in liquid nitrogen and homogenized to a fine powder in a motorized grinder (#RM100; Retsch, Haan, Germany) continuously cooled with liquid nitrogen. The cell powder was thawed on ice and ten volumes of extraction buffer (20 mM HEPES [pH 7.4], 1.0% Triton X-100, 0.5% sodium deoxycholate, 0.3% sodium N-lauroyl-sarcosine, 0.1 mM MgCl2, 1 mM DTT, 1:500 protease inhibitor cocktail [#P-8340; Sigma, St. Louis, Missouri, United States]) were added to cells and homogenized at 4 °C with a Polytron (Kinematica, Littau-Luzerne, Switzerland). The cell lysate was clarified by centrifugation (2,000 g for 5 min at 4 °C). The magnetic beads were added to the extract to a ratio of about 8 × 109 beads per g of cells. After incubation for 1 h at 4 °C, the beads were magnetically recovered. The beads were washed, resuspended in 50 μl of reaction buffer (according to the manufacturer's specifications), and Asp-N (#1420488; Roche, Basel, Switzerland) or Lys-C (#1420429; Roche) proteinase was added to give a weight ratio of 1:200 of proteinase to the tagged nup. After incubation at different time points at 37 °C, bead aliquots were removed and washed, and tagged fragments were eluted with 0.5 M NH4OH containing 0.5 mM EDTA. The eluant was vacuum-dried, resuspended in SDS-PAGE sample buffer, and separated on a 4%–12% bis-Tris gel (Invitrogen, Carlsbad, California, United States). Proteins were then either transferred electrophoretically to nitrocellulose or PVDF and probed with HRP-rabbit IgG (#011–0303-003; Jackson ImmunoResearch, West Grove, Pennsylvania, United States), or analyzed by amino-terminal Edman sequencing (Fernandez et al. 1994).
Figure S1: Model Score Versus Length
The graphs plot the assessment score of the model (Melo Z-score) (Melo et al. 2002) versus the model size, for the "non-MOULDER" models in Tables S2–S6. The red circles indicate the entries in Table 1 in the main text of the paper. Because the Z-score depends on the number of residues in the model, the smallest model with the highest Z-score was considered most significant.
(87 KB DOC).
Figure S2: Agreement between Predicted and Modeled Secondary Structure
The secondary structure predicted from sequence by PROF (Rost and Liu 2003) and PSI-Pred (McGuffin et al. 2003) is compared to the secondary structure observed in the three-dimensional models presented in Table S1 (“…” represents regions that are not modeled). The numbers above the predicted secondary structures correspond to the confidence score returned by the servers. Current secondary structure prediction methods based on multiple alignments correctly predict the secondary structure state for 70%–80% of residues (in a three-state prediction) (Eyrich et al. 2001). Since the random prediction would predict only about 30% of the residues correctly, the fact that our predictions match the assignments at 58%–87% level is highly suggestive, supporting our fold assignments. A representative example, Nup85, is shown here. For the visualization of all the Nups, see the additional information web page (http://salilab.org/damien/NPC/).
(47 KB DOC).
Protocol S1: List of Proteins Modeled as β-Propeller and α-Solenoid Domains in ModBase
(42 KB DOC).
Table S1: Modeling Results for Yeast Nup84 Complex Proteins I (yNup133)
(491 KB DOC).
Table S2: Modeling Results for Yeast Nup84 Complex Proteins II (yNup133)
(101 KB DOC).
Table S3: Modeling Results for Yeast Nup84 Complex Proteins III (yNup133)
(115 KB DOC).
Table S4: Modeling Results for Yeast Nup84 Complex Proteins IV (yNup133)
(132 KB DOC).
Table S5: Modeling Results for Yeast Nup84 Complex Proteins V (yNup133)
(124 KB DOC).
Table S6: Modeling Results for Yeast Nup84 Complex Proteins (yNup133)
(93 KB DOC).
Table S7: Modeling Results for Human and Plant Nup84 Complex Proteins (yNup133)
(144 KB DOC).
Accession Numbers
Uniprot (Apweiler et al. 2004) accession numbers (http://www.pir.uniprot.org) for proteins discussed in this paper are as follows. Yeast: ySeh1 (P53011), ySec13 (Q04491), yNup84 (P52891), yNup85 (P46673), yNup120 (P35729), yNup133 (P36161), and yNup145C (P49687). Human: vSec13 l (Q96EE3), vSec13R (P55735), vNup107 (P57740), vNup75 (Q9BW27), vNup160 (Q12769), vNup133 (Q8WUM0), and vNup96 (P52948).
Acknowledgments
We thank the following colleagues for their helpful contributions and discussions: Joe Fernandez and the Proteomics Resource Center of the Rockefeller University, Marc Marti-Renom, Joseph Mancias, Martine Cadene, Mark Field, Günter Blobel, John Aitchison, Kelli Mullin, John Kilmartin, Margaret Robinson, Chris Akey, Mallur Madhusudhan, Miklós Müller, Miguel Andrade, Fred Davis, and Robert Fletterick. We also acknowledge the support of the Rockefeller University, NIH GM062427, NIH/NCRR RR00862, NIH/NCI R33 CA89810, SUN, IBM, Intel, The Rita Allen and Sinsheimer Foundations, the Irma T. Hirschl Trust, and The Sandler Family Supporting Foundation.
Abbreviations
COPcoat protein
ERendoplasmic reticulum
NEnuclear envelope
NPCnuclear pore complex
nupnucleoporin
WDtryptophan/aspartic acid

Footnotes
Conflicts of interest. The authors have declared that no conflicts of interest exist.
Academic Editor: Greg Petsko, Brandeis University
Citation: Devos D, Dokudovskaya S, Alber F, Williams R, Chait BT, et al. (2004) Components of coated vesicles and nuclear pore complexes share a common molecular architecture. PLoS Biol 2(12): e380.
  • Amos LA, van den Ent F, Lowe J. Structural/functional homology between the bacterial and eukaryotic cytoskeletons. Curr Opin Cell Biol. 2004;16:24–31. [PubMed]
  • Andrade MA, Perez-Iratxeta C, Ponting CP. Protein repeats: Structures, functions, and evolution. J Struct Biol. 2001a;134:117–131. [PubMed]
  • Andrade MA, Petosa C, O'Donoghue SI, Muller CW, Bork P. Comparison of ARM and HEAT protein repeats. J Mol Biol. 2001b;309:1–18. [PubMed]
  • Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, et al. UniProt: The Universal Protein knowledgebase. Nucleic Acids Res. 2004. pp. D115–D119. [PMC free article] [PubMed]
  • Bednenko J, Cingolani G, Gerace L. Nucleocytoplasmic transport: Navigating the channel. Traffic. 2003;4:127–135. [PubMed]
  • Belden WJ, Barlowe C. Purification of functional Sec13p-Sec31p complex, a subunit of COPII coat. Methods Enzymol. 2001;329:438–443. [PubMed]
  • Blobel G. Intracellular protein topogenesis. Proc Natl Acad Sci U S A. 1980;77:1496–1500. [PubMed]
  • Boehm M, Bonifacino JS. Adaptins: The final recount. Mol Biol Cell. 2001;12:2907–2920. [PMC free article] [PubMed]
  • Boehmer T, Enninga J, Dales S, Blobel G, Zhong H. Depletion of a single nucleoporin, Nup107, prevents the assembly of a subset of nucleoporins into the nuclear pore complex. Proc Natl Acad Sci U S A. 2003;100:981–985. [PubMed]
  • Bonifacino JS, Lippincott-Schwartz J. Coat proteins: Shaping membrane transport. Nat Rev Mol Cell Biol. 2003;4:409–414. [PubMed]
  • Cingolani G, Petosa C, Weis K, Muller CW. Structure of importin-beta bound to the IBB domain of importin-alpha. Nature. 1999;399:221–229. [PubMed]
  • Collins BM, McCoy AJ, Kent HM, Evans PR, Owen DJ. Molecular architecture and functional model of the endocytic AP2 complex. Cell. 2002;109:523–535. [PubMed]
  • Cronshaw JM, Krutchinsky AN, Zhang W, Chait BT, Matunis MJ. Proteomic analysis of the mammalian nuclear pore complex. J Cell Biol. 2002;158:915–927. [PMC free article] [PubMed]
  • Dacks JB, Doolittle WF. Novel syntaxin gene sequences from Giardia, Trypanosoma and algae: Implications for the ancient evolution of the eukaryotic endomembrane system. J Cell Sci. 2002;115:1635–1642. [PubMed]
  • Dacks JB, Field MC. Eukaryotic cell evolution from a comparative genomic perspective: The endomembrane system. In: Hirt R, Horner D, editors. Organelles, genomes and eukaryote phylogeny: An evolutionary synthesis in the age of genomics. Boca Raton: CRC Press; 2004.
  • Denning DP, Patel SS, Uversky V, Fink AL, Rexach M. Disorder in the nuclear pore complex: the FG repeat regions of nucleoporins are natively unfolded. Proc Natl Acad Sci U S A. 2003;100:2450–2455. [PubMed]
  • Evans PR, Owen DJ. Endocytosis and vesicle trafficking. Curr Opin Struct Biol. 2002;12:814–821. [PubMed]
  • Eyrich VA, Marti-Renom MA, Przybylski D, Madhusudhan MS, Fiser A, et al. EVA: Continuous automatic evaluation of protein structure prediction servers. Bioinformatics. 2001;17:1242–1243. [PubMed]
  • Fahrenkrog B, Koser J, Aebi U. The nuclear pore complex: A jack of all trades? Trends Biochem Sci. 2004;29:175–182. [PubMed]
  • Fernandez J, Andrews L, Mische SM. An improved procedure for enzymatic digestion of polyvinylidene difluoride-bound proteins for internal sequence analysis. Anal Biochem. 1994;218:112–117. [PubMed]
  • Fontoura BM, Blobel G, Matunis MJ. A conserved biogenesis pathway for nucleoporins: Proteolytic processing of a 186-kilodalton precursor generates Nup98 and the novel nucleoporin, Nup96. J Cell Biol. 1999;144:1097–1112. [PMC free article] [PubMed]
  • Gavin AC, Bosche M, Krause R, Grandi P, Marzioch M, et al. Functional organization of the yeast proteome by systematic analysis of protein complexes. Nature. 2002;415:141–147. [PubMed]
  • Gough J, Karplus K, Hughey R, Chothia C. Assignment of homology to genome sequences using a library of hidden Markov models that represent all proteins of known structure. J Mol Biol. 2001;313:903–919. [PubMed]
  • Harel A, Orjalo AV, Vincent T, Lachish-Zalait A, Vasu S, et al. Removal of a single pore subcomplex results in vertebrate nuclei devoid of nuclear pores. Mol Cell. 2003;11:853–864. [PubMed]
  • Hirano T, Kinoshita N, Morikawa K, Yanagida M. Snap helix with knob and hole: Essential repeats in S. pombe nuclear protein nuc2+ Cell. 1990;60:319–328. [PubMed]
  • John B, Sali A. Comparative protein structure modeling by iterative alignment, model building and model assessment. Nucleic Acids Res. 2003;31:3982–3992. [PMC free article] [PubMed]
  • Kirchhausen T. Clathrin. Annu Rev Biochem. 2000a;69:699–727. [PubMed]
  • Kirchhausen T. Three ways to make a vesicle. Nat Rev Mol Cell Biol. 2000b;1:187–198. [PubMed]
  • Kirchhausen T, Harrison SC. Structural domains of clathrin heavy chains. J Cell Biol. 1984;99:1725–1734. [PMC free article] [PubMed]
  • Koh IY, Eyrich VA, Marti-Renom MA, Przybylski D, Madhusudhan MS, et al. EVA: Evaluation of protein structure prediction servers. Nucleic Acids Res. 2003;31:3311–3315. [PMC free article] [PubMed]
  • Lederkremer GZ, Cheng Y, Petre BM, Vogan E, Springer S, et al. Structure of the Sec23p/24p and Sec13p/31p complexes of COPII. Proc Natl Acad Sci U S A. 2001;98:10704–10709. [PubMed]
  • Lippincott-Schwartz J, Liu W. Membrane trafficking: Coat control by curvature. Nature. 2003;426:507–508. [PubMed]
  • Lo Conte L, Brenner SE, Hubbard TJ, Chothia C, Murzin AG. SCOP database in 2002: Refinements accommodate structural genomics. Nucleic Acids Res. 2002;30:264–267. [PMC free article] [PubMed]
  • Loiodice I, Alves A, Rabut G, Van Overbeek M, Ellenberg J, et al. The entire nup107–160 complex, including three new members, is targeted as one entity to kinetochores in mitosis. Mol Biol Cell. 2004;15:3333–3344. [PMC free article] [PubMed]
  • Lutzmann M, Kunze R, Buerer A, Aebi U, Hurt E. Modular self-assembly of a Y-shaped multiprotein complex from seven nucleoporins. EMBO J. 2002;21:387–397. [PubMed]
  • Malik HS, Eickbush TH, Goldfarb DS. Evolutionary specialization of the nuclear targeting apparatus. Proc Natl Acad Sci U S A. 1997;94:13738–13742. [PubMed]
  • Marti-Renom MA, Madhusudhan MS, Sali A. Alignment of protein sequences by their profiles. Protein Sci. 2004;13:1071–1087. [PubMed]
  • Marti-Renom MA, Stuart AC, Fiser A, Sanchez R, Melo F, et al. Comparative protein structure modeling of genes and genomes. Annu Rev Biophys Biomol Struct. 2000;29:291–325. [PubMed]
  • Maynard Smith J, Szathmâary E. Oxford: Oxford University Press; 1997. The major transitions in evolution; 360 pp.
  • McGuffin LJ, Jones DT. Improvement of the GenTHREADER method for genomic fold recognition. Bioinformatics. 2003;19:874–881. [PubMed]
  • McGuffin LJ, Bryson K, Jones DT. The PSIPRED protein structure prediction server. Bioinformatics. 2000;16:404–405. [PubMed]
  • Melo F, Sanchez R, Sali A. Statistical potentials for fold assessment. Protein Sci. 2002;11:430–448. [PubMed]
  • Neer EJ, Schmidt CJ, Nambudripad R, Smith TF. The ancient regulatory-protein family of WD-repeat proteins. Nature. 1994;371:297–300. [PubMed]
  • Notredame C, Higgins DG, Heringa J. T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000;302:205–217. [PubMed]
  • O'Donovan C, Martin MJ, Gattiker A, Gasteiger E, Bairoch A, et al. High-quality protein knowledge resource: SWISS-PROT and TrEMBL. Brief Bioinform. 2002;3:275–284. [PubMed]
  • Orengo CA, Michie AD, Jones S, Jones DT, Swindells MB, et al. CATH--A hierarchic classification of protein domain structures. Structure. 1997;5:1093–1108. [PubMed]
  • Pieper U, Eswar N, Braberg H, Madhusudhan MS, Davis FP, et al. MODBASE, a database of annotated comparative protein structure models, and associated resources. Nucleic Acids Res. 2004. pp. D217–D222. [PMC free article] [PubMed]
  • Rost B. PHD: Predicting one-dimensional protein structure by profile-based neural networks. Methods Enzymol. 1996;266:525–539. [PubMed]
  • Rost B, Liu J. The PredictProtein server. Nucleic Acids Res. 2003;31:3300–3304. [PMC free article] [PubMed]
  • Rout MP, Aitchison JD. The nuclear pore complex as a transport machine. J Biol Chem. 2001;276:16593–16596. [PubMed]
  • Rout MP, Aitchison JD, Suprapto A, Hjertaas K, Zhao Y, et al. The yeast nuclear pore complex: Composition, architecture, and transport mechanism. J Cell Biol. 2000;148:635–651. [PMC free article] [PubMed]
  • Rout MP, Aitchison JD, Magnasco MO, Chait BT. Virtual gating and nuclear transport: The hole picture. Trends Cell Biol. 2003;13:622–628. [PubMed]
  • Ryan KJ, Wente SR. Isolation and characterization of new Saccharomyces cerevisiae mutants perturbed in nuclear pore complex assembly. BMC Genet. 2002;3:17. [PMC free article] [PubMed]
  • Salama NR, Chuang JS, Schekman RW. Sec31 encodes an essential component of the COPII coat required for transport vesicle budding from the endoplasmic reticulum. Mol Biol Cell. 1997;8:205–217. [PMC free article] [PubMed]
  • Sali A, Blundell TL. Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993;234:779–815. [PubMed]
  • Saxena K, Gaitatzes C, Walsh MT, Eck M, Neer EJ, et al. Analysis of the physical properties and molecular modeling of Sec13: A WD repeat protein involved in vesicular traffic. Biochemistry. 1996;35:15215–15221. [PubMed]
  • Schledzewski K, Brinkmann H, Mendel RR. Phylogenetic analysis of components of the eukaryotic vesicle transport system reveals a common origin of adaptor protein complexes 1, 2, and 3 and the F subcomplex of the coatomer COPI. J Mol Evol. 1999;48:770–778. [PubMed]
  • Shi J, Blundell TL, Mizuguchi K. FUGUE: Sequence-structure homology recognition using environment-specific substitution tables and structure-dependent gap penalties. J Mol Biol. 2001;310:243–257. [PubMed]
  • Shugrue CA, Kolen ER, Peters H, Czernik A, Kaiser C, et al. Identification of the putative mammalian orthologue of Sec31P, a component of the COPII coat. J Cell Sci. 1999;112:4547–4556. [PubMed]
  • Siniossoglou S, Wimmer C, Rieger M, Doye V, Tekotte H, et al. A novel complex of nucleoporins, which includes Sec13p and a Sec13p homolog, is essential for normal nuclear pores. Cell. 1996;84:265–275. [PubMed]
  • Siniossoglou S, Lutzmann M, Santos-Rosa H, Leonard K, Mueller S, et al. Structure and assembly of the Nup84p complex. J Cell Biol. 2000;149:41–54. [PMC free article] [PubMed]
  • Sippl MJ. Recognition of errors in three-dimensional structures of proteins. Proteins. 1993;17:355–362. [PubMed]
  • Smith TF, Gaitatzes C, Saxena K, Neer EJ. The WD repeat: A common architecture for diverse functions. Trends Biochem Sci. 1999;24:181–185. [PubMed]
  • Suntharalingam M, Wente SR. Peering through the pore: Nuclear pore complex structure, assembly, and function. Dev Cell. 2003;4:775–789. [PubMed]
  • Teixeira MT, Dujon B, Fabre E. Genome-wide nuclear morphology screen identifies novel genes involved in nuclear architecture and gene-silencing in Saccharomyces cerevisiae. J Mol Biol. 2002;321:551–561. [PubMed]
  • ter Haar E, Musacchio A, Harrison SC, Kirchhausen T. Atomic structure of clathrin: A beta propeller terminal domain joins an alpha zigzag linker. Cell. 1998;95:563–573. [PubMed]
  • Walther TC, Alves A, Pickersgill H, Loiodice I, Hetzer M, et al. The conserved Nup107–160 complex is critical for nuclear pore complex assembly. Cell. 2003;113:195–206. [PubMed]
  • Yanai I, Wolf YI, Koonin EV. Evolution of gene fusions: Horizontal transfer versus independent events. Genome Biol. 2002;3:research0024.1–0024.13. [PMC free article] [PubMed]
  • Yu L, Gaitatzes C, Neer E, Smith TF. Thirty-plus functional families from a single motif. Protein Sci. 2000;9:2470–2476. [PubMed]
  • Zhou H, Zhou Y. Distance-scaled, finite ideal-gas reference state improves structure-derived potentials of mean force for structure selection and stability prediction. Protein Sci. 2002;11:2714–2726. [PubMed]
Articles from PLoS Biology are provided here courtesy of
Public Library of Science