Solution parameters modulating DNA binding specificity of the restriction endonuclease EcoRV

The DNA binding stringency of restriction endonucleases is crucial for their proper function. The X-ray structures of the specific and non-cognate complexes of the restriction nuclease EcoRV are considerably different suggesting significant differences in the hydration and binding free energies. Nonetheless, the majority of studies performed at pH 7.5, optimal for enzymatic activity, have found a < 10-fold difference between EcoRV binding constants to the specific and nonspecific sequences in the absence of divalent ions. We used a recently developed self-cleavage assay to measure EcoRV-DNA competitive binding and to evaluate the influence of water activity, pH and salt concentration on the binding stringency of the enzyme in the absence of divalent ions. We find the enzyme can readily distinguish specific and nonspecific sequences. The relative specific-nonspecific binding constant increases strongly with increasing neutral solute concentration and with decreasing pH. The difference in number of associated waters between specific and nonspecific DNA-EcoRV complexes is consistent with the differences in the crystal structures. Despite the large pH dependence of the sequence specificity, the osmotic pressure dependence indicates little change in structure with pH. The large osmotic pressure dependence means that measurement of protein-DNA specificity in dilute solution cannot be directly applied to binding in the crowded environment of the cell. In addition to divalent ions, water activity and pH are key parameters that strongly modulate binding specificity of EcoRV.     

Non-cognate enzyme-DNA complex: structural and kinetic analysis of EcoRV endonuclease bound to the EcoRI recognition site GAATTC

The crystal structure of EcoRV endonuclease bound to non-cognate DNA at 2.0 angstroms resolution shows that very small structural adaptations are sufficient to ensure the extreme sequence specificity characteristic of restriction enzymes. EcoRV bends its specific GATATC site sharply by 50 degrees into the major groove at the center TA step, generating unusual base-base interactions along each individual DNA strand. In the symmetric non-cognate complex bound to GAATTC, the center step bend is relaxed to avoid steric hindrance caused by the different placement of the exocyclic thymine methyl groups. The decreased base-pair unstacking in turn leads to small conformational rearrangements in the sugar-phosphate backbone, sufficient to destabilize binding of crucial divalent metal ions in the active site. A second crystal structure of EcoRV bound to the base-analog GAAUTC site shows that the 50 degrees center-step bend of the DNA is restored. However, while divalent metals bind at high occupancy in this structure, one metal ion shifts away from binding at the scissile DNA phosphate to a position near the 3'-adjacent phosphate group. This may explain why the 10(4)-fold attenuated cleavage efficiency toward GAATTC is reconstituted by less than tenfold toward GAAUTC. Examination of DNA binding and bending by equilibrium and stopped-flow florescence quenching and fluorescence resonance energy transfer (FRET) methods demonstrates that the capacity of EcoRV to bend the GAATTC non-cognate site is severely limited, but that full bending of GAAUTC is achieved at only a threefold reduced rate compared with the cognate complex. Together, the structural and biochemical data demonstrate the existence of distinct mechanisms for ensuring specificity at the bending and catalytic steps, respectively. The limited conformational rearrangements observed in the EcoRV non-cognate complex provide a sharp contrast to the extensive structural changes found in a non-cognate BamHI-DNA crystal structure, thus demonstrating a diversity of mechanisms by which restriction enzymes are able to achieve specificity.     

The crystal structure of EcoRV endonuclease and of its complexes with cognate and non-cognate DNA fragments.

The crystal structure of EcoRV endonuclease has been determined at 2.5 A resolution and that of its complexes with the cognate DNA decamer GGGATATCCC (recognition sequence underlined) and the non-cognate DNA octamer CGAGCTCG at 3.0 A resolution. Two octamer duplexes of the non-cognate DNA, stacked end-to-end, are bound to the dimeric enzyme in B-DNA-like conformations. The protein--DNA interactions of this complex are prototypic for non-specific DNA binding. In contrast, only one cognate decamer duplex is bound and deviates considerably from canonical B-form DNA. Most notably, a kink of approximately 50 degrees is observed at the central TA step with a concomitant compression of the major groove. Base-specific hydrogen bonds between the enzyme and the recognition base pairs occur exclusively in the major groove. These interactions appear highly co-operative as they are all made through one short surface loop comprising residues 182-186. Numerous contacts with the sugar phosphate backbone extending beyond the recognition sequence are observed in both types of complex. However, the total surface area buried on complex formation is > 1800 A2 larger in the case of cognate DNA binding. Two acidic side chains, Asp74 and Asp90, are close to the reactive phosphodiester group in the cognate complex and most probably provide oxygen ligands for binding the essential cofactor Mg2+. An important role is also indicated for Lys92, which together with the two acidic functions appears to be conserved in the otherwise unrelated structure of EcoRI endonuclease. The structural results give new insight into the physical basis of the remarkable sequence specificity of this enzyme.     

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.     

Specificity from the synapsis of DNA elements by the Sfi I endonuclease.

The synapsis of DNA sites is a prerequisite for the reactions of many proteins that act at specific DNA sequences. The requirement for synapsis was investigated by analysing the reactions of Sfi I, a tetrameric restriction enzyme that cleaves DNA only after interacting with two recognition sites. In the presence of Mg2+, oligonucleotide duplexes with the cognate recognition sequence were cleaved rapidly, with cooperative kinetics, while non-cognate duplexes were not cleaved. In the absence of Mg2+, the primary complex formed by Sfi I with cognate DNA contained two duplexes synapsed by the tetramer: a secondary complex containing one duplex was seen only at elevated Sfi I concentrations. In contrast, the principal complex with non-cognate DNA contained one duplex bound to Sfi I. Pairs of non-cognate duplexes, or one cognate and one non-cognate duplex, generally failed to form synaptic complexes. On adding Mg2+to complexes with cognate DNA, cleavage occurred much more rapidly in the synaptic complex than in the secondary complex. DNA synapsis thus acts to enhance the specificity of Sfi I for its recognition sequence, by demanding two cognate sites for a catalytically active complex and by excluding non-cognate sites from the synaptic complex.     

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.     

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.     

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.     

Crystallization of complexes of EcoRV endonuclease with cognate and non-cognate DNA fragments

Complexes of the type II restriction endonuclease EcoRV with a variety of short, selfcomplementary deoxyoligonucleotides have been crystallized. The best crystals diffract to about 2.7 A resolution and consist of 1:1 complexes between endonuclease dimers and duplexes of the cognate decamer GGGATATCCC containing the hexameric RV recognition sequence GATATC. Crystals with the non-cognate DNA octamer duplexes CGAGCTCG and CGAATTCG diffract to 3.0 and 3.5 A resolution, respectively, and contain two DNA duplexes per enzyme dimer.     

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.     

The elastic network model reveals a consistent picture on intrinsic functional dynamics of type II restriction endonucleases

The vibrational dynamics of various type II restriction endonucleases, in complex with cognate/non-cognate DNA and in the apo form, are investigated with the elastic network model in order to reveal common functional mechanisms in this enzyme family. Scissor-like and tong-like motions observed in the slowest modes of all enzymes and their complexes point to common DNA recognition and cleavage mechanisms. Normal mode analysis further points out that the scissor-like motion has an important role in differentiating between cognate and non-cognate sequences at the recognition site, thus implying its catalytic relevance. Flexible regions observed around the DNA-binding site of the enzyme usually concentrate on the highly conserved β-strands, especially after DNA binding. These β-strands may have a structurally stabilizing role in functional dynamics for target site recognition and cleavage. In addition, hot spot residues based on high-frequency modes reveal possible communication pathways between the two distant cleavage sites in the enzyme family. Some of these hot spots also exist on the shortest path between the catalytic sites and are highly conserved.     

Dynamic and Structural Changes in the Minimally Restructuring EcoRI Bound to a Minimally Mutated DNA Chain

The dynamics of a protein plays an important role in protein functionality. Here, we examine the differences in the dynamics of a minimally restructuring protein, EcoRI, when it is bound to its cognate DNA and to a noncognate sequence which differs by just a single basepair. Molecular dynamics simulations of the complexes and essential dynamics analyses reveal that the overall dynamics of the protein subunits change from a coordinated motion in the cognate complex to a scrambled motion in the noncognate complex. This dynamical difference extends to the protein-DNA interface where EcoRI tries to constrict the DNA in the cognate complex. In the noncognate complex, absence of the constricting motion of interfacial residues, overall change in backbone dynamics and structural relaxation of the arms enfolding the DNA leave the DNA less-kinked relative to the situation in the cognate complex, thus indicating that the protein is poised for linear diffusion along the DNA rather than for catalytic action. In a larger context, the results imply that the DNA sequences dictate protein dynamics and that when a protein chances upon the recognition sequence some of the key domains of the protein undergo dynamical changes that prepare the protein for eventual catalytic action.     

Two crystal forms of the restriction enzyme MspI-DNA complex show the same novel structure

The crystal structure of the Type IIP restriction endonuclease MspI bound to DNA containing its cognate recognition sequence has been determined in both monoclinic and orthorhombic space groups. Significantly, these two independent crystal forms present an identical structure of a novel monomer-DNA complex, suggesting a functional role for this novel enzyme-DNA complex. In both crystals, MspI interacts with the CCGG DNA recognition sequence as a monomer, using an asymmetric mode of recognition by two different structural motifs in a single polypeptide. In the crystallographic asymmetric unit, the two DNA molecules in the two MspI-DNA complexes appear to stack with each other forming an end-to-end pseudo-continuous 19-mer duplex. They are primarily B-form and no major bends or kinks are observed. For DNA recognition, most of the specific contacts between the enzyme and the DNA are preserved in the orthorhombic structure compared with the monoclinic structure. A cation is observed near the catalytic center in the monoclinic structure at a position homologous to the 74/45 metal site of EcoRV, and the orthorhombic structure also shows signs of this same cation. However, the coordination ligands of the metal are somewhat different from those of the 74/45 metal site of EcoRV. Combined with structural information from other solved structures of Type II restriction enzymes, the possible relationship between the structures of the enzymes and their cleavage behaviors is discussed.     

Fidelity of DNA sequence recognition by the SfiI restriction endonuclease is determined by communications between its two DNA-binding sites

The SfiI restriction endonuclease is a tetramer in which two subunits form a dimeric unit that contains one DNA binding cleft and the other two subunits contain a second cleft on the opposite side of the protein. Full activity requires both clefts to be filled with its recognition sequence: SfiI has low activity when bound to one site. The ability of SfiI to cleave non-cognate sites, one base pair different from the true site, was initially tested on substrates that lacked specific sites but which contained either one or multiple non-cognate sites. No cleavage of the DNA with one non-cognate site was detected, while a small fraction of the DNA with multiple sites was nicked. The alternative sequences were, however, cleaved in both strands, albeit at low levels, when the DNA also carried either a recognition site for SfiI or the termini generated by SfiI. Further tests employed a mutant of SfiI, altered at the dimer interface, which was known to be more active than wild-type SfiI when bound to a single site. This mutant similarly failed to cleave DNA with one non-cognate site, but cleaved the substrates with multiple non-cognate sites more readily than did the native enzyme. To cleave additional sites, SfiI thus needs to interact concurrently with either two non-cognate sites or one non-cognate and one cognate site (or the termini thereof), yet this arrangement is still restrained from cleaving the alternative site unless the communication pathway between the two DNA-binding clefts is disrupted.     

Dynamic and Structural Changes in the Minimally Restructuring EcoRI Bound to a Minimally Mutated DNA Chain

Abstract

The dynamics of a protein plays an important role in protein functionality. Here, we examine the differences in the dynamics of a minimally restructuring protein, EcoRI, when it is bound to its cognate DNA and to a noncognate sequence which differs by just a single basepair. Molecular dynamics simulations of the complexes and essential dynamics analyses reveal that the overall dynamics of the protein subunits change from a coordinated motion in the cognate complex to a scrambled motion in the noncognate complex. This dynamical difference extends to the protein-DNA interface where EcoRI tries to constrict the DNA in the cognate complex. In the noncognate complex, absence of the constricting motion of interfacial residues, overall change in backbone dynamics and structural relaxation of the arms enfolding the DNA leave the DNA less-kinked relative to the situation in the cognate complex, thus indicating that the protein is poised for linear diffusion along the DNA rather than for catalytic action. In a larger context, the results imply that the DNA sequences dictate protein dynamics and that when a protein chances upon the recognition sequence some of the key domains of the protein undergo dynamical changes that prepare the protein for eventual catalytic action.     

Mechanism of DNA Recognition by the Restriction Enzyme EcoRV

EcoRV, a restriction enzyme in Escherichia coli, destroys invading foreign DNA by cleaving it at the center step of a GATATC sequence. In the EcoRV-cognate DNA crystallographic complex, a sharp kink of 50° has been found at the center base-pair step (TA). Here, we examine the interplay between the intrinsic propensity of the cognate sequence to kink and the induction by the enzyme by performing all-atom molecular dynamics simulations of EcoRV unbound and interacting with three DNA sequences: the cognate sequence, GATATC (TA); the non-cognate sequence, GAATTC (AT); and with the cognate sequence methylated on the first adenine GA_CH3TATC (TA-CH₃). In the unbound EcoRV, the cleft between the two C-terminal subdomains is found to be open. Binding to AT narrows the cleft and forms a partially bound state. However, the intrinsic bending propensity of AT is insufficient to allow tight binding. In contrast, the cognate TA sequence is easier to bend, allowing specific, high-occupancy hydrogen bonds to form in the complex. The absence of cleavage for this methylated sequence is found to arise from the loss of specific hydrogen bonds between the first adenine of the recognition sequence and Asn185. On the basis of the results, we suggest a three-step recognition mechanism. In the first step, EcoRV, in an open conformation, binds to the DNA at a random sequence and slides along it. In the second step, when the two outer base pairs, GAxxTC, are recognized, the R loops of the protein become more ordered, forming strong hydrogen-bonding interactions, resulting in a partially bound EcoRV-DNA complex. In the third step, the flexibility of the center base pair is probed, and in the case of the full cognate sequence the DNA bends, the complex strengthens and the protein and DNA interact more closely, allowing cleavage.     

Specificity in protein-DNA interactions: energetic recognition by the (cytosine-C5)-methyltransferase from HhaI

Sequence-specific interactions between proteins and DNA are essential for a variety of biological functions. The (cytosine-C5)-methyltransferase from HhaI (M.HhaI) specifically modifies the second base in GCGC sequences, employing a base flipping mechanism to access the target base being chemically modified. The mechanism of sequence-specific recognition of M.HhaI is not evident based on crystallographic structures, leading to the suggestion that recognition is linked to the flipping event itself, a process that may be referred to as energetic recognition. Using computational methods, it is shown that the free energy barriers to flipping are significantly higher in non-cognate versus the cognate sequence, supporting the energetic recognition mechanism. Energetic recognition is imparted by two protein "selectivity filters" that function via a "web" of protein-DNA interactions in short-lived, high energy states present along the base flipping pathway. Other sequence-specific DNA binding proteins whose function involves significant distortion of DNA's conformation may use a similar recognition mechanism.