Specificity of molecular interactions in transient protein-protein interaction interfaces

An integrative approach unveils a distal encounter site for rPTPε and phospho-Src complex formation

The structure determination of protein tyrosine phosphatase (PTP): phospho-protein complexes, which is essential to understand how specificity is achieved at the amino acid level, remains a significant challenge for protein crystallography and cryoEM due to the transient nature of binding interactions. Using rPTPεD1 and phospho-SrcKD as a model system, we have established an integrative workflow to address this problem, by means of which we generate a protein:phospho-protein complex model using predetermined protein structures, SAXS and pTyr-tailored MD simulations. Our model reveals transient protein-protein interactions between rPTPεD1 and phospho-SrcKD and is supported by three independent experimental validations. Measurements of the association rate between rPTPεD1 and phospho-SrcKD showed that mutations on the rPTPεD1: SrcKD complex interface disrupts these transient interactions, resulting in a reduction in protein-protein association rate and, eventually, phosphatase activity. This integrative approach is applicable to other PTP: phospho-protein complexes and the characterization of transient protein-protein interface interactions.

Using diverse potentials and scoring functions for the development of improved machine-learned models for protein-ligand affinity and docking pose prediction

The advent of computational drug discovery holds the promise of significantly reducing the effort of experimentalists, along with monetary cost. More generally, predicting the binding of small organic molecules to biological macromolecules has far-reaching implications for a range of problems, including metabolomics. However, problems such as predicting the bound structure of a protein-ligand complex along with its affinity have proven to be an enormous challenge. In recent years, machine learning-based methods have proven to be more accurate than older methods, many based on simple linear regression. Nonetheless, there remains room for improvement, as these methods are often trained on a small set of features, with a single functional form for any given physical effect, and often with little mention of the rationale behind choosing one functional form over another. Moreover, it is not entirely clear why one machine learning method is favored over another. In this work, we endeavor to undertake a comprehensive effort towards developing high-accuracy, machine-learned scoring functions, systematically investigating the effects of machine learning method and choice of features, and, when possible, providing insights into the relevant physics using methods that assess feature importance. Here, we show synergism among disparate features, yielding adjusted R² with experimental binding affinities of up to 0.871 on an independent test set and enrichment for native bound structures of up to 0.913. When purely physical terms that model enthalpic and entropic effects are used in the training, we use feature importance assessments to probe the relevant physics and hopefully guide future investigators working on this and other computational chemistry problems.

Soft regions of protein surface are potent for stable dimer formation

By having knowledge about the characteristics of protein interaction interfaces, we will be able to manipulate protein complexes for therapies. Dimer state is considered as the primary alphabet of the most proteins' quaternary structure. The properties of binding interface between subunits and of noninterface region define the specificity and stability of the intended protein complex. Considering some topological properties and amino acids' affinity for binding in interfaces of protein dimers, we construct the interface-specific recurrence plots. The data obtained from recurrence quantitative analysis, and accessibility-related metrics help us to classify the protein dimers into four distinct classes. Some mechanical properties of binding interfaces are computed for each predefined class of the dimers. The computed mechanical characteristics of binding patch region are compared with those of nonbinding region of proteins. Our observations indicate that the mechanical properties of protein binding sites have a decisive impact on determining the dimer stability. We introduce a new concept in analyzing protein structure by considering mechanical properties of protein structure. We conclude that the interface region between subunits of stable dimers is usually mechanically softer than the interface of unstable protein dimers. AbbreviationsAABaverage affinity for bindingANManisotropic network modelAPCaffinity propagation clusteringASAaccessible surface areaCCDinter residues distanceCSCcomplex stability codeDMdistance matrixΔG_dissPISA-computed dissociation free energyGNMGaussian normal mode analysisNMAnormal mode analysisPBPprotein binding patchPISAproteins, interfaces, structures and assembliesrASArelative accessible area in respect to unfolded state of residuesRMrecurrence matrixrPrelative protrusionRPrecurrence plotRQArecurrence quantitative analysisSEMstandard error of meanCommunicated by Ramaswamy H. Sarma.

Exploring the potential of 3D Zernike descriptors and SVM for protein-protein interface prediction

Background

The correct determination of protein–protein interaction interfaces is important for understanding disease mechanisms and for rational drug design. To date, several computational methods for the prediction of protein interfaces have been developed, but the interface prediction problem is still not fully understood. Experimental evidence suggests that the location of binding sites is imprinted in the protein structure, but there are major differences among the interfaces of the various protein types: the characterising properties can vary a lot depending on the interaction type and function. The selection of an optimal set of features characterising the protein interface and the development of an effective method to represent and capture the complex protein recognition patterns are of paramount importance for this task.

Results

In this work we investigate the potential of a novel local surface descriptor based on 3D Zernike moments for the interface prediction task. Descriptors invariant to roto-translations are extracted from circular patches of the protein surface enriched with physico-chemical properties from the HQI8 amino acid index set, and are used as samples for a binary classification problem. Support Vector Machines are used as a classifier to distinguish interface local surface patches from non-interface ones. The proposed method was validated on 16 classes of proteins extracted from the Protein–Protein Docking Benchmark 5.0 and compared to other state-of-the-art protein interface predictors (SPPIDER, PrISE and NPS-HomPPI).

Conclusions

The 3D Zernike descriptors are able to capture the similarity among patterns of physico-chemical and biochemical properties mapped on the protein surface arising from the various spatial arrangements of the underlying residues, and their usage can be easily extended to other sets of amino acid properties. The results suggest that the choice of a proper set of features characterising the protein interface is crucial for the interface prediction task, and that optimality strongly depends on the class of proteins whose interface we want to characterise. We postulate that different protein classes should be treated separately and that it is necessary to identify an optimal set of features for each protein class.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2043-3) contains supplementary material, which is available to authorized users.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

What Glues a Homodimer Together: Systematic Analysis of the Stabilizing Effect of an Aromatic Hot Spot in the Protein-Protein Interface of the tRNA-Modifying Enzyme Tgt

Shigella bacteria constitute the causative agent of bacillary dysentery, an acute inflammatory disease causing the death of more than one million humans per year. A null mutation in the tgt gene encoding the tRNA-modifying enzyme tRNA-guanine transglycosylase (Tgt) was found to drastically decrease the pathogenicity of Shigella bacteria, suggesting the use of Tgt as putative target for selective antibiotics. The enzyme is only functionally active as a homodimer; thus, interference with the formation of its protein-protein interface is an attractive opportunity for therapeutic intervention. To better understand the driving forces responsible for the assembly, stability, and formation of the homodimer, we studied the properties of the residues that establish the dimer interface in detail. We performed site-directed mutagenesis and controlled shifts in the monomer/dimer equilibrium ratio in solution in a concentration-dependent manner by native mass spectrometry and used crystal structure analysis to elucidate the geometrical modulations resulting from mutational variations. The wild-type enzyme exhibits nearly exclusive dimer geometry. A patch of four aromatic amino acids, embedded into a ring of hydrophobic residues and further stabilized by a network of H-bonds, is essential for the stability of the dimer's contact. Accordingly, any perturbance in the constitution of this aromatic patch by nonaromatic residues reduces dimer stability significantly, with some of these exchanges resulting in a nearly exclusively monomeric state. Apart from the aromatic hot spot, the interface comprises an extended loop-helix motif that exhibits remarkable flexibility. In the destabilized mutated variants, the loop-helix motif adopts deviating conformations in the interface region, and a number of water molecules, penetrating into the interface, are observed.

Affinity- and specificity-enhancing mutations are frequent in multispecific interactions between TIMP2 and MMPs

Multispecific proteins play a major role in controlling various functions such as signaling, regulation of transcription/translation, and immune response. Hence, a thorough understanding of the atomic-level principles governing multispecific interactions is important not only for the advancement of basic science but also for applied research such as drug design. Here, we study evolution of an exemplary multispecific protein, a Tissue Inhibitor of Matrix Metalloproteinases 2 (TIMP2) that binds with comparable affinities to more than twenty-six members of the Matrix Metalloproteinase (MMP) and the related ADAMs families. We postulate that due to its multispecific nature, TIMP2 is not optimized to bind to any individual MMP type, but rather embodies a compromise required for interactions with all MMPs. To explore this hypothesis, we perform computational saturation mutagenesis of the TIMP2 binding interface and predict changes in free energy of binding to eight MMP targets. Computational results reveal the non-optimality of the TIMP2 binding interface for all studied proteins, identifying many affinity-enhancing mutations at multiple positions. Several TIMP2 point mutants predicted to enhance binding affinity and/or binding specificity towards MMP14 were selected for experimental verification. Experimental results show high abundance of affinity-enhancing mutations in TIMP2, with some point mutations producing more than ten-fold improvement in affinity to MMP14. Our computational and experimental results collaboratively demonstrate that the TIMP2 sequence lies far from the fitness maximum when interacting with its target enzymes. This non-optimality of the binding interface and high potential for improvement might characterize all proteins evolved for binding to multiple targets.

Improving protein complex classification accuracy using amino acid composition profile

Protein complex prediction approaches are based on the assumptions that complexes have dense protein-protein interactions and high functional similarity between their subunits. We investigated those assumptions by studying the subunits' interaction topology, sequence similarity and molecular function for human and yeast protein complexes. Inclusion of amino acids' physicochemical properties can provide better understanding of protein complex properties. Principal component analysis is carried out to determine the major features. Adopting amino acid composition profile information with the SVM classifier serves as an effective post-processing step for complexes classification. Improvement is based on primary sequence information only, which is easy to obtain.

Large Scale Chemical Cross-linking Mass Spectrometry Perspectives

The spectacular heterogeneity of a complex protein mixture from biological samples becomes even more difficult to tackle when one’s attention is shifted towards different protein complex topologies, transient interactions, or localization of PPIs. Meticulous protein-by-protein affinity pull-downs and yeast-two-hybrid screens are the two approaches currently used to decipher proteome-wide interaction networks. Another method is to employ chemical cross-linking, which gives not only identities of interactors, but could also provide information on the sites of interactions and interaction interfaces. Despite significant advances in mass spectrometry instrumentation over the last decade, mapping Protein-Protein Interactions (PPIs) using chemical cross-linking remains time consuming and requires substantial expertise, even in the simplest of systems. While robust methodologies and software exist for the analysis of binary PPIs and also for the single protein structure refinement using cross-linking-derived constraints, undertaking a proteome-wide cross-linking study is highly complex. Difficulties include i) identifying cross-linkers of the right length and selectivity that could capture interactions of interest; ii) enrichment of the cross-linked species; iii) identification and validation of the cross-linked peptides and cross-linked sites.

In this review we examine existing literature aimed at the large-scale protein cross-linking and discuss possible paths for improvement. We also discuss short-length cross-linkers of broad specificity such as formaldehyde and diazirine-based photo-cross-linkers. These cross-linkers could potentially capture many types of interactions, without strict requirement for a particular amino-acid to be present at a given protein-protein interface. How these shortlength, broad specificity cross-linkers be applied to proteome-wide studies? We will suggest specific advances in methodology, instrumentation and software that are needed to make such a leap.

Counterbalance of ligand- and self-coupled motions characterizes multispecificity of ubiquitin

Date hub proteins are a type of proteins that show multispecificity in a time-dependent manner. To understand dynamic aspects of such multispecificity we studied Ubiquitin as a typical example of a date hub protein. Here we analyzed 9 biologically relevant Ubiquitin-protein (ligand) heterodimer structures by using normal mode analysis based on an elastic network model. Our result showed that the self-coupled motion of Ubiquitin in the complex, rather than its ligand-coupled motion, is similar to the motion of Ubiquitin in the unbound condition. The ligand-coupled motions are correlated to the conformational change between the unbound and bound conditions of Ubiquitin. Moreover, ligand-coupled motions favor the formation of the bound states, due to its in-phase movements of the contacting atoms at the interface. The self-coupled motions at the interface indicated loss of conformational entropy due to binding. Therefore, such motions disfavor the formation of the bound state. We observed that the ligand-coupled motions are embedded in the motions of unbound Ubiquitin. In conclusion, multispecificity of Ubiquitin can be characterized by an intricate balance of the ligand- and self-coupled motions, both of which are embedded in the motions of the unbound form.

Structural and functional analysis of multi-interface domains

A multi-interface domain is a domain that can shape multiple and distinctive binding sites to contact with many other domains, forming a hub in domain-domain interaction networks. The functions played by the multiple interfaces are usually different, but there is no strict bijection between the functions and interfaces as some subsets of the interfaces play the same function. This work applies graph theory and algorithms to discover fingerprints for the multiple interfaces of a domain and to establish associations between the interfaces and functions, based on a huge set of multi-interface proteins from PDB. We found that about 40% of proteins have the multi-interface property, however the involved multi-interface domains account for only a tiny fraction (1.8%) of the total number of domains. The interfaces of these domains are distinguishable in terms of their fingerprints, indicating the functional specificity of the multiple interfaces in a domain. Furthermore, we observed that both cooperative and distinctive structural patterns, which will be useful for protein engineering, exist in the multiple interfaces of a domain.

Computer-aided antibody design

Recent clinical trials using antibodies with low toxicity and high efficiency have raised expectations for the development of next-generation protein therapeutics. However, the process of obtaining therapeutic antibodies remains time consuming and empirical. This review summarizes recent progresses in the field of computer-aided antibody development mainly focusing on antibody modeling, which is divided essentially into two parts: (i) modeling the antigen-binding site, also called the complementarity determining regions (CDRs), and (ii) predicting the relative orientations of the variable heavy (V(H)) and light (V(L)) chains. Among the six CDR loops, the greatest challenge is predicting the conformation of CDR-H3, which is the most important in antigen recognition. Further computational methods could be used in drug development based on crystal structures or homology models, including antibody-antigen dockings and energy calculations with approximate potential functions. These methods should guide experimental studies to improve the affinities and physicochemical properties of antibodies. Finally, several successful examples of in silico structure-based antibody designs are reviewed. We also briefly review structure-based antigen or immunogen design, with application to rational vaccine development.

A novel method for protein-protein interaction site prediction using phylogenetic substitution models

Protein-protein binding events mediate many critical biological functions in the cell. Typically, functionally important sites in proteins can be well identified by considering sequence conservation. However, protein-protein interaction sites exhibit higher sequence variation than other functional regions, such as catalytic sites of enzymes. Consequently, the mutational behavior leading to weak sequence conservation poses significant challenges to the protein-protein interaction site prediction. Here, we present a phylogenetic framework to capture critical sequence variations that favor the selection of residues essential for protein-protein binding. Through the comprehensive analysis of diverse protein families, we show that protein binding interfaces exhibit distinct amino acid substitution as compared with other surface residues. On the basis of this analysis, we have developed a novel method, BindML, which utilizes the substitution models to predict protein-protein binding sites of protein with unknown interacting partners. BindML estimates the likelihood that a phylogenetic tree of a local surface region in a query protein structure follows the substitution patterns of protein binding interface and nonbinding surfaces. BindML is shown to perform well compared to alternative methods for protein binding interface prediction. The methodology developed in this study is very versatile in the sense that it can be generally applied for predicting other types of functional sites, such as DNA, RNA, and membrane binding sites in proteins.

Macromolecular complexes in crystals and solutions

This paper presents a discussion of existing methods for the analysis of macromolecular interactions and complexes in crystal packing. Typical situations and conditions where wrong answers may be obtained in the course of ordinary procedures are presented and discussed. The more general question of what the relationship is between natural (in-solvent) and crystallized assemblies is discussed and researched. A computational analysis suggests that weak interactions with K(d) ≥ 100 µM have a considerable chance of being lost during the course of crystallization. In such instances, crystal packing misrepresents macromolecular complexes and interactions. For as many as 20% of protein dimers in the PDB the likelihood of misrepresentation is estimated to be higher than 50%. Given that weak macromolecular interactions play an important role in many biochemical processes, these results suggest that a complementary noncrystallographic study should be always conducted when inferring structural aspects of weakly bound complexes.

ReplicOpter: a replicate optimizer for flexible docking

We present a computationally efficient method for flexible refinement of docking predictions that reflects observed motions within a protein's structural class. Using structural homologs, we derive deformation models that capture likely motions. The models or "replicates" typically align along a rigid core, with a handful of flexible loops, linkers and tails. A few replicates can generate a much larger number of conformers, by exchanging each flexible region independently of the others. In this way, 10 replicates of a protein having 6 flexible regions can be used to generate a million conformations of a molecule. While this has obvious advantages in terms of sampling, the cost of assessing energies at every conformer is prohibitive, particularly when both molecules are flexible. Our approach addresses this combinatorial explosion, using key assumptions to compress the sampling by many orders of magnitude. ReplicOpter can perform hierarchical clustering from a list of rigid docking predictions and find nearby structures to any promising cluster representatives. These predicted complexes can then be refined and rescored. ReplicOpter's scoring function includes a Lennard-Jones potential softened using the Anderson-Chandler-Weeks decomposition, a desolvation term derived from the Atomic Contact Energy function, Coulombic electrostatics, hydrogen bonding, and terms to model pi-pi and pi-cation interactions. ReplicOpter has performed well on several recent CAPRI systems. We are presently benchmarking ReplicOpter on the complete docking benchmark set to fully establish its utility in refining rigid docking predictions and identifying near-native solutions.

Stacking and energetic contribution of aromatic islands at the binding interface of antibody proteins

Background

The enrichment and importance of some aromatic residues, such as Tyr and Trp, have been widely noticed at the binding interfaces of antibodies from many experimental and statistical results, some of which were even identified as “hot spots” contributing significantly greater to the binding affinity than other amino acids. However, how these aromatic residues influence the immune binding still deserves further investigation. A large-scale examination was done regarding the local spatial environment around the interfacial Tyr or Trp residues. Energetic contribution of these Tyr and Trp residues to the binding affinity was then studied regarding 82 representative antibody interfaces covering 509 immune complexes from the PDB database and IMGT/3Dstructure-DB.

Results

The connectivity analysis of interfacial residues showed that Tyr and Trp tended to cluster into the spatial Aromatic Islands (AI) rather than being distributed randomly at the antibody interfaces. Out of 82 antibody-antigen complexes, 72% (59) interfaces were found to contain AI with more than 3 aromatic residues. The statistical test against an empirical distribution indicated that the existence of AI was significant in about 60% representative antibody interfaces. Secondly, the loss of solvent accessible surface area (SASA) for side chains of aromatic residues between actually crowded state and independent state was nicely correlated with the AI size increasing in a linearly positive way which indicated that the aromatic side chains in AI tended to take a compact and ordered stacking conformation at the interfaces. Interestingly, the SASA loss of AI was also correlated roughly with the averaged gap of binding free energy between the theoretical and experimental data for immune complexes.

Conclusions

The results of our study revealed the wide existence and statistical significance of “Aromatic Island” (AI) composed of the spatially clustered Tyr and Trp residues at the antibody interfaces. The regular arrangement and stacking of aromatic side chains in AI could probably produce extra cooperative effects to the binding affinity which was firstly observed through the large-scale data analysis. The finding in this work not only provides insights into the functional role of aromatic residues in the antibody-antigen interaction, but also may facilitate the antibody engineering and potential clinical applications.

Knitting and untying the protein network: modulation of protein ensembles as a therapeutic strategy

Proteins constitute the working machinery and structural support of all organisms. In performing a given function, they must adopt highly specific structures that can change with their level of activity, often through the direct or indirect action of other proteins. Indeed, proteins typically function within an ensemble, rather than individually. Hence, they must be sufficiently flexible to interact with each other and execute diverse tasks. The discovery that errors within these groups can ultimately cause disease has led to a paradigm shift in drug discovery, from an emphasis on single protein targets to a holistic approach whereby entire ensembles are targeted.

Geometric characteristics of hydrogen bonds involving sulfur atoms in proteins

Sulfur atoms have been known to participate in hydrogen bonds (H-bonds) and these sulfur-containing H-bonds (SCHBs) are suggested to play important roles in certain biological processes. This study aims to comprehensively characterize all the SCHBs in 500 high-resolution protein structures (< or =1.8 A). We categorized SCHBs into six types according to donor/acceptor behaviors and used explicit hydrogen approach to distinguish SCHBs from those of nonhydrogen bonding interactions. It is revealed that sulfur atom is a very poor H-bond acceptor, but a moderately good H-bond donor. In alpha-helix, considerable SCHBs were found between the sulphydryl group of cysteine residue i and the carbonyl oxygen of residue i-4, and these SCHBs exert effects in stabilizing helices. Although for other SCHBs, they possess no specific secondary structural preference, their geometric characteristics in proteins and in free small compounds are significantly distinct, indicating the protein SCHBs are geometrically distorted. Interestingly, sulfur atom in the disulfide bond tends to form bifurcated H-bond whereas in cysteine-cysteine pairs prefer to form dual H-bond. These special H-bonds remarkably boost the interaction between H-bond donor and acceptor. By oxidation/reduction manner, the mutual transformation between the dual H-bonds and disulfide bonds for cysteine-cysteine pairs can accurately adjust the structural stability and biological function of proteins in different environments. Furthermore, few loose H-bonds were observed to form between the sulphydryl groups and aromatic rings, and in these cases the donor H is almost over against the rim rather than the center of the aromatic ring.

2D depiction of nonbonding interactions for protein complexes

A program called the 2D-GraLab is described for automatically generating schematic representation of nonbonding interactions across the protein binding interfaces. The input file of this program takes the standard PDB format, and the outputs are two-dimensional PostScript diagrams giving intuitive and informative description of the protein-protein interactions and their energetics properties, including hydrogen bond, salt bridge, van der Waals interaction, hydrophobic contact, pi-pi stacking, disulfide bond, desolvation effect, and loss of conformational entropy. To ensure these interaction information are determined accurately and reliably, methods and standalone programs employed in the 2D-GraLab are all widely used in the chemistry and biology community. The generated diagrams allow intuitive visualization of the interaction mode and binding specificity between two subunits in protein complexes, and by providing information on nonbonding energetics and geometric characteristics, the program offers the possibility of comparing different protein binding profiles in a detailed, objective, and quantitative manner. We expect that this 2D molecular graphics tool could be useful for the experimentalists and theoreticians interested in protein structure and protein engineering.

A feature-based approach to modeling protein-protein interaction hot spots

Identifying features that effectively represent the energetic contribution of an individual interface residue to the interactions between proteins remains problematic. Here, we present several new features and show that they are more effective than conventional features. By combining the proposed features with conventional features, we develop a predictive model for interaction hot spots. Initially, 54 multifaceted features, composed of different levels of information including structure, sequence and molecular interaction information, are quantified. Then, to identify the best subset of features for predicting hot spots, feature selection is performed using a decision tree. Based on the selected features, a predictive model for hot spots is created using support vector machine (SVM) and tested on an independent test set. Our model shows better overall predictive accuracy than previous methods such as the alanine scanning methods Robetta and FOLDEF, and the knowledge-based method KFC. Subsequent analysis yields several findings about hot spots. As expected, hot spots have a larger relative surface area burial and are more hydrophobic than other residues. Unexpectedly, however, residue conservation displays a rather complicated tendency depending on the types of protein complexes, indicating that this feature is not good for identifying hot spots. Of the selected features, the weighted atomic packing density, relative surface area burial and weighted hydrophobicity are the top 3, with the weighted atomic packing density proving to be the most effective feature for predicting hot spots. Notably, we find that hot spots are closely related to π–related interactions, especially π · · · π interactions.

Computational structural analysis: multiple proteins bound to DNA

Background

With increasing numbers of crystal structures of protein∶DNA and protein∶protein∶DNA complexes publically available, it is now possible to extract sufficient structural, physical-chemical and thermodynamic parameters to make general observations and predictions about their interactions. In particular, the properties of macromolecular assemblies of multiple proteins bound to DNA have not previously been investigated in detail.

Methodology/Principal Findings

We have performed computational structural analyses on macromolecular assemblies of multiple proteins bound to DNA using a variety of different computational tools: PISA; PROMOTIF; X3DNA; ReadOut; DDNA and DCOMPLEX. Additionally, we have developed and employed an algorithm for approximate collision detection and overlapping volume estimation of two macromolecules. An implementation of this algorithm is available at http://promoterplot.fmi.ch/Collision1/. The results obtained are compared with structural, physical-chemical and thermodynamic parameters from protein∶protein and single protein∶DNA complexes. Many of interface properties of multiple protein∶DNA complexes were found to be very similar to those observed in binary protein∶DNA and protein∶protein complexes. However, the conformational change of the DNA upon protein binding is significantly higher when multiple proteins bind to it than is observed when single proteins bind. The water mediated contacts are less important (found in less quantity) between the interfaces of components in ternary (protein∶protein∶DNA) complexes than in those of binary complexes (protein∶protein and protein∶DNA).The thermodynamic stability of ternary complexes is also higher than in the binary interactions. Greater specificity and affinity of multiple proteins binding to DNA in comparison with binary protein-DNA interactions were observed. However, protein-protein binding affinities are stronger in complexes without the presence of DNA.

Conclusions/Significance

Our results indicate that the interface properties: interface area; number of interface residues/atoms and hydrogen bonds; and the distribution of interface residues, hydrogen bonds, van der Walls contacts and secondary structure motifs are independent of whether or not a protein is in a binary or ternary complex with DNA. However, changes in the shape of the DNA reduce the off-rate of the proteins which greatly enhances the stability and specificity of ternary complexes compared to binary ones.

Characterization of protein-protein interfaces

We analyze the characteristics of protein-protein interfaces using the largest datasets available from the Protein Data Bank (PDB). We start with a comparison of interfaces with protein cores and non-interface surfaces. The results show that interfaces differ from protein cores and non-interface surfaces in residue composition, sequence entropy, and secondary structure. Since interfaces, protein cores, and non-interface surfaces have different solvent accessibilities, it is important to investigate whether the observed differences are due to the differences in solvent accessibility or differences in functionality. We separate out the effect of solvent accessibility by comparing interfaces with a set of residues having the same solvent accessibility as the interfaces. This strategy reveals residue distribution propensities that are not observable by comparing interfaces with protein cores and non-interface surfaces. Our conclusions are that there are larger numbers of hydrophobic residues, particularly aromatic residues, in interfaces, and the interactions apparently favored in interfaces include the opposite charge pairs and hydrophobic pairs. Surprisingly, Pro-Trp pairs are over represented in interfaces, presumably because of favorable geometries. The analysis is repeated using three datasets having different constraints on sequence similarity and structure quality. Consistent results are obtained across these datasets. We have also investigated separately the characteristics of heteromeric interfaces and homomeric interfaces.