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.     

Insights into protein-DNA interactions through structure network analysis

Protein-DNA interactions are crucial for many cellular processes. Now with the increased availability of structures of protein-DNA complexes, gaining deeper insights into the nature of protein-DNA interactions has become possible. Earlier, investigations have characterized the interface properties by considering pairwise interactions. However, the information communicated along the interfaces is rarely a pairwise phenomenon, and we feel that a global picture can be obtained by considering a protein-DNA complex as a network of noncovalently interacting systems. Furthermore, most of the earlier investigations have been carried out from the protein point of view (protein-centric), and the present network approach aims to combine both the protein-centric and the DNA-centric points of view. Part of the study involves the development of methodology to investigate protein-DNA graphs/networks with the development of key parameters. A network representation provides a holistic view of the interacting surface and has been reported here for the first time. The second part of the study involves the analyses of these graphs in terms of clusters of interacting residues and the identification of highly connected residues (hubs) along the protein-DNA interface. A predominance of deoxyribose-amino acid clusters in beta-sheet proteins, distinction of the interface clusters in helix-turn-helix, and the zipper-type proteins would not have been possible by conventional pairwise interaction analysis. Additionally, we propose a potential classification scheme for a set of protein-DNA complexes on the basis of the protein-DNA interface clusters. This provides a general idea of how the proteins interact with the different components of DNA in different complexes. Thus, we believe that the present graph-based method provides a deeper insight into the analysis of the protein-DNA recognition mechanisms by throwing more light on the nature and the specificity of these interactions.     

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.     

Binding of the estradiol-receptor complex to reconstituted nucleoacidic protein from calf uterus

Non-histone protein-DNA complexes with acceptor activity for estradiol-receptor complexes were reconstituted from fractionated calf uterine chromatin. Acceptor activity had tissue specificity with target tissue binding exceeding non-target tissue binding. The binding of estradiol-receptor complexes to acceptor sites was dependent on intact non-histone protein-DNA complexes, reconstituted select non-histone proteins, and protein equivalent: DNA reconstitution ratios. [³H]Estradiol-receptor complexes were bound to reconstituted non-histone protein-DNA complexes (i.e., nucleoacidic protein) with a high affinity and with a limited number of binding sites. Fractionation of uterine chromatin non-histone proteins identified two major sets of non-histone proteins which had acceptor activity when reconstituted with DNA. Thus, it seems possible to reconstitute nucleoacidic protein fractions with specific acceptor activity for the calf uterine estrogen receptor.     

An overview of the structures of protein-DNA complexes

On the basis of a structural analysis of 240 protein-DNA complexes contained in the Protein Data Bank (PDB), we have classified the DNA-binding proteins involved into eight different structural/functional groups, which are further classified into 54 structural families. Here we present this classification and review the functions, structures and binding interactions of these protein-DNA complexes.     

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.     

Design of a Fluorescent Electrophoretic Mobility Shift Assay Improved for the Quantitative and Multiple Analysis of Protein-DNA Complexes

We describe a protocol for the fluorescent electrophoretic mobility shift assay improved for the quantitative analysis of protein-DNA complexes. Fluorescent-labeled oligonucleotide probes incubated with nuclear proteins were followed by electrophoresis. The signals for protein-DNA complexes were measured and normalized with fluorescent-labeled marker using fragment analysis software. This assay proved reliable measurement and multiple detection of DNA binding proteins.     

Design of a fluorescent electrophoretic mobility shift assay improved for the quantitative and multiple analysis of protein-DNA complexes

We describe a protocol for the fluorescent electrophoretic mobility shift assay improved for the quantitative analysis of protein-DNA complexes. Fluorescent-labeled oligonucleotide probes incubated with nuclear proteins were followed by electrophoresis. The signals for protein-DNA complexes were measured and normalized with fluorescent-labeled marker using fragment analysis software. This assay proved reliable measurement and multiple detection of DNA binding proteins.     

Intermolecular and intramolecular readout mechanisms in protein-DNA recognition

Protein-DNA recognition plays an essential role in the regulation of gene expression. Regulatory proteins are known to recognize specific DNA sequences directly through atomic contacts (intermolecular readout) and/or indirectly through the conformational properties of the DNA (intramolecular readout). However, little is known about the respective contributions made by these so-called direct and indirect readout mechanisms. We addressed this question by making use of information extracted from a structural database containing many protein-DNA complexes. We quantified the specificity of intermolecular (direct) readout by statistical analysis of base-amino acid interactions within protein-DNA complexes. The specificity of the intramolecular (indirect) readout due to DNA was quantified by statistical analysis of the sequence-dependent DNA conformation. Systematic comparison of these specificities in a large number of protein-DNA complexes revealed that both intermolecular and intramolecular readouts contribute to the specificity of protein-DNA recognition, and that their relative contributions vary depending upon the protein-DNA complexes. We demonstrated that combination of the intermolecular and intramolecular energies derived from the statistical analyses lead to enhanced specificity, and that the combined energy could explain experimental data on binding affinity changes caused by base mutations. These results provided new insight into the relationship between specificity and structure in the process of protein-DNA recognition, which would lead to prediction of specific protein-DNA binding sites.     

Geometric analysis and comparison of protein-DNA interfaces: why is there no simple code for recognition?

Structural studies of protein-DNA complexes have shown that there are many distinct families of DNA-binding proteins, and have shown that there is no simple "code" describing side-chain/base interactions. However, systematic analysis and comparison of protein-DNA complexes has been complicated by the diversity of observed contacts, the sheer number of complexes currently available and the absence of any consistent method of comparison that retains detailed structural information about the protein-DNA interface. To address these problems, we have developed geometric methods for characterizing the local structural environment in which particular side-chain/base interactions are observed. In particular, we develop methods for analyzing and comparing spatial relationships at the protein-DNA interface. Our method involves attaching local coordinate systems to the DNA bases and to the C(alpha) atoms of the peptide backbone (these are relatively rigid structural units). We use these tools to consider how the position and orientation of the polypeptide backbone (with respect to the DNA) helps to determine what contacts are possible at any given position in a protein-DNA complex. Here, we focus on base contacts that are made in the major groove, and we use spatial relationships in analyzing: (i) the observed patterns of side-chain/base interactions; (ii) observed helix docking orientations; (iii) family/subfamily relationships among DNA-binding proteins; and (iv) broader questions about evolution, altered specificity mutants and the limits for the design of new DNA-binding proteins. Our analysis, which highlights differences in spatial relationships in different complexes and at different positions in a complex, helps explain why there is no simple, general code for protein-DNA recognition.     

Structure model of a complex between the factor for inversion stimulation (FIS) and DNA: Modeling protein-DNA complexes with dyad symmetry and known protein structures

A method is presented to predict overall conformations of protein-DNA complexes on the basis of the known three-dimensional structures of the proteins. The method is restricted to proteins with a common twofold symmetry axis, which show only minor conformational changes upon binding to DNA. The method uses a numerical finite difference solution of the linearized Poisson-Boltzmann equation and subsequent energy minimization cycles. Structural parameters—the rotation angle of the DNA relative to the protein around the common symmetry axis, the protein-DNA distance, and intermolecular hydrogen-bonding contacts—are presented for two test cases, DNA bound to CAP (catabolite gene activator protein) and to the Cro-repressor of bacteriophage 434. The DNA curvature in the starting model of the docking procedure was chosen as a smoothed approximation of the conformation found in the X-ray structures of these complexes. The method is further used to predict the unknown structure of the complex between the factor for inversion stimulation (FIS) and DNA, which is bent upon binding to FIS. In contrast to the test cases, the unknown curvature of the starting model is derived from a calibration of electrostatic precalculations for different proteins according to crystallographically observed DNA bending. The results of the modeling are in good accordance with the experimentally observed overall structure of protein-DNA complexes for the two test cases; for FIS, they correspond to several of the experimentally proposed protein-DNA contacts. © 1996 Wiley-Liss, Inc.     

Application of the NZ‐1 Fab as a crystallization chaperone for PA tag‐inserted target proteins

An antibody fragment that recognizes the tertiary structure of a target protein with high affinity can be utilized as a crystallization chaperone. Difficulties in establishing conformation‐specific antibodies, however, limit the applicability of antibody fragment‐assisted crystallization. Here, we attempted to establish an alternative method to promote the crystallization of target proteins using an already established anti‐tag antibody. The monoclonal antibody NZ‐1 recognizes the PA tag with an extremely high affinity. It was also established that the PA tag is accommodated in the antigen‐binding pocket in a bent conformation, compatible with an insertion into loop regions on the target. We, therefore, explored the application of NZ‐1 Fab as a crystallization chaperone that complexes with a target protein displaying a PA tag. Specifically, we inserted the PA tag into the β‐hairpins of the PDZ tandem fragment of a bacterial Site‐2 protease. We crystallized the PA‐inserted PDZ tandem mutants with the NZ‐1 Fab and solved the co‐crystal structure to analyze their interaction modes. Although the initial insertion designs produced only moderate‐resolution structures, eliminating the solvent‐accessible space between the NZ‐1 Fab and target PDZ tandem improved the diffraction qualities remarkably. Our results demonstrate that the NZ‐1‐PA system efficiently promotes crystallization of the target protein. The present work also suggests that β‐hairpins are suitable sites for the PA insertion because the PA tag contains a Pro‐Gly sequence with a propensity for a β‐turn conformation.     

Characterization of the control catabolite protein of gluconeogenic genes repressor by fluorescence cross-correlation spectroscopy and other biophysical approaches

Determination of the physical parameters underlying protein-DNA interactions is crucial for understanding the regulation of gene expression. In particular, knowledge of the stoichiometry of the complexes is a prerequisite to determining their energetics and functional molecular mechanisms. However, the experimental determination of protein-DNA complex stoichiometries remains challenging. We used fluorescence cross-correlation spectroscopy (FCCS) to investigate the interactions of the control catabolite protein of gluconeogenic genes, a key metabolic regulator in Gram-positive bacteria, with two oligonucleotides derived from its target operator sequences, gapB and pckA. According to our FCCS experiments, the stoichiometry of binding is twofold larger for the pckA target than for gapB. Correcting the FCCS data for protein self-association indicated that control catabolite protein of gluconeogenic genes forms dimeric complexes on the gapB target and tetrameric complexes on the pckA target. Analytical ultracentrifugation coupled with fluorescence anisotropy and hydrodynamic modeling allowed unambiguous confirmation of this result. The use of multiple complementary techniques to characterize these complexes should be employed wherever possible. However, there are cases in which analytical ultracentrifugation is precluded, due to protein stability, solubility, or availability, or, more obviously, when the studies are carried out in live cells. If information concerning the self-association of the protein is available, FCCS can be used for the direct and simultaneous determination of the affinity, cooperativity, and stoichiometry of protein-DNA complexes in a concentration range and conditions relevant to the regulation of these interactions.     

Fine structure and activity of discrete RAG-HMG complexes on V(D)J recombination signals

Two lymphoid cell-specific proteins, RAG-1 and RAG-2, initiate V(D)J recombination by introducing DNA breaks at recombination signal sequences (RSSs). Although the RAG proteins themselves bind and cleave DNA substrates containing either a 12-RSS or a 23-RSS, DNA-bending proteins HMG-1 and HMG-2 are known to promote these processes, particularly with 23-RSS substrates. Using in-gel cleavage assays and DNA footprinting techniques, I analyzed the catalytic activity and protein-DNA contacts in discrete 12-RSS and 23-RSS complexes containing the RAG proteins and either HMG-1 or HMG-2. I found that both the cleavage activity and the pattern of protein-DNA contacts in RAG-HMG complexes assembled on 12-RSS substrates closely resembled those obtained from analogous 12-RSS complexes lacking HMG protein. In contrast, 23-RSS complexes containing both RAG proteins and either HMG-1 or HMG-2 exhibited enhanced cleavage activity and displayed an altered distribution of cleavage products compared to 23-RSS complexes containing only RAG-1 and RAG-2. Moreover, HMG-dependent heptamer contacts in 23-RSS complexes were observed. The protein-DNA contacts in RAG-RSS-HMG complexes assembled on 12-RSS or 23-RSS substrates were strikingly similar at comparable positions, suggesting that the RAG proteins mediate HMG-dependent heptamer contacts in 23-RSS complexes. Results of ethylation interference experiments suggest that the HMG protein is positioned 5' of the nonamer in 23-RSS complexes, interacting largely with the side of the duplex opposite the one contacting the RAG proteins. Thus, HMG protein plays the dual role of bringing critical elements of the 23-RSS heptamer into the same phase as the 12-RSS to promote RAG binding and assisting in the catalysis of 23-RSS cleavage.     

A new analytical scale DNA affinity binding assay for analyses of specific protein-DNA interactions.

We describe a rapid analytical assay for identification of proteins binding to specific DNA sequences. The DAPSTER assay (DNA affinity preincubation specificity test of recognition assay) is a DNA affinity chromatography-based microassay that can discriminate between specific and nonspecific protein-DNA interactions. The assay is sensitive and can detect protein-DNA interactions and larger multicomponent complexes that can be missed by other analytical methods. Here we describe in detail the optimization and utilization of the DAPSTER assay to isolate AP-1 complexes and associated proteins in multimeric complexes bound to the AP-1 DNA element.     

Uracil interference, a rapid and general method for defining protein-DNA interactions involving the 5-methyl group of thymines: the GCN4-DNA complex.

We describe a novel uracil interference method for examining protein contacts with the 5-methyl group of thymines. The protein of interest is incubated with target DNA containing randomly distributed deoxyuracil substitutions that is generated by carrying out the polymerase chain reaction in the presence of a mixture of TTP and dUTP. After separating DNA-protein complexes away from unbound DNA, the locations of deoxyuracil residues that either do or do not interfere with DNA-binding are determined by cleavage with uracil-N-glycosylase followed by piperidine. Using this uracil interference assay, we show that the methyl groups of the four core thymines, but not the two peripheral thymines, of the optimal binding site (ATG-ACTCAT) are important for high affinity binding of GCN4. Similar, but not identical, results are obtained using KMnO4 interference, another method used for studying protein-DNA interactions involving thymine residues. These observations strongly suggest that GCN4 directly contacts the 5-methyl groups of the four core thymines that lie in the major groove of the target DNA. Besides providing specific structural information about protein-DNA complexes, uracil interference should also be useful for identifying DNA-binding proteins and their target sites in eukaryotic promoter regions.