Evolutionary conservation of DNA-contact residues in DNA-binding domains
BMC Bioinformatics volume 9, Article number: S3 (2008)
DNA-binding proteins are of utmost importance to gene regulation. The identification of DNA-binding domains is useful for understanding the regulation mechanisms of DNA-binding proteins. In this study, we proposed a method to determine whether a domain or a protein can has DNA binding capability by considering evolutionary conservation of DNA-binding residues.
Our method achieves high precision and recall for 66 families of DNA-binding domains, with a false positive rate less than 5% for 250 non-DNA-binding proteins. In addition, experimental results show that our method is able to identify the different DNA-binding behaviors of proteins in the same SCOP family based on the use of evolutionary conservation of DNA-contact residues.
This study shows the conservation of DNA-contact residues in DNA-binding domains. We conclude that the members in the same subfamily bind DNA specifically and the members in different subfamilies often recognize different DNA targets. Additionally, we observe the co-evolution of DNA-contact residues and interacting DNA base-pairs.
DNA-binding proteins play a key role in living organisms of many genetic activities such as transcription, recombination, DNA replication and repair. One or more domains of these proteins interact with DNA, and they offer the specificity for direct and indirect readout of DNA . To identify the DNA-binding domains is very important for understanding the regulation mechanisms.
Recently, rapidly increasing amount of protein-DNA complexes from X-ray crystallography and nuclear magnetic resonance (NMR) have enabled the use of structural-based approaches for identifying DNA-binding proteins. Most of the structural DNA-binding domains can be categorized into several classes according to their structures or binding types [2–4]. However, some DNA-binding domains can not be well categorized, and for some DNA-binding domains structural information is unavailable [3, 5]. Several studies used various computational approaches to predict potential DNA-binding proteins by using protein-DNA complexes structure features, such as the overall charges, electric moments, and shape of binding sites [6–12]. Since the charge and conformational complementarities of binding sites are essential for protein-DNA binding, these features provide a reasonable basis to identify DNA-binding proteins. Another trend is to consider the degree of conservation of residues [13–15]. Luscombe and Thornton  have studied 21 families of DNA-binding proteins and showed that those amino acids interacting with the DNA are better conserved than those not interacting with DNA. Stawiski et al.  found that electrostatic patches of DNA-binding proteins have a higher percentage of aromatic and positive residues. According to the general properties of 20 amino acids, they also showed that residues of the patch are conserved at property levels.
In this paper, we propose a structure-based threading method by considering evolutionary conservation of DNA-contact residues in DNA-binding domains to identify DNA-binding domains. We use BLOSUM62 , an evolutionary-based scoring matrix for amino acid substitutions, to measure the degree of conservation of binding residues. Our method can achieve high precision and recall for 66 families of DNA-binding domains, with a false positive rate less than 5% for 250 non-DNA-binding proteins.
Given a query domain, our method identified similar DNA-binding structures or homologous protein sequences from the template library. To evaluate the performance of our method, for each DNA-contact domain (D) in the template library we generated its corresponding positive and negative sets. The members in the positive set contain the domains similar to domain D based on SCOP, while domains in the negative set do not. By applying our method on these two sets, we found that the scores of the domains in the positive set are significantly higher than those of domains in the negative set. We further determined a threshold to achieve high precision and recall. Combining with the threshold, we applied our method on 66 known SCOP families of DNA-binding domains and 250 non-DNA-binding proteins to examine the performance.
Positive and negative set for each contact domain
We collected DNA-binding contact domains from SCOP database, the detail is described in Method. To remove redundant contact domains, domains with highly similar sequences (identity > 90%) are grouped using the NCBI software BLASTCLUST. In each group, the one with the maximal number of contact residues is chosen as the representative domain of a group. For a representative domain R, these protein domains in the same SCOP family are considered as the member of R according to SCOP95 (members whose similarity greater than 95% are excluded). Each member of R was aligned to R using the CE. We define a residue of R as misaligned if it is aligned to a gap. A family member is discarded if more than 20% contact residues of R are misaligned between R and this member. Family members that satisfy the above criteria are considered to be in the positive set. If there are less than five members in the positive set of R, the entire family of R is discarded. We finally yielded 66 representative domains with corresponding positive sets. For each R, we artificially generated 1000 domains to be the negative set. To do this, for each artificial domain, we replicate its residues from R. Then we randomly mutated the residue type of each contact residue of R.
Determining the threshold of similar DNA-binding function of a contact domain
For each representative domain R, each member in the positive and negative sets was scored by the method we developed. Ideally, the scores of domains in the positive set should be on average significantly higher than those of the negative set. We used the Kolmogorov-Smirnov (KS) test to examine the above criterion. The KS test is a nonparametric test to determine if two distributions differ significantly. According to our results, the scores are significantly different for the positive set and the negative set in most domains (97% of 66 sets have a p value less than 0.05).
Further, given a contact domain, we would like to determine a threshold for determining which domains have a similar DNA-binding function. For the two sets (positive and negative) of a representative domain, we separately transform all members' scores to z-scores by
where s is the score of a member, μ is the mean score of the these two sets, and δ is the standard deviation. Figures 1A and 1B show the precision (ratio of the number of retrieved true positive data to all retrieved data) and the recall (ratio of the number of retrieved true positive data to all true positive) with various z-score thresholds, respectively. As shown in Figure 1A, when we set the threshold greater than two, the precisions of using different thresholds are very similar (>90%). If we set the z-score threshold to one, only 60% of families are with high precision. The results imply that larger thresholds will yield higher precisions, but the benefit is limited when the threshold is larger than two. Oppositely, as shown in Figure 1B, larger thresholds will reduce the recall. According to these results, we take the z-score threshold as 2.0 and the domains with a z-score higher than the threshold will be considered as putative DNA-binding domains.
We further apply our method to 250 non-nucleic-acid binding (non-DNA-binding) proteins, which were initially studied by Hobohm and Sander  and further specified by Stawiski et al. . We align all non-redundant contact domains to those non-DNA-binding proteins using CE. Alignments whose z-scores (defined by CE) are greater than 3.7 with the misalign rate of contact residues less than 20% are chosen as non-DNA-binding domains. 177 non-DNA-binding domains pass the constraints among 250 proteins. We applied our method on these non-DNA-binding domains and transformed their scores to z-scores. Figure 2 shows the distribution of z-scores of non-DNA-binding domains. The scores approximately follow a normal distribution and the peak of the density occurred at Z = -1~0. Given a z-score threshold, the false positive rate is the ratio of number of domains whose z-score are beyond the threshold to the total non-DNA-binding domains. According to our previous analysis, we set the threshold to 2.0 and the false positive rate is less than 0.05. It shows that for non-DNA-binding domains, our method can recognize their non-binding with high accuracy.
Figure 3A shows an example, which is the ultrabithorax homeodomain (Ubx) from Drosophila melanogaster (PDB entry 1B8I-A) selected from 66 representative domains to described the characteristics of our method. The DNA is represented in green. 18 DNA-contact residues are presented as yellow stick and other residues are denoted as blue. The protein sequence is also presented and a contact residue is marked with an asterisk. For the alignment of the representative domain (1B8I-A) to the domains of its member, Figure 3B presents a nice case (PDB entry 1PUF-A), which is a homeobox protein hox-a9 from mouse . We found that the contact residues is highly conserved in the aligned amino acids of the two domains and our scoring method shows this high z-score (z-score = 11.92). On the other hand, if we align 1B8I-A to 250 non-DNA-binding proteins, our method is able to discard the similar protein structures whose contact residues are not conserved (z-score = 0.58). Figure 3C shows an example of aligning 1B8I-A to 1BOB, which is histone acetyltransferase hat1 from S. cerevisiae in complex with acetyl coenzyme .
The z-score of DNA-binding domains in the same SCOP family may be variable for several representative domains (Figure 4A). The 1PUF-A and 1O4X-A1 (Oct-1 POU homeodomains from Homo sapiens ) are the members of the 1B8I-A representative domain. The z-scores are 11.92 (1PUF-A) and 4.4 (1O4X-A1) when 1B8I-A was used as the query (Figure 4A). The z-scores indicated that the contact residues between 1PUF-A and 1B8I-A are more conserved than the ones between 1O4X-A1 and 1B8I-A on contact residues interacting to the bases of the core binding site in the DNA.
To investigate variation of contact residues of DNA-binding domain in the same SCOP family, we compared the bound DNA sequences of two DNA-binding domains by aligning the double-strand sequences to each other. 1B8I-A binds two DNA sequences (i.e. PDB entry 1B8I-C and 1B8I-D) and 1O4X-A1 binds another two DNA sequences (PDB entry 1O4X-C and 1O4X-D). First we generated four pairing alignments: 1B8I-C and 1O4X-C; 1B8I-C and 1O4X-D; 1B8I-D and 1O4X-C; and 1B8I-D and 1O4X-D. We do not allow any gap insertion when aligning a-pairing DNA sequences. The alignments are obtained by sliding two sequences against each other until the best match is found. The alignment with the maximum number of identical aligned pairs is chosen, and as a result the alignment between 1B8I-C and 1O4X-C is the one chosen (Figure 4C). Then we adjust the alignment of the other DNA strand pairs (i.e. 1B8I-D and 1O4X-D) according to this best alignment (1B8I-C and 1O4X-C).
Figures 4B and 4C show that the number of identical nucleotides between 1B8I-C and 1PUF-E (10) as well as 1B8I-D and 1PUF-D (10) is much higher than those of 1B8I-C and 1O4X-C (6) as well as 1B8I-D and 1O4X-D (5) for whole DNA sequences. At the same time, 11 identical contact nucleotides are obtained from the alignments of 1B8I-C and 1PUF-E as well as 1B8I-D and 1PUF-D; but two identical contact nucleotides are yielded from the alignments of 1B8I-C and 1O4X-C as well as 1B8I-D and 1O4X-D (the contact nucleotides are the nucleotides that interact with contact residues of protein). With respect to 1B8I-A, 1PUF-A and 1O4X-A1 are different not only in the DNA sequences they bind to but also in their DNA-binding sites. These results show that the members in the same SCOP family may have different DNA-binding models and that our method is able to detect the different Protein-DNA interactions based on the evolutionary conservation of DNA-contact residues.
We produced multiple protein sequence alignments of 13 homeodomains (Figure 5) selected from SCOP 1.71 using a multiple structure alignment tool, MUSTANG . These domains were ranked by z-scores calculated by using our scoring method and the sequence of 1B8I-A as the query. According to z-scores, these 13 domains can be roughly divided into two groups, including the Ubx-like homeodomain colored in blue (e.g. PDB entry 9ANT-A (12.77), 1AHD-P (12.19), and 1SAN (11.96)) and the Oct-1 POU homeodomain colored in red (e.g. PDB entry 1E3O-C1 (6.40), 1GT0-C1 (6.38), and 1O4X-A1 (4.40)). Figure 5 shows that all Ubx-like homeodomains are significantly more conserved than Oct-1 POU homeodomains on contact residues (green). The Ubx homeodomain binds together with the extradenticle homeodomain (Exd) to recognize four DNA bases (ATAA)  based on four residues that are Ile47, Gln50, Asn51, and Met54, locating at α 3 helix in the Ubx (gray columns in Figure 5). The z-scores of the domains are higher if they are conserved on these four residues, such as three antennapedia homeodomains and two homeobox protein hox. These results show that contact residues interacting with bases in the DNA sequences are often conserved. This result is consistent to previous results  in which the homeodomain family was considered as a multi-specific family that consists of some subfamilies. This work concluded that members in the same subfamily bind DNA specifically but the members in different subfamilies recognize different DNA targets. In summary, we demonstrated the conservation of DNA-contact residues in DNA-binding domains.
The contact residues of DNA-binding domains are useful in discriminating DNA-binding domains from non-DNA-binding domains in a novel protein sequence. Our method, which considers evolutionary conservation of DNA-binding residues, can achieve high precision and recall for 66 families of DNA-binding domains, with a false positive rate less than 5% for 250 non-DNA-binding proteins. In addition, our method is able to identify the different DNA-binding behaviors of proteins in the same SCOP family based on the evolutionary conservation of DNA-contact residues. We also discussed the mutation of contact residues of DNA-binding domains can possibly change the bound DNA sequences. It implies that the co-change of DNA-contact residues and their DNA-binding bases.
Figure 6 shows the flowchart of our proposed method. We quantitatively evaluated whether a given protein domain M has a similar DNA-binding function to a known crystal protein-DNA structure. For each crystal structure of protein-DNA complex in Protein Data Bank (PDB), we first identified the DNA-contact domain (D) using geometry information and domain definitions from Structure Classification of Proteins (SCOP, version 1.71) . The structures and sequences of both protein-DNA complexes and their DNA-contact domains were collected in the template library. For a given protein sequence/structure M, we used sequence/structural alignment tools to find the homologous DNA-contact domain D from the template library. Finally, we proposed a score method to evaluate the similarity between M and D based on the BLOSUM matrix. Detailed descriptions are as follows.
We first collected protein-DNA complexes from PDB and each complex should contain at least one protein chain and a double-strand DNA. As in Luscombe et al. , a complex was excluded if its DNA is single-stranded or the length of the DNA is less than 4 bases. For each protein-DNA complex, we then identify contact residues and contact domains of this protein. Contact residues, whose heavy atoms are within a distance (distance ≤ 4.5 Å) of any heavy atoms of the bound DNA, are considered as the core parts of the contact domain in a complex . For each protein-DNA complex, we identified its DNA-contact domains according to contact residues and the definition of the SCOP database. Each domain must have more than 5 contact residues and the number of residues of this protein is more than 50 to make sure that the contact between the protein and DNA was reasonably extensive. Finally, 230 contact DNA-binding domains were identified and collected in the template library.
Homologous proteins searching
For a given protein sequence/structure M, we found a homologous DNA-binding protein from the template library using alignment tools. If M is a 3D-structure, we used a structure alignment (i.e. CE ) to align M to all contact domains. The CE will return a Z score for each alignment representing the structure similarity of the two aligned structures. DNA-binding proteins are considered as homologous proteins of query M if CE Z scores of exceed 3.7 based on CE's statistical model. On the other hand, if M is a protein sequence, we used sequence alignment (i.e. FASTA [29–31]) to search the template library. Here, a DNA-binding protein is considered a homologous protein of M if the sequence identity exceeds 25% according to observations of previous studies [32–37].
For an alignment of the query domain (M) and a contact domain (D) that satisfies the above criterion, we calculate the alignment score for the aligned contact residues by using the BLOSUM62 matrix. The scoring method is defined as:
where CR is the set of the contact residues between D and M; d i and m i denote the corresponding ithcontact residue of D and M, respectively. Here, the score of a misaligned residue is -4 which is the smallest in the BLOSUM62 matrix.
Authors' contributionsYLC and HKT carried out the design of scoring functions and data set preparation, participated in experimental designs and drafted the manuscript. CYK provided the design of this study. YCC and YJH provided the domain knowledge and useful comments. JMY provided the original idea, participated in the design and coordination of this study and helped to draft the manuscript. All authors read and approved the final manuscript.
Michael Gromiha M, Siebers JG, Selvaraj S, Kono H, Sarai A: Intermolecular and intramolecular readout mechanisms in protein-DNA recognition. J Mol Biol 2004,337(2):285–294. 10.1016/j.jmb.2004.01.033
Vinson CR, Sigler PB, McKnight SL: Scissors-grip model for DNA recognition by a family of leucine zipper proteins. Science 1989,246(4932):911–916. 10.1126/science.2683088
Harrison SC: A structural taxonomy of DNA-binding domains. Nature 1991,353(6346):715–719. 10.1038/353715a0
Luscombe NM, Austin SE, Berman HM, Thornton JM: An overview of the structures of protein-DNA complexes. Genome Biol 2000,1(1):REVIEWS001. 10.1186/gb-2000-1-1-reviews001
Johnson PF, McKnight SL: Eukaryotic transcriptional regulatory proteins. Annu Rev Biochem 1989, 58: 799–839. 10.1146/annurev.bi.58.070189.004055
Ahmad S, Sarai A: Moment-based prediction of DNA-binding proteins. J Mol Biol 2004,341(1):65–71. 10.1016/j.jmb.2004.05.058
Ahmad S, Gromiha MM, Sarai A: Analysis and prediction of DNA-binding proteins and their binding residues based on composition, sequence and structural information. Bioinformatics 2004,20(4):477–486. 10.1093/bioinformatics/btg432
Tsuchiya Y, Kinoshita K, Nakamura H: Structure-based prediction of DNA-binding sites on proteins using the empirical preference of electrostatic potential and the shape of molecular surfaces. Proteins 2004,55(4):885–894. 10.1002/prot.20111
Bhardwaj N, Langlois RE, Zhao G, Lu H: Kernel-based machine learning protocol for predicting DNA-binding proteins. Nucleic Acids Res 2005,33(20):6486–6493. 10.1093/nar/gki949
Bhardwaj N, Lu H: Residue-level prediction of DNA-binding sites and its application on DNA-binding protein predictions. FEBS Lett 2007,581(5):1058–1066. 10.1016/j.febslet.2007.01.086
Yu X, Cao J, Cai Y, Shi T, Li Y: Predicting rRNA-, RNA-, and DNA-binding proteins from primary structure with support vector machines. J Theor Biol 2006,240(2):175–184. 10.1016/j.jtbi.2005.09.018
Szilagyi A, Skolnick J: Efficient prediction of nucleic acid binding function from low-resolution protein structures. J Mol Biol 2006,358(3):922–933. 10.1016/j.jmb.2006.02.053
Ahmad S, Sarai A: PSSM-based prediction of DNA binding sites in proteins. BMC Bioinformatics 2005, 6: 33. 10.1186/1471-2105-6-33
Kuznetsov IB, Gou Z, Li R, Hwang S: Using evolutionary and structural information to predict DNA-binding sites on DNA-binding proteins. Proteins 2006,64(1):19–27. 10.1002/prot.20977
Tjong H, Zhou HX: DISPLAR: an accurate method for predicting DNA-binding sites on protein surfaces. Nucleic Acids Res 2007,35(5):1465–1477. 10.1093/nar/gkm008
Luscombe NM, Thornton JM: Protein-DNA interactions: amino acid conservation and the effects of mutations on binding specificity. J Mol Biol 2002,320(5):991–1009. 10.1016/S0022-2836(02)00571-5
Stawiski EW, Gregoret LM, Mandel-Gutfreund Y: Annotating nucleic acid-binding function based on protein structure. J Mol Biol 2003,326(4):1065–1079. 10.1016/S0022-2836(03)00031-7
Henikoff S, Henikoff JG: Amino acid substitution matrices from protein blocks. Proc Natl Acad Sci USA 1992,89(22):10915–10919. 10.1073/pnas.89.22.10915
Hobohm U, Sander C: Enlarged representative set of protein structures. Protein Sci 1994,3(3):522–524.
Passner JM, Ryoo HD, Shen L, Mann RS, Aggarwal AK: Structure of a DNA-bound Ultrabithorax-Extradenticle homeodomain complex. Nature 1999,397(6721):714–719. 10.1038/17833
LaRonde-LeBlanc NA, Wolberger C: Structure of HoxA9 and Pbx1 bound to DNA: Hox hexapeptide and DNA recognition anterior to posterior. Genes Dev 2003,17(16):2060–2072. 10.1101/gad.1103303
Dutnall RN, Tafrov ST, Sternglanz R, Ramakrishnan V: Structure of the histone acetyltransferase Hat1: a paradigm for the GCN5-related N-acetyltransferase superfamily. Cell 1998,94(4):427–438. 10.1016/S0092-8674(00)81584-6
Williams DC Jr, Cai M, Clore GM: Molecular basis for synergistic transcriptional activation by Oct1 and Sox2 revealed from the solution structure of the 42-kDa Oct1.Sox2.Hoxb1-DNA ternary transcription factor complex. J Biol Chem 2004,279(2):1449–1457. 10.1074/jbc.M309790200
Konagurthu AS, Whisstock JC, Stuckey PJ, Lesk AM: MUSTANG: a multiple structural alignment algorithm. Proteins 2006,64(3):559–574. 10.1002/prot.20921
Murzin AG, Brenner SE, Hubbard T, Chothia C: SCOP: a structural classification of proteins database for the investigation of sequences and structures. J Mol Biol 1995,247(4):536–540. 10.1006/jmbi.1995.0159
Luscombe NM, Laskowski RA, Thornton JM: Amino acid-base interactions: a three-dimensional analysis of protein-DNA interactions at an atomic level. Nucleic Acids Res 2001,29(13):2860–2874. 10.1093/nar/29.13.2860
Morozov AV, Havranek JJ, Baker D, Siggia ED: Protein-DNA binding specificity predictions with structural models. Nucleic Acids Res 2005,33(18):5781–5798. 10.1093/nar/gki875
Holm L, Sander C: Protein structure comparison by alignment of distance matrices. J Mol Biol 1993,233(1):123–138. 10.1006/jmbi.1993.1489
Pearson WR, Lipman DJ: Improved tools for biological sequence comparison. Proc Natl Acad Sci USA 1988,85(8):2444–2448. 10.1073/pnas.85.8.2444
Pearson WR: Effective protein sequence comparison. Methods Enzymol 1996, 266: 227–258.
Pearson WR: Flexible sequence similarity searching with the FASTA3 program package. Methods Mol Biol 2000, 132: 185–219.
Smith TF: The art of matchmaking: sequence alignment methods and their structural implications. Structure 1999,7(1):R7-R12. 10.1016/S0969-2126(99)80003-3
Skolnick J, Fetrow JS: From genes to protein structure and function: novel applications of computational approaches in the genomic era. Trends Biotechnol 2000,18(1):34–39. 10.1016/S0167-7799(99)01398-0
Needleman SB, Wunsch CD: A general method applicable to the search for similarities in the amino acid sequence of two proteins. J Mol Biol 1970,48(3):443–453. 10.1016/0022-2836(70)90057-4
Smith TF, Waterman MS: Identification of common molecular subsequences. J Mol Biol 1981,147(1):195–197. 10.1016/0022-2836(81)90087-5
Karplus K, Barrett C, Hughey R: Hidden Markov models for detecting remote protein homologies. Bioinformatics 1998,14(10):846–856. 10.1093/bioinformatics/14.10.846
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990,215(3):403–410.
J.-M. Yang was supported by National Science Council and partial support of the ATU plan by MOE. Authors are grateful to both the hardware and software supports of the Structural Bioinformatics Core Facility at National Chiao Tung University.
This article has been published as part of BMC Bioinformatics Volume 9 Supplement 6, 2008: Symposium of Computations in Bioinformatics and Bioscience (SCBB07). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/9?issue=S6.
The authors declare that they have no competing interests.
Yao-Lin Chang, Huai-Kuang Tsai contributed equally to this work.
About this article
Cite this article
Chang, YL., Tsai, HK., Kao, CY. et al. Evolutionary conservation of DNA-contact residues in DNA-binding domains. BMC Bioinformatics 9 (Suppl 6), S3 (2008). https://doi.org/10.1186/1471-2105-9-S6-S3