Explore the hidden treasure in protein-protein interaction networks - an iterative model for predicting protein functions

Protein-protein interaction networks constructed by high throughput technologies provide opportunities for predicting protein functions. A lot of approaches and algorithms have been applied on PPI networks to predict functions of unannotated proteins over recent decades. However, most of existing algorithms and approaches do not consider unannotated proteins and their corresponding interactions in the prediction process. On the other hand, algorithms which make use of unannotated proteins have limited prediction performance. Moreover, current algorithms are usually one-off predictions. In this paper, we propose an iterative approach that utilizes unannotated proteins and their interactions in prediction. We conducted experiments to evaluate the performance and robustness of the proposed iterative approach. The iterative approach maximally improved the prediction performance by 50%-80% when there was a high proportion of unannotated neighborhood protein in the network. The iterative approach also showed robustness in various types of protein interaction network. Importantly, our iterative approach initially proposes an idea that iteratively incorporates the interaction information of unannotated proteins into the protein function prediction and can be applied on existing prediction algorithms to improve prediction performance.

Protein function prediction based on protein – protein interaction networks.

This research has expanded the knowledge in Bioinformatics and Data mining. It makes an influential contribution to the future research in this area.

Interaction of metal salts with cytoskeletal motor protein systems

Interactions of chemicals with the microtubular network of cells may lead to genotoxicity. Micronuclei (MN) might be caused by interaction of metals with tubulin and/or kinesin. The genotoxic effects of inorganic lead and mercury salts were studied using the MN assay and the CREST analysis in V79 Chinese hamster fibroblasts. Effects on the functional activity of motor protein systems were examined by measurement of tubulin assembly and kinesin-driven motility. Lead and mercury salts induced MN dose-dependently. The no-effect-concentration for MN induction was 1.1 μM PbCl2, 0.05 μM Pb(OAc)2 and 0.01 μM HgCl2. The in vitro results obtained for PbCl2 correspond to reported MN induction in workers occupationally exposed to lead, starting at 1.2 μM Hg(II) (Vaglenov et al., 2001, Environ. Health Perspect. 109, 295-298). The CREST Analysis indicate aneugenic effects of Pb(II) and aneugenic and additionally clastogenic effects of Hg(II). Lead (chloride, acetate, and nitrate) and mercury (chloride and nitrate) interfered dose-dependently with tubulin assembly in vitro. The no-effect-concentration for lead salts in this assay was 10 μM. Inhibition of tubulin assembly by mercury started at 2 μM. The gliding velocity of microtubules along immobilised kinesin molecules was affected by 25 μM Pb(NO3)2 and 0.1 μM HgCl2 in a dose-dependent manner. Our data support the hypothesis that lead and mercury genotoxicity may result, at least in part, via disturbance of chromosome segregation via interaction with cytoskeletal proteins.

Nonspecific Interaction between DNA and Protein allows for Cooperativity: A Case Study with Mycobacterium DNA Binding Protein

Different DNA-binding proteins have different interaction modes with DNA. Sequence-specific DNA protein interaction has been mostly associated with regulatory processes inside a cell, and as such extensive studies have been made. Adequate data is also available on nonspecific DNA protein interaction, as an intermediate to protein's search for its cognate partner. Multidomain nonspecific DNA protein interaction involving physical sequestering of DNA has often been implicated to regulate gene expression indirectly. However, data available on this type of interaction is limited. One such interaction is the binding of DNA with mycobacterium DNA binding proteins. We have used the Langmuir-Blodgett technique to evaluate for the first time the kinetics and thermodynamics of Mycobacterium smegmatis Dps 1 binding to DNA. By immobilizing one of the interacting partners, we have shown that, when a kinetic bottleneck is applied, the binding mechanism showed cooperative binding (n = 2.72) at lower temperatures, but the degree of cooperativity gradually reduces (n = 1.38) as the temperature was increased We have also compared the kinetics and thermodynamics of sequence-specific and nonspecific DNA protein interactions under the same set of conditions.

Interaction of the D-2short dopamine receptor with G proteins: analysis of receptor G protein selectivity

The human D-2short (D-2S) dopamine receptor has been expressed together with the G proteins Gi2 and Go in insect cells using the baculovirus system. Levels of receptor were determined using [H-3]spiperone binding. Levels of G protein heterotrimer were determined using quantitative Western blot and using [S-35]GTPgammaS saturation binding experiments. Levels of the receptor and G protein and the receptor/G protein ratio were similar in the two preparations. Stimulation of [S-35]GTPgammaS binding by a range of agonists occurred with higher relative efficacy and in some cases higher potency in the preparation expressing Go, indicating that interaction of the D-2S receptor is more efficient with this G protein. The effects of various G protein-selective agents on 10,11-dihydroxy-N-n-propylnorapomorphine ([H-3]NPA) binding were used to examine the receptor/G protein complex in the two preparations. Suramin inhibited [H-3]NPA binding with slightly higher potency in the Gi2 preparation, whereas GppNHp inhibited [H-3]NPA binding with greater potency (similar to6-fold) in the Go preparation. This may imply that the G protein is more readily activated in the D-2S/Go preparation. [H-3]Spiperone binding occurred with an increased B-max in the presence of suramin in the Go preparation but not in the Gi2 preparation, suggesting a higher affinity interaction between the free receptor and this G protein. It is concluded that the higher efficiency activation of Go by the D-2S receptor may be a function of higher affinity receptor/G protein interaction as well as a greater ability to activate the G protein. (C) 2003 Elsevier Science Inc. All rights reserved.

Stoichiometry and thermodynamics of the interaction between the C-terminus of human 90 kDa heat shock protein Hsp90 and the mitochondrial translocase of outer membrane Tom70

A large majority of the 1000-1500 proteins in the mitochondria are encoded by the nuclear genome, and therefore, they are translated in the cytosol in the form and contain signals to enable the import of proteins into the organelle. The TOM complex is the major translocase of the outer membrane responsible for preprotein translocation. It consists of a general import pore complex and two membrane import receptors, Tom20 and Tom70. Tom70 contains a characteristic TPR domain, which is a docking site for the Hsp70 and Hsp90 chaperones. These chaperones are involved in protecting cytosolic preproteins from aggregation and then in delivering them to the TOM complex. Although highly significant, many aspects of the interaction between Tom70 and Hsp90 are still uncertain. Thus, we used biophysical tools to study the interaction between the C-terminal domain of Hsp90 (C-Hsp90), which contains the EEVD motif that binds to TPR domains, and the cytosolic fragment of Tom70. The results indicate a stoichiometry of binding of one monomer of Tom70 per dimer of C-Hsp90 with a K(D) of 360 30 nM, and the stoichiometry and thermodynamic parameters obtained suggested that Tom70 presents a different mechanism of interaction with Hsp90 when compared with other TPR proteins investigated. (C) 2011 Elsevier Inc. All rights reserved.

Analysis of protein sequence and interaction data for candidate disease gene prediction

Linkage analysis is a successful procedure to associate diseases with specific genomic regions. These regions are often large, containing hundreds of genes, which make experimental methods employed to identify the disease gene arduous and expensive. We present two methods to prioritize candidates for further experimental study: Common Pathway Scanning (CPS) and Common Module Profiling (CMP). CPS is based on the assumption that common phenotypes are associated with dysfunction in proteins that participate in the same complex or pathway. CPS applies network data derived from protein–protein interaction (PPI) and pathway databases to identify relationships between genes. CMP identifies likely candidates using a domain-dependent sequence similarity approach, based on the hypothesis that disruption of genes of similar function will lead to the same phenotype. Both algorithms use two forms of input data: known disease genes or multiple disease loci. When using known disease genes as input, our combined methods have a sensitivity of 0.52 and a specificity of 0.97 and reduce the candidate list by 13-fold. Using multiple loci, our methods successfully identify disease genes for all benchmark diseases with a sensitivity of 0.84 and a specificity of 0.63. Our combined approach prioritizes good candidates and will accelerate the disease gene discovery process.

Physical and functional interaction between the dopamine transporter and the synaptic vesicle protein synaptogyrin-3.

Uptake through the dopamine transporter (DAT) represents the primary mechanism used to terminate dopaminergic transmission in brain. Although it is well known that dopamine (DA) taken up by the transporter is used to replenish synaptic vesicle stores for subsequent release, the molecular details of this mechanism are not completely understood. Here, we identified the synaptic vesicle protein synaptogyrin-3 as a DAT interacting protein using the split ubiquitin system. This interaction was confirmed through coimmunoprecipitation experiments using heterologous cell lines and mouse brain. DAT and synaptogyrin-3 colocalized at presynaptic terminals from mouse striatum. Using fluorescence resonance energy transfer microscopy, we show that both proteins interact in live neurons. Pull-down assays with GST (glutathione S-transferase) proteins revealed that the cytoplasmic N termini of both DAT and synaptogyrin-3 are sufficient for this interaction. Furthermore, the N terminus of DAT is capable of binding purified synaptic vesicles from brain tissue. Functional assays revealed that synaptogyrin-3 expression correlated with DAT activity in PC12 and MN9D cells, but not in the non-neuronal HEK-293 cells. These changes were not attributed to changes in transporter cell surface levels or to direct effect of the protein-protein interaction. Instead, the synaptogyrin-3 effect on DAT activity was abolished in the presence of the vesicular monoamine transporter-2 (VMAT2) inhibitor reserpine, suggesting a dependence on the vesicular DA storage system. Finally, we provide evidence for a biochemical complex involving DAT, synaptogyrin-3, and VMAT2. Collectively, our data identify a novel interaction between DAT and synaptogyrin-3 and suggest a physical and functional link between DAT and the vesicular DA system.

Unveiling Protein Functions through the Dynamics of the Interaction Network

Protein interaction networks have become a tool to study biological processes, either for predicting molecular functions or for designing proper new drugs to regulate the main biological interactions. Furthermore, such networks are known to be organized in sub-networks of proteins contributing to the same cellular function. However, the protein function prediction is not accurate and each protein has traditionally been assigned to only one function by the network formalism. By considering the network of the physical interactions between proteins of the yeast together with a manual and single functional classification scheme, we introduce a method able to reveal important information on protein function, at both micro- and macro-scale. In particular, the inspection of the properties of oscillatory dynamics on top of the protein interaction network leads to the identification of misclassification problems in protein function assignments, as well as to unveil correct identification of protein functions. We also demonstrate that our approach can give a network representation of the meta-organization of biological processes by unraveling the interactions between different functional classes

The dimer-dimer interaction surface of the replication terminator protein of Bacillus subtilis and termination of DNA replication.

The replication terminator protein (RTP) of Bacillus subtilis causes polar fork arrest at replication termini by sequence-specific interaction of two dimeric proteins with the terminus sequence. The crystal structure of the RTP protein has been solved, and the structure has already provide valuable clues regarding the structural basis of its function. However, it provides little information as to the surface of the protein involved in dimer-dimer interaction. Using site-directed mutagenesis, we have identified three sites on the protein that appear to mediate the dimer-dimer interaction. Crystallographic analysis of one of the mutant proteins (Y88F) showed that its structure is unaltered when compared to the wild-type protein. The locations of the three sites suggested a model for the dimer-dimer interaction that involves an association between two beta-ribbon motifs. This model is supported by a fourth mutation that was predicted to disrupt the interaction and was shown to do so. Biochemical analyses of these mutants provide compelling evidence that cooperative protein-protein interaction between two dimers of RTP is essential to impose polar blocks to the elongation of both DNA and RNA chains.

Modeling of protein signaling networks in clinical proteomics

Molecular interactions that underlie pathophysiological states are being elucidated using techniques that profile proteomicend points in cellular systems. Within the field of cancer research, protein interaction networks play pivotal roles in the establishment and maintenance of the hallmarks of malignancy, including cell division, invasion, and migration. Multiple complementary tools enable a multifaceted view of how signal protein pathway alterations contribute to pathophysiological states.One pivotal technique is signal pathway profiling of patient tissue specimens. This microanalysis technology provides a proteomic snapshot at one point in time of cells directly procured from the native context of a tumor micro environment. To study the adaptive patterns of signal pathway events over time, before and after experimental therapy, it is necessary to obtain biopsies from patients before, during, and after therapy. A complementary approach is the profiling of cultured cell lines with and without treatment. Cultured cell models provide the opportunity to study short-term signal changes occurring over minutes to hours. Through this type of system, the effects of particular pharmacological agents may be used to test the effects of signal pathway inhibition or activation on multiple endpoints within a pathway. The complexity of the data generated has necessitated the development of mathematical models for optimal interpretation of interrelated signaling pathways. In combination,clinical proteomic biopsy profiling, tissue culture proteomic profiling, and mathematical modeling synergistically enable a deeper understanding of how protein associations lead to disease states and present new insights into the design of therapeutic regimens.

Protein Molecular Surface Mapped at Different Geometrical Resolutions

Many areas of biochemistry and molecular biology, both fundamental and applications-orientated, require an accurate construction, representation and understanding of the protein molecular surface and its interaction with other, usually small, molecules. There are however many situations when the protein molecular surface gets in physical contact with larger objects, either biological, such as membranes, or artificial, such as nanoparticles. The contribution presents a methodology for describing and quantifying the molecular properties of proteins, by geometrical and physico-chemical mapping of the molecular surfaces, with several analytical relationships being proposed for molecular surface properties. The relevance of the molecular surface-derived properties has been demonstrated through the calculation of the statistical strength of the prediction of protein adsorption. It is expected that the extension of this methodology to other phenomena involving proteins near solid surfaces, in particular the protein interaction with nanoparticles, will result in important benefits in the understanding and design of protein-specific solid surfaces. © 2013 Nicolau et al.

Protein loaded mesoporous silica spheres as a controlled delivery platform

Background The adsorption of bovine serum albumin (BSA) onto mesoporous silica spheres (MPS) synthesized from silica colloids was studied employing real time in situ measurements. The stabilities of the BSA at different pH values, their isoelectric points and zeta potentials were determined in order to probe the interactions between the protein and the mesoporous silica. Results The pore size of MPS was designed for protein, and this, coupled with an in depth understanding of the physico-chemical characteristics of the protein and MPS has yielded a better binding capacity and delivery profile. The adsorption isotherm at pH 4.2 fitted the Langmuir model and displayed the highest adsorption capacity (71.43 mg mL-1 MPS). Furthermore, the delivery rates of BSA from the MPS under physiological conditions were shown to be dependent on the ionic strength of the buffer and protein loading concentration. Conclusion Economics and scale-up considerations of mesoporous material synthesized via destabilization of colloids by electrolyte indicate the scaleability and commercial viability of this technology as a delivery platform for biopharmaceutical applications.

LacO-LacI interaction in affinity adsorption of plasmid DNA

Current approaches for purifying plasmids from bacterial production systems exploit the physiochemical properties of nucleic acids in non-specific capture systems. In this study, an affinity system for plasmid DNA (pDNA) purification has been developed utilizing the interaction between the lac operon (lacO) sequence contained in the pDNA and a 64mer synthetic peptide representing the DNA-binding domain of the lac repressor protein, LacI. Two plasmids were evaluated, the native pUC19 and pUC19 with dual lacO3/lacOs operators (pUC19lacO3/lacOs), where the lacOs operator is perfectly symmetrical. The DNA-protein affinity interaction was evaluated by surface plasmon resonance using a Biacore system. The affinity capture of DNA in a chromatography system was evaluated using LacI peptide that had been immobilized to Streamline™ adsorbent. The KD-values for double stranded DNA (dsDNA) fragments containing lacO1 and lacO3 and lacOs and lacO3 were 5.7 ± 0.3 × 10 -11 M and 4.1 ± 0.2 × 10-11 M respectively, which compare favorably with literature reports of 5 × 10-10 - 1 × 10-9 M for native laCO1 and 1-1.2 × 10-10 M for lacO1 in a saline buffer. Densitometric analysis of the gel bands from the affinity chromatography run clearly showed a significant preference for capture of the supercoiled fraction from the feed pDNA sample. The results indicate the feasibility of the affinity approach for pDNA capture and purification using native protein-DNA interaction.

On the Specificity of Carbohydrate-Lectin Recognition The Interaction of a Lectin from Ricinus communis Beans with Simple Saccharides and Concanavalin A

A galactose-specific protein (RC1) isolated from Ricinus communis beans was found to give a precipitin reaction with concanavalin A. Its carbohydrate content amounted to 8–9% of the total protein and was found to be rich in mannose. The interaction of RC1 with galactose and lactose was measured in 0.05 M phosphate buffer containing 0.2 M NaCl (pH 6.8) by the method of conventional equilibrium dialysis. From the analysis of the binding data according to Scatchard method the association constant (Ka) at 5°C was calculated as 3.8 mM−1 and 1.2 mM−1 for lactose and galactose, respectively. In both cases the number of binding sites per molecule of RC1 with molecular weight of 120000 was found to be 2. From the temperature-dependent Ka values for the binding of lactose, the values of –5.7 kcal/mol and –4.3 cal × mol−1× K−1 were calculated for ΔH and ΔS, respectively. The addition of concanavalin A to RC1 or vice versa led to the formation of the insoluble complex RC1· ConA4 containing one molecule of RC1 and one molecule of tetrameric concanavalin A (ConA4) which could be dissociated upon addition of concanavalin A-specific sugars. The complex formation results in a time-dependent appearance of turbidity in the time range from 10s to 10 min. From the measurement of the time-dependent appearance and disappearance of the turbidity the formation (kf) and dissociation (kd) rate constants were calculated as 3 mM−1× s−1 and 0.07 ks−1 respectively. The ratio kf/kd (43μM −1), that corresponds to the association constant of complex RC1· ConA4, is higher than that of mannoside · ConA4 and thereby suggests that protein-protein interaction contributes significantly in stabilising glycoprotein · lectin complexes. The relevance of this finding to the understanding of the chemical specificities that are involved in a model cell-lectin interaction is discussed.