Protein coadaptation and the design of novel approaches to identify protein-protein interactions

Proteins rarely function in isolation but they form part of complex networks of interactions with other proteins within or among cells. The importance of a particular protein for cell viability is directly dependent upon the number of interactions where it participates and the function it performs: the larger the number of interactions of a protein the greater its functional importance is for the cell. With the advent of genome sequencing and "omics" technologies it became feasible conducting large-scale searches for protein interacting partners. Unfortunately, the accuracy of such analyses has been underwhelming owing to methodological limitations and to the inherent complexity of protein interactions. In addition to these experimental approaches, many computational methods have been developed to identify protein-protein interactions by assuming that interacting proteins coevolve resulting from the coadaptation dynamics between the amino acids of their interacting faces. We review the main technological advances made in the field of interactomics and discuss the feasibility of computational methods to identify protein-protein interactions based on the estimation of coevolution. As proof-of-concept, we present a classical case study: the interactions of cell surface proteins (receptors) and their ligands. Finally, we take this discussion one step forward to include interactions between organisms and species to understand the generation of biological complexity. Development of technologies for accurate detection of protein-protein interactions may shed light on processes that go from the fine-tuning of pathways and metabolic networks to the emergence of biological complexity.     

Systems approach to explore components and interactions in the presynapse

The application of proteomic techniques to neuroscientific research provides an opportunity for a greater understanding of nervous system structure and function. As increasing amounts of neuroproteomic data become available, it is necessary to formulate methods to integrate these data in a meaningful way to obtain a more comprehensive picture of neuronal subcompartments. Furthermore, computational methods can be used to make biologically relevant predictions from large proteomic data sets. Here, we applied an integrated proteomics and systems biology approach to characterize the presynaptic (PRE) nerve terminal. For this, we carried out proteomic analyses of presynaptically enriched fractions, and generated a PRE literature‐based protein–protein interaction network. We combined these with other proteomic analyses to generate a core list of 117 PRE proteins, and used graph theory‐inspired algorithms to predict 92 additional components and a PRE complex containing 17 proteins. Some of these predictions were validated experimentally, indicating that the computational analyses can identify novel proteins and complexes in a subcellular compartment. We conclude that the combination of techniques (proteomics, data integration, and computational analyses) used in this study are useful in obtaining a comprehensive understanding of functional components, especially low‐abundance entities and/or interactions in the PRE nerve terminal.     

Targeted purification development enabled by computational biophysical modeling

Chromatographic and non‐chromatographic purification of biopharmaceuticals depend on the interactions between protein molecules and a solid–liquid interface. These interactions are dominated by the protein–surface properties, which are a function of protein sequence, structure, and dynamics. In addition, protein–surface properties are critical for in vivo recognition and activation, thus, purification strategies should strive to preserve structural integrity and retain desired pharmacological efficacy. Other factors such as surface diffusion, pore diffusion, and film mass transfer can impact chromatographic separation and resin design. The key factors that impact non‐chromatographic separations (e.g., solubility, ligand affinity, charges and hydrophobic clusters, and molecular dynamics) are readily amenable to computational modeling and can enhance the understanding of protein chromatographic. Previously published studies have used computational methods such as quantitative structure–activity relationship (QSAR) or quantitative structure–property relationship (QSPR) to identify and rank order affinity ligands based on their potential to effectively bind and separate a desired biopharmaceutical from host cell protein (HCP) and other impurities. The challenge in the application of such an approach is to discern key yet subtle differences in ligands and proteins that influence biologics purification. Using a relatively small molecular weight protein (insulin), this research overcame limitations of previous modeling efforts by utilizing atomic level detail for the modeling of protein–ligand interactions, effectively leveraging and extending previous research on drug target discovery. These principles were applied to the purification of different commercially available insulin variants. The ability of these computational models to correlate directionally with empirical observation is demonstrated for several insulin systems over a range of purification challenges including resolution of subtle product variants (amino acid misincorporations). Broader application of this methodology in bioprocess development may enhance and speed the development of a robust purification platform. © 2014 American Institute of Chemical Engineers Biotechnol. Prog., 31:154–164, 2015     

Computational alanine scanning with linear scaling semiempirical quantum mechanical methods

Alanine scanning is a powerful experimental tool for understanding the key interactions in protein–protein interfaces. Linear scaling semiempirical quantum mechanical calculations are now sufficiently fast and robust to allow meaningful calculations on large systems such as proteins, RNA and DNA. In particular, they have proven useful in understanding protein–ligand interactions. Here we ask the question: can these linear scaling quantum mechanical methods developed for protein–ligand scoring be useful for computational alanine scanning? To answer this question, we assembled 15 protein–protein complexes with available crystal structures and sufficient alanine scanning data. In all, the data set contains ΔΔGs for 400 single point alanine mutations of these 15 complexes. We show that with only one adjusted parameter the quantum mechanics‐based methods outperform both buried accessible surface area and a potential of mean force and compare favorably to a variety of published empirical methods. Finally, we closely examined the outliers in the data set and discuss some of the challenges that arise from this examination. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Characterization of protein–glycolipid recognition at the membrane bilayer

A growing number of important molecular recognition events are being shown to involve the interactions between proteins and glycolipids. Glycolipids are molecules in which one or more monosaccharides are glycosidically linked to a lipid moiety. The lipid moiety is generally buried in the cell membrane or other bilayer, leaving the oligosaccharide moiety exposed but in close proximity to the bilayer surface. This presents a unique environment for protein–carbohydrate interactions, and studies to determine the influence of the bilayer on these phenomena are in their infancy. One important property of the bilayer is the ability to orient and cluster glycolipid species, as strong interactions in biological systems are often achieved through multivalency arising from the simultaneous association of two or more proteins and receptors. This is especially true of protein–carbohydrate binding because of the unusually low affinities that characterize the monovalent interactions. More recent studies have also shown that the composition of the lipid bilayer is a critical parameter in protein–glycolipid recognition. The fluidity of the bilayer allows for correct geometric positioning of the oligosaccharide head group relative to the binding sites on the protein. In addition, there are activity‐based and structural data demonstrating the impact of the bilayer microenvironment on the modulation of oligosaccharide presentation. The use of model membranes in biosensor‐based methods has supplied decisive evidence of the importance of the membrane in receptor presentation. These data can be correlated with three‐dimensional structural information from X‐ray <?tw=98%>crystallography, NMR, and molecular mechanics to provide insight into specific protein–carbohydrate inter‐­actions at the bilayer. Copyright © 1999 National Research Council Canada and John Wiley & Sons, Ltd.     

Computational Prediction of Protein–Protein Interactions

Recently a number of computational approaches have been developed for the prediction of protein–protein interactions. Complete genome sequencing projects have provided the vast amount of information needed for these analyses. These methods utilize the structural, genomic, and biological context of proteins and genes in complete genomes to predict protein interaction networks and functional linkages between proteins. Given that experimental techniques remain expensive, time-consuming, and labor-intensive, these methods represent an important advance in proteomics. Some of these approaches utilize sequence data alone to predict interactions, while others combine multiple computational and experimental datasets to accurately build protein interaction maps for complete genomes. These methods represent a complementary approach to current high-throughput projects whose aim is to delineate protein interaction maps in complete genomes. We will describe a number of computational protocols for protein interaction prediction based on the structural, genomic, and biological context of proteins in complete genomes, and detail methods for protein interaction network visualization and analysis.     

Systematic prediction of human membrane receptor interactions

Membrane receptor‐activated signal transduction pathways are integral to cellular functions and disease mechanisms in humans. Identification of the full set of proteins interacting with membrane receptors by high‐throughput experimental means is difficult because methods to directly identify protein interactions are largely not applicable to membrane proteins. Unlike prior approaches that attempted to predict the global human interactome, we used a computational strategy that only focused on discovering the interacting partners of human membrane receptors leading to improved results for these proteins. We predict specific interactions based on statistical integration of biological data containing highly informative direct and indirect evidences together with feedback from experts. The predicted membrane receptor interactome provides a system‐wide view, and generates new biological hypotheses regarding interactions between membrane receptors and other proteins. We have experimentally validated a number of these interactions. The results suggest that a framework of systematically integrating computational predictions, global analyses, biological experimentation and expert feedback is a feasible strategy to study the human membrane receptor interactome.     

Identifying responsive functional modules from protein-protein interaction network

Proteins interact with each other within a cell, and those interactions give rise to the biological function and dynamical behavior of cellular systems. Generally, the protein interactions are temporal, spatial, or condition dependent in a specific cell, where only a small part of interactions usually take place under certain conditions. Recently, although a large amount of protein interaction data have been collected by high-throughput technologies, the interactions are recorded or summarized under various or different conditions and therefore cannot be directly used to identify signaling pathways or active networks, which are believed to work in specific cells under specific conditions. However, protein interactions activated under specific conditions may give hints to the biological process underlying corresponding phenotypes. In particular, responsive functional modules consist of protein interactions activated under specific conditions can provide insight into the mechanism underlying biological systems, e.g. protein interaction subnetworks found for certain diseases rather than normal conditions may help to discover potential biomarkers. From computational viewpoint, identifying responsive functional modules can be formulated as an optimization problem. Therefore, efficient computational methods for extracting responsive functional modules are strongly demanded due to the NP-hard nature of such a combinatorial problem. In this review, we first report recent advances in development of computational methods for extracting responsive functional modules or active pathways from protein interaction network and microarray data. Then from computational aspect, we discuss remaining obstacles and perspectives for this attractive and challenging topic in the area of systems biology.     

Uncovering the In Vivo Proxisome Using Proximity‐Tagging Methods

The development of new approaches is critical to gain further insights into biological processes that cannot be obtained by existing methods or technologies. The detection of protein–protein interaction is often challenging, especially for weak and transient interactions or for membrane proteins. Over the last decade, several proximity‐tagging methodologies have been developed to explore protein interactions in living cells. Among those, the most efficient are based on protein partner modification, such as biotinylation or pupylation. Such technologies are based on engineered variants of enzymes like peroxidases or ligases that release reactive molecules, in the presence of specific substrates, that bind surrounding proteins. Fusing a protein of interest (POI) to these enzymes allows the definition of an unbiased “proxisome,” that is, all of the proteins in interaction or in close vicinity of the POI. Here, the different proximity‐labeling tools available are described and comprehensive comparison to discuss advantages and limitations is provided.     

Structure- and sequence-based function prediction for non-homologous proteins

The structural genomics projects have been accumulating an increasing number of protein structures, many of which remain functionally unknown. In parallel effort to experimental methods, computational methods are expected to make a significant contribution for functional elucidation of such proteins. However, conventional computational methods that transfer functions from homologous proteins do not help much for these uncharacterized protein structures because they do not have apparent structural or sequence similarity with the known proteins. Here, we briefly review two avenues of computational function prediction methods, i.e. structure-based methods and sequence-based methods. The focus is on our recent developments of local structure-based and sequence-based methods, which can effectively extract function information from distantly related proteins. Two structure-based methods, Pocket-Surfer and Patch-Surfer, identify similar known ligand binding sites for pocket regions in a query protein without using global protein fold similarity information. Two sequence-based methods, protein function prediction and extended similarity group, make use of weakly similar sequences that are conventionally discarded in homology based function annotation. Combined together with experimental methods we hope that computational methods will make leading contribution in functional elucidation of the protein structures.     

Templated chemistry for bioorganic synthesis and chemical biology

In light of the 2018 Max Bergmann Medal, this review discusses advancements on chemical biology–driven templated chemistry developed in the author's laboratories. The focused review introduces the template categories applied to orient functional units such as functional groups, chromophores, biomolecules, or ligands in space. Unimolecular templates applied in protein synthesis facilitate fragment coupling of unprotected peptides. Templating via bimolecular assemblies provides control over proximity relationships between functional units of two molecules. As an instructive example, the coiled coil peptide–templated labelling of receptor proteins on live cells will be shown. Termolecular assemblies provide the opportunity to put the proximity of functional units on two (bio)molecules under the control of a third party molecule. This allows the design of conditional bimolecular reactions. A notable example is DNA/RNA–triggered peptide synthesis. The last section shows how termolecular and multimolecular assemblies can be used to better characterize and understand multivalent protein‐ligand interactions.     

CoBRA: Containerized Bioinformatics Workflow for Reproducible ChIP/ATAC-seq Analysis

Chromatin immunoprecipitation sequencing (ChIP-seq) and the Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) have become essential technologies to effectively measure protein–DNA interactions and chromatin accessibility. However, there is a need for a scalable and reproducible pipeline that incorporates proper normalization between samples, correction of copy number variations, and integration of new downstream analysis tools. Here we present Containerized Bioinformatics workflow for Reproducible ChIP/ATAC-seq Analysis (CoBRA), a modularized computational workflow which quantifies ChIP-seq and ATAC-seq peak regions and performs unsupervised and supervised analyses. CoBRA provides a comprehensive state-of-the-art ChIP-seq and ATAC-seq analysis pipeline that can be used by scientists with limited computational experience. This enables researchers to gain rapid insight into protein–DNA interactions and chromatin accessibility through sample clustering, differential peak calling, motif enrichment, comparison of sites to a reference database, and pathway analysis. CoBRA is publicly available online at https://bitbucket.org/cfce/cobra     

Computational analysis of membrane proteins: genomic occurrence, structure prediction and helix interactions

We review recent computational advances in the study of membrane proteins, focusing on those that have at least one transmembrane helix. Transmembrane protein regions are, in many respects, easier to investigate computationally than experimentally, due to the uniformity of their structure and interactions (e.g. consisting predominately of nearly parallel helices packed together) on one hand and presenting the challenges of solubility on the other. We present the progress made on identifying and classifying membrane proteins into families, predicting their structure from amino-acid sequence patterns (using many different methods), and analyzing their interactions and packing The total result of this work allows us for the first time to begin to think about the membrane protein interactome, the set of all interactions between distinct transmembrane helices in the lipid bilayer.     

Directional and quantitative phosphorylation networks.

Directionality in protein signalling networks is due to modulated protein-protein interactions and is fundamental for proper signal progression and response to external and internal cues. This property is in part enabled by linear motifs embedding post-translational modification sites. These serve as recognition sites, guiding phosphorylation by kinases and subsequent binding of modular domains (e.g. SH2 and BRCT). Characterization of such modification-modulated interactions on a proteome-wide scale requires extensive computational and experimental analysis. Here, we review the latest advances in methods for unravelling phosphorylation-mediated cellular interaction networks. In particular, we will discuss how the combination of new quantitative mass-spectrometric technologies and computational algorithms together are enhancing mapping of these largely uncharted dynamic networks. By combining quantitative measurements of phosphorylation events with computational approaches, we argue that systems level models will help to decipher complex diseases through the ability to predict cellular systems trajectories.     

Protein–DNA/RNA interactions: Machine intelligence tools and approaches in the era of artificial intelligence and big data

With the development of artificial intelligence (AI) technologies and the availability of large amounts of biological data, computational methods for proteomics have undergone a developmental process from traditional machine learning to deep learning. This review focuses on computational approaches and tools for the prediction of protein – DNA/RNA interactions using machine intelligence techniques. We provide an overview of the development progress of computational methods and summarize the advantages and shortcomings of these methods. We further compiled applications in tasks related to the protein – DNA/RNA interactions, and pointed out possible future application trends. Moreover, biological sequence-digitizing representation strategies used in different types of computational methods are also summarized and discussed.