INTRODUCTION

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

It is often technically challenging to analyze extraordinarily large proteins because of their size. In particular, it is extremely difficult to clone the full-length cDNAs encoding such huge proteins and equally difficult to purify intact, large mRNAs and their proteins. For these reasons, the analysis of huge proteins proceeds relatively slowly. However, because such large proteins often are known to play central roles in nervous system development, brain function, and neurological disorders (Lieber et al. 1992; D'Arcangelo et al. 1995), we believe that extraordinarily large proteins are important players in biology, warranting in-depth study. Moreover, their exceptionally large molecular size most likely contributes to a unique biological function. We are therefore interested in the comprehensive characterization of these large proteins in the brain; and over the course of the past six years, we have conducted a cDNA-sequencing project focused on human brain cDNA clones that encode large proteins.

In order to identify unknown human genes overlooked using conventional methods, we have used the following two approaches. The first approach is based on random sampling of long cDNA clones that produce proteins larger than 50 kD in an in vitro transcription/translation system (Nagase et al. 2000). Thus far, we have determined >1900 cDNA sequences (named KIAA) that have been subsequently deposited into public databases. The data obtained to date indicate that: (1) the predicted functions of large proteins are mainly classified into cell communication/signaling (e.g., receptor, adhesion molecule, and GTPase-activating protein), cell structure/motility (e.g., cytoskeletal proteins, membrane skeletal proteins, and motor proteins), and nucleic acid management (e.g., DNA-binding protein; Nagase et al. 1997); (2) although difficult to identify using conventional methods, homologous genes are efficiently isolated by this approach; and (3) several large unknown proteins with unknown functions are predicted from our cDNA sequences. The second approach is to systematically identify genes encoding proteins with a particular domain of interest. The epidermal growth factor (EGF)-like domain is thought to play an important role in extracellular events, including cell adhesion and receptor-ligand interactions. Because the amino acid sequences of EGF-like domains in various known proteins have little in common, except for the conserved location of six cysteine residues, genes coding for EGF-like domain-containing proteins could not be systematically identified either by colony hybridization under relaxed conditions or by the PCR method using degenerate primers. We have therefore developed a computer-assisted method, called motif-trap screening, for this purpose and have succeeded in identifying many large proteins with multiple EGF-like motifs and/or cadherin domains (Nakayama et al. 1998; Nagase et al. 2000; Nakajima et al. 2001).

One of the main research themes of our group is the comprehensive understanding of large proteins in the brain. Our ultimate goal is to deeply understand interactions between transmembrane proteins (Nakayama et al. 1998; Nagase et al. 2000; Nakajima et al. 2001; Mitsui et al. 2002) and proteins localized in submembrane regions (Ohara et al. 1998, 1999, 2000; Yamakawa et al. 1999). This work is prompted by the idea that interactions between transmembrane proteins and proteins localized in submembrane regions are necessary for developmental processes, synaptic plasticity, and neuronal events dependent on electrical activity (e.g., Wyszynski et al. 1997; Lisman and McIntyre 2001). One particular protein complex in neurons capturing the spotlight is the postsynaptic density (PSD), which appears as an electron-dense complex upon electron microscopic examination. PSDs are comprised of many components, including calmodulin kinase, glutamate receptors, and PSD95 (Kennedy 2000). As of yet, however, little information exists regarding the nature of the high-molecular-weight proteins present in PSDs. Large proteins have multiple domains potentially capable of binding many kinds of partners. It is conceivable, therefore, that such large proteins could function as an intricate framework of assembly protein complexes.

Therefore, as part of an ongoing study of large-molecular-weight proteins, we have characterized interactions between large proteins. Two types of approaches were used to analyze genome-wide protein-protein interactions using the yeast two-hybrid system, namely, matrix compilation and library screening (Legrain and Selig 2000).

Many technical problems are encountered, however, when using standard methods to analyze large proteins with yeasts. Even when using either one-by-one or high-throughput matrix approaches (Uetz et al. 2000; Ito et al. 2001), stable activator domain-fused, full-length proteins must be produced reliably in yeast. It is rarely possible, however, to stably produce large proteins derived from other species in yeasts. Analysis of large proteins using library screening is also technically challenging. If conventional methods were to be used, most higher eukaryotic cDNA libraries would contain mainly highly redundant cDNAs and fewer long cDNAs (Ohara et al. 1997). Because commercially available expression libraries for yeast two-hybrid systems simply represent conversions of cDNA libraries produced using conventional methods, one cannot expect to find large amounts of long cDNAs in these commercially available libraries. Even if an expression library is prepared using random primers, it is doubtful that such libraries would contain cDNAs encoding the N-terminal regions of large proteins and contain enough cDNA fragments covering full-length, long cDNAs. For these reasons, we have constructed a normalized representative library composed of long cDNAs. This novel library is comprised of KIAA cDNAs identified in our laboratory. Using the yeast two-hybrid system and this unique expression library, we have identified many protein-protein interactions using the cytoplasmic domains of large transmembrane proteins as bait.

	RESULTS

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

Construction of a Functionally Classified Library Composed of Long cDNAs

Figure 1 illustrates our strategy for efficiently examining interactions between a large number of protein-protein pairs. When combinations of unclassified protein pairs are tested (Fig. 1A), interactions are detected with relatively low frequency. However, if protein pairs are matched according to similar function or cellular location derived from experimental results or bioinformatics, the probability of detecting interactions is much greater. Thus, in this study, we focused on interactions between proteins localized in the submembrane region and the cytoplasmic domain of transmembrane proteins. Because putative transcription factors in the nuclear and extracellular domains of the transmembrane protein are not located within the cytoplasm, to avoid false-positive interactions, these domains were removed from the prey library, as described in Methods. Proteins of unknown function and location were not removed from the prey library.

View larger version (30K): [in this window] [in a new window]   Figure 1   Diagram of our strategy for constructing a functionally classified library composed of long cDNAs. (A) Illustration of different combinations according to functionally classified order. Because there are ~30,000 human proteins, our ultimate aim would be to prepare a complete protein-protein interaction map by assessing interactions between ~30,000 × ~30,000 protein combinations. One dot represents the interaction between one vertical-ordered protein (i.e., left, ordinate) and one horizontal-ordered protein (i.e., left, abscissa). The set of possible interactions to be assayed is extremely large. If proteins are, instead, first categorized according to their location and function, as shown in the right panel, interaction dots are clustered to smaller regions, and thus the set of possible interactions to be investigated is reduced considerably. (B) Schematic diagram of library construction. Black and white bars represent open reading and untranslated regions, respectively. Triangle indicates the first methionine codon. After functional classification in silico, randomly fragmented cDNAs were introduced into vectors for yeast two-hybrid screening using the Gateway system.

To date, we have documented 1987 long cDNAs (named KIAA) that encode large proteins. The average length of the protein predicted from KIAA cDNAs is 872 amino acids. Of these KIAA proteins, 175 are >1500 amino acids in length. Because it is difficult to stably express large ORFs in yeast, we chose a strategy that enables the expression of short fragments of proteins containing each domain at equal ratios. Of the 1642 KIAA genes present during the time of library construction, 1087 KIAA genes for proteins localized in the submembrane region that can potentially interact with the cytoplasmic domain of transmembrane proteins were selected (Supplemental Table 1 online contains a list of these KIAA genes). To produce random fragments, these cDNAs were pooled at almost equal ratios and sonicated. cDNA fragments encoding 333-666 amino acids were isolated and introduced into entry vectors (Fig. 1B). These lengths were chosen because we reasoned from our experience that they could be expressed in yeast stably. Because distinction of their three reading frame and orientation cannot be achieved using this method, only one-sixth of the cDNA fragments produced AD fusion proteins with the correct reading frames for the prey proteins.

The total numbers of entry clones, expression clones in Escherichia coli, and expression clones in yeast implemented in this study were 3.5 × 10⁶, 7.3 × 10⁶, and 16 × 10⁶, respectively. The average sizes of the cDNA inserts, the ORFs, and the vectors were approximately 5 kb, 2.6 kb, and 3 kb, respectively. Therefore, 167 different parts of any particular protein per ORF were produced in an in-frame orientation. In fact, as shown in Figure 2, after two-hybrid screening, many independent clones containing various protein regions with core-binding domains were isolated. This indicates that this library is composed of many kinds of expressible cDNA fragments, probably covering all ORF regions.

View larger version (39K): [in this window] [in a new window]   Figure 2   Yeast two-hybrid screening and subsequent determination of interacting domains. The domain organization for bait (left) and prey (right) KIAA proteins are shown. Blue bars indicate independent clones (one bar represents one prey). The length of the blue bar indicates the sequence of the KIAA protein contained within an isolated prey. The interacting domains located in the bait and preys are indicated by pink (bait) and green (prey) brackets. In the case of KIAA1634, we estimated that KIAA1132 can bind two different domains: the first and second PDZ domain, as shown.

Selection of Cytoplasmic Domains of Extra-Large KIAA Transmembrane Proteins as Baits

For baits, we selected 71 KIAA cDNAs on the basis of three features: (1) Expression in the brain. Only KIAA cDNAs derived from human brains were used. The initial KIAA library contained both brain-derived and human KG-1 myeloblast-derived cDNAs. (2) Protein size. KIAA cDNAs encoding peptides of >1000 amino acids in length were used. (3) Protein location and type (i.e., transmembrane proteins). The hydrophobicity of each KIAA protein was predicted using the SOSUI program, and only proteins with more than one transmembrane region (not including the N-terminal signal sequence) were used.

To reduce self-activation-associated noise measured with the two-hybrid screening, an assay for self-activation was performed before screening. As shown in Table 1, 19% of the baits assayed displayed self-activation. This percentage is consistent with previous results obtained from comprehensive research of yeast proteomes (Uetz et al. 2000). Although it is known that acidic amphipathic domains are often responsible for unwanted transcriptional activation or noise (Ruden et al. 1991), to the best of our knowledge, bait self-activation has never been assessed in silico. Thus, in the present study, we provide a comprehensive list of bait amino acid sequences displaying self-activation activity. These data supply important basic information that can be used to create computer programs for predicting potential bait self-activation sequences (see Supplemental Table 2 online). Such a program would be very useful in the field of higher eukaryotic proteome research for raising the efficiency of screening by eliminating 19% of the noise caused by bait self-activation.

Yeast Two-Hybrid Screening Using Bait and Prey Vectors With the Gateway System

The Gateway system is a useful tool to easily transfer one or more cDNA fragments into various vectors, all the while maintaining orientation and reading frame in parallel reactions contained within one tube. We introduced this system into our vector for yeast two-hybrid screening because we plan to convert our representative cDNA library and the screened prey cDNA into other types of expression vectors (for an example, see Fig. 3 and Discussion). The ProQuest Gateway recombination system was made commercially available after we constructed the library described above using the standard Gateway recombination and yeast two-hybrid screening methods. The backbone of our yeast two-hybrid screening, however, is the MATCHMAKER GAL4-based two-hybrid system supplied by Clontech. Unlike the ProQuest system, both the yeast DNA-BD and AD vectors contain multicopy plasmids. Yeast AH109 containing the bait plasmid and its mating partner, Y187, containing the prey plasmid were mated together, and the resulting diploid yeasts were selected according to nutritional requirements for Ade, His, Leu, and Trp. Candidates for a two-hybrid interaction were retested with the same nutritional selective plates containing X--Gal for the colorimetric detection of the MEL1 reporter gene product, -galactosidase. These three reporters are driven by different promoters. This triple selection method using 3-amino-1,2,3-triazole (3-AT) to suppress leaky His3 expression is thought to be one of the most effective selection systems for yeast two-hybrid screening at the present time.

View larger version (28K): [in this window] [in a new window]   Figure 3   Reassessment of interactions after exchanging cDNAs in baits and preys via the Gateway system. Image of yeast two hybrids grown on selection media. After exchanging the plasmids contained in the original bait and prey using the Gateway system, new bait plasmids (KIAA0777, KIAA1256, KIAA1296, or KIAA1496), new prey plasmid (KIAA1368), and an empty vector (pACTGW-attR) were transformed into yeast AH109. The viability of the transformants was assessed on SD/Leu/-Trp and SD/Ade/His/Leu/Trp plates. The + symbol located in the right-hand column indicates that the MEL1 gene product, -galactosidase, is actively expressed. Blue-colored growth on nutritional selection plates in the presence of X--Gal indicates active expression of -galactosidase.

Yeast Two-Hybrid Screening and Confirmation of Protein Interactions

In the present study, 2 × 10⁸ CFU of prey library per bait was used; therefore, if the mating efficiency is >1.75% (normally it is >3%), we estimate that all possible protein-protein combinations (3.5 × 10⁶ independent clones) will be tested. Because a large number of clones were screened, several independent clones containing the same binding domain were isolated in many cases (Table 1; Fig. 2). When two or more clones containing the same binding domain were detected, at least in yeast, the interaction between two proteins is real and reproducible. An example of KIAA1368 is shown in Figure 2. Using its cytoplasmic domain as bait, 14 prey colonies were isolated, and their inserts were sequenced. Every sixth prey was either KIAA1296 or KIAA1494. Every other prey was KIAA0777 or KIAA1256. Although prey clones originating from KIAA1296 are independent clones with different sizes of cDNA fragments, all preys contained an identical region consisting of KIAA1296 protein. Thus, as shown in Figure 2, one advantage of screening randomly generated fragments is that functional interacting domains can be directly mapped alongside the screening experiment that identified the interacting protein. With regard to the interacting domains of the four preys described above, namely, the binding partners of KIAA1368/semaphorin 6A (Fig. 2), all of the binding regions of these four prey KIAA proteins contained the SH3 (Src homology 3) domain. Moreover, the amino acid sequence of the binding region of KIAA1296 is similar to that of KIAA0777 and KIAA1494 (59% and 38% sequence identity, respectively).

To confirm the yeast two-hybrid screening results, the cDNAs in the baits and preys used in the initial screening were exchanged, and the resulting plasmid DNAs were transformed into haploid yeast, AH109, to assess for protein-protein interactions. Using the Gateway system, bait and prey inserts can easily be exchanged. Figure 3 shows the result of retesting protein interactions after this exchange. Transformants bearing baits containing KIAA0777, KIAA1256, or KIAA1494 and empty prey vectors resulted in no growth on the selection plates. The transformants bearing bait containing KIAA1296 and empty prey vectors, however, displayed self-activation. On the other hand, transformants bearing baits containing KIAA0777, KIAA1256, or KIAA1494 and prey containing KIAA1368/semaphorin 6A exhibited growth on selection plates containing X--Gal. These colonies were blue, indicating that the MEL1 reporter, -galactosidase, was activated. Retests involving KIAA0777, KIAA1256, and KIAA1494 clones confirmed the interaction with KIAA1368/semaphorin 6A.

The bait-prey exchange experiments confirmed the results of the initial two-hybrid screening data and demonstrated that these initial results represented true interactions, even in cases in which only one clone (e.g., KIAA0077 and KIAA1256) was obtained for DNA sequencing. Moreover, bait-prey exchange experiments that examined interactions between KIAA1206 and KIAA0076, KIAA0380, KIAA0708, or KIAA1634 (data not shown) were also confirmed.

Remarkable Interaction Pairs Following Their Annotation and Predicted Functional Domain

As shown in Table 2, we have identified 84 protein-protein interactions (more detailed accounts of these interactions are available at http://www.kazusa.or.jp/huge/ppi/). Since our initial publications (Ohara et al. 1997; Nagase et al. 1999), some KIAA proteins have been identified and described by other groups using different criteria, and rodent KIAA homologs have been analyzed extensively. The function of most KIAA proteins, however, remains unknown. Our data presented here shed light on the potential functions of some of these KIAA proteins because specific binding sites on KIAA partners were identified. Furthermore, a number of notable interaction pairs were identified in our comprehensive analysis (Fig. 2). KIAA1368/semaphorin 6A, a transmembrane semaphorin that can act as a repulsive axon guidance cue, bound KIAA proteins containing multiple SH3 domains. Multiple SH3 domains are suggested to act as an adapter that mediates signal transduction (Tong et al. 2000; Li et al. 2001). KIAA1206/Plexin-B3, a neuronal receptor for semaphorin, bound proteins containing a PDZ domain and proteins that regulate small G proteins. Its binding partner, KIAA0380, contains both PDZ and RhoGEF domains encoding the guanine nucleotide exchange factor for Rho/Rac/Cdc42-like GTPases. It is worth noting that Vikis et al. (2000) reported that human plexin-B1 specifically interacts with active Rac in a ligand-dependent manner. Hu et al. (2001) also reported that, in Drosophila, plexin B mediated axon guidance by simultaneously inhibiting Rac activity and enhancing RhoA signaling. Moreover, one of its other binding partners, KIAA0076, contains a RhoGAP domain for GTPase activator proteins, although its scores for sequence family classification using Pfam search is low (E-value, 9.3). As shown in Table 2, KIAA0708, a homolog of KIAA0076, can also bind KIAA1206/Plexin-B3. Comparing the interacting regions of KIAA0076 to that of KIAA0708, we determined that these regions bore the same amino acid sequences, which may represent a new binding domain responsible for organizing protein complexes in the submembrane region. KIAA1132, an uncharacterized adhesion molecule, also binds proteins containing PDZ domains and KIAA1092/phospholipase C, which is thought to be involved in the production of second messengers.

We found that 16 KIAA proteins interact with the PDZ domain and target the C-terminal tail, as shown in Table 3A. KIAA0343, KIAA0816, KIAA1132, and KIAA1206 were designated as class I, and KIAA0578 and KIAA1763 were categorized as class II, according to the classification scheme of Tochio et al. (1999). Although the C terminus of KIAA1206, TDL, does not exactly match the class I consensus sequence (S/T-X-I/V), the properties of Leu are very similar to that of Ile; thus, KIAA1206 was designated as possessing a variant class I PDZ domain and was grouped into class I, according to Harris and Lim (2001). Because the C terminus of KIAA1216, HSM, was never described before (Harris and Lim 2001; Vaccaro et al. 2001), PDZ domains of KIAA0316 and KIAA0807 might be assigned to a new specificity class, which we temporarily designated as class new. Interactions between the PDZ domain and the target C terminus occur in vivo in many cases (Harris and Lim 2001). Although the interactions shown in Table 2 were not reported previously, our data strongly support that true interactions were detected in our screening. In addition, as shown in Table 3B, several interactions involving the SH3 domain were also identified.

As shown in Table 3, KIAA proteins that interact with the cytoplasmic domains of transmembrane proteins often contain distinctive domains (e.g., PDZ domain, SH3 domain, spectrin domain, or Pleckstrin homology domain) known to be localized within the submembrane region. These observations are also consistent with the target at which this study aims, and support our idea that the interactions observed in our two-hybrid screening were, indeed, true interactions. Taken together, many of the interactions we identified in this study are most likely to be involved in extra-large machinery composed of multiple protein components, acting as scaffolds that organize signaling complexes.

	DISCUSSION

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

Strategies for Efficient Testing of Protein-Protein Interactions

In the present study, we sought to efficiently examine interactions between large-molecular-weight proteins. Inspection of the draft DNA sequence of the human genome indicates that there are an estimated 30,000-40,000 genes that encode proteins (Lander et al. 2001). The ultimate aim would be to test all possible interaction combinations between proteins encoded by the predicted human ORFs (i.e., 0.9-1.6 × 10⁹ combinations). The vastness of such an undertaking can be better appreciated when comparing the total number of these combinations with the size of the human genome (i.e., 3 × 10⁹ bp). To obtain a higher ratio of interaction per tested combination, we attempted to avoid testing combinations with very low probabilities of interaction based on disjunctive locations and functions of proteins. Figure 1 summarizes this strategy. When combinations of unclassified protein pairs are tested (Fig. 1A), interactions are rarely detected. However, if protein pairs are matched according to similar function or location, the probability of detecting interactions is much greater. For example, one would expect that a protein localized in the nucleus binds other proteins localized in the nucleus. Similarly, a protein localized in the cytoplasm would likely interact with other proteins in the cytoplasm. One can be confident of a small probability of false positives if reporter genes of a yeast two-hybrid screening system are activated because of an interaction between a nuclear protein and an extracellular protein. A priori classification of potential protein-protein pairs based on known function or location will decrease the number of false positives by avoiding combinations that are unlikely to colocalize in vivo. There are exceptions, however, in which some proteins may move to different parts of the cell, where they may interact with other proteins. Our guiding principal, therefore, for this study was that only protein-protein combinations with a higher probability of interaction would be tested first.

Another problem encountered in examining interactions between large-molecular-weight proteins using the yeast two-hybrid system is that it is difficult to stably express large ORFs in yeast. To address this problem, we chose a strategy that enables the expression of short fragments of proteins containing each domain at equal ratios. cDNAs were pooled at almost equal ratios and sonicated to produce random fragments. Short cDNA fragments (333-666 amino acids) that could be stably expressed in yeast were isolated and introduced into entry vectors (Fig. 1B). Although one report describes a yeast expression vector (bearing in-frame cloned ORFs) that can be selected in enrichment media (Holz et al. 2001), selection with a yield of >60% for ORF clones with correct reading frames is still difficult to achieve. Moreover, bias arising from such preselection enrichment techniques cannot be avoided. Because prey clones containing cDNAs fused out-frame with the AD produce mostly short peptides that do not have binding activity, we chose normalization without selection even if effective clones are decreased to one-sixth of all clones. Instead of excluding out-of-frame clones, special attention was paid to maintaining an adequately large number of independent clones. By preparing the clones in this manner, we were assured that a maximum number of protein pairs could be assessed using two-hybrid screening.

The Success Rate for Detecting Interactions per Tested Combination

Previously, Uetz et al. (2000) obtained the success rate of one interaction per 30,000 combinations by systematic analysis of the yeast two-hybrid system for yeast proteins. Using a mammalian two-hybrid screening method to assess interactions between mouse proteins, Suzuki et al. (2001) detected approximately one interaction for every 80,000 combinations tested. In the present study, we detected 84 interactions per 52 baits tested in a library composed of 1087 KIAA cDNAs. Although we cannot compare our interaction values directly with those of Suzuki et al. (2001), because the two screening systems implemented were different, our ratio of interactions per combination tested was relatively higher. For one combination (i.e., one bait for every one KIAA cDNA in the library), our success rate was one interaction for every 673 combinations tested. There are several reasons for this higher ratio. First, because each protein was produced from a short piece of cDNA containing a binding domain, rather than a longer cDNA containing the full ORF, a larger amount of bait protein may have been expressed stably in the yeast. This decreases the number of false negatives inherent to the use of full ORFs in two-hybrid screening using the matrix approach. Second, in many cases, the cytoplasmic domain of transmembrane proteins binds two or more proteins (Kennedy 2000; Angst et al. 2001; Harris and Lim 2001). For this reason, the cytoplasmic domain is thought to be a good target for detecting interactions efficiently. Thus, our baits, which contain such cytoplasmic domains, tend to interact with many more proteins compared with the average interaction per protein. As a result, we detected many interactions at a high rate. Third, when the screening is performed with a functionally classified library, the likelihood of detecting true interactions is maximized by first removing the region encoding the nuclear proteins and extracellular proteins that extend into the extracellular matrix. Fourth, false positives cannot be excluded when using the yeast two-hybrid system to assess protein-protein interactions. However, because two or more clones containing overlapping cDNAs were isolated, our interaction data are consistent with previously described criteria from comprehensive protein-protein analyses (Uetz et al. 2000; Ito et al. 2001). Interaction analyses using other methodologies (e.g., BIAcore analysis, GST pull down, coimmunoprecipitation, and immunohistochemistry for localization) should be performed to confirm our yeast two-hybrid data. The Gateway system, the method we used to introduce the cDNA fragment into our vectors for two-hybrid screening, is also helpful to transfer the screened candidate cDNA fragment into other expression vectors for other methodologies in order to confirm and analyze extensively our yeast two-hybrid data.

The Unique Position Our Screening Occupies in the Field of Comprehensive Protein-Protein Interaction Research

Data resulting from our comprehensive screening of KIAA proteins comprise one step in the ongoing analyses of protein-protein interactions. Ultimately, our data, when taken together with protein-interaction data from other groups, will contribute to the complete mapping of all possible protein-protein interactions. Of course, our library contains only a limited number of long cDNAs; thus, interaction information about proteins encoded by genes not presently found in our library is unknown. Nevertheless, our library is inherently useful, and, unlike with other conventional libraries, we have a clear grasp of all the cDNAs of our library. When more cDNAs are made available as a result of further cDNA analysis, we will be able to construct new libraries containing cDNAs not contained in our present library. Because our libraries will never contain overlapping sequences, duplicate testing of protein-protein pairs derived from combinations from all of these libraries (i.e., the mixed matrix approach) will be avoided.

As an alternative approach to analyze large proteins, we have an ongoing project to construct another type of library specifically intended to analyze long cDNAs. We previously compiled a specialized cDNA library for the isolation of unidentified long cDNAs that encode large proteins (Ohara et al. 1997). The work was prompted by the fact that typical cDNA libraries generated using conventional procedures based on growth rate tend to bias the selection of plasmids containing short cDNAs as opposed to long cDNAs. To avoid this growth-rate bias, long cDNAs are now being selected following separation via agarose gel electrophoresis; this procedure enables the enrichment of long cDNAs of approximately the same size. Using this series of strictly size-fractionated cDNAs as starting material, we are able to construct another library using the same procedure as reported in this study, namely, random fragmentation and introduction to two-hybrid vectors. Although this library should contain a large variety of long cDNAs, one problem encountered is redundancy at the level of expression. Because of the distinctive features of each of these libraries, we suggest that one of these two libraries should be chosen to match the user's specific needs.

Why We Focus on Protein-Protein Interactions Between Large Proteins: The Putative Role of Large Proteins in Biological Multiprotein Complexes

One good example of a large protein that participates in biological multiprotein complexes is RNA polymerase II. Because its components and subunit assembly are well characterized (Woychik and Hampsey 2002), it serves as a good model of how large proteins function in multiprotein complexes. Eukaryotic RNA polymerase II, which catalyzes DNA-dependent synthesis of mRNA, is composed of 12 subunits. Its overall configuration is "jaw-shaped." Two of its largest subunits (1733 and 1224 amino acids in Saccharomyces cerevisiae) contribute to this shape (Allison et al. 1985; Goffeau et al. 1996) in which the largest subunit (Rpb1) forms most of the lower jaw, and the second largest subunit (Rpb2) forms most of the upper jaw. These two large subunits contain both catalytic domains and binding sites to the other subunits, and play a central role in RNA polymerase function, structure, and subunit assembly.

Another example of a large protein involved in multiprotein complexes is the polycomb group in Drosophila. These proteins form large multiprotein complexes composed of 30 proteins (Saurin et al. 2001). Polycomb repressive complex 1, which was identified via biochemical purification techniques, contains at least two large proteins, posterior sex combs (PSC, 1603 amino acids; Brunk et al. 1991) and polyhomeotic-proximal chromatin protein (PH, 1589 amino acids; DeCamillis et al. 1992). Genetic experimental results show that PSC and PH are essential in maintaining a repressed transcriptional state for homeotic genes; to date, however, their roles in complex assembly have not been described in detail.

Recently, Gavin et al. (2002) described 232 distinct multiprotein complexes in S. cerevisiae using tandem-affinity purification (TAP) and mass spectrometry techniques. To compare these complexes with the human equivalents, using TAP, they purified human Ccr4/Not and Bet3 complexes from transfected cell cultures generated by retrovirus-mediated gene transfer. It is worthy to note that these two complexes contained large proteins (i.e., KIAA1007, KIAA1194, and KIAA1012) previously identified by our laboratory.

Taken together, these examples clearly indicate that most biological multiprotein complexes contain one or more large proteins that are believed to be essential components, both functionally and structurally. These examples are fortuitous cases, however, because it is generally difficult to purify large proteins from higher eukaryotic cells.

Because large proteins may bind both large and small proteins, interactions between large and small proteins would need to be tested with the two-hybrid screening method in our extensive research. As an initial goal, we focused on especially large protein complexes that are detectable at the electron microscopic level with the aim of examining how these complexes assemble. Because large proteins possess multiple, potential binding domains, we speculate that large proteins bind other large proteins at these binding sites to create a framework of proteins that could be characterized as a superstructure containing many scaffolds (i.e., a scaffold of scaffolds).

Concluding Remarks

Soon after the introduction of the draft sequence of the human genome, many groups initiated comprehensive functional analyses of the genome; however, if special attention is not focused on collecting data about large proteins, most of these studies will most likely accumulate information primarily about small proteins and short cDNAs. Data from our study is one step to bridge this important and technically difficult gap by providing protein-protein interaction information about large proteins.

	METHODS

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

Construction of a Prey Library Composed of Long cDNAs (KIAA Clones)

In this study, we focused on binding of the cytoplasmic domain of a large, transmembrane protein with proteins localized in the submembrane region. To minimize false-positive identification of interactions with proteins believed not to be colocalized in vivo, the nuclear and extracellular proteins located in domains outside of the cytoplasm (i.e., putative transcription factors and extracellular domains of the transmembrane protein) were removed from the prey library. This was achieved by analyzing the amino acid sequence of KIAA genes using Pfam searches and the SOSUI program. KIAA genes containing motifs for transcription factors (e.g., Pfam ID, PF00010, PF00096, PF00170, PF00249, PF00628, PF00642, PF00643, and PF01530) and the predicted transmembrane region were omitted from the library used for this study. Plasmids totaling 1087 KIAA clones (Supplemental Table 1 online) were pooled at almost equal ratios and sonicated, and blunt ends were produced using mung bean nuclease and T4 DNA polymerase (TAKARA).

DNA fragments ranging from 1 kb to 2 kb were isolated electrophoretically in agarose gels (low melting point temperature), then cloned between the DraI and EcoRV sites of the pENTR 1A vector (Invitrogen). A pool of cDNA plasmids was recovered from 3.5 × 10⁶ primary transformants (ElectroMax DH10B cells; Invitrogen) and grown on agar plates containing 30 µg/mL kanamycin. To construct the destination plasmid pACTGW-attR, plasmid attR-pSP73 (Ohara and Temple 2001) was double-digested with HindIII/SacI, and an attR cassette containing the ccdB and chloramphenicol-resistance genes was isolated and ligated into the SmaI-XhoI site of pACT2 (Clontech) with blunt ends. A pool of cDNA plasmids was converted to expression clones via the Gateway LR reaction (Invitrogen) using pACTGW-attR. The reaction products were transformed into Escherichia coli cells. Plasmids were recovered from transformants (7.3 × 10⁶ ) on ampicillin-containing agar plates after growth at 30°C for 14 h. The cDNA plasmids were transformed into yeast, Y187 (MAT). Transformants (1.6 × 10⁷ ) were grown on synthetic medium SD/Leu plates and pooled. Aliquots from the pooled yeast transformant were used for mating in two-hybrid screening.

Preparation of Bait and Mating-Mediated Yeast Two-Hybrid Screening

To construct the destination plasmid pASGW-attR, plasmid attR-pSP73 was double-digested with HindIII/SacI, and an attR cassette containing the ccdB and chloramphenicol-resistance genes was isolated and ligated into the EcoRI-SalI site of pAS2-1 (Clontech) with blunt ends. Among the KIAA proteins, KIAA genes encoding transmembrane regions were selected using predictions from the SOSUI program. Using KIAA cDNA as a template, each cDNA fragment containing a cytoplasmic domain was amplified by PCR and oligomers GWF and GWR: 5'-GGGGACAAGTTTGTACAAAAAAGCAGGCTT G(X)₂₄-3' and 5'-GGGGACCACTTTGTACAAGAAAGCTGGG TTA(X)₂₄-3', respectively, where (X)₂₄ is the gene-specific sequence. (Supplemental Table 2 online describes these cDNAs) Following the BP reaction, the attB-PCR product was introduced into expression clones via the Gateway LR reaction using pASGW-attR and a donor plasmid (Invitrogen). These clones were introduced into E. coli DH5, and the inserts of all bait plasmids were confirmed by DNA sequencing. The bait plasmid was transformed into yeast AH109 (MATa), and transformants were grown on SD/Trp plates. Self-activation assays were performed on SD/Ade/His/Trp plates. Two transformants per each bait exhibiting no background growth on Ade/His were pooled, cultured, and mated with an aliquot of yeast prey containing the transformant library. Diploid cells were spread onto 9 cm × 13 cm 20-25 plates containing SD/Ade/His/Leu/Trp/15 mM 3-AT (3-amino-1,2,4-triazole). After 7-10 d of growth at 30°C, typically 30 colonies were replated onto SD/Ade/His/Leu/Trp/15 mM 3-AT plates containing X--gal. Blue colonies were selected and cultured in SD/Leu medium. Each plasmid was purified from the yeast and used as a template for PCR with oligomers OADN (5'-TCG ATG ATG AAG ATA CCC CAC-3') and OADC (5'-AAG AAA TTG AGA TGG TGC ACG-3'). DNA sequences of prey clones were determined to identify the interacting region of the KIAA protein. The resulting sequences were compared with the HUGE database (http://www.kazusa.or.jp/huge/) using BLAST2.

Exchanging of Bait and Prey Sequences and Confirmation of Binding Activity by Retest

To confirm the binding activity detected by the library screening of yeast two-hybrids, bait and prey sequences were exchanged using the Gateway reaction, and both were transfected into haploid AH109 yeasts. Baits were prepared by introducing the PCR product (from OADN and OADC oligomers) into expression clones via the Gateway BP and LR reaction and pASGW-attR and donor plasmids. Preys were similarly prepared by introducing the PCR product (from GWF and GWR oligomers) into expression clones using the Gateway BP and LR reaction and pACTGW-attR and donor plasmids. Gateway reaction products were transformed into E. coli DH5. After purification of the plasmids, the new bait and prey plasmids or the empty vector were transformed into yeast AH109; viability was checked on SD/Leu/Trp and SD/Ade/His/Leu/Trp plates.

	WEB SITE REFERENCES

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

http://www.kazusa.or.jp/huge/; HUGE database.

http://www.kazusa.or.jp/huge/ppi/; detailed accounts of 84 protein-protein interactions.