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.     

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.     

Rational design of DNA sequence-specific zinc fingers

We developed a rational scheme for designing DNA binding proteins. The scheme was applied for a zinc finger protein and the designed sequences were experimentally characterized with high DNA sequence specificity. Starting with the backbone of a known finger structure, we initially calculated amino acid sequences compatible with the expected structure and the secondary structures of the designed fingers were then experimentally confirmed. The DNA-binding function was added to the designed finger by reconsidering a section of the amino acid sequence and computationally selecting amino acids to have the lowest protein-DNA interaction energy for the target DNA sequences. Among the designed proteins, one had a gap between the lowest and second lowest protein-DNA interaction energies that was sufficient to give DNA sequence-specificity.     

Protein-DNA interactions: amino acid conservation and the effects of mutations on binding specificity

We investigate the conservation of amino acid residue sequences in 21 DNA-binding protein families and study the effects that mutations have on DNA-sequence recognition. The observations are best understood by assigning each protein family to one of three classes: (i) non-specific, where binding is independent of DNA sequence; (ii) highly specific, where binding is specific and all members of the family target the same DNA sequence; and (iii) multi-specific, where binding is also specific, but individual family members target different DNA sequences. Overall, protein residues in contact with the DNA are better conserved than the rest of the protein surface, but there is a complex underlying trend of conservation for individual residue positions. Amino acid residues that interact with the DNA backbone are well conserved across all protein families and provide a core of stabilising contacts for homologous protein-DNA complexes. In contrast, amino acid residues that interact with DNA bases have variable levels of conservation depending on the family classification. In non-specific families, base-contacting residues are well conserved and interactions are always found in the minor groove where there is little discrimination between base types. In highly specific families, base-contacting residues are highly conserved and allow member proteins to recognise the same target sequence. In multi-specific families, base-contacting residues undergo frequent mutations and enable different proteins to recognise distinct target sequences. Finally, we report that interactions with bases in the target sequence often follow (though not always) a universal code of amino acid-base recognition and the effects of amino acid mutations can be most easily understood for these interactions.     

Site-directed mutagenesis studies of tn5 transposase residues involved in synaptic complex formation

Transposition (the movement of discrete segments of DNA, resulting in rearrangement of genomic DNA) initiates when transposase forms a dimeric DNA-protein synaptic complex with transposon DNA end sequences. The synaptic complex is a prerequisite for catalytic reactions that occur during the transposition process. The transposase-DNA interactions involved in the synaptic complex have been of great interest. Here we undertook a study to verify the protein-DNA interactions that lead to synapsis in the Tn5 system. Specifically, we studied (i) Arg342, Glu344, and Asn348 and (ii) Ser438, Lys439, and Ser445, which, based on the previously published cocrystal structure of Tn5 transposase bound to a precleaved transposon end sequence, make cis and trans contacts with transposon end sequence DNA, respectively. By using genetic and biochemical assays, we showed that in all cases except one, each of these residues plays an important role in synaptic complex formation, as predicted by the cocrystal structure.     

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.     

Biochemical and genomic analysis of substrate recognition by the double-stranded RNA binding domain of yeast RNase III

Members of the RNase III family of double-stranded RNA (dsRNA) endonucleases are important enzymes of RNA metabolism in eukaryotic cells. Rnt1p is the only known member of the RNase III family of endonucleases in Saccharomyces cerevisiae. Previous studies have shown that Rnt1p cleaves dsRNA capped by a conserved AGNN tetraloop motif, which is a major determinant for Rnt1p binding and cleavage. The solution structure of the dsRNA-binding domain (dsRBD) of Rnt1p bound to a cognate RNA substrate revealed the structural basis for binding of the conserved tetraloop motif by alpha-helix 1 of the dsRBD. In this study, we have analyzed extensively the effects of mutations of helix 1 residues that contact the RNA. We show, using microarray analysis, that mutations of these amino acids induce substrate-specific processing defects in vivo. Cleavage kinetics and binding studies show that these mutations affect RNA cleavage and binding in vitro to different extents and suggest a function for some specific amino acids of the dsRBD in the catalytic positioning of the enzyme. Moreover, we show that 2'-hydroxyl groups of nucleotides of the tetraloop or adjacent base pairs predicted to interact with residues of alpha-helix 1 are important for Rnt1p cleavage in vitro. This study underscores the importance of a few amino acid contacts for positioning of a dsRBD onto its RNA target, and implicates the specific orientation of helix 1 on the RNA for proper positioning of the catalytic domain.     

Influence of DNA stiffness in protein-DNA recognition

Protein-DNA recognition plays an essential role in the regulation of gene expression. The protein-DNA binding specificity is based on direct atomic contacts between protein and DNA and/or the conformational properties of DNA. In this work, we have analyzed the influence of DNA stiffness (E) to the specificity of protein-DNA complexes. The average DNA stiffness parameters for several protein-DNA complexes have been computed using the structure based sequence dependent stiffness scale. The relationship between DNA stiffness and experimental protein-DNA binding specificity has been brought out. We have investigated the importance of DNA stiffness with the aid of experimental free energy changes (DeltaDeltaG) due to binding in several protein-DNA complexes, such as, ETS proteins, 434, lambda, Mnt and trp repressors, 434 cro protein, EcoRV endonuclease V and zinc fingers. We found a correlation in the range 0.65-0.97 between DeltaDeltaG and E in these examples. Further, we have qualitatively analyzed the effect of mutations in the target sequence of lambda repressor and we observed that the DNA stiffness could correctly identify 70% of the correct bases among the considered nine positions.     

Human herpesvirus 6B origin-binding protein: DNA-binding domain and consensus binding sequence.

We previously demonstrated by a DNA-binding assay that the human herpesvirus 6B (HHV-6B) replication origin has a structure similar to those of alphaherpesviruses, although the HHV-6B and herpes simplex virus type 1 (HSV-1) origin-binding proteins (OBPs) and origins are not interchangeable. Here we describe additional properties of the interaction between HHV-6B OBP and the HHV-6B origin. Competitive electrophoretic mobility shift assays (EMSAs) with DNA duplexes containing single-base alterations allowed deduction of a consensus DNA sequence for HHV-6B-specific OBP binding, YGWYCWCCY, where Y is T or C and W is T or A, while that for HSV-1-specific binding was reported to be YGYTCGCACT. By EMSA, the HHV-6B OBP DNA-binding domain was mapped to a segment containing amino acids 482 to 770. However, in Southwestern (protein-DNA) blotting, the region sufficient for the DNA binding encompassed only amino acids 657 to 770. Similarly, Southwestern blotting showed that amino acids 689 to 851 of HSV-1 OBP had HSV-1 origin-binding activity, although this region was insufficient for origin binding in the EMSA. Although the longer DNA-binding domains identified by EMSA have marginal overall homology among HHV-6B and alphaherpesvirus OBP homologs, the smaller regions sufficient for the binding observed by Southwestern blotting have significant similarity. From these results, we propose a hypothesis that the DNA-binding domain of herpesvirus OBPs consists of two subdomains, one containing a conserved motif that contacts DNA directly, and another, less well conserved, that may modulate either the conformation or accessibility of the binding domain.     

The M.EcoRV DNA-(adenine N6)-methyltransferase uses DNA bending for recognition of an expanded EcoDam recognition site

The M.EcoRV DNA methyltransferase recognizes GATATC sites. It is related to EcoDam, which methylates GATC sites. The DNA binding domain of M.EcoRV is similar to that of EcoDam suggesting a similar mechanism of DNA recognition. We show that amino acid residue Lys11 of M.EcoRV is involved in recognition of Gua1 and Arg128 contacts the Gua in base pair 6. These residues correspond to Lys9 and Arg124 in EcoDam, which recognize the Gua residues in both strands of the Dam recognition sequence, indicating that M.EcoRV and EcoDam make similar contacts to outermost base pairs of their recognition sequences and M.EcoRV recognizes its target site as an expanded GATC site. In contrast to EcoDam, M.EcoRV considerably bends the DNA (59+/-4 degrees) suggesting indirect readout of the AT-rich inner sequence. Recognition of an expanded target site by DNA bending is a new principle for changing DNA recognition specificity of proteins during molecular evolution. R128A is inefficient in DNA bending and binding, whereas K11A bends DNA with relaxed sequence specificity. These results suggest a temporal order of the formation of protein-DNA contacts in which the Gua6-Arg128 contact forms early followed by DNA bending and, finally, the formation of the Lys11-Gua1 contact.