- Research article
- Open Access
Comparison of tertiary structures of proteins in protein-protein complexes with unbound forms suggests prevalence of allostery in signalling proteins
© Swapna et al.; licensee BioMed Central Ltd. 2012
- Received: 27 November 2011
- Accepted: 5 April 2012
- Published: 3 May 2012
Most signalling and regulatory proteins participate in transient protein-protein interactions during biological processes. They usually serve as key regulators of various cellular processes and are often stable in both protein-bound and unbound forms. Availability of high-resolution structures of their unbound and bound forms provides an opportunity to understand the molecular mechanisms involved. In this work, we have addressed the question “What is the nature, extent, location and functional significance of structural changes which are associated with formation of protein-protein complexes?”
A database of 76 non-redundant sets of high resolution 3-D structures of protein-protein complexes, representing diverse functions, and corresponding unbound forms, has been used in this analysis. Structural changes associated with protein-protein complexation have been investigated using structural measures and Protein Blocks description. Our study highlights that significant structural rearrangement occurs on binding at the interface as well as at regions away from the interface to form a highly specific, stable and functional complex. Notably, predominantly unaltered interfaces interact mainly with interfaces undergoing substantial structural alterations, revealing the presence of at least one structural regulatory component in every complex.
Interestingly, about one-half of the number of complexes, comprising largely of signalling proteins, show substantial localized structural change at surfaces away from the interface. Normal mode analysis and available information on functions on some of these complexes suggests that many of these changes are allosteric. This change is largely manifest in the proteins whose interfaces are altered upon binding, implicating structural change as the possible trigger of allosteric effect. Although large-scale studies of allostery induced by small-molecule effectors are available in literature, this is, to our knowledge, the first study indicating the prevalence of allostery induced by protein effectors.
The enrichment of allosteric sites in signalling proteins, whose mutations commonly lead to diseases such as cancer, provides support for the usage of allosteric modulators in combating these diseases.
- Protein Data Bank
- Root Mean Square Deviation
- Normal Mode Analysis
- Interface Residue
- Wilcoxon Paired Test
Protein-protein interactions participate in myriad processes of the cell such as replication, transcription, translation, signal transduction, immune response, metabolism, membrane-associated processes and development (e.g., [1–4]). Protein-protein interactions offer an excellent way of combining its limited working parts, the proteins, to achieve large functional diversity using a limited genetic repertoire . Abnormal interactions between proteins within the cell or from pathogens cause many human diseases . Protein binding can also elicit an allosteric response. Allostery is an integral and pervasive mechanism employed by nature to modulate cellular processes [7–11]. It serves as a key mechanism for obtaining fine-tuned regulation in several cellular processes – from metabolic pathways, signalling systems  to gene regulation . Functional modulation is achieved either by enhancing (positive co-operativity) or decreasing (negative co-operativity) levels of function. The effect at target site can be varied, e.g., activation of catalysis, regulation of ligand-binding, control of complex formation.
Given their importance, several high-throughput interaction assays [14, 15], such as yeast two-hybrid and tandem affinity purification, have been developed to supplement the dataset of protein-protein interactions from low-throughput methods [16, 17]. However, such large-scale experimental methods suffer from high false-positive rates . The gold standard for protein-protein interactions is usually a dataset of complexes of interacting proteins solved using X-ray crystallography [19–21]. Although it is a much smaller and incomplete dataset in comparison to high-throughput protein-protein interaction datasets, it is reliable and enables mapping of interaction regions and structural changes which accompany interactions. Several derived databases provide protein-protein interaction datasets in various easy to –study and –use formats. SCOPPI , iPfam , SNAPPI-DB , 3D Complex , InterEvol  and ProtCID  are some of the available 3D structural databases of protein-protein complexes.
Protein-protein interactions can be classified into different kinds : homo-oligomers and hetero-oligomers; obligate and non-obligate complexes; permanent and transient complexes. Non-obligate complexes form an important class since they serve as key regulators in maintaining and regulating cellular homeostasis [29–31]. They are also valuable from the viewpoint of structural biology since both the unbound and bound forms can be crystallized owing to their stability. Several such structures have been solved by various groups and deposited in the Protein Data Bank (PDB) . An invaluable non-redundant dataset of structures of both the interacting partners solved in unbound and bound form has been collated, curated and updated by Weng and colleagues [33, 34]. The ComSin database provides a unique collection of structures of proteins solved in unbound and bound form, targeted towards disorder–order transitions .
Earlier studies of structures of protein-protein complexes using both the unbound and bound form of proteins reveal that proteins undergo changes in their structure upon binding. Betts and Sternberg  were the first to compare the bound and unbound forms using a dataset of 39 complexes. Martin et al. [37, 38] analyzed a dataset of 83 complexes in terms of local structural variations. The alterations in structure as a result of protein-protein interactions manifest either as a rigid-body shift of a segment or as a conformational change from one secondary structural form to another . The extent of conformational change observed at the interface upon binding prompted several studies to understand and predict these changes [40–42]. Such studies aim to improve protein-protein docking methods  and help in the accurate docking of protein-protein interactions, which can be used to understand the mechanism of functioning of the complex or design inhibitors.
In this work, we have used a curated and non-redundant dataset of 76 protein-protein complexes, solved using X-ray crystallography in high resolution in both unbound and bound form, to address questions about the nature, extent and location of structural changes upon binding. We noticed that, in addition to changes in the interface, possibly allosteric changes causing structural alteration occur in about half of the complexes, indicating a much higher prevalence of this phenomenon caused due to protein binding than appreciated before.
Proteins bound to other proteins undergo larger structural changes than unliganded proteins
Structural change observed in different forms of a protein could be due to experimental artifacts , intrinsic flexibility  or due to a biologically important external perturbation , such as ligand binding or post-translational modification. To differentiate structural changes potentially related to protein-protein interactions from those which are artefacts, we compared variations occurring in the dataset of protein-protein complexes with two control datasets (see Additional file 1: Table S1). The first control set (named Control – Rigid) consists of 50 structures, solved at a resolution ≤2.5 Å, of two fairly rigid and extensively studied proteins: bovine ribonuclease A and sperm whale myoglobin, and provides an indicator of co-ordinate uncertainties. The second control set (named Control – Monomer) consists of a non-redundant set of 95 clusters of structures of monomers, also solved at a resolution ≤2.5 Å, which serve as a heterogeneous set since this dataset contains both rigid and flexible proteins, thus serving as a control set for understanding intrinsic flexibility. The main dataset of our study named PPC (Protein-protein complexes) is an extensively curated dataset of non-obligatory proteins with their 3-D structures solved in both unbound and bound forms (Additional file 2: Table S2). It consists of 76 non-obligatory complexes representing members of diverse functions (25 enzyme-inhibitor, 11 antigen-antibody and 40 ‘other’ complexes, which largely comprises of signalling proteins). The number of proteins involved in the 76 complexes represent the major SCOP (Structural Classification of Proteins)  classes (all α - 32, all β - 84, α/β - 57, α + β - 37). Since the dataset has been pruned to exclude cases with a large percentage of missing residues at the interface, disordered proteins are under-represented. The complexes predominantly involve two-chain interactions and some interactions involving three chains, in which two chains are considered as a single entity (for example, in the case of light and heavy chains of antibody and in the case of Gβ-Gγ subunits in heterotrimeric G-proteins). The proteins constitute a mixture of single-domain and multi-domain members. Although some of the structures in the PPC dataset are solved at a resolution poorer than 2.5 Å, the highest resolution of the Control-Rigid and Control-Monomer datasets, the magnitude of structural changes captured across the datasets can be compared since 50/76 complexes of the PPC dataset were solved with a resolution ≤2.5 Å. The conclusions of comparison of various parameters capturing structural change in the different datasets, discussed below, remained unaltered when using either the ‘50’ or ‘76’ set of complexes (data not shown). The conclusions described below are for the entire dataset of 76 complexes.
In order to distinguish structural variations caused due to protein binding from those occurring due to crystallographic artifacts, the upper bound values corresponding to the Control-Rigid dataset were used as reference for the three parameters (see Additional file 3: Figure S1).
It is observed that the main protein-protein complex dataset comprises of complexes with varying range of interface area and size of proteins (see Additional file 2: Table S2). Therefore, dependence of the parameters capturing structural change for interface area and size of the protein (represented as length of the protein) were analyzed. The analysis indicates that there is slight dependence of RMSD, PBc and PBSSc for interface area buried by the complex whereas the parameters have negligible dependence on the lengths of proteins (see Additional file 4: Figure S2).
Independently, we also captured structural change using all-atom RMSD which includes consideration of sidechain atoms. Although there is expected variation in all-atom RMSDs for interacting residues (see Additional file 5a: Figure S3), there is no profound variation when the Cα RMSDs are compared with corresponding all-atom RMSD values (see Additional file 5b: Figure S3). Therefore, the present analysis is confined to Cα RMSD based comparison in this study.
Pre-made interfaces predominantly bind to structurally-altered interfaces
Usually, interfaces with an average Cα RMSD of ≥1.5 Å showed substantial changes at interface, exemplifying the concept of induced-fit hypothesis  for formation of protein-protein complexes (see Additional file 6B: Figure S4). 35 interfaces with average Cα RMSD of ≥1.5 Å are found. Of these predominantly altered interfaces, 10 are partners of pre-made interfaces, 4 are partners of like-wise induced-fit interfaces and the rest have values in between (see Additional file 6C,D: Figure S4).
Comparison of the structural change in terms of Cα RMSD and normalized PB substitution score can help in distinguishing cases of rigid body movements from conformational changes. Induced-fit interface regions with 0% PB change at interface can be considered to have rigid-body movements (see Additional file 7A: Figure S5). However, since PBs are very sensitive to backbone torsion angle changes, two very similar PBs will also be considered as PB changes. Therefore, normalized PB substitution score is a more pertinent metric to grade the local conformational change (see Additional file 7B: Figure S5).
Non-interacting regions away from the interface undergo substantial structural changes on binding
Features of proteins with substantial structural change in non-interacting regions and no/moderate change at interface
Normalized PB substitution score
(NIR_PBc – IR)
(NIR_PBc – IR)
PIP3 kinase (O)
HISF protein (O)
Son of Sevenless (O)
HPr kinase C-ter domain (E)
DH/PH domain of TRIO (O)
TGFbeta receptor (O)
Vitamin D binding protein (O)
Sporulation response factor B (O)
Ran GTPase (O)
HEW lysozyme (A)
Transferrin receptor ectodomain (O)
Anthrax toxin receptor (O)
Xylanase inhibitor (E)
Complement C3 (O)
Proteins with substantial structural change in non-interacting regions and interfacial regions
Differences from global RMSD
Normalized PB substitution score
IR* – Global RMSD
NIR_PBc* – Global RMSD
Arf1 GTPase (O)
Ras GTPase (O)
CDK2 kinase (E)
FC fragment of human IgG 1 (A)
Ran GTPase (O)
Cystein protease (E)
Rab21 GTPase (O)
CDC42 GTPase (O)
Rac GTPase (O)
Rac GTPase (O)
Glycoprotein IB-alpha (E)
Conformational changes occurring away from the interface are potentially allosteric: Literature-based, structure-based and normal mode analysis
Features of non-interacting regions with substantial conformational change upon protein-protein interaction
PDB SITE record/Crystal atlas
Proteins with almost no structural change at the interface
PIP3 kinase (O)
Poor resolution and many missing residues
HISF protein (O)
Catalytic site residues not present in NIR_PBc/IR
Out of 4 active site residues, one in IR
NIR_PBc is near IR and are rigid in unbound form, one of the symmetry-related molecules is situated ~6Å towards this region.
Son of Sevenless (O)
Missing residues near this region
HPr kinase C-ter domain (E)
Resolution of PDBs is 2.8 Å
DH/PH domain of TRIO (O)
Temperature factors for region under consideration abnormally high!
Not so high
4/12 important residues are present in NIR_PBc
Region involved in dimerization in unbound form
Important residues present in IR
TGFbeta receptor (O)
3/5 of active site residues are present in NIR_PBc
Vitamin D binding protein (O)
NIR_PBc values are close to the ones observed are IR and rigid in unbound form, one of the 4 symmetry-related molecules comes ~7Å towards this region. Another NIR_PBc value is far away from IR and is slightly mobile - the same symmetry-related molecule comes <5Å close to this region!
Ran GTPase (O)
Nh (nearby IR), H
Sporulation response factor B (O)
Missing residues near this region
Ran GTPase (O)
One important residue present in IR
HEW lysozyme (A)
Catalytic site residues not present in NIR_PBc/IR
Transferrin receptor ectodoma in (O)
Important residues not present in IR/NIR_PBc
Anthrax toxin receptor (O)
Missing residues near this region
Xylanase inhibitor (E)
Catalytic site residues not present in NIR_PBc/IR
Complement C3 (O)
Important residues are not present in NIR_PBc/IR
Proteins with substantial structural change at interface
Arf1 GTPase (O)
Ran GTPase (O)
Important residues are present in NIR_PBC/IR
We observed that most of the proteins with changes are signalling proteins (17/22 in Case 1: Conformational changes in non-interacting surface regions of proteins with invariant interfaces & 8/12 in Case 2: Conformational changes in non-interacting surface regions of proteins with altered interfaces). Therefore, 25/40 ‘other’ complexes (predominantly signalling proteins) from the PPC dataset show significant structural changes at distal sites, indicating prevalence of this phenomenon in signalling proteins. In contrast, only 7/25 and 3/11 complexes of the enzyme-inhibitor and antigen-antibody classes, respectively, show such changes.
Detailed information about the residue positions and the nature of conformational change observed in the residues possibly forming the target site for all examples of Case 1 and Case 2 are listed in Additional file 10: Table S3 & Additional file 11: Table S4, respectively.
Results of normal mode analysis for proteins with ‘non-interacting regions with PB change’
Region of interest
Cases with substantial change in non-interacting residues whereas interacting region has very less change
NIR_PBc far away (mainly 86-91)
NIR_PBc far away (264-274)
NIR_PBc far away (120-123)
NIR_PBc near interface (191-194)
NIR_PBc far away (314-325)
NIR_PBc near interface (250-258)
NIR_PBc far away (156-158)
NIR_PBc near interface (89-93)
NIR_PBc with interface (100-105)
NIR_PBc near interface (15-22)
NIR_PBc near interface (152-157)
Cases with substantial change in non-interacting residues near interacting region as well as in interacting region
NIR_PBc far away (145-149)
NIR_PBc near interface (32-29)
TolB – Peptidoglycan associated lipoprotein (Pal) complex. The proteins TolB and Pal constitute the complex used by Escherichia coli and other ‘group A colicins’ to penetrate and kill cells . TolB protein comprises of two domains – a smaller N-terminal domain and a larger C-terminal β-propeller domain which interacts with Pal protein. Large rigid-body motions and conformational changes were seen far away from the interface between the unbound and bound forms of the TolB protein (Figure 6A, encircled region). The biological relevance of these changes is supported by experiments which prove that Pal binding results in conformational changes being transmitted to N terminal α/β domain of TolB . Further, a recent study shows that TolA binds to the N-terminal region exhibiting structural change in TolB , indicating that the structural changes occurring upon Pal binding serves as an allosteric signal. Additional support comes from GNM-based normal mode analysis of the unbound form. The region of our interest (mainly residues 86–91), was seen to have intrinsic tendency for reconfiguration in the second most significant mode pertaining to local motions (mode 3, Figure 6A, purple coloured region). However, we identify that this region is not evolutionary conserved and most of the sites occurred in the least conserved bin of residues in the protein.
Complement C3 and Epstein-Barr virus receptor C2 complex. Complement component C3d binds to antigenic molecules. This binding helps in further amplification of B cell responses as a consequence of the simultaneous binding of antigen-bound C3d with complement receptor type 2 and binding to B cell receptor via bound antigen . Complement C3’s interaction with C2 receptor causes conformational changes, identified using PBs, at residues 264–274. No literature information specific for this region was available. However, GNM analysis of the unbound form indicated that the region near to interface and the region in the opposite side have intrinsic motion (Figure 6B, purple coloured regions). The region of our interest occurs opposite to the interface and is indicated to be partly flexible, implying that this motion could be biologically relevant (Figure 6B, green coloured segment in purple coloured region). However, this region is moderately conserved.
Ran GTPase and Regulator of chromosome condensation (RCC1) complex. Ran GTPase is a key component of G-protein signaling. It serves as a molecular switch which cycles between GDP- and GTP- bound states. It requires regulators for enhancing its low intrinsic hydrolysis and nucleotide dissociation rates. Guanine nucleotide exchange factors form the latter group which bind to G-proteins and induce rapid dissociation of bound GTP and hence enable fast activation to GTP bound form. The structure under consideration is a complex of Ran GTPase with the guanine nucleotide exchange factor RCC1 .
Ran GTPase, which adopts the P-loop containing nucleoside triphosphate hydrolases fold, contains two regions of interest, one near the interface (residues 41–44) and the other far away from the interface (residues 151–155) (Figure 6C, encircled regions). The interface also undergoes substantial changes on binding. The region of interest near the interface is seen to be intrinsically mobile in the second most important mode. The region of interest far away from the interface is seen to be a part of a mobile region in the fourth important mode. This region is very near to GDP-binding site in the unbound form. It appears that binding of Ran GTPase and RCC1 causes structural changes at distant sites (GDP-binding site) to bring about exchange of nucleotides. This case provided a clear example of signal transduction within the molecule to bring about the desired biochemical effect. Such sites can also probably be targeted by human intervention to prevent disease manifestations, such as cancer in Ran signalling pathway . The sampling of homologous sequences was not diverse enough to obtain a reliable answer about its evolutionary conservation.
Structural changes away from the interface observed largely in proteins with structurally altered interfaces
Since structural changes away from the interface appear to be common on protein binding (n = 34/76), we studied the structures of the proteins in the complex to understand if there are any common characteristics of the protein exhibiting these changes (target protein) versus the binding partner (effector protein).
When we analysed the type of interface (pre-made, induced-fit, other) in the target protein, we observe that pre-made interfaces constitute only 7/34 complexes undergoing structural changes away from the interface. Induced-fit and moderately induced-fit (‘other’ interfaces) constitute the bulk, accounting for 15/34 and 13/34 complexes, respectively. Allostery appears to be the most plausible explanation to connect protein binding with structural changes away from the interface, as discussed in the previous section. Therefore, most of the proteins comprising of pre-made interfaces (26/33) do not appear to be the target proteins involved in allosteric communication. On the contrary, they seem to serve as the effector molecules for transmitting the allosteric signal to the partner protein. Comparison of the global RMSDs of case 1 complexes indicates that the target proteins (Cα RMSD = 1.52 ± 0.71 Å) are significantly more flexible than their effector proteins (Cα RMSD of 0.86 ± 0.60 Å with a p value of 0.0018 for a Wilcoxon paired test). However, the disparity in the values could be partly caused due to the differences in lengths (effector - 183 ± 91, target - 389 ± 247, Wilcoxon paired test p value - 0.0006) since RMSD is dependent on the number of residues. For proteins of case 2 complexes, there is no significant variation in terms of their lengths (Effector proteins – 236 ± 105, Target proteins – 245 ± 88, Wilcoxon paired test p value – 0.73). However, the global RMSDs of the target proteins (1.65 ± 0.64) are slightly yet significantly higher than those of the effector proteins (1.04 ± 0.67) (p value – 0.0342, Wilcoxon paired test). Therefore, it appears that binding of an effector protein causes changes at the interface of the target protein that are propagated towards the distant allosteric site, providing credence to the views implicating flexibility in allosteric modulation .
Availability of bound and unbound structures of proteins provides an opportunity to address various questions regarding structural alterations occurring due to protein-protein interactions. Our study underlines that macromolecular liganded forms of proteins undergo larger structural alterations in terms of change in local conformation (captured using PBs) as well as atomic positions (captured using RMSD) compared to unliganded proteins (Figure 1). These changes are much larger than those observed due to random fluctuations characteristic of intrinsic flexibility or experimental artifacts (see Additional file 3: Figure S1).
Non-obligatory complexes occupy a niche position as key regulators of cellular homeostasis. Their specific and timely association and dissociation are crucial for bringing about required biological function. Spatial and temporal regulation of the interacting proteins is one of the ways of avoiding unsuitable complexation . The other mechanism could be the use of different conformations of binding sites, which provide favourable or unfavourable binding-competence to the partner.
Transformation of binding site structure into the active form can serve as a switch to ensure correct binding at the appropriate time. Our analysis of structural alterations provides credence to this view. Surprisingly, pre-made interfaces, which are structurally invariant upon binding, shows distribution of %PB changes similar to that observed for induced-fit interfaces. This indicates that there is some extent of conformational change in all interfaces; only the nature and magnitude varies (Figure 2B). Additionally, interface of the partner of pre-made interface is usually observed to undergo significant structural changes (Figure 2A). In essence, there are no ‘completely pre-made’ interfaces in non-obligatory complexes. Crucially, significant structural changes are observed at backbone level in most of the interfaces used in this study. It is well known that side-chains undergo large structural changes upon protein-protein complexation . Considered together, these results support the view that structural conformations by themselves can serve as a good mechanism to implement the required tight regulation. Lower magnitude of structural changes is generally observed to optimize complex formation, whereas larger magnitude of structural changes is observed to remove steric clashes.
Flexible regions are dynamic and can take up several distinct conformations, which can have specific functional relevance . Several studies have revealed that flexibility is localized to certain regions of protein structure and such dynamic sites are usually involved in both small and large molecular interaction [90, 91] and enzymatic catalysis. In our study, the interface regions of TolB - Pal complex (Figure 6A) and Complement C3 - Epstein-Barr virus receptor C2 complex (Figure 6B) are shown by normal mode analysis to be intrinsically mobile.
Studies show that thermodynamic entropy redistribution is a common outcome of protein-protein interaction, irrespective of the net change in entropy after complexation . Loss of entropy at interacting sites is many times accompanied by gain of entropy in other regions of surface. ‘Entropy-entropy compensation’ may be due to significant intermolecular motion between the interacting molecules, which recovers about half of the entropy lost due to rotational and translational components . This compensatory mechanism has been postulated to be the mechanism responsible for high-specificity binding of multiple ligands at the same region of a protein.
The region may be functionally relevant, for e.g. a ligand/macromolecule binding site, whose conformation is regulated by an allosteric mechanism. Since binding sites are observed to be a combination of flexible and rigid sites , the signal based on protein-protein complexation may alter the stability and facilitate conformational change at the functionally relevant distant region. The complex of Ran GTPase with its cognate guanine nucleotide exchange factor probably utilizes this mechanism since complexation helps in altering the accessibility to the ligand on Rho protein (Figure 6C).
Crystallization is known to induce substantially altered conformations . In our study, we ensure that this bias is accounted for (Table 3) and that the conformational changes observed are not due to such effects.
Trivial factors, such as missing residues near the region of interest (or) the region being near termini, could contribute to such changes . Since we ruled out complexes exhibiting such changes (Table 3), the changes observed have other biological origin.
In-depth analysis of several complexes using rigorous coarse-grained NMA and literature survey indicates that a fair proportion of structural changes upon protein-protein complexation are allosteric (Figure 6, Tables 4). Such communication is largely enriched in signalling proteins, which seems plausible considering the complex regulation of signal transduction pathways achieved using the interplay of several modular elements . The lesser frequency of occurrence of such changes in enzyme-inhibitor and antibody-antigen complexes is expected. In the case of the former, their interaction is usually the result of an allosteric modulation and in the latter, a very high-affinity complex is formed, which needs to be cleared.
The classical view of allostery is as a mechanism of effector binding causing functionally relevant conformational changes at a distant site . The salient features of the models involves two key attributes: the presence of two conformational states of the protein, one stabilised in the unbound state and the other favoured upon binding of the allosteric effector, and induction of structural change at the target site leading to functional modulation. However, studies in the last two decades have thrown new light on this phenomenon. The observations of allosteric modulation in the absence of conformational change  and the introduction of allosteric perturbation in non-allosteric proteins  have raised the viewpoint that all dynamic proteins are possibly allosteric . These studies indicate that proteins in their unbound states exist in several conformational sub-states, characterized by different population densities . Allosteric perturbation results in change in the relative populations of these conformers . Such studies resulted in a paradigm shift in the understanding of allostery from a structure-centric to a thermodynamics-centric phenomenon . Although newer studies on allostery indicate that change in dynamics also enables allosteric communication in many cases [94–96, 98, 99], in this study we have confined ourselves to the study of allostery in the classical sense, as communicated by structural changes. Surprisingly, allosteric communication established only via structural changes appears to be established in almost half of the complexes upon protein binding. Consideration of dynamics along with structural changes would most probably lead to uncovering of many more protein-induced allosteric changes. Therefore, our study suggests that protein-protein binding in the case of signalling complexes, is often likely to result in downstream effects. The smaller of the two proteins in a complex, usually comprising of an unaltered interface upon protein-protein complexation, appears to be the effector molecule in most cases. The binding event generally causes changes at the interface and concomitant structural changes at the target site.
Signalling proteins are key drug targets and the usage of allosteric modulators as drugs is gaining acceptance [100, 101]. In such a scenario, the understanding that most protein-protein interactions in signalling proteins are allosteric provides impetus for the design of allosteric modulators as drugs. Allosteric regulators provide certain advantages over traditional drugs, which are usually competitive inhibitors. Binding of an allosteric drug at a distant site provides reduced side-effects, saturability, modulation in the presence of true agonist etc. [84, 100]. We hope that knowledge of possible allosterically modified sites identified in the signalling complexes studied in our analysis (see Additional file 10: Table S3 & Additional file 11: Table S4) serves as a starting point for combating disease manifestations.
Comparison of bound and unbound structures of protein-protein complexes enables us to address various questions regarding structural alterations occurring due to interaction. Non-obligatory complexes occupy a niche position as key regulators of cellular homeostasis with appropriate and timely association and dissociation which are crucial for eliciting the necessary biological function. Structural alterations in most of the interfaces of these non-obligatory complexes support the view that conformational features by themselves can serve as a good mechanism to implement the required tight regulation.
The interface is the most altered region in the entire protein structure upon protein-protein binding, as expected. The modifications are largely conformational in nature. In the rare case of one the partners remaining unaltered, the other partner is usually observed to undergo significant structural modification, thereby supporting the ‘induced fit hypothesis’  more than the ‘lock and key hypothesis’ .
The observation of a substantial proportion of instances with significant structural changes in non-interacting regions away from the interface implies that the binding is likely to result in downstream effects. In-depth analysis of several complexes using rigorous coarse-grained NMA and literature survey indicates that these changes have functional relevance, with most of them being allosteric. The observation of allostery-like structural changes in about half of the transient complexes suggests this phenomenon is much more prevalent in signalling complexes than appreciated before. It also appears that the reversible nature of protein-protein association and dissociation, characteristic of transient complexes, affords nature with an attractive means to bring about allostery which is generally a reversible process.
Rigid-proteins dataset (Control dataset 1): A dataset of 50 independently determined structures of two rigid proteins (see Additional file 1: Table S1), bovine ribonuclease (32 structures) and sperm whale myoglobin (18 structures), were taken from Rashin et.al . Values calculated from this dataset for different parameters are used as thresholds to account for positional coordinate uncertainty.
Monomeric-proteins dataset (Control dataset 2): To get a general idea about the flexibility in atomic positions for a random dataset, the PDB was mined for crystal structures of proteins with the following criteria: a single chain is present in the asymmetric unit and biological unit; crystallographic resolution of the structure should be 2.5 Å or better and the structure should not contain DNA, RNA, DNA-RNA hybrid, or other ligands bound to the protein. These molecules were clustered at a sequence identity of 95% and length coverage of 100% using BLASTCLUST (http://www.csc.fi/english/research/sciences/bioscience/programs/blast/blastclust). Finally, the clusters were refined to contain only one entry for each PubmedID per cluster, which ensures that mutants are not considered, to arrive at a dataset containing 95 clusters (see Additional file 2: Table S2) of 319 independently solved protein structures.
Protein-protein complex (PPC) dataset
The set of curated non-obligatory protein-protein interaction complexes solved in both unbound and bound form is taken from Benchmark 3.0 dataset . The set was further pruned using PISA  and PDB biological unit information to exclude cases containing different non-biological oligomeric forms of a protein in the unbound and bound forms (eg. X-X in unbound form and X-Y in bound form) and bound to other small ligands or peptides. All antibody-antigen complexes in the original dataset in which only the bound structure of the antibody was solved were discarded since the corresponding unbound form was not available. The final dataset consists of 76 non-obligatory complexes (see Additional file 2: Table S2). The resolution of these entries is 3.5 Å or better. Proteins in every interacting pair in the dataset is non-redundant at the level of SCOP family . Although a much larger dataset can be compiled if only one of the interacting proteins is available in unbound and bound form, such a dataset was not used since our objective is to compare the changes occurring in both the proteins upon complexation.
Although our dataset is intended to contain entries of identical proteins or protein domains available in both protein-bound and free forms, practically there could be some differences in the length and region of known 3-D structures in the bound and free forms. However the overwhelming majority of the same protein available in bound and free forms have >90% sequence identity (see Additional file 2: Table S2) indicating that the bound and unbound forms are almost the same. In all the cases with % sequence identity less than 90%, it is observed that the aligned region is identical or contains very few substitutions. Further, of the 3 cases showing large length variation between the bound and unbound forms (PDB codes: 1gcq, 1qa9, 1e6j) only 1e6j features in our analysis of cases showing structural changes away from the interface. So, it appears that the analysis is robust to length variations between bound and unbound forms of a protein.
As mentioned before the dataset used in the present analysis was derived from the robust list of protein-protein complexes proposed by Weng and coworkers  in their protein-protein docking benchmark version 3.0. In this dataset the authors have carefully avoided the complexes with significant extent of disordered regions. Indeed in the dataset used in the current analysis none of the complex structures used has any disordered residue at the protein-protein interfaces. This could be ensured on the basis of information on missing residues given in the PDB file, by checking the distance between Cα atoms of putative adjacent residues and by checking for the presence of all the expected atoms in a residue.
Identification of interfacial residues
If the distance between any two atoms of residues from the two proteins is less than sum of their van der Waals radii + 0.5 Å, the two residues are considered to be in the interface . The van der Waals radii were taken from the literature .
Classification of residues based on solvent accessibility
The residues in a structure are classified on the basis of their residue surface accessibility (RSA) which is calculated using NACCESS [104, 105]. This parameter provides a normalized measure of the accessible surface area of any residue in the protein, calculated with respect to the extended form of the residue, using the NACCESS program. The cut-offs employed are: ≤5% RSA (buried residues) and ≥10% RSA (surface residues). The 5% cut-off was adopted from , who optimized and used it to define residues buried in monomeric proteins. Buried, surface, and interface residues constitute ~25%, 75% and 10-20% of the residues in a protein, respectively.
Quantification of structural change
Structural change is estimated for a given residue in unbound and bound forms. A sequence alignment of the unbound and bound forms performed using CLUSTALW  provides the residue equivalences. Structural change is captured using two measures: RMSD and Protein Blocks. Structural change is classically captured by means of root mean square deviation (RMSD), where RMSD is calculated as follows: for i ranging from residue 1 to n of the dataset and d is the distance between N pairs of equivalent atoms. Two measures of RMSD have been employed: Cα RMSD and all-atom RMSD, based on deviation between the Cα positions of the same residue in unbound and bound forms in the former and between all-atoms of the same residue in unbound and bound forms for the latter. Deviation in side chain positions are generally expected  whereas large backbone changes are comparatively uncommon. Therefore, the deviation between the Cα positions of the same residue in unbound and bound form is used as an indicator of structural change mainly. The changes are captured at structural level and averaged out for the entire protein or a set of residues in a protein (for e.g. interface residues) and the averaged measures are used in the analysis. Small yet significant changes in local conformation of a protein can be captured using Protein Blocks . The three dimensional structural information in the bound and unbound forms is represented in a one-dimensional form using Protein Blocks (PBs). They consist of 16 structural prototypes, each of which approximates the backbone of a five-residue peptide. Given a 3D structure, each overlapping sequence of 5-residue fragments is associated with its closest PB. The sequence of PBs is annotated in the sequence alignment obtained using CLUSTALW. Two parameters are calculated using this measure. The first parameter indicates the presence of conformational change and is calculated as % changes in PBs between unbound and bound form (PBc). The second parameter indicates the magnitude of observed change and is calculated using PB substitution score (PBSSc) for the equivalent residues.
‘Pre-made’ versus ‘induced-fit’ interfaces
An interface with ≥0.5 Å Cα RMSD difference between the bound and unbound forms is classified as ‘pre-made’ interface whereas an interface with ≥1.5 Å Cα RMSD difference between the bound and unbound forms is classified as an ‘induced-fit’ interface. However, there are some interfaces with lower difference in terms of magnitude but with substantial difference at the interface in comparison to the rest of the surface residues (≥10% RSA). This cut-off was chosen since 90% of the interface residues have an RSA equal to or greater than this value in the unbound form. A normalization-based metric was used to identify induced-fit interfaces exhibiting smaller structural changes as where ‘CαRMSDI’ indicates the average Cα RMSD difference between bound and unbound form for interface, and ‘CαRMSDROS’ indicates the average Cα RMSD difference between bound and unbound form for the rest of the surface.
For example, a value of two indicates a doubled change in magnitude of the interface with respect to the rest of the surface. This value is used as a cut-off to identify substantial changes localized to the interface.
Identification of proteins with substantial structural change in non-interacting regions
To identify cases where the interface is largely invariant/moderately altered:
Criterion 1: Average Cα RMSD of ‘non-interacting residues with PB change’ should be ≥1 Å than the average Cα RMSD of the interface residues. For this comparison, the individual segments under consideration, were superimposed using SUPER (B.S.Neela, unpublished).
Criterion 2: Normalized PB substitution score of ‘non-interacting residues with PB change’ should be ≤ −2 than the normalized PB substitution score of interacting residues.
To identify cases where there are large changes at interface:
Criterion 1: Average Cα RMSD of both ‘non-interacting residues with PB change’ and interacting region should be ≥2 Å.
Criterion 2: Average Cα RMSD of both ‘non-interacting residues with PB change’ and interacting region should be > > global Cα RMSD.
Analysis of B-factors
The B-factor (temperature factor/atomic displacement factor) of an atom reflects the degree of isotropic smearing of electron density around its center . A low B-factor indicates small uncertainty in the position of an atom. A high B-factor can be caused by different factors: high thermal fluctuations, alternate conformation of an atom, and domain motion, to name a few.
To ascertain the flexibility/rigidity of a particular residue in a structure, its normalized backbone B-factor was considered . Normalization with respect to all the other residues provides an idea of increase/decrease in flexibility on a standard scale. Only surface residues (≥10% RSA) were considered for the normalization since all interacting residues and non-interacting surface residues form the crux of this study. The three most N-terminal and C-terminal surface residues were excluded since their B-factors are usually high and can affect the ‘mean’ of the values. B-factors of only backbone atoms were considered as we are studying backbone changes and also since side chain are generally more flexible than backbone atoms. The normalized B-factor per residue (B i,N ) was computed as where B i is the B-factor of residue i, <B i > is the mean B-factor of the protein surface residues and σ Bi is the s.d. for the same.
Residues with backbone B-factors ≥3, ≥2, and <1 standard deviations from the mean backbone B-factors for surface residues can be considered to have ‘very high’, ‘high’ and ‘low’ flexibility, respectively.
Identification of regions of protein structure with intrinsic collective motions
GNM-based NMA of the unbound form of a protein was undertaken to identify intrinsic collective motions of the molecule. In this model, the biomolecule is modelled as a harmonic oscillator with every residue represented as a single site, connected by springs to nearby residues . The oGNM web server  calculates low frequency normal modes for the unbound structure based on GNM. In GNM, the motions are isotropic by definition, thereby predicting only regions exhibiting intrinsic motion and magnitude of change. The directionality of motion cannot be predicted using GNM models. The server constructs the elastic network model of the structure by considering each of the Cα atoms as a node and identifying all interacting nodes using a distance cut-off of 10 Å. The six most low frequency modes were analyzed to check whether any of the ‘non-interacting regions with PB change’ far away from the interface show probable biologically relevant intrinsic motion.
Determining the extent of conservation of non-interacting residues with conformational changes
xBased on the assumption that evolutionary conservation of a site in a protein family is an indicator of its functional relevance and/or structural integrity, the degree of conservation of all sites in a protein family was calculated using the Jensen-Shannon divergence measure. This metric operates on the premise that most sites in a protein family are not under any evolutionary pressure and hence have a distribution similar to background amino acid distribution. Sites under evolutionary pressure, such as functional or stabilizing sites, show amino acid distribution significantly different from the background distribution.
Homologous sequences for every protein in our PPC dataset were identified by a search employing PSI-BLAST  against the UNIREF90  database at an e-value cutoff of 0.0001 for 3 iterations. Further, only sequences with ≥30% identity were considered. A multiple sequence alignment (MSA) of the query sequence with only the aligned regions of the homologous sequences was generated using CLUSTALW. The conservation scores for every site in the MSA was calculated using Jensen-Shannon divergence measure . The sites with top 30% conservation scores are considered to be well conserved  and sites with bottom 30% conservation scores are considered to be poorly conserved.
Generation of symmetry-related molecules using PyMOL
Symmetry-related molecules were generated for the bound and unbound molecule using PyMOL . The crystal packing after generation of symmetry-related molecules was checked to ascertain if any crystal packing could cause the observed structural changes in a complex.
This research is supported by Indo-French collaborative grant (CEFIPRA/IFCPAR number 3903-E) between AdB and NS. Part of this project is also supported for NS by Department of Biotechnology (DBT). Part of this was supported for AdB by grants from the Ministry of Research, University Paris Diderot – Paris 7, National Institute for Blood Transfusion (INTS) and the Institute for Health and Medical Research (INSERM). LSS and SM are supported by CEFIPRA/IFCPAR and DBT, respectively. The authors are grateful to Dr. Manoj Tyagi for providing the program to perform PB assignments and to Dr. Agnel P. Joseph for providing the PB substitution matrix. The authors express their immense gratitude to Professor Ruth Nussinov for key suggestions.
- Reichmann D, Rahat O, Cohen M, Neuvirth H, Schreiber G: The molecular architecture of protein-protein binding sites. Curr Opin Struct Biol 2007, 17(1):67–76.PubMedView ArticleGoogle Scholar
- Levy ED, Pereira-Leal JB: Evolution and dynamics of protein interactions and networks. Curr Opin Struct Biol 2008, 18(3):349–357.PubMedView ArticleGoogle Scholar
- Vidal M, Cusick ME, Barabasi AL: Interactome networks and human disease. Cell 2011, 144(6):986–998.PubMed CentralPubMedView ArticleGoogle Scholar
- Schreiber G, Keating AE: Protein binding specificity versus promiscuity. Curr Opin Struct Biol 2011, 21(1):50–61.PubMed CentralPubMedView ArticleGoogle Scholar
- Janin J, Wodak SJ: Protein modules and protein-protein interaction. Introduction. Adv Protein Chem 2002, 61: 1–8.PubMedView ArticleGoogle Scholar
- Ryan DP, Matthews JM: Protein-protein interactions in human disease. Curr Opin Struct Biol 2005, 15(4):441–446.PubMedView ArticleGoogle Scholar
- Swain JF, Gierasch LM: The changing landscape of protein allostery. Curr Opin Struct Biol 2006, 16(1):102–108.PubMedView ArticleGoogle Scholar
- Cui Q, Karplus M: Allostery and cooperativity revisited. Protein Sci 2008, 17(8):1295–1307.PubMed CentralPubMedView ArticleGoogle Scholar
- Laskowski RA, Gerick F, Thornton JM: The structural basis of allosteric regulation in proteins. FEBS Lett 2009, 583(11):1692–1698.PubMedView ArticleGoogle Scholar
- Tsai CJ, Del Sol A, Nussinov R: Protein allostery, signal transmission and dynamics: a classification scheme of allosteric mechanisms. Mol Biosyst 2009, 5(3):207–216.PubMed CentralPubMedView ArticleGoogle Scholar
- Tzeng SR, Kalodimos CG: Protein dynamics and allostery: an NMR view. Curr Opin Struct Biol 2011, 21(1):62–67.PubMedView ArticleGoogle Scholar
- Changeux JP, Edelstein SJ: Allosteric mechanisms of signal transduction. Science 2005, 308(5727):1424–1428.PubMedView ArticleGoogle Scholar
- Beckett D: Regulating transcription regulators via allostery and flexibility. Proc Natl Acad Sci U S A 2009, 106(52):22035–22036.PubMed CentralPubMedView ArticleGoogle Scholar
- Drewes G, Bouwmeester T: Global approaches to protein-protein interactions. Curr Opin Cell Biol 2003, 15(2):199–205.PubMedView ArticleGoogle Scholar
- Piehler J: New methodologies for measuring protein interactions in vivo and in vitro. Curr Opin Struct Biol 2005, 15(1):4–14.PubMedView ArticleGoogle Scholar
- Phizicky EM, Fields S: Protein-protein interactions: methods for detection and analysis. Microbiol Rev 1995, 59(1):94–123.PubMed CentralPubMedGoogle Scholar
- Lakey JH, Raggett EM: Measuring protein-protein interactions. Curr Opin Struct Biol 1998, 8(1):119–123.PubMedView ArticleGoogle Scholar
- von Mering C, Krause R, Snel B, Cornell M, Oliver SG, Fields S, Bork P: Comparative assessment of large-scale data sets of protein-protein interactions. Nature 2002, 417(6887):399–403.PubMedView ArticleGoogle Scholar
- Henrick K, Feng Z, Bluhm WF, Dimitropoulos D, Doreleijers JF, Dutta S, Flippen-Anderson JL, Ionides J, Kamada C, Krissinel E, et al.: Remediation of the protein data bank archive. Nucleic Acids Res 2008, 36(Database issue):D426-D433.PubMed CentralPubMedGoogle Scholar
- Janin J, Rodier F, Chakrabarti P, Bahadur RP: Macromolecular recognition in the Protein Data Bank. Acta Crystallogr D: Biol Crystallogr 2007, 63(Pt 1):1–8.View ArticleGoogle Scholar
- Bordner AJ, Gorin AA: Comprehensive inventory of protein complexes in the Protein Data Bank from consistent classification of interfaces. BMC Bioinf 2008, 9: 234.View ArticleGoogle Scholar
- Kim WK, Henschel A, Winter C, Schroeder M: The many faces of protein-protein interactions: A compendium of interface geometry. PLoS Comput Biol 2006, 2(9):e124.PubMed CentralPubMedView ArticleGoogle Scholar
- Finn RD, Marshall M, Bateman A: iPfam: visualization of protein-protein interactions in PDB at domain and amino acid resolutions. Bioinformatics 2005, 21(3):410–412.PubMedView ArticleGoogle Scholar
- Jefferson ER, Walsh TP, Roberts TJ, Barton GJ: SNAPPI-DB: a database and API of Structures, iNterfaces and Alignments for Protein-Protein Interactions. Nucleic Acids Res 2007, 35(Database issue):D580-D589.PubMed CentralPubMedView ArticleGoogle Scholar
- Levy ED, Pereira-Leal JB, Chothia C, Teichmann SA: 3D complex: a structural classification of protein complexes. PLoS Comput Biol 2006, 2(11):e155.PubMed CentralPubMedView ArticleGoogle Scholar
- Faure G, Andreani J, Guerois R: InterEvol database: exploring the structure and evolution of protein complex interfaces. Nucleic Acids Res 2011, 40(Database issue):D847-D856.PubMed CentralPubMedGoogle Scholar
- Xu Q, Dunbrack RL: The protein common interface database (ProtCID)--a comprehensive database of interactions of homologous proteins in multiple crystal forms. Nucleic Acids Res 2011, 39(Database issue):D761-D770.PubMed CentralPubMedView ArticleGoogle Scholar
- Nooren IM, Thornton JM: Diversity of protein-protein interactions. EMBO J 2003, 22(14):3486–3492.PubMed CentralPubMedView ArticleGoogle Scholar
- Ansari S, Helms V: Statistical analysis of predominantly transient protein-protein interfaces. Proteins 2005, 61(2):344–355.PubMedView ArticleGoogle Scholar
- Mintseris J, Weng Z: Structure, function, and evolution of transient and obligate protein-protein interactions. Proc Natl Acad Sci U S A 2005, 102(31):10930–10935.PubMed CentralPubMedView ArticleGoogle Scholar
- Perkins JR, Diboun I, Dessailly BH, Lees JG, Orengo C: Transient protein-protein interactions: structural, functional, and network properties. Structure 2010, 18(10):1233–1243.PubMedView ArticleGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res 2000, 28(1):235–242.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen R, Mintseris J, Janin J, Weng Z: A protein-protein docking benchmark. Proteins 2003, 52(1):88–91.PubMedView ArticleGoogle Scholar
- Hwang H, Pierce B, Mintseris J, Janin J, Weng Z: Protein-protein docking benchmark version 3.0. Proteins 2008, 73(3):705–709.PubMed CentralPubMedView ArticleGoogle Scholar
- Lobanov MY, Shoemaker BA, Garbuzynskiy SO, Fong JH, Panchenko AR, Galzitskaya OV: ComSin: database of protein structures in bound (complex) and unbound (single) states in relation to their intrinsic disorder. Nucleic Acids Res 2010, 38(Database issue):D283-D287.PubMed CentralPubMedView ArticleGoogle Scholar
- Betts MJ, Sternberg MJ: An analysis of conformational changes on protein-protein association: implications for predictive docking. Protein Eng 1999, 12(4):271–283.PubMedView ArticleGoogle Scholar
- Martin J, Regad L, Lecornet H, Camproux AC: Structural deformation upon protein-protein interaction: a structural alphabet approach. BMC Struct Biol 2008, 8: 12.PubMed CentralPubMedView ArticleGoogle Scholar
- Martin J, Regad L, Etchebest C, Camproux AC: Taking advantage of local structure descriptors to analyze interresidue contacts in protein structures and protein complexes. Proteins 2008, 73(3):672–689.PubMedView ArticleGoogle Scholar
- Agarwal A, Dinesh DC, Srinivasan N, de Brevern AG: Characterization of conformational patterns in active and inactive forms of kinases using protein blocks approach. In Computational Intelligence and Pattern Analysis in Biological Informatics. Edited by: Maulik U, Bandyopadhyay S, Wang JTL. Wiley, New Jersey; 2010:169–187.View ArticleGoogle Scholar
- Smith GR, Sternberg MJ, Bates PA: The relationship between the flexibility of proteins and their conformational states on forming protein-protein complexes with an application to protein-protein docking. J Mol Biol 2005, 347(5):1077–1101.PubMedView ArticleGoogle Scholar
- Grunberg R, Leckner J, Nilges M: Complementarity of structure ensembles in protein-protein binding. Structure 2004, 12(12):2125–2136.PubMedView ArticleGoogle Scholar
- Grant BJ, Gorfe AA, McCammon JA: Large conformational changes in proteins: signaling and other functions. Curr Opin Struct Biol 2010, 20(2):142–147.PubMed CentralPubMedView ArticleGoogle Scholar
- Zacharias M: Accounting for conformational changes during protein-protein docking. Curr Opin Struct Biol 2010, 20(2):180–186.PubMedView ArticleGoogle Scholar
- Rashin AA, Rashin AH, Jernigan RL: Protein flexibility: coordinate uncertainties and interpretation of structural differences. Acta Crystallogr D: Biol Crystallogr 2009, 65(Pt 11):1140–1161.View ArticleGoogle Scholar
- Teilum K, Olsen JG, Kragelund BB: Functional aspects of protein flexibility. Cell Mol Life Sci 2009, 66(14):2231–2247.PubMedView ArticleGoogle Scholar
- Juritz EI, Alberti SF, Parisi GD: PCDB: a database of protein conformational diversity. Nucleic Acids Res 2011, 39(Database issue):D475-D479.PubMed CentralPubMedView ArticleGoogle Scholar
- Murzin AG, Brenner SE, Hubbard T, Chothia C: SCOP: a structural classification of proteins database for the investigation of sequences and structures. J Mol Biol 1995, 247(4):536–540.PubMedGoogle Scholar
- de Brevern AG, Etchebest C, Hazout S: Bayesian probabilistic approach for predicting backbone structures in terms of protein blocks. Proteins 2000, 41(3):271–287.PubMedView ArticleGoogle Scholar
- Joseph AP, Agarwal G, Mahajan S, Gelly J-C, Swapna LS, Offmann B, Cadet F, Bornot A, Tyagi M, Valadié H, et al.: A short survey on Protein Blocks. Biophys Rev 2010, 2(3):137–145.PubMed CentralPubMedView ArticleGoogle Scholar
- Joseph AP, Srinivasan N, de Brevern AG: Improvement of protein structure comparison using a structural alphabet. Biochimie 2011, 93(9):1434–1445.PubMedView ArticleGoogle Scholar
- Fischer E: Einfluss der configuration auf die wirkung der enzyme. Ber Dtsch Chem Ges 1894, 27: 2985–2993.View ArticleGoogle Scholar
- Koshland DE: Application of a Theory of Enzyme Specificity to Protein Synthesis. Proc Natl Acad Sci U S A 1958, 44(2):98–104.PubMed CentralPubMedView ArticleGoogle Scholar
- Keskin O, Tsai CJ, Wolfson H, Nussinov R: A new, structurally nonredundant, diverse data set of protein-protein interfaces and its implications. Protein Sci 2004, 13(4):1043–1055.PubMed CentralPubMedView ArticleGoogle Scholar
- Jones S, Thornton JM: Protein-protein interactions: a review of protein dimer structures. Prog Biophys Mol Biol 1995, 63(1):31–65.PubMedView ArticleGoogle Scholar
- Porter CT, Bartlett GJ, Thornton JM: The Catalytic Site Atlas: a resource of catalytic sites and residues identified in enzymes using structural data. Nucleic Acids Res 2004, 32(Database issue):D129-D133.PubMed CentralPubMedView ArticleGoogle Scholar
- Pacold ME, Suire S, Perisic O, Lara-Gonzalez S, Davis CT, Walker EH, Hawkins PT, Stephens L, Eccleston JF, Williams RL: Crystal structure and functional analysis of Ras binding to its effector phosphoinositide 3-kinase gamma. Cell 2000, 103(6):931–943.PubMedView ArticleGoogle Scholar
- Douangamath A, Walker M, Beismann-Driemeyer S, Vega-Fernandez MC, Sterner R, Wilmanns M: Structural evidence for ammonia tunneling across the (beta alpha)(8) barrel of the imidazole glycerol phosphate synthase bienzyme complex. Structure 2002, 10(2):185–193.PubMedView ArticleGoogle Scholar
- Cheetham JC, Smith DM, Aoki KH, Stevenson JL, Hoeffel TJ, Syed RS, Egrie J, Harvey TS: NMR structure of human erythropoietin and a comparison with its receptor bound conformation. Nat Struct Biol 1998, 5(10):861–866.PubMedView ArticleGoogle Scholar
- Johnston SC, Larsen CN, Cook WJ, Wilkinson KD, Hill CP: Crystal structure of a deubiquitinating enzyme (human UCH-L3) at 1.8 A resolution. EMBO J 1997, 16(13):3787–3796.PubMed CentralPubMedView ArticleGoogle Scholar
- Freedman TS, Sondermann H, Friedland GD, Kortemme T, Bar-Sagi D, Marqusee S, Kuriyan J: A Ras-induced conformational switch in the Ras activator Son of sevenless. Proc Natl Acad Sci U S A 2006, 103(45):16692–16697.PubMed CentralPubMedView ArticleGoogle Scholar
- Hart PJ, Deep S, Taylor AB, Shu Z, Hinck CS, Hinck AP: Crystal structure of the human TbetaR2 ectodomain–TGF-beta3 complex. Nat Struct Biol 2002, 9(3):203–208.PubMedGoogle Scholar
- Fieulaine S, Morera S, Poncet S, Mijakovic I, Galinier A, Janin J, Deutscher J, Nessler S: X-ray structure of a bifunctional protein kinase in complex with its protein substrate HPr. Proc Natl Acad Sci U S A 2002, 99(21):13437–13441.PubMed CentralPubMedView ArticleGoogle Scholar
- Chhatriwala MK, Betts L, Worthylake DK, Sondek J: The DH and PH domains of Trio coordinately engage Rho GTPases for their efficient activation. J Mol Biol 2007, 368(5):1307–1320.PubMed CentralPubMedView ArticleGoogle Scholar
- Dementiev A, Simonovic M, Volz K, Gettins PG: Canonical inhibitor-like interactions explain reactivity of alpha1-proteinase inhibitor Pittsburgh and antithrombin with proteinases. J Biol Chem 2003, 278(39):37881–37887.PubMedView ArticleGoogle Scholar
- Huse M, Muir TW, Xu L, Chen YG, Kuriyan J, Massague J: The TGF beta receptor activation process: an inhibitor- to substrate-binding switch. Mol Cell 2001, 8(3):671–682.PubMedView ArticleGoogle Scholar
- Zapf J, Sen U: Madhusudan, Hoch JA, Varughese KI: A transient interaction between two phosphorelay proteins trapped in a crystal lattice reveals the mechanism of molecular recognition and phosphotransfer in signal transduction. Structure 2000, 8(8):851–862.PubMedView ArticleGoogle Scholar
- Stewart M, Kent HM, McCoy AJ: Structural basis for molecular recognition between nuclear transport factor 2 (NTF2) and the GDP-bound form of the Ras-family GTPase Ran. J Mol Biol 1998, 277(3):635–646.PubMedView ArticleGoogle Scholar
- Bennett MJ, Lebron JA, Bjorkman PJ: Crystal structure of the hereditary haemochromatosis protein HFE complexed with transferrin receptor. Nature 2000, 403(6765):46–53.PubMedView ArticleGoogle Scholar
- Sondermann P, Huber R, Oosthuizen V, Jacob U: The 3.2-A crystal structure of the human IgG1 Fc fragment-Fc gammaRIII complex. Nature 2000, 406(6793):267–273.PubMedView ArticleGoogle Scholar
- Delprato A, Lambright DG: Structural basis for Rab GTPase activation by VPS9 domain exchange factors. Nat Struct Mol Biol 2007, 14(5):406–412.PubMed CentralPubMedView ArticleGoogle Scholar
- Bahar I, Lezon TR, Yang LW, Eyal E: Global dynamics of proteins: bridging between structure and function. Annu Rev Biophys 2010, 39: 23–42.PubMed CentralPubMedView ArticleGoogle Scholar
- Ma J: Usefulness and limitations of normal mode analysis in modeling dynamics of biomolecular complexes. Structure 2005, 13(3):373–380.PubMedView ArticleGoogle Scholar
- Bahar I, Rader AJ: Coarse-grained normal mode analysis in structural biology. Curr Opin Struct Biol 2005, 15(5):586–592.PubMed CentralPubMedView ArticleGoogle Scholar
- Tobi D, Bahar I: Structural changes involved in protein binding correlate with intrinsic motions of proteins in the unbound state. Proc Natl Acad Sci U S A 2005, 102(52):18908–18913.PubMed CentralPubMedView ArticleGoogle Scholar
- Yang LW, Rader AJ, Liu X, Jursa CJ, Chen SC, Karimi HA, Bahar I: oGNM: online computation of structural dynamics using the Gaussian Network Model. Nucleic Acids Res 2006, 34(Web Server issue):W24-W31.PubMed CentralPubMedView ArticleGoogle Scholar
- Chennubhotla C, Bahar I: Signal propagation in proteins and relation to equilibrium fluctuations. PLoS Comput Biol 2007, 3(9):1716–1726.PubMedView ArticleGoogle Scholar
- Delano WL: The PyMOL Molecular Graphics System. DeLano Scientific, San Carlos, CA, USA; 2002.Google Scholar
- Abergel C, Bouveret E, Claverie JM, Brown K, Rigal A, Lazdunski C, Benedetti H: Structure of the Escherichia coli TolB protein determined by MAD methods at 1.95 A resolution. Structure 1999, 7(10):1291–1300.PubMedView ArticleGoogle Scholar
- Bonsor DA, Grishkovskaya I, Dodson EJ, Kleanthous C: Molecular mimicry enables competitive recruitment by a natively disordered protein. J Am Chem Soc 2007, 129(15):4800–4807.PubMedView ArticleGoogle Scholar
- Bonsor DA, Hecht O, Vankemmelbeke M, Sharma A, Krachler AM, Housden NG, Lilly KJ, James R, Moore GR, Kleanthous C: Allosteric beta-propeller signalling in TolB and its manipulation by translocating colicins. EMBO J 2009, 28(18):2846–2857.PubMed CentralPubMedView ArticleGoogle Scholar
- Szakonyi G, Guthridge JM, Li D, Young K, Holers VM, Chen XS: Structure of complement receptor 2 in complex with its C3d ligand. Science 2001, 292(5522):1725–1728.PubMedView ArticleGoogle Scholar
- Renault L, Kuhlmann J, Henkel A, Wittinghofer A: Structural basis for guanine nucleotide exchange on Ran by the regulator of chromosome condensation (RCC1). Cell 2001, 105(2):245–255.PubMedView ArticleGoogle Scholar
- Sanderson HS, Clarke PR: Cell biology: Ran, mitosis and the cancer connection. Curr Biol 2006, 16(12):R466-R468.PubMedView ArticleGoogle Scholar
- del Sol A, Tsai CJ, Ma B, Nussinov R: The origin of allosteric functional modulation: multiple pre-existing pathways. Structure 2009, 17(8):1042–1050.PubMed CentralPubMedView ArticleGoogle Scholar
- Stein A, Pache RA, Bernado P, Pons M, Aloy P: Dynamic interactions of proteins in complex networks: a more structured view. FEBS J 2009, 276(19):5390–5405.PubMedView ArticleGoogle Scholar
- Ruvinsky AM, Kirys T, Tuzikov AV, Vakser IA: Side-chain conformational changes upon Protein-Protein Association. J Mol Biol 2011, 408(2):356–365.PubMed CentralPubMedView ArticleGoogle Scholar
- Fayos R, Melacini G, Newlon MG, Burns L, Scott JD, Jennings PA: Induction of flexibility through protein-protein interactions. J Biol Chem 2003, 278(20):18581–18587.PubMedView ArticleGoogle Scholar
- Yu L, Zhu CX, Tse-Dinh YC, Fesik SW: Backbone dynamics of the C-terminal domain of Escherichia coli topoisomerase I in the absence and presence of single-stranded DNA. Biochemistry 1996, 35(30):9661–9666.PubMedView ArticleGoogle Scholar
- Arumugam S, Gao G, Patton BL, Semenchenko V, Brew K, Van Doren SR: Increased backbone mobility in beta-barrel enhances entropy gain driving binding of N-TIMP-1 to MMP-3. J Mol Biol 2003, 327(3):719–734.PubMedView ArticleGoogle Scholar
- Luque I, Freire E: Structural stability of binding sites: consequences for binding affinity and allosteric effects. Proteins 2000, 41(S4):63–71.View ArticleGoogle Scholar
- Teague SJ: Implications of protein flexibility for drug discovery. Nat Rev Drug Discov 2003, 2(7):527–541.PubMedView ArticleGoogle Scholar
- Grunberg R, Nilges M, Leckner J: Flexibility and conformational entropy in protein-protein binding. Structure 2006, 14(4):683–693.PubMedView ArticleGoogle Scholar
- Faure G, Bornot A, de Brevern AG: Analysis of protein contacts into Protein Units. Biochimie 2009, 91(7):876–887.PubMedView ArticleGoogle Scholar
- Popovych N, Sun S, Ebright RH, Kalodimos CG: Dynamically driven protein allostery. Nat Struct Mol Biol 2006, 13(9):831–838.PubMed CentralPubMedView ArticleGoogle Scholar
- Ikeda Y, Taniguchi N, Noguchi T: Dominant negative role of the glutamic acid residue conserved in the pyruvate kinase M(1) isozyme in the heterotropic allosteric effect involving fructose-1,6-bisphosphate. J Biol Chem 2000, 275(13):9150–9156.PubMedView ArticleGoogle Scholar
- Gunasekaran K, Ma B, Nussinov R: Is allostery an intrinsic property of all dynamic proteins? Proteins 2004, 57(3):433–443.PubMedView ArticleGoogle Scholar
- Lange OF, Lakomek NA, Fares C, Schroder GF, Walter KF, Becker S, Meiler J, Grubmuller H, Griesinger C, de Groot BL: Recognition dynamics up to microseconds revealed from an RDC-derived ubiquitin ensemble in solution. Science 2008, 320(5882):1471–1475.PubMedView ArticleGoogle Scholar
- Tsai CJ, del Sol A, Nussinov R: Allostery: absence of a change in shape does not imply that allostery is not at play. J Mol Biol 2008, 378(1):1–11.PubMed CentralPubMedView ArticleGoogle Scholar
- Mitternacht S, Berezovsky IN: Binding leverage as a molecular basis for allosteric regulation. PLoS Comput Biol 2011, 7(9):e1002148.PubMed CentralPubMedView ArticleGoogle Scholar
- Kar G, Keskin O, Gursoy A, Nussinov R: Allostery and population shift in drug discovery. Curr Opin Pharmacol 2010, 10(6):715–722.PubMedView ArticleGoogle Scholar
- Acuner Ozbabacan SE, Gursoy A, Keskin O, Nussinov R: Conformational ensembles, signal transduction and residue hot spots: application to drug discovery. Curr Opin Drug Discov Devel 2010, 13(5):527–537.PubMedGoogle Scholar
- Krissinel E, Henrick K: Inference of macromolecular assemblies from crystalline state. J Mol Biol 2007, 372(3):774–797.PubMedView ArticleGoogle Scholar
- Chothia C: Structural invariants in protein folding. Nature 1975, 254(5498):304–308.PubMedView ArticleGoogle Scholar
- Lee B, Richards FM: The interpretation of protein structures: estimation of static accessibility. J Mol Biol 1971, 55(3):379–400.PubMedView ArticleGoogle Scholar
- Hubbard SJ, Campbell SF, Thornton JM: Molecular recognition. Conformational analysis of limited proteolytic sites and serine proteinase protein inhibitors. J Mol Biol 1991, 220(2):507–530.PubMedView ArticleGoogle Scholar
- Miller S, Lesk AM, Janin J, Chothia C: The accessible surface area and stability of oligomeric proteins. Nature 1987, 328(6133):834–836.PubMedView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994, 22(22):4673–4680.PubMed CentralPubMedView ArticleGoogle Scholar
- Parthasarathy S, Murthy MR: Analysis of temperature factor distribution in high-resolution protein structures. Protein Sci 1997, 6(12):2561–2567.PubMed CentralPubMedView ArticleGoogle Scholar
- Yuan Z, Zhao J, Wang ZX: Flexibility analysis of enzyme active sites by crystallographic temperature factors. Protein Eng 2003, 16(2):109–114.PubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25(17):3389–3402.PubMed CentralPubMedView ArticleGoogle Scholar
- Suzek BE, Huang H, McGarvey P, Mazumder R, Wu CH: UniRef: comprehensive and non-redundant UniProt reference clusters. Bioinformatics 2007, 23(10):1282–1288.PubMedView ArticleGoogle Scholar
- Capra JA, Singh M: Predicting functionally important residues from sequence conservation. Bioinformatics 2007, 23(15):1875–1882.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.