Specific endonuclease BbvBI from Bacillus brevis]

A new restriction endonuclease BbvBI free from contaminating nonspecific nucleases and phosphatases was isolated from the Bacillus brevis cells. The enzyme was purified by fractionating the sonicated cell-free extract in a two-phase PEG/dextran system and subsequent chromatographies on DEAE-sepharose, blue sepharose and heparin sepharose. The endonuclease BbvBI displayed the maximal activity at 45 degrees C, pH between 8.0 and 8.5, MgCl2 concentration in the range of 5-10 mM and at the low ionic strength. It is shown that the enzyme cleaves the sequence G'GYPC'C, with the preferential cleavage of GGTACC and GGCACC sites as compared with GGTGCC and GGCGCC. Thus, the restriction endonuclease BbvBI is a true isoschizomer of nuclease BanI.     

Effects of functional group changes in the EcoRV recognition site on the cleavage reaction catalyzed by the endonuclease

Oligodeoxynucleotides have been prepared which contain changes in the functional group pattern present in the EcoRV recognition site d(GATATC). These modifications involve the deletion of specific functional groups or the reversal of the relative positions of functional groups within the canonical six base pair recognition site. The duplex stability of these modified oligodeoxynucleotides has been assessed by determining the thermodynamic parameters characterizing helix formation. Steady-state kinetic parameters have been used to characterize the interaction of the modified oligodeoxynucleotides with the EcoRV endonuclease. The enzyme is very sensitive to the deletion of either of the adenine amino or thymine methyl groups, or the reversal of the relative positions of the adenine amino group and thymine carboxy group which form an interstrand hydrogen bond in the major groove of the B-DNA helix. Conversely, deletion of the guanine amino group had only minimal effects upon the measured kinetic parameters. Deletion of the exocyclic amino group from the "inner" dA-dT base pair resulted in the fragment which interacted with the enzyme on the basis of observed inhibition experiments but was not cleaved. The results suggest that the endonuclease interacts with its recognition sequence via contacts in the major groove of the B-DNA helix and that both hydrogen bonding to the adenine amino groups and also hydrophobic interactions with the thymine methyl groups are involved.     

Molecular mechanics of the interactions of spermine with DNA: DNA bending as a result of ligand binding.

We used energy minimization of a molecular mechanical force field to evaluate spermine interactions with B-form DNA oligomers with either alternating purine/pyrimidine or homopolymeric sequences. Four different positions for spermine docking--within, along, and bridging the minor groove and bridging the major groove--were assessed for each sequence. Interaction at the major groove of alternating purine/pyrimidine sequences appears to be the most favorable of all models assessed, and are associated with significant bending of DNA. Interactions at the major groove of homopolymers were less favorable than those of heteropolymers and showed little or no bending. Interactions with the minor groove were most favorable for spermine positioned near the base of the groove, and became less favorable as spermine was moved toward the top of the groove. Association along the phosphate backbone alone was the least favorable of the interactions.     

A comparison of the structural requirements for DNA cleavage by the isoschizomers HaeIII, BspRI and BsuRI

We have investigated the structural requirements for DNA cleavage by the isoschizomers HaeIII, BspRI and BsuRI which recognize the sequence -d(GGCC)-. For this purpose decadeoxynucleotides were synthesized by the solid-phase phosphotriester method and purified by high-performance liquid chromatography. The kinetics of cleavage of these oligodeoxynucleotides were determined for the three isoschizomers with the following results. The sequence adjacent to the recognition site strongly influences the rate of cleavage. The preference is qualitatively the same for all three enzymes: AGGCCT greater than TGGCCA greater than GGGCCC approximately equal to CGGCCG, and follows the thermal stability of the different decanucleotides. Substitutions within the recognition site, namely dI for dG and dU for dC, affect the rate of cleavage differently for the three enzymes. The results can be rationalized in terms of an interaction of HaeIII with the major and minor groove of the DNA, of BspRI mainly with the minor groove and of BsuRI with the major groove of DNA. It is obvious from our data that the mechanism of recognition of the same site is different for the three isoschizomers.     

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

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

Molecular dynamics of spermine-DNA interactions: sequence specificity and DNA bending for a simple ligand.

We used molecular dynamics to model interactions between the physiologically important polyamine spermine and two B-DNA oligomers, the homopolymer (dG)10-(dC)10 and the heteropolymer (dGdC)5-(dGdC)5. Water and counterions were included in the simulation. Starting coordinates for spermine-DNA complexes were structures obtained by molecular mechanics modeling of spermine with the two oligomers; in these models, spermine binding induced a bend in the heteropolymer but not in the homopolymer. During approximately 40 psec of molecular dynamics simulation, spermine moves away from the floor of the major groove and interacts nospecifically with d(G)10-d(C)10. In contrast, a spermine-induced bend in the helix of (dGdC)5-(dGdC)5 is maintained throughout the simulation and spermine remains closely associated with the major groove. These results provide further evidence that the binding of spermine to nucleic acids can be sequence specific and that bending of alternating purine-pyrimidine sequences may be a physiologically important result of spermine binding.     

DNA curvature in native and modified EcoRI recognition sites and possible influence upon the endonuclease cleavage reaction

The ligation of a decadeoxynucleotide containing the EcoRI recognition site forms a series of multimers which appear to be curved based on observed anomalous gel migration in polyacrylamide gels. The degree of DNA curvature present in the recognition sequence, based upon the observed migration anomaly, can be altered by modifications to the purine functional groups at the 2- and 6-positions. Deletion of the guanine 2-amino group, occurring in the minor groove of the B-DNA helix, is most effective in increasing the observed DNA curvature. Conversely, the displacement of an amino group from the major groove to the minor groove eliminates curvature. DNA curvature is also modulated by the exocyclic group at the purine 6-position with decreasing curvature observed when changing the amino group to a carbonyl or proton substituent. Differences in the kinetic parameters characterizing the cleavage reaction by the endonuclease for many of the modified sequences are the result of modifications of functional groups in the major groove, which are likely to contact the endonuclease during catalysis. However, with two examples, significant decreases in the observed specificity constant (kcat/Km), characterizing the protein-nucleic acid interaction, cannot be easily explained in terms of such functional group contacts. It is more likely in these cases that the functional group modifications affect the efficiency of the endonuclease-DNA interaction by modulation of the structure of the double-stranded DNA helix. With both examples, modifications have been made to minor groove substituents. The extent of DNA curvature is increased significantly for one and decreased for the other, compared with that observed for the native recognition site. The results suggest that curvature of the DNA helix axis is an intrinsic property of the d(GAATTC) sequence which helps to optimize the protein-nucleic acid interactions observed for the EcoRI restriction endonuclease.     

Effects of functional group changes in the EcoRI recognition site on the cleavage reaction catalyzed by the endonuclease

Oligodeoxynucleotides have been prepared that contain changes in the functional group pattern present in the EcoRI recognition site. These changes involve "functional group deletions", "functional group reversals", and "displaced functional groups". Steady-state kinetic parameters have been used to characterize the interaction of these modified recognition sites with the EcoRI endonuclease. Changes in the functional group pattern have varying effects upon the cleavage reaction. Both the exocyclic amino groups of the two adenine residues and the methyl groups of the thymine residues appear to interact with the endonuclease quite differently. In both cases efficient catalysis was observed when these functional groups were present at the "outer" dA-dT base pair. Selectivity was decreased by over an order of magnitude largely via increases in Km when these functional groups were deleted. Similar modifications at the "inner" dA-dT base pair did not alter the kinetic parameters significantly from those observed with the native sequence. Addition of an amino group to the minor groove at the outer dA-dT base pair resulted in a modified recognition site that interacted with the enzyme, on the basis of observed competitive inhibition kinetics, but was not cleaved.     

DNA sequence-dependent positioning of the linker histone in a nucleosome: A single-pair FRET study

Linker histones (LHs) bind to nucleosomes with their globular domain (gH) positioned in either an on- or an off-dyad binding mode. Here, we study the effect of the linker DNA (L-DNA) sequence on the binding of a full-length LH, Xenopus laevis H1.0b, to a Widom 601 nucleosome core particle (NCP) flanked by two 40 bp long L-DNA arms, by single-pair FRET spectroscopy. We varied the sequence of the 11 bp of L-DNA adjoining the NCP on either side, making the sequence either A-tract, purely GC, or mixed with 64% AT. The labeled gH consistently exhibited higher FRET efficiency with the labeled L-DNA containing the A-tract than that with the pure-GC stretch, even when the stretches were swapped. However, it did not exhibit higher FRET efficiency with the L-DNA containing 64% AT-rich mixed DNA when compared to the pure-GC stretch. We explain our observations with a model that shows that the gH binds on dyad and that two arginines mediate recognition of the A-tract via its characteristically narrow minor groove. To investigate whether this on-dyad minor groove-based recognition was distinct from previously identified off-dyad major groove-based recognition, a nucleosome was designed with A-tracts on both the L-DNA arms. One A-tract was complementary to thymine and the other to deoxyuridine. The major groove of the thymine-tract was lined with methyl groups that were absent from the major groove of the deoxyuridine tract. The gH exhibited similar FRET for both these A-tracts, suggesting that it does not interact with the thymine methyl groups exposed on the major groove. Our observations thus complement previous studies that suggest that different LH isoforms may employ different ways of recognizing AT-rich DNA and A-tracts. This adaptability may enable the LH to universally compact scaffold-associated regions and constitutive heterochromatin, which are rich in such sequences.     

A conserved major groove antideterminant for Saccharomyces cerevisiae RNase III recognition

Rnt1p, the only known Saccharomyces cerevisiae RNase III double-stranded RNA endonuclease, plays important roles in the processing of precursors of ribosomal RNAs and small nuclear and nucleolar RNAs and in the surveillance of unspliced pre-mRNAs. Specificity of cleavage by Rnt1p relies on the presence of RNA tetraloop structures with the consensus sequence AGNN at the top of the target dsRNA. The sequences of 79 fungal RNase III substrates were inspected to identify additional conserved sequence elements or antideterminants that may contribute to Rnt1p recognition of the double-stranded RNA. Surprisingly, U-A sequences at the base pair adjacent to the conserved terminal tetraloop (closing base pair) were found to be absent from all but one inspected sequence. Analysis of chemically modified variants of the closing base pair showed that the presence of exocyclic groups in the major groove of purines 3' to the last nucleotide of the tetraloop inhibits Rnt1p cleavage without strongly inhibiting Rnt1p binding. We propose that these groups interfere with the recognition of the RNA substrate by the catalytic domain of Rnt1p. These results identify exocyclic groups of purines in the major groove downstream of the tetraloop as a major antideterminant in S. cerevisiae RNase III activity, and suggest a rationale for their apparent counter selection in RNA processing sites.     

Oligonucleotide duplexes containing inosine, 7-deazainosine, tubercidin, nebularine and 7-deazanebularine as substrates for restriction endonucleases HindII, SalI and TaqI.

Synthetic hexadecanucleotide duplexes containing a single purine nucleotide analogue in the recognition sites of the restriction endonucleases HindII, SalI and TaqI were used to investigate the restriction site determinants required by these enzymes for sequence recognition and phosphodiester bond cleavage. The enzymes were, in general, unaffected by changes introduced into the minor groove of the helix. SalI was found to be inhibited by the major groove modifications introduced into the fourth position of its recognition sequence GTCGAC. HindII and TaqI were, by contrast, able to cleave the sites containing the analogues at this position. TaqI and, to a lesser extent, HindII could also be shown to tolerate "mismatch analogues" at this site.     

Ordered water structure around a B-DNA dodecamer: A quantitative study

The crystal structure of the double-helical B-DNA dodecamer of sequence C-G-C-G-A-A-T-T-C-G-C-G has been solved and refined independently in three forms: (1) the parent sequence at room temperature; (2) the same sequence at 16 K; and (3) the 9-bromo variant C-G-C-G-A-A-T-T^BrC-G-C-G at 7 °C in 60% (v/v) 2-methyl-2.4-pentanediol. The latter two structures show extensive hydration along the phosphate backbone, a feature that was invisible in the native structure because of high temperature factors (indicating thermal or static disorder) of the backbone atoms. Sixty-five solvent peaks are associated with the phosphate backbone, or an average of three per phosphate group. Nineteen other molecules form a first shell of hydration to base edge N and O atoms within the major groove, and 36 more are found in upper hydration layers. The latter tend to occur in strings or clusters spanning the major groove from one phosphate group to another. A single spermine molecule also spans the major groove. In the minor groove, the zig-zag spine of hydration that we believe to be principally responsible for stabilizing the B form of DNA is found in all three structures. Upper level hydration in the minor groove is relatively sparse, and consists mainly of strings of water molecules extending across the groove, with few contacts to the spine below. Sugar O-1′ atoms are closely associated with water molecules, but these are chiefly molecules in the spine, so the association may reflect the geometry of the minor groove rather than any intrinsic attraction of O-1′ atoms for hydration. The phosphate O-3′ and O-5′ atoms within the backbone chain are least hydrated of all, although no physical or steric impediment seems to exist that would deny access to these oxygen atoms by water molecules.     

Conformational characteristics and correlations in crystal structures of nucleic acid oligonucleotides: evidence for sub-states

Sugar phosphate backbone conformations are a structural element inextricably involved in a complete understanding of specific recognition nucleic acid ligand interactions, from early stage discrimination of the correct target to complexation per se, including any structural adaptation on binding. The collective results of high resolution DNA, RNA and protein/DNA crystal structures provide an opportunity for an improved and enhanced statistical analysis of standard and unusual sugar-phosphate backbone conformations together with corresponding dinucleotide sequence effects as a basis for further exploration of conformational effects on binding. In this study, we have analyzed the conformations of all relevant crystal structures in the nucleic acids data base, determined the frequency distribution of all possible epsilon, zeta, alpha, beta and gamma backbone angle arrangements within four nucleic acid categories (A-RNA and A-DNA, free and bound B-DNA) and explored the relationships between backbone angles, sugar puckers and selected helical parameters. The trends in the correlations are found to be similar regardless of the nucleic acid category. It is interesting that specific structural effects exhibited by the different unusual backbone sub-states are in some cases contravariant. Certain alpha/gamma changes are accompanied by C3' endo (north) sugars, small twist angles and positive values of base pair roll, and favor a displacement of nucleotide bases towards the minor groove compared to that of canonical B form structures. Unusual epsilon/zeta combinations occur with C2' (south) sugars, high twist angles, negative values of base pair roll, and base displacements towards the major groove. Furthermore, any unusual backbone correlates with a reduced dispersion of equilibrium structural parameters of the whole double helix, as evidenced by the reduced standard deviations of almost all conformational parameters. Finally, a strong sequence effect is displayed in the free oligomers, but reduced somewhat in the ligand bound forms. The most variable steps are GpA and CpA, and, to a lesser extent, their partners TpC and TpG. The results provide a basis for considering if the variable and non-variable steps within a biological active sequence precisely determine morphological structural features as the curvature direction, the groove depth, and the accessibility of base pair for non covalent associations.     

O6-methylguanine in place of guanine causes asymmetric single-strand cleavage of DNA by some restriction enzymes

The interactions of restriction enzymes with their cognate DNA recognition sequences present a model for protein-DNA interactions. We have investigated the effect of O6-methylguanine on restriction enzyme cleavage of DNA; O6-methylguanine is a carcinogenic lesion and a structural analogue of the biological restriction inhibitor N6-methyladenine. O6-Methylguanine was synthesized into oligonucleotides at unique positions. The oligonucleotides were purified and analyzed by high-pressure liquid chromatography to assure that, within the limits of our detection, O6-methylguanine was the only modified base present. These oligonucleotides were annealed with their complement so that cytosine, and in one case thymine, opposed O6-methylguanine. DNA cleavage by restriction enzymes that recognize a unique DNA sequence, HpaII, HhaI, HinPI, NaeI, NarI, PvuII, and XhoI, was inhibited by a single O6-methylguanine in place of guanine (adenine for PvuII) within the appropriate recognition sequences. However, only the modified strand was nicked by HpaII, NaeI, and XhoI with O6-methylguanine at certain positions, indicating asymmetric strand cleavage. For all the restriction enzymes studied but AhaII, BanI, and NarI, lack of double- or single-strand cleavage correlated with inability of the O6-methylguanine-containing recognition sequence to measurably bind enzyme. None of the restriction enzymes studied were inhibited by O6-methylguanine outside their cognate recognition sequences.     

Synthesis of fluorescent microgonotropens (FMGTs) and their interactions with dsDNA

A new class of microgonotropen compounds (FIMGTs), which fluoresce upon binding to dsDNA, is introduced. The FMGTs consist of a minor groove binding moiety based upon Hoescht 33258 covalently attached to a polyamine chain capable of interacting with the phosphodiester backbone of dsDNA. The interactions of FMGTs with dsDNA were investigated by fluorescence and UV spectroscopy. Several different dsDNA oligomers were studied to determine the effect of binding site sequence on stoichiometric and binding affinity. The FMGTs were found to bind a dsDNA oligomer that contained the sequence 5'-AATTT-3' with FMGT:dsDNA stoichiometrics equal to 2:1 or 3:1. Hoechst 33258 bound the same dsDNA oligomer with a 1:1 stoichiometry. The second and third order equilibrium constants for complexation were determined to be Log(K1K2) = 17.9 M(-2) and Log(K1K2K3) = 26.1 M(-3), respectively, for two of strongest binding FMGTs. From thermal melting experiments deltaTm for Hoechst 33258 was determined to be 10 degrees C while the deltaTm values for FMGTs ranged from 20-26 degrees C indicating the greater stability of the latter.     

Insight Derived from Molecular Dynamics Simulation into Substrate-Induced Changes in Protein Motions of Proteinase K

Abstract

Because of the significant industrial, agricultural and biotechnological importance of serine protease proteinase K, it has been extensively investigated using experimental approaches such as X-ray crystallography, site-directed mutagenesis and kinetic measurement. However, detailed aspects of enzymatic mechanism such as substrate binding, release and relevant regulation remain unstudied. Molecular dynamics (MD) simulations of the proteinase K alone and in complex with the peptide substrate AAPA were performed to investigate the effect of substrate binding on the dynamics/molecular motions of proteinase K. The results indicate that during simulations the substrate-complexed proteinase K adopt a more compact and stable conformation than the substrate-free form. Further essential dynamics (ED) analysis reveals that the major internal motions are confined within a subspace of very small dimension. Upon substrate binding, the overall flexibility of the protease is reduced; and the noticeable displacements are observed not only in substrate-binding regions but also in regions opposite the substrate-binding groove/pockets. The dynamic pockets caused by the large concerted motions are proposed to be linked to the substrate recognition, binding, orientation and product release; and the significant displacements in regions opposite the binding groove/pockets are considered to play a role in modulating the dynamics of enzyme-substrate interaction. Our simulation results complement the biochemical and structural studies, highlighting the dynamic mechanism of the functional properties of proteinase K.