Prediction of Protein-Protein Binding Affinities from Unbound Protein Structures.

Proteins are the workhorses of cells to carry out sophisticated and complex cellular processes. Such processes require a coordinated and regulated interactions between proteins that are both time and location specific. The strength, or binding affinity, of protein-protein interactions ranges between the micro- and the nanomolar association constant, often dictating the molecular mechanisms underlying the interaction and the longevity of the complex, i.e., transient or permanent. In consequence, there is a need to quantify the strength of protein-protein interactions for biological, biomedical, and biotechnological applications. While experimental methods are labor intensive and costly, computational ones are useful tools to predict the affinity of protein-protein interactions. In this chapter, we review the methods developed by us to address this question. We briefly present two methods to comprehend the structure of the protein complex derived by either comparative modeling or docking. Then we introduce BADOCK, a method to predict the binding energy without requiring the structure of the protein complex, thus overcoming one of the major limitations of structure-based methods for the prediction of binding affinity. BADOCK utilizes the structure of unbound proteins and the protein docking sampling space to predict protein-protein binding affinities. We present step-by-step protocols to utilize these methods, describing the inputs and potential pitfalls as well as their respective strengths and limitations.

SPServer: split-statistical potentials for the analysis of protein structures and protein-protein interactions.

BACKGROUND: Statistical potentials, also named knowledge-based potentials, are scoring functions derived from empirical data that can be used to evaluate the quality of protein folds and protein-protein interaction (PPI) structures. In previous works we decomposed the statistical potentials in different terms, named Split-Statistical Potentials, accounting for the type of amino acid pairs, their hydrophobicity, solvent accessibility and type of secondary structure. These potentials have been successfully used to identify near-native structures in protein structure prediction, rank protein docking poses, and predict PPI binding affinities. RESULTS: Here, we present the SPServer, a web server that applies the Split-Statistical Potentials to analyze protein folds and protein interfaces. SPServer provides global scores as well as residue/residue-pair profiles presented as score plots and maps. This level of detail allows users to: (1) identify potentially problematic regions on protein structures; (2) identify disrupting amino acid pairs in protein interfaces; and (3) compare and analyze the quality of tertiary and quaternary structural models. CONCLUSIONS: While there are many web servers that provide scoring functions to assess the quality of either protein folds or PPI structures, SPServer integrates both aspects in a unique easy-to-use web server. Moreover, the server permits to locally assess the quality of the structures and interfaces at a residue level and provides tools to compare the local assessment between structures. SERVER ADDRESS: https://sbi.upf.edu/spserver/ .

Galaxy InteractoMIX: An Integrated Computational Platform for the Study of Protein-Protein Interaction Data.

Protein interactions play a crucial role among the different functions of a cell and are central to our understanding of cellular processes both in health and disease. Here we present Galaxy InteractoMIX (http://galaxy.interactomix.com), a platform composed of 13 different computational tools each addressing specific aspects of the study of protein-protein interactions, ranging from large-scale cross-species protein-wide interactomes to atomic resolution level of protein complexes. Galaxy InteractoMIX provides an intuitive interface where users can retrieve consolidated interactomics data distributed across several databases or uncover links between diseases and genes by analyzing the interactomes underlying these diseases. The platform makes possible large-scale prediction and curation protein interactions using the conservation of motifs, interology, or presence or absence of key sequence signatures. The range of structure-based tools includes modeling and analysis of protein complexes, delineation of interfaces and the modeling of peptides acting as inhibitors of protein-protein interactions. Galaxy InteractoMIX includes a range of ready-to-use workflows to run complex analyses requiring minimal intervention by users. The potential range of applications of the platform covers different aspects of life science, biomedicine, biotechnology and drug discovery where protein associations are studied.

Using collections of structural models to predict changes of binding affinity caused by mutations in protein-protein interactions.

Protein-protein interactions (PPIs) in all the molecular aspects that take place both inside and outside cells. However, determining experimentally the structure and affinity of PPIs is expensive and time consuming. Therefore, the development of computational tools, as a complement to experimental methods, is fundamental. Here, we present a computational suite: MODPIN, to model and predict the changes of binding affinity of PPIs. In this approach we use homology modeling to derive the structures of PPIs and score them using state-of-the-art scoring functions. We explore the conformational space of PPIs by generating not a single structural model but a collection of structural models with different conformations based on several templates. We apply the approach to predict the changes in free energy upon mutations and splicing variants of large datasets of PPIs to statistically quantify the quality and accuracy of the predictions. As an example, we use MODPIN to study the effect of mutations in the interaction between colicin endonuclease 9 and colicin endonuclease 2 immune protein from Escherichia coli. Finally, we have compared our results with other state-of-art methods.

On the prediction of DNA-binding preferences of C2H2-ZF domains using structural models: application on human CTCF.

Cis2-His2 zinc finger (C2H2-ZF) proteins are the largest family of transcription factors in human and higher metazoans. To date, the DNA-binding preferences of many members of this family remain unknown. We have developed a computational method to predict their DNA-binding preferences. We have computed theoretical position weight matrices (PWMs) of proteins composed by C2H2-ZF domains, with the only requirement of an input structure. We have predicted more than two-third of a single zinc-finger domain binding site for about 70% variants of Zif268, a classical member of this family. We have successfully matched between 60 and 90% of the binding-site motif of examples of proteins composed by three C2H2-ZF domains in JASPAR, a standard database of PWMs. The tests are used as a proof of the capacity to scan a DNA fragment and find the potential binding sites of transcription-factors formed by C2H2-ZF domains. As an example, we have tested the approach to predict the DNA-binding preferences of the human chromatin binding factor CTCF. We offer a server to model the structure of a zinc-finger protein and predict its PWM.

Zip6 Transporter Is an Essential Component of the Lymphocyte Activation Machinery.

Zinc deficiency causes immune dysfunction. In T lymphocytes, hypozincemia promotes thymus atrophy, polarization imbalance, and altered cytokine production. Zinc supplementation is commonly used to boost immune function to prevent infectious diseases in at-risk populations. However, the molecular players involved in zinc homeostasis in lymphocytes are poorly understood. In this paper, we wanted to determine the identity of the transporter responsible for zinc entry into lymphocytes. First, in human Jurkat cells, we characterized the effect of zinc on proliferation and activation and found that zinc supplementation enhances activation when T lymphocytes are stimulated using anti-CD3/anti-CD28 Abs. We show that zinc entry depends on specific pathways to correctly tune the NFAT, NF-κB, and AP-1 activation cascades. Second, we used various human and murine models to characterize the zinc transporter family, Zip, during T cell activation and found that Zip6 was strongly upregulated early during activation. Therefore, we generated a Jurkat Zip6 knockout (KO) line to study how the absence of this transporter affects lymphocyte physiology. We found that although Zip6KO cells showed no altered zinc transport or proliferation under basal conditions, under activation, these KO cells showed deficient zinc transport and a drastically impaired activation program. Our work shows that zinc entry into activated lymphocytes depends on Zip6 and that this transporter is essential for the correct function of the cellular activation machinery.