Templates are available to model nearly all complexes of structurally characterized proteins.

Traditional approaches to protein-protein docking sample the binding modes with no regard to similar experimentally determined structures (templates) of protein-protein complexes. Emerging template-based docking approaches utilize such similar complexes to determine the docking predictions. The docking problem assumes the knowledge of the participating proteins' structures. Thus, it provides the possibility of aligning the structures of the proteins and the template complexes. The progress in the development of template-based docking and the vast experience in template-based modeling of individual proteins show that, generally, such approaches are more reliable than the free modeling. The key aspect of this modeling paradigm is the availability of the templates. The current common perception is that due to the difficulties in experimental structure determination of protein-protein complexes, the pool of docking templates is insignificant, and thus a broad application of template-based docking is possible only at some future time. The results of our large scale, systematic study show that, surprisingly, in spite of the limited number of protein-protein complexes in the Protein Data Bank, docking templates can be found for complexes representing almost all the known protein-protein interactions, provided the components themselves have a known structure or can be homology-built. About one-third of the templates are of good quality when they are compared to experimental structures in test sets extracted from the Protein Data Bank and would be useful starting points in modeling the complexes. This finding dramatically expands our ability to model protein interactions, and has far-reaching implications for the protein docking field in general.

The targets of CAPRI rounds 20-27.

Eight CAPRI prediction rounds with a total of 15 targets were held in the years 2010-2012. Only five of the targets were protein assemblies comparable with those of earlier CAPRI rounds. In one target, the solvent positions at the interface had to be predicted; another was a protein-polysaccharide complex. The remainders were designed complexes issued from protein engineering experiments, and the prediction concerned either their structure or the binding affinity of the designed ligand. Affinity prediction was a new experiment in CAPRI, and a challenge for its participants. It pushed the community into developing novel procedures and score functions that will improve the performance of docking methods, help designing binders, and yield better structure-based estimates of the binding free energy of natural assemblies.

A survey of hemoglobin quaternary structures.

We perform an analysis of the quaternary structure and dimer/dimer interface in the crystal structures of 165 human hemoglobin tetramers; 112 are in the T, 17 the R, 14 the Y (or R2) state; 11 are high-affinity T state mutants, and 11 may either be intermediates between the states, or off the allosteric transition pathway. The tertiary structure is fixed within each state, in spite of the different ligands, mutations, and chemical modifications present in individual entries. The geometry of the tetramer assembly is essentially the same in all the R or the Y state entries; it is slightly different in high salt and low salt crystals of T state hemoglobins. The dimer/dimer interface differs in terms of size, chemical composition and polar interactions, between the states. It is loosely packed, like crystal packing contacts or the subunit interface of weakly associated homodimers, and unlike most oligomeric proteins, which have close-packed interfaces. The loose packing is most obvious in the liganded forms, where the tetramer is known to dissociate at low concentration. We identify cavities that contribute to the loose packing of the α1ß2 and α2ß1 contacts. Two pairs of cavities occur recurrently in both the T and the R state tetramers. They may contribute to the allosteric mechanism by facilitating the subunit movements and the tertiary structure changes that accompany the transition from T to R to Y.

The targets of CAPRI Rounds 13-19.

Seven rounds of CAPRI predictions with a total of 14 targets were held in the period June 2007-November 2009. In addition to protease/inhibitor complexes and complexes with G-proteins, some of the targets displayed novel features presenting new challenges to the predictors: a complex with RNA, a leucine zipper to be built ab initio, and a human-designed protein. Nine targets were unbound or required model building; the other five had a bound component. Thirteen were assessed, and the results show that the predictor and scorer groups submitted two- or three-star (medium or high quality) models for eight of them, one-star (acceptable) models for three but failed on the unbound RNA complex and another unbound target.

Side-chain rotamer transitions at protein-protein interfaces.

We compare the changes in side chain conformations that accompany the formation of protein-protein complexes, in residues forming either the interface or the remainder of the solvent-accessible surface of the proteins in the Docking Benchmark 3.0. We find that the interface residues undergo significantly more changes than other surface residues, and these changes are more likely to convert them from a high-energy torsion angle state to a lower-energy one than the reverse. Moreover, in both the unbound proteins and the complexes, the interface residues are more frequently found to be in a high-energy torsion angle state than the noninterface residues. As these differences exist before the binding step, they may be relevant to specificity and help in identifying binding sites for docking predictions.

Protein-protein docking benchmark version 4.0.

We updated our protein-protein docking benchmark to include complexes that became available since our previous release. As before, we only considered high-resolution complex structures that are nonredundant at the family-family pair level, for which the X-ray or NMR unbound structures of the constituent proteins are also available. Benchmark 4.0 adds 52 new complexes to the 124 cases of Benchmark 3.0, representing an increase of 42%. Thus, benchmark 4.0 provides 176 unbound-unbound cases that can be used for protein-protein docking method development and assessment. Seventeen of the newly added cases are enzyme-inhibitor complexes, and we found no new antigen-antibody complexes. Classifying the new cases according to expected difficulty for protein-protein docking algorithms gives 33 rigid body cases, 11 cases of medium difficulty, and 8 cases that are difficult. Benchmark 4.0 listings and processed structure files are publicly accessible at http://zlab.umassmed.edu/benchmark/.

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.

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.

The third CAPRI assessment meeting Toronto, Canada, April 20-21, 2007.

CAPRI is a community-wide experiment to test protein-protein docking methods in blind predictions. The Toronto meeting assessed structure predictions made from 2005-2007 on nine target protein-protein complexes or homodimers, and reported new developments in functions used to score predicted interactions, in treatment of conformational flexibility, and in taking nonstructural information into account in the predictions.

Residue conservation in viral capsid assembly.

To evaluate the evolutionary constraints placed on viral proteins by the structure and assembly of the capsid, we calculate Shannon entropies in the aligned sequences of 45 polypeptide chains in 32 icosahedral viruses, and relate these entropies to the residue location in the three-dimensional structure of the capsids. Three categories of residues have entropies lower than the chain average implying that they are better conserved than average: residues that are buried within a subunit (the protein core), residues that contain atoms buried at an interface between subunits (the interface core), and residues that contribute to several such interfaces. The interface core is also conserved in homomeric proteins and in transient protein-protein complexes, which have only one interface whereas capsids have many. In capsids, the subunit interfaces implicate most of the polypeptide chain: on average, 66% of the capsid residues are at an interface, 34% at more than one, and 47% at the interface core. Nevertheless, we observe that the degree of residue conservation can vary widely between interfaces within a capsid and between regions within an interface. The interfaces and regions of interfaces that show a low sequence variability are likely to play major roles in the self-assembly of the capsid, with implications on its mechanism that we discuss taking adeno-associated virus as an example.

Protein-protein docking benchmark version 3.0.

We present version 3.0 of our publicly available protein-protein docking benchmark. This update includes 40 new test cases, representing a 48% increase from Benchmark 2.0. For all of the new cases, the crystal structures of both binding partners are available. As with Benchmark 2.0, Structural Classification of Proteins (Murzin et al., J Mol Biol 1995;247:536-540) was used to remove redundant test cases. The 124 unbound-unbound test cases in Benchmark 3.0 are classified into 88 rigid-body cases, 19 medium-difficulty cases, and 17 difficult cases, based on the degree of conformational change at the interface upon complex formation. In addition to providing the community with more test cases for evaluating docking methods, the expansion of Benchmark 3.0 will facilitate the development of new algorithms that require a large number of training examples. Benchmark 3.0 is available to the public at http://zlab.bu.edu/benchmark.

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.

The targets of CAPRI rounds 6-12.

Six protein-protein complexes and two homodimeric proteins involved in a variety of biological processes were offered as targets to CAPRI by crystallographers in Rounds 6-12. CAPRI predictor groups had to predict their structure by docking the free proteins, which they did with a degree of success that depended largely on the amplitude of the conformation changes. In one case at least, the prediction pointed to alternative possibilities of interactions in the crystal of a complex, showing that docking methods have value even when there is an experimental structure.

Cloning, production, and purification of proteins for a medium-scale structural genomics project.

The South-Paris Yeast Structural Genomics Pilot Project (http://www.genomics.eu.org) aims at systematically expressing, purifying, and determining the three-dimensional structures of Saccharomyces cerevisiae proteins. We have already cloned 240 yeast open reading frames in the Escherichia coli pET system. Eighty-two percent of the targets can be expressed in E. coli, and 61% yield soluble protein. We have currently purified 58 proteins. Twelve X-ray structures have been solved, six are in progress, and six other proteins gave crystals. In this chapter, we present the general experimental flowchart applied for this project. One of the main difficulties encountered in this pilot project was the low solubility of a great number of target proteins. We have developed parallel strategies to recover these proteins from inclusion bodies, including refolding, coexpression with chaperones, and an in vitro expression system. A limited proteolysis protocol, developed to localize flexible regions in proteins that could hinder crystallization, is also described.

Genome-wide studies of protein-protein interaction.

Recent large-scale studies of protein complexes in yeast have demonstrated that the wide majority of proteins exist in the cell as parts of multicomponent assemblies, mostly novel and of unknown function. The structural and functional analysis of these complexes should be a priority for structural biologists in coming years. In silico methods such as docking simulations, which may contribute to this analysis, are being tested in the CAPRI community-wide experiment, which assesses blind predictions of the structure of protein-protein complexes.

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.

Revisiting the Voronoi description of protein-protein interfaces.

We developed a model of macromolecular interfaces based on the Voronoi diagram and the related alpha-complex, and we tested its properties on a set of 96 protein-protein complexes taken from the Protein Data Bank. The Voronoi model provides a natural definition of the interfaces, and it yields values of the number of interface atoms and of the interface area that have excellent correlation coefficients with those of the classical model based on solvent accessibility. Nevertheless, some atoms that do not lose solvent accessibility are part of the interface defined by the Voronoi model. The Voronoi model provides robust definitions of the curvature and of the connectivity of the interfaces, and leads to estimates of these features that generally agree with other approaches. Our implementation of the model allows an analysis of protein-water contacts that highlights the role of structural water molecules at protein-protein interfaces.

Crystal structure of the yeast His6 enzyme suggests a reaction mechanism.

The Saccharomyces cerevisiae His6 gene codes for the enzyme phosphoribosyl-5-amino-1-phosphoribosyl-4-imidazolecarboxamide isomerase, catalyzing the fourth step in histidine biosynthesis. To get an insight into the structure and function of this enzyme, we determined its X-ray structure at a resolution of 1.30 A using the anomalous diffraction signal of the protein's sulphur atoms at 1.77 A wavelength. His6 folds in an (alpha/beta)8 barrel similar to HisA, which performs the same function in bacteria and archaea. We found a citrate molecule from the buffer bound in a pocket near the expected position of the active site and used it to model the open form of the substrate (phosphoribulosyl moiety), which is a reaction intermediate. This model enables us to identify catalytic residues and to propose a reaction mechanism where two aspartates act as acid/base catalysts: Asp134 as a proton donor for ring opening, and Asp9 as a proton acceptor and donor during enolization of the aminoaldose. Asp9 is conserved in yeast His6 and bacterial or archaeal HisA sequences, and Asp134 has equivalents in both HisA and TrpF, but they occur at a different position in the protein sequence.

Crystal structure of the YDR533c S. cerevisiae protein, a class II member of the Hsp31 family.

The ORF YDR533c from Saccharomyces cerevisiae codes for a 25.5 kDa protein of unknown biochemical function. Transcriptome analysis of yeast has shown that this gene is activated in response to various stress conditions together with proteins belonging to the heat shock family. In order to clarify its biochemical function, we determined the crystal structure of YDR533c to 1.85 A resolution by the single anomalous diffraction method. The protein possesses an alpha/beta hydrolase fold and a putative Cys-His-Glu catalytic triad common to a large enzyme family containing proteases, amidotransferases, lipases, and esterases. The protein has strong structural resemblance with the E. coli Hsp31 protein and the intracellular protease I from Pyrococcus horikoshii, which are considered class I and class III members of the Hsp31 family, respectively. Detailed structural analysis strongly suggests that the YDR533c protein crystal structure is the first one of a class II member of the Hsp31 family.