Type II restriction endonuclease R.Eco29kI is a member of the GIY-YIG nuclease superfamily
© Ibryashkina et al; licensee BioMed Central Ltd. 2007
Received: 12 February 2007
Accepted: 12 July 2007
Published: 12 July 2007
The majority of experimentally determined crystal structures of Type II restriction endonucleases (REases) exhibit a common PD-(D/E)XK fold. Crystal structures have been also determined for single representatives of two other folds: PLD (R.BfiI) and half-pipe (R.PabI), and bioinformatics analyses supported by mutagenesis suggested that some REases belong to the HNH fold. Our previous bioinformatic analysis suggested that REase R.Eco29kI shares sequence similarities with one more unrelated nuclease superfamily, GIY-YIG, however so far no experimental data were available to support this prediction. The determination of a crystal structure of the GIY-YIG domain of homing endonuclease I-TevI provided a template for modeling of R.Eco29kI and prompted us to validate the model experimentally.
Using protein fold-recognition methods we generated a new alignment between R.Eco29kI and I-TevI, which suggested a reassignment of one of the putative catalytic residues. A theoretical model of R.Eco29kI was constructed to illustrate its predicted three-dimensional fold and organization of the active site, comprising amino acid residues Y49, Y76, R104, H108, E142, and N154. A series of mutants was constructed to generate amino acid substitutions of selected residues (Y49A, R104A, H108F, E142A and N154L) and the mutant proteins were examined for their ability to bind the DNA containing the Eco29kI site 5'-CCGCGG-3' and to catalyze the cleavage reaction. Experimental data reveal that residues Y49, R104, E142, H108, and N154 are important for the nuclease activity of R.Eco29kI, while H108 and N154 are also important for specific DNA binding by this enzyme.
Substitutions of residues Y49, R104, H108, E142 and N154 predicted by the model to be a part of the active site lead to mutant proteins with strong defects in the REase activity. These results are in very good agreement with the structural model presented in this work and with our prediction that R.Eco29kI belongs to the GIY-YIG superfamily of nucleases. Our study provides the first experimental evidence for a Type IIP REase that does not belong to the PD-(D/E)XK or HNH superfamilies of nucleases, and is instead a member of the unrelated GIY-YIG superfamily.
Type II restriction endonucleases (REases) are one of the largest groups of biochemically characterized enzymes . They usually recognize a short (4–8 bp) palindromic sequence of double-stranded DNA and in the presence of Mg2+ catalyze the hydrolysis of phosphodiester bonds at precise positions within or close to this sequence, leaving "blunt" ends or "sticky" (5' or 3') overhangs. REases that do not fit this definition or exhibit certain structural and functional peculiarities have been classified into several subtypes . To date (May 2007), crystal structures of 24 Type II REases have been solved (for the continuously updated list see the REBASE database ). All but one share a common catalytic domain and a weakly conserved bipartite catalytic motif PD-(D/E)XK, suggesting that they are evolutionarily related . The same core was also found in a number of other nucleases, now grouped together with the afore-mentioned REases into the "PD-(D/E)XK superfamily" . Nonetheless, hundreds of known REase sequences typically exhibit little sequence similarity to one another and to other proteins, and it is not necessarily obvious that the available crystal structures are representative of all REases.
Indeed, based on sequence analyses it was predicted that REases may belong not only to the PD-(D/E)XK superfamily, but also to completely unrelated superfamilies of nucleases: PLD , HNH [6, 7], or GIY-YIG . So far, biochemical and crystallographic analyses confirmed that the catalytic domain of R.BfiI is a member of the Nuc/PLD superfamily . More recently, another study first predicted that R.PabI may exhibit a different three-dimensional fold  and later confirmed it experimentally, revealing a new fold termed 'half-pipe' . Besides, biochemical analyses have supported the prediction that the catalytic domain of R.KpnI belongs to the HNH superfamily . Thus, there is experimental evidence that REases may belong to four different folds/superfamilies: PD-(D/E)XK, PLD, half-pipe, and HNH. However, compelling experimental support is still missing for the prediction that REases may also belong to the GIY-YIG superfamily .
Our initial bioinformatics analyses  revealed that REase R.Eco29kI (an enzyme from the Type II restriction-modification system Eco29kI) shares sequence similarities with the catalytic domain of GIY-YIG endonucleases, such as a homing endonuclease I-TevI. At that time, no crystal structure of any GIY-YIG superfamily member was available, therefore no three-dimensional model of the R.Eco29kI structure could be obtained to guide the analysis of sequence-structure-function relationships. We have carried out a Monte-Carlo simulation of the catalytic domain of the homing endonuclease I-TevI using spatial restraints derived from the published sparse NMR data  to generate the first, preliminary three-dimensional model of the GIY-YIG domain . Subsequent crystallographic resolution of the I-TevI catalytic domain structure confirmed that the computational prediction of the three-dimensional fold of the GIY-YIG domain was correct . Subsequently, structures of other GIY-YIG superfamily members were solved, including the 3' nuclease domain of UvrC by X-ray crystallography  the putative nuclease domain of a bacterial Slx-1 homolog EF2693 by NMR , and a comparative analysis of the GIY-YIG structures , providing robust templates for modeling other members of the superfamily. In this work, we present a revised sequence alignment between R.Eco29kI and the experimentally determined structure of a GIY-YIG domain, the first three-dimensional model of R.Eco29kI, and experimental validation of predicted active site residues. This is the first experimental study of R.Eco29kI that demonstrates importance of particular residues for its REase activity.
Molecular modeling of R.Eco29kI
DNA binding properties of R.Eco29kI and mutant proteins
Mutants with amino acid substitutions (Y49A, R104A, H108F, E142A and N154L) were constructed. In accordance with the structural model of R.Eco29kI, these mutations were to affect the residues localized in the active enzyme center. We also constructed a mutant protein with an alanine substitution of R86, i.e. a residue predicted to be located in the active site in our previous analysis  but regarded as not important functionally, according to the current model. Using a simple purification procedure under native conditions on a Ni-NTA column, the wt and mutant proteins were isolated to 90% purity. First of all, it was necessary to ascertain that the introduced mutations did not result in perturbation of protein conformation. Ability of the mutants to specifically bind DNA containing the Eco29kI site (i.e. the 5'-CCGCGG-3' sequence recognized specifically by the R.Eco29kI REase) would demonstrate that the mutations do not perturb the overall folding of the protein, and only change the local environment of the active site.
DNA cleavage and binding activities of wt R.Eco29kI and its mutants
DNA cleavage activity [% of wt]
standard bimolecular binding model
cooperative binding model
16 ± 5
299 ± 50
2.1 ± 0.2
49 ± 5
529 ± 50
1.9 ± 0.3
21 ± 5
344 ± 19
2.0 ± 0.4
43 ± 18
1911 ± 356
1.8 ± 0.2
114 ± 25
14285 ± 2000
2.0 ± 0.1
45 ± 11
970 ± 200
1.8 ± 0.4
50 ± 15
3100 ± 700
2.0 ± 0.3
DNA cleavage activity of wt R.Eco29kI and mutant proteins
The activity was measured from the intensity pattern of DNA cleavage products. The activity of the mutants was calculated as percentage of the wt R.Eco29kI cleavage activity taking into account the dilution of wt enzyme and reaction incubation time. The summarized data are presented in Table 1. To test mutant proteins for specific nicking activity the reaction products of cleavage of DNA substrate with a single Eco29kI site were analyzed under denaturing conditions (Figure 5(c), lanes 2–8). Mutants Y49A, H108F, and E142A had neither double-stranded DNA (ds DNA) cleavage nor nicking activity. Comparison of the profiles of R104A and N154L hydrolysis reaction products under native and denaturing conditions demonstrates a concerted weakening of ds DNA cleavage and nicking activities.
Most of Type II REases characterized to date were found to belong to the PD-(D/E)XK superfamily of nucleases , although theoretical predictions have been made that some REases belong to different, unrelated superfamilies. To date, predictions have been experimentally confirmed for REases from three non-PD-(D/E)XK superfamilies: R.BfiI belongs to the PLD superfamily , R.PabI exhibits a new 'half-pipe' fold  and R.KpnI and a few other enzymes belong to the HNH superfamily . Here, we present a refined alignment and the first structural model of the R.Eco29kI structure, and the first set of experimental data that support the prediction that REases may belong also to the GIY-YIG superfamily of nucleases.
According to the model, the active center of R.Eco29kI comprises amino acid residues Y49, Y76, R104, H108, E142, and N154. Alanine substitutions of the corresponding residues in I-TevI, the archetypal member of the GIY-YIG superfamily, have been shown to be critical for nuclease activity of the enzyme . The analysis of the crystal structure of the catalytic domain I-TevI confirmed the importance of these residues for the catalytic activity of the enzyme, having demonstrated that the putative catalytic residues are located on a shallow concave surface, and E75 (E142 of R.Eco29kI) was identified as a divalent cation-binding site . We constructed mutants of R.Eco29kI with amino acid substitutions Y49A, R104A, H108F, E142A and N154L, and R86A. Proteins with Ala and Phe substitutions of Y76 and H108A and N154A were not obtained. We did not find plasmids with required substitutions possibly because of toxicity of the resulting mutant proteins. We selected Phe as a replacement for His and Leu as a replacement for Asn because they are of nearly the same size and are chemically inert. The Ala substitution of R86 residue, predicted in this work to be functionally irrelevant, has no effect on the catalysis or binding as it should. The data on the specificity and mechanism of binding show that, although the introduced substitutions change the efficiency of binding of mutant protein forms, the spatial structure of constructed mutant forms of R.Eco29kI undergoes no dramatic changes.
Substitutions Y49A, H108F, and E142A completely eliminate the nuclease activity of R.Eco29kI just as I-TevI mutants Y6A and E75A. Further, substitutions R104A and N154L greatly reduce the R.Eco29kI catalytic activity, but do not abolish it completely. It is noteworthy that the residual nuclease activity was also found in the mutants of I-TevI, in which R27 and N90 (homolog's of R104 and N154 in R.EcoR29kI), were replaced by Ala [13, 14]. The efficiency of DNA-binding of Y49A, R104A and E142A mutants proved to be slightly reduced, whereas in N154L and in particular in H108F it is more significantly compromised. These data show that all the above-mentioned residues are to some extent involved in DNA-binding, but in the case of Y49, R104, and E142 the effect of mutations concerns mostly catalysis, while H108 and N154 are probably mostly important for DNA binding (as of today, it is difficult to determine the contribution of these residues to catalysis, the efficiency of which obviously depends on binding). It should be emphasized that residues R27 and E75 in I-TevI, homologous to residues R104 and E142 of R.Eco29kI, are of critical importance for catalysis . All these results are in very good agreement with the structural model presented in this work and with our prediction that R.Eco29kI, unlike all other REases studied to date, belongs to the GIY-YIG superfamily of nucleases . Importantly, modeling combined with mutagenesis allowed us to verify and refine the previous prediction of the active site in R.Eco29kI enzyme, which was based solely on the sequence considerations , demonstrating the utility of the methods for protein structure prediction. We would like to emphasize that our analysis relies primarily on the detection of similarities between R.Eco29kI and I-TevI, and does not represent an attempt to identify the functionally important residues of R.Eco29kI ab initio. There is no direct biochemical and biophysical data on the precise role of individual residues in the catalytic mechanism of DNA cleavage by I-TevI, therefore we refrain from making excessive speculations about the exact roles of the corresponding residues in R.EcoR29kI. On the other hand, any refinement of the I-TevI mechanism of action and the studies on the function of its residues may be reflected in appropriate adjustments of the interpretation of the R.Eco29kI model presented here.
It is noteworthy that at such a low level of similarity as between R.Eco29kI and I-TevI, one can expect some local structural divergence to occur between protein structures. In particular, secondary structure elements that are not directly hydrogen-bonded are prone to mutual translations or rotations induced by repacking of side chains in the protein core. Such shifts of homologous helices have been indeed observed between the remotely related GIY-YIG nucleases I-TevI and UvrC . Therefore, the root mean square deviation of the presented R.Eco29kI model from the true (currently unknown) structure of this protein can turn out to be substantial, in the range of similarity between I-TevI and UvrC, i.e. 2.9 Å for the 74 structurally superimposable core residues. Nonetheless, we are confident that the mutual orientation of major secondary structure elements that constitute the GIY-YIG fold, and the functionally important residues of R.Eco29kI studied in this work, is very accurate.
Extreme problems with the aggregation of R.Eco29kI at high concentrations (data not shown) hampered the attempts to obtain crystals or concentrated solutions suitable for NMR analyses, therefore the theoretical model of R.Eco29kI structure developed and validated in this work will serve as a starting point for computational docking and experimental analyses aimed at the identification of residues important for specific recognition of DNA by R.Eco29kI and understanding of the cleavage mechanism and architecture of the protein-DNA complex.
Extreme problems with the aggregation of R.Eco29kI at high concentrations preclude the experimental characterization of its structure by X-ray crystallography or NMR. Based on bioinformatics analyses, we predicted that R.EcoR29kI is a member of the GIY-YIG superfamily of nucleases and thus, it is unrelated to all other REases studied to date. A theoretical model of the R.Eco29kI structure was constructed using a combination of the protein fold-recognition approach to detect structural templates for homology modeling and de novo folding of insertions and terminal extensions unique to R.Eco29kI and not present in other protein structures. Substitutions of residues Y49, R104, H108, E142 and N154 predicted by the model to be a part of the active site lead to mutant proteins with strong defects in the catalytic activity. The mutants retain the ability to bind the substrate DNA, suggesting that the effect of substitutions is limited to the active site. The experimental data support the model based on the revised alignment that identified R104 as the catalytic residue and suggested R86 not to be important. Our study provides the first experimental evidence for a Type IIP REase that does not belong to the PD-(D/E)XK or HNH superfamilies of nucleases, and is instead a member of the unrelated GIY-YIG superfamily.
Protein structure prediction
Secondary structure prediction and fold-recognition analysis of R.Eco29kI were carried out via the GeneSilico metaserver gateway . Homology modeling of the catalytic core was carried out using the "FRankenstein's monster" approach (see  for a detailed description). Briefly, alternative sequence alignments between R.Eco29kI and template structures obtained from various FR servers with significant scores (all members of the GIY-YIG superfamily) were used to build preliminary models with MODELLER . The preliminary models were scored by VERIFY3D  and the best-scoring parts from all models were merged to form a "hybrid" model.
The "hybrid" model obtained was used as a starting point for folding simulations of the complete sequence using ROSETTA . The homology-modeled core of R.Eco29kI was completely "frozen" and the search of conformational space for the variable regions was restricted by the choice of fragments from known crystal structures that were compatible with the sequence and predicted secondary structure of R.Eco29kI. The final simulation of R.Eco29kI folding was conducted by the CABS method, which uses a reduced lattice representation of the protein chain . Simulations in CABS were guided by the long-range spatial restraints derived from the well-scored parts of the ROSETTA models and short-range restraints derived from the predicted secondary structure in the regions absent or poorly scored in the FR model. The "decoy" models generated during the simulation were clustered using the HCPM method , and the most representative member of the largest cluster (in terms of the root mean square deviation, RMSD) was selected as the final model.
Bacterial strains, plasmids and DNA
Wild-type eco29kIR gene was cloned in vector plasmid pQE-30 (Apr) (Qiagene, USA). The DNA of pQE30 plasmid with wt eco29kIR gene was used as a template for PCR amplification using mutagenic oligonucleotide primers listed below (mutant variants were generated using the QuickChange Site-Directed Mutagenesis Kit by Stratagene): Y49A: 5'-GGTGCTGGGG TGGCT GCTCT TTACT-3'; E142A: 5'-GATTAGTACA GTTGCG GCCG CTCTTA-3'; R104A: 5'-CTATCTAATGCA ATTAGAGA ACATGG-3'; H108F: 5'-CTAATAGAAT TAGAGAATTC GGCCGAAATA TAGC-3'; N154L: 5'-GCCTTTGTGG CTG ACCGTTG TTGATG-3'; N154A: 5'-GCCTTTGTGG GCT ACCGTTG TTGATG-3'; R86A: 5'-GGTTGGGCG C AGTCTAGAAT TAGTG-3'; Y76A: 5'-AACCTTCCTA TTGCT GTTGG CAAGGC-3'; Y76F: 5'-AACCTTCCTA TTTTC GTTGG CAAGGC-3'. Plasmids bearing the mutant variants of eco29kIR gene were sequenced to make sure that the substitution of the selected codon was introduced. Strain M15 [pREP4] (Qiagene, USA)? with plasmid p29Cm (Cmr) bearing eco29kIM gene  was used as a host for transformation of plasmids containing the wt and mutant eco29kIR genes.
Purification of wt and mutant R.Eco29kI
One-step purification of wt and mutant R.Eco29kI variants was performed under native conditions by affinity chromatography on Ni-NTA resin according to the handbook (Qiagen, USA). The protein was more than 90% homogeneous according to SDS-PAGE analysis. Protein concentration was determined by measuring OD280 on a Shimadzu UV-1601 spectrophotometer (Japan) with an extinction coefficient of 38120 M-1 cm-1 (the value was calculated using the PCGene software package).
Electrophoretic mobility shift assay (EMSA)
We studied the DNA binding activity of wt and mutant enzymes by EMSA with the PvuII-PstI fragments of pUC128 plasmid, containing a single Eco29kI site in the 150 bp DNA fragment, in the absence of Mg2+ ions. The binding reaction was performed in assay buffer AB (10 mM Tris-HCl, pH 7.5, 50 mM NaCl) with 250 nM of enzymes for 15 min at room temperature. Each sample was loaded onto 6% polyacrylamide gel. Products were analyzed under non-denaturing conditions.
Fluorescence anisotropy measurements
All fluorescence polarization measurements were performed using the thermo-jacketed Beacon 2000 Fluorescence Polarization System. As a substrate, we used an oligonucleotide duplex with the recognition sequence of R.Eco29kI conjugated to a 5(6)-carboxyfluorescein (FAM) via a six-carbon spacer at the 5' terminus (5'-FAM-TTTTGGTACCGCGG CCGCAAGCTT-3').
The 24 bp double-stranded oligonucleotide was prepared by annealing equimolar amounts of each strand. The DNA fragment (1 nM) was added to each diluted enzyme in AB buffer containing BSA (0.1 mg/mL) and the binding reaction mixture was incubated at room temperature for 20 min. Fluorescence polarization in each probe was measured in the static mode (measurement time 5 sec). The equilibrium binding data were analyzed using the non-linear regression according to two enzyme-DNA binding models: the standard bimolecular binding model (M1) and the cooperative binding model (M2):
E + S ↔ ES, (1)
nE + S ↔ EnS. (2)
In the cooperative binding model M2 , the apparent binding constant Kd app was calculated according to:mP = mP o + (mP max - mP 0 ) *E 0 n *(1/ K d app )/(1 + E 0 n *(1/ K d app )),
where mP – fluorescence polarization (FP) of the DNA-enzyme complex,
mP0 – minimum polarization (i.e., polarization of free DNA prior to the addition of protein),
mPmax – final polarization (i.e., polarization of DNA totally bound to protein),
S0 – initial concentration of the DNA duplex,
E0 – initial concentration of protein,
Kd – dissociation constant in M1,
n – magnitude of cooperativity,
Kdapp – apparent binding constant in M2.
Kd and Kdapp were calculated as described Baskunov et al.
DNA cleavage assay
The double-stranded (ds) DNA cleavage activity and nicking activities of wt R.Eco29kI and the mutants were determined. The ds DNA cleavage activity was assayed using the bacteriophage φ80vir DNA  and 200 bp DNA fragment with the Eco29kI site. The fragment was amplified using specific primers for pUC128: M13/pUC reverse sequencing primer (-26) and M13/pUC sequencing primer (+9). Only one strand of specific DNA was 5'-labeled with [γ-32 P] ATP using T4 polynucleotide kinase (Fermentas).
The reactions were carried out in 10 μl reaction buffer RB (10 mM Tris-HCl pH 7.5, 10 mM MgCl2, 50 mM NaCl) containing 1 μg of phage φ80vir DNA or 0.01 μM of DNA fragment. Mutant proteins and serially diluted wild type were incubated with DNA substrate at 37°C. The aliquots were removed after fixed time intervals until the reaction rate plateaued. The cleavage reaction products of φ80vir DNA were analyzed by 0.8% agarose gel electrophoresis. The cleavage reaction products of 200 bp DNA fragment were applied onto 5% polyacrylamide gel under non-denaturing conditions. Then the gel was autoradiographed.
The deoxyribonuclease activity was measured by quantifying the intensity of bands corresponding to DNA cleavage products in the PAGE analysis. Activity of mutants was calculated as percentage of the wt R.Eco29kI cleavage activity taking into account the dilution of wt enzyme and reaction incubation time . The DNA nicking activity of the proteins was checked by hydrolysis of the same 200 bp DNA fragment with the Eco29kI site. In this case, aliquots were applied onto 8% polyacrylamide gel containing 7M urea. After electrophoresis, the gel was soaked in 10% acetic acid, then was dried and autoradiographed.
We thank Dr. S. A. Eremin (Moscow University) for advice and help with Beacon 2000 Fluorescence Polarization System. The work at the laboratory of JMB was supported by the NIH (Fogarty International Center grant R03 TW007163-01). The work at the laboratory of AK was supported by KBN (grant PBZ-KBN-088/P04/2003). The work at IBPM RAS was supported by the Russian Foundation for Basic Research (Project 04-04-48784-a; M.Z.) and the NIH (Fogarty International Research Collaboration Award grant PA02-057).
- Pingoud AM: Restriction endonucleases. Springer-Verlag, Berlin, Heidelberg; 2004.View Article
- Roberts RJ, Vincze T, Posfai J, Macelis D: REBASE: restriction enzymes and methyltransferases. Nucleic Acids Res 2003, 31: 418–420. 10.1093/nar/gkg069PubMed CentralView ArticlePubMed
- Pingoud A, Fuxreiter M, Pingoud V, Wende W: Type II restriction endonucleases: structure and mechanism. Cell Mol Life Sci 2005, 62: 685–707. 10.1007/s00018-004-4513-1View ArticlePubMed
- Kosinski J, Feder M, Bujnicki JM: The PD-(D/E)XK superfamily revisited: identification of new members among proteins involved in DNA metabolism and functional predictions for domains of (hitherto) unknown function. BMC Bioinformatics 2005, 12(6):172. 10.1186/1471-2105-6-172View Article
- Sapranauskas R, Sasnauskas G, Lagunavicius A, Vilkaitis G, Lubys A, Siksnys V: Novel subtype of type IIs restriction enzymes. BfiI endonuclease exhibits similarities to the EDTAresistant nuclease Nuc of Salmonella typhimurium . J Biol Chem 2000, 275: 30878–30885. 10.1074/jbc.M003350200View ArticlePubMed
- Aravind L, Makarova KS, Koonin EV: SURVEY AND SUMMARY: holliday junction resolvases and related nucleases: identification of new families, phyletic distribution and evolutionary trajectories. Nucleic Acids Res 2000, 28: 3417–3432. 10.1093/nar/28.18.3417PubMed CentralView ArticlePubMed
- Bujnicki JM, Radlinska M, Rychlewski L: Polyphyletic evolution of type II restriction enzymes revisited: two independent sources of second-hand folds revealed. Trends Biochem Sci 2001, 26: 9–11. 10.1016/S0968-0004(00)01690-XView ArticlePubMed
- Grazulis S, Manakova E, Roessle M, Bochtler M, Tamulaitiene G, Huber R, Siksnys V: Structure of the metal-independent restriction enzyme BfiI reveals fusion of a specific DNA-binding domain with a nonspecific nuclease. Proc Natl Acad Sci USA 2005, 102: 15797–15802. 10.1073/pnas.0507949102PubMed CentralView ArticlePubMed
- Ishikawa K, Watanabe M, Kuroita T, Uchiyama I, Bujnicki JM, Kawakami B, Tanokura M, Kobayashi I: Discovery of a novel restriction endonuclease by genome comparison and application of a wheat-germ-based cell-free translation assay: PabI (5'GTA/C) from the hyperthermophilic archaeon Pyrococcus abyssi. Nucleic Acids Res 33(13):e112. 2005, Jul 21 10.1093/nar/gni113
- Miyazono K, Watanabe M, Kosinski J, Ishikawa K, Kamo M, Sawasaki T, Nagata K, Bujnicki JM, Endo Y, Tanokura M, Kobayashi I: Novel protein fold discovered in the PabI family of restriction enzymes. Nucleic Acids Res 2007, 35(6):1908–1918. Epub 2007 Mar 1 10.1093/nar/gkm091PubMed CentralView ArticlePubMed
- Saravanan M, Bujnicki JM, Cymerman IA, Rao DN, Nagaraja V: Type II restriction endonuclease R.KpnI is a member of the HNH nuclease superfamily. Nucleic Acids Res 2004, 32: 6129–6135. 10.1093/nar/gkh951PubMed CentralView ArticlePubMed
- Dunin-Horkawicz S, Feder M, Bujnicki JM: Phylogenomic analysis of the GIY-YIG nuclease superfamily. BMC Genomics 7(1):98. 2006, Apr 28 10.1186/1471-2164-7-98PubMed CentralView ArticlePubMed
- Kowalski JC, Belfort M, Stapleton MA, Holpert M, Dansereau JT, Pietrokovski S, Baxter SM, Derbyshire V: Configuration of the catalytic GIY-YIG domain of intron endonuclease I-TevI: coincidence of computational and molecular findings. Nucleic Acids Res 1999, 27: 2115–2125. 10.1093/nar/27.10.2115PubMed CentralView ArticlePubMed
- Van Roey P, Meehan L, Kowalski JC, Belfort M, Derbyshire V: Catalytic domain structure and hypothesis for function of GIY-YIG intron endonuclease I-TevI. Nat Struct Biol 2002, 9: 806–811.PubMed
- Truglio JJ, Rhau B, Croteau DL, Wang L, Skorvaga M, Karakas E, Dellavecchia MJ, Wang H, Van Houten B, Kisker C: Structural insights into the first incision reaction during nucleotide excision repair. EMBO J 2005, 24: 885–894. 10.1038/sj.emboj.7600568PubMed CentralView ArticlePubMed
- Snyder DA, Chen Y, Denissova NG, Acton T, Aramini JM, Ciano M, Karlin R, Liu J, Manor P, Rajan PA, Rossi P, Swapna GV, Xiao R, Rost B, Hunt J, Montelione GT: Comparisons of NMR spectral quality and success in crystallization demonstrate that NMR and X-ray crystallography are complementary methods for small protein structure determination. J Am Chem Soc 2005, 127(47):16505–16511. 10.1021/ja053564hView ArticlePubMed
- Bujnicki JM, Rotkiewicz P, Kolinski A, Rychlewski L: Three-dimensional modeling of the I-TevI homing endonuclease catalytic domain, a GIY-YIG superfamily member, using NMR restraints and Monte Carlo dynamics. Protein Eng 2001, 14: 717–721. 10.1093/protein/14.10.717View ArticlePubMed
- Kosinski J, Cymerman IA, Feder M, Kurowski MA, Sasin JM, Bujnicki JM: A "FRankenstein's monster" approach to comparative modeling: merging the finest fragments of Fold-Recognition models and iterative model refinement aided by 3D structure evaluation. Proteins 2003, 53(Suppl 6):369–379. 10.1002/prot.10545View ArticlePubMed
- Kolinski A, Bujnicki JM: Generalized protein structure prediction based on combination of fold-recognition with de novo folding and evaluation of models. Proteins 2005, 61(Suppl 7):84–90. 10.1002/prot.20723View ArticlePubMed
- Kosinski J, Gajda MJ, Cymerman IA, Kurowski MA, Pawlowski M, Boniecki M, Obarska A, Papaj G, Sroczynska-Obuchowicz P, Tkaczuk KL, Sniezynska P, Sasin JM, Augustyn A, Bujnicki JM, Feder M: FRankenstein becomes a cyborg: the automatic recombination and realignment of Fold-Recognition models in CASP6. Proteins 2005, 61(Suppl 7):106–113. 10.1002/prot.20726View ArticlePubMed
- GeneSilico model repository[ftp://genesilico.pl/iamb/models/]
- Pingoud A, Jeltsch A: Recognition and cleavage of DNA by type-II restriction endonucleases. Eur J Biochem 1997, 246: 1–22. 10.1111/j.1432-1033.1997.t01-6-00001.xView ArticlePubMed
- Bujnicki JM: Molecular phylogenetics of restriction endonucleases. In Restriction Endonucleases: V 14. Edited by: Pingoud A. Springer-Verlag, Berlin, Heidelberg; 2004:63–93.View Article
- Cymerman IA, Obarska A, Skowronek KJ, Lubys A, Bujnicki JM: Identification of a new subfamily of HNH nucleases and experimental characterization of a representative member, HphI restriction endonuclease. Proteins 65(4):867–876. 2006, Dec 1 10.1002/prot.21156View ArticlePubMed
- Kurowski MA, Bujnicki JM: GeneSilico protein structure prediction meta-server. Nucleic Acids Res 2003, 31: 3305–3307. 10.1093/nar/gkg557PubMed CentralView ArticlePubMed
- Fiser A, Sali A: Modeller: generation and refinement of homology-based protein structure models. Methods Enzymol 2003, 374: 461–491.View ArticlePubMed
- Luthy R, Bowie JU, Eisenberg D: Assessment of protein models with three-dimensional profiles. Nature 1992, 356: 83–85. 10.1038/356083a0View ArticlePubMed
- Simons KT, Kooperberg C, Huang E, Baker D: Assembly of protein tertiary structures from fragments with similar local sequences using simulated annealing and Bayesian scoring functions. J Mol Biol 1997, 268: 209–225. 10.1006/jmbi.1997.0959View ArticlePubMed
- Kolinski A: Protein modeling and structure prediction with a reduced representation. Acta Biochim Pol 2004, 51: 349–371.PubMed
- Gront D, Kolinski A: HCPM--program for hierarchical clustering of protein models. Bioinformatics 2005, 21: 3179–3180. 10.1093/bioinformatics/bti450View ArticlePubMed
- Zakharova MV, Beletskaya IV, Bolovin DV, Yurkova TV, Semenova LM, Solonin AS: Structural plasmid evolution as a result of coupled recombinations at bom and cer sites. Mol Genet Genomics 2003, 270: 415–419. 10.1007/s00438-003-0933-0View ArticlePubMed
- Reid SL, Parry D, Liu H-H, Connolly BA: Binding and recognition of GATATC target sequences by the EcoRV restriction endonuclease: a study using fluorescent oligonucleotides and fluorescence polarization. Biochemistry 2001, 40: 2484–2494. 10.1021/bi001956pView ArticlePubMed
- Baskunov VB, Subach FV, Kolbanovskiy A, Kolbanovskiy M, Eremin SA, Johnson F, Bonala R, Geacintov NE, Gromova ES: Effects of benzo[a]pyrene-deoxyguanosine lesions on DNA methylation catalyzed by EcoRII DNA methyltransferase and on DNA cleavage effected by EcoRII restriction endonuclease. Biochemistry 2005, 44: 1054–1066. 10.1021/bi048130yView ArticlePubMed
- Ackers GK, Johnson AD, Shea MA: Quantitative model for gene regulation by lambda phage repressor. Proc Natl Acad Sci U S A 1982, 79: 1129–1133. 10.1073/pnas.79.4.1129PubMed CentralView ArticlePubMed
- Miller JH: Experiments in Molecular Genetics. Cold Spring Harbor Laboratory; 1972.
- Holtz JK, Topal MD: Location of putative binding and catalytic sites of NaeI by random mutagenesis. J Biol Chem 1994, 269(44):27286–27290.PubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.