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.     

DBD-Hunter: a knowledge-based method for the prediction of DNA-protein interactions

The structures of DNA-protein complexes have illuminated the diversity of DNA-protein binding mechanisms shown by different protein families. This lack of generality could pose a great challenge for predicting DNA-protein interactions. To address this issue, we have developed a knowledge-based method, DNA-binding Domain Hunter (DBD-Hunter), for identifying DNA-binding proteins and associated binding sites. The method combines structural comparison and the evaluation of a statistical potential, which we derive to describe interactions between DNA base pairs and protein residues. We demonstrate that DBD-Hunter is an accurate method for predicting DNA-binding function of proteins, and that DNA-binding protein residues can be reliably inferred from the corresponding templates if identified. In benchmark tests on approximately 4000 proteins, our method achieved an accuracy of 98% and a precision of 84%, which significantly outperforms three previous methods. We further validate the method on DNA-binding protein structures determined in DNA-free (apo) state. We show that the accuracy of our method is only slightly affected on apo-structures compared to the performance on holo-structures cocrystallized with DNA. Finally, we apply the method to approximately 1700 structural genomics targets and predict that 37 targets with previously unknown function are likely to be DNA-binding proteins. DBD-Hunter is freely available at http://cssb.biology.gatech.edu/skolnick/webservice/DBD-Hunter/.     

DNA deformation energy as an indirect recognition mechanism in protein-DNA interactions

Proteins that bind to specific locations in genomic DNA control many basic cellular functions. Proteins detect their binding sites using both direct and indirect recognition mechanisms. Deformation energy, which models the energy required to bend DNA from its native shape to its shape when bound to a protein, has been shown to be an indirect recognition mechanism for one particular protein, integration host factor (IHF). This work extends the analysis of deformation to two other DNA-binding proteins, CRP and SRF, and two endonucleases, I-Crel and I-Ppol. Known binding sites for all five proteins showed statistically significant differences in mean deformation energy as compared to random sequences. Binding sites for the three DNA-binding proteins and one of the endonucleases had mean deformation energies lower than random sequences. Binding sites for I-Ppol had mean deformation energy higher than random sequences. Classifiers that were trained using the deformation energy at each base pair step showed good cross-validated accuracy when classifying unseen sequences as binders or nonbinders. These results support DNA deformation energy as an indirect recognition mechanism across a wider range of DNA-binding proteins. Deformation energy may also have a predictive capacity for the underlying catalytic mechanism of DNA-binding enzymes     

From Nonspecific DNA–Protein Encounter Complexes to the Prediction of DNA–Protein Interactions

DNA–protein interactions are involved in many essential biological activities. Because there is no simple mapping code between DNA base pairs and protein amino acids, the prediction of DNA–protein interactions is a challenging problem. Here, we present a novel computational approach for predicting DNA-binding protein residues and DNA–protein interaction modes without knowing its specific DNA target sequence. Given the structure of a DNA-binding protein, the method first generates an ensemble of complex structures obtained by rigid-body docking with a nonspecific canonical B-DNA. Representative models are subsequently selected through clustering and ranking by their DNA–protein interfacial energy. Analysis of these encounter complex models suggests that the recognition sites for specific DNA binding are usually favorable interaction sites for the nonspecific DNA probe and that nonspecific DNA–protein interaction modes exhibit some similarity to specific DNA–protein binding modes. Although the method requires as input the knowledge that the protein binds DNA, in benchmark tests, it achieves better performance in identifying DNA-binding sites than three previously established methods, which are based on sophisticated machine-learning techniques. We further apply our method to protein structures predicted through modeling and demonstrate that our method performs satisfactorily on protein models whose root-mean-square Cα deviation from native is up to 5 Å from their native structures. This study provides valuable structural insights into how a specific DNA-binding protein interacts with a nonspecific DNA sequence. The similarity between the specific DNA–protein interaction mode and nonspecific interaction modes may reflect an important sampling step in search of its specific DNA targets by a DNA-binding protein.     

DISPLAR: an accurate method for predicting DNA-binding sites on protein surfaces

Structural and physical properties of DNA provide important constraints on the binding sites formed on surfaces of DNA-targeting proteins. Characteristics of such binding sites may form the basis for predicting DNA-binding sites from the structures of proteins alone. Such an approach has been successfully developed for predicting protein–protein interface. Here this approach is adapted for predicting DNA-binding sites. We used a representative set of 264 protein–DNA complexes from the Protein Data Bank to analyze characteristics and to train and test a neural network predictor of DNA-binding sites. The input to the predictor consisted of PSI-blast sequence profiles and solvent accessibilities of each surface residue and 14 of its closest neighboring residues. Predicted DNA-contacting residues cover 60% of actual DNA-contacting residues and have an accuracy of 76%. This method significantly outperforms previous attempts of DNA-binding site predictions. Its application to the prion protein yielded a DNA-binding site that is consistent with recent NMR chemical shift perturbation data, suggesting that it can complement experimental techniques in characterizing protein–DNA interfaces.     

A nuclear DNA-binding protein expressed during early stages of B cell differentiation interacts with diverse segments within and 3' of the Ig H chain gene cluster.

We have used electrophoretic mobility shift assays (EMSA) to detect B cell lineage-specific nuclear proteins that bind to diverse segments within and 3' of the Ig H chain gene cluster. DNA binding sites include sequences 5' of each of the following C region genes: mu, gamma 1, gamma 2a, epsilon, and alpha. For the most part, these binding sites lie 5' of CH-associated tandem repeats. Binding sites for the same B cell lineage-specific proteins have also been defined in the region 3' of C alpha, close to a recently described B cell-specific enhancer element. Cross-competition of EMSA indicates that the B cell lineage-specific nucleoprotein is indistinguishable from those described previously by others: S alpha-BP and BSAP. Because of the diverse sequences recognized by this protein, we term it NF-HB, B-lineage-specific nuclear factor that binds to Ig H gene segments. EMSA using segments 5' of S gamma 2a (5'S gamma 2a-176) and 3' of C alpha (3' alpha-88) shows multiple binding complexes, two of which are B cell lineage specific. The B cell-specific complex with fastest mobility contains only NF-HB, and the one with slowest mobility contains NF-HB together with a ubiquitous DNA-binding protein(s). The ubiquitous binding protein is different for 5' S gamma 2a-176 and for 3' alpha-88, representing the formation of protein-NF-HB complexes specific for these particular Ig DNA regions. Spleen cells show a single band upon EMSA with either 5'S gamma 2a-176 or 3' alpha-88. Upon LPS stimulation, additional binding complexes of slower mobility were formed resulting in a pattern comparable to those detected in pro-B, pre-B, and B cell lines. We hypothesize that NF-HB may promote physical interactions between the 3' alpha-enhancer and segments of the Ig H gene cluster.     

Enhancement of DNA replication by transcription factors NFI and NFIII/Oct-1 depends critically on the positions of their binding sites in the adenovirus origin of replication

The origin of DNA replication of many human adenoviruses is composed of a highly conserved core origin and an auxiliary region, containing the binding sites for NFI and NFIII/Oct-1. We examined enhancement of DNA replication in vitro by the purified functional DNA-binding domains of NFI (NFI-BD) and NFIII/Oct-1 (the POU domain), using origins in which the positions of the binding sites for these proteins were transposed. Insertion or deletion of two or three base pairs between the core origin and the NFI binding site resulted in a 3-5-fold decrease of stimulation, whereas larger insertions gradually reduced the stimulation further. Mutants in which the NFI binding site was separated approximately one or two helical turns from the core origin by AT-rich sequences could still be stimulated by NFI. In contrast, insertion of two or more base pairs between the NFI and NFIII/Oct-1 binding sites abolished stimulation by NFIII/Oct-1 almost completely. Furthermore, stimulation by this protein was lost when the Ad2 NFIII/Oct-1 binding site was transposed to a position closer to the core origin, destroying the NFI binding site. This shows that the position of the NFIII/Oct-1 binding site is essential for stimulation. Models to explain these position-dependent effects on stimulation are discussed.     

Facilitated diffusion of a DNA binding protein on chromatin.

Facilitated diffusion accounts for the rapid rate of association of many bacterial DNA binding proteins with specific DNA sequences in vitro. In this mechanism the proteins bind at random to non-specific sites on the DAN and diffuse (by 'sliding' or 'hopping') along the DNA chain until they arrive at their specific functional sites. We have investigated whether such a mechanism can operate in chromatin by using a bacterial DNA binding protein, Escherichia coli RNA polymerase, that depends on linear diffusion to locate initiation sites on DNA. We have measured the competition between chromatin and its free DNA for the formation of initiation complexes. Only the short linker segments exposed by the removal of histone H1 are available for interaction with the polymerase, but the sparsely distributed promoter sites on the linker DNA of such a polynucleosome chain are located at the same rate as those on DNA. We conclude that the polymerase is free to migrate between the separate linker DNA segments of a polynucleosome chain to reach a promoter site. This chain thus permits the 'hopping' of proteins between neighboring linker segments in their search for a target site on the accessible DNA.     

Cellular protein interactions with herpes simplex virus type 1 oriS.

The herpes simplex virus type 1 (HSV-1) origin of DNA replication, oriS, contains an AT-rich region and three highly homologous sequences, sites I, II, and III, identified as binding sites for the HSV-1 origin-binding protein (OBP). In the present study, interactions between specific oriS DNA sequences and proteins in uninfected cell extracts were characterized. The formation of one predominant protein-DNA complex, M, was demonstrated in gel shift assays following incubation of uninfected cell extracts with site I DNA. The cellular protein(s) that comprises complex M has been designated origin factor I (OF-I). The OF-I binding site was shown to partially overlap the OBP binding site within site I. Complexes with mobilities indistinguishable from that of complex M also formed with site II and III DNAs in gel shift assays. oriS-containing plasmid DNA mutated in the OF-I binding site exhibited reduced replication efficiency in transient assays, demonstrating a role for this site in oriS function. The OF-I binding site is highly homologous to binding sites for the cellular CCAAT DNA-binding proteins. The binding site for the CCAAT protein CP2 was found to compete for OF-I binding to site I DNA. These studies support a model involving the participation of cellular proteins in the initiation of HSV-1 DNA synthesis at oriS.     

The euryarchaeota, nature's medium for engineering of single-stranded DNA-binding proteins

The architecture of single-stranded DNA-binding proteins, which play key roles in DNA metabolism, is based on different combinations of the oligonucleotide/oligosaccharide binding (OB) fold. Whereas the polypeptide serving this function in bacteria contains one OB fold, the eukaryotic functional homolog comprises a complex of three proteins, each harboring at least one OB fold. Here we show that unlike these groups of organisms, the Euryarchaeota has exploited the potential in the OB fold to re-invent single-stranded DNA-binding proteins many times. However, the most common form is a protein with two OB folds and one zinc finger domain. We created several deletion mutants of this protein based on its conserved motifs, and from these structures functional chimeras were synthesized, supporting the hypothesis that gene duplication and recombination could lead to novel functional forms of single-stranded DNA-binding proteins. Biophysical studies showed that the orthologs of the two OB fold/one zinc finger replication protein A in Methanosarcina acetivorans and Methanopyrus kandleri exhibit two binding modes, wrapping and stretching of DNA. However, the ortholog in Ferroplasma acidarmanus possessed only the stretching mode. Most interestingly, a second single-stranded DNA-binding protein, FacRPA2, in this archaeon exhibited the wrapping mode. Domain analysis of this protein, which contains a single OB fold, showed that its architecture is similar to the functional homologs thought to be unique to the Crenarchaeotes. Most unexpectedly, genes coding for similar proteins were found in the genomes of eukaryotes, including humans. Although the diversity shown by archaeal single-stranded DNA-binding proteins is unparalleled, the presence of their simplest form in many organisms across all domains of life is of greater evolutionary consequence.     

Twelve species of the nucleoid-associated protein from Escherichia coli. Sequence recognition specificity and DNA binding affinity.

The genome of Escherichia coli is composed of a single molecule of circular DNA with the length of about 47,000 kilobase pairs, which is associated with about 10 major DNA-binding proteins, altogether forming the nucleoid. We expressed and purified 12 species of the DNA-binding protein, i.e. CbpA (curved DNA-binding protein A), CbpB or Rob (curved DNA-binding protein B or right arm of the replication origin binding protein), DnaA (DNA-binding protein A), Dps (DNA-binding protein from starved cells), Fis (factor for inversion stimulation), Hfq (host factor for phage Q(beta)), H-NS (histone-like nucleoid structuring protein), HU (heat-unstable nucleoid protein), IciA (inhibitor of chromosome initiation A), IHF (integration host factor), Lrp (leucine-responsive regulatory protein), and StpA (suppressor of td(-) phenotype A). The sequence specificity of DNA binding was determined for all the purified nucleoid proteins using gel-mobility shift assays. Five proteins (CbpB, DnaA, Fis, IHF, and Lrp) were found to bind to specific DNA sequences, while the remaining seven proteins (CbpA, Dps, Hfq, H-NS, HU, IciA, and StpA) showed apparently sequence-nonspecific DNA binding activities. Four proteins, CbpA, Hfq, H-NS, and IciA, showed the binding preference for the curved DNA. From the apparent dissociation constant (K(d)) determined using the sequence-specific or nonspecific DNA probes, the order of DNA binding affinity were determined to be: HU > IHF > Lrp > CbpB(Rob) > Fis > H-NS > StpA > CbpA > IciA > Hfq/Dps, ranging from 25 nM (HU binding to the non-curved DNA) to 250 nM (Hfq binding to the non-curved DNA), under the assay conditions employed.     

Zinc fingers and other domains cooperate in binding of Drosophila sry beta and delta proteins at specific chromosomal sites.

The closely related Drosophila serendipity (sry) beta and delta zinc finger proteins display consensus in vitro DNA recognition sequences differing by 4 of 13 nucleotide positions and bind in vivo to distinct sets of sites on polytene chromosomes. We compared the pattern of in vivo chromosomal binding of deleted forms of the sry delta protein fused to beta-galactosidase and expressed in Drosophila transgenic lines. Results show that the carboxy-terminal DNA-binding finger domain is required and sufficient for binding at specific chromosomal sites but that this binding does not nearly reproduce the wild-type pattern. An NH2-terminal domain of the sry delta protein is essential to its specificity of in vivo interaction with chromatin. In vitro and in vivo experiments using reciprocal finger swap between the sry beta and delta proteins suggest that the in vivo specificity is dependent on selective protein-protein contacts at defined chromosomal sites, in addition to DNA specific recognition.