Sequence Recognition of DNA by Protein-Induced Conformational Transitions

The binding of proteins to specific sequences of DNA is an important feature of virtually all DNA transactions. Proteins recognize specific DNA sequences using both direct readout (sensing types and positions of DNA functional groups) and indirect readout (sensing DNA conformation and deformability). Previously we showed that the P22 c2 repressor N-terminal domain (P22R NTD) forces the central non-contacted 5'-ATAT-3' sequence of the DNA operator into the B′ state, a state known to affect DNA hydration, rigidity and bending. Usually the B′ state, with a narrow minor groove and a spine of hydration, is reserved for A-tract DNA (TpA steps disrupt A-tracts). Here, we have co-crystallized P22R NTD with an operator containing a central 5′-ACGT-3′ sequence in the non-contacted region. C·G base pairs have not previously been observed in the B′ state and are thought to prevent it. However, P22R NTD induces a narrow minor groove and a spine of hydration to 5'-ACGT-3'. We observe that C·G base pairs have distinctive destabilizing and disordering effects on the spine of hydration. It appears that the reduced stability of the spine results in a higher energy cost for the B to B′ transition. The differential effect of DNA sequence on the barrier to this transition allows the protein to sense the non-contacted DNA sequence.     

Monovalent cations regulate DNA sequence recognition by 434 repressor

The bacteriophage 434 repressor distinguishes between its six naturally occurring binding sites using indirect readout. In indirect readout, sequence-dependent differences in the structure and flexibility of non-contacted bases in a protein's DNA-binding site modulate the affinity of DNA for protein. The conformation and flexibility of a DNA sequence can be influenced by the interaction of the DNA bases or backbone with solution components. We examined the effect of changing the cation-type present in solution on the stability and structure of 434 repressor complexes with wild-type and mutant OR1 and OR3, binding sites that differ in their contacted and non-contacted base sequences. We find that the affinity of repressor for OR1, but not for OR3, depends remarkably on the type and concentration of monovalent cation. Moreover, the formation of a stable, specific repressor-OR1 complex requires the presence of monovalent cations; however, repressor-OR3 complex formation has no such requirement. Changing monovalent cation type alters the ability of repressor to protect OR1, but not OR3, from *OH radical cleavage. Altering the relative length of the poly(dA) x poly(dT) tract in the non-contacted regions of the OR1 and OR3 can reverse the cation sensitivity of repressor's affinities for these two sites. Taken together these findings show that cation-dependent alterations in DNA structure underlies indirect readout of DNA sequence by 434 repressor and perhaps other proteins.     

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.     

Recognition of DNA structure by 434 repressor.

In complexes of bacteriophage 434 binding sites with 434 repressor the central 4 bp of the 14 bp site are not contacted by the protein, although changes in these bases alter binding site affinity for the repressor. Our previous data suggested that the ability of the non-contacted central bases to be overtwisted in repressor-DNA complexes governs affinity of the binding site for 434 repressor. This idea was tested by examining the affinity of two central sequence variant 434 binding sites for 434 repressor as a function of binding site average twist. The 434 repressor preferred the relatively overwound binding site to the two more underwound forms. The greatest affinity enhancement resulting from increasing twist was observed with a binding site that is relatively underwound and more resistant to twisting deformation. Consistent with the idea that 434 repressor overtwists its binding site upon DNA binding, we show that 434 repressor is capable of binding to sites bearing a single base insertion in their center (a 15mer), but binds poorly to binding sites bearing central base deletions (12mer and 13mer). The N-terminal dimer interface plays a large role in determining 434 repressor central base preferences. Mutations in this interface eliminate central base discrimination and/or site size preferences. These mutations also lead to changes in the size of the repressor footprint on the various sized DNA sites that are consistent with their binding characteristics.     

P22 c2 repressor-operator complex: mechanisms of direct and indirect readout

The P22 c2 repressor protein (P22R) binds to DNA sequence-specifically and helps to direct the temperate lambdoid bacteriophage P22 to the lysogenic developmental pathway. We describe the 1.6 A X-ray structure of the N-terminal domain (NTD) of P22R in a complex with a DNA fragment containing the synthetic operator sequence [d(ATTTAAGATATCTTAAAT)]2. This operator has an A-T base pair at position 9L and a T-A base pair at position 9R and is termed DNA9T. Direct readout: nondirectional van der Waals interactions between protein and DNA appear to confer sequence-specificity. The structure of the P22R NTD-DNA9T complex suggests that sequence-specificity arises substantially from lock-and-key interaction of a valine with a complementary binding cleft on the major groove surface of DNA9T. The cleft is formed by four methyl groups on sequential base pairs of 5'-TTAA-3'. The valine cleft is intrinsic to the DNA sequence and does not arise from protein-induced DNA conformational changes. Protein-DNA hydrogen bonding plays a secondary role in specificity. Indirect readout: it is known that the noncontacted bases in the center of the complex are important determinants of affinity. The protein induces a transition of the noncontacted region from B-DNA to B'-DNA. The B' state is characterized by a narrow minor groove and a zigzag spine of hydration. The free energy of the transition from B- to B'-DNA is known to depend on the sequence. Thus, the observed DNA conformation and hydration allows for the formulation of a predictive model of the indirect readout phenomenon.     

Specific minor groove solvation is a crucial determinant of DNA binding site recognition

The DNA sequence preferences of nearly all sequence specific DNA binding proteins are influenced by the identities of bases that are not directly contacted by protein. Discrimination between non-contacted base sequences is commonly based on the differential abilities of DNA sequences to allow narrowing of the DNA minor groove. However, the factors that govern the propensity of minor groove narrowing are not completely understood. Here we show that the differential abilities of various DNA sequences to support formation of a highly ordered and stable minor groove solvation network are a key determinant of non-contacted base recognition by a sequence-specific binding protein. In addition, disrupting the solvent network in the non-contacted region of the binding site alters the protein''s ability to recognize contacted base sequences at positions 5–6 bases away. This observation suggests that DNA solvent interactions link contacted and non-contacted base recognition by the protein.     

Solution measurement of DNA curvature in papillomavirus E2 binding sites

‘Indirect readout’ refers to the proposal that proteins can recognize the intrinsic three-dimensional shape or flexibility of a DNA binding sequence apart from direct protein contact with DNA base pairs. The differing affinities of human papillomavirus (HPV) E2 proteins for different E2 binding sites have been proposed to reflect indirect readout. DNA bending has been observed in X-ray structures of E2 protein–DNA complexes. X-ray structures of three different E2 DNA binding sites revealed differences in intrinsic curvature. DNA sites with intrinsic curvature in the direction of protein-induced bending were bound more tightly by E2 proteins, supporting the indirect readout model. We now report solution measurements of intrinsic DNA curvature for three E2 binding sites using a sensitive electrophoretic phasing assay. Measured E2 site curvature agrees well the predictions of a dinucleotide model and supports an indirect readout hypothesis for DNA recognition by HPV E2.     

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.     

Recognition of DNA by single-chain derivatives of the phage 434 repressor: high affinity binding depends on both the contacted and non-contacted base pairs.

Single-chain derivatives of the phage 434 repressor, termed single-chain repressors, contain covalently dimerized DNA-binding domains (DBD) which are connected with a peptide linker in a head-to-tail arrangement. The prototype RR69 contains two wild-type DBDs, while RR*69 contains a wild-type and an engineered DBD. In this latter domain, the DNA- contacting amino acids of thealpha3 helix of the 434 repressor are replaced by the corresponding residues of the related P22 repressor. We have used binding site selection, targeted mutagenesis and binding affinity studies to define the optimum DNA recognition sequence for these single-chain proteins. It is shown that RR69 recognizes DNA sequences containing the consensus boxes of the 434 operators in a palindromic arrangement, and that RR*69 optimally binds to non-palindromic sequences containing a 434 operator box and a TTAA box of which the latter is present in most P22 operators. The spacing of these boxes, as in the 434 operators, is 6 bp. The DNA-binding of both single-chain repressors, similar to that of the 434 repressor, is influenced indirectly by the sequence of the non-contacted, spacer region. Thus, high affinity binding is dependent on both direct and indirect recognition. Nonetheless, the single-chain framework can accommodate certain substitutions to obtain altered DNA-binding specificity and RR*69 represents an example for the combination of altered direct and unchanged indirect readout mechanisms.     

Sequence-dependent DNA deformability studied using molecular dynamics simulations

Proteins recognize specific DNA sequences not only through direct contact between amino acids and bases, but also indirectly based on the sequence-dependent conformation and deformability of the DNA (indirect readout). We used molecular dynamics simulations to analyze the sequence-dependent DNA conformations of all 136 possible tetrameric sequences sandwiched between CGCG sequences. The deformability of dimeric steps obtained by the simulations is consistent with that by the crystal structures. The simulation results further showed that the conformation and deformability of the tetramers can highly depend on the flanking base pairs. The conformations of xATx tetramers show the most rigidity and are not affected by the flanking base pairs and the xYRx show by contrast the greatest flexibility and change their conformations depending on the base pairs at both ends, suggesting tetramers with the same central dimer can show different deformabilities. These results suggest that analysis of dimeric steps alone may overlook some conformational features of DNA and provide insight into the mechanism of indirect readout during protein–DNA recognition. Moreover, the sequence dependence of DNA conformation and deformability may be used to estimate the contribution of indirect readout to the specificity of protein–DNA recognition as well as nucleosome positioning and large-scale behavior of nucleic acids.     

A code governing specific binding of regulatory proteins to DNA and structure of stereospecific sites of regulatory proteins]

A model is proposed for the structure of stereospecific sites in regulatory proteins. On its basis a possible code is suggested that governs the binding of regulatory proteins at specific control sites on DNA. Stereospecific sites of regulatory proteins are assumed to contain pairs of antiparallel polypeptide chain segments which form a right-hand twisted antiparallel beta-sheet, with single-stranded regions at the ends of the beta-structure. The model predicts that binding reaction between a regulatory protein and double-helical DNA is a cooperative phenomenon and is accompanied by significant structural alteration at the stereospecific site of the protein. Half of hydrogen bonds normally existing in beta-structure are broken upon complex formation with DNA and a new set of hydrogen bonds is formed between polypeptide amide groups and DNA base pairs. In a stereospecific site, one chain (t-chain) is attached through hydrogen bonds to the carbonyl oxygens of pyramides and N3 adenines lying in one DNA strand, while the second polypeptide chain (g chain) is hydrogen bonded to the 2-amino groups of guanine residues lying in the opposite DNA strand. The amide groups serve as specific reaction sites being hydrogen bond acceptors in g-chain and hydrogen bond donors in t-chain. The single-stranded portions of t- and g-chains lying in neighbouring subunits of regulatory protein interact with each other forming deformed beta-sheets. The recognition of regulatory sequences by proteins is based on the structural complementarity between stereospecific sites of regulatory proteins and base pairs sequences at the control sites. An essential feature of these sequences is the asymmetrical distribution of guanine residues between the two DNA strands. The code predicts that there are six fundamental amino acid residues (serine, threonine, asparagine, histidine, glutamine and cysteine) whose sequence in stereospecific site determines the base pair sequence to which a given regulatory protein would bind preferentially. The code states a correspondence between four amino acid residues at the stereospecific site of regulatory protein with the two residues being in t- and g-segments, respectively, and AT(GC) base pair at the control site. It is thus possible to determine which amino acid residues in the repressor and which base pairs in the operator DNA are involved in specific interactions with each other, as exemplified by lac repressor binding to lac operator.     

Determinants of DNA mismatch recognition within the polymerase domain of the Klenow fragment.

The Klenow fragment of Escherichia coli DNA polymerase I catalyzes template-directed synthesis of DNA and uses a separate 3'-5' exonuclease activity to edit misincorporated bases. The polymerase and exonuclease activities are contained in separate structural domains. In this study, nine Klenow fragment derivatives containing mutations within the polymerase domain were examined for their interaction with model primer-template duplexes. The partitioning of the DNA primer terminus between the polymerase and 3'-5' exonuclease active sites of the mutant proteins was assessed by time-resolved fluorescence anisotropy, utilizing a dansyl fluorophore attached to the DNA. Mutation of N845 or R668 disrupted favorable interactions between the Klenow fragment and a duplex containing a matched terminal base pair but had little effect when the terminus was mismatched. Thus, N845 and R668 are required for recognition of correct terminal base pairs in the DNA substrate. Mutation of N675, R835, R836, or R841 resulted in tighter polymerase site binding of DNA, suggesting that the side chains of these residues induce strain in the DNA and/or protein backbone. A double mutant (N675A/R841A) showed an even greater polymerase site partitioning than was displayed by either single mutation, indicating that such strain is additive. In both groups of mutant proteins, the ability to discriminate between duplexes containing matched or mismatched base pairs was impaired. In contrast, mutation of K758 or Q849 had no effect on partitioning relative to wild type, regardless of DNA mismatch character. These results demonstrate that DNA mismatch recognition is dependent on specific amino acid residues within the polymerase domain and is not governed solely by thermodynamic differences between correct and mismatched base pairs. Moreover, this study suggests a mechanism whereby the Klenow fragment is able to recognize polymerase errors following a misincorporation event, leading to their eventual removal by the 3'-5' exonuclease activity.     

The structure of SgrAI bound to DNA; recognition of an 8 base pair target

The three-dimensional X-ray crystal structure of the ‘rare cutting’ type II restriction endonuclease SgrAI bound to cognate DNA is presented. SgrAI forms a dimer bound to one duplex of DNA. Two Ca²⁺ bind in the enzyme active site, with one ion at the interface between the protein and DNA, and the second bound distal from the DNA. These sites are differentially occupied by Mn²⁺, with strong binding at the protein–DNA interface, but only partial occupancy of the distal site. The DNA remains uncleaved in the structures from crystals grown in the presence of either divalent cation. The structure of the dimer of SgrAI is similar to those of Cfr10I, Bse634I and NgoMIV, however no tetrameric structure of SgrAI is observed. DNA contacts to the central CCGG base pairs of the SgrAI canonical target sequence (CR|CCGGYG, | marks the site of cleavage) are found to be very similar to those in the NgoMIV/DNA structure (target sequence G|CCGGC). Specificity at the degenerate YR base pairs of the SgrAI sequence may occur via indirect readout using DNA distortion. Recognition of the outer GC base pairs occurs through a single contact to the G from an arginine side chain located in a region unique to SgrAI.     

Sequence-specific interaction of the Salmonella Hin recombinase in both major and minor grooves of DNA.

The Hin recombinase of Salmonella catalyzes a site-specific recombination event which leads to flagellar phase variation. Starting with a fully symmetrical recombination site, hixC, a set of 40 recombination sites which vary by pairs of single base substitutions was constructed. This set was incorporated into the Salmonella-specific bacteriophage P22 based challenge phage selection and used to define the DNA sequence determinants for the binding of Hin to DNA in vivo. The critical sequence-specific contacts between a Hin monomer and a 13 bp hix half-site are at two T:A base pairs in the major groove of the DNA which are separated by one base pair, and two consecutive A:T contacts in the minor groove. The base substitutions in the major groove recognition portion which were defective in binding Hin still retained residual binding capability in vivo, while the base pair substitutions affecting the minor groove recognition region lost all in vivo binding. Using in vitro binding assays, Hin was found to bind to hix symmetrical sites with A:T base pairs or I:C base pairs in the minor groove recognition sequences, but not to G:C base pairs. In separate in vitro binding assays, Hin was equally defective in binding to either a G:C or a I:C contact in a major groove recognition sequence. Results from in vitro binding assays to hix sites in which 3-deazaadenine was substituted for adenine are consistent with Hin making a specific contact to either the N3 of adenine or O2 of thymine in the minor groove within the hix recombination site on each symmetric half-site. These results taken with the results of previous studies on the DNA binding domain of Hin suggest a sequence-specific minor groove DNA binding motif.     

EcoRV restriction endonuclease binds all DNA sequences with equal affinity

In the presence of MgCl2, the EcoRV restriction endonuclease cleaves its recognition sequence on DNA at least a million times more readily than any other sequence. In this study, the binding of the EcoRV restriction enzyme to DNA was examined in the absence of Mg2+. With each DNA fragment tested, several DNA-protein complexes were detected by electrophoresis through polyacrylamide. No differences were observed between isogenic DNA molecules that either contained or lacked the EcoRV recognition site. The number of complexes with each fragment varied with the length of the DNA. Three complexes were formed with a DNA molecule of 55 base pairs, corresponding to the DNA bound to 1, 2, or 3 molecules of the protein, while greater than 15 complexes were formed with a DNA of 381 base pairs. A new method was developed to analyze the binding of a protein to multiple sites on DNA. The method showed that the EcoRV enzyme binds to all DNA sequences, including the EcoRV recognition site, with the same equilibrium constant, though two molecules of the protein bind preferentially to adjacent sites on the DNA in a cooperative fashion. All of the complexes with a substrate that contained the EcoRV site dissociated upon addition of competitor DNA, but when the competitor was mixed with MgCl2, a fraction of the substrate was cleaved at the EcoRV site. The fraction cleaved was due mainly to the translocation of the enzyme from nonspecific sites on the DNA to the specific site.     

Human alkyladenine DNA glycosylase employs a processive search for DNA damage

DNA repair proteins conduct a genome-wide search to detect and repair sites of DNA damage wherever they occur. Human alkyladenine DNA glycosylase (AAG) is responsible for recognizing a variety of base lesions, including alkylated and deaminated purines, and initiating their repair via the base excision repair pathway. We have investigated the mechanism by which AAG locates sites of damage using an oligonucleotide substrate containing two sites of DNA damage. This substrate was designed so that AAG randomly binds to either of the two lesions. AAG-catalyzed base excision creates a repair intermediate, and the subsequent partitioning between dissociation and diffusion to the second site can be quantified from the rates of formation of the different products. Our results demonstrate that AAG has the ability to slide for short distances along DNA at physiological salt concentrations. The processivity of AAG decreases with increasing ionic strength to become fully distributive at high ionic strengths, suggesting that electrostatic interactions between the negatively charged DNA and the positively charged DNA binding surface are important for nonspecific DNA binding. Although the amino terminus of the protein is dispensable for glycosylase activity at a single site, we find that deletion of the 80 amino-terminal amino acids significantly decreases the processivity of AAG. These observations support the idea that diffusion on undamaged DNA contributes to the search for sites of DNA damage.     

Distamycin and penta-N-methylpyrrolecarboxamide binding sites on native DNA. A comparison of methidiumpropyl-EDTA-Fe(II) footprinting and DNA affinity cleaving

Using two direct methods we have studied the binding locations and site sizes of distamycin and penta-N-methylpyrrolecarboxamide on three DNA restriction fragments from pBR322 plasmid. We find that methidiumpropyl-EDTA.Fe(II) footprinting and DNA affinity cleaving methods report common binding locations and site sizes for the tri- and pentapeptides bound to heterogeneous DNA. The tripeptide distamycin binds 5-base-pair sites with a preference for poly(dA).poly(dT) regions. The pentapeptide binds 6-7-base-pair sites with a preference for poly(dA).poly(dT) regions. These results are consistent with distamycin binding as an isogeometric helix to the minor groove of DNA with the four carboxamide N-H's hydrogen bonding five A + T base pairs. The data supports a model where each of the carboxamide N-H's can hydrogen bond to two bases, either O(2) of thymine or N(3) of adenine, located on adjacent base pairs on opposite strands of the helix. In most (but not all) cases the tri- and pentapeptide can adopt two orientations at each A + T rich binding site.