The interactome: predicting the protein-protein interactions in cells

Spastin and alsin protein interactome analyses begin to reveal key canonical pathways and suggest novel druggable targets

Developing effective and long-term treatment strategies for rare and complex neurodegenerative diseases is challenging. One of the major roadblocks is the extensive heterogeneity among patients. This hinders understanding the underlying disease-causing mechanisms and building solutions that have implications for a broad spectrum of patients. One potential solution is to develop personalized medicine approaches based on strategies that target the most prevalent cellular events that are perturbed in patients. Especially in patients with a known genetic mutation, it may be possible to understand how these mutations contribute to problems that lead to neurodegeneration. Protein-protein interaction analyses offer great advantages for revealing how proteins interact, which cellular events are primarily involved in these interactions, and how they become affected when key genes are mutated in patients. This line of investigation also suggests novel druggable targets for patients with different mutations. Here, we focus on alsin and spastin, two proteins that are identified as "causative" for amyotrophic lateral sclerosis and hereditary spastic paraplegia, respectively, when mutated. Our review analyzes the protein interactome for alsin and spastin, the canonical pathways that are primarily important for each protein domain, as well as compounds that are either Food and Drug Administration-approved or are in active clinical trials concerning the affected cellular pathways. This line of research begins to pave the way for personalized medicine approaches that are desperately needed for rare neurodegenerative diseases that are complex and heterogeneous.

Dual fluorescence reporter mice for Ccl3 transcription, translation, and intercellular communication

Chemokines guide immune cells during their response against pathogens and tumors. Various techniques exist to determine chemokine production, but none to identify cells that directly sense chemokines in vivo. We have generated CCL3-EASER (ErAse, SEnd, Receive) mice that simultaneously report for Ccl3 transcription and translation, allow identifying Ccl3-sensing cells, and permit inducible deletion of Ccl3-producing cells. We infected these mice with murine cytomegalovirus (mCMV), where Ccl3 and NK cells are critical defense mediators. We found that NK cells transcribed Ccl3 already in homeostasis, but Ccl3 translation required type I interferon signaling in infected organs during early infection. NK cells were both the principal Ccl3 producers and sensors of Ccl3, indicating auto/paracrine communication that amplified NK cell response, and this was essential for the early defense against mCMV. CCL3-EASER mice represent the prototype of a new class of dual fluorescence reporter mice for analyzing cellular communication via chemokines, which may be applied also to other chemokines and disease models.

Encoding the space of protein-protein binding interfaces by artificial intelligence

The physical interactions between proteins are largely determined by the structural properties at their binding interfaces. It was found that the binding interfaces in distinctive protein complexes are highly similar. The structural properties underlying different binding interfaces could be further captured by artificial intelligence. In order to test this hypothesis, we broke protein-protein binding interfaces into pairs of interacting fragments. We employed a generative model to encode these interface fragment pairs in a low-dimensional latent space. After training, new conformations of interface fragment pairs were generated. We found that, by only using a small number of interface fragment pairs that were generated by artificial intelligence, we were able to guide the assembly of protein complexes into their native conformations. These results demonstrate that the conformational space of fragment pairs at protein-protein binding interfaces is highly degenerate. Features in this degenerate space can be well characterized by artificial intelligence. In summary, our machine learning method will be potentially useful to search for and predict the conformations of unknown protein-protein interactions.

Woven bone formation and mineralization by rat mesenchymal stromal cells imply increased expression of the intermediate filament desmin

Background

Disordered and hypomineralized woven bone formation by dysfunctional mesenchymal stromal cells (MSCs) characterize delayed fracture healing and endocrine -metabolic bone disorders like fibrous dysplasia and Paget disease of bone. To shed light on molecular players in osteoblast differentiation, woven bone formation, and mineralization by MSCs we looked at the intermediate filament desmin (DES) during the skeletogenic commitment of rat bone marrow MSCs (rBMSCs), where its bone-related action remains elusive.

Results

Monolayer cultures of immunophenotypically- and morphologically - characterized, adult male rBMSCs showed co-localization of desmin (DES) with vimentin, F-actin, and runx2 in all cell morphotypes, each contributing to sparse and dense colonies. Proteomic analysis of these cells revealed a topologically-relevant interactome, focused on cytoskeletal and related enzymes//chaperone/signalling molecules linking DES to runx2 and alkaline phosphatase (ALP). Osteogenic differentiation led to mineralized woven bone nodules confined to dense colonies, significantly smaller and more circular with respect to controls. It significantly increased also colony-forming efficiency and the number of DES-immunoreactive dense colonies, and immunostaining of co-localized DES/runx-2 and DES/ALP. These data confirmed pre-osteoblastic and osteoblastic differentiation, woven bone formation, and mineralization, supporting DES as a player in the molecular pathway leading to the osteogenic fate of rBMSCs.

Conclusion

Immunocytochemical and morphometric studies coupled with proteomic and bioinformatic analysis support the concept that DES may act as an upstream signal for the skeletogenic commitment of rBMSCs. Thus, we suggest that altered metabolism of osteoblasts, woven bone, and mineralization by dysfunctional BMSCs might early be revealed by changes in DES expression//levels. Non-union fractures and endocrine - metabolic bone disorders like fibrous dysplasia and Paget disease of bone might take advantage of this molecular evidence for their early diagnosis and follow-up.

Electrostatic complementarity at the interface drives transient protein-protein interactions

Understanding the mechanisms driving bio-molecules binding and determining the resulting complexes' stability is fundamental for the prediction of binding regions, which is the starting point for drug-ability and design. Characteristics like the preferentially hydrophobic composition of the binding interfaces, the role of van der Waals interactions, and the consequent shape complementarity between the interacting molecular surfaces are well established. However, no consensus has yet been reached on the role of electrostatic. Here, we perform extensive analyses on a large dataset of protein complexes for which both experimental binding affinity and pH data were available. Probing the amino acid composition, the disposition of the charges, and the electrostatic potential they generated on the protein molecular surfaces, we found that (i) although different classes of dimers do not present marked differences in the amino acid composition and charges disposition in the binding region, (ii) homodimers with identical binding region show higher electrostatic compatibility with respect to both homodimers with non-identical binding region and heterodimers. Interestingly, (iii) shape and electrostatic complementarity, for patches defined on short-range interactions, behave oppositely when one stratifies the complexes by their binding affinity: complexes with higher binding affinity present high values of shape complementarity (the role of the Lennard-Jones potential predominates) while electrostatic tends to be randomly distributed. Conversely, complexes with low values of binding affinity exploit Coulombic complementarity to acquire specificity, suggesting that electrostatic complementarity may play a greater role in transient (or less stable) complexes. In light of these results, (iv) we provide a novel, fast, and efficient method, based on the 2D Zernike polynomial formalism, to measure electrostatic complementarity without the need of knowing the complex structure. Expanding the electrostatic potential on a basis of 2D orthogonal polynomials, we can discriminate between transient and permanent protein complexes with an AUC of the ROC of [Formula: see text] 0.8. Ultimately, our work helps shedding light on the non-trivial relationship between the hydrophobic and electrostatic contributions in the binding interfaces, thus favoring the development of new predictive methods for binding affinity characterization.

Selected Approaches to Disrupting Protein-Protein Interactions within the MAPK/RAS Pathway

Within the MAPK/RAS pathway, there exists a plethora of protein-protein interactions (PPIs). For many years, scientists have focused efforts on drugging KRAS and its effectors in hopes to provide much needed therapies for patients with KRAS-mutant driven cancers. In this review, we focus on recent strategies to inhibit RAS-signaling via disrupting PPIs associated with SOS1, RAF, PDEδ, Grb2, and RAS.

Computational analyses of the interactome between TNF and TNFR superfamilies

Proteins in the tumor necrosis factor (TNF) superfamily (TNFSF) regulate diverse cellular processes by interacting with their receptors in the TNF receptor (TNFR) superfamily (TNFRSF). Ligands and receptors in these two superfamilies form a complicated network of interactions, in which the same ligand can bind to different receptors and the same receptor can be shared by different ligands. In order to study these interactions on a systematic level, a TNFSF-TNFRSF interactome was constructed in this study by searching the database which consists of both experimentally measured and computationally predicted protein-protein interactions (PPIs). The interactome contains a total number of 194 interactions between 18 TNFSF ligands and 29 TNFRSF receptors in human. We modeled the structure for each ligand-receptor interaction in the network. Their binding affinities were further computationally estimated based on modeled structures. Our computational outputs, which are all publicly accessible, serve as a valuable addition to the currently limited experimental resources to study TNF-mediated cell signaling.

Data Incompleteness May form a Hard-to-Overcome Barrier to Decoding Life’s Mechanism

Simple Summary

The influence of data incompleteness on the correctness of conclusions about the structure and functions of the objects under study is widely discussed in the literature. It was noted that even a small percentage of missing data can lead to incorrect conclusions and imperfect knowledge. In particular, incompleteness can lead to critical errors in the qualitative and quantitative assessments of interactions in biological systems and a distorted understanding of the functioning mechanisms of living systems. In this brief review, we attempt to demonstrate the extent of this incompleteness in functional information about living systems using the best-studied examples. We suggest that this incompleteness may form seemingly insurmountable barriers in deciphering the mechanisms of the functioning of complex systems with unpredictable properties arising from the interaction of the system components.

Abstract

In this brief review, we attempt to demonstrate that the incompleteness of data, as well as the intrinsic heterogeneity of biological systems, may form very strong and possibly insurmountable barriers for researchers trying to decipher the mechanisms of the functioning of live systems. We illustrate this challenge using the two most studied organisms: E. coli, with 34.6% genes lacking experimental evidence of function, and C. elegans, with identified proteins for approximately 50% of its genes. Another striking example is an artificial unicellular entity named JCVI-syn3.0, with a minimal set of genes. A total of 31.5% of the genes of JCVI-syn3.0 cannot be ascribed a specific biological function. The human interactome mapping project identified only 5–10% of all protein interactions in humans. In addition, most of the available data are static snapshots, and it is barely possible to generate realistic models of the dynamic processes within cells. Moreover, the existing interactomes reflect the de facto interaction but not its functional result, which is an unpredictable emerging property. Perhaps the completeness of molecular data on any living organism is beyond our reach and represents an unsolvable problem in biology.

SIN-3 functions through multi-protein interaction to regulate apoptosis, autophagy, and longevity in Caenorhabditis elegans

SIN3/HDAC is a multi-protein complex that acts as a regulatory unit and functions as a co-repressor/co-activator and a general transcription factor. SIN3 acts as a scaffold in the complex, binding directly to HDAC1/2 and other proteins and plays crucial roles in regulating apoptosis, differentiation, cell proliferation, development, and cell cycle. However, its exact mechanism of action remains elusive. Using the Caenorhabditis elegans (C. elegans) model, we can surpass the challenges posed by the functional redundancy of SIN3 isoforms. In this regard, we have previously demonstrated the role of SIN-3 in uncoupling autophagy and longevity in C. elegans. In order to understand the mechanism of action of SIN3 in these processes, we carried out a comparative analysis of the SIN3 protein interactome from model organisms of different phyla. We identified conserved, expanded, and contracted gene classes. The C. elegans SIN-3 interactome -revealed the presence of  well-known proteins, such as DAF-16, SIR-2.1, SGK-1, and AKT-1/2, involved in autophagy, apoptosis, and longevity. Overall, our analyses propose  potential mechanisms by which SIN3 participates in multiple biological processes and their conservation across species and identifies candidate genes for further experimental analysis.

EXCESP: A Structure-Based Online Database for Extracellular Interactome of Cell Surface Proteins in Humans

The interactions between ectodomains of cell surface proteins are vital players in many important cellular processes, such as regulating immune responses, coordinating cell differentiation, and shaping neural plasticity. However, while the construction of a large-scale protein interactome has been greatly facilitated by the development of high-throughput experimental techniques, little progress has been made to support the discovery of extracellular interactome for cell surface proteins. Harnessed by the recent advances in computational modeling of protein-protein interactions, here we present a structure-based online database for the extracellular interactome of cell surface proteins in humans, called EXCESP. The database contains both experimentally determined and computationally predicted interactions among all type-I transmembrane proteins in humans. All structural models for these interactions and their binding affinities were further computationally modeled. Moreover, information such as expression levels of each protein in different cell types and its relation to various signaling pathways from other online resources has also been integrated into the database. In summary, the database serves as a valuable addition to the existing online resources for the study of cell surface proteins. It can contribute to the understanding of the functions of cell surface proteins in the era of systems biology.

Classification of protein-protein association rates based on biophysical informatics

BACKGROUND

Proteins form various complexes to carry out their versatile functions in cells. The dynamic properties of protein complex formation are mainly characterized by the association rates which measures how fast these complexes can be formed. It was experimentally observed that the association rates span an extremely wide range with over ten orders of magnitudes. Identification of association rates within this spectrum for specific protein complexes is therefore essential for us to understand their functional roles.

RESULTS

To tackle this problem, we integrate physics-based coarse-grained simulations into a neural-network-based classification model to estimate the range of association rates for protein complexes in a large-scale benchmark set. The cross-validation results show that, when an optimal threshold was selected, we can reach the best performance with specificity, precision, sensitivity and overall accuracy all higher than 70%. The quality of our cross-validation data has also been testified by further statistical analysis. Additionally, given an independent testing set, we can successfully predict the group of association rates for eight protein complexes out of ten. Finally, the analysis of failed cases suggests the future implementation of conformational dynamics into simulation can further improve model.

CONCLUSIONS

In summary, this study demonstrated that a new modeling framework that combines biophysical simulations with bioinformatics approaches is able to identify protein-protein interactions with low association rates from those with higher association rates. This method thereby can serve as a useful addition to a collection of existing experimental approaches that measure biomolecular recognition.

Regulatory role of local tissue signal Del-1 in cancer and inflammation: a review

Developmental endothelial locus-1 (Del-1) is a secretory, multifunctional domain protein. It can bind to integrins and phosphatidylserine. As a local tissue signal, it plays a regulatory role in the cancer microenvironment and inflammation. Del-1 has destructive effects in most cancers and is associated with the progression and invasion of some cancers. In contrast, Del-1 also plays a protective role in inflammation. Del-1 regulates inflammation by regulating the generation of neutrophils in bone marrow, inhibiting the recruitment and migration of neutrophils and accelerating the clearance of neutrophils by macrophages. Del-1 and IL-17 are reciprocally regulated, and their balance maintains immune system homeostasis. Del-1 is expected to become a new therapeutic target for inflammatory disorders such as multiple sclerosis.

Computational Assessment of Protein-Protein Binding Affinity by Reverse Engineering the Energetics in Protein Complexes

The cellular functions of proteins are maintained by forming diverse complexes. The stability of these complexes is quantified by the measurement of binding affinity, and mutations that alter the binding affinity can cause various diseases such as cancer and diabetes. As a result, accurate estimation of the binding stability and the effects of mutations on changes of binding affinity is a crucial step to understanding the biological functions of proteins and their dysfunctional consequences. It has been hypothesized that the stability of a protein complex is dependent not only on the residues at its binding interface by pairwise interactions but also on all other remaining residues that do not appear at the binding interface. Here, we computationally reconstruct the binding affinity by decomposing it into the contributions of interfacial residues and other non-interfacial residues in a protein complex. We further assume that the contributions of both interfacial and non-interfacial residues to the binding affinity depend on their local structural environments such as solvent-accessible surfaces and secondary structural types. The weights of all corresponding parameters are optimized by Monte-Carlo simulations. After cross-validation against a large-scale dataset, we show that the model not only shows a strong correlation between the absolute values of the experimental and calculated binding affinities, but can also be an effective approach to predict the relative changes of binding affinity from mutations. Moreover, we have found that the optimized weights of many parameters can capture the first-principle chemical and physical features of molecular recognition, therefore reversely engineering the energetics of protein complexes. These results suggest that our method can serve as a useful addition to current computational approaches for predicting binding affinity and understanding the molecular mechanism of protein–protein interactions.

2D Zernike polynomial expansion: Finding the protein-protein binding regions

We present a method for efficiently and effectively assessing whether and where two proteins can interact with each other to form a complex. This is still largely an open problem, even for those relatively few cases where the 3D structure of both proteins is known. In fact, even if much of the information about the interaction is encoded in the chemical and geometric features of the structures, the set of possible contact patches and of their relative orientations are too large to be computationally affordable in a reasonable time, thus preventing the compilation of reliable interactome. Our method is able to rapidly and quantitatively measure the geometrical shape complementarity between interacting proteins, comparing their molecular iso-electron density surfaces expanding the surface patches in term of 2D Zernike polynomials. We first test the method against the real binding region of a large dataset of known protein complexes, reaching a success rate of 0.72. We then apply the method for the blind recognition of binding sites, identifying the real region of interaction in about 60% of the analyzed cases. Finally, we investigate how the efficiency in finding the right binding region depends on the surface roughness as a function of the expansion order.

Understanding the Impacts of Conformational Dynamics on the Regulation of Protein-Protein Association by a Multiscale Simulation Method

Complexes formed among diverse proteins carry out versatile functions in nearly all physiological processes. Association rates which measure how fast proteins form various complexes are of fundamental importance to characterize their functions. The association rates are not only determined by the energetic features at binding interfaces of a protein complex but also influenced by the intrinsic conformational dynamics of each protein in the complex. Unfortunately, how this conformational effect regulates protein association has never been calibrated on a systematic level. To tackle this problem, we developed a multiscale strategy to incorporate the information on protein conformational variations from Langevin dynamic simulations into a kinetic Monte Carlo algorithm of protein-protein association. By systematically testing this approach against a large-scale benchmark set, we found the association of a protein complex with a relatively rigid structure tends to be reduced by its conformational fluctuations. With specific examples, we further show that higher degrees of structural flexibility in various protein complexes can facilitate the searching and formation of intermolecular interactions and thereby accelerate their associations. In general, the integration of conformational dynamics can improve the correlation between experimentally measured association rates and computationally derived association probabilities. Finally, we analyzed the statistical distributions of different secondary structural types on protein-protein binding interfaces and their preference to the change of association rates. Our study, to the best of our knowledge, is the first computational method that systematically estimates the impacts of protein conformational dynamics on protein-protein association. It throws lights on the molecular mechanisms of how protein-protein recognition is kinetically modulated.

Using Coarse-Grained Simulations to Characterize the Mechanisms of Protein-Protein Association

The formation of functionally versatile protein complexes underlies almost every biological process. The estimation of how fast these complexes can be formed has broad implications for unravelling the mechanism of biomolecular recognition. This kinetic property is traditionally quantified by association rates, which can be measured through various experimental techniques. To complement these time-consuming and labor-intensive approaches, we developed a coarse-grained simulation approach to study the physical processes of protein–protein association. We systematically calibrated our simulation method against a large-scale benchmark set. By combining a physics-based force field with a statistically-derived potential in the simulation, we found that the association rates of more than 80% of protein complexes can be correctly predicted within one order of magnitude relative to their experimental measurements. We further showed that a mixture of force fields derived from complementary sources was able to describe the process of protein–protein association with mechanistic details. For instance, we show that association of a protein complex contains multiple steps in which proteins continuously search their local binding orientations and form non-native-like intermediates through repeated dissociation and re-association. Moreover, with an ensemble of loosely bound encounter complexes observed around their native conformation, we suggest that the transition states of protein–protein association could be highly diverse on the structural level. Our study also supports the idea in which the association of a protein complex is driven by a “funnel-like” energy landscape. In summary, these results shed light on our understanding of how protein–protein recognition is kinetically modulated, and our coarse-grained simulation approach can serve as a useful addition to the existing experimental approaches that measure protein–protein association rates.

TCF 4 tumor suppressor: a molecular target in the prognosis of sporadic colorectal cancer in humans

Background

A huge array of function is played by the Wnt/β-catenin signaling pathway in development by balancing gene expression through the modulation of cell-specific DNA binding downstream effectors such as T-cell factor/lymphoid enhancer factor (TCF/LEF). The β-catenin/TCF-4 complex is a central regulatory switch for differentiation and proliferation of intestinal cells (both normal and malignant). Thus, in the present study we evaluated each of 60 cases of sporadic adenocarcinoma, alongside adjoining and normal mucosa specimens of colorectum in humans, for mutation and expression analysis of the gene coding for TCF-4 protein.

Methods

DNA sequencing following PCR amplification and SSCP analysis (single strand conformation polymorphism) was employed to detect TCF-4 gene mutations in the case of exon 1. Quantitative real-time (qRT) PCR, immunohistochemistry (IHC), confocal microscopy and western blot analysis were used to detect TCF-4 gene/protein expression.

Results

Sequencing analysis confirmed 5/60 patients with a point mutation in exon 1 of the TCF-4 gene in tumor samples. mRNA expression using qRT-PCR showed approximately 83% decreased TCF-4 mRNA expression in tumor tissue and adjoining mucosa compared to normal mucosa. Similarly, a significant decrease in protein expression using IHC showed decreased TCF-4 protein expression in tumor tissue and adjoining mucosa compared to normal mucosa, which also corresponds to some important clinicopathological factors, including disease metastasis and tumor grade. Mutational alterations and downregulation of TCF-4 mRNA and hence decreased expression of TCF-4 protein in tumors suggest its involvement in the pathogenesis of CRC.

Conclusions

A remarkable decrease in TCF-4 mRNA and protein expression was detected in tumorous and adjoining tissues compared to normal mucosa. Hence the alterations in genomic architecture along with downregulation of TCF-4 mRNA and decreased expression of TCF-4 protein in tumors, which is in accordance with clinical features, suggest its involvement in the pathogenesis of CRC. Thus, deregulation and collaboration of TCF-4 with CRC could be a concrete and distinctive feature in the prognosis of the disease at an early stage of development.

Computational studies of protein-protein dissociation by statistical potential and coarse-grained simulations: a case study on interactions between colicin E9 endonuclease and immunity proteins

Proteins carry out their diverse functions in cells by forming interactions with each other. The dynamics of these interactions are quantified by the measurement of association and dissociation rate constants. Relative to the efforts made to model the association of biomolecules, little has been studied to understand the principles of protein complex dissociation. Using the interaction between colicin E9 endonucleases and immunity proteins as a test system, here we develop a coarse-grained simulation method to explore the dissociation mechanisms of protein complexes. The interactions between proteins in the complex are described by the knowledge-based potential that was constructed by the statistics from available protein complexes in the structural database. Our study provides the supportive evidences to the dual recognition mechanism for the specificity of binding between E9 DNase and immunity proteins, in which the conserved residues of helix III of Im2 and Im9 proteins act as the anchor for binding, while the sequence variations in helix II make positive or negative contributions to specificity. Beyond that, we further suggest that this binding specificity is rooted in the process of complex dissociation instead of association. While we increased the flexibility of protein complexes, we further found that they are less prone to dissociation, suggesting that conformational fluctuations of protein complexes play important functional roles in regulating their binding and dissociation. Our studies therefore bring new insights to the molecule mechanisms of protein-protein interactions, while the method can serve as a new addition to a suite of existing computational tools for the simulations of protein complexes.

Integrating Structural Information to Study the Dynamics of Protein-Protein Interactions in Cells

The information of how two proteins interact is embedded in the atomic details of their binding interfaces. These interactions, spatial-temporally coordinating each other as a network in a variable cytoplasmic environment, dominate almost all biological functions. A feasible and reliable computational model is highly demanded to realistically simulate these cellular processes and unravel the complexities beneath them. We therefore present a multiscale framework that integrates simulations on two different scales. The higher-resolution model incorporates structural information of proteins and energetics of their binding, while the lower-resolution model uses a highly simplified representation of proteins to capture the long-time-scale dynamics of a system with multiple proteins. Through a systematic benchmark test and two practical applications of biomolecular systems with specific cellular functions, we demonstrated that this method could be a powerful approach to understand molecular mechanisms of dynamic interactions between biomolecules and their functional impacts with high computational efficiency.

Predicting Protein-protein Association Rates using Coarse-grained Simulation and Machine Learning

Protein–protein interactions dominate all major biological processes in living cells. We have developed a new Monte Carlo-based simulation algorithm to study the kinetic process of protein association. We tested our method on a previously used large benchmark set of 49 protein complexes. The predicted rate was overestimated in the benchmark test compared to the experimental results for a group of protein complexes. We hypothesized that this resulted from molecular flexibility at the interface regions of the interacting proteins. After applying a machine learning algorithm with input variables that accounted for both the conformational flexibility and the energetic factor of binding, we successfully identified most of the protein complexes with overestimated association rates and improved our final prediction by using a cross-validation test. This method was then applied to a new independent test set and resulted in a similar prediction accuracy to that obtained using the training set. It has been thought that diffusion-limited protein association is dominated by long-range interactions. Our results provide strong evidence that the conformational flexibility also plays an important role in regulating protein association. Our studies provide new insights into the mechanism of protein association and offer a computationally efficient tool for predicting its rate.

Understand protein functions by comparing the similarity of local structural environments

The three-dimensional structures of proteins play an essential role in regulating binding between proteins and their partners, offering a direct relationship between structures and functions of proteins. It is widely accepted that the function of a protein can be determined if its structure is similar to other proteins whose functions are known. However, it is also observed that proteins with similar global structures do not necessarily correspond to the same function, while proteins with very different folds can share similar functions. This indicates that function similarity is originated from the local structural information of proteins instead of their global shapes. We assume that proteins with similar local environments prefer binding to similar types of molecular targets. In order to testify this assumption, we designed a new structural indicator to define the similarity of local environment between residues in different proteins. This indicator was further used to calculate the probability that a given residue binds to a specific type of structural neighbors, including DNA, RNA, small molecules and proteins. After applying the method to a large-scale non-redundant database of proteins, we show that the positive signal of binding probability calculated from the local structural indicator is statistically meaningful. In summary, our studies suggested that the local environment of residues in a protein is a good indicator to recognize specific binding partners of the protein. The new method could be a potential addition to a suite of existing template-based approaches for protein function prediction.

A computational interactome and functional annotation for the human proteome

We present a database, PrePPI (Predicting Protein-Protein Interactions), of more than 1.35 million predicted protein-protein interactions (PPIs). Of these at least 127,000 are expected to constitute direct physical interactions although the actual number may be much larger (~500,000). The current PrePPI, which contains predicted interactions for about 85% of the human proteome, is related to an earlier version but is based on additional sources of interaction evidence and is far larger in scope. The use of structural relationships allows PrePPI to infer numerous previously unreported interactions. PrePPI has been subjected to a series of validation tests including reproducing known interactions, recapitulating multi-protein complexes, analysis of disease associated SNPs, and identifying functional relationships between interacting proteins. We show, using Gene Set Enrichment Analysis (GSEA), that predicted interaction partners can be used to annotate a protein's function. We provide annotations for most human proteins, including many annotated as having unknown function.

Decomposing the space of protein quaternary structures with the interface fragment pair library

BACKGROUND

The physical interactions between proteins constitute the basis of protein quaternary structures. They dominate many biological processes in living cells. Deciphering the structural features of interacting proteins is essential to understand their cellular functions. Similar to the space of protein tertiary structures in which discrete patterns are clearly observed on fold or sub-fold motif levels, it has been found that the space of protein quaternary structures is highly degenerate due to the packing of compact secondary structure elements at interfaces. Therefore, it is necessary to further decompose the protein quaternary structural space into a more local representation.

RESULTS

Here we constructed an interface fragment pair library from the current structure database of protein complexes. After structural-based clustering, we found that more than 90% of these interface fragment pairs can be represented by a limited number of highly abundant motifs. These motifs were further used to guide complex assembly. A large-scale benchmark test shows that the native-like binding is highly likely in the structural ensemble of modeled protein complexes that were built through the library.

CONCLUSIONS

Our study therefore presents supportive evidences that the space of protein quaternary structures can be represented by the combination of a small set of secondary-structure-based packing at binding interfaces. Finally, after future improvements such as adding sequence profiles, we expect this new library will be useful to predict structures of unknown protein-protein interactions.

Plant protein interactomes

Protein-protein interactions are a critical element of biological systems, and the analysis of interaction partners can provide valuable hints about unknown functions of a protein. In recent years, several large-scale protein interaction studies have begun to unravel the complex networks through which plant proteins exert their functions. Two major classes of experimental approaches are used for protein interaction mapping: analysis of direct interactions using binary methods such as yeast two-hybrid or split ubiquitin, and analysis of protein complexes through affinity purification followed by mass spectrometry. In addition, bioinformatics predictions can suggest interactions that have evaded detection by other methods or those of proteins that have not been investigated. Here we review the major approaches to construct, analyze, use, and carry out quality control on plant protein interactome networks. We present experimental and computational approaches for large-scale mapping, methods for validation or smaller-scale functional studies, important bioinformatics resources, and findings from recently published large-scale plant interactome network maps.

The effect of hyperosmolar stimuli and cyclophosphamide on the culture of normal rat urothelial cells in vitro

Highly concentrated urine may induce a harmful effect on the urinary bladder. Therefore, we considered osmolarity of the urine as a basic pathomechanism of mucosal damage. The influence of both cyclophosphamide (CYP) and hyperosmolar stimuli (HS) on the urothelium are not well described. The purpose was to evaluate the effect of CYP and HS on rat urothelial cultured cells (RUCC). 15 Wistar rats were used for RUCC preparation. RUCC were exposed to HS (2080 and 3222 mOsm/l NaCl) for 15 min and CYP (1 mg/ml) for 4 hrs. APC-labelled annexin V was used to quantitatively determine the percentage of apoptotic cells and propidium iodide (PI) as a standard flow cytometric viability probe to distinguish necrotic cells from viable ones. Annexin V-APC (+), annexin V-APC and PI (+), and PI (+) cells were analysed as apoptotic, dead, and necrotic cells, respectively. The results were presented in percentage values. The flow cytometric analysis was done on a FACSCalibur Flow Cytometer using Cell-Quest software. Treatment with 2080 and 3222 mOsm/l HS resulted in 23.7 ± 3.9% and 26.0 ± 1.5% apoptotic cells, respectively, 14.3 ± 1.4% and 19.4 ± 2.7% necrotic cells, respectively and 60.5 ± 1.4% and 48.6 ± 5.3% dead cells, respectively. The effect of CYP on RUCC was similar to the effect of HS. After CYP the apoptotic and necrotic cells were 23.1 ± 0.3% and 17.9 ± 7.4%, respectively. The percentage of dead cells was 57.7 ± 10.8%. CYP and HS induced apoptosis and necrosis in RUCC. 3222 mOsm/l HS had the most harmful effect based on the percentage of necrotic and apoptotic cells.

Network approaches to the functional analysis of microbial proteins

Large amounts of detailed biological data have been generated over the past few decades. Much of these data is freely available in over 1000 online databases; an enticing, but frustrating resource for microbiologists interested in a systems-level view of the structure and function of microbial cells. The frustration engendered by the need to trawl manually through hundreds of databases in order to accumulate information about a gene, protein, pathway, or organism of interest can be alleviated by the use of computational data integration to generated network views of the system of interest. Biological networks can be constructed from a single type of data, such as protein-protein binding information, or from data generated by multiple experimental approaches. In an integrated network, nodes usually represent genes or gene products, while edges represent some form of interaction between the nodes. Edges between nodes may be weighted to represent the probability that the edge exists in vivo. Networks may also be enriched with ontological annotations, facilitating both visual browsing and computational analysis via web service interfaces. In this review, we describe the construction, analysis of both single-data source and integrated networks, and their application to the inference of protein function in microbes.

Interactome networks and human disease

Complex biological systems and cellular networks may underlie most genotype to phenotype relationships. Here, we review basic concepts in network biology, discussing different types of interactome networks and the insights that can come from analyzing them. We elaborate on why interactome networks are important to consider in biology, how they can be mapped and integrated with each other, what global properties are starting to emerge from interactome network models, and how these properties may relate to human disease.

GIDMP: good protein-protein interaction data metamining practice

Studying the interactome is one of the exciting frontiers of proteomics, as shown lately at the recent bioinformatics conferences (for example ISMB 2010, or ECCB 2010). Distribution of data is facilitated by a large number of databases. Metamining databases have been created in order to allow researchers access to several databases in one search, but there are serious difficulties for end users to evaluate the metamining effort. Therefore we suggest a new standard, "Good Interaction Data Metamining Practice" (GIDMP), which could be easily automated and requires only very minor inclusion of statistical data on each database homepage. Widespread adoption of the GIDMP standard would provide users with: a standardized way to evaluate the statistics provided by each metamining database, thus enhancing the end-user experience; a stable contact point for each database, allowing the smooth transition of statistics; a fully automated system, enhancing time- and cost-effectiveness. The proposed information can be presented as a few hidden lines of text on the source database www page, and a constantly updated table for a metamining database included in the source/credits web page.

PPI_SVM: prediction of protein-protein interactions using machine learning, domain-domain affinities and frequency tables

Protein-protein interactions (PPI) control most of the biological processes in a living cell. In order to fully understand protein functions, a knowledge of protein-protein interactions is necessary. Prediction of PPI is challenging, especially when the three-dimensional structure of interacting partners is not known. Recently, a novel prediction method was proposed by exploiting physical interactions of constituent domains. We propose here a novel knowledge-based prediction method, namely PPI_SVM, which predicts interactions between two protein sequences by exploiting their domain information. We trained a two-class support vector machine on the benchmarking set of pairs of interacting proteins extracted from the Database of Interacting Proteins (DIP). The method considers all possible combinations of constituent domains between two protein sequences, unlike most of the existing approaches. Moreover, it deals with both single-domain proteins and multi domain proteins; therefore it can be applied to the whole proteome in high-throughput studies. Our machine learning classifier, following a brainstorming approach, achieves accuracy of 86%, with specificity of 95%, and sensitivity of 75%, which are better results than most previous methods that sacrifice recall values in order to boost the overall precision. Our method has on average better sensitivity combined with good selectivity on the benchmarking dataset. The PPI_SVM source code, train/test datasets and supplementary files are available freely in the public domain at: http://code.google.com/p/cmater-bioinfo/.

PRIN: a predicted rice interactome network

Background

Protein-protein interactions play a fundamental role in elucidating the molecular mechanisms of biomolecular function, signal transductions and metabolic pathways of living organisms. Although high-throughput technologies such as yeast two-hybrid system and affinity purification followed by mass spectrometry are widely used in model organisms, the progress of protein-protein interactions detection in plants is rather slow. With this motivation, our work presents a computational approach to predict protein-protein interactions in Oryza sativa.

Results

To better understand the interactions of proteins in Oryza sativa, we have developed PRIN, a Predicted Rice Interactome Network. Protein-protein interaction data of PRIN are based on the interologs of six model organisms where large-scale protein-protein interaction experiments have been applied: yeast (Saccharomyces cerevisiae), worm (Caenorhabditis elegans), fruit fly (Drosophila melanogaster), human (Homo sapiens), Escherichia coli K12 and Arabidopsis thaliana. With certain quality controls, altogether we obtained 76,585 non-redundant rice protein interaction pairs among 5,049 rice proteins. Further analysis showed that the topology properties of predicted rice protein interaction network are more similar to yeast than to the other 5 organisms. This may not be surprising as the interologs based on yeast contribute nearly 74% of total interactions. In addition, GO annotation, subcellular localization information and gene expression data are also mapped to our network for validation. Finally, a user-friendly web interface was developed to offer convenient database search and network visualization.

Conclusions

PRIN is the first well annotated protein interaction database for the important model plant Oryza sativa. It has greatly extended the current available protein-protein interaction data of rice with a computational approach, which will certainly provide further insights into rice functional genomics and systems biology.

PRIN is available online at http://bis.zju.edu.cn/prin/.

Computational approaches for protein function prediction: a combined strategy from multiple sequence alignment to molecular docking-based virtual screening

The functional characterization of proteins represents a daily challenge for biochemical, medical and computational sciences. Although finally proved on the bench, the function of a protein can be successfully predicted by computational approaches that drive the further experimental assays. Current methods for comparative modeling allow the construction of accurate 3D models for proteins of unknown structure, provided that a crystal structure of a homologous protein is available. Binding regions can be proposed by using binding site predictors, data inferred from homologous crystal structures, and data provided from a careful interpretation of the multiple sequence alignment of the investigated protein and its homologs. Once the location of a binding site has been proposed, chemical ligands that have a high likelihood of binding can be identified by using ligand docking and structure-based virtual screening of chemical libraries. Most docking algorithms allow building a list sorted by energy of the lowest energy docking configuration for each ligand of the library. In this review the state-of-the-art of computational approaches in 3D protein comparative modeling and in the study of protein-ligand interactions is provided. Furthermore a possible combined/concerted multistep strategy for protein function prediction, based on multiple sequence alignment, comparative modeling, binding region prediction, and structure-based virtual screening of chemical libraries, is described by using suitable examples. As practical examples, Abl-kinase molecular modeling studies, HPV-E6 protein multiple sequence alignment analysis, and some other model docking-based characterization reports are briefly described to highlight the importance of computational approaches in protein function prediction.

Computer applications for prediction of protein-protein interactions and rational drug design

In recent years, protein-protein interactions are becoming the object of increasing attention in many different fields, such as structural biology, molecular biology, systems biology, and drug discovery. From a structural biology perspective, it would be desirable to integrate current efforts into the structural proteomics programs. Given that experimental determination of many protein-protein complex structures is highly challenging, and in the context of current high-performance computational capabilities, different computer tools are being developed to help in this task. Among them, computational docking aims to predict the structure of a protein-protein complex starting from the atomic coordinates of its individual components, and in recent years, a growing number of docking approaches are being reported with increased predictive capabilities. The improvement of speed and accuracy of these docking methods, together with the modeling of the interaction networks that regulate the most critical processes in a living organism, will be essential for computational proteomics. The ultimate goal is the rational design of drugs capable of specifically inhibiting or modifying protein-protein interactions of therapeutic significance. While rational design of protein-protein interaction inhibitors is at its very early stage, the first results are promising.

The era of 'omics unlimited

Determining the primary sequences of informational macromolecules is no longer a limiting factor for our ability to completely understand the biological functioning of cells and organisms. Similarly, our understanding of transcriptional regulation (transcriptomics) has been greatly enhanced by the availability of microarrays. Our next hurdle is to learn the biochemical functions of all the gene products (proteomics) and the totality of all the interactions among them (interactomics). Using traditional biochemical methods, this will take a very long time. More efficient methods are needed to address these questions, or at least to suggest possible candidates for further testing. High-resolution imaging using molecule-specific tags will reveal details of cellular architecture that are expected to provide additional insights and clues about the interactions and functions of many gene products. Computer modeling of macromolecular structures and functional systems will be of key importance. We present here a brief historical and futuristic perspective of genomics and some of its other ‘omics offshoots in the post-genomic era.