Do water molecules mediate protein-DNA recognition?

A comprehensive analysis of interfacial water molecules in the structures of 109 unique protein-DNA complexes is presented together with a new view on their role in protein-DNA recognition. Location of interfacial water molecules as reported in the crystal structures and as emerging from a series of molecular dynamics studies on protein-DNA complexes with explicit solvent and counterions, was analyzed based on their acceptor, donor hydrogen bond relationships with the atoms and residues of the macromolecules, electrostatic field calculations and packing density considerations. Water molecules for the purpose of this study have been categorized into four classes: viz. (I) those that contact both the protein and the DNA simultaneously and thus mediate recognition directly; (II) those that contact either the protein or the DNA exclusively via hydrogen bonds solvating each solute separately; (III) those that contact the hydrophobic groups in either the protein or the DNA; and, lastly (IV) those that contact another water molecule. Of the 17,963 crystallographic water molecules under examination, about 6% belong to class I and 76% belong to class II. About three-fourths of class I and class II water molecules are exclusively associated with hydrogen bond acceptor atoms of both protein and DNA. Noting that DNA is polyanionic, it is significant that a majority of the crystallographically observed water molecules as well as those from molecular dynamics simulations should be involved in facilitating binding by screening unfavorable electrostatics. Less than 2% of the reported water molecules occur between hydrogen bond donor atoms of protein and acceptor atoms of DNA. These represent cases where protein atoms cannot reach out to DNA to make favorable hydrogen bond interactions due to packing/structural restrictions and interfacial water molecules provide an extension to side-chains to accomplish hydrogen bonding.     

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.     

Structure-based prediction of DNA target sites by regulatory proteins

Regulatory proteins play a critical role in controlling complex spatial and temporal patterns of gene expression in higher organism, by recognizing multiple DNA sequences and regulating multiple target genes. Increasing amounts of structural data on the protein-DNA complex provides clues for the mechanism of target recognition by regulatory proteins. The analyses of the propensities of base-amino acid interactions observed in those structural data show that there is no one-to-one correspondence in the interaction, but clear preferences exist. On the other hand, the analysis of spatial distribution of amino acids around bases shows that even those amino acids with strong base preference such as Arg with G are distributed in a wide space around bases. Thus, amino acids with many different geometries can form a similar type of interaction with bases. The redundancy and structural flexibility in the interaction suggest that there are no simple rules in the sequence recognition, and its prediction is not straightforward. However, the spatial distributions of amino acids around bases indicate a possibility that the structural data can be used to derive empirical interaction potentials between amino acids and bases. Such information extracted from structural databases has been successfully used to predict amino acid sequences that fold into particular protein structures. We surmised that the structures of protein-DNA complexes could be used to predict DNA target sites for regulatory proteins, because determining DNA sequences that bind to a particular protein structure should be similar to finding amino acid sequences that fold into a particular structure. Here we demonstrate that the structural data can be used to predict DNA target sequences for regulatory proteins. Pairwise potentials that determine the interaction between bases and amino acids were empirically derived from the structural data. These potentials were then used to examine the compatibility between DNA sequences and the protein-DNA complex structure in a combinatorial "threading" procedure. We applied this strategy to the structures of protein-DNA complexes to predict DNA binding sites recognized by regulatory proteins. To test the applicability of this method in target-site prediction, we examined the effects of cognate and noncognate binding, cooperative binding, and DNA deformation on the binding specificity, and predicted binding sites in real promoters and compared with experimental data. These results show that target binding sites for several regulatory proteins are successfully predicted, and our data suggest that this method can serve as a powerful tool for predicting multiple target sites and target genes for regulatory proteins.     

Rearrangement of side-chains in a Zif268 mutant highlights the complexities of zinc finger-DNA recognition.

Structural and biochemical studies of Cys(2)His(2) zinc finger proteins initially led several groups to propose a "recognition code" involving a simple set of rules relating key amino acid residues in the zinc finger protein to bases in its DNA site. One recent study from our group, involving geometric analysis of protein-DNA interactions, has discussed limitations of this idea and has shown how the spatial relationship between the polypeptide backbone and the DNA helps to determine what contacts are possible at any given position in a protein-DNA complex. Here we report a study of a zinc finger variant that highlights yet another source of complexity inherent in protein-DNA recognition. In particular, we find that mutations can cause key side-chains to rearrange at the protein-DNA interface without fundamental changes in the spatial relationship between the polypeptide backbone and the DNA. This is clear from a simple analysis of the binding site preferences and co-crystal structures for the Asp20-->Ala point mutant of Zif268. This point mutation in finger one changes the specificity of the protein from GCG TGG GCG to GCG TGG GC(G/T), and we have solved crystal structures of the D20A mutant bound to both types of sites. The structure of the D20A mutant bound to the GCG site reveals that contacts from key residues in the recognition helix are coupled in complex ways. The structure of the complex with the GCT site also shows an important new water molecule at the protein-DNA interface. These side-chain/side-chain interactions, and resultant changes in hydration at the interface, affect binding specificity in ways that cannot be predicted either from a simple recognition code or from analysis of spatial relationships at the protein-DNA interface. Accurate computer modeling of protein-DNA interfaces remains a challenging problem and will require systematic strategies for modeling side-chain rearrangements and change in hydration.     

Crystal structure of paromomycin docked into the eubacterial ribosomal decoding A site

BACKGROUND: Aminoglycoside antibiotics interfere with translation in both gram-positive and gram-negative bacteria by binding to the tRNA decoding A site of the 16S ribosomal RNA. RESULTS: Crystals of complexes between oligoribonucleotides incorporating the sequence of the ribosomal A site of Escherichia coli and the aminoglycoside paromomycin have been solved at 2.5 A resolution. Each RNA fragment contains two A sites inserted between Watson-Crick pairs. The paromomycin molecules interact in an enlarged deep groove created by two bulging and one unpaired adenines. In both sites, hydroxyl and ammonium side chains of the antibiotic form 13 direct hydrogen bonds to bases and backbone atoms of the A site. In the best-defined site, 8 water molecules mediate 12 other hydrogen bonds between the RNA and the antibiotics. Ring I of paromomycin stacks over base G1491 and forms pseudo-Watson-Crick contacts with A1408. Both the hydroxyl group and one ammonium group of ring II form direct and water-mediated hydrogen bonds to the U1495oU1406 pair. The bulging conformation of the two adenines A1492 and A1493 is stabilized by hydrogen bonds between phosphate oxygens and atoms of rings I and II. The hydrophilic sites of the bulging A1492 and A1493 contact the shallow groove of G=C pairs in a symmetrical complex. CONCLUSIONS: Water molecules participate in the binding specificity by exploiting the antibiotic hydration shell and the typical RNA water hydration patterns. The observed contacts rationalize the protection, mutation, and resistance data. The crystal packing mimics the intermolecular contacts induced by aminoglycoside binding in the ribosome.     

Insights into DNA solvation found in protein-DNA structures

The proteins that bind double-helical DNA present various microenvironments that sense and/or induce signals in the genetic material. The high-resolution structures of protein-DNA complexes reveal the nature of both the microenvironments and the conformational responses in DNA and protein. Complex networks of interactions within the structures somehow tie the protein and DNA together and induce the observed spatial forms. Here we show how the cumulative buildup of amino acid atoms around the sugars, phosphates, and bases in different protein-DNA complexes produces a binding cloud around the double helix and how different types of atoms fill that cloud. Rather than focusing on the principles of molecular binding and recognition suggested by the arrangements of amino acids and nucleotides in the macromolecular complexes, we consider the proteins in contact with DNA as organized solvents. We describe differences in the mix of atoms that come in closest contact with DNA, subtle sequence-dependent features in the microenvironment of the sugar-phosphate backbone, a direct link between the localized buildup of ionic species and the electrostatic potential surfaces of the DNA bases, and sites of atomic buildup above and below the basepair planes that transmit the unique features of the base environments along the chain backbone. The inferences about solvation that can be drawn from the survey provide new stimuli for improvement of nucleic acid force fields and fresh ideas for exploration of the properties of DNA in solution.     

The role of residue 50 and hydration water molecules in homeodomain DNA recognition

We conducted molecular dynamics simulations on several wild-type and mutant homeodomain-DNA complexes to investigate the role of residue 50 in homeodomain-DNA interaction and the behavior of interfacial hydration water. Our results suggest that this residue interacts more favorably with its consensus sequence and thus plays a considerable role in DNA recognition. However, residue 50 was not responsible for DNA recognition alone. Other residues in the vicinity could interact with residue 50 in cooperation upon DNA binding. We also found the lifetime for some water in the protein-DNA interface can be as high as nanoseconds and that a few well-conserved sites for water-mediated hydrogen bonds from protein to DNA are occupied by high-mobility hydrating waters.     

NUCPLOT: a program to generate schematic diagrams of protein-nucleic acid interactions.

Proteins that bind to DNA are found in all areas of genetic activity within the cell. To help understand how these proteins perform their various functions, it is useful to analyse which residues are involved in binding to the DNA and how they interact with the bases and sugar-phosphate backbone of nucleic acids. Here we describe a program called NUCPLOT which can automatically identify these interactions from the 3D atomic coordinates of the complex from a PDB file and generate a plot that shows all the interactions in a schematic manner. The program produces a PostScript output file representing hydrogen, van der Waals and covalent bonds between the protein and the DNA. The resulting diagram is both clear and simple and allows immediate identification of important interactions within the structure. It also facilitates comparison of binding found in different structures. NUCPLOT is a completely automatic program, which can be used for any protein-DNA complex and will also work for certain protein-RNA structures.     

Hydrogen and hydration in proteins

Neutron diffraction provides an experimental method of directly locating hydrogen atoms in proteins. High-resolution neutron diffractometers dedicated to biological macromolecules (BIX-type diffractometer) have been constructed at the Japan Atomic Energy Research Institute and they have been used in the 1.5A-resolution crystal structure analyses of several proteins. Interesting topics relevant to hydrogen and hydration in proteins, such as (1) the detailed geometry of hydrogen bonds; (2) information regarding hydrogen/deuterium exchange behavior; (3) the acidities of certain H atoms; (4) the role of hydrogen atoms in enzyme mechanisms and thermostability; (5) the location methyl hydrogen atoms; and (6) dynamical behavior of hydration structures that include H positions have been extracted from these structural results. In addition, a method for the systematic growth of large single crystals based on phase diagrams has been introduced and will be briefly described in this article.     

Modeling helix-turn-helix protein-induced DNA bending with knowledge-based distance restraints.

A crucial element of many gene functions is protein-induced DNA bending. Computer-generated models of such bending have generally been derived by using a presumed bending angle for DNA. Here we describe a knowledge-based docking strategy for modeling the structure of bent DNA recognized by a major groove-inserting alpha-helix of proteins with a helix-turn-helix (HTH) motif. The method encompasses a series of molecular mechanics and dynamics simulations and incorporates two experimentally derived distance restraints: one between the recognition helix and DNA, the other between respective sites of protein and DNA involved in chemical modification-enabled nuclease scissions. During simulation, a DNA initially placed at a distance was "steered" by these restraints to dock with the binding protein and bends. Three prototype systems of dimerized HTH DNA binding were examined: the catabolite gene activator protein (CAP), the phage 434 repressor (Rep), and the factor for inversion stimulation (Fis). For CAP-DNA and Rep-DNA, the root mean square differences between model and x-ray structures in nonhydrogen atoms of the DNA core domain were 2.5 A and 1.6 A, respectively. An experimental structure of Fis-DNA is not yet available, but the predicted asymmetrical bending and the bending angle agree with results from a recent biochemical analysis.     

Quantitative parameters for amino acid-base interaction: implications for prediction of protein-DNA binding sites.

Inspection of the amino acid-base interactions in protein-DNA complexes is essential to the understanding of specific recognition of DNA target sites by regulatory proteins. The accumulation of information on protein-DNA co-crystals challenges the derivation of quantitative parameters for amino acid-base interaction based on these data. Here we use the coordinates of 53 solved protein-DNA complexes to extract all non-homologous pairs of amino acid-base that are in close contact, including hydrogen bonds and hydrophobic interactions. By comparing the frequency distribution of the different pairs to a theoretical distribution and calculating the log odds, a quantitative measure that expresses the likelihood of interaction for each pair of amino acid-base could be extracted. A score that reflects the compatibility between a protein and its DNA target can be calculated by summing up the individual measures of the pairs of amino acid-base involved in the complex, assuming additivity in their contributions to binding. This score enables ranking of different DNA binding sites given a protein binding site and vice versa and can be used in molecular design protocols. We demonstrate its validity by comparing the predictions using this score with experimental binding results of sequence variants of zif268 zinc fingers and their DNA binding sites.     

Protein-DNA interactions: A structural analysis

A detailed analysis of the DNA-binding sites of 26 proteins is presented using data from the Nucleic Acid Database (NDB) and the Protein Data Bank (PDB). Chemical and physical properties of the protein-DNA interface, such as polarity, size, shape, and packing, were analysed. The DNA-binding sites shared common features, comprising many discontinuous sequence segments forming hydrophilic surfaces capable of direct and water-mediated hydrogen bonds. These interface sites were compared to those of protein-protein binding sites, revealing them to be more polar, with many more intermolecular hydrogen bonds and buried water molecules than the protein-protein interface sites. By looking at the number and positioning of protein residue-DNA base interactions in a series of interaction footprints, three modes of DNA binding were identified (single-headed, double-headed and enveloping). Six of the eight enzymes in the data set bound in the enveloping mode, with the protein presenting a large interface area effectively wrapped around the DNA.A comparison of structural parameters of the DNA revealed that some values for the bound DNA (including twist, slide and roll) were intermediate of those observed for the unbound B-DNA and A-DNA. The distortion of bound DNA was evaluated by calculating a root-mean-square deviation on fitting to a canonical B-DNA structure. Major distortions were commonly caused by specific kinks in the DNA sequence, some resulting in the overall bending of the helix. The helix bending affected the dimensions of the grooves in the DNA, allowing the binding of protein elements that would otherwise be unable to make contact. From this structural analysis a preliminary set of rules that govern the bending of the DNA in protein-DNA complexes, are proposed.     

Hydration heat capacity of nucleic acid constituents determined from the random network model

The heat capacities of hydration (dCp) of the five nucleic acid bases A, G, C, T, and U, the sugars ribose and deoxyribose, and the phosphate backbone were determined using Monte Carlo simulations and the random network model. Solute-induced changes in the mean length and root mean square angle of hydrogen bonds between hydration shell waters were used to compute dCp for these solutes. For all solutes the dCp is significantly more positive than predicted from accessible surface area (ASA) models of heat capacity. In ASA models, nitrogen, oxygen, and phosphorus atoms are considered as uniformly polar, therefore making a negative contribution to dCp. However, the simulations show that many of these polar atoms are hydrated by water whose hydrogen bonds are less distorted than in bulk, leading to a positive dCp. This is in contrast to the effect of polar groups seen previously in small molecules and amino acids, which increase the water H-bond distortion, giving negative dCp contributions. Our results imply that dCp accompanying DNA dehydration in DNA-ligand and DNA-protein binding reactions may be significantly more negative than previously believed and that dehydration is a significant contributor to the large decrease in heat capacity seen in experiments.     

Binding of the estrogen receptor to DNA. The role of waters.

Molecular dynamics simulations are carried out to investigate the binding of the estrogen receptor, a member of the nuclear hormone receptor family, to specific and non-specific DNA. Two systems have been simulated, each based on the crystallographic structure of a complex of a dimer of the estrogen receptor DNA binding domain with DNA. One structure includes the dimer and a consensus segment of DNA, ds(CCAGGTCACAGTGACCTGG); the other structure includes the dimer and a nonconsensus segment of DNA, ds(CCAGAACACAGTGACCTGG). The simulations involve an atomic model of the protein-DNA complex, counterions, and a sphere of explicit water with a radius of 45 A. The molecular dynamics package NAMD was used to obtain 100 ps of dynamics for each system with complete long-range electrostatic interactions. Analysis of the simulations revealed differences in the protein-DNA interactions for consensus and nonconsensus sequences, a bending and unwinding of the DNA, a slight rearrangement of several amino acid side chains, and inclusion of water molecules at the protein-DNA interface region. Our results indicate that binding specificity and stability is conferred by a network of direct and water mediated protein-DNA hydrogen bonds. For the consensus sequence, the network involves three water molecules, residues Glu-25, Lys-28, Lys-32, Arg-33, and bases of the DNA. The binding differs for the nonconsensus DNA sequence in which case the fluctuating network of hydrogen bonds allows water molecules to enter the protein-DNA interface. We conclude that water plays a role in furnishing DNA binding specificity to nuclear hormone receptors.     

The role of the minor groove substituents in indirect readout of DNA sequence by 434 repressor

The sequence of non-contacted bases at the center of the 434 repressor binding site affects the strength of the repressor-DNA complex by influencing the structure and flexibility of DNA (Koudelka, G. B., and Carlson, P. (1992) Nature 355, 89-91). We synthesized 434 repressor binding sites that differ in their central sequence base composition to test the importance of minor groove substituents and/or the number of base pair hydrogen bonds between these base pairs on DNA structure and strength of the repressor-DNA complex. We show here that the number of base pair H-bonds between the central bases apparently has no role in determining the relative affinity of a DNA site for repressor. Instead we find that the affinity of DNA for repressor depends on the absence or presence the N2-NH(2) group on the purine bases at the binding site center. The N2-NH(2) group on bases at the center of the 434 binding site appears to destabilize 434 repressor-DNA complexes by decreasing the intimacy of the specific repressor-DNA contacts, while increasing the reliance on protein contacts to the DNA phosphate backbone. Thus, the presence of an N2-NH(2) group on the purines at the center of a binding site globally alters the precise conformation of the protein-DNA interface.     

Hydration of the DNA bases is local.

Ordered hydration sites were determined for the nucleotide bases in B-type conformations using the crystal structure data on 14 B-DNA decamer structures. A method of density representation was extended so that positions, occupancies, and distributions of the hydration sites were predicted around a B-DNA double helix by a method analogous to crystallographic refinement. The predicted hydration sites correctly reproduce the main features of hydration around the B-DNA dodecamer. In contrast to the previous observations, the newly available crystal data show the same extent of hydration of guanine and adenine, and of cytosine and thymine.