Interaction of DNA with clusters of amino acids in proteins

Protein-DNA interactions facilitate the fundamental functions of living cells and are universal in all living organisms. Several investigations have been carried out, essentially identifying pairs of interactions between the amino acid residues in proteins and the bases in DNA. In the present study, we have detected the recognition motifs that may constitute a cluster of spatially interacting residues in proteins, which interact with the bases of DNA. Graph spectral algorithm has been used to detect side chain clusters comprising Arg, Lys, Asn, Gln and aromatic residues from proteins interacting with DNA. We find that the interaction of proteins with DNA is through clusters in about half of the proteins in the dataset and through individual residues in the rest. Furthermore, inspection of the clusters has revealed additional interactions in a few cases, which have not been reported earlier. The geometry of the interaction between the DNA base and the protein residue is quantified by the distance d and the angle theta. These parameters have been identified for the cation-pi/H-bond stair motif that was reported earlier. Among the Arg, Lys, Asn and Gln residues, the range of (d, theta) values of the interacting Arg clearly falls into the cation-pi and the hydrogen bond interactions of the 'cation-pi/H-bond' stair motif. Analysis of the cluster composition reveals that the Arg residue is predominant than the Lys, Asn and Gln residues. The clusters are classified into Type I and Type II based on the presence or absence of aromatic residues (Phe, Tyr) in them. Residue conservation in these clusters has been examined. Apart from the conserved residues identified previously, a few more residues mainly Phe, Tyr and Arg have also been identified as conserved and interactive with the DNA. Interestingly, a few residues that are parts of interacting clusters and do not interact directly with the DNA have also been conserved. This emphasizes the importance of recognizing the protein side chain cluster motifs interacting with the DNA, which could serve as signatures of protein-DNA recognition in the families of DNA binding proteins.     

Probing the energetic and structural role of amino acid/nucleobase cation-pi interactions in protein-ligand complexes

X-ray structures of proteins bound to ligand molecules containing a nucleic acid base were systematically searched for cation-pi interactions between the base and a positively charged or partially charged side chain group located above it, using geometric criteria. Such interactions were found in 38% of the complexes and are thus even more frequent than pi-pi stacking interactions. They are moreover well conserved in families of related proteins. The overwhelming majority of cation-pi contacts involve Ade bases, as these constitute by far the most frequent ligand building block; Arg-Ade is the most frequent cation-pi pair. Ab initio energy calculations at MP2 level were performed on all recorded pairs. Though cation-pi interactions involving the net positive charge carried by Arg or Lys side chains are the most favorable energetically, those involving the partial positive charge of Asn and Gln side chain amino groups (sometimes referred to as amino-pi interactions) are favorable too, owing to the electron correlation energy contribution. Chains of cation-pi interactions with a nucleobase bound simultaneously to two charged groups or a charged group sandwiched between two aromatic moieties are found in several complexes. The systematic association of these motifs with specific ligand molecules in unrelated protein sequences raises the question of their role in protein-ligand structure, stability, and recognition.     

Structural analysis of cation-pi interactions in DNA binding proteins

Cation-pi interactions play an important role in the stability of protein structures. In this work, we have analyzed the influence of cation-pi interactions in DNA binding proteins. We observed cation-pi interactions in 45 out of 62 DNA binding proteins and there is no significant correlation between the number of amino acid residues and number of cation-pi interactions. These interactions are mainly formed by long-range contacts, and the role of short and medium-range contacts is minimal. The preference of Arg is higher than Lys to form cation-pi interactions. The pair-wise cation-pi interaction energy between aromatic and positively charged residues shows that Arg-Tyr energy is the strongest among the possible six pairs. The structural analysis of cation-pi interaction forming residues shows that Lys, Trp, and Tyr prefer to be in the binding site of protein-DNA complexes. Further, the accessible surface areas of cation-pi interaction forming cationic residues are significantly less than that of other residues. The preference of cation-pi interaction forming residues in different secondary structures shows that Lys prefers to be in strand and Phe prefers to be in turn regions. The results obtained in the present study will be useful in understanding the contribution of cation-pi interactions to the stability and specificity of protein-DNA complexes.     

Sequence-specific high mobility group box factors recognize 10-12-base pair minor groove motifs

Sequence-specific high mobility group (HMG) box factors bind and bend DNA via interactions in the minor groove. Three-dimensional NMR analyses have provided the structural basis for this interaction. The cognate HMG domain DNA motif is generally believed to span 6-8 bases. However, alignment of promoter elements controlled by the yeast genes ste11 and Rox1 has indicated strict conservation of a larger DNA motif. By site selection, we identify a highly specific 12-base pair motif for Ste11, AGAACAAAGAAA. Similarly, we show that Tcf1, MatMc, and Sox4 bind unique, highly specific DNA motifs of 12, 12, and 10 base pairs, respectively. Footprinting with a deletion mutant of Ste11 reveals a novel interaction between the 3' base pairs of the extended DNA motif and amino acids C-terminal to the HMG domain. The sequence-specific interaction of Ste11 with these 3' base pairs contributes significantly to binding and bending of the DNA motif.     

Stabilization of bound polycyclic aromatic hydrocarbons by a pi-cation interaction

Proteins can use aromatic side-chains to stabilize bound cationic ligands through cation-pi interactions. Here, we report the first example of the reciprocal process, termed pi-cation, in which a cationic protein side-chain stabilizes a neutral aromatic ligand. Site-directed mutagenesis revealed that an arginine side-chain located in the deep binding pocket of a monoclonal antibody (4D5) is essential for binding the neutral polynuclear aromatic hydrocarbon benzo[a]pyrene. This Arg was very likely selected for in the primary response, further underscoring the importance of the pi-cation interaction for ligand binding, which should be considered in protein analysis and design when ligands include aromatic groups.     

Classification of pseudo pairs between nucleotide bases and amino acids by analysis of nucleotide-protein complexes

Nucleotide bases are recognized by amino acid residues in a variety of DNA/RNA binding and nucleotide binding proteins. In this study, a total of 446 crystal structures of nucleotide-protein complexes are analyzed manually and pseudo pairs together with single and bifurcated hydrogen bonds observed between bases and amino acids are classified and annotated. Only 5 of the 20 usual amino acid residues, Asn, Gln, Asp, Glu and Arg, are able to orient in a coplanar fashion in order to form pseudo pairs with nucleotide bases through two hydrogen bonds. The peptide backbone can also form pseudo pairs with nucleotide bases and presents a strong bias for binding to the adenine base. The Watson-Crick side of the nucleotide bases is the major interaction edge participating in such pseudo pairs. Pseudo pairs between the Watson-Crick edge of guanine and Asp are frequently observed. The Hoogsteen edge of the purine bases is a good discriminatory element in recognition of nucleotide bases by protein side chains through the pseudo pairing: the Hoogsteen edge of adenine is recognized by various amino acids while the Hoogsteen edge of guanine is only recognized by Arg. The sugar edge is rarely recognized by either the side-chain or peptide backbone of amino acid residues.     

Structure-based analysis of protein-RNA interactions using the program ENTANGLE

Until recently, drawing general conclusions about RNA recognition by proteins has been hindered by the paucity of high-resolution structures. We have analyzed 45 PDB entries of protein-RNA complexes to explore the underlying chemical principles governing both specific and non-sequence specific binding. To facilitate the analysis, we have constructed a database of interactions using ENTANGLE, a JAVA-based program that uses available structural models in their PDB format and searches for appropriate hydrogen bonding, stacking, electrostatic, hydrophobic and van der Waals interactions. The resulting database of interactions reveals correlations that suggest the basis for the discrimination of RNA from DNA and for base-specific recognition. The data illustrate both major and minor interaction strategies employed by families of proteins such as tRNA synthetases, ribosomal proteins, or RNA recognition motifs with their RNA targets. Perhaps most surprisingly, specific RNA recognition appears to be mediated largely by interactions of amide and carbonyl groups in the protein backbone with the edge of the RNA base. In cases where a base accepts a proton, the dominant amino acid donor is arginine, whereas in cases where the base donates a proton, the predominant acceptor is the backbone carbonyl group, not a side-chain group. This is in marked contrast to DNA-protein interactions, which are governed predominantly by amino acid side-chain interactions with functional groups that are presented in the accessible major groove. RNA recognition often proceeds through loops, bulges, kinks and other irregular structures that permit use of all the RNA functional groups and this is seen throughout the protein-RNA interaction database.     

Contribution of cation-pi interactions to the stability of protein-DNA complexes

Cation-pi interactions between an aromatic ring and a positive charge located above it have proven to be important in protein structures and biomolecule associations. Here, the role of these interactions at the interface of protein-DNA complexes is investigated, by means of ab initio quantum mechanics energy calculations and X-ray structure analyses. Ab initio energy calculations indicate that Na ions and DNA bases can form stable cation-pi complexes, whose binding strength strongly depends on the type of base, on the position of the Na ion, and whether the base is isolated or included in a double-stranded B-DNA. A survey of protein-DNA complex structures using appropriate geometrical criteria revealed cation-pi interactions in 71% of the complexes. More than half of the cation-pi pairs involve arginine residues, about one-third asparagine or glutamine residues that only carry a partial charge, and one-seventh lysine residues. The most frequently observed pair, which is also the most stable as monitored by ab initio energy calculations, is arginine- guanine. Arginine-adenine interactions are also favorable in general, although to a lesser extent, whereas those with thymine and cytosine are not. Our calculations show that the major contribution to cation-pi interactions with DNA bases is of electrostatic nature. These interactions often occur concomitantly with hydrogen bonds with adjacent bases; their strength is estimated to be from three to four times lower than that of hydrogen bonds. Finally, the role of cation-pi interactions in the stability and specificity of protein-DNA complexes is discussed.     

Cation-pi interactions at non-redundant protein-RNA interfaces

Cation-pi interactions have proved to be important in proteins and protein-ligand complexes. Here, cation-pi interactions are analyzed for 282 non-redundant protein-RNA interfaces. The statistical results show that this kind of interactions exists in 65% of the interfaces. The four RNA bases are ranked as Gua > Ade > Ura > Cyt according to their propensity to participate in cation-pi interactions. The corresponding ranking for the involved amino acid residues is: Arg > Lys > Asn > Gln. The same trends are obtained based on the empirical energy calculation. The Arg-Gua pairs have the greatest stability and are also most frequently observed. The number of cation-pi pairs involving unpaired bases is 2.5 times as many as those involving paired bases. Hence, cation-pi interactions show sequence and structural specificities. For the bicyclic bases, Gua and Ade, their 5-atom rings participate in cation-pi interactions somewhat more than the 6-atom rings, with percentages of 54 and 46%, respectively, which is due to the higher cation-pi participation proportion (63%) of 5-atom rings in the paired bases. These results give a general view of cation-pi interactions at protein-RNA interfaces and are helpful in understanding the specific recognition between protein and RNA.     

Cations Form Sequence Selective Motifs within DNA Grooves via a Combination of Cation-Pi and Ion-Dipole/Hydrogen Bond Interactions

The fine conformational subtleties of DNA structure modulate many fundamental cellular processes including gene activation/repression, cellular division, and DNA repair. Most of these cellular processes rely on the conformational heterogeneity of specific DNA sequences. Factors including those structural characteristics inherent in the particular base sequence as well as those induced through interaction with solvent components combine to produce fine DNA structural variation including helical flexibility and conformation. Cation-pi interactions between solvent cations or their first hydration shell waters and the faces of DNA bases form sequence selectively and contribute to DNA structural heterogeneity. In this paper, we detect and characterize the binding patterns found in cation-pi interactions between solvent cations and DNA bases in a set of high resolution x-ray crystal structures. Specifically, we found that monovalent cations (Tl⁺) and the polarized first hydration shell waters of divalent cations (Mg²⁺, Ca²⁺) form cation-pi interactions with DNA bases stabilizing unstacked conformations. When these cation-pi interactions are combined with electrostatic interactions a pattern of specific binding motifs is formed within the grooves.     

The role of lysine 55 in determining the specificity of the purine repressor for its operators through minor groove interactions.

The interaction of the dimeric Escherichia coli purine repressor (PurR) with its cognate sequences leads to a 45 degrees to 50 degrees kink at a central CpG base step towards the major groove, as dyad-related leucine side-chains interdigitate between these bases from the minor groove. The resulting broadening of the minor groove increases the accessibility of the six central base-pairs towards minor groove interactions with residues from PurR. It has been shown that lysine 55 of PurR makes a direct contact with the adenine base (Ade8) directly 5' to the central CpG base-pair step in the high-affinity purF operator sequence. We have investigated the importance of this interaction in the specificity and affinity of wild-type PurR (WT) for its operators and we have studied a mutant of PurR in which Lys55 is replaced with alanine (K55A). Complexes of WT and K55A with duplex DNA containing pur operator sequences varied at position 8 were investigated crystallographically, and binding studies were performed using fluorescence anisotropy. The structures of the protein-DNA complexes reveal a relatively unperturbed global conformation regardless of the identity of the base-pair at position 8 or residue 55. In all structures the combination of higher resolution and a palindromic purF operator site allowed several new PurR.DNA interactions to be observed, including contacts by Thr15, Thr16 and His20. The side-chain of Lys55 makes productive, though varying, interactions with the adenine, thymine or cytosine base at position 8 that result in equilibrium dissociation constants of 2.6 nM, 10 nM and 35 nM, respectively. However, the bulk of the lysine side-chain apparently blocks high-affinity binding of operators with guanine at position 8 (Kd620 nM). Also, the high-affinity binding conformation appears blocked, as crystals of WT bound to DNA with guanine at position 8 could not be grown. In complexes containing K55A, the alanine side-chain is too far removed to engage in van der Waals interactions with the operator, and, with the loss of the general electrostatic interaction between the phosphate backbone and the ammonium group of lysine, K55A binds each operator weakly. However, the mutation leads to a swap of specificity of PurR for the base at position 8, with K55A exhibiting a twofold preference for guanine over adenine. In addition to defining the role of Lys55 in PurR minor groove binding, these studies provide structural insight into the minor groove binding specificities of other LacI/GalR family members that have either alanine (e.g. LacI, GalR, CcpA) or a basic residue (e.g. RafR, ScrR, RbtR) at the comparable position.     

Improved modeling of side-chain--base interactions and plasticity in protein--DNA interface design

Combinatorial sequence optimization for protein design requires libraries of discrete side-chain conformations. The discreteness of these libraries is problematic, particularly for long, polar side chains, since favorable interactions can be missed. Previously, an approach to loop remodeling where protein backbone movement is directed by side-chain rotamers predicted to form interactions previously observed in native complexes (termed "motifs") was described. Here, we show how such motif libraries can be incorporated into combinatorial sequence optimization protocols and improve native complex recapitulation. Guided by the motif rotamer searches, we made improvements to the underlying energy function, increasing recapitulation of native interactions. To further test the methods, we carried out a comprehensive experimental scan of amino acid preferences in the I-AniI protein-DNA interface and found that many positions tolerated multiple amino acids. This sequence plasticity is not observed in the computational results because of the fixed-backbone approximation of the model. We improved modeling of this diversity by introducing DNA flexibility and reducing the convergence of the simulated annealing algorithm that drives the design process. In addition to serving as a benchmark, this extensive experimental data set provides insight into the types of interactions essential to maintain the function of this potential gene therapy reagent.     

Exploring the environmental preference of weak interactions in (alpha/beta)8 barrel proteins

The environmental preference for the occurrence of noncanonical hydrogen bonding and cation-pi interactions, in a data set containing 71 nonredundant (alpha/beta)(8) barrel proteins, with respect to amino acid type, secondary structure, solvent accessibility, and stabilizing residues has been performed. Our analysis reveals some important findings, which include (a) higher contribution of weak interactions mediated by main-chain atoms irrespective of the amino acids involved; (b) domination of the aromatic amino acids among interactions involving side-chain atoms; (c) involvement of strands as the principal secondary structural unit, accommodating cross strand ion pair interaction and clustering of aromatic amino acid residues; (d) significant contribution to weak interactions occur in the solvent exposed areas of the protein; (e) majority of the interactions involve long-range contacts; (f) the preference of Arg is higher than Lys to form cation-pi interaction; and (g) probability of theoretically predicted stabilizing amino acid residues involved in weak interaction is higher for polar amino acids such as Trp, Glu, and Gln. On the whole, the present study reveals that the weak interactions contribute to the global stability of (alpha/beta)(8) TIM-barrel proteins in an environment-specific manner, which can possibly be exploited for protein engineering applications.     

Novel FXXFF and FXXMF motifs in androgen receptor cofactors mediate high affinity and specific interactions with the ligand-binding domain

Upon hormone binding, a hydrophobic coactivator binding groove is induced in the androgen receptor (AR) ligand-binding domain (LBD). This groove serves as high affinity docking site for alpha-helical FXXLF motifs present in the AR N-terminal domain and in AR cofactors. Study of the amino acid requirements at position +4 of the AR FXXLF motif revealed that most amino acid substitutions strongly reduced or completely abrogated AR LBD interaction. Strong interactions were still observed following substitution of Leu+4 by Phe or Met residues. Leu+4 to Met or Phe substitutions in the FXXLF motifs of AR cofactors ARA54 and ARA70 were also compatible with strong AR LBD binding. Like the corresponding FXXLF motifs, interactions of FXXFF and FXXMF variants of AR and ARA54 motifs were AR specific, whereas variants of the less AR-selective ARA70 motif displayed increased AR specificity. A survey of currently known AR-binding proteins revealed the presence of an FXXFF motif in gelsolin and an FXXMF motif in PAK6. In vivo fluorescence resonance energy transfer and functional protein-protein interaction assays showed direct, efficient, and specific interactions of both motifs with AR LBD. Mutation of these motifs abrogated interaction of gelsolin and PAK6 proteins with AR. In conclusion, we have demonstrated strong interaction of FXXFF and FXXMF motifs to the AR coactivator binding groove, thereby mediating specific binding of a subgroup of cofactors to the AR LBD.     

Adenine recognition: a motif present in ATP-, CoA-, NAD-, NADP-, and FAD-dependent proteins.

Adenosine triphosphate (ATP) plays an essential role in energy transfer within the cell. In the form of NAD, adenine participates in multiple redox reactions. Phosphorylation and ATP-hydrolysis reactions have key roles in signal transduction and regulation of many proteins, especially enzymes. In each cell, proteins with many different functions use adenine and its derivatives as ligands; adenine, of course, is present in DNA and RNA. We show that an adenine binding motif, which differs according to the backbone chain direction of a loop that binds adenine (and in one variant by the participation of an aspartate side-chain), is common to many proteins; it was found from an analysis of all adenylate-containing protein structures from the Protein Data Bank. Indeed, 224 protein-ligand complexes (86 different proteins) from a total of 645 protein structure files bind ATP, CoA, NAD, NADP, FAD, or other adenine-containing ligands, and use the same structural elements to recognize adenine, regardless of whether the ligand is a coenzyme, cofactor, substrate, or an allosteric effector. The common adenine-binding motif shown in this study is simple to construct. It uses only (1) backbone polar interactions that are not dependent on the protein sequence or particular properties of amino acid side-chains, and (2) nonspecific hydrophobic interactions. This is probably why so many different proteins with different functions use this motif to bind an adenylate-containing ligand. The adenylate-binding motif reported is present in "ancient proteins" common to all living organisms, suggesting that adenine-containing ligands and the common motif for binding them were exploited very early in evolution. The geometry of adenine binding by this motif mimics almost exactly the geometry of adenine base-pairing seen in DNA and RNA.     

Cation-pi and amino-acceptor interactions between hydrated metal cations and DNA bases. A quantum-chemical view

Cation-pi interactions between cytosine and hexahydrated cations have been characterized using ab initio method with inclusion of electron correlation effects, assuming idealized and crystal geometries of the interacting species. Hydrated metal cations can interact with nucleobases in a cation-pi manner. The stabilization energy of such complexes would be large and comparable to the one for cation-pi complex with benzene. Further, polarized water molecules belonging to the hydration shell of the cation are capable to form a strong hydrogen bond interaction with the nitrogen lone electron pair of the amino groups of bases and enforce a pronounced sp3 pyramidalization of the nucleobase amino groups. However, in contrast to the benzene-cation complexes, the cation-pi configurations are highly unstable for a nucleobase since the conventional in plane binding of hydrated cations to the acceptor sites on the nucleobase is strongly preferred. Thus, a cation-pi interaction with a nucleobase can occur only if the position of the cation is locked above the nucleobase plane by another strong interaction. This indeed can occur in biopolymers and may have an effect on the local DNA architecture. Nevertheless, nucleobases have no intrinsic propensity to form cation-pi interactions.     

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.