Dissecting subunit interfaces in homodimeric proteins

Crystal structures of barley thioredoxin h isoforms HvTrxh1 and HvTrxh2 reveal features involved in protein recognition and possibly in discriminating the isoform specificity

H-type thioredoxins (Trxs) constitute a particularly large Trx sub-group in higher plants. Here, the crystal structures are determined for the two barley Trx h isoforms, HvTrxh1 and HvTrxh2, in the partially radiation-reduced state to resolutions of 1.7 A, and for HvTrxh2 in the oxidized state to 2.0 A. The two Trxs have a sequence identity of 51% and highly similar fold and active-site architecture. Interestingly, the four independent molecules in the crystals of HvTrxh1 form two relatively large and essentially identical protein-protein interfaces. In each interface, a loop segment of one HvTrxh1 molecule is positioned along a shallow hydrophobic groove at the primary nucleophile Cys40 of another HvTrxh1 molecule. The association mode can serve as a model for the target protein recognition by Trx, as it brings the Met82 Cgamma atom (gamma position as a disulfide sulfur) of the bound loop segment in the proximity of the Cys40 thiol. The interaction involves three characteristic backbone-backbone hydrogen bonds in an antiparallel beta-sheet-like arrangement, similar to the arrangement observed in the structure of an engineered, covalently bound complex between Trx and a substrate protein, as reported by Maeda et al. in an earlier paper. The occurrence of an intermolecular salt bridge between Glu80 of the bound loop segment and Arg101 near the hydrophobic groove suggests that charge complementarity plays a role in the specificity of Trx. In HvTrxh2, isoleucine corresponds to this arginine, which emphasizes the potential for specificity differences between the coexisting barley Trx isoforms.

Automated technologies and novel techniques to accelerate protein crystallography for structural genomics

The sequence infrastructure that has arisen through large-scale genomic projects dedicated to protein analysis, has provided a wealth of information and brought together scientists and institutions from all over the world. As a consequence, the development of novel technologies and methodologies in proteomics research is helping to unravel the biochemical and physiological mechanisms of complex multivariate diseases at both a functional and molecular level. In the late sixties, when X-ray crystallography had just been established, the idea of determining protein structure on an almost universal basis was akin to an impossible dream or a miracle. Yet only forty years after, automated protein structure determination platforms have been established. The widespread use of robotics in protein crystallography has had a huge impact at every stage of the pipeline from protein cloning, over-expression, purification, crystallization, data collection, structure solution, refinement, validation and data management- all of which have become more or less automated with minimal human intervention necessary. Here, recent advances in protein crystal structure analysis in the context of structural genomics will be discussed. In addition, this review aims to give an overview of recent developments in high throughput instrumentation, and technologies and strategies to accelerate protein structure/function analysis.

Dissecting protein-RNA recognition sites

We analyze the protein–RNA interfaces in 81 transient binary complexes taken from the Protein Data Bank. Those with tRNA or duplex RNA are larger than with single-stranded RNA, and comparable in size to protein–DNA interfaces. The protein side bears a strong positive electrostatic potential and resembles protein–DNA interfaces in its amino acid composition. On the RNA side, the phosphate contributes less, and the sugar much more, to the interaction than in protein–DNA complexes. On average, protein–RNA interfaces contain 20 hydrogen bonds, 7 that involve the phosphates, 5 the sugar 2′OH, and 6 the bases, and 32 water molecules. The average H-bond density per unit buried surface area is less with tRNA or single-stranded RNA than with duplex RNA. The atomic packing is also less compact in interfaces with tRNA. On the protein side, the main chain NH and Arg/Lys side chains account for nearly half of all H-bonds to RNA; the main chain CO and side chain acceptor groups, for a quarter. The 2′OH is a major player in protein–RNA recognition, and shape complementarity an important determinant, whereas electrostatics and direct base–protein interactions play a lesser part than in protein–DNA recognition.

PiQSi: protein quaternary structure investigation

PiQSi facilitates the manual investigation of the quaternary structure of protein complexes in the Protein Data Bank (PDB). Users can browse and obtain an overview of the quaternary structure information of a given protein together with its evolutionary relatives, which helps in the determination of the biological quaternary state. I have used this framework to annotate over 10,000 structures from the PDB Biological Unit and corrected the quaternary state of approximately 15% of them. A benchmark shows that the annotations are of high quality and stresses the need for manual curation, in particular for ambiguous cases such as proteins in equilibrium between two quaternary states. The approximately 10,000 annotations already in the database can be used to improve the accuracy of analyses on protein structure or to benchmark methods that predict protein quaternary structure. In addition, PiQSi incorporates a community-based curation system, which I hope will allow us to reach an accurate and complete description of the biological quaternary state of proteins in PDB. PiQSi is accessible at http://www.PiQSi.org/.

The origin of protein interactions and allostery in colocalization

Two fundamental principles can account for how regulated networks of interacting proteins originated in cells. These are the law of mass action, which holds that the binding of one molecule to another increases with concentration, and the fact that the colocalization of molecules vastly increases their local concentrations. It follows that colocalization can amplify the effect on one protein of random mutations in another protein and can therefore, through natural selection, lead to interactions between proteins and to a startling variety of complex allosteric controls. It also follows that allostery is common and that homologous proteins can have different allosteric mechanisms. Thus, the regulated protein networks of organisms seem to be the inevitable consequence of natural selection operating under physical laws.

DiMoVo: a Voronoi tessellation-based method for discriminating crystallographic and biological protein-protein interactions

MOTIVATION

Knowledge of the oligomeric state of a protein is often essential for understanding its function and mechanism. Within a protein crystal, each protein monomer is in contact with many others, forming many small interfaces and a few larger ones that are biologically significant if the protein is a homodimer in solution, but not if the protein is monomeric. Telling such 'crystal dimers' from real ones remains a difficult task.

RESULTS

It has already been demonstrated that the interfaces of native and non-native protein-protein complexes can be distinguished using a combination of parameters computed with a method on the Voronoi tessellation. We show in this article that the same parameters highlight significant differences between the interfaces of biological and crystal dimers. Using these parameters as descriptors in machine learning methods leads to accurate classification of specific and non-specific protein-protein interfaces.

AVAILABILITY

Software is available at http://fifi.ibbmc.u-psud.fr/DiMoVo.

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.

PPIDD: an extraction and visualisation method of biological protein-protein interfaces

Today, the information for generating reliable protein-protein complex datasets is not directly accessible from PDB structures. Moreover, in X-ray protein structures, different types of contacts can be observed between proteins: contacts in homodimers or inside heterocomplexes considered to be specific, and contacts induced by crystallogenesis processes, considered to be non-specific. However, none of the databases giving access to protein-protein complexes allows the crystallographic interfaces to be distinguished from the biological interfaces. For this reason we developed PPIDD (Protein-Protein Interface Description Database), an innovative tool, which allows the extraction and visualisation of biological protein-protein interfaces from an annotated subset of crystallographic structures of proteins. This tool is focused on the description of protein-protein interfaces corresponding to well-identified classes of protein assemblies. It permits the representation of any of these protein-protein assemblies (duplex) and their interfaces as well as the export of the corresponding molecular structures under a flexible format, which is an extension of the PDBML. Moreover, PPIDD facilitates the construction of subsets of interfaces presenting user-specified common characteristics, to enhance the understanding of the determinants of specific protein-protein interactions.

Energetics of protein homodimerization: Effects of water sequestering on the formation of β-lactoglobulin dimer

Transient protein-protein interactions are functionally relevant as a control mechanism in a variety of biological processes. Analysis of the 3D structure of protein-protein complexes indicates that water molecules trapped at the interface are very common; however, their role in the stability and specificity of protein homodimer interactions has been not addressed yet. To provide new insights into the energetic bases that govern the formation of highly hydrated interfaces, the dissociation process of bovine beta lg variant A at a neutral pH was characterized here thermodynamically by conducting dilution experiments with an isothermal titration calorimeter. Association was enthalpically driven throughout the temperature range spanned. DeltaH and deltaC(p) were significantly more negative than estimates based on surface area changes, suggesting the occurrence of effects additional to the dehydration of the contact surfaces between subunits. Near-UV CD spectra proved to be independent of protein concentration, indicating a rigid body-like association. Furthermore, the process proved not to be coupled to significant changes in the protonation state of ionizable groups or counterion exchange. In contrast, both osmotic stress experiments and a computational analysis of the dimer's 3D structure indicated that a large number of water molecules are incorporated into the interface upon association. Numerical estimates considering the contributions of interface area desolvation and water immobilization accounted satisfactorily for the experimental deltaC(p). Thus, our study highlights the importance of explicitly considering the effects of water sequestering to perform a proper quantitative analysis of the formation of homodimers with highly hydrated interfaces.

Protein crystallization in drug design: towards a rational approach

X-ray crystallography is the method of choice for the detailed characterization of stereochemistry of interactions of drug leads and potential chemotherapeutics with their protein targets. The resulting atomic models allow for rational enhancement of the lead properties and consequently for the design of high-affinity inhibitors. However, a major bottleneck of the technique is the requirement for the protein and its complexes to yield high quality single crystals. Furthermore, it is highly desirable that such crystals diffract to high resolution, preferably ≥ 1.2 Å, revealing the structures in atomic detail. Unfortunately, only a small portion of proteins readily crystallize in that fashion. New approaches are being developed to circumvent this problem. One proposed option includes rational protein surface engineering to systematically improve the crystallizability of the protein. This is accomplished by creating surface patches readily mediating weak, but specific, intermolecular interactions that take on the role of crystal contacts during nucleation and crystal growth phase.

Structural adaptation to low temperatures − analysis of the subunit interface of oligomeric psychrophilic enzymes

Enzymes from psychrophiles show higher catalytic efficiency in the 0-20 degrees C temperature range and often lower thermostability in comparison with meso/thermophilic homologs. Physical and chemical characterization of these enzymes is currently underway in order to understand the molecular basis of cold adaptation. Psychrophilic enzymes are often characterized by higher flexibility, which allows for better interaction with substrates, and by a lower activation energy requirement in comparison with meso/thermophilic counterparts. In their tertiary structure, psychrophilic enzymes present fewer stabilizing interactions, longer and more hydrophilic loops, higher glycine content, and lower proline and arginine content. In this study, a comparative analysis of the structural characteristics of the interfaces between oligomeric psychrophilic enzyme subunits was carried out. Crystallographic structures of oligomeric psychrophilic enzymes, and their meso/thermophilic homologs belonging to five different protein families, were retrieved from the Protein Data Bank. The following structural parameters were calculated: overall and core interface area, characterization of polar/apolar contributions to the interface, hydrophobic contact area, quantity of ion pairs and hydrogen bonds between monomers, internal area and total volume of non-solvent-exposed cavities at the interface, and average packing of interface residues. These properties were compared to those of meso/thermophilic enzymes. The results were analyzed using Student's t-test. The most significant differences between psychrophilic and mesophilic proteins were found in the number of ion pairs and hydrogen bonds, and in the apolarity of their subunit interface. Interestingly, the number of ion pairs at the interface shows an opposite adaptation to those occurring at the monomer core and surface.

Secondary structure based analysis and classification of biological interfaces: identification of binding motifs in protein–protein interactions

MOTIVATION

The increasing amount of data on protein-protein interaction needs to be rationalized for deriving guidelines for the alteration or design of an interface between two proteins.

RESULTS

We present a detaild structural analysis and comparison of homo- versus heterodimeric protein-protein interfaces. Regular secondary structures (helices and strands) are the main components of the former, whereas non-regular structures (turns, loops, etc.) frequently mediate interactions in the latter. Interface helices get longer with increasing interface area, but only in heterocomplexes. On average, the homodimers have longer helical segments and prominent helix-helix pairs. There is a surprising distinction in the relative orientation of interface helices, with a tendency for aligned packing in homodimers and a clear preference for packing at 90 degrees in heterodimers. Arg and the aromatic residues have a higher preference to occur in all secondary structural elements (SSEs) in the interface. Based on the dominant SSE, the interfaces have been grouped into four classes: alpha, beta, alphabeta and non-regular. Identity between protein and interface classes is the maximum for alpha proteins, but rather mediocre for the other protein classes. The interface classes of the two chains forming a heterodimer are often dissimilar. Eleven binding motifs can capture the prominent architectural features of most of the interfaces.

PIER: protein interface recognition for structural proteomics

Recent advances in structural proteomics call for development of fast and reliable automatic methods for prediction of functional surfaces of proteins with known three-dimensional structure, including binding sites for known and unknown protein partners as well as oligomerization interfaces. Despite significant progress the problem is still far from being solved. Most existing methods rely, at least partially, on evolutionary information from multiple sequence alignments projected on protein surface. The common drawback of such methods is their limited applicability to the proteins with a sparse set of sequential homologs, as well as inability to detect interfaces in evolutionary variable regions. In this study, the authors developed an improved method for predicting interfaces from a single protein structure, which is based on local statistical properties of the protein surface derived at the level of atomic groups. The proposed Protein IntErface Recognition (PIER) method achieved the overall precision of 60% at the recall threshold of 50% at the residue level on a diverse benchmark of 490 homodimeric, 62 heterodimeric, and 196 transient interfaces (compared with 25% precision at 50% recall expected from random residue function assignment). For 70% of proteins in the benchmark, the binding patch residues were successfully detected with precision exceeding 50% at 50% recall. The calculation only took seconds for an average 300-residue protein. The authors demonstrated that adding the evolutionary conservation signal only marginally influenced the overall prediction performance on the benchmark; moreover, for certain classes of proteins, using this signal actually resulted in a deteriorated prediction. Thorough benchmarking using other datasets from literature showed that PIER yielded improved performance as compared with several alignment-free or alignment-dependent predictions. The accuracy, efficiency, and dependence on structure alone make PIER a suitable tool for automated high-throughput annotation of protein structures emerging from structural proteomics projects.

Inference of macromolecular assemblies from crystalline state

We discuss basic physical-chemical principles underlying the formation of stable macromolecular complexes, which in many cases are likely to be the biological units performing a certain physiological function. We also consider available theoretical approaches to the calculation of macromolecular affinity and entropy of complexation. The latter is shown to play an important role and make a major effect on complex size and symmetry. We develop a new method, based on chemical thermodynamics, for automatic detection of macromolecular assemblies in the Protein Data Bank (PDB) entries that are the results of X-ray diffraction experiments. As found, biological units may be recovered at 80-90% success rate, which makes X-ray crystallography an important source of experimental data on macromolecular complexes and protein-protein interactions. The method is implemented as a public WWW service.

Character and evolution of protein-protein interfaces

Protein-protein interactions create the macromolecular assemblies and sequential signaling pathways essential for cell function. Their number far exceeds the number of proteins themselves and their experimental characterization, while improving, remains relatively slow. For these reasons, novel computational methods have important roles to play in understanding the physical basis of protein interactions, and in constraining the molecular basis of their specificity. This paper discusses methods based on multiple sequence alignments of protein homologues and phylogenetic trees.

Crosstalk between the subunits of the homodimeric enzyme triosephosphate isomerase

Homodimeric triosephosphate isomerase (TIM) from Trypanosoma cruzi (TcTIM) and T. brucei (TbTIM) are markedly similar in amino acid sequence and three-dimensional structure. In their dimer interfaces, each monomer has a Cys15 that is surrounded by loop3 of the adjoining subunit. Perturbation of Cys15 by methylmethane thiosulfonate (MMTS) induces abolition of catalysis and structural changes. In the two TIMs, the structural arrangements of their Cys15 are almost identical. Nevertheless, the susceptibility of TcTIM to MMTS is nearly 100-fold higher than in TbTIM. To ascertain the extent to which the characteristics of the interface Cys depend on the dynamics of its own monomer or on those of the adjacent monomer, we studied MMTS action on mutants of TcTIM that had the interface residues of TbTIM, and hybrids that have only one interfacial Cys15 (C15ATcTIM-wild type TbTIM). We found that the solvent exposure of the interfacial Cys depends predominantly on the characteristics of the adjoining monomer. The maximal inhibition of activity induced by perturbation of the sole interface Cys in the C15ATcTIM-TbTIM hybrid is around 60%. Hybrids formed with C15ATcTIM monomers and catalytically inert TbTIM monomers (E168DTbTIM) were also studied. Their activity drops by nearly 50% when the only interfacial Cys is perturbed. These results in conjunction with those on C15ATcTIM-wild type TbTIM hybrid indicate that about half of the activity of each monomer depends on the integrity of each of the two Cys15-loop3 portions of the interface. This could be another reason of why TIM is an obligatory dimer.

Peptide segments in protein-protein interfaces

An important component of functional genomics involves the understanding of protein association. The interfaces resulting from protein-protein interactions - (i) specific, as represented by the homodimeric quaternary structures and the complexes formed by two independently occurring protein components, and (ii) non-specific, as observed in the crystal lattice of monomeric proteins - have been analysed on the basis of the length and the number of peptide segments. In 1000 A2 of the interface area, contributed by a polypeptide chain, there would be 3.4 segments in homodimers, 5.6 in complexes and 6.3 in crystal contacts. Concomitantly, the segments are the longest (with 8.7 interface residues) in homodimers. Core segments (likely to contribute more towards binding) are more in number in homodimers (1.7) than in crystal contacts (0.5), and this number can be used as one of the parameters to distinguish between the two types of interfaces. Dominant segments involved in specific interactions, along with their secondary structural features, are enumerated.

Comprehensive statistical analysis of residues interaction specificity at protein-protein interfaces

We calculated interchain contacts on the atomic level for nonredundant set of 4602 protein-protein interfaces using an unbiased Voronoi-Delaune tessellation method, and made 20x20 residue contact matrixes both for homodimers and heterocomplexes. The area of contacts and the distance distribution for these contacts were calculated on both the residue and the atomic levels. We analyzed residue area distribution and showed the existence of two types of interresidue contacts: stochastic and specific. We also derived formulas describing the distribution of contact area for stochastic and specific interactions in parametric form. Maximum pairing preference index was found for Cys-Cys contacts and for oppositely charged interactions. A significant difference in residue contacts was observed between homodimers and heterocomplexes. Interfaces in homodimers were enriched with contacts between residues of the same type due to the effects of structure symmetry.

A Dissection of the Protein–Protein Interfaces in Icosahedral Virus Capsids

We selected 49 icosahedral virus capsids whose crystal structures are reported in the Protein Data Bank. They belong to the T=1, T=3, pseudo T=3 and other lattice types. We identified in them 779 unique interfaces between pairs of subunits, all repeated by icosahedral symmetry. We analyzed the geometric and physical chemical properties of these interfaces and compared with interfaces in protein-protein complexes and homodimeric proteins, and with crystal packing contacts. The capsids contain one to 16 subunits implicated in three to 66 unique interfaces. Each subunit loses 40-60% of its accessible surface in contacts with an average of 8.5 neighbors. Many of the interfaces are very large with a buried surface area (BSA) that can exceed 10,000 A(2), yet 39% are small with a BSA<800 A(2) comparable to crystal packing contacts. Pairwise capsid interfaces overlap, so that one-third of the residues are part of more than one interface. Those with a BSA>800 A(2) resemble homodimer interfaces in their chemical composition. Relative to the protein surface, they are non-polar, enriched in aliphatic residues and depleted of charged residues, but not of neutral polar residues. They contain one H-bond per about 200 A(2) BSA. Small capsid interfaces (BSA<800 A(2)) are only slightly more polar. They have a similar amino acid composition, but they bury fewer atoms and contain fewer H-bonds for their size. Geometric parameters that estimate the quality of the atomic packing suggest that the small capsid interfaces are loosely packed like crystal packing contacts, whereas the larger interfaces are close-packed as in protein-protein complexes and homodimers. We discuss implications of these findings on the mechanism of capsid assembly, assuming that the larger interfaces form first to yield stable oligomeric species (capsomeres), and that medium-size interfaces allow the stepwise addition of capsomeres to build larger intermediates.

X-ray crystallographic analysis of the sulfur carrier protein SoxY from Chlorobium limicola f. thiosulfatophilum reveals a tetrameric structure

Dissimilatory oxidation of thiosulfate in the green sulfur bacterium Chlorobium limicola f. thiosulfatophilum is carried out by the ubiquitous sulfur-oxidizing (Sox) multi-enzyme system. In this system, SoxY plays a key role, functioning as the sulfur substrate-binding protein that offers its sulfur substrate, which is covalently bound to a conserved C-terminal cysteine, to another oxidizing Sox enzyme. Here, we report the crystal structures of a stand-alone SoxY protein of C. limicola f. thiosulfatophilum, solved at 2.15 A and 2.40 A resolution using X-ray diffraction data collected at 100 K and room temperature, respectively. The structure reveals a monomeric Ig-like protein, with an N-terminal alpha-helix, that oligomerizes into a tetramer via conserved contact regions between the monomers. The tetramer can be described as a dimer of dimers that exhibits one large hydrophobic contact region in each dimer and two small hydrophilic interface patches in the tetramer. At the tetramer interface patch, two conserved redox-active C-terminal cysteines form an intersubunit disulfide bridge. Intriguingly, SoxY exhibits a dimer/tetramer equilibrium that is dependent on the redox state of the cysteines and on the type of sulfur substrate component bound to them. Taken together, the dimer/tetramer equilibrium, the specific interactions between the subunits in the tetramer, and the significant conservation level of the interfaces strongly indicate that these SoxY oligomers are biologically relevant.

Structure-based secondary structure-independent approach to design protein ligands: Application to the design of Kv1.2 potassium channel blockers

We have developed a structure-based approach to the design of protein ligands. This approach is based on the transfer of a functional binding motif of amino acids, often referred as to the "hot spot", on a host protein able to reproduce the functional topology of these residues. The scaffolds were identified by a systematic in silico search in the Protein Data Bank for proteins possessing a group of residues in a topology similar to that adopted by the functional motif in a reference ligand of known 3D structure. In contrast to previously reported studies, this search is independent of the particular secondary structure supporting the functional motif. To take into account the global properties of the host protein, two additional criteria were taken into account in the selection process: (1) Only those scaffolds sterically compatible with the positioning of the functional motif as observed in a reference complex model were retained. (2) Host proteins displaying electrostatic potentials, in the region of the transferred functional motif, similar to that of the reference ligand were selected. This approach was applied to the development of protein ligands of the Kv1.2 channel using BgK, a small protein isolated from the sea anemone Bunodosoma granulifera, as the reference ligand. Four proteins obtained by this approach were produced for experimental evaluation. The X-ray structure of one of these proteins was determined to check for similarity of the transferred functional motif with the structure it adopts in the reference ligand. Three of these protein ligands bind the Kv1.2 channel with inhibition constants of 0.5, 1.5, and 1.6 microM. Several mutants of these designed protein ligands gave binding results consistent with the presumed binding mode. These results show that protein ligands can be designed by transferring a binding motif on a protein host selected to reproduce the functional topology of this motif, irrespective to the secondary structure supporting the functional motif, if the host protein possesses steric and electrostatic properties compatible with the binding to the target. This result opens the way to the design of protein ligands by taking advantage of the considerable structural repertoire of the Protein Data Bank.

Interaction geometry involving planar groups in protein-protein interfaces

The geometry of interactions of planar residues is nonrandom in protein tertiary structures and gives rise to conventional, as well as nonconventional (X--H...pi, X--H...O, where X = C, N, or O) hydrogen bonds. Whether a similar geometry is maintained when the interaction is across the protein-protein interface is addressed here. The relative geometries of interactions involving planar residues, and the percentage of contacts giving rise to different types of hydrogen bonds are quite similar in protein structures and the biological interfaces formed by protein chains in homodimers and protein-protein heterocomplexes--thus pointing to the similarity of chemical interactions that occurs during protein folding and binding. However, the percentage is considerably smaller in the nonspecific and nonphysiological interfaces that are formed in crystal lattices of monomeric proteins. The C--H...O interaction linking the aromatic and the peptide groups is quite common in protein structures as well as the three types of interfaces. However, as the interfaces formed by crystal contacts are depleted in aromatic residues, the weaker hydrogen bond interactions would contribute less toward their stability.

From protein-protein complexes to interactomics

Protein-protein interactions (or PPIs) are key elements for the normal functioning of a living cell. A large description of the protein interactomics field is given in this review where different aspects will be discussed. We first give an introduction of the different large scale experimental approaches from yeast two-hybrid to mass spectrometry used to discover PPIs and build protein interaction maps. Single PPI validation techniques such as co-immunoprecipitation or fluorescence methods are then presented as they are more and more integrated in global PPI discovery strategy. Data from different experimental sets are compared and an assessment of the different large scale technologies is presented. Bioinformatics tools can also predict with a good accuracy PPIs in silico, PPIs databases are now numerous and topological analysis has led to interesting insights into the nature of network connection. Finally, PPI, as an association of two proteins, has been structurally characterized for many protein complexes and is largely discussed throughout existing examples. The results obtained so far already provide the biologist with a large set of structured data from which knowledge on pathways and associated protein function can be extracted.

Macromolecular recognition in the Protein Data Bank

Crystal structures deposited in the Protein Data Bank illustrate the diversity of biological macromolecular recognition: transient interactions in protein-protein and protein-DNA complexes and permanent assemblies in homodimeric proteins. The geometric and physical chemical properties of the macromolecular interfaces that may govern the stability and specificity of recognition are explored in complexes and homodimers compared with crystal-packing interactions. It is found that crystal-packing interfaces are usually much smaller; they bury fewer atoms and are less tightly packed than in specific assemblies. Standard-size interfaces burying 1200-2000 A2 of protein surface occur in protease-inhibitor and antigen-antibody complexes that assemble with little or no conformation changes. Short-lived electron-transfer complexes have small interfaces; the larger size of the interfaces observed in complexes involved in signal transduction and homodimers correlates with the presence of conformation changes, often implicated in biological function. Results of the CAPRI (critical assessment of predicted interactions) blind prediction experiment show that docking algorithms efficiently and accurately predict the mode of assembly of proteins that do not change conformation when they associate. They perform less well in the presence of large conformation changes and the experiment stimulates the development of novel procedures that can handle such changes.

Theoretical model of the three-dimensional structure of a disease resistance gene homolog encoding resistance protein in Vigna mungo

Plant disease resistance (R) genes, the key players of innate immunity system in plants encode 'R' proteins. 'R' protein recognizes product of avirulance gene from the pathogen and activate downstream signaling responses leading to disease resistance. No three dimensional (3D) structural information of any 'R' proteins is available as yet. We have reported a 'R' gene homolog, the 'VMYR1', encoding 'R' protein in Vigna mungo. Here, we describe the homology modeling of the 'VMYR1' protein. The model was created by using the 3D structure of an ATP-binding cassette transporter protein from Vibrio cholerae as a template. The strategy for homology modeling was based on the high structural conservation in the superfamily of P-loop containing nucleoside triphosphate hydrolase in which target and template proteins belong. This is the first report of theoretical model structure of any 'R' proteins.

Structural similarity enhances interaction propensity of proteins

We study statistical properties of interacting protein-like surfaces and predict two strong, related effects: (i) statistically enhanced self-attraction of proteins; (ii) statistically enhanced attraction of proteins with similar structures. The effects originate in the fact that the probability to find a pattern self-match between two identical, even randomly organized interacting protein surfaces is always higher compared with the probability for a pattern match between two different, promiscuous protein surfaces. This theoretical finding explains statistical prevalence of homodimers in protein-protein interaction networks reported earlier. Further, our findings are confirmed by the analysis of curated database of protein complexes that showed highly statistically significant overrepresentation of dimers formed by structurally similar proteins with highly divergent sequences ("superfamily heterodimers"). We suggest that promiscuous homodimeric interactions pose strong competitive interactions for heterodimers evolved from homodimers. Such evolutionary bottleneck is overcome using the negative design evolutionary pressure applied against promiscuous homodimer formation. This is achieved through the formation of highly specific contacts formed by charged residues as demonstrated both in model and real superfamily heterodimers.

The many faces of protein-protein interactions: A compendium of interface geometry

A systematic classification of protein-protein interfaces is a valuable resource for understanding the principles of molecular recognition and for modelling protein complexes. Here, we present a classification of domain interfaces according to their geometry. Our new algorithm uses a hybrid approach of both sequential and structural features. The accuracy is evaluated on a hand-curated dataset of 416 interfaces. Our hybrid procedure achieves 83% precision and 95% recall, which improves the earlier sequence-based method by 5% on both terms. We classify virtually all domain interfaces of known structure, which results in nearly 6,000 distinct types of interfaces. In 40% of the cases, the interacting domain families associate in multiple orientations, suggesting that all the possible binding orientations need to be explored for modelling multidomain proteins and protein complexes. In general, hub proteins are shown to use distinct surface regions (multiple faces) for interactions with different partners. Our classification provides a convenient framework to query genuine gene fusion, which conserves binding orientation in both fused and separate forms. The result suggests that the binding orientations are not conserved in at least one-third of the gene fusion cases detected by a conventional sequence similarity search. We show that any evolutionary analysis on interfaces can be skewed by multiple binding orientations and multiple interaction partners. The taxonomic distribution of interface types suggests that ancient interfaces common to the three major kingdoms of life are enriched by symmetric homodimers. The classification results are online at http://www.scoppi.org.

Intramolecular surface contacts contain information about protein-protein interface regions

MOTIVATION

Some amino acids clearly show preferences over others in protein-protein interfaces. These preferences, or so-called interface propensities can be used for a priori interface prediction. We investigated whether the prediction accuracy could be improved by considering not single but pairs of residues in an interface. Here we present the first systematic analysis of intramolecular surface contacts in interface prediction.

RESULTS

We show that preferences do exist for contacts within and around an interface region within one molecule: specific pairs of amino acids are more often occurring than others. Using intramolecular contact propensities in a blind test, higher average scores were assigned to interface residues than to non-interface residues. This effect persisted as small but significant when the contact propensities were corrected to eliminate the influence of single amino acid interface propensity. This indicates that intramolecular contact propensities may replace interface propensities in protein-protein interface prediction.

AVAILABILITY

The source code is available on request from the authors.

ProFace: a server for the analysis of the physicochemical features of protein-protein interfaces

Background

Molecular recognition is all pervasive in biology. Protein molecules are involved in enzyme regulation, immune response, signal transduction, oligomer assembly, etc. Delineation of physical and chemical features of the interface formed by protein-protein association would allow us to better understand protein interaction networks on one hand, and to design molecules that can engage a given interface and thereby control protein function on the other hand.

Results

ProFace is a suite of programs that uses a file, containing atomic coordinates of a multi-chain molecule, as input and analyzes the interface between any two or more subunits. The interface residues are shown segregated into spatial patches (if such a clustering is possible based on an input threshold distance) and/or core and rim regions. A number of physicochemical parameters defining the interface is tabulated. Among the different output files, one contains the list of interacting residues across the interface. Results can be used to infer if a particular interface belongs to a homodimeric molecule.

Conclusion

A web-server, ProFace (available at ) has been developed for dissecting protein-protein interfaces and deriving various physicochemical parameters.

The origins and evolution of functional modules: lessons from protein complexes

Modularity is an attribute of a system that can be decomposed into a set of cohesive entities that are loosely coupled. Many cellular networks can be decomposed into functional modules-each functionally separable from the other modules. The protein complexes in physical protein interaction networks are a good example of this, and here we focus on their origins and evolution. We investigate the emergence of protein complexes and physical interactions between proteins by duplication, and review other mechanisms. We dissect the dataset of protein complexes of known three-dimensional structure, and show that roughly 90% of these complexes contain contacts between identical proteins within the same complex. Proteins that are shared across different complexes occur frequently, and they tend to be essential genes more often than members of a single protein complex. We also provide a perspective on the evolutionary mechanisms driving the growth of other modular cellular networks such as transcriptional regulatory and metabolic networks.

Drugs Targeting Protein–Protein Interactions

Most biological processes involve permanent and nonpermanent interactions between different proteins, and many protein complexes play a key role in various human diseases. Therefore, molecules that prevent the formation of these protein complexes could be valuable new therapeutic agents to treat these diseases. Protein interfaces have not evolved to bind low-molecular-weight molecules, as is the case with enzyme catalytic sites. It is therefore difficult to identify small compounds that inhibit protein-protein interactions. However, there is considerable diversity in the structure of protein interfaces, some of which may be more attractive than others for medicinal chemistry. One of the main challenges in drug discovery is to identify these interfaces and to exploit their properties to make marketable drugs. Herein, the properties of protein interfaces are discussed in light of their use as drug targets.

Primary and secondary interactions between CK2alpha and CK2beta lead to ring-like structures in the crystals of the CK2 holoenzyme

Protein kinase CK2 predominantly exists as a heterotetrameric holoenyzme consisting of two catalytic subunits (CK2alpha) and two non-catalytic subunits (CK2beta). Early investigations which we review here had revealed the presence of two types of contacts between CK2alpha and CK2beta: a primary interaction responsible for the stability of the CK2 holoenzyme and stimulatory for the catalytic activity, and a secondary interaction which is inhibitory and in which the acidic loop of CK2beta associates with the basic stretch and the (p+1)-loop of CK2alpha. At the end of the last decade both types of interactions were assumed to occur within the same tetrameric complex. The CK2 holoenyzme structure, however, suggested that the secondary interactions must happen between different CK2 tetramers. Such a behaviour should lead to higher-ordered aggregates consistent with several previous reports about a distinct aggregation propensity of CK2. We demonstrate here that in the CK2 holoenzyme crystals contacts between different CK2 tetramers exists which provide structural details of the secondary CK2alpha/CK2beta interactions. These mainly ionic interactions lead to trimeric rings of CK2 holoenzymes in the crystal. In these rings each CK2 tetramer possesses one CK2alpha subunit open for substrate binding and another one whose active site is blocked by a secondary contact with CK2beta from a neighbouring tetramer. This observation fits to previous findings that salt-sensitive ring-like aggregates of CK2 holoenzymes can exist which possess significant catalytic activity. Furthermore it suggests that earlier ideas about a regulatory role of the enzyme's aggregation propensity may be worth to be revitalised.

Structural basis of reduction-dependent activation of human cystatin F

Cystatins are important natural cysteine protease inhibitors targeting primarily papain-like cysteine proteases, including cathepsins and parasitic proteases like cruzipain, but also mammalian asparaginyl endopeptidase. Mammalian cystatin F, which is expressed almost exclusively in hematopoietic cells and accumulates in lysosome-like organelles, has been implicated in the regulation of antigen presentation and other immune processes. It is an unusual cystatin superfamily member with a redox-regulated activation mechanism and a restricted specificity profile. We describe the 2.1A crystal structure of human cystatin F in its dimeric "off" state. The two monomers interact in a fashion not seen before for cystatins or cystatin-like proteins that is crucially dependent on an unusual intermolecular disulfide bridge, suggesting how reduction leads to monomer formation and activation. Strikingly, core sugars for one of the two N-linked glycosylation sites of cystatin F are well ordered, and their conformation and interactions with the protein indicate that this unique feature of cystatin F may modulate its inhibitory properties, in particular its reduced affinity toward asparaginyl endopeptidase compared with other cystatins.

Crystal structure of Bacillus subtilis TrmB, the tRNA (m7G46) methyltransferase

The structure of Bacillus subtilis TrmB (BsTrmB), the tRNA (m⁷G46) methyltransferase, was determined at a resolution of 2.1 Å. This is the first structure of a member of the TrmB family to be determined by X-ray crystallography. It reveals a unique variant of the Rossmann-fold methyltransferase (RFM) structure, with the N-terminal helix folded on the opposite site of the catalytic domain. The architecture of the active site and a computational docking model of BsTrmB in complex with the methyl group donor S-adenosyl-l-methionine and the tRNA substrate provide an explanation for results from mutagenesis studies of an orthologous enzyme from Escherichia coli (EcTrmB). However, unlike EcTrmB, BsTrmB is shown here to be dimeric both in the crystal and in solution. The dimer interface has a hydrophobic core and buries a potassium ion and five water molecules. The evolutionary analysis of the putative interface residues in the TrmB family suggests that homodimerization may be a specific feature of TrmBs from Bacilli, which may represent an early stage of evolution to an obligatory dimer.

Identification of hot spot residues at protein-protein interface

It is known that binding free energy of protein-protein interaction is mainly contributed by hot spot (high energy) interface residues. Here, we investigate the characteristics of hot spots by examining inter-atomic sidechain-sidechain interactions using a dataset of 296 alanine-mutated interface residues. Results show that hot spots participate in strong and energetically favorable sidechain-sidechain interactions. Subsequently, we describe a novel, yet simple 'hot spot' prediction model with an accuracy that is similar to many available approaches. The model is also shown to efficiently distinguish specific protein-protein interactions from non-specific interactions.

Survey of the geometric association of domain-domain interfaces

Considering the limited success of the most sophisticated docking methods available and the amount of computation required for systematic docking, cataloging all the known interfaces may be an alternative basis for the prediction of protein tertiary and quaternary structures. We classify domain interfaces according to the geometry of domain-domain association. By applying a simple and efficient method called "interface tag clustering," more than 4,000 distinct types of domain interfaces are collected from Protein Quaternary Structure Server and Protein Data Bank. Given a pair of interacting domains, we define "face" as the set of interacting residues in each single domain and the pair of interacting faces as an "interface." We investigate how the geometry of interfaces relates to a network of interacting protein families, such as how many different binding orientations are possible between two families or whether a family uses distinct surfaces or the same surface when the family has diverse interaction partners from various families. We show there are, on average, 1.2-1.9 different types of interfaces between interacting domains and a significant number of family pairs associate in multiple orientations. In general, a family tends to use distinct faces for each partner when the family has diverse interaction partners. Each face is highly specific to its interaction partner and the binding orientation. The relative positions of interface residues are generally well conserved within the same type of interface even between remote homologs. The classification result is available at http://www.biotec.tu-dresden.de/~wkim/supplement.

Assessing methods for identifying pair-wise atomic contacts across binding interfaces

An essential step in understanding the molecular basis of protein-protein interactions is the accurate identification of inter-protein contacts. We evaluate a number of common methods used in analyzing protein-protein interfaces: a Voronoi polyhedra-based approach, changes in solvent accessible surface area (DeltaSASA) and various radial cutoffs (closest atom, Cbeta, and centroid). First, we compared the Voronoi polyhedra-based analysis to the DeltaSASA and show that using Voronoi polyhedra finds knob-in-hole contacts. To assess the accuracy between the Voronoi polyhedra-based approach and the various radial cutoff methods, two sets of data were used: a small set of 75 experimental mutants and a larger one of 592 structures of protein-protein interfaces. In an assessment using the small set, the Voronoi polyhedra-based methods, a solvent accessible surface area method, and the closest atom radial method identified 100% of the direct contacts defined by mutagenesis data, but only the Voronoi polyhedra-based method found no false positives. The other radial methods were not able to find all of the direct contacts even using a cutoff of 9A. With the larger set of structures, we compared the overall number contacts using the Voronoi polyhedra-based method as a standard. All the radial methods using a 6-A cutoff identified more interactions, but these putative contacts included many false positives as well as missed many false negatives. While radial cutoffs are quicker to calculate as well as to implement, this result highlights why radial cutoff methods do not have the proper resolution to detail the non-homogeneous packing within protein interfaces, and suggests an inappropriate bias in pair-wise contact potentials. Of the radial cutoff methods, using the closest atom approach exhibits the best approximation to the more intensive Voronoi calculation. Our version of the Voronoi polyhedra-based method QContacts is available at .

Energetics of outer membrane phospholipase A (OMPLA) dimerization

Outer membrane phospholipase A (OMPLA) is a widely conserved transmembrane enzyme found in Gram-negative bacteria, and it is implicated in the virulence of a number of pathogenic organisms. The regulation of the protein's phospholipase activity is not well understood despite the existence of a number of high resolution structures. Previous biochemical studies have demonstrated that dimerization of OMPLA is a prerequisite for its phospholipase activity, and it has been shown in vitro that this dimerization is dependent on calcium and substrate binding. Therefore, to fully understand the regulation of OMPLA, it is necessary to understand the stability of the protein dimer and the extent to which it is influenced by its effector molecules. We have used sedimentation equilibrium analytical ultracentrifugation to dissect the energetics of Escherichia coli OMPLA dimerization in detergent micelles. We find that calcium contributes relatively little stability to the dimer, while interactions with the substrate acyl chain are the predominant force in stabilizing the dimeric conformation of the enzyme. The resulting thermodynamic cycle suggests that interactions between effector molecules are additive. These energetic measurements not only provide insight into the activation of OMPLA, but they also represent the first quantitative investigation of the association energetics of a transmembrane beta-barrel. This thermodynamic study allows us to begin to address the differences between protein-protein interfaces in transmembrane proteins with a helical fold to those of a beta-barrel fold and to more fully understand the forces involved in membrane protein interactions.

Optimal docking area: a new method for predicting protein-protein interaction sites

Understanding energetics and mechanism of protein-protein association remains one of the biggest theoretical problems in structural biology. It is assumed that desolvation must play an essential role during the association process, and indeed protein-protein interfaces in obligate complexes have been found to be highly hydrophobic. However, the identification of protein interaction sites from surface analysis of proteins involved in non-obligate protein-protein complexes is more challenging. Here we present Optimal Docking Area (ODA), a new fast and accurate method of analyzing a protein surface in search of areas with favorable energy change when buried upon protein-protein association. The method identifies continuous surface patches with optimal docking desolvation energy based on atomic solvation parameters adjusted for protein-protein docking. The procedure has been validated on the unbound structures of a total of 66 non-homologous proteins involved in non-obligate protein-protein hetero-complexes of known structure. Optimal docking areas with significant low-docking surface energy were found in around half of the proteins. The 'ODA hot spots' detected in X-ray unbound structures were correctly located in the known protein-protein binding sites in 80% of the cases. The role of these low-surface-energy areas during complex formation is discussed. Burial of these regions during protein-protein association may favor the complexed configurations with near-native interfaces but otherwise arbitrary orientations, thus driving the formation of an encounter complex. The patch prediction procedure is freely accessible at http://www.molsoft.com/oda and can be easily scaled up for predictions in structural proteomics.

Insights into the quaternary association of proteins through structure graphs: a case study of lectins

The unique three-dimensional structure of both monomeric and oligomeric proteins is encoded in their sequence. The biological functions of proteins are dependent on their tertiary and quaternary structures, and hence it is important to understand the determinants of quaternary association in proteins. Although a large number of investigations have been carried out in this direction, the underlying principles of protein oligomerization are yet to be completely understood. Recently, new insights into this problem have been gained from the analysis of structure graphs of proteins belonging to the legume lectin family. The legume lectins are an interesting family of proteins with very similar tertiary structures but varied quaternary structures. Hence they have become a very good model with which to analyse the role of primary structures in determining the modes of quaternary association. The present review summarizes the results of a legume lectin study as well as those obtained from a similar analysis carried out here on the animal lectins, namely galectins, pentraxins, calnexin, calreticulin and rhesus rotavirus Vp4 sialic-acid-binding domain. The lectin structure graphs have been used to obtain clusters of non-covalently interacting amino acid residues at the intersubunit interfaces. The present study, performed along with traditional sequence alignment methods, has provided the signature sequence motifs for different kinds of quaternary association seen in lectins. Furthermore, the network representation of the lectin oligomers has enabled us to detect the residues which make extensive interactions ('hubs') across the oligomeric interfaces that can be targetted for interface-destabilizing mutations. The present review also provides an overview of the methodology involved in representing oligomeric protein structures as connected networks of amino acid residues. Further, it illustrates the potential of such a representation in elucidating the structural determinants of protein-protein association in general and will be of significance to protein chemists and structural biologists.

X-ray structures of free and leupeptin-complexed human alphaI-tryptase mutants: indication for an alpha-->beta-tryptase transition

Tryptases alpha and beta are trypsin-like serine proteinases expressed in large amounts by mast cells. Beta-tryptase is a tetramer that has enzymatic activity, but requires heparin binding to maintain functional and structural stability, whereas alpha-tryptase has little, if any, enzymatic activity but is a stable tetramer in the absence of heparin. As shown previously, these differences can be mainly attributed to the different conformations of the 214-220 segment. Interestingly, the replacement of Asp216 by Gly, which is present in beta-tryptase, results in enzymatically active but less stable alpha-tryptase mutants. We have solved the crystal structures of both the single (D216G) and the double (K192Q/D216G) mutant forms of recombinant human alphaI-tryptase in complex with the peptide inhibitor leupeptin, as well as the structure of the non-inhibited single mutant. The inhibited mutants exhibited an open functional substrate binding site, while in the absence of an inhibitor, the open (beta-tryptase-like) and the closed (alpha-tryptase-like) conformations were present simultaneously. This shows that both forms are in a two-state equilibrium, which is influenced by the residues in the vicinity of the active site and by inhibitor/substrate binding. Novel insights regarding the observed stability differences as well as a potential proteolytic activity of wild-type alpha-tryptase, which may possess a cryptic active site, are discussed.

Oligomeric protein structure networks: insights into protein-protein interactions

BACKGROUND

Protein-protein association is essential for a variety of cellular processes and hence a large number of investigations are being carried out to understand the principles of protein-protein interactions. In this study, oligomeric protein structures are viewed from a network perspective to obtain new insights into protein association. Structure graphs of proteins have been constructed from a non-redundant set of protein oligomer crystal structures by considering amino acid residues as nodes and the edges are based on the strength of the non-covalent interactions between the residues. The analysis of such networks has been carried out in terms of amino acid clusters and hubs (highly connected residues) with special emphasis to protein interfaces.

RESULTS

A variety of interactions such as hydrogen bond, salt bridges, aromatic and hydrophobic interactions, which occur at the interfaces are identified in a consolidated manner as amino acid clusters at the interface, from this study. Moreover, the characterization of the highly connected hub-forming residues at the interfaces and their comparison with the hubs from the non-interface regions and the non-hubs in the interface regions show that there is a predominance of charged interactions at the interfaces. Further, strong and weak interfaces are identified on the basis of the interaction strength between amino acid residues and the sizes of the interface clusters, which also show that many protein interfaces are stronger than their monomeric protein cores. The interface strengths evaluated based on the interface clusters and hubs also correlate well with experimentally determined dissociation constants for known complexes. Finally, the interface hubs identified using the present method correlate very well with experimentally determined hotspots in the interfaces of protein complexes obtained from the Alanine Scanning Energetics database (ASEdb). A few predictions of interface hot spots have also been made based on the results obtained from this analysis, which await experimental verification.

CONCLUSION

The construction and analysis of oligomeric protein structure networks and their comparison with monomeric protein structure networks provide insights into protein association. Further, the interface hubs identified using the present method can be effective targets for interface de-stabilizing mutations. We believe this analysis will significantly enhance our knowledge of the principles behind protein association and also aid in protein design.

Homodimeric hexaprenyl pyrophosphate synthase from the thermoacidophilic crenarchaeon Sulfolobus solfataricus displays asymmetric subunit structures

Hexaprenyl pyrophosphate synthase (HexPPs) from Sulfolobus solfataricus catalyzes the synthesis of trans-C(30)-hexaprenyl pyrophosphate (HexPP) by reacting two isopentenyl pyrophosphate molecules with one geranylgeranyl pyrophosphate. The crystal structure of the homodimeric C(30)-HexPPs resembles those of other trans-prenyltransferases, including farnesyl pyrophosphate synthase (FPPs) and octaprenyl pyrophosphate synthase (OPPs). In both subunits, 10 core helices are arranged about a central active site cavity. Leu164 in the middle of the cavity controls the product chain length. Two protein conformers are observed in the S. solfataricus HexPPs structure, and the major difference between them occurs in the flexible region of residues 84 to 100. Several helices (alphaI, alphaJ, alphaK, and part of alphaH) and the associated loops have high-temperature factors in one monomer, which may be related to the domain motion that controls the entrance to the active site. Different side chain conformations of Trp136 in two HexPPs subunits result in weaker hydrophobic interactions at the dimer interface, in contrast to the symmetric pi-pi stacking interactions of aromatic side chains found in FPPs and OPPs. Finally, the three-conformer switched model may explain the catalytic process for HexPPs.

Conservation and relative importance of residues across protein-protein interfaces

A core region surrounded by a rim characterizes biological interfaces. We ascertain the importance of the core by showing the sequence entropies of the residues comprising the core to be smaller than those in the rim. Such a distinction is not seen in the 2-fold-related, nonphysiological interfaces formed in crystal lattices of monomeric proteins, thereby providing a procedure for characterizing the oligomeric state from crystal structures of protein molecules. This method is better than those that rely on the comparison of the sequence entropies in the interface and the rest of the protein surface, especially in cases where the surface harbors additional binding sites. To a good approximation there is a correlation between the accessible surface area lost because of complexation and DeltaDeltaG values obtained through alanine-scanning mutagenesis (26-38 cal per A(2) of the surface buried) for residues located in the core, a relationship that is not discernable for rim residues. If, however, a residue participates in hydrogen bonding across the interface, the extent of stabilization is 52 cal/mol per 1 A(2) of the nonpolar surface area buried by the residue. As opposed to an amino acid classification used earlier, an environment-based grouping of residues yields a better discrimination in the sequence entropy between the core and the rim.

Analysis of the quaternary structure of the MutL C-terminal domain

The dimeric DNA mismatch repair protein MutL has a key function in communicating mismatch recognition by MutS to downstream repair processes. Dimerization of MutL is mediated by the C-terminal domain, while activity of the protein is modulated by the ATP-dependent dimerization of the highly conserved N-terminal domain. Recently, a crystal structure analysis of the Escherichia coli MutL C-terminal dimerization domain has been reported and a model for the biological dimer was proposed. In this model, dimerization is mediated by the internal (In) subdomain comprising residues 475-569. Here, we report a computational analysis of all protein interfaces observed in the crystal structure and suggest that the biological dimer interface is formed by a hydrophobic surface patch of the external (Ex) subdomain (residues 432-474 and 570-615). Moreover, sequence analysis revealed that this surface patch is conserved among the MutL proteins. To test this hypothesis, single and double-cysteine variants of MutL were generated and tested for their ability to be cross-linked with chemical cross-linkers of various size. Finally, deletion of the C-terminal residues 605-615 abolished homodimerization. The biochemical data are fully compatible with a revised model for the biological dimer, which has important implications for understanding the heterodimerization of eukaryotic MutL homologues, modeling the MutL holoenzyme and predicting protein-protein interaction sites.

A protein domain interaction interface database: InterPare

BACKGROUND

Most proteins function by interacting with other molecules. Their interaction interfaces are highly conserved throughout evolution to avoid undesirable interactions that lead to fatal disorders in cells. Rational drug discovery includes computational methods to identify the interaction sites of lead compounds to the target molecules. Identifying and classifying protein interaction interfaces on a large scale can help researchers discover drug targets more efficiently.

DESCRIPTION

We introduce a large-scale protein domain interaction interface database called InterPare http://interpare.net. It contains both inter-chain (between chains) interfaces and intra-chain (within chain) interfaces. InterPare uses three methods to detect interfaces: 1) the geometric distance method for checking the distance between atoms that belong to different domains, 2) Accessible Surface Area (ASA), a method for detecting the buried region of a protein that is detached from a solvent when forming multimers or complexes, and 3) the Voronoi diagram, a computational geometry method that uses a mathematical definition of interface regions. InterPare includes visualization tools to display protein interior, surface, and interaction interfaces. It also provides statistics such as the amino acid propensities of queried protein according to its interior, surface, and interface region. The atom coordinates that belong to interface, surface, and interior regions can be downloaded from the website.

CONCLUSION

InterPare is an open and public database server for protein interaction interface information. It contains the large-scale interface data for proteins whose 3D-structures are known. As of November 2004, there were 10,583 (Geometric distance), 10,431 (ASA), and 11,010 (Voronoi diagram) entries in the Protein Data Bank (PDB) containing interfaces, according to the above three methods. In the case of the geometric distance method, there are 31,620 inter-chain domain-domain interaction interfaces and 12,758 intra-chain domain-domain interfaces.

Adaptation of a fast Fourier transform-based docking algorithm for protein design

Designing proteins with novel protein/protein binding properties can be achieved by combining the tools that have been developed independently for protein docking and protein design. We describe here the sequence-independent generation of protein dimer orientations by protein docking for use as scaffolds in protein sequence design algorithms. To dock monomers into sequence-independent dimer conformations, we use a reduced representation in which the side chains are approximated by spheres with atomic radii derived from known C2 symmetry-related homodimers. The interfaces of C2-related homodimers are usually more hydrophobic and protein core-like than the interfaces of heterodimers; we parameterize the radii for docking against this feature to capture and recreate the spatial characteristics of a hydrophobic interface. A fast Fourier transform-based geometric recognition algorithm is used for docking the reduced representation protein models. The resulting docking algorithm successfully predicted the wild-type homodimer orientations in 65 out of 121 dimer test cases. The success rate increases to approximately 70% for the subset of molecules with large surface area burial in the interface relative to their chain length. Forty-five of the predictions exhibited less than 1 A C(alpha) RMSD compared to the native X-ray structures. The reduced protein representation therefore appears to be a reasonable approximation and can be used to position protein backbones in plausible orientations for homodimer design.

Interresidue Contacts in Proteins and Protein−Protein Interfaces and Their Use in Characterizing the Homodimeric Interface

The environment of amino acid residues in protein tertiary structures and three types of interfaces formed by protein-protein association--in complexes, homodimers, and crystal lattices of monomeric proteins--has been analyzed in terms of the propensity values of the 20 amino acid residues to be in contact with a given residue. On the basis of the similarity of the environment, twenty residues can be divided into nine classes, which may correspond to a set of reduced amino acid alphabet. There is no appreciable change in the environment in going from the tertiary structure to the interface, those participating in the crystal contacts showing the maximum deviation. Contacts between identical residues are very prominent in homodimers and crystal dimers and arise due to 2-fold related association of residues lining the axis of rotation. These two types of interfaces, representing specific and nonspecific associations, are characterized by the types of residues that partake in "self-contacts"--most notably Leu in the former and Glu in the latter. The relative preference of residues to be involved in "self-contacts" can be used to develop a scoring function to identify homodimeric proteins from crystal structures. Thirty-four percent of such residues are fully conserved among homologous proteins in the homodimer dataset, as opposed to only 20% in crystal dimers. Results point to Leu being the stickiest of all amino acid residues, hence its widespread use in motifs, such as leucine zippers.