Structures Of Protein-protein Complexes Research Articles

Protein–protein docking: Progress in CAPRI rounds 6–12 using a combination of methods: The introduction of steered solvated molecular dynamics

In recent rounds of CAPRI, the Bii group has employed a combination of techniques for the prediction of the structure of protein-protein complexes. We currently use third-party software for rigid-body and semiflexible docking (MolFit, 3D-Dock, RosettaDock), and our own steered molecular dynamics (SMD) technique for flexible refinement. SMD has also been found to be useful for discriminating near-native from false positive docking decoys. In addition to this, a variety of sources of information, including multiple descriptors of interface quality combined with a QSAR-like technique, published biological information, and continuum electrostatics calculations, are also used in the assessment of candidate complexes. We shall concentrate on results for CAPRI rounds 9-11 (targets 24-27). In these rounds, the Bii group has been successful in submitting a medium quality model for each of CAPRI targets 25 and 26, and a model of acceptable quality for target 27.

Protein–protein docking in CAPRI using ATTRACT to account for global and local flexibility

A reduced protein model combined with a systematic docking approach has been employed to predict protein-protein complex structures in CAPRI rounds 6-11. The docking approach termed ATTRACT is based on energy minimization in translational and rotational degrees of freedom of one protein with respect to the second protein starting from many thousand initial protein partner placements. It also allows for approximate inclusion of global flexibility of protein partners during systematic docking by conformational relaxation of the partner proteins in precalculated soft collective backbone degrees of freedom. We have submitted models for six targets, achieved acceptable docking solutions for two targets, and predicted >20% correct contacts for five targets. Possible improvements of the docking approach in particular at the scoring and refinement steps are discussed.

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.

FastContact: a free energy scoring tool for protein-protein complex structures

‘FastContact’ is a server that estimates the direct electrostatic and desolvation interaction free energy between two proteins in units of kcal/mol. Users submit two proteins in PDB format, and the output is emailed back to the user in three files: one output file, and the two processed proteins. Besides the electrostatic and desolvation free energy, the server reports residue contact free energies that rapidly highlight the hotspots of the interaction and evaluates the van der Waals interaction using CHARMm. Response time is ∼1 min. The server has been successfully tested and validated, scoring refined complex structures and blind sets of docking decoys, as well as proven useful predicting protein interactions. ‘FastContact’ offers unique capabilities from biophysical insights to scoring and identifying important contacts.

Trp/Met/Phe Hot Spots in Protein-Protein Interactions: Potential Targets in Drug Design

Protein-protein interactions are crucial to biological functions. Consequently, designing drugs to control protein-protein interactions is receiving increasing attention. Protein structures can associate in different ways. Analysis of the structures of protein-protein complexes using amino acid sequence order-independent multiple structural comparison algorithms, led us to conclude that the amino acids Trp, Met, and Phe are important for protein-protein interactions. Hence, in principle, drug design targeting the Trp/Met/Phe should modulate protein functions effectively. Several clusters of the Trp/Met/Phe residues are involved in the p53 protein-protein interactions. The best example in this regard is the Phe19/Trp23 of p53, which binds to transcriptional factors and to the MDM2 protein. In the HIV related proteins, the Trp/Met/Phe residues have roles in the dimerization of the transcriptase (p51/p66) and in cell-fusion processes, including the gp120-CD4 interaction and the gp41 six-helix bundle formation. Trp/Met/Phe residues are preferred in 'normal' functional protein-protein interactions and they also appear to be exploited in amyloid formation, especially the phenylalanine. Comparison of binding propensity and amyloid formation preference reveals that apart from Lysine, Isoleucine is the least structurally conserved in protein binding sites and has a high propensity in sequences forming amyloids. Thus, this may suggest that nature tends to avoid Ile conservation in protein-protein interaction to avoid amyloid formation. In this regards, Trp/Met/Phe as well as Ile may be targeted to modulate protein-protein interaction.

PROTCOM: searchable database of protein complexes enhanced with domain-domain structures

The database of protein complexes (PROTCOM) is a compilation of known 3D structures of protein–protein complexes enriched with artificially created domain–domain structures using the available entries in the Protein Data Bank. The domain–domain structures are generated by parsing single chain structures into loosely connected domains and are important features of the database. The database () could be used for benchmarking purposes of the docking and other algorithms for predicting 3D structures of protein–protein complexes. The database can be utilized as a template database in the homology or threading methods for modeling the 3D structures of unknown protein–protein complexes. PROTCOM provides the scientific community with an integrated set of tools for browsing, searching, visualizing and downloading a pool of protein complexes. The user is given the option to select a subset of entries using a combination of up to 10 different criteria. As on July 2006 the database contains 1770 entries, each of which consists of the known 3D structures and additional relevant information that can be displayed either in text-only or in visual mode.

PIPER: An FFT‐based protein docking program with pairwise potentials

The Fast Fourier Transform (FFT) correlation approach to protein-protein docking can evaluate the energies of billions of docked conformations on a grid if the energy is described in the form of a correlation function. Here, this restriction is removed, and the approach is efficiently used with pairwise interaction potentials that substantially improve the docking results. The basic idea is approximating the interaction matrix by its eigenvectors corresponding to the few dominant eigenvalues, resulting in an energy expression written as the sum of a few correlation functions, and solving the problem by repeated FFT calculations. In addition to describing how the method is implemented, we present a novel class of structure-based pairwise intermolecular potentials. The DARS (Decoys As the Reference State) potentials are extracted from structures of protein-protein complexes and use large sets of docked conformations as decoys to derive atom pair distributions in the reference state. The current version of the DARS potential works well for enzyme-inhibitor complexes. With the new FFT-based program, DARS provides much better docking results than the earlier approaches, in many cases generating 50% more near-native docked conformations. Although the potential is far from optimal for antibody-antigen pairs, the results are still slightly better than those given by an earlier FFT method. The docking program PIPER is freely available for noncommercial applications.

RESIDUE CONSERVATION INFORMATION FOR GENERATING NEAR-NATIVE STRUCTURES IN PROTEIN–PROTEIN DOCKING

Protein-protein docking algorithms typically generate large numbers of possible complex structures with only a few of them resembling the native structure. Recently (Duan et al., Protein Sci, 14:316-218, 2005), it was observed that the surface density of conserved residue positions is high at the interface regions of interacting protein surfaces, except for antibody-antigen complexes, where a lesser number of conserved positions than average is observed at the interface regions. Using this observation, we identified putative interacting regions on the surface of interacting partners and significantly improved docking results by assigning top ranks to near-native complex structures. In this paper, we combine the residue conservation information with a widely used shape complementarity algorithm to generate candidate complex structures with a higher percentage of near-native structures (hits). What is new in this work is that the conservation information is used early in the generation stage and not only in the ranking stage of the docking algorithm. This results in a significantly larger number of generated hits and an improved predictive ability in identifying the native structure of protein-protein complexes. We report on results from 48 well-characterized protein complexes, which have enough residue conservation information from the same 59 benchmark complexes used in our previous work. We compute conservation indices of residue positions on the surfaces of interacting proteins using available homologous sequences from UNIPROT and calculate the solvent accessible surface area. We combine this information with shape-complementarity scores to generate candidate protein-protein complex structures. When compared with pure shape-complementarity algorithms, performed by FTDock, our method results in significantly more hits, with the improvement being over 100% in many instances. We demonstrate that residue conservation information is useful not only in refinement and scoring of docking solutions, but also helpful in enrichment of near-native-structures during the generation of candidate geometries of complex structures.

Lanthanide Labeling Offers Fast NMR Approach to 3D Structure Determinations of Protein−Protein Complexes

A novel nuclear magnetic resonance (NMR) strategy based on labeling with lanthanides achieves rapid determinations of accurate three-dimensional (3D) structures of protein-protein complexes. The method employs pseudocontact shifts (PCS) induced by a site-specifically bound lanthanide ion to anchor the coordinate system of the magnetic susceptibility tensor in the molecular frames of the two molecules. Simple superposition of the tensors detected in the two protein molecules brings them together in a 3D model of the protein-protein complex. The method is demonstrated with the 30 kDa complex between two subunits of Escherichia coli polymerase III, comprising the N-terminal domain of the exonuclease subunit epsilon and the subunit theta. The 3D structures of the individual molecules were docked based on a limited number of PCS observed in 2D 15N-heteronuclear single quantum coherence spectra. Degeneracies in the mutual orientation of the protein structures were resolved by the use of two different lanthanide ions, Dy3+ and Er3+.

Prediction of interface residues in protein–protein complexes by a consensus neural network method: Test against NMR data

The number of structures of protein-protein complexes deposited to the Protein Data Bank is growing rapidly. These structures embed important information for predicting structures of new protein complexes. This motivated us to develop the PPISP method for predicting interface residues in protein-protein complexes. In PPISP, sequence profiles and solvent accessibility of spatially neighboring surface residues were used as input to a neural network. The network was trained on native interface residues collected from the Protein Data Bank. The prediction accuracy at the time was 70% with 47% coverage of native interface residues. Now we have extensively improved PPISP. The training set now consisted of 1156 nonhomologous protein chains. Test on a set of 100 nonhomologous protein chains showed that the prediction accuracy is now increased to 80% with 51% coverage. To solve the problem of over-prediction and under-prediction associated with individual neural network models, we developed a consensus method that combines predictions from multiple models with different levels of accuracy and coverage. Applied on a benchmark set of 68 proteins for protein-protein docking, the consensus approach outperformed the best individual models by 3-8 percentage points in accuracy. To demonstrate the predictive power of cons-PPISP, eight complex-forming proteins with interfaces characterized by NMR were tested. These proteins are nonhomologous to the training set and have a total of 144 interface residues identified by chemical shift perturbation. cons-PPISP predicted 174 interface residues with 69% accuracy and 47% coverage and promises to complement experimental techniques in characterizing protein-protein interfaces. .

SiteEngines: recognition and comparison of binding sites and protein-protein interfaces

Protein surface regions with similar physicochemical properties and shapes may perform similar functions and bind similar binding partners. Here we present two web servers and software packages for recognition of the similarity of binding sites and interfaces. Both methods recognize local geometrical and physicochemical similarity, which can be present even in the absence of overall sequence or fold similarity. The first method, SiteEngine (), receives as an input two protein structures and searches the complete surface of one protein for regions similar to the binding site of the other. The second, Interface-to-Interface (I2I)-SiteEngine (), compares protein–protein interfaces, which are regions of interaction between two protein molecules. It receives as an input two structures of protein–protein complexes, extracts the interfaces and finds the three-dimensional transformation that maximizes the similarity between two pairs of interacting binding sites. The output of both servers consists of a superimposition in PDB file format and a list of physicochemical properties shared by the compared entities. The methods are highly efficient and the freely available software packages are suitable for large-scale database searches of the entire PDB.

ATTRACT: Protein–protein docking in CAPRI using a reduced protein model

Protein-protein complex structures have been predicted for CAPRI Rounds 3 and 5 using a reduced protein model. Proteins are represented by up to 3 pseudoatoms per amino acid. The docking approach termed ATTRACT is based on energy minimization in translational and rotational degrees of freedom of one protein with respect to another protein. The reduced protein model allows one to perform systematic docking minimization of many thousand start structures in reasonable computer time. Flexibility of critical surface side-chains can be accounted for by a multiple conformational copy approach. The multicopy approach allows simultaneous adjustment of side-chain conformations and optimization of translational and rotational degrees of freedom of one protein with respect to the partner during docking. For 3 (Targets 8, 14, and 19) out of 5 CAPRI targets, the approach resulted in predictions in close agreement with experiment [root-mean-square deviation (RMSD) of backbone atoms within 10 A of the protein-protein interface < 1.8 A]. The comparison of predicted and experimental structures of the CAPRI targets indicates that besides local conformational changes (e.g., changes in side-chain conformations), global conformational changes of the protein backbone can be critical for complex formation. These conformational changes not accounted for during docking are a likely reason for the unrealistic predictions in 2 cases (Targets 9 and 18).

Assessment of CAPRI predictions in rounds 3–5 shows progress in docking procedures

The current status of docking procedures for predicting protein-protein interactions starting from their three-dimensional (3D) structure is reassessed by evaluating blind predictions, performed during 2003-2004 as part of Rounds 3-5 of the community-wide experiment on Critical Assessment of PRedicted Interactions (CAPRI). Ten newly determined structures of protein-protein complexes were used as targets for these rounds. They comprised 2 enzyme-inhibitor complexes, 2 antigen-antibody complexes, 2 complexes involved in cellular signaling, 2 homo-oligomers, and a complex between 2 components of the bacterial cellulosome. For most targets, the predictors were given the experimental structures of 1 unbound and 1 bound component, with the latter in a random orientation. For some, the structure of the free component was derived from that of a related protein, requiring the use of homology modeling. In some of the targets, significant differences in conformation were displayed between the bound and unbound components, representing a major challenge for the docking procedures. For 1 target, predictions could not go to completion. In total, 1866 predictions submitted by 30 groups were evaluated. Over one-third of these groups applied completely novel docking algorithms and scoring functions, with several of them specifically addressing the challenge of dealing with side-chain and backbone flexibility. The quality of the predicted interactions was evaluated by comparison to the experimental structures of the targets, made available for the evaluation, using the well-agreed-upon criteria used previously. Twenty-four groups, which for the first time included an automatic Web server, produced predictions ranking from acceptable to highly accurate for all targets, including those where the structures of the bound and unbound forms differed substantially. These results and a brief survey of the methods used by participants of CAPRI Rounds 3-5 suggest that genuine progress in the performance of docking methods is being achieved, with CAPRI acting as the catalyst.

Assessing predictions of protein–protein interaction: The CAPRI experiment

The Critical Assessment of PRedicted Interactions (CAPRI) experiment was designed in 2000 to test protein docking algorithms in blind predictions of the structure of protein-protein complexes. In four years, 17 complexes offered by crystallographers as targets prior to publication, have been subjected to structure prediction by docking their two components. Models of these complexes were submitted by predictor groups and assessed by comparing their geometry to the X-ray structure and by evaluating the quality of the prediction of the regions of interaction and of the pair wise residue contacts. Prediction was successful on 12 of the 17 targets, most of the failures being due to large conformation changes that the algorithms could not cope with. Progress in the prediction quality observed in four years indicates that the experiment is a powerful incentive to develop new procedures that allow for flexibility during docking and incorporate nonstructural information. We therefore call upon structural biologists who study protein-protein complexes to provide targets for further rounds of CAPRI predictions.

Topology of small-world networks of protein-protein complex structures

The majority of real examples of small-world networks exhibit a power law distribution of edges among the nodes, therefore not fitting into the wiring model proposed by Watts and Strogatz. However, protein structures can be modeled as small-world networks, with a distribution of the number of links decaying exponentially as in the case of this wiring model. We approach the protein-protein interaction mechanism by viewing it as a particular rewiring occurring in the system of two small-world networks represented by the monomers, where a re-arrangement of links takes place upon dimerization leaving the small-world character in the dimer network. Due to this rewiring, the most central residues at the complex interfaces tend to form clusters, which are not homogenously distributed. We show that these highly central residues are strongly correlated with the presence of hot spots of binding free energy. ao-mesa@fujirebio.co.jp http://www.fujirebio.co.jp/support/index.php (under construction).

A Brownian Dynamics Study: The Effect of a Membrane Environment on an Electron Transfer System

During the past few years, three-dimensional crystal structures of many of the important integral membrane proteins responsible for the bioenergetic processes of photosynthesis and respiration have been determined. Moreover, a few crystal structures of protein-protein complexes have become available that characterize the interaction between those membrane proteins and the electron carrier protein cytochrome c. Here, we address the association kinetics for binding of cytochrome c to cytochrome c oxidase (COX) from Paracoccus denitrificans by Brownian dynamics simulations. The effects of ionic strength and protein mutations were studied for two different cytochrome c species: the positively charged, dipolar horse heart cytochrome c and the negatively charged physiological electron transfer partner cytochrome c 552. We studied association toward “naked” COX and toward membrane-embedded COX where the membrane is represented as an uncharged DPPC bilayer modeled in atomistic detail. For the nonnatural association toward “naked” COX, the association rates are >100 times larger for horse heart cytochrome c than for cytochrome c 552. Interestingly, the presence of the lipid bilayer leads to a dramatic decrease of the association rate of horse heart cytochrome c, but slightly enhances association of cytochrome c 552, leading to very similar association rates of both proteins to membrane-embedded COX. This finding from computational modeling studies may reflect the optimization of surface patches and of the total net charge on electron transfer pairs in nature.

Identification of the protein kinase A regulatory RIalpha-catalytic subunit interface by amide H/2H exchange and protein docking.

An important goal after structural genomics is to build up the structures of higher-order protein-protein complexes from structures of the individual subunits. Often structures of higher order complexes are difficult to obtain by crystallography. We have used an alternative approach in which the structures of the individual catalytic (C) subunit and RIalpha regulatory (R) subunit of PKA were first subjected to computational docking, and the top 100,000 solutions were subsequently filtered based on amide hydrogen/deuterium (H/2H) exchange interface protection data. The resulting set of filtered solutions forms an ensemble of structures in which, besides the inhibitor peptide binding site, a flat interface between the C-terminal lobe of the C-subunit and the A- and B-helices of RIalpha is uniquely identified. This holoenzyme structure satisfies all previous experimental data on the complex and allows prediction of new contacts between the two subunits.

Determination of molecular alignment tensors without backbone resonance assignment: Aid to rapid analysis of protein-protein interactions.

Based on high-resolution structures of the free molecules accurate determination of structures of protein complexes by NMR spectroscopy is possible using residual dipolar couplings. In order, however, to be able to apply these methods, protein backbone resonances have to be assigned first. This NMR assignment process is particularly difficult and time consuming for protein sizes above 20 kDa. Here we show that, when NMR resonances belonging to a specific amino acid type are selected either by amino acid specific labeling, by their characteristic Calpha/Cbeta chemical shifts or by dedicated NMR experiments, molecular alignment tensors of proteins up to 80 kDa can be determined without prior backbone resonance assignment. This offers the opportunity to greatly accelerate determination of three-dimensional structures of protein-protein and protein-ligand complexes, and validation of multimeric states of proteins. Moreover, exhaustive back-calculation can be performed using only 1DNH couplings. Therefore, it avoids expensive 13C-labeling and it gives access to orientational information for large proteins that strongly aggregate at concentrations above 50 microM, i.e., experimental conditions where 3D triple resonance experiments are not sensitive enough to allow backbone resonance assignment.

Structures of protein-protein complexes are docked using only NMR restraints from residual dipolar coupling and chemical shift perturbations.

NMR structures of protein-protein and protein-ligand complexes rely heavily on intermolecular NOEs. Recent work has shown that if no significant conformational changes occur upon complex formation residual dipolar coupling can replace most of the NOE restraints in protein-protein complexes, while restraints derived from chemical shift perturbations can largely replace intermolecular NOEs in protein-ligand structures. By combining restraints from chemical shift perturbations with orientation restraints derived from measurements of residual dipolar couplings, we show that the structure of the EIN-HPr complex can be calculated without NOE restraints. The final structure, built from the crystal structures of EIN and HPr in their uncomplexed form and docked only with NMR restraints, places HPr within 2.5 A of the position determined from the mean NMR structure of the complex.

Synchrotron Protein Footprinting: A Technique to Investigate Protein-Protein Interactions

Traditional approaches for macromolecular structure elucidation, including NMR, crystallography and cryo-EM have made significant progress in defining the structures of protein-protein complexes. A substantial number of macromolecular structures, however, have not been examined with atomic detail due to sample size and heterogeneity, or resolution limitations of the technique; therefore, the general applicability of each method is greatly reduced. Synchrotron footprinting attempts to bridge the gap in these methods by monitoring changes in accessible surface areas of discrete macromolecular moieties. As evidenced by our previous studies on RNA folding and DNA-protein interactions, the three-dimensional structure is probed by examining the reactions of these moieties with hydroxyl radicals generated by synchrotron X-rays. Here we report the application of synchrotron foot- printing to the investigation of protein-protein interactions, as the novel technique has been utilized to successfully map the contact sites of gelsolin segment-1 in the gelsolin segment 1/actin complex. Footprinting results demonstrate that phenylalanine 104, located on the actin binding helix of gelsolin segment 1, is protected from hydroxyl radical modification in the presence of actin. This change in reactivity results from the specific protection of gel- solin segment-1, consistent with the substantial decrease in solvent accessibility of F104 upon actin binding, as calculated from the crystal structural of the gelsolin segment 1/actin complex. The results presented here establish synchrotron footprinting as a broadly applicable method to probe structural features of macromolecular complexes that are not amenable to conventional approaches.