Protein interaction networks--more than mere modules

It is widely believed that the modular organization of cellular function is reflected in a modular structure of molecular networks. A common view is that a “module” in a network is a cohesively linked group of nodes, densely connected internally and sparsely interacting with the rest of the network. Many algorithms try to identify functional modules in protein-interaction networks (PIN) by searching for such cohesive groups of proteins. Here, we present an alternative approach independent of any prior definition of what actually constitutes a “module”. In a self-consistent manner, proteins are grouped into “functional roles” if they interact in similar ways with other proteins according to their functional roles. Such grouping may well result in cohesive modules again, but only if the network structure actually supports this. We applied our method to the PIN from the Human Protein Reference Database (HPRD) and found that a representation of the network in terms of cohesive modules, at least on a global scale, does not optimally represent the network's structure because it focuses on finding independent groups of proteins. In contrast, a decomposition into functional roles is able to depict the structure much better as it also takes into account the interdependencies between roles and even allows groupings based on the absence of interactions between proteins in the same functional role. This, for example, is the case for transmembrane proteins, which could never be recognized as a cohesive group of nodes in a PIN. When mapping experimental methods onto the groups, we identified profound differences in the coverage suggesting that our method is able to capture experimental bias in the data, too. For example yeast-two-hybrid data were highly overrepresented in one particular group. Thus, there is more structure in protein-interaction networks than cohesive modules alone and we believe this finding can significantly improve automated function prediction algorithms.

Cellular function is widely believed to be organized in a modular fashion. On all scales and at all levels of complexity, relatively independent sub-units perform relatively independent sub-tasks. This functional modularity must be reflected in the topology of molecular networks. But how a functional module should be represented in an interaction network is an open question. On a small scale, one can identify a protein-complex as a module in protein-interaction networks (PIN), i.e., modules are understood as densely linked (interacting) groups of proteins, that are only sparsely interacting with the rest of the network. In this contribution, we show that extrapolating this concept of cohesively linked clusters of proteins as modules to the scale of the entire PIN inevitably misses important and functionally relevant structure inherent in the network. As an alternative, we introduce a novel way of decomposing a network into functional roles and show that this represents network structure and function more efficiently. This finding should have a profound impact on all module assisted methods of protein function prediction and should shed new light on how functional modules can be represented in molecular interaction networks in general.

Activator Gcn4 employs multiple segments of Med15/Gal11, including the KIX domain, to recruit mediator to target genes in vivo

Mediator is a multisubunit coactivator required for initiation by RNA polymerase II. The Mediator tail subdomain, containing Med15/Gal11, is a target of the activator Gcn4 in vivo, critical for recruitment of native Mediator or the Mediator tail subdomain present in sin4Delta cells. Although several Gal11 segments were previously shown to bind Gcn4 in vitro, the importance of these interactions for recruitment of Mediator and transcriptional activation by Gcn4 in cells was unknown. We show that interaction of Gcn4 with the Mediator tail in vitro and recruitment of this subcomplex and intact Mediator to the ARG1 promoter in vivo involve additive contributions from three different segments in the N terminus of Gal11. These include the KIX domain, which is a critical target of other activators, and a region that shares a conserved motif (B-box) with mammalian coactivator SRC-1, and we establish that B-box is a critical determinant of Mediator recruitment by Gcn4. We further demonstrate that Gcn4 binds to the Gal11 KIX domain directly and, by NMR chemical shift analysis combined with mutational studies, we identify the likely binding site for Gcn4 on the KIX surface. Gcn4 is distinctive in relying on comparable contributions from multiple segments of Gal11 for efficient recruitment of Mediator in vivo.

Molecular basis for association of PIPKI gamma-p90 with clathrin adaptor AP-2

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P(2)) is an essential determinant in clathrin-mediated endocytosis (CME). In mammals three type I phosphatidylinositol-4-phosphate 5-kinase (PIPK) enzymes are expressed, with the I gamma-p90 isoform being highly expressed in the brain where it regulates synaptic vesicle (SV) exo-/endocytosis at nerve terminals. How precisely PI(4,5)P(2) metabolism is controlled spatially and temporally is still uncertain, but recent data indicate that direct interactions between type I PIPK and components of the endocytic machinery, in particular the AP-2 adaptor complex, are involved. Here we demonstrated that PIPKI gamma-p90 associates with both the mu and beta2 subunits of AP-2 via multiple sites. Crystallographic data show that a peptide derived from the splice insert of the human PIPKI gamma-p90 tail binds to a cognate recognition site on the sandwich subdomain of the beta2 appendage. Partly overlapping aromatic and hydrophobic residues within the same peptide also can engage the C-terminal sorting signal binding domain of AP-2mu, thereby potentially competing with the sorting of conventional YXXØ motif-containing cargo. Biochemical and structure-based mutagenesis analysis revealed that association of the tail domain of PIPKI gamma-p90 with AP-2 involves both of these sites. Accordingly the ability of overexpressed PIPKI gamma tail to impair endocytosis of SVs in primary neurons largely depends on its association with AP-2 beta and AP-2mu. Our data also suggest that interactions between AP-2 and the tail domain of PIPKI gamma-p90 may serve to regulate complex formation and enzymatic activity. We postulate a model according to which multiple interactions between PIPKI gamma-p90 and AP-2 lead to spatiotemporally controlled PI(4,5)P(2) synthesis during clathrin-mediated SV endocytosis.

Identifying the functional part of heparin-binding protein (HBP) as a monocyte stimulator and the novel role of monocytes as HBP producers

Heparin-binding protein (HBP), an evolutionary ancient and biologically highly important molecule in inflammation, is an inactive serine protease due to mutations in the catalytic triad. The histidine (position 41) in the conserved sequence TAAHC is mutated to serine and this sequence (TAASC) plays a crucial role when HBP binds to monocytes. We synthesized a 20-44 HBP peptide, cyclicized by a sulphur bridge, which encompasses this amino acid and functions as full-length HBP. Using a human monocyte cell line, we have shown that lipopolysaccharide (LPS)-triggered secretion of IL-6 is enhanced up to 10-fold when full-length HBP or the peptide are present in low-to-moderate concentrations. A monoclonal antibody neutralizing HBP also neutralizes the peptide, indicating that the ligand for the HBP receptor is located near serine in position 41 on the HBP surface. A 'back mutated' 20-44 peptide (serine→histidine) has some, but not significant, stimulatory effect on monocytes. Normally, HBP production and release is ascribed to neutrophil granulocytes, but here we find that also monocytes secrete HBP when stimulated with LPS. Furthermore, a small amount of HBP can be demonstrated when monocytes are incubated in medium alone. Our efforts to identify a suggested HBP receptor on monocytes has failed so far.

Involvement of three glutamine tracts in human androgen receptor transactivation

The androgen receptor (AR) possesses a polymorphic polyglutamine tract (polyQ), whose length is inversely correlated with its transcriptional activity. Here, we investigated whether 6 and 5 repetitive glutamine tracts (Q6 and Q5, respectively) in the N-terminal domain of AR also have effects on AR transactivation. In a reporter gene assay using two-tandem repeats of an androgen response element, deletion of glutamine tracts significantly increased AR transactivation in the following order: wild-type<a single deletion of polyQ or Q5<double deletion of polyQ and Q6<double deletion of polyQ and Q5<triple deletion. Deletion of polyQ alone or combined deletion of polyQ and Q5 from an AR mutant lacking the ligand-binding domain, which is constitutively active due to activation function-1, increased AR transactivation. However, the glutamine tracts had no influence on activation function-1 activity, suggesting that the glutamine tracts modulate the binding of AR to DNA. Q5, like polyQ, was found to be involved in the interaction between the NH(2)- and COOH-terminal regions of AR (N-C interaction). These results indicate that the inhibitory effects of polyQ and Q5 on AR transactivation are the due, at least in part, to their negative regulation of N-C interaction.

H2O2-dependent translocation of TCTP into the nucleus enables its interaction with VDR in human keratinocytes: TCTP as a further module in calcitriol signalling

Translationally controlled tumour protein (TCTP) is an evolutionarily highly conserved molecule implicated in many processes related to cell cycle progression, proliferation and growth, to the protection against harmful conditions including apoptosis and to the human allergic response. We are showing here that after application of mild oxidative stress, human TCTP relocates from the cytoplasm to the nuclei of HaCaT keratinocytes where it directly associates with the ligand-binding domain of endogenous vitamin D(3) receptor (VDR) through its helical domain 2 (AA 71-132). Interestingly, the latter harbours a putative nuclear hormone receptor coregulatory LxxLL-like motif which seems to be involved in the interaction. Moreover, we demonstrate that VDR transcriptionally induces the expression of TCTP by binding to a previously unknown VDR response element within the TCTP promotor. Conversely, ectopically overexpressed TCTP downregulates the amount of VDR on both mRNA as well as protein level. These data, to conclude, suggest a kind of feedback regulation between TCTP and VDR to regulate a variety of (Ca(2+) dependent) cellular effects and in this way further underscore the physiological relevance of this novel protein-protein interaction.

Structure of protein interaction networks and their implications on drug design

Protein-protein interaction networks (PINs) are rich sources of information that enable the network properties of biological systems to be understood. A study of the topological and statistical properties of budding yeast and human PINs revealed that they are scale-rich and configured as highly optimized tolerance (HOT) networks that are similar to the router-level topology of the Internet. This is different from claims that such networks are scale-free and configured through simple preferential-attachment processes. Further analysis revealed that there are extensive interconnections among middle-degree nodes that form the backbone of the networks. Degree distributions of essential genes, synthetic lethal genes, synthetic sick genes, and human drug-target genes indicate that there are advantageous drug targets among nodes with middle- to low-degree nodes. Such network properties provide the rationale for combinatorial drugs that target less prominent nodes to increase synergetic efficacy and create fewer side effects.

Genome-wide data on interactions between proteins are now available, and networks of protein interactions are the keys to understanding diseases and finding accurate drug targets. This study revealed that the architectural properties of the backbones of protein interaction networks (PINs) were similar to those of the Internet router-level topology by using statistical analyses of genome-wide budding yeast and human PINs. This type of network is known as a highly optimized tolerance (HOT) network that is robust against failures in its components and that ensures high levels of communication. Moreover, we also found that a large number of the most successful drug-target proteins are on the backbone of the human PIN. We made a list of proteins on the backbone of the human PIN, which may help drug companies to search more efficiently for new drug targets.

RUNX3 directly interacts with intracellular domain of Notch1 and suppresses Notch signaling in hepatocellular carcinoma cells

RUNX3 takes a strong suppressive effect in many tumors including hepatocellular carcinoma (HCC). HES-1, a downstream target of Notch signaling, is shown to be decreased in human HCC cell line SMMC7721 with RUNX3 gene transfection. Since Notch signaling is oncogenic in HCC, RUNX3 might exert its inhibitory effect in HCC partly through the suppression on Notch signaling. To investigate the possible mechanism of the down-regulation of HES-1 by RUNX3, we performed Western blot and reporter assay and found that RUNX3 suppressed intracellular domain of Notch1 (ICN1)-mediated transactivation of Notch signaling while it did not alter the expression of ICN1 and recombination signal binding protein-J kappa (RBP-J) in SMMC7721 cells. Besides, confocal microscopy, co-immunoprecipitation and GST pull-down assays showed that RUNX3 could co-localize with ICN1 and RBP-J, forming a complex with these two molecules in nucleus of SMMC7721 cells by its direct interaction with ICN1. Furthermore, RUNX3 was recruited to RBP-J recognition motif of HES-1 promoter, which was identified by chromatin immunoprecipitation (ChIP) and electrophoretic mobility shift assay (EMSA). Taken together, these findings indicate that RUNX3 suppresses Notch signaling in HCC SMMC7721 cells by its interaction with ICN1 and thus recruitment to the RBP-J recognition motif of downstream genes of Notch signaling.

Interaction of heparin and heparin-derived oligosaccharides with synthetic peptide analogues of the heparin-binding domain of heparin/heparan sulfate-interacting protein

BACKGROUND

Although protamine is effective as an antidote of heparin, there is a need to replace protamine due to its side effects. HIP peptide has been reported to neutralize the anticoagulant activity of heparin. The interaction of HIP analog peptides with heparin and heparin-derived oligosaccharides is investigated in this paper.

METHODS

Seven analogues of the heparin-binding domain of heparin/heparan sulfate-interacting protein (HIP) were synthesized, and their interaction with heparin was characterized by heparin affinity chromatography, isothermal titration calorimetry, and NMR.

RESULTS

NMR results indicate the imidazolium groups of the His side chains of histidine-containing Hip analog peptide interact site-specifically with heparin at pH 5.5. Heparin has identical affinities for HIP analog peptides of opposite chirality. Analysis by counterion condensation theory indicates the peptide AC-SRPKAKAKAKAKDQTK-NH2 makes on average approximately 3 ionic interactions with heparin that result in displacement of approximately 2 Na+ ions, and ionic interactions account for approximately 46% of the binding free energy at a Na+ concentration of 0.15 M.

CONCLUSIONS

The affinity of heparin for the peptides is strongly dependent on the nature of the cationic side chains and pH. The thermodynamic parameters measured for the interaction of HIP peptide analogs with heparin are strongly dependent on the peptide sequence and pH.

GENERAL SIGNIFICANCE

The information obtained in this research will be of use in the design of new agents for neutralization of the anticoagulant activity of heparin. The site-specific binding of protonated histidine side chains to heparin provides a molecular-level explanation for the pH-dependent binding of beta-amyloid peptides by heparin and heparan sulfate proteoglycan and may have implications for amyloid formation.

Down syndrome candidate region-1 protein interacts with Tollip and positively modulates interleukin-1 receptor-mediated signaling

BACKGROUND

The Down syndrome candidate region-1 gene (DSCR1, also known as RCAN1) is situated close to the Down Syndrome Critical Region (DSCR), which contains genes responsible for many features of Down syndrome. DSCR1 modulates calcineurin phosphatase activity, though its functional role is incompletely understood.

METHODS

Here we investigated the role of DSCR1-1S isoform in IL-1 receptor (IL-1R)-mediated signaling by analyzing interaction between DSCR1-1S and the IL-1R pathway components Tollip, IRAK-1, and TRAF6.

RESULTS

Co-immunoprecipitation analyses of HEK293 cells revealed that DSCR1-1S interacted with Tollip, an IRAK-1 inhibitor, leading to the dissociation of IRAK-1 from Tollip. Similarly, both DSCR1-1S and Tollip interacted with TRAF6, with DSCR1 reducing interaction between Tollip and TRAF6. DSCR1-1S also stimulated IL-1R-mediated signaling pathways, TAK1 activation, NF-kappaB transactivation, and IL-8 production, all downstream consequences of IL-1R activation.

GENERAL SIGNIFICANCE

Together, these results suggest that DSCR1-1S isoform positively modulates IL-1R-mediated signaling pathways by regulating Tollip/IRAK-1/TRAF6 complex formation.

Multiple core homeodomain binding motifs differentially contribute to transcriptional activity of the murine gonadotropin-releasing hormone receptor gene promoter

Multiple homeodomain (Hbox) proteins have been shown to organize expression of key markers of gonadotropes. Nine putative Hbox-binding sites, characterized by the homeospecific TAAT motif, are located within the proximal 600 bp of the murine GnRHR promoter. Homeoproteins bind separate Hbox sites within this promoter, supporting basal- and endocrine-directed transcription. The function of the most proximal sites (Hbox1 and Hbox2) in the murine GnRHR is unknown; thus, understanding of the global contribution of homeospecific TAAT sites to promoter function is incomplete. Site-directed mutagenesis revealed that loss of Hbox2 reduced promoter activity in a cell-specific manner, having no effect in alphaT3-1 cells but reducing promoter function in LbetaT2 cells, another gonadotrope-derived cell line representing a later developmental stage. In contrast, eliminating Hbox1 reduced basal activity in both lines. This region displayed specific binding to homeoprotein Oct-1. Mutagenesis of a previously identified Oct-1-binding site in concert with Hbox1 led to further reduction in activity. We suggest that the two most proximal homeodomain-binding sites in the murine GnRHR promoter may regulate the promoter in a developmentally dependent fashion and that Oct-1 acts at multiple but distinct TAAT sites to support basal transcription.

A novel androgen receptor amino terminal region reveals two classes of amino/carboxyl interaction-deficient variants with divergent capacity to activate responsive sites in chromatin

The androgen receptor (AR) is an important signaling molecule in multiple tissues, yet its mode of action and cell-specific activities remain enigmatic. AR function has been best studied in the prostate, in which it is essential for growth and homeostasis of the normal organ as well as each stage of cancer development. Investigation of mechanisms responsible for continued AR action that evolve during prostate cancer progression or after hormonal management of the disease have been instructive in defining AR signaling pathways. In the current paper, we use sequence similarity and the collocation of somatic mutations in prostate cancer to define residues 501-535 of the AR amino-terminal domain as an important mediator of receptor function. Specifically, the 501-535 region is required for optimal interaction of the amino-terminal domain with both the p160 coactivator, nuclear receptor coactivator-2, and the AR-ligand binding domain in the amino/carboxyl (N/C) interaction. The N/C interaction is decreased by deletion of the 501-535 region but is distinct from deletion of the (23)FQNLF(27) peptide in that it does not affect the capacity of the AR to activate transcription from a chromatin integrated reporter or recruitment of the receptor to androgen-responsive loci in vivo. Collectively, we have been able to outline two classes of N/C-deficient AR variant that are divergent in their capacity to act in a chromatin context, thereby further defining the interplay between N/C interaction and coregulator recruitment via multiple receptor domains. These mechanisms are likely to be key determinants of the cell and promoter specific activities of the AR.

Gating of the CFTR Cl- channel by ATP-driven nucleotide-binding domain dimerisation

The cystic fibrosis transmembrane conductance regulator (CFTR) plays a fundamental role in fluid and electrolyte transport across epithelial tissues. Based on its structure, function and regulation, CFTR is an ATP-binding cassette (ABC) transporter. These transporters are assembled from two membrane-spanning domains (MSDs) and two nucleotide-binding domains (NBDs). In the vast majority of ABC transporters, the NBDs form a common engine that utilises the energy of ATP hydrolysis to pump a wide spectrum of substrates through diverse transmembrane pathways formed by the MSDs. By contrast, in CFTR the MSDs form a pathway for passive anion flow that is gated by cycles of ATP binding and hydrolysis by the NBDs. Here, we consider how the interaction of ATP with two ATP-binding sites, formed by the NBDs, powers conformational changes in CFTR structure to gate the channel pore. We explore how conserved sequences from both NBDs form ATP-binding sites at the interface of an NBD dimer and highlight the distinct roles that each binding site plays during the gating cycle. Knowledge of how ATP gates the CFTR Cl- channel is critical for understanding CFTR's physiological role, its malfunction in disease and the mechanism of action of small molecules that modulate CFTR channel gating.

Talin phosphorylation by Cdk5 regulates Smurf1-mediated talin head ubiquitylation and cell migration

Cell migration is a dynamic process that requires temporal and spatial regulation of integrin activation and focal adhesion assembly/disassembly. Talin, an actin and beta-integrin tail-binding protein, is essential for integrin activation and focal adhesion formation. Calpain-mediated cleavage of talin has a key role in focal adhesion turnover; however, the talin head domain, one of the two cleavage products, stimulates integrin activation, localizes to focal adhesions and maintains cell edge protrusions, suggesting that other steps, downstream of talin proteolysis, are required for focal adhesion disassembly. Here we show that talin head binds Smurf1, an E3 ubiquitin ligase involved in cell polarity and migration, more tightly than full-length talin does and that this interaction leads to talin head ubiquitylation and degradation. We found that talin head is a substrate for Cdk5, a cyclin-dependent protein kinase that is essential for cell migration, synaptic transmission and cancer metastasis. Cdk5 phosphorylated talin head at Ser 425, inhibiting its binding to Smurf1, thus preventing talin head ubiquitylation and degradation. Expression of the mutant tal(S425A), which resists Cdk5 phosphorylation thereby increasing its susceptibility to Smurf1-mediated ubiqitylation, resulted in extensive focal adhesion turnover and inhibited cell migration. Thus, talin head produced by calpain-induced cleavage of talin is degraded through Smurf1-mediated ubiquitylation; moreover, phosphorylation by Cdk5 regulates the binding of Smurf1 to talin head, controlling talin head turnover, adhesion stability and ultimately, cell migration.

ECM1 interacts with fibulin-3 and the beta 3 chain of laminin 332 through its serum albumin subdomain-like 2 domain

The extracellular matrix protein 1 (ECM1) is an 85 kDa secreted glycoprotein, comprising four variants and playing a pivotal role in endochondral bone formation, angiogenesis, and tumour biology. A computational model for the three-dimensional structure of ECM1a was determined to identify the potential and/or concealed region(s) for binding with candidate partners in human skin. Multiple alignments for the secondary structure of ECM1a and b revealed similarity with serum albumin. The N-terminal domain of ECM1a consists mainly of alpha-helices (alphaD1), while the remaining three domains, namely serum albumin subdomain-like (SASDL) domains 2-4, were topologically comparable with the subdomain of the third serum albumin domain. Yeast-two-hybrid screening of a human foreskin cDNA library using both full-length ECM1a and the hot spot region for ECM1 gene mutations in lipoid proteinosis, an autosomal recessive genodermatosis (complete SASDL2 and the linker to SASDL3: aa177-aa361), as bait, isolated seven extracellular proteins. The site-specific interaction of ECM1a with two of these candidate binders, laminin 332 beta-3 chain and fibulin-3, was confirmed by in vitro and in vivo co-immunoprecipitation experiments. Immunohistologically both binders co-localized with ECM1 in human skin. Together, ECM1 is a multifunctional binding core and/or a scaffolding protein interacting with a variety of extracellular and structural proteins, contributing to the maintenance of skin integrity and homeostasis. Hence, disruption of the ECM1 function may cause the failure of multi-communication among the surrounding skin interstitial molecules, as seen in lipoid proteinosis pathology.

[An infrared imaging system for detecting electrophoretic mobility shift of DNA-protein complexes]

OBJECTIVE

To establish a new non-radioactive method for electrophoretic mobility shift assay (EMSA) to investigate the binding between glucocorticoid induced leucine zipper (GILZ) and peroxisome proliferator-activated receptor-gamma 2 (PPARgamma2) promoter oligonucleotides.

METHODS

GILZ protein prepared by prokaryotic expression was linked to PPARgamma2 promoter oligonucleotides end-labeled with IRDye 800 infrared dye. The DNA-protein complex was separated with non-denatured polyacrylamide gel and scanned with the Odyssey. Infrared Imaging System.

RESULTS

One lane of DNA-protein complex was clearly presented, and the signal intensity increased along with the increment of the protein load.

CONCLUSION

This infrared imaging system can be used for EMSA for detecting the DNA-protein complex with high sensitivity efficiency and allows easy operation.

Negative cooperativity in H2 relaxin binding to a dimeric relaxin family peptide receptor 1

H2 relaxin, a member of the insulin superfamily, binds to the G-protein-coupled receptor RXFP1 (relaxin family peptide 1), a receptor that belongs to the leucine-rich repeat (LRR)-containing subgroup (LGRs) of class A GPCRs. We recently demonstrated negative cooperativity in INSL3 binding to RXFP2 and showed that this subgroup of GPCRs functions as constitutive dimers. In this work, we investigated whether the binding of H2 relaxin to RXFP1 also shows negative cooperativity, and whether this receptor functions as a dimer using BRET(2). Both binding and dissociation were temperature dependent, and the pH optimum for binding was pH 7.0. Our results showed that RXFP1 is a constitutive dimer with negative cooperativity in ligand binding, that dimerization occurs through the 7TM domain, and that the ectodomain has a stabilizing effect on this interaction. Dimerization and negative cooperativity appear to be general properties of LGRs involved in reproduction as well as other GPCRs.

The major histocompatibility complex class Ib molecule HLA-E at the interface between innate and adaptive immunity

The non-classical major histocompatibility complex (MHC) class I molecule human leucocyte antigen (HLA)-E is the least polymorphic of all the MHC class I molecules and acts as a ligand for receptors of both the innate and the adaptive immune systems. The recognition of self-peptides complexed to HLA-E by the CD94-NKG2A receptor expressed by natural killer (NK) cells represents a crucial checkpoint for immune surveillance by NK cells. However, HLA-E can also be recognised by the T-cell receptor expressed by alphabeta CD8 T cells and therefore can play a role in the adaptive immune response to invading pathogens. The recent resolution of HLA-E in complex with both innate and adaptive ligands has provided insight into the dual role of this molecule in immunity.

Organization of physical interactomes as uncovered by network schemas

Large-scale protein-protein interaction networks provide new opportunities for understanding cellular organization and functioning. We introduce network schemas to elucidate shared mechanisms within interactomes. Network schemas specify descriptions of proteins and the topology of interactions among them. We develop algorithms for systematically uncovering recurring, over-represented schemas in physical interaction networks. We apply our methods to the S. cerevisiae interactome, focusing on schemas consisting of proteins described via sequence motifs and molecular function annotations and interacting with one another in one of four basic network topologies. We identify hundreds of recurring and over-represented network schemas of various complexity, and demonstrate via graph-theoretic representations how more complex schemas are organized in terms of their lower-order constituents. The uncovered schemas span a wide range of cellular activities, with many signaling and transport related higher-order schemas. We establish the functional importance of the schemas by showing that they correspond to functionally cohesive sets of proteins, are enriched in the frequency with which they have instances in the H. sapiens interactome, and are useful for predicting protein function. Our findings suggest that network schemas are a powerful paradigm for organizing, interrogating, and annotating cellular networks.

Intramolecular cohesion of coils mediated by phenylalanine--glycine motifs in the natively unfolded domain of a nucleoporin

The nuclear pore complex (NPC) provides the sole aqueous conduit for macromolecular exchange between the nucleus and the cytoplasm of cells. Its diffusion conduit contains a size-selective gate formed by a family of NPC proteins that feature large, natively unfolded domains with phenylalanine–glycine repeats (FG domains). These domains of nucleoporins play key roles in establishing the NPC permeability barrier, but little is known about their dynamic structure. Here we used molecular modeling and biophysical techniques to characterize the dynamic ensemble of structures of a representative FG domain from the yeast nucleoporin Nup116. The results showed that its FG motifs function as intramolecular cohesion elements that impart order to the FG domain and compact its ensemble of structures into native premolten globular configurations. At the NPC, the FG motifs of nucleoporins may exert this cohesive effect intermolecularly as well as intramolecularly to form a malleable yet cohesive quaternary structure composed of highly flexible polypeptide chains. Dynamic shifts in the equilibrium or competition between intra- and intermolecular FG motif interactions could facilitate the rapid and reversible structural transitions at the NPC conduit needed to accommodate passing karyopherin–cargo complexes of various shapes and sizes while simultaneously maintaining a size-selective gate against protein diffusion.

The nuclear pore complex is a molecular filter that gates macromolecular exchange between the cytoplasm and the nucleoplasm of cells. It contains a size-selective diffusion barrier at its center composed of proteins named FG nucleoporins. These nucleoporins feature large, structurally disordered domains that are highly decorated with phenylalanine–glycine (FG) sequence motifs. The dynamic structure of these disordered FG domains excludes them from classical structural biology analyses such as X-ray crystallography; thus, new approaches are needed to characterize their shape. Here computational and biophysical approaches were used to elucidate the ensemble of structures adopted by the FG domain of a nucleoporin. The analyses showed that the FG motifs function as intramolecular cohesion elements that compact the shape of the FG domain, forcing it to adopt loosely knit globular configurations that are constantly reconfiguring. Within the nuclear pore complex, dozens of these nucleoporin FG domains may stack as loosely knit globules forming a porous sieve that gates molecular diffusion by size exclusion.

Classifying RNA-binding proteins based on electrostatic properties

Protein structure can provide new insight into the biological function of a protein and can enable the design of better experiments to learn its biological roles. Moreover, deciphering the interactions of a protein with other molecules can contribute to the understanding of the protein's function within cellular processes. In this study, we apply a machine learning approach for classifying RNA-binding proteins based on their three-dimensional structures. The method is based on characterizing unique properties of electrostatic patches on the protein surface. Using an ensemble of general protein features and specific properties extracted from the electrostatic patches, we have trained a support vector machine (SVM) to distinguish RNA-binding proteins from other positively charged proteins that do not bind nucleic acids. Specifically, the method was applied on proteins possessing the RNA recognition motif (RRM) and successfully classified RNA-binding proteins from RRM domains involved in protein–protein interactions. Overall the method achieves 88% accuracy in classifying RNA-binding proteins, yet it cannot distinguish RNA from DNA binding proteins. Nevertheless, by applying a multiclass SVM approach we were able to classify the RNA-binding proteins based on their RNA targets, specifically, whether they bind a ribosomal RNA (rRNA), a transfer RNA (tRNA), or messenger RNA (mRNA). Finally, we present here an innovative approach that does not rely on sequence or structural homology and could be applied to identify novel RNA-binding proteins with unique folds and/or binding motifs.

Gene expression in all living organisms is regulated by a complex set of events at both transcriptional and posttranscriptional levels. RNA-binding proteins play a key role in posttranscriptional events including splicing, stability, transport, and translation. Nowadays, there is increasing evidence that many other cellular processes may be mediated by RNA. Identifying new proteins involved in interaction with RNA is thus essential to unraveling the cellular processes in which these interactions are involved. In the current study we present a successful computational approach for classifying RNA-binding proteins and distinguishing them from other proteins based on structural and electrostatic properties. We test the method on a unique protein domain, the RNA recognition motif (RRM), which mediates both RNA and protein interactions. We show that we can discriminate RNA-binding RRMs from protein-binding RRMs. Further, we demonstrate that we can classify known RNA-binding proteins based on their RNA target (mRNA, rRNA, or tRNA). Our method does not rely on any kind of evolutionary information and thus can be applied to identify RNA-binding proteins with novel modes of RNA recognition.

Identification of residues in the receptor-binding domain (RBD) of the spike protein of human coronavirus NL63 that are critical for the RBD-ACE2 receptor interaction

Human coronavirus NL63 (NL63), a member of the group I coronaviruses, may cause acute respiratory diseases in young children and immunocompromised adults. Like severe acute respiratory syndrome coronavirus (SARS-CoV), NL63 also employs the human angiotensin-converting enzyme 2 (hACE2) receptor for cellular entry. To identify residues in the spike protein of NL63 that are important for hACE2 binding, this study first generated a series of S1-truncated variants, examined their associations with the hACE2 receptor and subsequently mapped a minimal receptor-binding domain (RBD) that consisted of 141 residues (aa 476-616) towards the C terminus of the S1 domain. The data also demonstrated that the NL63 RBD bound to hACE2 more efficiently than its full-length counterpart and had a binding efficiency comparable to the S1 or RBD of SARS-CoV. A further series of RBD variants was generated using site-directed mutagenesis and random mutant library screening assays, and identified 15 residues (C497, Y498, V499, C500, K501, R518, R530, V531, G534, G537, D538, S540, E582, W585 and T591) that appeared to be critical for the RBD-hACE2 association. These critical residues clustered in three separate regions (designated RI, RII and RIII) inside the RBD, which may represent three receptor-binding sites. These results may help to delineate the molecular interactions between the S protein of NL63 and the hACE2 receptor, and may also enhance our understanding of the pathogenesis of NL63 and SARS-CoV.

Impact of heterodimerization on intracellular localization of the ecdysteroid receptor (EcR)

Initially, nuclear import of the ecdysteroid receptor (EcR) in vertebrate cells (CHO-K1 and COS-7) does not afford a heterodimerization partner. Later on, EcR is retained in the nucleus only in the presence of a heterodimerization partner. Ultraspiracle (Usp) is more efficient compared to its vertebrate orthologue RXR and leads to an exclusively nuclear localization of EcR even in the absence of ligand. The DNA binding domain of the heterodimerization partner is important for retainment of EcR in the nucleus as shown by Usp4 (Usp(R130C)), which has lost its DNA binding capability. The C-terminal end of Usp (Usp(Delta205-508)) encompassing the C-terminal part of the D-domain and the E- and F-domains are essential for retainment of EcR in the nucleus. Nuclear localization is further influenced by cell-specific factors, since hormone and heterodimerization stabilizes the EcR protein in a cell-specific way.

Diverse and conserved roles of CLE peptides

The function of plant CLAVATA3 (CLV3)/ENDOSPERM SURROUNDING REGION (ESR) (CLE) peptides in shoot meristem differentiation has been expanded in recent years to implicate roles in root growth and vascular development among different CLE family members. Recent evidence suggests that nematode pathogens within plant roots secrete ligand mimics of plant CLE peptides to modify selected host cells into multinucleate feeding sites. This discovery demonstrated an unprecedented adaptation of an animal gene product to functionally mimic a plant peptide involved in cellular signaling for parasitic benefit. This review highlights the diverse and conserved role of CLE peptides in these different contexts.

The interaction domains of the plant Myc-like bHLH transcription factors can regulate the transactivation strength

The N-terminal region of the plant Myc-like basic helix-loop-helix transcription factors (bHLH TFs) contains two domains. Approximately, 190 amino acids at the N-terminus comprise an interaction domain, a.k.a. Myb-interacting-region (MIR) for its primary function of interacting with Myb-like TFs. Following, the interaction domain is an activation (or acidic) domain responsible for transactivation. We have previously discovered that a lysine to methionine substitution (K157M) in the interaction domain of Myc-RP of Perilla frutescens leads to a 50-fold increase in transactivation activity. The result suggests that mutations in the interaction domain affect transactivation. The highly conserved nature of this lysine residue in many Myc-like bHLH TFs prompted us to explore the functional importance of this residue within the TF family and the influence of the interaction domain on the activation domain in transactivation. We found that the replacement of the equivalent lysine with methionine significantly affects the transactivation activities of two other Myc-RP homologues, Delila from snapdragon and Lc from maize. In addition to methionine, substitution with several other amino acids at this position has positive effects on transcriptional activity. A neighboring conserved alanine residue (A159 in Myc-RP, A161 in Delila and A172 in Lc) also affects transactivation. Substitution of this alanine residue to an aspartic acid abolished transactivation of both Myc-RP and Delila and severely reduced transactivation of Lc. Ectopic expression of a Myc-RP K157M mutant in transgenic tobacco resulted in increased anthocyanin accumulation compared to plants expressing the wild-type gene. Our study reveals the potential cooperation between functional domains of the bHLH TFs.

Smad3 and Pitx2 cooperate in stimulation of FSHbeta gene transcription

Activin is a member of the TGFbeta superfamily of growth and differentiation factors that control a variety of cellular and physiological functions. The canonical intracellular pathway of this ligand is well established and involves Smad signaling molecules. The tissue- and cell-specificity of activin action is achieved by Smad interaction with various transcriptional co-factors in the nucleus. In the reproductive axis, activin induces biosynthesis and secretion of follicle stimulating hormone (FSH) through transcriptional control of FSHbeta-subunit. Whereas it has been well demonstrated that this regulation is mediated by Smad pathway, the molecular mechanisms underlying gonadotrope-specific expression of the FSHbeta gene are not fully understood. Previously, we have identified Pitx2 as a pituitary-expressed transcription factor involved in activin-dependent induction of the FSHbeta promoter. Present data demonstrate that Pitx2 is not only sufficient, but also necessary for FSHbeta gene transcription, as a siRNA-mediated downregulation of Pitx2 protein expression abrogates both Smad3- and activin-mediated stimulation of the FSHbeta promoter. In addition, downregulation of Smad3 protein expression has a significant effect on Pitx2-dependent stimulation of the FSHbeta promoter, suggesting that cooperation between these factors is necessary for full transcriptional activation of the FSHbeta promoter. Furthermore, we show that Pitx2/Smad protein complexes assemble and can be co-immunoprecipitated. This interaction is mediated through the homeodomain of Pitx2 and is important for stimulation of FSHbeta gene transcription. Overall, these data contribute to the emerging molecular mechanism underlying both basal and activin-dependent FSHbeta gene regulation.

The division inhibitor EzrA contains a seven-residue patch required for maintaining the dynamic nature of the medial FtsZ ring

The essential cytoskeletal protein FtsZ assembles into a ring-like structure at the nascent division site and serves as a scaffold for the assembly of the prokaryotic division machinery. We previously characterized EzrA as an inhibitor of FtsZ assembly in Bacillus subtilis. EzrA interacts directly with FtsZ to prevent aberrant FtsZ assembly and cytokinesis at cell poles. EzrA also concentrates at the cytokinetic ring in an FtsZ-dependent manner, although its precise role at this position is not known. Here, we identified a conserved patch of amino acids in the EzrA C terminus that is essential for localization to the FtsZ ring. Mutations in this patch (designated the "QNR patch") abolish EzrA localization to midcell but do not significantly affect EzrA's ability to inhibit FtsZ assembly at cell poles. ezrA QNR patch mutant cells exhibit stabilized FtsZ assembly at midcell and are significantly longer than wild-type cells, despite lacking extra FtsZ rings. These results indicate that EzrA has two distinct activities in vivo: (i) preventing aberrant FtsZ ring formation at cell poles through inhibition of de novo FtsZ assembly and (ii) maintaining proper FtsZ assembly dynamics within the medial FtsZ ring, thereby rendering it sensitive to the factors responsible for coordinating cell growth and cell division.

Polyelectrostatic interactions of disordered ligands suggest a physical basis for ultrasensitivity

Regulation of biological processes often involves phosphorylation of intrinsically disordered protein regions, thereby modulating protein interactions. Initiation of DNA replication in yeast requires elimination of the cyclin-dependent kinase inhibitor Sic1 via the SCF(Cdc4) ubiquitin ligase. Intriguingly, the substrate adapter subunit Cdc4 binds to Sic1 only after phosphorylation of a minimum of any six of the nine cyclin-dependent kinase sites on Sic1. To investigate the physical basis of this ultrasensitive interaction, we consider a mean-field statistical mechanical model for the electrostatic interactions between a single receptor site and a conformationally disordered polyvalent ligand. The formulation treats phosphorylation sites as negative contributions to the total charge of the ligand and addresses its interplay with the strength of the favorable ligand-receptor contact. Our model predicts a threshold number of phosphorylation sites for receptor-ligand binding, suggesting that ultrasensitivity in the Sic1-Cdc4 system may be driven at least in part by cumulative electrostatic interactions. This hypothesis is supported by experimental affinities of Cdc4 for Sic1 fragments with different total charges. Thus, polyelectrostatic interactions may provide a simple yet powerful framework for understanding the modulation of protein interactions by multiple phosphorylation sites in disordered protein regions.

[XBP-1 interacts with estrogen receptor alpha (ERalpha)]

Estrogen receptor alpha (ERalpha) has been a primary target of treatment as well as a prognostic indicator for breast cancer. The level of human X-box binding protein 1 (XBP-1) mRNA was related with that of ERalpha in breast tumors and was over-expressed in some breast tumors. These previous studies suggested that XBP-1 may interact with ERalpha. XBP-1 has two isoforms, XBP-1S and XBP-1U, as the result of unique splicing. GST pull-down assay showed that both XBP-1S and XBP-1U bound to ERalpha in vitro. The binding of XBP-1S to ERalpha was stronger than that of XBP-1U to ERalpha. Co-immunoprecipitation revealed that the binding was in a ligand-independent manner. XBP-1S and XBP-1U interacted with the region of ERalpha that contains a DNA-binding domain. The ERalpha-interacting regions on XBP-1S and XBP-1U have been mapped to two regions, the N-terminal basic region leucine zipper domain (bzip) and the C-terminal activation domain. These findings suggest that XBP-1S and XBP-1U may participate in ERalpha signaling pathway through the mediation of ERalpha.