DNA distortion and specificity in a sequence-specific endonuclease

Five new structures of the Q138F HincII enzyme bound to a total of three different DNA sequences and three different metal ions (Ca²⁺, Mg²⁺, and Mn²⁺) are presented. While previous structures were produced from soaking Ca²⁺ into preformed Q138F HincII/DNA crystals, the new structures are derived from cocrystallization with Ca²⁺, Mg²⁺, or Mn²⁺. The Mn²⁺-bound structure provides the first view of a product complex of Q138F HincII with cleaved DNA. Binding studies and a crystal structure show how Ca²⁺ allows trapping of a Q138F HincII complex with noncognate DNA in a catalytically incompetent conformation. Many Q138F HincII/DNA structures show asymmetry, despite the binding of a symmetric substrate by a symmetric enzyme. The various complexes are fit into a model describing the different conformations of the DNA-bound enzyme and show how DNA conformational energetics determine DNA-cleavage rates by the Q138F HincII enzyme.     

Ca2+ binding in the active site of HincII: implications for the catalytic mechanism

The 2.8 A crystal structure of the type II restriction endonuclease HincII bound to Ca(2+) and cognate DNA containing GTCGAC is presented. The DNA is uncleaved, and one calcium ion is bound per active site, in a position previously described as site I in the related blunt cutting type II restriction endonuclease EcoRV [Horton, N. C., Newberry, K. J., and Perona, J. J. (1998) Proc. Natl. Acad. Sci. U.S.A. 95 (23), 13489-13494], as well as that found in other related enzymes. Unlike the site I metal in EcoRV, but similar to that of PvuII, NgoMIV, BamHI, BglII, and BglI, the observed calcium cation is directly ligated to the pro-S(p) oxygen of the scissile phosphate. A calcium ion-ligated water molecule is well positioned to act as the nucleophile in the phosphodiester bond cleavage reaction, and is within hydrogen bonding distance of the conserved active site lysine (Lys 129), as well as the pro-R(p) oxygen of the phosphate group 3' of the scissile phosphate, suggesting possible roles for these groups in the catalytic mechanism. Kinetic data consistent with an important role for the 3'-phosphate group in DNA cleavage by HincII are presented. The previously observed sodium ion [Horton, N. C., Dorner, L. F., and Perona, J. J. (2002) Nat. Struct. Biol. 9, 42-47] persists in the active sites of the Ca(2+)-bound structure; however, kinetic data show little effect on the single-turnover rate of DNA cleavage in the absence of Na(+) ions.     

Mechanistic insights from the structures of HincII bound to cognate DNA cleaved from addition of Mg2+ and Mn2+

The three-dimensional X-ray crystal structures of HincII bound to cognate DNA containing GTCGAC and Mn(2+) or Mg(2+), at 2.50A and 2.95A resolution, respectively, are presented. In both structures, the DNA is found cleaved, and the positions of the active-site groups, cleaved phosphate group, and 3' oxygen atom of the leaving group are in very similar positions. Two highly occupied Mn(2+) positions are found in each active site of the four crystallographically independent subunit copies in the HincII/DNA/Mn(2+) structure. The manganese ion closest to the previously identified single Ca(2+) position of HincII is shifted 1.7A and has lost direct ligation to the active-site aspartate residue, Asp127. A Mn(2+)-ligated water molecule in a position analogous to that seen in the HincII/DNA/Ca(2+) structure, and proposed to be the attacking nucleophile, is beyond hydrogen bonding distance from the active-site lysine residue, Lys129, but remains within hydrogen bonding distance from the proRp oxygen atom of the phosphate group 3' to the scissile phosphate group. In addition, the position of the cleaved phosphate group is on the opposite side of the axis connecting the two metal ions relative to that found in the BamHI/product DNA/Mn(2+) structure. Mechanistic implications are discussed, and a model for the two-metal-ion mechanism of DNA cleavage by HincII is proposed.     

Sequence selectivity and degeneracy of a restriction endonuclease mediated by DNA intercalation.

The crystal structure of the HincII restriction endonuclease-DNA complex shows that degenerate specificity for blunt-ended cleavage at GTPyPuAC sequences arises from indirect readout of conformational preferences at the center pyrimidine-purine step. Protein-induced distortion of the DNA is accomplished by intercalation of glutamine side chains into the major groove on either side of the recognition site, generating bending by either tilt or roll at three distinct loci. The intercalated side chains propagate a concerted shift of all six target-site base pairs toward the minor groove, producing an unusual cross-strand purine stacking at the center pyrimidine-purine step. Comparison of the HincII and EcoRV cocrystal structures suggests that sequence-dependent differences in base-stacking free energies are a crucial underlying factor mediating protein recognition by indirect readout.     

Rearrangement of side-chains in a Zif268 mutant highlights the complexities of zinc finger-DNA recognition.

Structural and biochemical studies of Cys(2)His(2) zinc finger proteins initially led several groups to propose a "recognition code" involving a simple set of rules relating key amino acid residues in the zinc finger protein to bases in its DNA site. One recent study from our group, involving geometric analysis of protein-DNA interactions, has discussed limitations of this idea and has shown how the spatial relationship between the polypeptide backbone and the DNA helps to determine what contacts are possible at any given position in a protein-DNA complex. Here we report a study of a zinc finger variant that highlights yet another source of complexity inherent in protein-DNA recognition. In particular, we find that mutations can cause key side-chains to rearrange at the protein-DNA interface without fundamental changes in the spatial relationship between the polypeptide backbone and the DNA. This is clear from a simple analysis of the binding site preferences and co-crystal structures for the Asp20-->Ala point mutant of Zif268. This point mutation in finger one changes the specificity of the protein from GCG TGG GCG to GCG TGG GC(G/T), and we have solved crystal structures of the D20A mutant bound to both types of sites. The structure of the D20A mutant bound to the GCG site reveals that contacts from key residues in the recognition helix are coupled in complex ways. The structure of the complex with the GCT site also shows an important new water molecule at the protein-DNA interface. These side-chain/side-chain interactions, and resultant changes in hydration at the interface, affect binding specificity in ways that cannot be predicted either from a simple recognition code or from analysis of spatial relationships at the protein-DNA interface. Accurate computer modeling of protein-DNA interfaces remains a challenging problem and will require systematic strategies for modeling side-chain rearrangements and change in hydration.     

Understanding the immutability of restriction enzymes: crystal structure of BglII and its DNA substrate at 1.5 A resolution

Restriction endonucleases are remarkably resilient to alterations in their DNA binding specificity. To understand the basis of this immutability, we have determined the crystal structure of endonuclease BglII bound to its recognition sequence (AGATCT), at 1. 5 A resolution. We compare the structure of BglII to endonuclease BamHI, which recognizes a closely related DNA site (GGATCC). We show that both enzymes share a similar alpha/beta core, but in BglII, the core is augmented by a beta-sandwich domain that encircles the DNA to provide extra specificity. Remarkably, the DNA is contorted differently in the two structures, leading to different protein-DNA contacts for even the common base pairs. Furthermore, the BglII active site contains a glutamine in place of the glutamate at the general base position in BamHI, and only a single metal is found coordinated to the putative nucleophilic water and the phosphate oxygens. This surprising diversity in structures shows that different strategies can be successful in achieving site-specific recognition and catalysis in restriction endonucleases.     

Specific DNA recognition by the Antp homeodomain: MD simulations of specific and nonspecific complexes

Four molecular dynamics simulation trajectories of complexes between the wild-type or a mutant Antennapedia homeodomain and 2 DNA sequences were generated in order to probe the mechanisms governing the specificity of DNA recognition. The starting point was published affinity measurements showing that a single protein mutation combined with a replacement of 2 base pairs yields a new high-affinity complex, whereas the other combinations, with changes on only 1 macromolecule, exhibited lower affinity. The simulations of the 4 complexes yielded fluctuating networks of interaction. On average, these networks differ significantly, explaining the switch of affinity caused by the alterations in the macromolecules. The network of mostly hydrogen-bonding interactions involving several water molecules, which was suggested both by X-ray and NMR structures of the wild-type homeodomain and its DNA operator sequence, could be reproduced in the trajectory. More interestingly, the high-affinity complex with alterations in both the protein and the DNA yielded again a dynamic but very tight network of intermolecular interactions, however, attributing a significantly stronger role to direct hydrophobic interactions at the expense of water bridges. The other 2 homeodomain-DNA complexes, with only 1 molecule altered, show on average over the trajectories a clearly reduced number of protein-DNA interactions. The observations from these simulations suggest specific experiments and thus close the circle formed by biochemical, structural, and computational studies. The shift from a water-dominated to a more "dry" interface may prove important in the design of proteins binding DNA in a specific manner.     

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.     

Engineered extrahelical base destabilization enhances sequence discrimination of DNA methyltransferase M.HhaI

Improved sequence specificity of the DNA cytosine methyltransferase HhaI was achieved by disrupting interactions at a hydrophobic interface between the active site of the enzyme and a highly conserved flexible loop. Transient fluorescence experiments show that mutations disrupting this interface destabilize the positioning of the extrahelical, "flipped" cytosine base within the active site. The ternary crystal structure of the F124A M.HhaI bound to cognate DNA and the cofactor analogue S-adenosyl-l-homocysteine shows an increase in cavity volume between the flexible loop and the core of the enzyme. This cavity disrupts the interface between the loop and the active site, thereby destabilizing the extrahelical target base. The favored partitioning of the base-flipped enzyme-DNA complex back to the base-stacked intermediate results in the mutant enzyme discriminating better than the wild-type enzyme against non-cognate sites. Building upon the concepts of kinetic proofreading and our understanding of M.HhaI, we describe how a 16-fold specificity enhancement achieved with a double mutation at the loop/active site interface is acquired through destabilization of intermediates prior to methyltransfer rather than disruption of direct interactions between the enzyme and the substrate for M.HhaI.     

Structure and calcium-binding studies of a recoverin mutant (E85Q) in an allosteric intermediate state

Recoverin, a member of the EF-hand superfamily, serves as a calcium sensor in retinal rod cells. A myristoyl or related fatty acyl group covalently attached to the N-terminus of recoverin facilitates the binding of recoverin to retinal disk membranes by a mechanism known as the Ca2+-myristoyl switch. Previous structural studies revealed that the myristoyl group of recoverin is sequestered inside the protein core in the absence of calcium. The cooperative binding of two calcium ions to the second and third EF-hands (EF-2 and EF-3) of recoverin leads to the extrusion of the fatty acid. Here we present nuclear magnetic resonance (NMR), fluorescence, and calcium-binding studies of a myristoylated recoverin mutant (myr-E85Q) designed to abolish high-affinity calcium binding to EF-2 and thereby trap the myristoylated protein with calcium bound solely to EF-3. Equilibrium calcium-binding studies confirm that only one Ca2+ binds to myr-E85Q under the conditions of this study with a dissociation constant of 100 microM. Fluorescence and NMR spectra of the Ca2+-free myr-E85Q are identical to those of Ca2+-free wild type, indicating that the E85Q mutation does not alter the stability and structure of the Ca2+-free protein. In contrast, the fluorescence and NMR spectra of half-saturated myr-E85Q (one bound Ca2+) look different from those of Ca2+-saturated wild type (two bound Ca2+), suggesting that half-saturated myr-E85Q may represent a structural intermediate. We report here the three-dimensional structure of Ca2+-bound myr-E85Q as determined by NMR spectroscopy. The N-terminal myristoyl group of Ca2+-bound myr-E85Q is sequestered within a hydrophobic cavity lined by many aromatic residues (F23, W31, Y53, F56, F83, and Y86) resembling that of Ca2+-free recoverin. The structure of Ca2+-bound myr-E85Q in the N-terminal region (residues 2-90) is similar to that of Ca2+-free recoverin, whereas the C-terminal region (residues 100-202) is more similar to that of Ca2+-bound wild type. Hence, the structure of Ca2+-bound myr-E85Q represents a hybrid between the structures of recoverin with zero and two Ca2+ bound. The binding of Ca2+ to EF-3 leads to local structural changes within the EF-hand that alter the domain interface and cause a 45 degrees swiveling of the N- and C-terminal domains, resulting in a partial unclamping of the myristoyl group. We propose that Ca2+-bound myr-E85Q may represent a stable intermediate state in the kinetic mechanism of the calcium-myristoyl switch.     

Contribution of cation-pi interactions to the stability of protein-DNA complexes

Cation-pi interactions between an aromatic ring and a positive charge located above it have proven to be important in protein structures and biomolecule associations. Here, the role of these interactions at the interface of protein-DNA complexes is investigated, by means of ab initio quantum mechanics energy calculations and X-ray structure analyses. Ab initio energy calculations indicate that Na ions and DNA bases can form stable cation-pi complexes, whose binding strength strongly depends on the type of base, on the position of the Na ion, and whether the base is isolated or included in a double-stranded B-DNA. A survey of protein-DNA complex structures using appropriate geometrical criteria revealed cation-pi interactions in 71% of the complexes. More than half of the cation-pi pairs involve arginine residues, about one-third asparagine or glutamine residues that only carry a partial charge, and one-seventh lysine residues. The most frequently observed pair, which is also the most stable as monitored by ab initio energy calculations, is arginine- guanine. Arginine-adenine interactions are also favorable in general, although to a lesser extent, whereas those with thymine and cytosine are not. Our calculations show that the major contribution to cation-pi interactions with DNA bases is of electrostatic nature. These interactions often occur concomitantly with hydrogen bonds with adjacent bases; their strength is estimated to be from three to four times lower than that of hydrogen bonds. Finally, the role of cation-pi interactions in the stability and specificity of protein-DNA complexes is discussed.     

Structural characterization of human 8-oxoguanine DNA glycosylase variants bearing active site mutations

The human 8-oxoguanine DNA glycosylase (hOGG1) protein is responsible for initiating base excision DNA repair of the endogenous mutagen 8-oxoguanine. Like nearly all DNA glycosylases, hOGG1 extrudes its substrate from the DNA helix and inserts it into an extrahelical enzyme active site pocket lined with residues that participate in lesion recognition and catalysis. Structural analysis has been performed on mutant versions of hOGG1 having changes in catalytic residues but not on variants having altered 7,8-dihydro-8-oxoguanine (oxoG) contact residues. Here we report high resolution structural analysis of such recognition variants. We found that Ala substitution at residues that contact the phosphate 5' to the lesion (H270A mutation) and its Watson-Crick face (Q315A mutation) simply removed key functionality from the contact interface but otherwise had no effect on structure. Ala substitution at the only residue making an oxoG-specific contact (G42A mutation) introduced torsional stress into the DNA contact surface of hOGG1, but this was overcome by local interactions within the folded protein, indicating that this oxoG recognition motif is "hardwired." Introduction of a side chain intended to sterically obstruct the active site pocket (Q315F mutation) led to two different structures, one of which (Q315F(*149)) has the oxoG lesion in an exosite flanking the active site and the other of which (Q315F(*292)) has the oxoG inserted nearly completely into the lesion recognition pocket. The latter structure offers a view of the latest stage in the base extrusion pathway yet observed, and its lack of catalytic activity demonstrates that the transition state for displacement of the lesion base is geometrically demanding.     

A new method for evaluating the specificity of indirect readout in protein-DNA recognition

Proteins recognize a specific DNA sequence not only through direct contact (direct readout) with base pairs but also through sequence-dependent conformation and/or flexibility of DNA (indirect readout). However, it is difficult to assess the contribution of indirect readout to the sequence specificity. What is needed is a straightforward method for quantifying its contributions to specificity. Using Bayesian statistics, we derived the probability of a particular sequence for a given DNA structure from the trajectories of molecular dynamics (MD) simulations of DNAs containing all possible tetramer sequences. Then, we quantified the specificity of indirect readout based on the information entropy associated with the probability. We tested this method with known structures of protein-DNA complexes. This method enabled us to correctly predict those regions where experiments suggested the involvement of indirect readout. The results also indicated new regions where the indirect readout mechanism makes major contributions to the recognition. The present method can be used to estimate the contribution of indirect readout without approximations to the distributions in the conformational ensembles of DNA, and would serve as a powerful tool to study the mechanism of protein-DNA recognition.     

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.     

The crystal structure of the gene targeting homing endonuclease I-SceI reveals the origins of its target site specificity

The I-SceI homing endonuclease enhances gene targeting by introducing double-strand breaks at specific chromosomal loci, thereby increasing the recombination frequency. Here, we report the crystal structure of the enzyme complexed to its DNA substrate and Ca(2+) determined at 2.25A resolution. The structure shows the prototypical beta-saddle of LAGLIDADG homing endonucleases that is contributed by two pseudo-symmetric domains. The high specificity of I-SceI is explained by the large number of protein-DNA contacts, many that are made by a long beta-hairpin loop that reaches into the major groove of the DNA. The DNA minor groove is compressed at the catalytic center, bringing the two scissile phosphodiester bonds into close proximity. The protein-Ca(2+)-DNA structure shows the protein bound to its DNA substrate in a pre-reactive state that is defined by the presence of two asymmetric active sites, one of which appears poised to first cleave the DNA bottom strand.