Cloning of NruI and Sbo13I restriction and modification sstems in E. coli and amino acid sequence comparison of M.NruI and M.Sbo13I with other amino-methyltransferases
© Benner et al; licensee BioMed Central Ltd. 2010
Received: 23 April 2010
Accepted: 24 May 2010
Published: 24 May 2010
NruI and Sbo13I are restriction enzyme isoschizomers with the same recognition sequence 5' TCG↓CGA 3' (cleavage as indicated↓). Here we report the cloning of NruI and Sbo13I restriction-modification (R-M) systems in E. coli. The NruI restriction endonuclease gene (nruIR) was cloned by PCR and inverse PCR using primers designed from the N-terminal amino acid sequence. The NruI methylase gene (nruIM) was derived by inverse PCR walking.
The amino acid sequences of NruI endonuclease and methylase are very similar to the Sbo13I R-M system which has been cloned and expressed in E. coli by phage selection of a plasmid DNA library. Dot blot analysis using rabbit polyclonal antibodies to N6mA- or N4mC-modified DNA indicated that M.NruI is possibly a N6mA-type amino-methyltransferase that most likely modifies the external A in the 5' TCGCGA 3' sequence. M.Sbo13I, however, is implicated as a probable N4mC-type methylase since plasmid carrying sbo13IM gene is not restricted by Mrr endonuclease and Sbo13I digestion is not blocked by Dam methylation of the overlapping site. The amino acid sequence of M.NruI and M.Sbo13I did not show significant sequence similarity to many known amino-methyltransferases in the α, β, and γ groups, except to a few putative methylases in sequenced microbial genomes.
The order of the conserved amino acid motifs (blocks) in M.NruI/M.Sbo13I is similar to the γ. group amino-methyltranferases, but with two distinct features: In motif IV, the sequence is DPPY instead of NPPY; there are two additional conserved motifs, IVa and Xa as extension of motifs IV and X, in this family of enzymes. We propose that M.NruI and M.Sbo13I form a subgroup in the γ group of amino-methyltransferases.
Among the four types of restriction-modification (R-M) systems discovered from microbial sources based on subunit complexity, ATP/GTP requirement, and methylation-dependency, the Type II restriction endonucleases (REases) are useful tools in cleaving DNA into specific fragments for gene cloning and analysis . Type II restriction endonuclease genes are usually accompanied by companion methylase genes encoding methylases that modify the same target sites to avoid self-destruction of genomic DNA or extra chromosomal DNA . In bacterial warfare, foreign or phage DNAs are unmodified by the host resident methylases and therefore subjected to restriction by the companion endonuclease, while "self" DNA is partially or fully modified and thus resistant to Type II endonucelase attack. Over 3500 R-M systems have been found in nature with approximately 300 unique specificities containing 4-8 bp recognition sequences . There are three major types of base modification in bacteria and archaea: 5mC-methylation of the cytosine pyrimidine ring carbon producing 5-methylcytosine, N4-methylcytosine (N4mC), and N6-methyladenine (N6mA) [4–6]. In 5mC methylases such as M.HhaI, there are ten conserved amino acid motifs (blocks) that are arranged in the order of I to X (circular permutation of motifs IX and X in 5mC methylases has also been observed previously) [7, 8]. Among the N4mC and N6mA methylases, there are three major groups of methylases, i.e, α, β, and γ, based on the order of the amino acid motifs involved in S-adenosyl-L-methionine (AdoMet, methyl donor) binding (motifs X-I-II-III), catalytic function (motifs IV-V-VI-VII-VIII), and DNA target recognizing (TRD) [9, 2]. In the α group of amino-methyltransferases, the AdoMet binding region precedes the TRD and the catalytic region (X-I-II-III-TRD-IV-V-VI-VII-VIII, Motif IV = DPPY). In the β group, the catalytic region is arranged before the TRD and AdoMet binding region (IV-V-VI-VII-VIII-TRD-X-I-II-III, motif IV = DPPY or SPPY). In the γ group, the conserved motifs are arranged in the order of AdoMet binding region, catalytic region, and TRD (X-I-II-III-IV-V-VI-VII-VIII-TRD, motif IV = NPPY).
The NruI restriction endonuclease (REase or R) and methyltransferase (methylase or M) are enzymes isolated from the bacterium Nocardia rubra (Comb D.G., Schildkraut I., Greenough L. unpublished results cited in ). The NruI endonuclease binds to the symmetric sequence 5'-TCGCGA-3' in double-stranded DNA (dsDNA) and cleaves the DNA between the G and C in both strands (5'-TCG↓CGA-3'), thus producing DNA fragments with blunt ends. Sbo13I endonuclease is an NruI isoschizomer found in the strain Shigella boydii C13 (strain # NCTC 9361) . The goal of this work was cloning of the NruI R-M system in E. coli. During the initial cloning attempt, the methylase gene selection strategy  was not successful in cloning the nruIM gene, probably as a result of poor expression of nruIM in E. coli. Therefore, we sequenced the NruI endonuclease protein and obtained partial N-terminal amino acid sequence, which was then used to design degenerate primers for inverse PCR amplification of the coding sequence. Additional PCR reactions were carried out to amplify the entire NruI R-M system and the nruIR gene was successfully cloned in E. coli. The Sbo13I R-M system was cloned in E. coli by phage selection from plasmid expression libraries. When the M.NruI and M.Sbo13I amino acid sequences were searched against other known amino-methyltransferases in protein database in BlastP analysis , very little sequence similarity was detected, except five putative methylases. We propose that M.NruI and M.Sbo13I should be included in the γ group of amino-methyltransferases with two distinct features: motif IV is DPPY instead of NPPY, and two additional amino acid motifs IVa and Xa are also present, which might be involved in DNA target recognition or catalytic activity.
Strains, plasmid vectors, enzymes, primers, and genomic DNA preparation
E. coli strain ER2683 is an MM294 derivative (F'ProA+B+ lacI q ΔlacZM15 miniTn10 (KanR) fhuA2 Δ(lacI-lacA)200 glnV44 e14- rfbD1? relA1? endA1 spoT1? thi-1 Δ(mcrC-mrr)114::IS10) (strain provided by Lise Raleigh, NEB). Plasmid vectors pUC19, pACYC184, restriction enzymes, modification enzymes, Taq and Vent® DNA polymerases were produced at New England Biolabs, Inc (NEB). Degenerate inverse PCR primers were designed from the N-terminal NruI amino acid sequence (MGFLAD XDLSYDEINE LLTDN). The following primers were made by the NEB Organic Synthesis Division for inverse PCR amplification of the nruIR coding sequence:
Pho RTCNGCYAARAANCCCA (R = A/G; N = A/T/C/G; Y = T/C; Pho = 5'Phosphate)
(P270, mixed primers corresponding to aa MGFLAD)
Pho RTCNGCNAGTAANCCCA (P271, mixed primers corresponding to aa MGFLAD)
Pho TAYGAYGARATHAAYGA (P272, mixed primers corresponding to aa YDEINE).
Two-step gradient inverse PCR was carried out using primer pair P270 + P272 or P271 + P272 (DNA denaturing temperature at 95°C for 30 seconds, annealing temperature gradient at 37°C to 55°C for 30 seconds with 0.5°C increase in every cycle, 72°C extension for 2 min, 4 units of Taq DNA polymerase, 200 to 400 ng of endonuclease-digested and self-ligated genomic DNA template for 35 thermocycles).
Genomic DNA was prepared from 10 g of Nocardia rubra cells by a modified procedure based on the references [13, 14]. Cell paste was resuspended in 35 ml of 0.1 M Tris-HCl, pH 7, 0.1 M EDTA. Cell lysis was carried out by addition of 25 ml of 2 mg/ml fresh lysozyme in 0.1 M Tris-HCl, 0.1 M EDTA, pH 7.6, with incubation at 37°C for 1 hour. Protease K was added to the cell lysate at 0.1 mg/ml, with incubation at 37°C for 1 hour. To further facilitate cell lysis, SDS and sarcosyl solutions were added to 0.1% and 0.9% final concentration, respectively, with incubation at 55°C for 1 hour. The lysed cells were subjected to Phenol-CHCl3 extraction (3 ×), and CHCl3 extraction (2×). The nucleic acids were dialyzed in 4 L of TE buffer (10 mM Tris-HCl, pH 7.5, 0.1 mM EDTA) at 4°C overnight (2×). RNA was removed by digestion with RNaseA at 37°C for 1 hour. The quality of genomic DNA was analyzed on a 0.8% agarose gel. The size of the genomic DNA was estimated to be larger than 10 kb.
E. coli cell extracts were prepared as described previously . Native NruI endonuclease was purified from Nocardia rubra cell lysate by chromatography through following columns: heparin Hyper-D, Source Q, Affigel blue, ceramic HTP and Mono Q columns. N-terminus amino acid sequences were determined by the Edman degradation method . NruI endonuclease activity was carried out in a 1× NruI buffer: 100 mM KCl, 50 mM Tri-HCl, pH 7.7, 10 mM MgCl2, in a total volume of 50 μl at 37°C for 1 h on DNA substrate. Plasmid DNA and inverse PCR products were sequenced using the BigDye terminator cycle sequencing kit (Applied Biosystems). Dot blot analysis using N6mA or N4mC antibodies was carried out as described previously .
Multiple sequence alignment of M.NruI and M. Sbo13I, and methylase homologs were performed using the PROMALS3D web server . The amino acid sequences of NruI and Sbo13I and putative isoschizomers were analyzed using the Clustal W web server . MAFFT (v 6.0), using the minimum linkage method, was used to calculate distances from the PROMALS3D alignment and an unrooted distance tree was constructed for the M.NruI/M.Sbo13I family of N6mA methylases . Cluster analysis of multiple sequences was carried out using CLANS to visualize the clustering (formation of closely related protein families) among the M.NruI/M.Sbo13I group and known N6mA methyltransferases. CLANS is a Java utility tool based on the Fruchterman-Reingold graph layout algorithm .
Results and Discussion
Attempt to clone NruI endonuclease gene (nruIR) and methylase gene (nruIM) by methylase selection method
At first, the methylase selection method [11, 13] was used to clone the nruIM gene, but without success. It's known that M.NruI exists in the native strain because the genomic DNA from Nocardia rubra is resistant to NruI digestion (data not shown). However, the cloning of the NruI methylase gene (nruIM) proved to be extremely difficult. Dozens of genomic DNA libraries were constructed with different vectors and various enzyme digested genomic DNA fragments. The plasmid DNA libraries were subjected to NruI digestion and the digested DNA was then transferred into a McrBC- Mrr- E. coli host. But no true resistant clones were ever recovered (data not shown). The negative results may be due to the low expression of the methylase gene in E. coli and poor modification of NruI sites on the plasmid. It's known that Nocardia has a very high GC contents (the GC contents of Nocardia farcinica IFM10152 are 70%) , and genes from Nocardia may not be expressed efficiently in E. coli.
Purification of native NruI and N-terminal amino acid sequencing
Native NruI was purified from Nocardia rubra cell lysates by chromatography through the following columns: heparin Hyper-D, Source Q, Affigel blue, ceramic HTP and Mono Q columns. The starting crude extracts had 5.6 × 106 units, and the final purified enzyme contained 1.5 × 106 units. The enzyme preparation was only functionally pure, free of DNA, RNA and other nonspecific nucleases (data not shown). The purified NruI endonuclease was subjected to SDS-PAGE and more than ten protein bands were identified (data not shown). The four bands between 25 kDa and 47.5 kDa were selected and their N-terminus amino acid sequences were determined by Edman degradation method (most Type IIP restriction enzymes are in the range of 20 to 50 kDa in molecular masses, see REBASE). Only the protein band close to 25 kDa produced a protein sequence: MGFLADXDLSYDEINELLTDN (X, unidentified amino acid residue). The rest of the protein bands did not yield meaningful data (data not shown).
Design of degenerate inverse PCR primers and amplification of nruIR gene by inverse PCR
Genomic DNA was prepared from Nocardia rubra by phenol-chloroform extractions. One μg of genomic DNA was digested with individual restriction enzymes with 4-6 recognition sequences. Digested DNA was then purified and self-ligated at a low DNA concentration. The ligated DNA was further purified by Qiagen spin columns. The final volume was 50 μl and 10 μl was used as the template for degenerate inverse PCR.
The remaining DNA of the NruI R-M system was obtained by further inverse PCRs. The NruI R-M system consists of R and M genes running in the opposite direction and the two genes were separated by a large gap (326-bp non-coding sequence) (Figure 1). The nruIM gene is 804 bp and encodes a polypeptide of 267 amino acids, and the nruIR gene is 660 bp coding for a polypeptide of 219 amino acids (Figure 1). The N-terminus of the translated nruIR is MGFLEDWDLSYDEINELTDN (Figure 2), which closely matches the N-terminal aa sequence derived from the native NruI endonuclease (MGFLADXDLSYDEINELLTDN). Another new finding was that when the amino acid sequence of M.NruI was compared to other known amino-methyltransferases in GenBank and REBASE by blastx/blastp, M.NruI only shows sequence similarity to M.Sbo13I and a few other putative methylases. M.NruI/M.Sbo13I and a few others share a unique amino acid sequence similarity among the DNA amino-methyltransferases. This high degree of amino acid sequence similarity was not found in other N4C or N6A methylases (more detailed sequence comparison below).
BlastP analysis using NruI as a query also identified two additional putative endonucleases Cma ORF1763 (GenBank ID: ZP_04834997) from Corynebacterium matruchotii ATCC 14266 (annotated as a helix-turn-helix XRE-family DNA binding protein) and RspRS ORF4066 found in the sequenced genome of Roseiflexus sp. RS-1 (GenBank accession number CP000686, coordinate 5095066-5095614, complement strand sequence). Figure 2 shows the Clustal W multiple sequence alignment of NruI, Sbo13I, Cma ORF1763 (ZP_04834997), and RspRS ORF4066R. The predicted catalytic sites in NruI, Sbo13I and Cma ORF1763 (ZP_04834997) are PD-X21-EVK or DD-X21-EVK. In RspRS ORF4066R, however, the predicted catalytic residues are DD-X20-QIK (Mrr endonuclease-like catalytic site) . The endonuclease activity of Cma ORF1763 (ZP_04834997) and RspRS ORF4066R remains to be examined by experimentation. There are putative companion methylases (ZP_04834998; YP_001278360 = RspRS_ORF4066M, complete genome coordinate 5094219 to 5095064) adjacent to Cma ORF1763 (ZP_04834997) and RspRS ORF4066R respectively that show significant sequence similarity to M.NruI/M.Sbo13I (see below). The gene organization of the two putative NruI isoschizomers is shown in Figure 1. Similar to NruI and Sbo13I, the homologous R-M genes are transcribed in opposite directions, but lacking a large non-coding region between the R-M genes.
Cloning of Sbo13I Restriction modification system in E. coli
Sbo13I restriction endonuclease and modification methylase were derived from Shigella boydii C13 . Sbo13I and NruI recognize the same DNA sequence TCG/CGA and cleave at the same position. Previous results indicated that the Sbo13I R-M system was encoded by a resident plasmid (data not shown). Plasmid DNA of Shigella boydii C13 was prepared by the cesium chloride-ethidium bromide equilibrium density centrifugation method. One isolate (A5) was used to map the 5.4 kb restriction system-containing plasmid for single restriction sites. The purified A5 isolate DNA was then digested with ClaI, HindIII and NsiI, respectively and ligated into pUC19 or pBR322 (CIP-treated with compatible ends). The ligated DNA was used to transform E. coli strain RR1 and selected on LB agar Amp plates. AmpR transformants were picked into 200 μl of LB with antibiotic into micro-titer plates. Master plates were prepared by stamping on LB with antibiotic. Replica plates were prepared by stamping onto four levels of T7 phage: 109, 107, 105, and 103 phage/plate. Individual colonies which survived at all levels of phage infection were considered to be phage-resistant and hence likely to carry an active restriction system.
A number of plasmids were found to carry ClaI or NsiI fragments of approximately 5.4 kb in length and these plasmids were resistant to Sbo13I digestion. These plasmids were subsequently shown to carry the complete Sbo13I R-M system. The recombinant plasmid pBLSboC13M8.1 which carries the gene encoding the Sbo13I restriction endonuclease and methylase was transferred into E. coli strain RR1 by transformation. Approximately, 105 units of Sbo13I REase per gram of wet cell paste were produced from the expression strain (data not shown). The sbo13IR and sbo13IM genes were sequenced by the Sanger sequencing method (data not shown). The gene organization of the Sbo13I R-M system is shown in Figure 1. Similar to the NruI R-M system, the sbo13IR and sbo13IM genes are transcribed in opposite directions. In addition, there is a controller gene sbo13IC in front of the sbo13IR gene, which may regulate the transcription of the R gene, by analogy to the PvuII and BclI controller proteins [25, 26].
The NruI endonuclease was expressed in pre-modified E. coli hosts and the recombinant NruI endonuclease was purified (ZZ, SYX, unpublished results).
Dot blot analysis to determine N6mA modification for M.NruI and M.Sbo13I
Sequence comparison of M.NruI and M. Sbo13I
Comparison of M.NruI and M.Sbo13I amino acid sequences indicates that they are significantly similar to each other (49% aa sequence identity). Furthermore, these enzymes may form a new family (cluster) of methyltransferases based on their lack of similarity to many other known methyltransferases available in REBASE. Alternatively, they may belong to a subgroup of already established α, β and γ groups of N6mA methyltransferases (Malone et al. 1995).
The other members of this family (cluster) were identified by querying M.NruI and M.Sbo13I protein sequences against non-redundant protein sequences (NR) and REBASE databases using BlastP searches. Significant hits with an expectation value <1E-20 and sequence length of more than 200 amino acids were considered potential M.NruI and M.Sbo13I homologs (putative isoschizomer). M.NruI and M.Sbo13I enzymes share amino acid sequence similarities with proteins YP_001278360 (RspRS_4066 from Roseiflexus sp. RS-1), ZP_01906951 (PPSIR1_36012 from Plesiocystis pacifica SIR-1), ZP_03711821 (CORMATOL_02672 from Corynebacterium matruchotii), YP_001432817 (Rcas_2728 from Roseiflexus castenholzii DSM 13941). The GenBank protein and gene IDs are described as following: ZP_01906951 (gi_149918461, PPSIR1_36012 from Plesiocystis pacifica SIR-1), ZP_03711821 (gi_225022629, CORMATOL_02672 from Corynebacterium matruchotii), ZP_04834998 (gi_252123721, M.Cma ORF1763, ATCC_14266 from Corynebacterium matruchotii), and M.Rca13941ORF2728P (DSM 13941 from Roseiflexus castenholzii).
Motif I contains a core G-loop (Gly-X-Gly) tripeptide that may bind the methionine moiety of AdoMet. X indicates any amino acid.
Motif II contains a charged residue at the end of a predicted β strand which may interact with ribose hydroxyls of AdoMet. It is followed by a bulky hydrophobic side-chain (Ile or Leu) that may make van der Waals contact with the adenine of AdoMet.
Motif IV contains a P-loop, DPPY motif, which may form the active site, along with motifs V to VIII. Superficially, the DPPY motif suggests that M.NruI and M.Sbo13I may be related to the α group of N6mA methyltransferases. However, the phylogenetic distance analysis indicates otherwise.
Motif V contains the consensus Asp-Leu-Tyr-X-X-Ala-(Gly/Ser) which differs from the standard motif V (Asn/Asp-Leu-Tyr-X-X-Phe-(Leu/Val/Ile) defined in the literature for the γ group of N6mA methyltransferases. The first three amino acids are the same as standard amino acids, but the last two amino acids are different.
Motif VI starts with Gly and ends with a hydrophobic amino acid (Val) which is very similar to the standard motif defined in literature.
Motif VII is highly hydrophobic which is consistent with the standard motif. This hydrophobic region may interact with the target DNA adenine.
Motif X is located at the N-terminus of the primary sequence and contains conserved hydrophobic residues.
Both NruI and Sbo13I R-M systems have been cloned in E. coli. BlastP analysis of proteins in GenBank indicated that two ORFs are probable isoschizomers of NruI/SboI or enzymes with similar recognition sequence. BlastP analysis of proteins in GenBank using the M.NruI amino acid sequence as a query revealed five additional putative amino-methylases that share similar motif organization. The M.NruI/M.Sbo13I family enzymes appear to belong to the γ group of amino-methyltransferases with two distinct features: Motif IV contains DPPY instead of NPPY sequence; two additional motifs (IVa and Xa) are also present, which may be involved in target recognition and AdoMet binding. Despite the extensive amino acid sequence similarity between M.NruI and M.Sbo13I, M.NruI is most likely an N6mA amino-methyltransferase and M.Sbo13I is likely an N4mC methyltransferase.
We thank Rich Roberts and Lise Raleigh for critical comments; Geoff Wilson and Bill Jack for helpful discussions; Brian Anton and Janos Posfai for help with amino acid motif analysis; Don Comb and Jim Ellard for financial support and encouragement. This work was supported by New England Biolabs, Inc. Request for NruI clone should be directed to firstname.lastname@example.org; and request for Sbo13I clone should be addressed to email@example.com.
- Roberts RJ, Belfort M, Bestor T, Bhagwat AS, Bickle TA, Bitinaite J, Blumenthal RM, Degtyarev S, Dryden DT, Dybvig K, et al: A nomenclature for restriction enzymes, DNA methyltransferases, homing endonucleases and their genes. Nucleic Acids Res. 2003, 31 (7): 1805-1812. 10.1093/nar/gkg274.PubMed CentralPubMedView ArticleGoogle Scholar
- Wilson GG, Murray NE: Restriction and Modification Systems. Ann Rev Genet. 1991, 25: 585-627. 10.1146/annurev.ge.25.120191.003101.PubMedView ArticleGoogle Scholar
- Roberts RJ, Vincze T, Posfai J, Macelis D: REBASE--a database for DNA restriction and modification: enzymes, genes and genomes. Nucleic Acids Res. 2009, D234-236. 38 Database
- Cheng X: Structure and function of DNA methyltransferases. Annu Rev Biophys Biomol Struct. 1995, 24: 293-318. 10.1146/annurev.bb.24.060195.001453.PubMedView ArticleGoogle Scholar
- Vanyushin BF, Belozersky AN, Kokurina NA, Kadirova DX: 5-methylcytosine and 6-methylamino-purine in bacterial DNA. Nature. 1968, 218 (5146): 1066-1067. 10.1038/2181066a0.PubMedView ArticleGoogle Scholar
- Janulaitis A, Klimasauskas S, Petrusyte M, Butkus V: Cytosine modification in DNA by BcnI methylase yields N4-methylcytosine. FEBS Lett. 1983, 161 (1): 131-134. 10.1016/0014-5793(83)80745-5.PubMedView ArticleGoogle Scholar
- Posfai J, Bhagwat AS, Posfai G, Roberts RJ: Predictive motifs derived from cytosine methyltransferases. Nucleic Acids Res. 1989, 17: 2421-2435. 10.1093/nar/17.7.2421.PubMed CentralPubMedView ArticleGoogle Scholar
- Xu Sy, Xiao J-p, Posfai J, Maunus R, Benner J: Cloning of the BssHII restriction-modification system in Escherichia coli: BssHII methyltransferase contains circularly permuted cytosine-5 methyltransferase motifs. Nucleic Acids Res. 1997, 25: 3991-3994. 10.1093/nar/25.20.3991.PubMed CentralPubMedView ArticleGoogle Scholar
- Malone T, Blumenthal RM, Cheng X: Structure-guided analysis reveals nine sequence motifs conserved among DNA amino-methyl-transferases, and suggests a catalytic mechanism for these enzymes. J Mol Biol. 1995, 253: 618-632. 10.1006/jmbi.1995.0577.PubMedView ArticleGoogle Scholar
- Mise K, Nakajima K, Terakado N, Ishidate M: Production of restriction endonucleases using multicopy Hsd plasmids occurring naturally in pathogenic Escherichia coli and Shigella boydii. Gene. 1986, 44 (1): 165-169. 10.1016/0378-1119(86)90058-2.PubMedView ArticleGoogle Scholar
- Szomolanyi I, Kiss A, Venetianer P: Cloning the modification methylase gene of Bacillus sphaericus R in Escherichia coli. Gene. 1980, 10: 219-225. 10.1016/0378-1119(80)90051-7.PubMedView ArticleGoogle Scholar
- Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL: BLAST+: architecture and applications. BMC Bioinformatics. 2009, 10: 421-10.1186/1471-2105-10-421.PubMed CentralPubMedView ArticleGoogle Scholar
- Lunnen KD, Barsomian JM, Camp RR, Card CO, Chen SZ, Croft R, Looney MC, Meda MM, Moran LS, Nwankwo DO, et al: Cloning Type II restriction and modification genes. Gene. 1988, 74: 25-32. 10.1016/0378-1119(88)90242-9.PubMedView ArticleGoogle Scholar
- Xu Sy, Xiao J-p, Ettwiller L, Holden M, Aliotta J, Poh CL, Dalton M, Robinson DP, Petronzio TR, Moran L, et al: Cloning and expression of the ApaLI, NspI, NspHI, SacI, ScaI, and SapI restriction-modification systems in Escherichia coli. Mol Gen Genet. 1998, 260: 226-231. 10.1007/s004380050890.PubMedView ArticleGoogle Scholar
- Xu S-y, Schildkraut I: Isolation of BamHI variants with reduced cleavage activities. J Biol Chem. 1991, 266: 4425-4429.PubMedGoogle Scholar
- Edman P: Sequence determination. Mol Biol Biochem Biophys. 1970, 8: 211-255.PubMedGoogle Scholar
- Kong H, Lin LF, Porter N, Stickel S, Byrd D, Posfai J, Roberts RJ: Functional analysis of putative restriction-modification system genes in the Helicobacter pylori J99 genome. Nucleic Acids Res. 2000, 28 (17): 3216-3223. 10.1093/nar/28.17.3216.PubMed CentralPubMedView ArticleGoogle Scholar
- Pei J, Kim BH, Grishin NV: PROMALS3D: a tool for multiple protein sequence and structure alignments. Nucleic Acids Res. 2008, 36 (7): 2295-2300. 10.1093/nar/gkn072.PubMed CentralPubMedView ArticleGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23 (21): 2947-2948. 10.1093/bioinformatics/btm404.PubMedView ArticleGoogle Scholar
- Katoh K, Asimenos G, Toh H: Multiple alignment of DNA sequences with MAFFT. Methods Mol Biol. 2009, 537: 39-64. full_text.PubMedView ArticleGoogle Scholar
- Frickey T, Lupas A: CLANS: a Java application for visualizing protein families based on pairwise similarity. Bioinformatics. 2004, 20 (18): 3702-3704. 10.1093/bioinformatics/bth444.PubMedView ArticleGoogle Scholar
- Ishikawa J, Yamashita A, Mikami Y, Hoshino Y, Kurita H, Hotta K, Shiba T, Hattori M: The complete genomic sequence of Nocardia farcinica IFM 10152. Proc Natl Acad Sci USA. 2004, 101 (41): 14925-14930. 10.1073/pnas.0406410101.PubMed CentralPubMedView ArticleGoogle Scholar
- Wilson GG: Cloned restriction-modification systems - a review. Gene. 1988, 74: 281-289. 10.1016/0378-1119(88)90304-6.PubMedView ArticleGoogle Scholar
- Bujnicki JM, Rychlewski L: Identification of a PD-(D/E)XK-like domain with a novel configuration of the endonuclease active site in the methyl-directed restriction enzyme Mrr and its homologs. Gene. 2001, 267 (2): 183-191. 10.1016/S0378-1119(01)00405-X.PubMedView ArticleGoogle Scholar
- Knowle D, Lintner RE, Touma YM, Blumenthal RM: Nature of the Promoter Activated by C.PvuII, an Unusual Regulatory Protein Conserved among Restriction-Modification Systems. J Bacteriol. 2005, 187 (2): 488-497. 10.1128/JB.187.2.488-497.2005.PubMed CentralPubMedView ArticleGoogle Scholar
- Sawaya MR, Zhu Z, Mersha F, Chan SH, Dabur R, Xu SY, Balendiran GK: Crystal structure of the restriction-modification system control element C.Bcll and mapping of its binding site. Structure. 2005, 13 (12): 1837-1847. 10.1016/j.str.2005.08.017.PubMedView ArticleGoogle Scholar
- Brooks JE, Nathan PD, Landry D, Sznyter LA, Waite-Rees P, Ives CL, Moran LS, Slatko BE, Benner JS: Characterization of the cloned BamHI restriction modification system: its nucleotide sequence, properties of the methylase, and expression in heterologous hosts. Nucleic Acids Res. 1991, 19: 841-850. 10.1093/nar/19.4.841.PubMed CentralPubMedView ArticleGoogle Scholar
- Samuelson JC, Xu SY: Directed evolution of restriction endonuclease BstYI to achieve increased substrate specificity. J Mol Biol. 2002, 319 (3): 673-683. 10.1016/S0022-2836(02)00343-1.PubMedView ArticleGoogle Scholar
- Patel Y, Nelson M, McClelland M: Methylation at overlapping dam (Gm6ATC) sites does not block cleavage by the NruI (TCGCGA) isoschizomer restriction endonucleases AmaI, SalDI and Sbo13I. Nucleic Acids Res. 1989, 17 (9): 3613-10.1093/nar/17.9.3613.PubMed CentralPubMedView ArticleGoogle Scholar
- Degtyarev S, Prikhod'ko EA, Rechkunova NI, Prikhod'ko GG, Krasnykh VN: Biochemical characterization of VspI methyltransferase. Gene. 1995, 157 (1-2): 65-66. 10.1016/0378-1119(94)00710-A.PubMedView ArticleGoogle Scholar