A "solvated rotamer" approach to modeling water-mediated hydrogen bonds at protein-protein interfaces

Water-mediated hydrogen bonds play critical roles at protein-protein and protein-nucleic acid interfaces, and the interactions formed by discrete water molecules cannot be captured using continuum solvent models. We describe a simple model for the energetics of water-mediated hydrogen bonds, and show that, together with knowledge of the positions of buried water molecules observed in X-ray crystal structures, the model improves the prediction of free-energy changes upon mutation at protein-protein interfaces, and the recovery of native amino acid sequences in protein interface design calculations. We then describe a "solvated rotamer" approach to efficiently predict the positions of water molecules, at protein-protein interfaces and in monomeric proteins, that is compatible with widely used rotamer-based side-chain packing and protein design algorithms. Finally, we examine the extent to which the predicted water molecules can be used to improve prediction of amino acid identities and protein-protein interface stability, and discuss avenues for overcoming current limitations of the approach.     

Kinase inhibitors and the case for CH...O hydrogen bonds in protein-ligand binding

Although the hydrogen bond is known to be an important mediator of intermolecular interactions, there has yet to be an analysis of the role of CH...O hydrogen bonds in protein-ligand complexes. In this work, we present evidence for such nonstandard hydrogen bonds from a survey of aromatic ligands in 184 kinase crystal structures and 358 high-resolution structures from the Protein Data Bank. CH groups adjacent to the positively charged nitrogen of nicotinamide exhibit geometric preferences strongly suggestive of hydrogen bonding interactions, as do heterocyclic CH groups in kinase ligands, while other aromatic CH groups do not exhibit these characteristics. Ab initio calculations reveal a considerable range of CH...O hydrogen bonding potentials among different aromatic ring systems, with nicotinamide and heterocycles preferred in kinase inhibitors showing particularly favorable interactions. These results provide compelling evidence for the existence of CH...O hydrogen bonds in protein-ligand interactions, as well as information on the relative strength of various aromatic CH donors. Such knowledge will be of considerable value in protein modeling, ligand design, and structure-activity analysis.     

Hydrogen bonding isosteres: bimolecular carboxylic acid and amine-N-oxide interactions mediated via CH...O hydrogen bonds

The X-ray structures of cocrystals between 2,2'-dipyridyl-N,N'-dioxide (1) with fumaric acid (2), itaconic acid (3), succinic acid (4), and oxalic acid (5) were solved to determine if concurrent CH...O interactions were capable of orienting the bimolecular association of the two molecules. Cocrystals 1.2, 1.3 and 1.4 produce cyclic hydrogen bonded motifs employing pair-wise OH...O and CH...O hydrogen bonds, whereas cocrystal 1.5 forms analogous OH...O hydrogen bonds with a different set of intermolecular CH...O hydrogen bonds. Evidence of cocrystal formation was also observed for these complexes by differential scanning calorimetry and FT-IR spectroscopy. The structures of 1.2, 1.3 and 1.4 demonstrate the potential of the pair-wise OH...O and CH...O hydrogen bonding interactions and serve to illustrate their use as hydrogen bonding isosteres in crystal engineering, molecular recognition, and drug design.     

Investigations on unconventional hydrogen bonds in RNA binding proteins: the role of CH...O=C interactions

We have investigated the roles played by C-H...O=C interactions in RNA binding proteins. There was an average of 78 CH...O=C interactions per protein and also there was an average of one significant CH...O=C interactions for every 6 residues in the 59 RNA binding proteins studied. Main chain-Main chain (MM) CH...O=C interactions are the predominant type of interactions in RNA binding proteins. The donor atom contribution to CH...O=C interactions was mainly from aliphatic residues. The acceptor atom contribution for MM CH...O=C interactions was mainly from Val, Phe, Leu, Ile, Arg and Ala. The secondary structure preference analysis of CH...O=C interacting residues showed that, Arg, Gln, Glu and Tyr preferred to be in helix, while Ala, Asp, Cys, Gly, Ile, Leu, Lys, Met, Phe, Trp and Val preferred to be in strand conformation. Most of the CH...O=C interacting polar amino acid residues were solvent exposed while, majority of the CH...O=C interacting non polar residues were excluded from the solvent. Long and medium-range CH...O=C interactions are the predominant type of interactions in RNA binding proteins. More than 50% of CH...O=C interacting residues had a higher conservation score. Significant percentage of CH...O=C interacting residues had one or more stabilization centers. Sixty-six percent of the theoretically predicted stabilizing residues were also involved in CH...O=C interactions and hence these residues may also contribute additional stability to RNA binding proteins.     

N-H...O, O-H...O, and C-H...O hydrogen bonds in protein-ligand complexes: strong and weak interactions in molecular recognition

The characteristics of N-H...O, O-H...O, and C-H...O hydrogen bonds are examined in a group of 28 high-resolution crystal structures of protein-ligand complexes from the Protein Data Bank and compared with interactions found in small-molecule crystal structures from the Cambridge Structural Database. It is found that both strong and weak hydrogen bonds are involved in ligand binding. Because of the prevalence of multifurcation, the restrictive geometrical criteria set up for hydrogen bonds in small-molecule crystal structures may need to be relaxed in macromolecular structures. For example, there are definite deviations from linearity for the hydrogen bonds in protein-ligand complexes. The formation of C-H...O hydrogen bonds is influenced by the activation of the C(alpha)-H atoms and by the flexibility of the side-chain atoms. In contrast to small-molecule structures, anticooperative geometries are common in the macromolecular structures studied here, and there is a gradual lengthening as the extent of furcation increases. C-H...O bonds formed by Gly, Phe, and Tyr residues are noteworthy. The numbers of hydrogen bond donors and acceptors agree with Lipinski's "rule of five" that predicts drug-like properties. Hydrogen bonds formed by water are also seen to be relevant in ligand binding. Ligand C-H...O(w) interactions are abundant when compared to N-H...O(w) and O-H...O(w). This suggests that ligands prefer to use their stronger hydrogen bond capabilities for use with the protein residues, leaving the weaker interactions to bind with water. In summary, the interplay between strong and weak interactions in ligand binding possibly leads to a satisfactory enthalpy-entropy balance. The implications of these results to crystallographic refinement and molecular dynamics software are discussed.     

Non-classical hydrogen bonds in interleukins: the role of C-H...O interactions

The role of classical hydrogen bonds in the structural stability of biological macro-molecules is well understood. In the present study, we explore the influence of C-H...O interactions in relation to other environmental preferences in interleukins. Main chain-main chain interactions are predominant. Pro residues might stabilize helices and strands by C-H...O H-bonds in interleukins. Majority of the C-H...O interacting residues were solvent exposed. 62% of C-H...O interactions was long-range interactions. The results presented in this study might be useful for structural stability studies in interleukins.     

Side-chain conformational entropy at protein-protein interfaces

Protein-protein interactions are the key to many biological processes. How proteins selectively and correctly associate with their required protein partner(s) is still unclear. Previous studies of this "protein-docking problem" have found that shape complementarity is a major determinant of interaction, but the detailed balance of energy contributions to association remains unclear. This study estimates side-chain conformational entropy (per unit solvent accessible area) for various protein surface regions, using a self-consistent mean field calculation of rotamer probabilities. Interfacial surface regions were less flexible than the rest of the protein surface for calculations with monomers extracted from homodimer datasets in 21 of 25 cases, and in 8 of 9 for the large protomer from heterodimer datasets. In surface patch analysis, based on side-chain conformational entropy, 68% of true interfaces were ranked top for the homodimer set and 66% for the large protomer/heterodimer set. The results indicate that addition of a side-chain entropic term could significantly improve empirical calculations of protein-protein association.     

Alignment of non-covalent interactions at protein-protein interfaces

Background

The study and comparison of protein-protein interfaces is essential for the understanding of the mechanisms of interaction between proteins. While there are many methods for comparing protein structures and protein binding sites, so far no methods have been reported for comparing the geometry of non-covalent interactions occurring at protein-protein interfaces.

Methodology/Principal Findings

Here we present a method for aligning non-covalent interactions between different protein-protein interfaces. The method aligns the vector representations of van der Waals interactions and hydrogen bonds based on their geometry. The method has been applied to a dataset which comprises a variety of protein-protein interfaces. The alignments are consistent to a large extent with the results obtained using two other complementary approaches. In addition, we apply the method to three examples of protein mimicry. The method successfully aligns respective interfaces and allows for recognizing conserved interface regions.

Conclusions/Significance

The Galinter method has been validated in the comparison of interfaces in which homologous subunits are involved, including cases of mimicry. The method is also applicable to comparing interfaces involving non-peptidic compounds. Galinter assists users in identifying local interface regions with similar patterns of non-covalent interactions. This is particularly relevant to the investigation of the molecular basis of interaction mimicry.     

C-H...O hydrogen bonds in the nuclear receptor RARgamma--a potential tool for drug selectivity

Hydrogen bonds between polarized atoms play a crucial role in protein interactions and are often used in drug design, which usually neglects the potential of C-H...O hydrogen bonds. The 1.4 A resolution crystal structure of the ligand binding domain of the retinoic acid receptor RARgamma complexed with the retinoid SR11254 reveals several types of C-H...O hydrogen bonds. A striking example is the hydroxyl group of the ligand that acts as an H bond donor and acceptor, leading to a synergy between classical and C-H...O hydrogen bonds. This interaction introduces both specificity and affinity within the hydrophobic ligand pocket. The similarity of intraprotein and protein-ligand C-H...O interactions suggests that such bonds should be considered in rational drug design approaches.     

The occurrence of C--H...O hydrogen bonds in alpha-helices and helix termini in globular proteins

A comprehensive database analysis of C--H...O hydrogen bonds in 3124 alpha-helices and their corresponding helix termini has been carried out from a nonredundant data set of high-resolution globular protein structures resolved at better than 2.0 A in order to investigate their role in the helix, the important protein secondary structural element. The possible occurrence of 5 --> 1 C--H...O hydrogen bond between the ith residue CH group and (i - 4)th residue C==O with C...O < or = 3.8 A is studied, considering as potential donors the main-chain Calpha and the side-chain carbon atoms Cbeta, Cgamma, Cdelta and Cepsilon. Similar analysis has been carried out for 4 --> 1 C--H...O hydrogen bonds, since the C--H...O hydrogen bonds found in helices are predominantly of type 5 --> 1 or 4 --> 1. A total of 17,367 (9310 of type 5 --> 1 and 8057 of type 4 --> 1) C--H...O hydrogen bonds are found to satisfy the selected criteria. The average stereochemical parameters for the data set suggest that the observed C--H...O hydrogen bonds are attractive interactions. Our analysis reveals that the Cgamma and Cbeta hydrogen atom(s) are frequently involved in such hydrogen bonds. A marked preference is noticed for aliphatic beta-branched residue Ile to participate in 5 --> 1 C--H...O hydrogen bonds involving methylene Cgamma 1 atom as donor in alpha-helices. This may be an enthalpic compensation for the greater loss of side-chain conformational entropy for beta-branched amino acids due to the constraint on side-chain torsion angle, namely, chi1, when they occur in helices. The preference of amino acids for 4 --> 1 C--H...O hydrogen bonds is found to be more for Asp, Cys, and for aromatic residues Trp, Phe, and His. Interestingly, overall propensity for C--H...O hydrogen bonds shows that a majority of the helix favoring residues such as Met, Glu, Arg, Lys, Leu, and Gln, which also have large side-chains, prefer to be involved in such types of weak attractive interactions in helices. The amino acid side-chains that participate in C--H...O interactions are found to shield the acceptor carbonyl oxygen atom from the solvent. In addition, C--H...O hydrogen bonds are present along with helix stabilizing salt bridges. A novel helix terminating interaction motif, X-Gly with Gly at C(cap) position having 5 --> 1 Calpha--H...O, and a chain reversal structural motif having 1 --> 5 Calpha-H...O have been identified and discussed. Our analysis highlights that a multitude of local C--H...O hydrogen bonds formed by a variety of amino acid side-chains and Calpha hydrogen atoms occur in helices and more so at the helix termini. It may be surmised that the main-chain Calpha and the side-chain CH that participate in C--H...O hydrogen bonds collectively augment the cohesive energy and thereby contribute together with the classical N--H...O hydrogen bonds and other interactions to the overall stability of helix and therefore of proteins.     

Residue frequencies and pairing preferences at protein-protein interfaces

We used a nonredundant set of 621 protein-protein interfaces of known high-resolution structure to derive residue composition and residue-residue contact preferences. The residue composition at the interfaces, in entire proteins and in whole genomes correlates well, indicating the statistical strength of the data set. Differences between amino acid distributions were observed for interfaces with buried surface area of less than 1,000 A(2) versus interfaces with area of more than 5,000 A(2). Hydrophobic residues were abundant in large interfaces while polar residues were more abundant in small interfaces. The largest residue-residue preferences at the interface were recorded for interactions between pairs of large hydrophobic residues, such as Trp and Leu, and the smallest preferences for pairs of small residues, such as Gly and Ala. On average, contacts between pairs of hydrophobic and polar residues were unfavorable, and the charged residues tended to pair subject to charge complementarity, in agreement with previous reports. A bootstrap procedure, lacking from previous studies, was used for error estimation. It showed that the statistical errors in the set of pairing preferences are generally small; the average standard error is approximately 0.2, i.e., about 8% of the average value of the pairwise index (2.9). However, for a few pairs (e.g., Ser-Ser and Glu-Asp) the standard error is larger in magnitude than the pairing index, which makes it impossible to tell whether contact formation is favorable or unfavorable. The results are interpreted using physicochemical factors and their implications for the energetics of complex formation and for protein docking are discussed. Proteins 2001;43:89-102.     

Hydration of protein-protein interfaces

We present an analysis of the water molecules immobilized at the protein-protein interfaces of 115 homodimeric proteins and 46 protein-protein complexes, and compare them with 173 large crystal packing interfaces representing nonspecific interactions. With an average of 15 waters per 1000 A2 of interface area, the crystal packing interfaces are more hydrated than the specific interfaces of homodimers and complexes, which have 10-11 waters per 1000 A2, reflecting the more hydrophilic composition of crystal packing interfaces. Very different patterns of hydration are observed: Water molecules may form a ring around interfaces that remain "dry," or they may permeate "wet" interfaces. A majority of the specific interfaces are dry and most of the crystal packing interfaces are wet, but counterexamples exist in both categories. Water molecules at interfaces form hydrogen bonds with protein groups, with a preference for the main-chain carbonyl and the charged side-chains of Glu, Asp, and Arg. These interactions are essentially the same in specific and nonspecific interfaces, and very similar to those observed elsewhere on the protein surface. Water-mediated polar interactions are as abundant at the interfaces as direct protein-protein hydrogen bonds, and they may contribute to the stability of the assembly.     

Functional role of C(alpha)-H...O hydrogen bonds between transmembrane alpha-helices in photosystem I

The crystal structure at 2.5A resolution of the membrane-intrinsic, homotrimeric photosystem I (PSI) isolated from the thermophilic cyanobacterium Synechococcus elongatus shows that each monomer is composed of 12 protein subunits of which nine are embedded in the membrane and feature a total of 34 transmembrane alpha-helices (TMH). Hence, PSI provides an ideal case to study "conventional" and C(alpha)-H...O hydrogen bonds between TMH engaged in intra- and intersubunit interactions. Of the total of 75 C(alpha)-H...O hydrogen bonds between TMHs, 72 are intrasubunit and only three are intersubunit. The two largest subunits PsaA and PsaB are each folded into 11 TMHs showing 29 and 24 intrasubunit C(alpha)-H...O hydrogen bonds, respectively, that are not distributed randomly but many of them flank chlorophyll a (Chl a) co-ordinating amino acids, suggesting stabilisation of these structural segments. As major constituent of the trimerisation domain, subunit PsaL is located next to the 3-fold axis relating the three monomers of PSI. PsaL features a unique number of 19 intrasubunit C(alpha)-H...O hydrogen bonds that connect two of its three TMHs but there are no intersubunit C(alpha)-H...O hydrogen bonds between the three PsaL. Of the three intersubunit C(alpha)-H...O hydrogen bonds, two are formed between PsaA and PsaB and one between PsaB and PsaM. The large number of 75 C(alpha)-H...O hydrogen bonds contrasts the 49 conventional hydrogen bonds, indicating that the former and van der Waals contacts determine association and orientation of TMHs in PSI.     

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.     

Cation-pi interactions in protein-protein interfaces

Arginine is an abundant residue in protein-protein interfaces. The importance of this residue relates to the versatility of its side chain in intermolecular interactions. Different classes of protein-protein interfaces were surveyed for cation-pi interactions. Approximately half of the protein complexes and one-third of the homodimers analyzed were found to contain at least one intermolecular cation-pi pair. Interactions between arginine and tyrosine were found to be the most abundant. The electrostatic interaction energy was calculated to be approximately 3 kcal/mol, on average. A distance-based search of guanidinium:aromatic interactions was also performed using the Macromolecular Structure Database (MSD). This search revealed that half of the guanidinium:aromatic pairs pack in a coplanar manner. Furthermore, it was found that the cationic group of the cation-pi pair is frequently involved in intermolecular hydrogen bonds. In this manner the arginine side chain can participate in multiple interactions, providing a mechanism for inter-protein specificity. Thus, the cation-pi interaction is established as an important contributor to protein-protein interfaces.     

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.     

Anchored design of protein-protein interfaces

Background

Few existing protein-protein interface design methods allow for extensive backbone rearrangements during the design process. There is also a dichotomy between redesign methods, which take advantage of the native interface, and de novo methods, which produce novel binders.

Methodology

Here, we propose a new method for designing novel protein reagents that combines advantages of redesign and de novo methods and allows for extensive backbone motion. This method requires a bound structure of a target and one of its natural binding partners. A key interaction in this interface, the anchor, is computationally grafted out of the partner and into a surface loop on the design scaffold. The design scaffold''s surface is then redesigned with backbone flexibility to create a new binding partner for the target. Careful choice of a scaffold will bring experimentally desirable characteristics into the new complex. The use of an anchor both expedites the design process and ensures that binding proceeds against a known location on the target. The use of surface loops on the scaffold allows for flexible-backbone redesign to properly search conformational space.

Conclusions and Significance

This protocol was implemented within the Rosetta3 software suite. To demonstrate and evaluate this protocol, we have developed a benchmarking set of structures from the PDB with loop-mediated interfaces. This protocol can recover the correct loop-mediated interface in 15 out of 16 tested structures, using only a single residue as an anchor.     

Hydrophobic folding units at protein-protein interfaces: implications to protein folding and to protein-protein association.

A hydrophobic folding unit cutting algorithm, originally developed for dissecting single-chain proteins, has been applied to a dataset of dissimilar two-chain protein-protein interfaces. Rather than consider each individual chain separately, the two-chain complex has been treated as a single chain. The two-chain parsing results presented in this work show hydrophobicity to be a critical attribute of two-state versus three-state protein-protein complexes. The hydrophobic folding units at the interfaces of two-state complexes suggest that the cooperative nature of the two-chain protein folding is the outcome of the hydrophobic effect, similar to its being the driving force in a single-chain folding. In analogy to the protein-folding process, the two-chain, two-state model complex may correspond to the formation of compact, hydrophobic nuclei. On the other hand, the three-state model complex involves binding of already folded monomers, similar to the association of the hydrophobic folding units within a single chain. The similarity between folding entities in protein cores and in two-state protein-protein interfaces, despite the absence of some chain connectivities in the latter, indicates that chain linkage does not necessarily affect the native conformation. This further substantiates the notion that tertiary, non-local interactions play a critical role in protein folding. These compact, hydrophobic, two-chain folding units, derived from structurally dissimilar protein-protein interfaces, provide a rich set of data useful in investigations of the role played by chain connectivity and by tertiary interactions in studies of binding and of folding. Since they are composed of non-contiguous pieces of protein backbones, they may also aid in defining folding nuclei.     

The contribution of C alpha-H...O hydrogen bonds to membrane protein stability depends on the position of the amide

Structural analyses of membrane proteins reveal a large number of C(alpha)-H...O contacts between transmembrane helices, presumed to be hydrogen bonds. Recent experiments produced conflicting results for the contribution of such hydrogen bonds to membrane protein stability. An FTIR study estimated an energy of -0.88 kcal/mol for the G79-C(alpha)-H...I76-O hydrogen bond in glycophorin A, whereas a mutagenesis study showed that the A51-C(alpha)-H...T24-O(gamma) hydrogen bond does not stabilize bacteriorhodopsin. Here, we reconcile these results using molecular mechanics calculations and an implicit membrane model (IMM1). With explicit hydrogen atoms, the potential energy of the G79-C(alpha)-H...I76-O interaction in GpA ranges from -0.54 to -0.9 kcal/mol and its contribution to stability (effective energy) from -0.49 to -0.83 kcal/mol, depending on the structural model used. The average values of these quantities in GpA-like motifs are similar. In bR, the corresponding numbers for the A51-C(alpha)-H...T24-O(gamma) interaction are +0.15 and +0.32 kcal/mol. The difference results from the different arrangement of the interacting groups and specifically the position of the acceptor with respect to the C(alpha) and N atoms. This conclusion likely applies to soluble proteins as well.     

New measures for estimating surface complementarity and packing at protein-protein interfaces

A number of methods exist that use different approaches to assess geometric properties like the surface complementarity and atom packing at the protein-protein interface. We have developed two new and conceptually different measures using the Delaunay tessellation and interface slice selection to compute the surface complementarity and atom packing at the protein-protein interface in a straightforward manner. Our measures show a strong correlation among themselves and with other existing measures, and can be calculated in a highly time-efficient manner. The measures are discriminative for evaluating biological, as well as non-biological protein-protein contacts, especially from large protein complexes and large-scale structural studies (http://pallab.serc.iisc.ernet.in/nip_nsc).