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

Sliding and jumping of single EcoRV restriction enzymes on non-cognate DNA

The restriction endonuclease EcoRV can rapidly locate a short recognition site within long non-cognate DNA using 'facilitated diffusion'. This process has long been attributed to a sliding mechanism, in which the enzyme first binds to the DNA via nonspecific interaction and then moves along the DNA by 1D diffusion. Recent studies, however, provided evidence that 3D translocations (hopping/jumping) also help EcoRV to locate its target site. Here we report the first direct observation of sliding and jumping of individual EcoRV molecules along nonspecific DNA. Using fluorescence microscopy, we could distinguish between a slow 1D diffusion of the enzyme and a fast translocation mechanism that was demonstrated to stem from 3D jumps. Salt effects on both sliding and jumping were investigated, and we developed numerical simulations to account for both the jump frequency and the jump length distribution. We deduced from our study the 1D diffusion coefficient of EcoRV, and we estimated the number of jumps occurring during an interaction event with nonspecific DNA. Our results substantiate that sliding alternates with hopping/jumping during the facilitated diffusion of EcoRV and, furthermore, set up a framework for the investigation of target site location by other DNA-binding proteins.

DNA conformations and their sequence preferences

The geometry of the phosphodiester backbone was analyzed for 7739 dinucleotides from 447 selected crystal structures of naked and complexed DNA. Ten torsion angles of a near-dinucleotide unit have been studied by combining Fourier averaging and clustering. Besides the known variants of the A-, B- and Z-DNA forms, we have also identified combined A + B backbone-deformed conformers, e.g. with alpha/gamma switches, and a few conformers with a syn orientation of bases occurring e.g. in G-quadruplex structures. A plethora of A- and B-like conformers show a close relationship between the A- and B-form double helices. A comparison of the populations of the conformers occurring in naked and complexed DNA has revealed a significant broadening of the DNA conformational space in the complexes, but the conformers still remain within the limits defined by the A- and B- forms. Possible sequence preferences, important for sequence-dependent recognition, have been assessed for the main A and B conformers by means of statistical goodness-of-fit tests. The structural properties of the backbone in quadruplexes, junctions and histone-core particles are discussed in further detail.

Modulation of Escherichia coli DNA methyltransferase activity by biologically derived GATC-flanking sequences

Escherichia coli DNA adenine methyltransferase (EcoDam) methylates the N-6 position of the adenine in the sequence 5'-GATC-3' and plays vital roles in gene regulation, mismatch repair, and DNA replication. It remains unclear how the small number of critical GATC sites involved in the regulation of replication and gene expression are differentially methylated, whereas the approximately 20,000 GATCs important for mismatch repair and dispersed throughout the genome are extensively methylated. Our prior work, limited to the pap regulon, showed that methylation efficiency is controlled by sequences immediately flanking the GATC sites. We extend these studies to include GATC sites involved in diverse gene regulatory and DNA replication pathways as well as sites previously shown to undergo differential in vivo methylation but whose function remains to be assigned. EcoDam shows no change in affinity with variations in flanking sequences derived from these sources, but methylation kinetics varied 12-fold. A-tracts immediately adjacent to the GATC site contribute significantly to these differences in methylation kinetics. Interestingly, only when the poly(A) is located 5' of the GATC are the changes in methylation kinetics revealed. Preferential methylation is obscured when two GATC sites are positioned on the same DNA molecule, unless both sites are surrounded by large amounts of nonspecific DNA. Thus, facilitated diffusion and sequences immediately flanking target sites contribute to higher order specificity for EcoDam; we suggest that the diverse biological roles of the enzyme are in part regulated by these two factors, which may be important for other enzymes that sequence-specifically modify DNA.

Cleavage of mispaired heteroduplex DNA substrates by numerous restriction enzymes

The utility of restriction endonucleases as a tool in molecular biology is in large part due to the high degree of specificity with which they cleave well-characterized DNA recognition sequences. The specificity of restriction endonucleases is not absolute, yet many commonly used assays of biological phenomena and contemporary molecular biology techniques rely on the premise that restriction enzymes will cleave only perfect cognate recognition sites. In vitro, mispaired heteroduplex DNAs are commonly formed, especially subsequent to polymerase chain reaction amplification. We investigated a panel of restriction endonucleases to determine their ability to cleave mispaired heteroduplex DNA substrates. Two straightforward, non-radioactive assays are used to evaluate mispaired heteroduplex DNA cleavage: a PCR amplification method and an oligonucleotide-based assay. These assays demonstrated that most restriction endonucleases are capable of site-specific double-strand cleavage with heteroduplex mispaired DNA substrates, however, certain mispaired substrates do effectively abrogate cleavage to undetectable levels. These data are consistent with mispaired substrate cleavage previously reported for Eco RI and, importantly, extend our knowledge of mispaired heteroduplex substrate cleavage to 13 additional enzymes.

Structural and evolutionary classification of Type II restriction enzymes based on theoretical and experimental analyses

For a very long time, Type II restriction enzymes (REases) have been a paradigm of ORFans: proteins with no detectable similarity to each other and to any other protein in the database, despite common cellular and biochemical function. Crystallographic analyses published until January 2008 provided high-resolution structures for only 28 of 1637 Type II REase sequences available in the Restriction Enzyme database (REBASE). Among these structures, all but two possess catalytic domains with the common PD-(D/E)XK nuclease fold. Two structures are unrelated to the others: R.BfiI exhibits the phospholipase D (PLD) fold, while R.PabI has a new fold termed 'half-pipe'. Thus far, bioinformatic studies supported by site-directed mutagenesis have extended the number of tentatively assigned REase folds to five (now including also GIY-YIG and HNH folds identified earlier in homing endonucleases) and provided structural predictions for dozens of REase sequences without experimentally solved structures. Here, we present a comprehensive study of all Type II REase sequences available in REBASE together with their homologs detectable in the nonredundant and environmental samples databases at the NCBI. We present the summary and critical evaluation of structural assignments and predictions reported earlier, new classification of all REase sequences into families, domain architecture analysis and new predictions of three-dimensional folds. Among 289 experimentally characterized (not putative) Type II REases, whose apparently full-length sequences are available in REBASE, we assign 199 (69%) to contain the PD-(D/E)XK domain. The HNH domain is the second most common, with 24 (8%) members. When putative REases are taken into account, the fraction of PD-(D/E)XK and HNH folds changes to 48% and 30%, respectively. Fifty-six characterized (and 521 predicted) REases remain unassigned to any of the five REase folds identified so far, and may exhibit new architectures. These enzymes are proposed as the most interesting targets for structure determination by high-resolution experimental methods. Our analysis provides the first comprehensive map of sequence-structure relationships among Type II REases and will help to focus the efforts of structural and functional genomics of this large and biotechnologically important class of enzymes.

A spectroscopic method to determine the activity of the restriction endonuclease EcoRV that involves a single reaction

A one-step protocol is presented to determine the activity of EcoRV as a model of restriction enzymes. The protocol involved a molecular beacon as DNA substrate, with the target sequence recognized by EcoRV in the stem. EcoRV cleaved the stem forming two fragments, one of which contained the fluorophore and quencher, initially bound by 3 bp. This shorter fragment rapidly dissociated at 37 °C, causing an increase of fluorescence intensity that was used to gauge the reaction kinetics. The reaction can be described using the Michaelis-Menten mechanism, and the kinetic parameters obtained were compared with literature values involving other protocols.

On-gel fluorescent visualization and the site identification of S-nitrosylated proteins

Mounting evidence indicates that S-nitrosylation of critical cysteine residues in a protein represents a common feature of protein function regulation and cell signaling. However, the progress in studying the exact role of S-nitrosylation has been hampered by a lack of rapid and accurate methods for the detection of these S-nitrosylated proteins and the exact modification sites. In this article, we report a fluorescence-based method in which the S-nitrosylated cysteines are converted into 7-amino-4-methylcoumarin-3-acetic acid (AMCA) fluorophore-labeled cysteines-termed the AMCA switch method. The labeled proteins are then analyzed by nonreducing SDS-PAGE, and the S-nitrosylated proteins can be readily detected as brilliant blue bands after the activation of ultraviolet light. Furthermore, the sites of modification can be determined by liquid chromatography-tandem mass spectrometry (LC-MS/MS) after in-gel tryptic digestion of the fluorescent band, and the recognizable AMCA tag in the MS spectra ensures the accurate site identification of the nitrosocysteines. Therefore, our method offers some apparent advantages by direct visualization of on-gel image of S-nitrosylated proteins, shorter experiment time by skipping the anti-biotin immunoblotting step, and elimination of the potential interference of endogenous biotinylated proteins. Based on this method, we detected the S-nitrosylation and the modified site in bovine serum albumin and gankyrin after in vitro S-nitrosylation. Overall, our results indicate that the AMCA switch method is a fast and accurate method to identify the S-nitrosylated protein and the modification sites.

Unexpected A-form formation of 4'-thioDNA in solution, revealed by NMR, and the implications as to the mechanism of nuclease resistance

Fully modified 4′-thioDNA, an oligonucleotide only comprising 2′-deoxy-4′-thionucleosides, exhibited resistance to an endonuclease, in addition to preferable hybridization with RNA. Therefore, 4′-thioDNA is promising for application as a functional oligonucleotide. Fully modified 4′-thioDNA was found to behave like an RNA molecule, but no details of its structure beyond the results of circular dichroism analysis are available. Here, we have determined the structure of fully modified 4′-thioDNA with the sequence of d(CGCGAATTCGCG) by NMR. Most sugars take on the C3′-endo conformation. The major groove is narrow and deep, while the minor groove is wide and shallow. Thus, fully modified 4′-thioDNA takes on the A-form characteristic of RNA, both locally and globally. The only structure reported for 4′-thioDNA showed that partially modified 4′-thioDNA that contained some 2′-deoxy-4′-thionucleosides took on the B-form in the crystalline form. We have determined the structure of 4′-thioDNA in solution for the first time, and demonstrated unexpected differences between the two structures. The origin of the formation of the A-form is discussed. The remarkable biochemical properties reported for fully modified 4′-thioDNA, including nuclease-resistance, are rationalized in the light of the elucidated structure.

A switch in the mechanism of communication between the two DNA-binding sites in the SfiI restriction endonuclease

While many Type II restriction enzymes are dimers with a single DNA-binding cleft between the subunits, SfiI is a tetramer of identical subunits. Two of its subunits (a dimeric unit) create one DNA-binding cleft, and the other two create a second cleft on the opposite side of the protein. The two clefts bind specific DNA cooperatively to give a complex of SfiI with two recognition sites. This complex is responsible for essentially all of the DNA-cleavage reactions by SfiI: virtually none is due to the complex with one site. The communication between the DNA-binding clefts was examined by disrupting one of the very few polar interactions in the otherwise hydrophobic interface between the dimeric units: a tyrosine hydroxyl was removed by mutation to phenylalanine. The mutant protein remained tetrameric in solution and could bind two DNA sites. But instead of being activated by binding two sites, like wild-type SfiI, it showed maximal activity when bound to a single site and had a lower activity when bound to two sites. This interaction across the dimer interface thus enforces in wild-type SfiI a cooperative transition between inactive and active states in both dimers, but without this interaction as in the mutant protein, a single dimer can undergo the transition to give a stable intermediate with one inactive dimer and one active dimer.

Use of plasmon coupling to reveal the dynamics of DNA bending and cleavage by single EcoRV restriction enzymes

Pairs of Au nanoparticles have recently been proposed as "plasmon rulers" based on the dependence of their light scattering on the interparticle distance. Preliminary work has suggested that plasmon rulers can be used to measure and monitor dynamic distance changes over the 1- to 100-nm length scale in biology. Here, we substantiate that plasmon rulers can be used to measure dynamical biophysical processes by applying the ruler to a system that has been investigated extensively by using ensemble kinetic measurements: the cleavage of DNA by the restriction enzyme EcoRV. Temporal resolutions of up to 240 Hz were obtained, and the end-to-end extension of up to 1,000 individual dsDNA enzyme substrates could be simultaneously monitored for hours. The kinetic parameters extracted from our single-molecule cleavage trajectories agree well with values obtained in bulk through other methods and confirm well known features of the cleavage process, such as DNA bending before cleavage. Previously unreported dynamical information is revealed as well, for instance, the degree of softening of the DNA just before cleavage. The unlimited lifetime, high temporal resolution, and high signal/noise ratio make the plasmon ruler a unique tool for studying macromolecular assemblies and conformational changes at the single-molecule level.

Structural basis for sequence-dependent DNA cleavage by nonspecific endonucleases

Nonspecific endonucleases hydrolyze DNA without sequence specificity but with sequence preference, however the structural basis for cleavage preference remains elusive. We show here that the nonspecific endonuclease ColE7 cleaves DNA with a preference for making nicks after (at 3'O-side) thymine bases but the periplasmic nuclease Vvn cleaves DNA more evenly with little sequence preference. The crystal structure of the 'preferred complex' of the nuclease domain of ColE7 bound to an 18 bp DNA with a thymine before the scissile phosphate had a more distorted DNA phosphate backbone than the backbones in the non-preferred complexes, so that the scissile phosphate was compositionally closer to the endonuclease active site resulting in more efficient DNA cleavage. On the other hand, in the crystal structure of Vvn in complex with a 16 bp DNA, the DNA phosphate backbone was similar and not distorted in comparison with that of a previously reported complex of Vvn with a different DNA sequence. Taken together these results suggest a general structural basis for the sequence-dependent DNA cleavage catalyzed by nonspecific endonucleases, indicating that nonspecific nucleases could induce DNA to deform to distinctive levels depending on the local sequence leading to different cleavage rates along the DNA chain.

Direct observation of DNA bending/unbending kinetics in complex with DNA-bending protein IHF

Regulation of gene expression involves formation of specific protein-DNA complexes in which the DNA is often bent or sharply kinked. Kinetics measurements of DNA bending when in complex with the protein are essential for understanding the molecular mechanism that leads to precise recognition of specific DNA-binding sites. Previous kinetics measurements on several DNA-bending proteins used stopped-flow techniques that have limited time resolution of few milliseconds. Here we use a nanosecond laser temperature-jump apparatus to probe, with submillisecond time resolution, the kinetics of bending/unbending of a DNA substrate bound to integration host factor (IHF), an architectural protein from Escherichia coli. The kinetics are monitored with time-resolved FRET, with the DNA substrates end-labeled with a FRET pair. The temperature-jump measurements, in combination with stopped-flow measurements, demonstrate that the binding of IHF to its cognate DNA site involves an intermediate state with straight or, possibly, partially bent DNA. The DNA bending rates range from approximately 2 ms(-1) at approximately 37 degrees C to approximately 40 ms(-1) at approximately 10 degrees C and correspond to an activation energy of approximately 14 +/- 3 kcal/mol. These rates and activation energy are similar to those of a single A:T base pair opening inside duplex DNA. Thus, our results suggest that spontaneous thermal disruption in base-paring, nucleated at an A:T site, may be sufficient to overcome the free energy barrier needed to partially bend/kink DNA before forming a tight complex with IHF.

Positively charged C-terminal subdomains of EcoRV endonuclease: contributions to DNA binding, bending, and cleavage

The carboxy-terminal subdomains of the homodimeric EcoRV restriction endonuclease each bear a net charge of +4 and are positioned on the inner concave surface of the 50 degree DNA bend that is induced by the enzyme. A complete kinetic and structural analysis of a truncated EcoRV mutant lacking these domains was performed to assess the importance of this diffuse charge in facilitating DNA binding, bending, and cleavage. At the level of formation of an enzyme-DNA complex, the association rate for the dimeric mutant enzyme was sharply decreased by 10(3)-fold, while the equilibrium dissociation constant was weakened by nearly 10(6)-fold compared with that of wild-type EcoRV. Thus, the C-terminal subdomains strongly stabilize the enzyme-DNA ground-state complex in which the DNA is known to be bent. Further, the extent of DNA bending as observed by fluorescence resonance energy transfer was also significantly decreased. The crystal structure of the truncated enzyme bound to DNA and calcium ions at 2.4 A resolution reveals that the global fold is preserved and suggests that a divalent metal ion crucial to catalysis is destabilized in the active site. This may explain the 100-fold decrease in the rate of metal-dependent phosphoryl transfer observed for the mutant. These results show that diffuse positive charge associated with the C-terminal subdomains of EcoRV plays a key role in DNA association, bending, and cleavage.

NMR structural and kinetic characterization of a homeodomain diffusing and hopping on nonspecific DNA

Nonspecific protein-DNA interactions are inherently dynamic and involve both diffusion of the protein along the DNA and hopping of the protein from one DNA molecule or segment to another. Understanding how gene regulatory proteins interact nonspecifically with DNA in terms of both structure and dynamics is challenging because the experimental observables are an ensemble average of many rapidly exchanging states. By using a variety of NMR spectroscopic techniques, including relaxation analysis, paramagnetic relaxation enhancement, and residual dipolar couplings, we have characterized structural and kinetic aspects of the interaction of the HoxD9 homeodomain with a nonspecific, 24-bp DNA duplex in a system in which the protein is not constrained to any particular site. The data reveal that HoxD9 binds to nonspecific DNA with the same binding mode and orientation as that observed in the specific complex. The mobility, however, of Arg side-chains contacting the DNA is increased in the nonspecific complex relative to the specific one. The kinetics of intermolecular translocation between two different nonspecific DNA molecules have also been analyzed and reveal that at high DNA concentrations (such as those present in vivo) direct transfer from one nonspecific complex to another nonspecific DNA molecule occurs without going through the intermediary of free protein. This finding provides a simple mechanism for accelerating the target search in vivo for the specific site in a sea of nonspecific sites by permitting more effective sampling of available DNA sites as the protein jumps from one segment to another.

Dynamics of single DNA looping and cleavage by Sau3AI and effect of tension applied to the DNA

Looping and cleavage of single DNA molecules by the two-site restriction endonuclease Sau3AI were measured with optical tweezers. A DNA template containing many recognition sites was used, permitting loop sizes from approximately 10 to 10,000 basepairs. At high enzyme concentration, cleavage events were detected within 5 s and nearly all molecules were cleaved within 5 min. Activity decreased approximately 10-fold as the DNA tension was increased from 0.03 to 0.7 pN. Substituting Ca(2+) for Mg(2+) blocked cleavage, permitting measurement of stable loops. At low tension, the initial rates of cleavage and looping were similar (approximately 0.025 s(-1) at 0.1 pN), suggesting that looping is rate limiting. Short loops formed more rapidly than long loops. The optimum size decreased from approximately 250 to 45 basepairs and the average number of loops (in 1 min) from 4.2 to 0.75 as tension was increased from 0.03 to 0.7 pN. No looping was detected at 5 pN. These findings are in qualitative agreement with recent theoretical predictions considering only DNA mechanics, but we observed weaker suppression with tension and smaller loop sizes. Our results suggest that the span and elasticity of the protein complex, nesting of loops, and protein-induced DNA bending and wrapping play an important role.

Sequence-dependent enhancement of hydrolytic deamination of cytosines in DNA by the restriction enzyme PspGI

Hydrolytic deamination of cytosines in DNA creates uracil and, if unrepaired, these lesions result in C to T mutations. We have suggested previously that a possible way in which cells may prevent or reduce this chemical reaction is through the binding of proteins to DNA. We use a genetic reversion assay to show that a restriction enzyme, PspGI, protects cytosines within its cognate site, 5'-CCWGG (W is A or T), against deamination under conditions where no DNA cleavage can occur. It decreases the rate of cytosine deamination to uracil by 7-fold. However, the same protein dramatically increases the rate of deaminations within the site 5'-CCSGG (S is C or G) by approximately 15-fold. Furthermore, a similar increase in cytosine deaminations is also seen with a catalytically inactive mutant of the enzyme showing that endonucleolytic ability of the protein is dispensable for its mutagenic action. The sequences of the mutants generated in the presence of PspGI show that only one of the cytosines in CCSGG is predominantly converted to thymine. Our results are consistent with PspGI 'sensitizing' the cytosine in the central base pair in CCSGG for deamination. Remarkably, PspGI sensitizes this base for damage despite its inability to form stable complexes at CCSGG sites. These results can be explained if the enzyme has a transient interaction with this sequence during which it flips the central cytosine out of the helix. This prediction was validated by modeling the structure of PspGI-DNA complex based on the structure of the related enzyme Ecl18kI which is known to cause base-flipping.

Tension-dependent DNA cleavage by restriction endonucleases: two-site enzymes are "switched off" at low force

DNA looping occurs in many important protein-DNA interactions, including those regulating replication, transcription, and recombination. Recent theoretical studies predict that tension of only a few piconewtons acting on DNA would almost completely inhibit DNA looping. Here, we study restriction endonucleases that require interaction at two separated sites for efficient cleavage. Using optical tweezers we measured the dependence of cleavage activity on DNA tension with 15 known or suspected two-site enzymes (BfiI, BpmI, BsgI, BspMI, Cfr9I, Cfr10I, Eco57I, EcoRII, FokI, HpaII, MboII, NarI, SacII, Sau3AI, and SgrAI) and six one-site enzymes (BamHI, EcoRI, EcoRV, HaeIII, HindIII, and DNaseI). All of the one-site enzymes were virtually unaffected by 5 pN of tension, whereas all of the two-site enzymes were completely inhibited. These enzymes thus constitute a remarkable example of a tension sensing "molecular switch." A detailed study of one enzyme, Sau3AI, indicated that the activity decreased exponentially with tension and the decrease was approximately 10-fold at 0.7 pN. At higher forces (approximately 20-40 pN) cleavage by the one-site enzymes EcoRV and HaeIII was partly inhibited and cleavage by HindIII was enhanced, whereas BamHI, EcoRI, and DNaseI were largely unaffected. These findings correlate with structural data showing that EcoRV bends DNA sharply, whereas BamHI, EcoRI, and DNaseI do not. Thus, DNA-directed enzyme activity involving either DNA looping or bending can be modulated by tension, a mechanism that could facilitate mechanosensory transduction in vivo.

DNA looping by two-site restriction endonucleases: heterogeneous probability distributions for loop size and unbinding force

Proteins interacting at multiple sites on DNA via looping play an important role in many fundamental biochemical processes. Restriction endonucleases that must bind at two recognition sites for efficient activity are a useful model system for studying such interactions. Here we used single DNA manipulation to study sixteen known or suspected two-site endonucleases. In eleven cases (BpmI, BsgI, BspMI, Cfr10I, Eco57I, EcoRII, FokI, HpaII, NarI, Sau3AI and SgrAI) we found that substitution of Ca2+ for Mg2+ blocked cleavage and enabled us to observe stable DNA looping. Forced disruption of these loops allowed us to measure the frequency of looping and probability distributions for loop size and unbinding force for each enzyme. In four cases we observed bimodal unbinding force distributions, indicating conformational heterogeneity and/or complex binding energy landscapes. Measured unlooping events ranged in size from 7 to 7500 bp and the most probable size ranged from less than 75 bp to nearly 500 bp, depending on the enzyme. In most cases the size distributions were in much closer agreement with theoretical models that postulate sharp DNA kinking than with classical models of DNA elasticity. Our findings indicate that DNA looping is highly variable depending on the specific protein and does not depend solely on the mechanical properties of DNA.

Changes in DNA bending and flexing due to tethered cations detected by fluorescence resonance energy transfer

Local DNA deformation arises from an interplay among sequence-related base stacking, intrastrand phosphate repulsion, and counterion and water distribution, which is further complicated by the approach and binding of a protein. The role of electrostatics in this complex chemistry was investigated using tethered cationic groups that mimic proximate side chains. A DNA duplex was modified with one or two centrally located deoxyuracils substituted at the 5-position with either a flexible 3-aminopropyl group or a rigid 3-aminopropyn-1-yl group. End-to-end helical distances and duplex flexibility were obtained from measurements of the time-resolved Förster resonance energy transfer between 5′- and 3′-linked dye pairs. A novel analysis utilized the first and second moments of the G(t) function, which encompasses only the energy transfer process. Duplex flexibility is altered by the presence of even a single positive charge. In contrast, the mean 5′–3′ distance is significantly altered by the introduction of two adjacently tethered cations into the double helix but not by a single cation: two adjacent aminopropyl groups decrease the 5′–3′ distance while neighboring aminopropynyl groups lengthen the helix.

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.

A view of consecutive binding events from structures of tetrameric endonuclease SfiI bound to DNA

Many reactions in cells proceed via the sequestration of two DNA molecules in a synaptic complex. SfiI is a member of a growing family of restriction enzymes that can bind and cleave two DNA sites simultaneously. We present here the structures of tetrameric SfiI in complex with cognate DNA. The structures reveal two different binding states of SfiI: one with both DNA-binding sites fully occupied and the other with fully and partially occupied sites. These two states provide details on how SfiI recognizes and cleaves its target DNA sites, and gives insight into sequential binding events. The SfiI recognition sequence (GGCCNNNN[downward arrow]NGGCC) is a subset of the recognition sequence of BglI (GCCNNNN[downward arrow]NGGC), and both enzymes cleave their target DNAs to leave 3-base 3' overhangs. We show that even though SfiI is a tetramer and BglI is a dimer, and there is little sequence similarity between the two enzymes, their modes of DNA recognition are unusually similar.

Interfacial water as a "hydration fingerprint" in the noncognate complex of BamHI

The molecular code of specific DNA recognition by proteins as a paradigm in molecular biology remains an unsolved puzzle primarily because of the subtle interplay between direct protein-DNA interaction and the indirect contribution from water and ions. Transformation of the nonspecific, low affinity complex to a specific, high affinity complex is accompanied by the release of interfacial water molecules. To provide insight into the conversion from the loose to the tight form, we characterized the structure and energetics of water at the protein-DNA interface of the BamHI complex with a noncognate sequence and in the specific complex. The fully hydrated models were produced with Grand Canonical Monte Carlo simulations. Proximity analysis shows that water distributions exhibit sequence dependent variations in both complexes and, in particular, in the noncognate complex they discriminate between the correct and the star site. Variations in water distributions control the number of water molecules released from a given sequence upon transformation from the loose to the tight complex as well as the local entropy contribution to the binding free energy. We propose that interfacial waters can serve as a "hydration fingerprint" of a given DNA sequence.

Common patterns in type II restriction enzyme binding sites

Restriction enzymes are among the best studied examples of DNA binding proteins. In order to find general patterns in DNA recognition sites, which may reflect important properties of protein–DNA interaction, we analyse the binding sites of all known type II restriction endonucleases. We find a significantly enhanced GC content and discuss three explanations for this phenomenon. Moreover, we study patterns of nucleotide order in recognition sites. Our analysis reveals a striking accumulation of adjacent purines (R) or pyrimidines (Y). We discuss three possible reasons: RR/YY dinucleotides are characterized by (i) stronger H-bond donor and acceptor clusters, (ii) specific geometrical properties and (iii) a low stacking energy. These features make RR/YY steps particularly accessible for specific protein–DNA interactions. Finally, we show that the recognition sites of type II restriction enzymes are underrepresented in host genomes and in phage genomes.

DNA-tension dependence of restriction enzyme activity reveals mechanochemical properties of the reaction pathway

Type II restriction endonucleases protect bacteria against phage infections by cleaving recognition sites on foreign double-stranded DNA (dsDNA) with extraordinary specificity. This capability arises primarily from large conformational changes in enzyme and/or DNA upon target sequence recognition. In order to elucidate the connection between the mechanics and the chemistry of DNA recognition and cleavage, we used a single-molecule approach to measure rate changes in the reaction pathway of EcoRV and BamHI as a function of DNA tension. We show that the induced-fit rate of EcoRV is strongly reduced by such tension. In contrast, BamHI is found to be insensitive, providing evidence that both substrate binding and hydrolysis are not influenced by this force. Based on these results, we propose a mechanochemical model of induced-fit reactions on DNA, allowing determination of induced-fit rates and DNA bend angles. Finally, for both enzymes a strongly decreased association rate is obtained on stretched DNA, presumably due to the absence of intradomain dissociation/re-association between non-specific sites (jumping). The obtained results should apply to many other DNA-associated proteins.

Structure of HinP1I endonuclease reveals a striking similarity to the monomeric restriction enzyme MspI

HinP1I, a type II restriction endonuclease, recognizes and cleaves a palindromic tetranucleotide sequence (G↓CGC) in double-stranded DNA, producing 2 nt 5′ overhanging ends. Here, we report the structure of HinP1I crystallized as one protein monomer in the crystallographic asymmetric unit. HinP1I displays an elongated shape, with a conserved catalytic core domain containing an active-site motif of SDX₁₈QXK and a putative DNA-binding domain. Without significant sequence homology, HinP1I displays striking structural similarity to MspI, an endonuclease that cleaves a similar palindromic DNA sequence (C↓CGG) and binds to that sequence crystallographically as a monomer. Almost all the structural elements of MspI can be matched in HinP1I, including both the DNA recognition and catalytic elements. Examining the protein–protein interactions in the crystal lattice, HinP1I could be dimerized through two helices located on the opposite side of the protein to the active site, generating a molecule with two active sites and two DNA-binding surfaces opposite one another on the outer surfaces of the dimer. A possible functional link between this unusual dimerization mode and the tetrameric restriction enzymes is discussed.

A homology model of restriction endonuclease SfiI in complex with DNA

Background

Restriction enzymes (REases) are commercial reagents commonly used in recombinant DNA technologies. They are attractive models for studying protein-DNA interactions and valuable targets for protein engineering. They are, however, extremely divergent: the amino acid sequence of a typical REase usually shows no detectable similarities to any other proteins, with rare exceptions of other REases that recognize identical or very similar sequences. From structural analyses and bioinformatics studies it has been learned that some REases belong to at least four unrelated and structurally distinct superfamilies of nucleases, PD-DxK, PLD, HNH, and GIY-YIG. Hence, they are extremely hard targets for structure prediction and homology-based inference of sequence-function relationships and the great majority of REases remain structurally and evolutionarily unclassified.

Results

SfiI is a REase which recognizes the interrupted palindromic sequence 5'GGCCNNNN^NGGCC3' and generates 3 nt long 3' overhangs upon cleavage. SfiI is an archetypal Type IIF enzyme, which functions as a tetramer and cleaves two copies of the recognition site in a concerted manner. Its sequence shows no similarity to other proteins and nothing is known about the localization of its active site or residues important for oligomerization. Using the threading approach for protein fold-recognition, we identified a remote relationship between SfiI and BglI, a dimeric Type IIP restriction enzyme from the PD-DxK superfamily of nucleases, which recognizes the 5'GCCNNNN^NGGC3' sequence and whose structure in complex with the substrate DNA is available. We constructed a homology model of SfiI in complex with its target sequence and used it to predict residues important for dimerization, tetramerization, DNA binding and catalysis.

Conclusions

The bioinformatics analysis suggest that SfiI, a Type IIF enzyme, is more closely related to BglI, an "orthodox" Type IIP restriction enzyme, than to any other REase, including other Type IIF REases with known structures, such as NgoMIV. NgoMIV and BglI belong to two different, very remotely related branches of the PD-DxK superfamily: the α-class (EcoRI-like), and the β-class (EcoRV-like), respectively. Thus, our analysis provides evidence that the ability to tetramerize and cut the two DNA sequences in a concerted manner was developed independently at least two times in the evolution of the PD-DxK superfamily of REases. The model of SfiI will also serve as a convenient platform for further experimental analyses.

Kinetics of protein-DNA interaction: facilitated target location in sequence-dependent potential

Recognition and binding of specific sites on DNA by proteins is central for many cellular functions such as transcription, replication, and recombination. In the process of recognition, a protein rapidly searches for its specific site on a long DNA molecule and then strongly binds this site. Here we aim to find a mechanism that can provide both a fast search (1-10 s) and high stability of the specific protein-DNA complex (Kd=10(-15)-10(-8) M). Earlier studies have suggested that rapid search involves sliding of the protein along the DNA. Here we consider sliding as a one-dimensional diffusion in a sequence-dependent rough energy landscape. We demonstrate that, despite the landscape's roughness, rapid search can be achieved if one-dimensional sliding is accompanied by three-dimensional diffusion. We estimate the range of the specific and nonspecific DNA-binding energy required for rapid search and suggest experiments that can test our mechanism. We show that optimal search requires a protein to spend half of its time sliding along the DNA and the other half diffusing in three dimensions. We also establish that, paradoxically, realistic energy functions cannot provide both rapid search and strong binding of a rigid protein. To reconcile these two fundamental requirements we propose a search-and-fold mechanism that involves the coupling of protein binding and partial protein folding. The proposed mechanism has several important biological implications for search in the presence of other proteins and nucleosomes, simultaneous search by several proteins, etc. The proposed mechanism also provides a new framework for interpretation of experimental and structural data on protein-DNA interactions.

How do site-specific DNA-binding proteins find their targets?

Essentially all the biological functions of DNA depend on site-specific DNA-binding proteins finding their targets, and therefore 'searching' through megabases of non-target DNA. In this article, we review current understanding of how this sequence searching is done. We review how simple diffusion through solution may be unable to account for the rapid rates of association observed in experiments on some model systems, primarily the Lac repressor. We then present a simplified version of the 'facilitated diffusion' model of Berg, Winter and von Hippel, showing how non-specific DNA-protein interactions may account for accelerated targeting, by permitting the protein to sample many binding sites per DNA encounter. We discuss the 1-dimensional 'sliding' motion of protein along non-specific DNA, often proposed to be the mechanism of this multiple site sampling, and we discuss the role of short-range diffusive 'hopping' motions. We then derive the optimal range of sliding for a few physical situations, including simple models of chromosomes in vivo, showing that a sliding range of approximately 100 bp before dissociation optimizes targeting in vivo. Going beyond first-order binding kinetics, we discuss how processivity, the interaction of a protein with two or more targets on the same DNA, can reveal the extent of sliding and we review recent experiments studying processivity using the restriction enzyme EcoRV. Finally, we discuss how single molecule techniques might be used to study the dynamics of DNA site-specific targeting of proteins.

Mycobacterium tuberculosis RecA intein, a LAGLIDADG homing endonuclease, displays Mn(2+) and DNA-dependent ATPase activity

Mycobacterium tuberculosis RecA intein (PI-MtuI), a LAGLIDADG homing endonuclease, displays dual target specificity in response to alternative cofactors. While both ATP and Mn(2+) were required for optimal cleavage of an inteinless recA allele (hereafter referred to as cognate DNA), Mg(2+) alone was sufficient for cleavage of ectopic DNA sites. In this study, we have explored the ability of PI-MtuI to catalyze ATP hydrolysis in the presence of alternative metal ion cofactors and DNA substrates. Our results indicate that PI-MtuI displays maximum ATPase activity in the presence of cognate but not ectopic DNA. Kinetic analysis revealed that Mn(2+) was able to stimulate PI-MtuI catalyzed ATP hydrolysis, whereas Mg(2+) failed to do so. Using UV crosslinking, limited proteolysis and amino acid sequence analysis, we show that (32)P-labeled ATP was bound to a 14 kDa peptide containing the putative Walker A motif. Furthermore, the limited proteolysis approach disclosed that cognate DNA was able to induce structural changes in PI-MtuI. Mutation of the presumptive metal ion-binding ligands (Asp122 and Asp222) in the LAGLIDADG motifs of PI-MtuI impaired its affinity for ATP, thus resulting in a reduction in or loss of its endonuclease activity. Together, these results suggest that PI-MtuI is a (cognate) DNA- and Mn(2+)-dependent ATPase, unique from the LAGLIDADG family of homing endonucleases, and implies a possible role for ATP hydrolysis in the recognition and/or cleavage of homing site DNA sequence.

Energetic and structural considerations for the mechanism of protein sliding along DNA in the nonspecific BamHI-DNA complex

The molecular mechanism by which DNA-binding proteins find their specific binding sites is still unclear. To gain insights into structural and energetic elements of this mechanism, we used the crystal structure of the nonspecific BamHI-DNA complex as a template to study the dominant electrostatic interaction in the nonspecific association of protein with DNA, and the possible sliding pathways that could be sustained by such an interaction. Based on calculations using the nonlinear Poisson-Boltzmann method and Brownian dynamics, a model is proposed for the initial nonspecific binding of BamHI to B-form DNA that differs from that seen in the crystal structure of the nonspecific complex. The model is electrostatically favorable, and the salt dependence as well as other thermodynamic parameters calculated for this model are in good agreement with experimental results. Several residues in BamHI are identified for their important contribution to the energy in the nonspecific binding model, and specific mutagenesis experiments are proposed to test the model on this basis. We show that a favorable sliding pathway of the protein along DNA is helical.

Protein motion from non-specific to specific DNA by three-dimensional routes aided by supercoiling

DNA-binding proteins are generally thought to locate their target sites by first associating with the DNA at random and then translocating to the specific site by one-dimensional (1D) diffusion along the DNA. We report here that non-specific DNA conveys proteins to their target sites just as well when held near the target by catenation as when co-linear with the target. Hence, contrary to the prevalent view, proteins move from random to specific sites primarily by three-dimensional (3D) rather than 1D pathways, by multiple dissociation/re-association events within a single DNA molecule. We also uncover a role for DNA supercoiling in target-site location. Proteins find their sites more readily in supercoiled than in relaxed DNA, again indicating 3D rather than 1D routes.

Baculovirus alkaline nuclease possesses a 5'-->3' exonuclease activity and associates with the DNA-binding protein LEF-3

Alkaline nuclease (AN) of the Autographa californica multiple-capsid nucleopolyhedrovirus (AcMNPV) (open reading frame 133) was expressed in recombinant baculovirus as a His(6)-tagged fusion and purified by sequential chromatography on Ni-NTA-agarose, DEAE-Toyopearl, and heparin-Sepharose. At all stages of purification, AcMNPV AN was found to copurify with a 44-kDa polypeptide which was identified as the baculovirus single-stranded DNA (ssDNA)-binding (SSB) protein, LEF-3. Sedimentation analysis in glycerol gradients of highly purified samples suggested that AN and LEF-3 are associated in a complex (designated *AN/L3), predominantly as heterodimers, although oligomeric forms containing both proteins were evident. In reactions with single- or double-stranded 62-mer oligonucleotides that were labeled with (32)P at the 5' or 3' ends, *AN/L3 carried out exonucleolytic hydrolysis of both substrates exclusively in a 5'-->3' direction. Saturation of ssDNA with an excess of LEF-3 prior to the addition of *AN/L3 resulted in a marked decrease in the rate of ssDNA hydrolysis. This suggests that excess LEF-3 may protect ssDNA from digestion by a AN-LEF-3 complex, thus regulating its activity in infected cells. The association of baculovirus AN with the viral SSB LEF-3 and the 5'-->3' exonuclease activity of this complex suggests that AN and LEF-3 may participate in homologous recombination of the baculovirus genome in a manner similar to that of exonuclease (Redalpha) and DNA-binding protein (Redbeta) of the Red-mediated homologous recombination system of bacteriophage lambda.

Protein stability indicates divergent evolution of PD-(D/E)XK type II restriction endonucleases

Type II restriction endonucleases recognize 4-8 base-pair-long DNA sequences and catalyze their cleavage with remarkable specificity. Crystal structures of the PD-(DE)XK superfamily revealed a common alpha/beta core motif and similar active site. In contrast, these enzymes show little sequence similarity and use different strategies to interact with their substrate DNA. The intriguing question is whether this enzyme family could have evolved from a common origin. In our present work, protein structure stability elements were analyzed and compared in three parts of PD-(DE)XK type II restriction endonucleases: (1) core motif, (2) active-site residues, and (3) residues playing role in DNA recognition. High correlation was found between the active-site residues and those stabilization factors that contribute to preventing structural decay. DNA recognition sites were also observed to participate in stabilization centers. It indicates that recognition motifs and active sites in PD-(DE)XK type II restriction endonucleases should have been evolutionary more conserved than other parts of the structure. Based on this observation it is proposed that PD-(DE)XK type II restriction endonucleases have developed from a common ancestor with divergent evolution.

Involvement of conserved histidine, lysine and tyrosine residues in the mechanism of DNA cleavage by the caspase-3 activated DNase CAD

The caspase-activated DNase (CAD) is involved in DNA degradation during apoptosis. Chemical modification of murine CAD with the lysine-specific reagent 2,4,6-trinitrobenzenesulphonic acid and the tyrosine-specific reagent N-acetylimidazole leads to inactivation of the nuclease, indicating that lysine and tyrosine residues are important for DNA cleavage by this enzyme. The presence of DNA or the inhibitor ICAD-L protects the enzyme from modification. Amino acid substitution in murine CAD of lysines and tyrosines conserved in CADs from five different species leads to variants with little if any catalytic activity, but unaltered DNA binding (K155Q, K301Q, K310Q, Y247F), with the exception of Y170F, which retains wild-type activity. Similarly, as observed for the previously characterised H242N, H263N, H308N and H313N variants, the newly introduced His-->Asp/Glu or Arg exchanges lead to variants with <1% of wild-type activity, with two exceptions: H313R shows wild-type activity, and H308D at pH 5.0 exhibits approximately 5% of wild-type activity at this pH. Y170F and H313R produce a specific pattern of fragments, different from wild-type CAD, which degrades DNA non-specifically. The recombinant nuclease variants produced in Escherichia coli were tested for their ability to form nucleolytically active oligomers. They did not show any significant deviation from the wild-type enzyme. Based on these and published data possible roles of the amino acid residues under investigation are discussed.

Crystal structure of the Bse634I restriction endonuclease: comparison of two enzymes recognizing the same DNA sequence

Crystal structures of Type II restriction endonucleases demonstrate a conserved common core and active site residues but diverse structural elements involved in DNA sequence discrimination. Comparative structural analysis of restriction enzymes recognizing the same nucleotide sequence might therefore contribute to our understanding of the structural diversity of specificity determinants within restriction enzymes. We have solved the crystal structure of the Bacillus stearothermophilus restriction endonuclease Bse634I by the multiple isomorphous replacement technique to 2.17 A resolution. Bse634I is an isoschisomer of the Cfr10I restriction enzyme whose crystal structure has been reported previously. Comparative structural analysis of the first pair of isoschisomeric enzymes revealed conserved structural determinants of sequence recognition and catalysis. However, conformations of the N-terminal subdomains differed between Bse634I/Cfr10I, suggesting a rigid body movement that might couple DNA recognition and catalysis. Structural similarities extend to the quaternary structure level: crystal contacts suggest that Bse634I similarly to Cfr10I is arranged as a tetramer. Kinetic analysis reveals that Bse634I is able to interact simultaneously with two recognition sites supporting the tetrameric architecture of the protein. Thus, restriction enzymes Bse634I, Cfr10I and NgoMIV, recognizing overlapping nucleotide sequences, exhibit a conserved tetrameric architecture that is of functional importance.

Restriction enzymes increase efficiencies of illegitimate DNA integration but decrease homologous integration in mammalian cells

Mammalian cells repair DNA double-strand breaks by illegitimate end-joining or by homologous recombination. We investigated the effects of restriction enzymes on illegitimate and homologous DNA integration in mammalian cells. A plasmid containing the neo(R) expression cassette, which confers G418 resistance, was used to select for illegitimate integration events in CHO wild-type and xrcc5 mutant cells. Co-transfection with the restriction enzymes BamHI, BglII, EcoRI and KpnI increased the efficiency of linearized plasmid integration up to 5-fold in CHO cells. In contrast, the restriction enzymes did not increase the integration efficiency in xrcc5 mutant cells. Effects of restriction enzymes on illegitimate and homologous integration were also studied in mouse embryonic stem (ES) cells using a plasmid containing the neo(R) gene flanked by exon 3 of HPRT: The enzymes BamHI, BglII and EcoRI increased the illegitimate integration efficiency of transforming DNA several-fold, similar to the results for CHO cells. However, all three enzymes decreased the absolute frequency of homologous integration approximately 2-fold, and the percentage of homologous integration decreased >10-fold. This suggests that random DNA breaks attract illegitimate recombination (IR) events that compete with homology search.

The structural basis of damaged DNA recognition and endonucleolytic cleavage for very short patch repair endonuclease

Endonucleases in DNA repair must be able to recognize damaged DNA as well as cleave the phosphodiester backbone. These functional prerequisites are manifested in very short patch repair (Vsr) endonuclease through a common endonuclease topology that has been tailored for recognition of TG mismatches. Structural and biochemical comparison with type II restriction enzymes illustrates how Vsr resembles these endonucleases in overall topology but also how Vsr diverges in terms of the detailed catalytic mechanism. A histidine and two metal-water clusters catalyze the phosphodiester cleavage. The mode of DNA damage recognition is also unique to Vsr. All other structurally characterized DNA damage-binding enzymes employ a nucleotide flipping mechanism for substrate recognition and for catalysis. Vsr, on the other hand, recognizes the TG mismatch as a wobble base pair and penetrates the DNA with three aromatic residues on one side of the mismatch. Thus, Vsr endonuclease provides important counterpoints in our understanding of endonucleolytic mechanisms and of damaged DNA recognition.

Homing endonucleases: structural and functional insight into the catalysts of intron/intein mobility

Homing endonucleases confer mobility to their host intervening sequence, either an intron or intein, by catalyzing a highly specific double-strand break in a cognate allele lacking the intervening sequence. These proteins are characterized by their ability to bind long DNA target sites (14-40 bp) and their tolerance of minor sequence changes in these sites. A wealth of biochemical and structural data has been generated for these enzymes over the past few years. Herein we review our current understanding of homing endonucleases, including their diversity and evolution, DNA-binding and catalytic mechanisms, and attempts to engineer them to bind novel DNA substrates.

Structure and function of type II restriction endonucleases

More than 3000 type II restriction endonucleases have been discovered. They recognize short, usually palindromic, sequences of 4-8 bp and, in the presence of Mg(2+), cleave the DNA within or in close proximity to the recognition sequence. The orthodox type II enzymes are homodimers which recognize palindromic sites. Depending on particular features subtypes are classified. All structures of restriction enzymes show a common structural core comprising four beta-strands and one alpha-helix. Furthermore, two families of enzymes can be distinguished which are structurally very similar (EcoRI-like enzymes and EcoRV-like enzymes). Like other DNA binding proteins, restriction enzymes are capable of non-specific DNA binding, which is the prerequisite for efficient target site location by facilitated diffusion. Non-specific binding usually does not involve interactions with the bases but only with the DNA backbone. In contrast, specific binding is characterized by an intimate interplay between direct (interaction with the bases) and indirect (interaction with the backbone) readout. Typically approximately 15-20 hydrogen bonds are formed between a dimeric restriction enzyme and the bases of the recognition sequence, in addition to numerous van der Waals contacts to the bases and hydrogen bonds to the backbone, which may also be water mediated. The recognition process triggers large conformational changes of the enzyme and the DNA, which lead to the activation of the catalytic centers. In many restriction enzymes the catalytic centers, one in each subunit, are represented by the PD. D/EXK motif, in which the two carboxylates are responsible for Mg(2+) binding, the essential cofactor for the great majority of enzymes. The precise mechanism of cleavage has not yet been established for any enzyme, the main uncertainty concerns the number of Mg(2+) ions directly involved in cleavage. Cleavage in the two strands usually occurs in a concerted fashion and leads to inversion of configuration at the phosphorus. The products of the reaction are DNA fragments with a 3'-OH and a 5'-phosphate.

SfiI endonuclease activity is strongly influenced by the non-specific sequence in the middle of its recognition site

The SfiI endonuclease cleaves DNA at the sequence GGCCNNNN NGGCC, where N is any base and downward arrow is the point of cleavage. Proteins that recognise discontinuous sequences in DNA can be affected by the unspecified sequence between the specified base pairs of the target site. To examine whether this applies to SFII, a series of DNA duplexes were made with identical sequences apart from discrete variations in the 5 bp spacer. The rates at which SFII cleaved each duplex were measured under steady-state conditions: the steady-state rates were determined by the DNA cleavage step in the reaction pathway. SFII cleaved some of these substrates at faster rates than other substrates. For example, the change in spacer sequence from AACAA to AAACA caused a 70-fold increase in reaction rate. In general, the extrapolated values for k(cat) and K(m) were both higher on substrates with inflexible spacers than those with flexible structures. The dinucleotide at the site of cleavage was largely immaterial. SFII activity is thus highly dependent on conformational variations in the spacer DNA.

Calculation of weak protein-protein interactions: the pH dependence of the second virial coefficient

Interactions between proteins are often sufficiently weak that their study through the use of conventional structural techniques becomes problematic. Of the few techniques capable of providing experimental measures of weak protein-protein interactions, perhaps the most useful is the second virial coefficient, B(22), which quantifies a protein solution's deviations from ideal behavior. It has long been known that B(22) can in principle be computed, but only very recently has it been demonstrated that such calculations can be performed using protein models of true atomic detail (Biophys. J. 1998, 75:2469-2477). The work reported here extends these previous efforts in an attempt to develop a transferable energetic model capable of reproducing the experimental trends obtained for two different proteins over a range of pH and ionic strengths. We describe protein-protein interaction energies by a combination of three separate terms: (i) an electrostatic interaction term based on the use of effective charges, (ii) a term describing the electrostatic desolvation that occurs when charged groups are buried by an approaching protein partner, and (iii) a solvent-accessible surface area term that is used to describe contributions from van der Waals and hydrophobic interactions. The magnitude of the third term is governed by an adjustable, empirical parameter, gamma, that is altered to optimize agreement between calculated and experimental values of B(22). The model is applied separately to the proteins lysozyme and chymotrypsinogen, yielding optimal values of gamma that are almost identical. There are, however, clear difficulties in reproducing B(22) values at the extremes of pH. Explicit calculation of the protonation states of ionizable amino acids in the 200 most energetically favorable protein-protein structures suggest that these difficulties are due to a neglect of the protonation state changes that can accompany complexation. Proper reproduction of the pH dependence of B(22) will, therefore, almost certainly require that account be taken of these protonation state changes. Despite this problem, the fact that almost identical gamma values are obtained from two different proteins suggests that the basic energetic formulation used here, which can be evaluated very rapidly, might find use in dynamical simulations of weak protein-protein interactions at intermediate pH values.

Sequence-specific recognition and cleavage of telomeric repeat (TTAGG)(n) by endonuclease of non-long terminal repeat retrotransposon TRAS1

The telomere of the silkworm Bombyx mori consists of (TTAGG/CCTAA)(n) repeats and harbors a large number of telomeric repeat-specific non-long terminal repeat retrotransposons, such as TRAS1 and SART1. To understand how these retrotransposons recognize and integrate into the telomeric repeat in a sequence-specific manner, we expressed the apurinic-apryrimidinic endonuclease-like endonuclease domain of TRAS1 (TRAS1 EN), which is supposed to digest the target DNA, and characterized its enzymatic properties. Purified TRAS1 EN could generate specific nicks on both strands of the telomeric repeat sequence between T and A of the (TTAGG)(n) strand (bottom strand) and between C and T of the (CCTAA)(n) strand (top strand). These sites are consistent with insertion sites expected from the genomic structure of boundary regions of TRAS1. Time course studies of nicking activities on both strands revealed that the cleavages on the bottom strand preceded those on the top strand, supporting the target-primed reverse transcription model. TRAS1 EN could cleave the telomeric repeats specifically even if it was flanked by longer tracts of nontelomeric sequence, indicating that the target site specificity of the TRAS1 element was mainly determined by its EN domain. Based on mutation analyses, TRAS1 EN recognizes less than 10 bp around the initial cleavage site (upstream 7 bp and downstream 3 bp), and the GTTAG sequence especially is essential for the cleavage reaction on the bottom strand (5'. TTAGGTT downward arrow AGG. 3'). TRAS1 EN, the first identified endonuclease digesting telomeric repeats, may be used as a genetic tool to shorten the telomere in insects and some other organisms.

Backbone dynamics of sequence specific recognition and binding by the yeast Pho4 bHLH domain probed by NMR

Backbone dynamics of the basic/helix-loop-helix domain of Pho4 from Saccharomyces cerevisae have been probed by NMR techniques, in the absence of DNA, nonspecifically bound to DNA and bound to cognate DNA. Alpha proton chemical shift indices and nuclear Overhauser effect patterns were used to elucidate the secondary structure in these states. These secondary structures are compared to the co-crystal complex of Pho4 bound to a cognate DNA sequence (Shimizu T. Toumoto A, Ihara K, Shimizu M, Kyogou Y, Ogawa N, Oshima Y, Hakoshima T, 1997, EMBO J 15: 4689-4697). The dynamic information provides insight into the nature of this DNA binding domain as it progresses from free in solution to a specifically bound DNA complex. Relative to the unbound form, we show that formation of either the nonspecific and cognate DNA bound complexes involves a large change in conformation and backbone dynamics of the basic region. The nonspecific and cognate complexes, however, have nearly identical secondary structure and backbone dynamics. We also present evidence for conformational flexibility at a highly conserved glutamate basic region residue. These results are discussed in relation to the mechanism of sequence specific recognition and binding.

One- and three-dimensional pathways for proteins to reach specific DNA sites

Proteins that interact with specific DNA sites bind to DNA at random and then translocate to the target site. This may occur by one-dimensional diffusion along the DNA, or through three-dimensional space via multiple dissociation/re-associations. To distinguish these routes, reactions of the ECO:RV endonuclease were studied on substrates with two ECO:RV sites separated by varied distances. The fraction of encounters between the DNA and the protein that resulted in the processive cleavage of both sites decreased as the length of intervening DNA was increased, but not in the manner demanded for one-dimensional diffusion. The variation in processivity with inter-site spacing shows instead that protein moves from one site to another through three-dimensional space, by successive dissociation/re-associations, though each re-association to a new site is followed by a search of the DNA immediately adjacent to that site. Although DNA-binding proteins are usually thought to find their target sites by one-dimensional pathways, three-dimensional routes may be more common than previously anticipated.

SURVEY AND SUMMARY: holliday junction resolvases and related nucleases: identification of new families, phyletic distribution and evolutionary trajectories

Holliday junction resolvases (HJRs) are key enzymes of DNA recombination. A detailed computer analysis of the structural and evolutionary relationships of HJRs and related nucleases suggests that the HJR function has evolved independently from at least four distinct structural folds, namely RNase H, endonuclease, endonuclease VII-colicin E and RusA. The endonuclease fold, whose structural prototypes are the phage lambda exonuclease, the very short patch repair nuclease (Vsr) and type II restriction enzymes, is shown to encompass by far a greater diversity of nucleases than previously suspected. This fold unifies archaeal HJRs, repair nucleases such as RecB and Vsr, restriction enzymes and a variety of predicted nucleases whose specific activities remain to be determined. Within the RNase H fold a new family of predicted HJRs, which is nearly ubiquitous in bacteria, was discovered, in addition to the previously characterized RuvC family. The proteins of this family, typified by Escherichia coli YqgF, are likely to function as an alternative to RuvC in most bacteria, but could be the principal HJRs in low-GC Gram-positive bacteria and AQUIFEX: Endonuclease VII of phage T4 is shown to serve as a structural template for many nucleases, including MCR:A and other type II restriction enzymes. Together with colicin E7, endonuclease VII defines a distinct metal-dependent nuclease fold. As a result of this analysis, the principal HJRs are now known or confidently predicted for all bacteria and archaea whose genomes have been completely sequenced, with many species encoding multiple potential HJRs. Horizontal gene transfer, lineage-specific gene loss and gene family expansion, and non-orthologous gene displacement seem to have been major forces in the evolution of HJRs and related nucleases. A remarkable case of displacement is seen in the Lyme disease spirochete Borrelia burgdorferi, which does not possess any of the typical HJRs, but instead encodes, in its chromosome and each of the linear plasmids, members of the lambda exonuclease family predicted to function as HJRs. The diversity of HJRs and related nucleases in bacteria and archaea contrasts with their near absence in eukaryotes. The few detected eukaryotic representatives of the endonuclease fold and the RNase H fold have probably been acquired from bacteria via horizontal gene transfer. The identity of the principal HJR(s) involved in recombination in eukaryotes remains uncertain; this function could be performed by topoisomerase IB or by a novel, so far undetected, class of enzymes. Likely HJRs and related nucleases were identified in the genomes of numerous bacterial and eukaryotic DNA viruses. Gene flow between viral and cellular genomes has probably played a major role in the evolution of this class of enzymes. This analysis resulted in the prediction of numerous previously unnoticed nucleases, some of which are likely to be new restriction enzymes.

Recognition of native DNA methylation by the PvuII restriction endonuclease

Recognizing the methylation status of specific DNA sequences is central to the function of many systems in eukaryotes and prokaryotes. Restriction-modification systems have to distinguish between 'self' and 'non-self' DNA and depend on the inability of restriction endonucleases to cleave their DNA substrates when the DNA is appropriately methylated. These endonucleases thus provide a model system for studying the recognition of DNA methylation by proteins. We have characterized the interaction of R.PVU:II with DNA containing the physiologically relevant N4-methylcytosine modification. R.PVU:II binds (N4m)C-modified DNA and cleaves it very slowly. Methylated strands in hemimethylated duplexes were cleaved at a higher rate than in fully methylated duplexes, in parallel with a higher binding affinity for hemimethylated DNA. The co-crystal structures of R.PVU:II-DNA, together with a mutagenesis study, have implicated specific amino acids in recognition of the methylatable base; one of these is His84. We report that replacing His84 with Ala reduced the rate of cleavage of unmodified DNA but, in contrast, slightly increased the cleavage of (N4m)C-modified DNA.

Crystal structure of NaeI-an evolutionary bridge between DNA endonuclease and topoisomerase

NAE:I is transformed from DNA endonuclease to DNA topoisomerase and recombinase by a single amino acid substitution. The crystal structure of NAE:I was solved at 2.3 A resolution and shows that NAE:I is a dimeric molecule with two domains per monomer. Each domain contains one potential DNA recognition motif corresponding to either endonuclease or topoisomerase activity. The N-terminal domain core folds like the other type II restriction endonucleases as well as lambda-exonuclease and the DNA repair enzymes MutH and Vsr, implying a common evolutionary origin and catalytic mechanism. The C-terminal domain contains a catabolite activator protein (CAP) motif present in many DNA-binding proteins, including the type IA and type II topoisomerases. Thus, the NAE:I structure implies that DNA processing enzymes evolved from a few common ancestors. NAE:I may be an evolutionary bridge between endonuclease and DNA processing enzymes.

Crystallographic snapshots along a protein-induced DNA-bending pathway

Two new high-resolution cocrystal structures of EcoRV endonuclease bound to DNA show that a large variation in DNA-bending angles is sampled in the ground state binary complex. Together with previous structures, these data reveal a contiguous series of protein conformational states delineating a specific trajectory for the induced-fit pathway. Rotation of the DNA-binding domains, together with movements of two symmetry-related helices binding in the minor groove, causes base unstacking at a key base-pair step and propagates structural changes that assemble the active sites. These structures suggest a complex mechanism for DNA bending that depends on forces generated by interacting protein segments, and on selective neutralization of phosphate charges along the inner face of the bent double helix.