Amino acid–base interactions: a three-dimensional analysis of protein–DNA interactions at an atomic level

To assess whether there are universal rules that govern amino acid–base recognition, we investigate hydrogen bonds, van der Waals contacts and water-mediated bonds in 129 protein–DNA complex structures. DNA–backbone interactions are the most numerous, providing stability rather than specificity. For base interactions, there are significant base–amino acid type correlations, which can be rationalised by considering the stereochemistry of protein side chains and the base edges exposed in the DNA structure. Nearly two-thirds of the direct read-out of DNA sequences involves complex networks of hydrogen bonds, which enhance specificity. Two-thirds of all protein–DNA interactions comprise van der Waals contacts, compared to about one-sixth each of hydrogen and water-mediated bonds. This highlights the central importance of these contacts for complex formation, which have previously been relegated to a secondary role. Although common, water-mediated bonds are usually non-specific, acting as space-fillers at the protein–DNA interface. In conclusion, the majority of amino acid–base interactions observed follow general principles that apply across all protein–DNA complexes, although there are individual exceptions. Therefore, we distinguish between interactions whose specificities are ‘universal’ and ‘context-dependent’. An interactive Web-based atlas of side chain–base contacts provides access to the collected data, including analyses and visualisation of the three-dimensional geometry of the interactions.     

Protein-RNA interactions: structural analysis and functional classes

A data set of 89 protein-RNA complexes has been extracted from the Protein Data Bank, and the nucleic acid recognition sites characterized through direct contacts, accessible surface area, and secondary structure motifs. The differences between RNA recognition sites that bind to RNAs in functional classes has also been analyzed. Analysis of the complete data set revealed that van der Waals interactions are more numerous than hydrogen bonds and the contacts made to the nucleic acid backbone occur more frequently than specific contacts to nucleotide bases. Of the base-specific contacts that were observed, contacts to guanine and adenine occurred most frequently. The most favored amino acid-nucleotide pairings observed were lysine-phosphate, tyrosine-uracil, arginine-phosphate, phenylalanine-adenine and tryptophan-guanine. The amino acid propensities showed that positively charged and polar residues were favored as expected, but also so were tryptophan and glycine. The propensities calculated for the functional classes showed trends similar to those observed for the complete data set. However, the analysis of hydrogen bond and van der Waal contacts showed that in general proteins complexed with messenger RNA, transfer RNA and viral RNA have more base specific contacts and less backbone contacts than expected, while proteins complexed with ribosomal RNA have less base-specific contacts than the expected. Hence, whilst the types of amino acids involved in the interfaces are similar, the distribution of specific contacts is dependent upon the functional class of the RNA bound.     

Structural Insights into the Dual Strategy of Recognition by Peptidoglycan Recognition Protein,PGRP-S: Structure of the Ternary Complex of PGRP-S with Lipopolysaccharide and Stearic Acid

Peptidoglycan recognition proteins (PGRPs) are part of the innate immune system. The 19 kDa Short PGRP (PGRP-S) is one of the four mammalian PGRPs. The concentration of PGRP-S in camel (CPGRP-S) has been shown to increase considerably during mastitis. The structure of CPGRP-S consists of four protein molecules designated as A, B, C and D forming stable intermolecular contacts, A–B and C–D. The A–B and C–D interfaces are located on the opposite sides of the same monomer leading to the the formation of a linear chain with alternating A–B and C–D contacts. Two ligand binding sites, one at C–D contact and another at A–B contact have been observed. CPGRP-S binds to the components of bacterial cell wall molecules such as lipopolysaccharide (LPS), lipoteichoic acid (LTA), and peptidoglycan (PGN) from both Gram-positive and Gram-negative bacteria. It also binds to fatty acids including mycolic acid of the Mycobacterium tuberculosis (Mtb). Previous structural studies of binary complexes of CPGRP-S with LPS and stearic acid (SA) have shown that LPS binds to CPGRP-S at C–D contact (Site-1) while SA binds to it at the A–B contact (Site-2). The binding studies using surface plasmon resonance showed that LPS and SA bound to CPGRP-S in the presence of each other. The structure determination of the ternary complex showed that LPS and SA bound to CPGRP-S at Site-1 and Site-2 respectively. LPS formed 13 hydrogen bonds and 159 van der Waals contacts (distances ≤4.2 Å) while SA formed 56 van der Waals contacts. The ELISA test showed that increased levels of productions of pro-inflammatory cytokines TNF-α and IFN-γ due to LPS and SA decreased considerably upon the addition of CPGRP-S.     

Incorporation of thio-pseudoisocytosine into triplex-forming peptide nucleic acids for enhanced recognition of RNA duplexes

Peptide nucleic acids (PNAs) have been developed for applications in biotechnology and therapeutics. There is great potential in the development of chemically modified PNAs or other triplex-forming ligands that selectively bind to RNA duplexes, but not single-stranded regions, at near-physiological conditions. Here, we report on a convenient synthesis route to a modified PNA monomer, thio-pseudoisocytosine (L), and binding studies of PNAs incorporating the monomer L. Thermal melting and gel electrophoresis studies reveal that L-incorporated 8-mer PNAs have superior affinity and specificity in recognizing the duplex region of a model RNA hairpin to form a pyrimidine motif major-groove RNA₂–PNA triplex, without appreciable binding to single-stranded regions to form an RNA–PNA duplex or, via strand invasion, forming an RNA–PNA₂ triplex at near-physiological buffer condition. In addition, an L-incorporated 8-mer PNA shows essentially no binding to single-stranded or double-stranded DNA. Furthermore, an L-modified 6-mer PNA, but not pseudoisocytosine (J) modified or unmodified PNA, binds to the HIV-1 programmed −1 ribosomal frameshift stimulatory RNA hairpin at near-physiological buffer conditions. The stabilization of an RNA₂–PNA triplex by L modification is facilitated by enhanced van der Waals contacts, base stacking, hydrogen bonding and reduced dehydration energy. The destabilization of RNA–PNA and DNA–PNA duplexes by L modification is due to the steric clash and loss of two hydrogen bonds in a Watson–Crick-like G–L pair. An RNA₂–PNA triplex is significantly more stable than a DNA₂–PNA triplex, probably because the RNA duplex major groove provides geometry compatibility and favorable backbone–backbone interactions with PNA. Thus, L-modified triplex-forming PNAs may be utilized for sequence-specifically targeting duplex regions in RNAs for biological and therapeutic applications.     

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)₂ and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.     

Probing binding hot spots at protein–RNA recognition sites

We use evolutionary conservation derived from structure alignment of polypeptide sequences along with structural and physicochemical attributes of protein–RNA interfaces to probe the binding hot spots at protein–RNA recognition sites. We find that the degree of conservation varies across the RNA binding proteins; some evolve rapidly compared to others. Additionally, irrespective of the structural class of the complexes, residues at the RNA binding sites are evolutionary better conserved than those at the solvent exposed surfaces. For recognitions involving duplex RNA, residues interacting with the major groove are better conserved than those interacting with the minor groove. We identify multi-interface residues participating simultaneously in protein–protein and protein–RNA interfaces in complexes where more than one polypeptide is involved in RNA recognition, and show that they are better conserved compared to any other RNA binding residues. We find that the residues at water preservation site are better conserved than those at hydrated or at dehydrated sites. Finally, we develop a Random Forests model using structural and physicochemical attributes for predicting binding hot spots. The model accurately predicts 80% of the instances of experimental ΔΔG values in a particular class, and provides a stepping-stone towards the engineering of protein–RNA recognition sites with desired affinity.     

Conformational plasticity of RNA for target recognition as revealed by the 2.15?? crystal structure of a human IgG–aptamer complex

Aptamers are short single-stranded nucleic acids with high affinity to target molecules and are applicable to therapeutics and diagnostics. Regardless of an increasing number of reported aptamers, the structural basis of the interaction of RNA aptamer with proteins is poorly understood. Here, we determined the 2.15 Å crystal structure of the Fc fragment of human IgG1 (hFc1) complexed with an anti-Fc RNA aptamer. The aptamer adopts a characteristic structure fit to hFc1 that is stabilized by a calcium ion, and the binding activity of the aptamer can be controlled many times by calcium chelation and addition. Importantly, the aptamer–hFc1 interaction involves mainly van der Waals contacts and hydrogen bonds rather than electrostatic forces, in contrast to other known aptamer–protein complexes. Moreover, the aptamer–hFc1 interaction involves human IgG-specific amino acids, rendering the aptamer specific to human IgGs, and not crossreactive to other species IgGs. Hence, the aptamer is a potent alternative for protein A affinity purification of Fc-fusion proteins and therapeutic antibodies. These results demonstrate, from a structural viewpoint, that conformational plasticity and selectivity of an RNA aptamer is achieved by multiple interactions other than electrostatic forces, which is applicable to many protein targets of low or no affinity to nucleic acids.     

Protein and drug interactions in the minor groove of DNA

Interactions between proteins, drugs, water and B-DNA minor groove have been analyzed in crystal structures of 60 protein–DNA and 14 drug–DNA complexes. It was found that only purine N3, pyrimidine O2, guanine N2 and deoxyribose O4′ are involved in the interactions, and that contacts to N3 and O2 are most frequent and more polar than contacts to O4′. Many protein contacts are mediated by water, possibly to increase the DNA effective surface. Fewer water-mediated contacts are observed in drug complexes. The distributions of ligands around N3 are significantly more compact than around O2, and distributions of water molecules are the most compact. Distributions around O4′ are more diffuse than for the base atoms but most distributions still have just one binding site. Ligands bind to N3 and O2 atoms in analogous positions, and simultaneous binding to N3 and N2 in guanines is extremely rare. Contacts with two consecutive nucleotides are much more frequent than base–sugar contacts within one nucleotide. The probable reason for this is the large energy of deformation of hydrogen bonds for the one nucleotide motif. Contacts of Arg, the most frequent amino acid ligand, are stereochemically indistinguishable from the binding of the remaining amino acids except asparagine (Asn) and phenylalanine (Phe). Asn and Phe bind in distinct ways, mostly to a deformed DNA, as in the complexes of TATA-box binding proteins. DNA deformation concentrates on dinucleotide regions with a distinct deformation of the δ and backbone torsion angles for the Asn and δ, , ζ and χ for the Phe-contacted regions.     

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.     

PRI-Modeler: extracting RNA structural elements from PDB files of protein-RNA complexes

A complete understanding of protein and RNA structures and their interactions is important for determining the binding sites in protein-RNA complexes. Computational approaches exist for identifying secondary structural elements in proteins from atomic coordinates. However, similar methods have not been developed for RNA, due in part to the very limited structural data so far available. We have developed a set of algorithms for extracting and visualizing secondary and tertiary structures of RNA and for analyzing protein-RNA complexes. These algorithms have been implemented in a web-based program called PRI-Modeler (protein-RNA interaction modeler). Given one or more protein data bank files of protein-RNA complexes, PRI-Modeler analyzes the conformation of the RNA, calculates the hydrogen bond (H bond) and van der Waals interactions between amino acids and nucleotides, extracts secondary and tertiary RNA structure elements, and identifies the patterns of interactions between the proteins and RNAs. This paper presents PRI-Modeler and its application to the hydrogen bond and van der Waals interactions in the most representative set of protein-RNA complexes. The analysis reveals several interesting interaction patterns at various levels. The information provided by PRI-Modeler should prove useful for determining the binding sites in protein-RNA complexes. PRI-Modeler is accessible at http://wilab.inha.ac.kr/primodeler/, and supplementary materials are available in the analysis results section at http://wilab.inha.ac.kr/primodeler/.     

Recognition of DNA sequences by the repressor of bacteriophage 434

The structure of a complex between the DNA-binding domain of phage 434 repressor and a 14 base-pair synthetic DNA operator reveals the molecular interactions important for sequence-specific recognition. A set of contacts with DNA backbone, notably involving hydrogen bonds between peptide-NH groups and DNA phosphates, position the repressor and fix the DNA configuration. Direct interactions between amino acid side chains and DNA bases involve nonpolar van der Waals contacts as well as hydrogen bonds. The structures of the repressor domain and of the 434 cro protein are extremely similar. There appear to be no major conformational changes in the proteins when they bind to DNA.     

Structural basis of replication origin recognition by the DnaA protein

Escherichia coli DnaA binds to 9 bp sequences (DnaA boxes) in the replication origin, oriC, to form a complex initiating chromosomal DNA replication. In the present study, we determined the crystal structure of its DNA-binding domain (domain IV) complexed with a DnaA box at 2.1 Å resolution. DnaA domain IV contains a helix–turn–helix motif for DNA binding. One helix and a loop of the helix– turn–helix motif are inserted into the major groove and 5 bp (3′ two-thirds of the DnaA box sequence) are recognized through base-specific hydrogen bonds and van der Waals contacts with the C5-methyl groups of thymines. In the minor groove, Arg399, located in the loop adjacent to the motif, recognizes three more base pairs (5′ one-third of the DnaA box sequence) by base-specific hydrogen bonds. DNA bending by ~28° was also observed in the complex. These base-specific interactions explain how DnaA exhibits higher affinity for the strong DnaA boxes (R1, R2 and R4) than the weak DnaA boxes (R3 and M) in the replication origin.     

A new hydrogen-bonding potential for the design of protein-RNA interactions predicts specific contacts and discriminates decoys

RNA-binding proteins play many essential roles in the regulation of gene expression in the cell. Despite the significant increase in the number of structures for RNA–protein complexes in the last few years, the molecular basis of specificity remains unclear even for the best-studied protein families. We have developed a distance and orientation-dependent hydrogen-bonding potential based on the statistical analysis of hydrogen-bonding geometries that are observed in high-resolution crystal structures of protein–DNA and protein–RNA complexes. We observe very strong geometrical preferences that reflect significant energetic constraints on the relative placement of hydrogen-bonding atom pairs at protein–nucleic acid interfaces. A scoring function based on the hydrogen-bonding potential discriminates native protein–RNA structures from incorrectly docked decoys with remarkable predictive power. By incorporating the new hydrogen-bonding potential into a physical model of protein–RNA interfaces with full atom representation, we were able to recover native amino acids at protein–RNA interfaces.     

Theoretical study of molecular recognition by Hoechst 33258 derivatives

The factors responsible for the binding of Hoechst 33258 with DNA residues have been investigated in this work using the AM1 method. First and foremost, it is found that, although all crystal structure determinations indicate a preference for binding at AT rich sites, the hydrogen bond strength is actually greater for complexes with cytosine and guanine. From this, it has been inferred that other factors such as electrostatic, van der Waals interactions and nonbonded contacts with the walls of the minor groove have a strong role to play in the binding process. The hydrogen bond is found to be stronger for complexation with the thymine O2 than with the adenine N3, in line with experimental observations. Combined QM/MM studies on the drug complexed with the Dickerson-Drew dodecamer reveal that binding induces structural changes in both the ligand as well as DNA. Electron donating substituents at the para position in the phenyl ring of Hoechst 33258 lead to stronger binding with DNA. A correlation with the octanol/water partition coefficients points to the importance of hydrophobic and electrostatic interactions.     

NMR investigations of protein-carbohydrate interactions binding studies and refined three-dimensional solution structure of the complex between the B domain of wheat germ agglutinin and N,N', N"-triacetylchitotriose.

The specific interaction of the isolated B domain of wheat germ agglutinin (WGA-B) with N,N',N"-triacetylchitotriose has been analyzed by 1H-NMR spectroscopy. The association constants for the binding of WGA-B to this trisaccharide have been determined from both 1H-NMR titration experiments and microcalorimetry methods. Entropy and enthalpy of binding have been obtained. The driving force for the binding process is provided by a negative DeltaH which is partially compensated by negative DeltaS. These negative signs indicate that hydrogen bonding and van der Waals forces are the major interactions stabilizing the complex. NOESY NMR experiments in water solution provided 327 protein proton-proton distance constraints. All the experimental constraints were used in a refinement protocol including restrained molecular dynamics in order to determine the refined solution conformation of this protein/carbohydrate complex. With regard to the NMR structure of the free protein, no important changes in the protein NOEs were observed, indicating that carbohydrate-induced conformational changes are small. The average backbone rmsd of the 35 refined structures was 1.05 A, while the heavy atom rmsd was 2.10 A. Focusing on the bound ligand, two different orientations of the trisaccharide within WGA-B binding site are possible. It can be deduced that both hydrogen bonds and van der Waals contacts confer stability to both complexes. A comparison of the three-dimensional structure of WGA-B in solution to that reported in the solid state and to those deduced for hevein and pseudohevein in solution has also been performed.     

Insights into negative modulation of E. coli replication initiation from the structure of SeqA-hemimethylated DNA complex

The SeqA protein binds clusters of fully methylated or hemimethylated GATC sequences at oriC and negatively modulates the initiation of DNA replication. We find that SeqA can be proteolytically cleaved into an N-terminal multimerization and a C-terminal DNA-binding domain and have determined the crystal structure of the C-terminal domain in complex with a hemimethylated GATC site. SeqA makes direct hydrogen bonds and van der Waals contacts with the hemimethylated A-T base pair in addition to interactions with the surrounding bases and DNA backbone. The tetrameric protein-DNA complex found in the crystal suggests that SeqA binds multiple GATC sites on separate DNA duplexes, altering the overall DNA topology and sequestering oriC from replication initiation.     

The Mu repressor-DNA complex contains an immobilized 'wing' within the minor groove

We have determined the solution structure of the complex between the 'winged-helix' enhancer binding domain of the Mu repressor protein and its cognate DNA site. The structure reveals an unusual use for the 'wing' which becomes immobilized upon DNA binding where it makes intermolecular hydrogen bond contacts deep within the minor groove. Although the wing is mobile in the absence of DNA, it partially negates the large entropic penalty associated with its burial by maintaining a small degree of structural order in the DNA-free state. Extensive contacts are also formed between the recognition helix and the DNA, which reads the major groove of a highly conserved region of the binding site through a single base-specific hydrogen bond and van der Waals contacts.     

Structural and Functional Characterization of an Anesthetic Binding Site in the Second Cysteine-Rich Domain of Protein Kinase Cδ?

Elucidating the principles governing anesthetic-protein interactions requires structural determinations at high resolutions not yet achieved with ion channels. Protein kinase C (PKC) activity is modulated by general anesthetics. We solved the structure of the phorbol-binding domain (C1B) of PKCδ complexed with an ether (methoxymethylcycloprane) and with an alcohol (cyclopropylmethanol) at 1.36-Å resolution. The cyclopropane rings of both agents displace a single water molecule in a surface pocket adjacent to the phorbol-binding site, making van der Waals contacts with the backbone and/or side chains of residues Asn-237 to Ser-240. Surprisingly, two water molecules anchored in a hydrogen-bonded chain between Thr-242 and Lys-260 impart elasticity to one side of the binding pocket. The cyclopropane ring takes part in π-acceptor hydrogen bonds with the amide of Met-239. There is a crucial hydrogen bond between the oxygen atoms of the anesthetics and the hydroxyl of Tyr-236. A Tyr-236-Phe mutation results in loss of binding. Thus, both van der Waals interactions and hydrogen-bonding are essential for binding to occur. Ethanol failed to bind because it is too short to benefit from both interactions. Cyclopropylmethanol inhibited phorbol-ester-induced PKCδ activity, but failed to do so in PKCδ containing the Tyr-236-Phe mutation.