A general model for preferential hetero-oligomerization of LIN-2/7 domains: mechanism underlying directed assembly of supramolecular signaling complexes

LIN-2/7 (L27) domains are protein interaction modules that preferentially hetero-oligomerize, a property critical for their function in directing specific assembly of supramolecular signaling complexes at synapses and other polarized cell-cell junctions. We have solved the solution structure of the heterodimer composed of the L27 domains from LIN-2 and LIN-7. Comparison of this structure with other L27 domain structures has allowed us to formulate a general model for why most L27 domains form an obligate heterodimer complex. L27 domains can be divided in two types (A and B), with each heterodimer comprising an A/B pair. We have identified two keystone positions that play a central role in discrimination. The residues at these positions are energetically acceptable in the context of an A/B heterodimer, but would lead to packing defects or electrostatic repulsion in the context of A/A and B/B homodimers. As predicted by the model, mutations of keystone residues stabilize normally strongly disfavored homodimers. Thus, L27 domains are specifically optimized to avoid homodimeric interactions.     

Identification of multiple binding partners for the amino-terminal domain of synapse-associated protein 97

Multiprotein complexes mediate static and dynamic functions to establish and maintain cell polarity in both epithelial cells and neurons. Membrane-associated guanylate kinase (MAGUK) proteins are thought to be scaffolding molecules in these processes and bind multiple proteins via their obligate postsynaptic density (PSD)-95/Disc Large/Zona Occludens-1, Src homology 3, and guanylate kinase-like domains. Subsets of MAGUK proteins have additional protein-protein interaction domains. An additional domain we identified in SAP97 called the MAGUK recruitment (MRE) domain binds the LIN-2,7 amino-terminal (L27N) domain of mLIN-2/CASK, a MAGUK known to bind mLIN-7. Here we show that SAP97 binds two other mLIN-7 binding MAGUK proteins. One of these MAGUK proteins, DLG3, coimmunoprecipitates with SAP97 in lysates from rat brain and transfected Madin-Darby canine kidney cells. This interaction requires the MRE domain of SAP97 and surprisingly, both the L27N and L27 carboxyl-terminal (L27C) domains of DLG3. We also demonstrate that SAP97 can interact with the MAGUK protein, DLG2, but not the highly related protein, PALS2. The ability of SAP97 to interact with multiple MAGUK proteins is likely to be important for the targeting of specific protein complexes in polarized cells.     

Structural basis for L27 domain-mediated assembly of signaling and cell polarity complexes

L27 is a protein-binding domain that can assemble essential proteins for signaling and cell polarity into complexes by interacting in a heterodimeric manner. One of these protein complexes is the PATJ/PALS1/Crumbs tripartite complex, which is crucial for the establishment and maintenance of cell polarity. To reveal the structural basis underlining the obligate heterodimerization, we have determined the crystal structure of the PALS1-L27N/PATJ-L27 heterodimer complex. Each L27 domain is composed of three helices. The two L27 domains heterodimerize by building a compact structure consisting of a four-helix bundle formed by the first two helices of each L27 domain and one coiled-coil formed by the third helix of each domain. The large hydrophobic packing interactions contributed by all the helices of both L27 domains predominantly drive the heterodimer formation, which is likely to be a general feature of L27 domains. Combined with mutational studies, we can begin to understand the structural basis for the specificity of L27 binding pairs. Our results provide unique insights into L27 domain heterodimer complex, which is critical for cell polarization.     

PDZ-mediated interactions retain the epithelial GABA transporter on the basolateral surface of polarized epithelial cells

The PDZ target motifs located in the C-terminal end of many receptors and ion channels mediate protein-protein interactions by binding to specific PDZ-containing proteins. These interactions are involved in the localization of surface proteins on specialized membrane domains of neuronal and epithelial cells. However, the molecular mechanism responsible for this PDZ protein-dependent polarized localization is still unclear. This study first demonstrated that the epithelial gamma-aminobutyric acid (GABA) transporter (BGT-1) contains a PDZ target motif that mediates the interaction with the PDZ protein LIN-7 in Madin-Darby canine kidney (MDCK) cells, and then investigated the role of this interaction in the basolateral localization of the transporter. It was found that although the transporters from which the PDZ target motif was deleted were still targeted to the basolateral surface, they were not retained but internalized in an endosomal recycling compartment. Furthermore, an interfering BGT peptide determined the intracellular relocation of the native transporter. These data indicate that interactions with PDZ proteins determine the polarized surface localization of target proteins by means of retention and not targeting mechanisms. PDZ proteins may, therefore, act as a sort of membrane protein sorting machinery which, by recognizing retention signals (the PDZ target sequences), prevents protein internalization.     

LIN-10 can promote LET-23 EGFR signaling and trafficking independently of LIN-2 and LIN-7

During Caenorhabditis elegans larval development, an inductive signal mediated by the LET-23 EGFR (epidermal growth factor receptor), specifies three of six vulva precursor cells (VPCs) to adopt vulval cell fates. An evolutionarily conserved complex consisting of PDZ domain-containing scaffold proteins LIN-2 (CASK), LIN-7 (Lin7 or Veli), and LIN-10 (APBA1 or Mint1) (LIN-2/7/10) mediates basolateral LET-23 EGFR localization in the VPCs to permit signal transmission and development of the vulva. We recently found that the LIN-2/7/10 complex likely forms at Golgi ministacks; however, the mechanism through which the complex targets the receptor to the basolateral membrane remains unknown. Here we found that overexpression of LIN-10 or LIN-7 can compensate for loss of their complex components by promoting LET-23 EGFR signaling through previously unknown complex-independent and receptor-dependent pathways. In particular, LIN-10 can independently promote basolateral LET-23 EGFR localization, and its complex-independent function uniquely requires its PDZ domains that also regulate its localization to Golgi. These studies point to a novel complex-independent function for LIN-7 and LIN-10 that broadens our understanding of how this complex regulates targeted sorting of membrane proteins.     

Linking up and interacting with BRCT domains

BRCT domains are present in an ever expanding family of proteins that includes many DNA repair and checkpoint proteins. The most prominent member of the BRCT family is BRCA1, mutations in which are responsible for a high proportion of breast and ovarian cancers. BRCT domains act as protein–protein interaction modules and facilitate the formation of hetero- and homo-oligomers. The domains occur either singly or in pairs, with up to eight domains in a single protein. When in pairs the domains are separated by a short inter-BRCT linker. Numerous crystal structures have been determined for BRCT domains from a range of different proteins, which indicate that the overall structure of the BRCT domains is generally well conserved. In contrast, the positions and structures of the linker regions are more varied, as are the roles of the linkers. Here, we describe the protein–protein interactions involving three different inter-BRCT linker regions, those of DNA ligase IV (LigIV), Schizosaccharomyces pombe Crb2 and human 53BP1.     

Emergence of scale-free distribution in protein–protein interaction networks based on random selection of interacting domain pairs

Recent analyses of biological and artificial networks have revealed a common network architecture, called scale-free topology. The origin of the scale-free topology has been explained by using growth and preferential attachment mechanisms. In a cell, proteins are the most important carriers of function, and are composed of domains as elemental units responsible for the physical interaction between protein pairs. Here, we propose a model for protein–protein interaction networks that reveals the emergence of two possible topologies. We show that depending on the number of randomly selected interacting domain pairs, the connectivity distribution follows either a scale-free distribution, even in the absence of the preferential attachment, or a normal distribution. This new approach only requires an evolutionary model of proteins (nodes) but not for the interactions (edges). The edges are added by means of random interaction of domain pairs. As a result, this model offers a new mechanistic explanation for understanding complex networks with a direct biological interpretation because only protein structures and their functions evolved through genetic modifications of amino acid sequences. These findings are supported by numerical simulations as well as experimental data.     

Identification of critical domains and putative partners for the<Emphasis Type="Italic"> Caenorhabditis elegans</Emphasis> spindle component LIN-5

Successful cell division requires proper assembly, placement and functioning of the spindle apparatus that segregates the chromosomes. The Caenorhabditis elegans gene lin-5 encodes a novel coiled-coil component of the spindle required for spindle positioning and chromosome segregation. To gain further insights into lin-5 function, we screened for dominant suppressors of the partial loss-of-function phenotype associated with the mutation lin-5(ev571ts), and isolated 68 suppressing mutations. Eight out of the ten suppressors sequenced contained intragenic missense mutations immediately upstream of the lesion in lin-5(ev571ts). These probably help to stabilize protein-protein interactions mediated by the coiled-coil domain. This domain was found to be required for binding to several putative LIN-5 interacting (LFI) proteins identified in yeast two-hybrid screens. Interestingly, interaction with the coiled-coil protein LFI-1 was specifically reduced by the lin-5(ev571ts) mutation and restored by a representative intragenic suppressor mutation. Immunostaining experiments showed that LIN-5 and LFI-1 may co-localize around the kinetochore microtubules during metaphase, indicating potential interaction in vivo. The coiled-coil domain of LIN-5 was also found to mediate homodimerization, while the C-terminal region of LIN-5 was sufficient for interaction with GPR-1, a recently identified component of a LIN-5 spindle-regulatory complex. A single amino-acid substitution in the N-terminal region of LIN-5, encoded by the e1457 allele, abolished all LIN-5 interactions. Taken together, our results indicate that the spindle functions of LIN-5 depend on interactions with multiple protein partners, and that these interactions are mediated through several different domains of LIN-5.Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by C. P. Hollenberg     

Statistical analysis of domains in interacting protein pairs

MOTIVATION: Several methods have recently been developed to analyse large-scale sets of physical interactions between proteins in terms of physical contacts between the constituent domains, often with a view to predicting new pairwise interactions. Our aim is to combine genomic interaction data, in which domain-domain contacts are not explicitly reported, with the domain-level structure of individual proteins, in order to learn about the structure of interacting protein pairs. Our approach is driven by the need to assess the evidence for physical contacts between domains in a statistically rigorous way. RESULTS: We develop a statistical approach that assigns p-values to pairs of domain superfamilies, measuring the strength of evidence within a set of protein interactions that domains from these superfamilies form contacts. A set of p-values is calculated for SCOP superfamily pairs, based on a pooled data set of interactions from yeast. These p-values can be used to predict which domains come into contact in an interacting protein pair. This predictive scheme is tested against protein complexes in the Protein Quaternary Structure (PQS) database, and is used to predict domain-domain contacts within 705 interacting protein pairs taken from our pooled data set.     

The origins of specificity in polyketide synthase protein interactions

Polyketides, a diverse group of heteropolymers with antibiotic and antitumor properties, are assembled in bacteria by multiprotein chains of modular polyketide synthase (PKS) proteins. Specific protein-protein interactions determine the order of proteins within a multiprotein chain, and thereby the order in which chemically distinct monomers are added to the growing polyketide product. Here we investigate the evolutionary and molecular origins of protein interaction specificity. We focus on the short, conserved N- and C-terminal docking domains that mediate interactions between modular PKS proteins. Our computational analysis, which combines protein sequence data with experimental protein interaction data, reveals a hierarchical interaction specificity code. PKS docking domains are descended from a single ancestral interacting pair, but have split into three phylogenetic classes that are mutually noninteracting. Specificity within one such compatibility class is determined by a few key residues, which can be used to define compatibility subclasses. We identify these residues using a novel, highly sensitive co-evolution detection algorithm called CRoSS (correlated residues of statistical significance). The residue pairs selected by CRoSS are involved in direct physical interactions in a docked-domain NMR structure. A single PKS system can use docking domain pairs from multiple classes, as well as domain pairs from multiple subclasses of any given class. The termini of individual proteins are frequently shuffled, but docking domain pairs straddling two interacting proteins are linked as an evolutionary module. The hierarchical and modular organization of the specificity code is intimately related to the processes by which bacteria generate new PKS pathways.     

aPKC phosphorylates NuMA-related LIN-5 to position the mitotic spindle during asymmetric division

The position of the mitotic spindle controls the plane of cell cleavage and determines whether polarized cells divide symmetrically or asymmetrically. In animals, an evolutionarily conserved pathway of LIN-5 (homologues: Mud and NuMA), GPR-1/2 (homologues: Pins, LGN, AGS-3) and Gα mediates spindle positioning, and acts downstream of the conserved PAR-3-PAR-6-aPKC polarity complex. However, molecular interactions between polarity proteins and LIN-5-GPR-Gα remain to be identified. Here we describe a quantitative mass spectrometry approach for in vivo identification of protein kinase substrates. Applying this strategy to Caenorhabditis elegans embryos, we found that depletion of the polarity kinase PKC-3 results in markedly decreased levels of phosphorylation of a cluster of four LIN-5 serine residues. These residues are directly phosphorylated by PKC-3 in vitro. Phospho-LIN-5 co-localizes with PKC-3 at the anterior cell cortex and temporally coincides with a switch from anterior- to posterior-directed spindle movements in the one-cell embryo. LIN-5 mutations that prevent phosphorylation increase the extent of anterior-directed spindle movements, whereas phosphomimetic mutations decrease spindle migration. Our results indicate that anterior-located PKC-3 inhibits cortical microtubule pulling forces through direct phosphorylation of LIN-5. This molecular interaction between polarity and spindle-positioning proteins may be used broadly in cell cleavage plane determination.     

H3K9me2/3 binding of the MBT domain protein LIN-61 is essential for Caenorhabditis elegans vulva development

MBT domain proteins are involved in developmental processes and tumorigenesis. In vitro binding and mutagenesis studies have shown that individual MBT domains within clustered MBT repeat regions bind mono- and dimethylated histone lysine residues with little to no sequence specificity but discriminate against the tri- and unmethylated states. However, the exact function of promiscuous histone methyl-lysine binding in the biology of MBT domain proteins has not been elucidated. Here, we show that the Caenorhabditis elegans four MBT domain protein LIN-61, in contrast to other MBT repeat factors, specifically interacts with histone H3 when methylated on lysine 9, displaying a strong preference for di- and trimethylated states (H3K9me2/3). Although the fourth MBT repeat is implicated in this interaction, H3K9me2/3 binding minimally requires MBT repeats two to four. Further, mutagenesis of residues conserved with other methyl-lysine binding MBT regions in the fourth MBT repeat does not abolish interaction, implicating a distinct binding mode. In vivo, H3K9me2/3 interaction of LIN-61 is required for C. elegans vulva development within the synMuvB pathway. Mutant LIN-61 proteins deficient in H3K9me2/3 binding fail to rescue lin-61 synMuvB function. Also, previously identified point mutant synMuvB alleles are deficient in H3K9me2/3 interaction although these target residues that are outside of the fourth MBT repeat. Interestingly, lin-61 genetically interacts with two other synMuvB genes, hpl-2, an HP1 homologous H3K9me2/3 binding factor, and met-2, a SETDB1 homologous H3K9 methyl transferase (H3K9MT), in determining C. elegans vulva development and fertility. Besides identifying the first sequence specific and di-/trimethylation binding MBT domain protein, our studies imply complex multi-domain regulation of ligand interaction of MBT domains. Our results also introduce a mechanistic link between LIN-61 function and biology, and they establish interplay of the H3K9me2/3 binding proteins, LIN-61 and HPL-2, as well as the H3K9MT MET-2 in distinct developmental pathways.     

Structure of an L27 domain heterotrimer from cell polarity complex Patj/Pals1/Mals2 reveals mutually independent L27 domain assembly mode

The assembly of supramolecular complexes in multidomain scaffold proteins is crucial for the control of cell polarity. The scaffold protein of protein associated with Lin-7 1 (Pals1) forms a complex with two other scaffold proteins, Pals-associated tight junction protein (Patj) and mammalian homolog-2 of Lin-7 (Mals2), through its tandem Lin-2 and Lin-7 (L27) domains to regulate apical-basal polarity. Here, we report the crystal structure of a 4-L27 domain-containing heterotrimer derived from the tripartite complex Patj/Pals1/Mals2. The heterotrimer consists of two cognate pairs of heterodimeric L27 domains with similar conformations. Structural analysis and biochemical data further show that the dimers assemble mutually independently. Additionally, such mutually independent assembly of the two heterodimers can be observed in another tripartite complex, Disks large homolog 1 (DLG1)/calcium-calmodulin-dependent serine protein kinase (CASK)/Mals2. Our results reveal a novel mechanism for tandem L27 domain-mediated, supramolecular complex assembly with a mutually independent mode.     

Conservation of protein-protein interactions - lessons from ascomycota

Interacting proteins from Saccharomyces cerevisiae are evolutionarily conserved and their likelihood of having an ortholog in other ascomycota species correlates with the number of interaction partners. Moreover, interacting proteins show a clear preference to be conserved as a pair, indicating that nature maintains selection pressure on the interaction links between proteins. The conservation of interacting protein pairs between different organisms does not exhibit any bias with respect to protein functional roles.     

Co-evolutionary analysis of domains in interacting proteins reveals insights into domain-domain interactions mediating protein-protein interactions

Recent advances in functional genomics have helped generate large-scale high-throughput protein interaction data. Such networks, though extremely valuable towards molecular level understanding of cells, do not provide any direct information about the regions (domains) in the proteins that mediate the interaction. Here, we performed co-evolutionary analysis of domains in interacting proteins in order to understand the degree of co-evolution of interacting and non-interacting domains. Using a combination of sequence and structural analysis, we analyzed protein-protein interactions in F1-ATPase, Sec23p/Sec24p, DNA-directed RNA polymerase and nuclear pore complexes, and found that interacting domain pair(s) for a given interaction exhibits higher level of co-evolution than the non-interacting domain pairs. Motivated by this finding, we developed a computational method to test the generality of the observed trend, and to predict large-scale domain-domain interactions. Given a protein-protein interaction, the proposed method predicts the domain pair(s) that is most likely to mediate the protein interaction. We applied this method on the yeast interactome to predict domain-domain interactions, and used known domain-domain interactions found in PDB crystal structures to validate our predictions. Our results show that the prediction accuracy of the proposed method is statistically significant. Comparison of our prediction results with those from two other methods reveals that only a fraction of predictions are shared by all the three methods, indicating that the proposed method can detect known interactions missed by other methods. We believe that the proposed method can be used with other methods to help identify previously unrecognized domain-domain interactions on a genome scale, and could potentially help reduce the search space for identifying interaction sites.     

The binding affinities of proteins interacting with the PDZ domain of PICK1

Protein interacting with C kinase (PICK1) is well conserved throughout evolution and plays a critical role in synaptic plasticity by regulating the trafficking and posttranslational modification of its interacting proteins. PICK1 contains a single PSD95/DlgA/Zo-1 (PDZ) protein-protein interaction domain, which is promiscuous and shown to interact with over 60 proteins, most of which play roles in neuronal function. Several reports have suggested the role of PICK1 in disorders such as epilepsy, pain, brain trauma and stroke, drug abuse and dependence, schizophrenia and psychosis. Importantly, lead compounds that block PICK1 interactions are also now becoming available. Here, a new modeling approach was developed to investigate binding affinities of PDZ interactions. Using these methods, the binding affinities of all major PICK1 interacting proteins are reported and the effects of PICK1 mutations on these interactions are described. These modeling methods have important implications in defining the binding properties of proteins interacting with PICK1 as well as the general structural requirements of PDZ interactions. The study also provides modeling methods to support in the drug design of ligands for PDZ domains, which may further aid in development of the family of PDZ domains as a drug target.     

Prediction of protein-protein interactions using protein signature profiling

Protein domains are conserved and functionally independent structures that play an important role in interactions among related proteins. Domain-domain interactions have been recently used to predict protein-protein interactions (PPI). In general, the interaction probability of a pair of domains is scored using a trained scoring function. Satisfying a threshold, the protein pairs carrying those domains are regarded as "interacting". In this study, the signature contents of proteins were utilized to predict PPI pairs in Saccharomyces cerevisiae, Caenorhabditis elegans, and Homo sapiens. Similarity between protein signature patterns was scored and PPI predictions were drawn based on the binary similarity scoring function. Results show that the true positive rate of prediction by the proposed approach is approximately 32% higher than that using the maximum likelihood estimation method when compared with a test set, resulting in 22% increase in the area under the receiver operating characteristic (ROC) curve. When proteins containing one or two signatures were removed, the sensitivity of the predicted PPI pairs increased significantly. The predicted PPI pairs are on average 11 times more likely to interact than the random selection at a confidence level of 0.95, and on average 4 times better than those predicted by either phylogenetic profiling or gene expression profiling.     

A direct interaction between the carboxyl-terminal region of CDC5L and the WD40 domain of PLRG1 is essential for pre-mRNA splicing.

The human proteins CDC5L (hCDC5) and PLRG1 are both highly conserved components of a multiprotein complex that is a subunit of the spliceosome. The respective homologues in yeast of both proteins are also associated with a sub-spliceosomal multiprotein complex that has been shown to be important for pre-mRNA splicing. We show that these two human proteins are associated in vivo and will interact directly in vitro. The regions containing the interacting domains in both proteins have been identified. Our results indicate that the carboxyl-terminal region of CDC5L and the WD40 domain of PLRG1 are essential for direct interaction between both proteins. By using a bacterially expressed mutant protein, containing the PLRG1 interacting domain in CDC5L, we show that the CDC5L-PLRG1 interaction in HeLa nuclear extract can be disrupted causing pre-mRNA splicing to be inhibited. Thus, a direct interaction between the CDC5L protein and PLRG1 in the CDC5L complex is essential for pre-mRNA splicing progression.     

Identification of peptide inhibitors of pre-mRNA splicing derived from the essential interaction domains of CDC5L and PLRG1

CDC5L and PLRG1 are both spliceosomal proteins that are highly conserved across species. They have both been shown to be part of sub- spliceosomal protein complexes that are essential for pre-mRNA splicing in yeast and humans. CDC5L and PLRG1 interact directly in vitro. This interaction is mediated by WD40 regions in PLRG1 and the C-terminal domain of CDC5L. In order to determine whether this interaction is important for the splicing mechanism, we have designed peptides corresponding to highly conserved sequences in the interaction domains of both proteins. These peptides were used in in vitro splicing experiments as competitors to the cognate sequences in the endogenous proteins. Certain peptides derived from the binding domains of both proteins were found to inhibit in vitro splicing. This splicing inhibition could be prevented by preincubating the peptides with the corresponding partner protein that had been expressed in Escherichia coli. The results from this study indicate that the interaction between CDC5L and PLRG1 is essential for pre-mRNA splicing and further demonstrate that small peptides can be used as effective splicing inhibitors.     

Ribosomal proteins L7 and L8 function in concert with six A3 assembly factors to propagate assembly of domains I and II of 25S rRNA in yeast 60S ribosomal subunits

Ribosome biogenesis is a complex multistep process that involves alternating steps of folding and processing of pre-rRNAs in concert with assembly of ribosomal proteins. Recently, there has been increased interest in the roles of ribosomal proteins in eukaryotic ribosome biogenesis in vivo, focusing primarily on their function in pre-rRNA processing. However, much less is known about participation of ribosomal proteins in the formation and rearrangement of preribosomal particles as they mature to functional subunits. We have studied ribosomal proteins L7 and L8, which are required for the same early steps in pre-rRNA processing during assembly of 60S subunits but are located in different domains within ribosomes. Depletion of either leads to defects in processing of 27SA(3) to 27SB pre-rRNA and turnover of pre-rRNAs destined for large ribosomal subunits. A specific subset of proteins is diminished from these residual assembly intermediates: six assembly factors required for processing of 27SA(3) pre-rRNA and four ribosomal proteins bound to domain I of 25S and 5.8S rRNAs surrounding the polypeptide exit tunnel. In addition, specific sets of ribosomal proteins are affected in each mutant: In the absence of L7, proteins bound to domain II, L6, L14, L20, and L33 are greatly diminished, while proteins L13, L15, and L36 that bind to domain I are affected in the absence of L8. Thus, L7 and L8 might establish RNP structures within assembling ribosomes necessary for the stable association and function of the A(3) assembly factors and for proper assembly of the neighborhoods containing domains I and II.