Resilience of protein-protein interaction networks as determined by their large-scale topological features

The relationship between the structure and function of biological networks constitutes a fundamental issue in systems biology. Particularly, the structure of protein-protein interaction networks is related to important biological functions. In this work, we investigated how such a resilience is determined by the large scale features of the respective networks. Four species are taken into account, namely yeast Saccharomyces cerevisiae, worm Caenorhabditis elegans, fly Drosophila melanogaster and Homo sapiens. We adopted two entropy-related measurements (degree entropy and dynamic entropy) in order to quantify the overall degree of robustness of these networks. We verified that while they exhibit similar structural variations under random node removal, they differ significantly when subjected to intentional attacks (hub removal). As a matter of fact, more complex species tended to exhibit more robust networks. More specifically, we quantified how six important measurements of the networks topology (namely clustering coefficient, average degree of neighbors, average shortest path length, diameter, assortativity coefficient, and slope of the power law degree distribution) correlated with the two entropy measurements. Our results revealed that the fraction of hubs and the average neighbor degree contribute significantly for the resilience of networks. In addition, the topological analysis of the removed hubs indicated that the presence of alternative paths between the proteins connected to hubs tend to reinforce resilience. The performed analysis helps to understand how resilience is underlain in networks and can be applied to the development of protein network models.

High functional coherence in k-partite protein cliques of protein interaction networks

We introduce a new topological concept called k-partite protein cliques to study protein interaction (PPI) networks. In particular, we examine functional coherence of proteins in k-partite protein cliques. A k-partite protein clique is a k-partite maximal clique comprising two or more nonoverlapping protein subsets between any two of which full interactions are exhibited. In the detection of PPI’s k-partite maximal cliques, we propose to transform PPI networks into induced K-partite graphs with proteins as vertices where edges only exist among the graph’s partites. Then, we present a k-partite maximal clique mining (MaCMik) algorithm to enumerate k-partite maximal cliques from K-partite graphs. Our MaCMik algorithm is applied to a yeast PPI network. We observe that there does exist interesting and unusually high functional coherence in k-partite protein cliques—most proteins in k-partite protein cliques, especially those in the same partites, share the same functions. Therefore, the idea of k-partite protein cliques suggests a novel approach to characterizing PPI networks, and may help function prediction for unknown proteins.

Topological aspects of the human protein interaction network

It is currently widely accepted that the understanding of complex cell functions depends on an integrated network theoretical approach and not on an isolated view of the different molecular agents. Aim of this thesis was the examination of topological properties that mirror known biological aspects by depicting the human protein network with methods from graph- and network theory. The presented network is a partial human interactome of 9222 proteins and 36324 interactions, consisting of single interactions reliably extracted from peer-reviewed scientific publications. In general, one can focus on intra- or intermodular characteristics, where a functional module is defined as "a discrete entity whose function is separable from those of other modules". It is found that the presented human network is also scale-free and hierarchically organised, as shown for yeast networks before. The interactome also exhibits proteins with high betweenness and low connectivity which are biologically analyzed and interpreted here as shuttling proteins between organelles (e.g. ER to Golgi, internal ER protein translocation, peroxisomal import, nuclear pores import/export) for the first time. As an optimisation for finding proteins that connect modules, a new method is developed here based on proteins located between highly clustered regions, rather than regarding highly connected regions. As a proof of principle, the Mediator complex is found in first place, the prime example for a connector complex. Focusing on intramodular aspects, the measurement of k-clique communities discriminates overlapping modules very well. Twenty of the largest identified modules are analysed in detail and annotated to known biological structures (e.g. proteasome, the NFκB-, TGF-β complex). Additionally, two large and highly interconnected modules for signal transducer and transcription factor proteins are revealed, separated by known shuttling proteins. These proteins yield also the highest number of redundant shortcuts (by calculating the skeleton), exhibit the highest numbers of interactions and might constitute highly interconnected but spatially separated rich-clubs either for signal transduction or for transcription factors. This design principle allows manifold regulatory events for signal transduction and enables a high diversity of transcription events in the nucleus by a limited set of proteins. Altogether, biological aspects are mirrored by pure topological features, leading to a new view and to new methods that assist the annotation of proteins to biological functions, structures and subcellular localisations. As the human protein network is one of the most complex networks at all, these results will be fruitful for other fields of network theory and will help understanding complex network functions in general.

NETWORK TOPOLOGY IN HUMAN PROTEIN INTERACTION DATA PREDICTS FUNCTIONAL ASSOCIATION

High-throughput assays, such as yeast two-hybrid system, have generated a huge amount of protein-protein interaction (PPI) data in the past decade. This tremendously increases the need for developing reliable methods to systematically and automatically suggest protein functions and relationships between them. With the available PPI data, it is now possible to study the functions and relationships in the context of a large-scale network. To data, several network-based schemes have been provided to effectively annotate protein functions on a large scale. However, due to those inherent noises in high-throughput data generation, new methods and algorithms should be developed to increase the reliability of functional annotations. Previous work in a yeast PPI network (Samanta and Liang, 2003) has shown that the local connection topology, particularly for two proteins sharing an unusually large number of neighbors, can predict functional associations between proteins, and hence suggest their functions. One advantage of the work is that their algorithm is not sensitive to noises (false positives) in high-throughput PPI data. In this study, we improved their prediction scheme by developing a new algorithm and new methods which we applied on a human PPI network to make a genome-wide functional inference. We used the new algorithm to measure and reduce the influence of hub proteins on detecting functionally associated proteins. We used the annotations of the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) as independent and unbiased benchmarks to evaluate our algorithms and methods within the human PPI network. We showed that, compared with the previous work from Samanta and Liang, our algorithm and methods developed in this study improved the overall quality of functional inferences for human proteins. By applying the algorithms to the human PPI network, we obtained 4,233 significant functional associations among 1,754 proteins. Further comparisons of their KEGG and GO annotations allowed us to assign 466 KEGG pathway annotations to 274 proteins and 123 GO annotations to 114 proteins with estimated false discovery rates of <21% for KEGG and <30% for GO. We clustered 1,729 proteins by their functional associations and made pathway analysis to identify several subclusters that are highly enriched in certain signaling pathways. Particularly, we performed a detailed analysis on a subcluster enriched in the transforming growth factor β signaling pathway (P<10-50) which is important in cell proliferation and tumorigenesis. Analysis of another four subclusters also suggested potential new players in six signaling pathways worthy of further experimental investigations. Our study gives clear insight into the common neighbor-based prediction scheme and provides a reliable method for large-scale functional annotations in this post-genomic era.

Protein-protein interaction between phototaxis receptor sensory Rhodopsin I and its transducer HtrI

The molecular complex containing the seven transmembrane helix photoreceptor S&barbelow;ensory R&barbelow;hodopsin I&barbelow; (SRI) and transducer protein HtrI (H&barbelow;alobacterial Transducer for SRI&barbelow;) mediates color-sensitive phototaxis responses in the archaeon Halobacterium salinarum. Orange light causes an attractant response by a one-photon reaction and white light (orange + UV light) a repellent response by a two-photon reaction. Three aspects of SRI-HtrI structure/function and the signal transduction pathway were explored. First, the coupling of HtrI to the photoactive site of SRI was analyzed by mutagenesis and kinetic spectroscopy. Second, SRI-HtrI mutations and suppressors were selected and characterized to elucidate the color-sensing mechanism. Third, the signal relay through the transducer-bound histidine kinase was analyzed using an in vitro reconstitution system with known and newly identified taxis components. ^ Twenty-one mutations on HtrI were introduced by site-directed mutagenesis. Several replacements of charged residues perturbed the photochemical kinetics of SRI which led to the finding of a cluster of residues at the membrane/cytoplasm interface in HtrI electrostatically coupled to the photoactive site of SRI. We found by laser-flash kinetic spectroscopy that the transducer and these residues have specific effects on the light-induced proton transfer between the retinal chromophore and the protein. ^ One of the mutations showed an unusual mutant phenotype we called “inverted” signaling, in which the cell produces a repellent response to normally attractant light. Therefore, this mutant (E56Q of HtrI) had lost the color-discrimination by the SRI-HtrI complex. We used suppressor analysis to better understand the phenotype. Certain suppressors resulted in return of attractant responses to orange light but with inversion of the normally repellent response to white light to an attractant response. To explain this and other results, we formulated the Conformational Shuttling model in which the HtrI-SRI complex is poised in a metastable equilibrium of two conformations shifted in opposite directions by orange and white light. We tested this model by behavioral analysis (computerized cell tracking and motion study) of double mutants of inverting and suppressing mutations and the results confirmed the equilibrium-shift explanation. ^ We developed an in vitro system for measuring the effect of purified transducer on the histidine-kinase CheAH that controls the flagellar motor switch. The rate of kinase autophosphorylation was stimulated >2 fold in the reconstitution of the complete signal transduction system from purified components from H. salinarum. The in vitro assay also showed that the kinase activity was reduced in the absence and in the presence of high levels of linker protein CheWH. (Abstract shortened by UMI.) ^

Inhibitor of cytochrome c messenger RNA -protein interaction in stimulated rat skeletal muscle

The purpose of this work was to examine the possible mechanisms for the regulation of cytochrome c gene expression in response to increased contractile activity in rat skeletal muscle. The working hypothesis was that increased contractile activity enhances cytochrome c gene expression through a cis-element. A 110% increase in cytochrome c mRNA concentration was observed in tibialis anterior (TA) muscle after 9 days of chronic stimulation. Similar difference (120%) exists between soleus (SO) muscle of higher contractile activity and white vastus lateralis (WV) muscle of lower contractile activity. These results suggest that the endogenous cytochrome c gene expression is regulated by contractile activity. Cytochrome c-reporter genes were injected into skeletal muscles to identify the cis-element that is responsible for the regulation. Although the data was inconclusive, part of it suggested the importance of the 3$\sp\prime$-untranslated region (3$\sp\prime$-UTR) in mediating the response to increased contractile activity.^ RNA gel mobility shift (GMSA) and ultraviolet (UV) cross-linking assays revealed specific RNA-protein interaction in a 50-nucleotide region of the 3$\sp\prime$-UTR in unstimulated TA muscle. Computer analysis predicted a stem-loop structure of 17 nucleotides, which provides a structural basis for RNA-protein interaction. These 17 nucleotides are 100% conserved among rat, mouse and human cytochrome c genes and their 13 pseudogenes, suggesting a functional role for this region. The RNA-protein interaction was significantly less in highly active SO muscle than in inactive WV muscle and was dramatically decreased in stimulated TA muscle due to a protein inhibitor(s) associated with ribosome. It is possible that cytochrome c mRNAs undergoing translation are subject to a compartmentalized regulatory influence.^ The conclusion from these results is that increases in contractile activity induce or activate a protein inhibitor(s) associated with ribosome in rat skeletal muscle. The inhibitor decreases RNA-protein interaction in the 3$\sp\prime$-UTR of cytochrome c mRNA, which may result in increased mRNA stability and/or translation. ^

Evidence for a voltage-dependent enhancement of neurotransmitter release mediated via the synaptic protein interaction site of N-type Ca2+ channels

Secretion of neurotransmitters is initiated by voltage-gated calcium influx through presynaptic, voltage-gated N-type calcium channels. These channels interact with the SNARE proteins, which are core components of the exocytosis process, via the synaptic protein interaction (synprint) site in the intracellular loop connecting domains II and III of their α1B subunit. Interruption of this interaction by competing synprint peptides inhibits fast, synchronous transmitter release. Here we identify a voltage-dependent, but calcium-independent, enhancement of transmitter release that is elicited by trains of action potentials in the presence of a hyperosmotic extracellular concentration of sucrose. This enhancement of transmitter release requires interaction of SNARE proteins with the synprint site. Our results provide evidence for a voltage-dependent signal that is transmitted by protein–protein interactions from the N-type calcium channel to the SNARE proteins and enhances neurotransmitter release by altering SNARE protein function.

A transcriptional activating region with two contrasting modes of protein interaction

A C-terminal segment of the yeast activator Gal4 manifests two functions: When tethered to DNA, it elicits gene activation, and it binds the inhibitor Gal80. Here we examine the effects on these two functions of cysteine and proline substitutions. We find that, although certain cysteine substitutions diminish interaction with Gal80, those substitutions have little effect on the activating function in vivo and interaction with TATA box-binding protein (TBP) in vitro. Proline substitutions introduced near residues critical for Gal80 binding abolish that interaction but once again have no effect on the activating function. Crosslinking experiments show that a defined position in the activating peptide is in close proximity to TBP and Gal80 in the two separate reactions and show that binding of the inhibitor blocks binding to TBP. Thus, the same stretch of amino acids are involved in two quite different protein–protein interactions: binding to Gal80, which depends on a precise sequence and the formation of a defined secondary structure, or interactions with the transcriptional machinery in vivo, which are not impaired by perturbations of either sequence or structure.

Genetic characterization of a mammalian protein-protein interaction domain by using a yeast reverse two-hybrid system.

Many biological processes rely upon protein-protein interactions. Hence, detailed analysis of these interactions is critical for their understanding. Due to the complexities involved, genetic approaches are often needed. In yeast and phage, genetic characterizations of protein complexes are possible. However, in multicellular organisms, such characterizations are limited by the lack of powerful selection systems. Herein we describe genetic selections that allow single amino acid changes that disrupt protein-protein interactions to be selected from large libraries of randomly generated mutant alleles. The strategy, based on a yeast reverse two-hybrid system, involves a first-step negative selection for mutations that affect interaction, followed by a second-step positive selection for a subset of these mutations that maintain expression of full-length protein (two-step selection). We have selected such mutations in the transcription factor E2F1 that affect its ability to heterodimerize with DP1. The mutations obtained identified a putative helix in the marked box, a region conserved among E2F family members, as an important determinant for interaction. This two-step selection procedure can be used to characterize any interaction domain that can be tested in the two-hybrid system.

Amino-terminal protein-protein interaction motif (POZ-domain) is responsible for activities of the promyelocytic leukemia zinc finger-retinoic acid receptor-alpha fusion protein.

Promyelocytic leukemia zinc finger-retinoic acid receptor a (PLZF-RARalpha), a fusion receptor generated as a result of a variant t(11;17) chromosomal translocation that occurs in a small subset of acute promyelocytic leukemia (APL) patients, has been shown to display a dominant-negative effect against the wild-type RARalpha/retinoid X receptor alpha (RXRalpha). We now show that its N-terminal region (called the POZ-domain), which mediates protein-protein interaction as well as specific nuclear localization of the wild-type PLZF and chimeric PLZF-RARalpha proteins, is primarily responsible for this activity. To further investigate the mechanisms of PLZF-RARalpha action, we have also studied its ligand-receptor, protein-protein, and protein-DNA interaction properties and compared them with those of the promyelocytic leukemia gene (PML)-RARalpha, which is expressed in the majority of APLs as a result of t(15;17) translocation. PLZF-RARalpha and PML-RARalpha have essentially the same ligand-binding affinities and can bind in vitro to retinoic acid response elements (RAREs) as homodimers or heterodimers with RXRalpha. PLZF-RARalpha homodimerization and heterodimerization with RXRalpha were primarily mediated by the POZ-domain and RARalpha sequence, respectively. Despite having identical RARalpha sequences, PLZF-RARalpha and PML-RARalpha homodimers recognized with different affinities distinct RAREs. Furthermore, PLZF-RARalpha could heterodimerize in vitro with the wild-type PLZF, suggesting that it may play a role in leukemogenesis by antagonizing actions of not only the retinoid receptors but also the wild-type PLZF and possibly other POZ-domain-containing regulators. These different protein-protein interactions and the target gene specificities of PLZF-RARalpha and PML-RARalpha may underlie, at least in part, the apparent resistance of APL with t(11;17) to differentiation effects of all-trans-retinoic acid.

Protein-protein interaction: a genetic selection for compensating mutations at the barnase-barstar interface.

Barnase and barstar are trivial names of the extracellular RNase and its intracellular inhibitor produced by Bacillus amyloliquefaciens. Inhibition involves the formation of a very tight one-to-one complex of the two proteins. With the crystallographic solution of the structure of the barnase-barstar complex and the development of methods for measuring the free energy of binding, the pair can be used to study protein-protein recognition in detail. In this report, we describe the isolation of suppressor mutations in barstar that compensate for the loss in interaction energy caused by a mutation in barnase. Our suppressor search is based on in vivo selection for barstar variants that are able to protect host cells against the RNAse activity of those barnase mutants not properly inhibited by wild-type barstar. This approach utilizes a plasmid system in which barnase expression is tightly controlled to keep the mutant barnase gene silent. When expression of barnase is turned on, failure to form a complex between the mutant barnase and barstar has a lethal effect on host cells unless overcome by substitution of the wild-type barstar by a functional suppressor derivative. A set of barstar suppressors has been identified for barnase mutants with substitutions in two amino acid positions (residues 102 and 59), which are critically involved in both RNase activity and barstar binding. The mutations selected as suppressors could not have been predicted on the basis of the known protein structures. The single barstar mutation with the highest information content for inhibition of barnase (H102K) has the substitution Y30W. The reduction in binding caused by the R59E mutation in barnase can be partly reversed by changing Glu-76 of barstar, which forms a salt bridge with the Arg-59 in the wild-type complex, to arginine, thus completing an interchange of the two charges.

A role for the Msx-1 homeodomain in transcriptional regulation: residues in the N-terminal arm mediate TATA binding protein interaction and transcriptional repression.

In a previous study we showed that the murine homeodomain protein Msx-1 is a potent transcriptional repressor and that this activity is independent of its DNA binding function. The implication of these findings is that repression by Msx-1 is mediated through its association with certain protein factors rather than through its interaction with DNA recognition sites, which prompted investigation of the relevant protein factors. Here we show that Msx-1 interacts directly with the TATA binding protein (TBP) but not with several other general transcription factors. This interaction is mediated by the Msx-1 homeodomain, specifically through residues in the N-terminal arm. These same N-terminal arm residues are required for repression by Msx-1, suggesting a functional relationship between TBP association and transcriptional repression. This is further supported by the observation that addition of excess TBP blocks the repressor action of Msx-1 in in vitro transcription assays. Finally, DNA binding activity is separable from both TBP interaction and repression, which further shows that these other activities of the Msx-1 homeodomain are distinct. Therefore, these findings define a role for the Msx-1 homeodomain, particularly the N-terminal arm residues in protein-protein interaction and transcriptional repression, and implicate a more complex role overall for homeodomains in transcriptional regulation.

Direct protein interaction underlies gene-for-gene specificity and coevolution of the flax resistance genes and flax rust avirulence genes

Plant resistance proteins (R proteins) recognize corresponding pathogen avirulence (Avr) proteins either indirectly through detection of changes in their host protein targets or through direct R-Avr protein interaction. Although indirect recognition imposes selection against Avr effector function, pathogen effector molecules recognized through direct interaction may overcome resistance through sequence diversification rather than loss of function. Here we show that the flax rust fungus AvrLS67 genes, whose products are recognized by the L5, L6, and L7 R proteins of flax, are highly diverse, with 12 sequence variants identified from six rust strains. Seven AvrL567 variants derived from Avr alleles induce necrotic responses when expressed in flax plants containing corresponding resistance genes (R genes), whereas five variants from avr alleles do not. Differences in recognition specificity between AvA567 variants and evidence for diversifying selection acting on these genes suggest they have been involved in a gene-specific arms race with the corresponding flax R genes. Yeast two-hybrid assays indicate that recognition is based on direct R-Avr protein interaction and recapitulate the interaction specificity observed in planta. Biochemical analysis of Escherichia coli-produced AvrL567 proteins shows that variants that escape recognition nevertheless maintain a conserved structure and stability, suggesting that the amino acid sequence differences directly affect the R-Avr protein interaction. We suggest that direct recognition associated with high genetic diversity at corresponding R and Avr gene loci represents an alternative outcome of plant-pathogen coevolution to indirect recognition associated with simple balanced polymorphisms for functional and nonfunctional R and Avr genes.

Tear protein interaction with hydrogel contact lenses

The design and synthesis of biomaterials covers a growing number of biomedical applications. The use of biomaterials in biological environment is associated with a number of problems, the most important of which is biocompatabUity. If the implanted biomaterial is not compatible with the environment, it will be rejected by the biological site. This may be manifested in many ways depending on the environment in which it is used. Adsorption of proteins takes place almost instantaneously when a biomaterial comes into contact with most biological fluids. The eye is a unique body site for the study of protein interactions with biomaterials, because of its ease of access and deceptive complexity of the tears. The use of contact lenses for either vision correction and cosmetic reasons or as a route for the controlled drug delivery, has significantly increased in recent years. It is relatively easy to introduce a contact lens Into the tear fluid and remove after a few minutes without surgery or trauma to the patient. A range of analytical techniques were used and developed to measure the proteins absorbed to some existing commercial contact lens materials and also to novel hydrogels synthesised within the research group. Analysis of the identity and quantity of proteins absorbed to biomaterials revealed the importance of many factors on the absorption process. The effect of biomaterial structure, protein nature in terms of size. shape and charge and pH of the environment on the absorption process were examined in order to determine the relative up-take of tear proteins. This study showed that both lysozyme and lactoferrin penetrate the lens matrix of ionic materials. Measurement of the mobility and activity of the protein deposited into the surface and within the matrix of ionic lens materials demonstrated that the mobility is pH dependent and, within the experimental errors, the biological activity of lysozyme remained unchanged after adsorption and desorption. The study on the effect of different monomers copolymerised with hydroxyethyl methacrylate (HEMA) on the protein up-take showed that monomers producing a positive charge on the copolymer can reduce the spoilation with lysozyme. The studies were extended to real cases in order to compare the patient dependent factors. The in-vivo studies showed that the spoilation is patient dependent as well as other factors. Studies on the extrinsic factors such as dye used in colour lenses showed that the addition of colourant affects protein absorption and, in one case, its effect is beneficial to the wearer as it reduces the quantity of the protein absorbed.