- Research article
- Open Access
A genetic dissection of the LlaJI restriction cassette reveals insights on a novel bacteriophage resistance system
© O'Driscoll et al; licensee BioMed Central Ltd. 2006
- Received: 21 March 2006
- Accepted: 28 April 2006
- Published: 28 April 2006
Restriction/modification systems provide the dual function of protecting host DNA against restriction by methylation of appropriate bases within their recognition sequences, and restriction of foreign invading un-methylated DNA, such as promiscuous plasmids or infecting bacteriphage. The plasmid-encoded LlaJI restriction/modification system from Lactococcus lactis recognizes an asymmetric, complementary DNA sequence, consisting of 5'GACGC'3 in one strand and 5'GCGTC'3 in the other and provides a prodigious barrier to bacteriophage infection. LlaJI is comprised of four similarly oriented genes, encoding two 5mC-MTases (M1.LlaJI and M2.LlaJI) and two subunits responsible for restriction activity (R1.LlaJI and R2.LlaJI). Here we employ a detailed genetic analysis of the LlaJI restriction determinants in an attempt to characterize mechanistic features of this unusual hetero-oligomeric endonuclease.
Detailed bioinformatics analysis confirmed the presence of a conserved GTP binding and hydrolysis domain within the C-terminal half of the R1.LlaJI amino acid sequence whilst the N-terminal half appeared to be entirely unique. This domain architecture was homologous with that of the "B" subunit of the GTP-dependent, methyl-specific McrBC endonuclease from E.coli K-12. R1.LlaJI did not appear to contain a catalytic centre, whereas this conserved motif; PD....D/EXK, was clearly identified within the amino acid sequence for R2.LlaJI. Both R1.LlaJI and R2.LlaJI were found to be absolutely required for detectable LlaJI activity in vivo. The LlaJI restriction subunits were purified and examined in vitro, which allowed the assignment of R1.LlaJI as the sole specificity determining subunit, whilst R2.LlaJI is believed to mediate DNA cleavage.
The hetero-subunit structure of LlaJI, wherein one subunit mediates DNA binding whilst the other subunit is predicted to catalyze strand hydrolysis distinguishes LlaJI from previously characterized restriction-modification systems. Furthermore, this distinction is accentuated by the fact that whilst LlaJI behaves as a conventional Type IIA system in vivo, in that it restricts un-methylated DNA, it resembles the Type IV McrBC endonuclease, an enzyme specific for methylated DNA. A number of similar restriction determinants were identified in the database and it is likely LlaJI together with these homologous systems, comprise a new subtype of the Type II class incorporating features of Type II and Type IV systems.
- High Molecular Weight Complex
- GTPase Domain
- Phage Resistance
- Intrinsic GTPase Activity
- Lactococcal Strain
The diverse structural and functional activities described for restriction endonucleases (REases) is reflected by their ubiquitous abundance among prokaryotic organisms . These enzymes are typically components of restriction/modification (R/M) systems where they primarily serve as defensive instruments against bacteriophage infection and other invading "foreign" DNA.
Four distinct groups of REases are known, with the "Type II" group being the most richly represented . Members of this group, of which there are eleven subtypes, generally recognize relatively short DNA sequences and specifically cut within or reasonably close to such sequences in a Mg+2 dependent manner.
Although generally lacking significant amino acid sequence homology, structural analysis of many of these REases has revealed in some instances a conserved core, consisting of a four/five-stranded mixed β-sheet flanked by α-helices, reviewed in . This core spatially localizes two carboxylates (present in the amino acids D and D/E) which are believed to be involved in divalent metal ion binding, and a lysine residue, possibly involved in positioning and/or activation of an attacking nucleophile . Together it is believed that these residues comprise the REase catalytic site, the PD...D/EXK motif, responsible for phosphodiester bond hydrolysis in the target DNA molecule.
For double-stranded scission this hydrolysis event must occur on both strands of the target DNA thereby requiring two catalytic sites , although single-site exceptions do exist [12, 13]. Consequently, the active forms of many Type II REases are homodimers (or even tetramers), with each identical subunit providing a single catalytic site for individual strand hydrolysis [9, 14]. Of course given the multifarious nature of the Type II group, it is not surprising that unusual hetero-multimeric REases, comprised of unique subunits, consisting of the associated products of different genes, have been described [1, 3, 15–18].
One such characterized enzyme, BbvCI, has been shown to be heterodimeric, with each individual subunit independently responsible for cleavage of separate strands of its asymmetric recognition sequence [15, 19]. Similarly, the Bpu10I REase requires two different subunits for its activity however, in this instance, a direct interaction between the subunits has yet to be completely established .
Previously we have described the cloning and transcriptional regulation of LlaJI R/M system [1, 20]. The LlaJI gene cluster, encoded by the lactococcal plasmid pNP40, consists of two 5mC MTase-encoding genes followed by two REase-encoding genes which mediate a prodigious barrier to bacteriophage infection . The specificity of the system was determined to be asymmetric but complementary, comprised of 5'GACGC'3 in one strand and 5'GCGTC'3 in the other . Perhaps the most striking feature of this R/M system is its homology to the methyl-specific REase McrBC. This highly unusual enzyme, which lacks an associated MTase, has been characterized extensively [10, 21–31]. Its orphan nature, unusual specificity and cofactor requirements define McrBC, together with its identified homologues, as a distinct REase group .
The mcrBC locus encodes two major products, McrB (53 kDa) and McrC (39 kDa), which mediate the restriction phenotype . McrBC recognizes two methylated/hemimethylated RmC sites (separated by 40 bp to up to 3000 bp) and cleaves at random positions between the two [24, 25, 32]. The optimal stoichiometry for cleavage in vitro has been suggested as 3–5:1, McrB:McrC (although a 1:1 ratio has been reported by another group ), with the cleavage reaction strictly dependent upon GTP and Mg+2 [25, 33].
The McrB subunit harbours the specific DNA binding activity (which resides within its N-terminus), in addition to GTP binding and hydrolysis activities [22, 30, 34]. The McrC subunit contains the catalytic centre for cleavage . Following binding to its methylated target site, McrBC facilitates a GTP-dependent translocation of the DNA, akin to the mechanisms of other NTP-dependent REases of the Type I and III groups, with cleavage ensuing upon stalling of the translocating complex [23, 35]. Furthermore, it has been shown that McrB forms high molecular weight complexes – heptameric rings, in the presence of GTP (although independent of GTP hydrolysis). These complexes can further multimerise to tetradecamers, structures stabilized by the binding of McrC – which constitutes the "active" form of the endonuclease .
In this report, we execute a detailed genetic and phenotypic assessment of the McrBC-like LlaJI endonuclease which serves to portray the complexity of this unique enzyme.
Sequence analysis of LlaJI
This domain has been investigated biochemically in the case of the R1.LlaJI homologue, the McrB subunit of the McrBC endonuclease – and has revealed interesting variations from the archetypal GTPase family (see discussion) .
Additional conserved modules (of identity) can be discerned within R1.LlaJI and its homologues; however, at this point no functional significance can be assigned to any of these regions (Fig 1B). It was suggested for McrB, that some of these conserved modules, at the time identified using completely different protein sets, may contribute in some manner to the diverged functionality of the GTP-binding domain . Hence, it is likely that, as for McrB, R1.LlaJI represents yet another member of a distinct subclass of the classical GTPase family.
For most of the identified R1.LlaJI homologues, significant alignment was only achieved with portions of the sequences corresponding to the C-terminal half of each, with the N-terminal portions apparently representing non-conserved sequence. The two exceptions to this, namely LlaIB and BsuMIB (Fig 1A and 1B), which are smaller (337 and 343 aa's, respectively) than R1.LlaJI (585 aa), aligned completely along the full length of their respective sequences. Interestingly, both appear to be encoded by "restriction cassettes" comprised of three ORF's as opposed to the two observed for the remaining homologues (Fig 1A).
Database searches with R2.LlaJI (48.8 kDa) revealed a set of similar proteins that are encoded by genes located adjacent to those that had been identified as encoding homologues of R1.LlaJI (see Fig 1A and 1C). 1A multiple alignment of R2.LlaJI with its homologues allowed the identification of a putative catalytic centre for DNA cleavage; PD...D/ExK. This assignment is supported by inclusion of the McrC subunit of McrBC in the alignment, whose catalytic centre has essentially been confirmed . As for R1.LlaJI above, R2.LlaJI shares significant, but scattered modules of identity with other proteins, but again the functional significance of this similarity remains purely speculative.
These data cumulatively suggest that LlaJI, and indeed other restriction systems exhibit structural and therefore likely mechanistic similarities with that of the McrBC methyl-specific endonuclease (see discussion).
Phenotypic and genetic analysis of LlaJI
Previously, we have described the bacteriophage resistance phenotype conferred by the LlaJI operon . To fully explore the contribution of the individual LlaJI restriction "subunits" to this phage resistance activity, we adopted a genetic subcloning approach. To this end, the llaJIR1 and llaJIR2 genes were amplified individually, and in tandem, and placed under the control of the strong nisin-inducible promoter of plasmid pNZ8020 to facilitate their over-expression.
Bacteria, bacteriophage and plasmids used in this study.
Strains and Plasmids
Reference or source
Cloning host, KanR.
F'proA+B+ lacIq Δ lacZM15 miniTn10 (KanR) fhuA2 Δ(lacI-lacA)200glnV44e14-rfbD1? relA1? endA1 spoT1?thi-1 Δ(mcrC-mrr)114::IS10.
Expression host for pQE 60 derivatives.
Plasmid-free derivative of L.Lactis supsp. lactis 712
MG1363, pepN::nisRK, nisin inducible overexpression host.
Partially cured derivative of NCK203. Host strain for phage Ul36
E.coli temperate bacteriophage.
Small isometric headed phage for MG1614.
Prolate headed phage for MG1614.
P335 lytic phage for NCK690 KP1.
Low-copy number E.coli-L.Lactis shuttle vector for cloning, TcR.
pPTPi derivative containing complete 6.2 kb LlaJI operon from pNP40.
pPTPi derivative containing 3.0 kb fragment from pNP40 encompassing promoter, llaJIM1 and llaJIM2 genes.
High-copy number E.coli-L.Lactis overexpression vector, PnisA, CmR.
pNZ8020 derivative containing llaJIR1 under control of PnisA.
pNZ8020 derivative containing llaJIR2 under control of PnisA.
pNZ8020 derivative containing llaJIR1 and llaJIR2 under control of PnisA.
High-copy number E.coli cloning vector, AmpR.
pUC19 derivative containing llaJIR1.
pUC19 derivative containing llaJIR2.
pUC19 derivative containing llaJIR1 and llaJIR2.
Medium-copy number expression plasmid
pMal-c2x derivative containing llaJIR2.
High-number 6xhis tag expression plasmid.
pQE60 derivative containing llaJIR1.
pQE60 derivative containing llaJIR2.
Low-copy number LacI expressing pQE60 companion
LlaJI mediated bacteriophage resistance to the lactococcal phage Sk1.
Plaque size after 16 hrs (mm)
Plaque size after 16 hrs (mm)
3.3 × 10-3
2.6 × 10-6
1–2 and pinpoint
1 × 10-5
5 × 10-6
1–2 and pinpoint
The fact that this optimal level of resistance may be achieved at 30°C by mere over-expression of R1.LlaJI and R2.LlaJI suggests that the temperature sensitivity of the LlaJI system previously reported, may be attributable to the intracellular concentration of active LlaJI at the respective temperatures (modulated by the LlaJI natural transcriptional signals), rather than to the temperature "activity" spectrum of assembled functional LlaJI, as had previously been speculated .
As stated, the results obtained with the hosts expressing R1.LlaJI and R2.LlaJI individually (on plasmids pNZ-R1, and pNZ-R2, respectively) indicated that neither subunit alone was capable of restricting the growth of phage. However, despite the lack of significant decrease in EOP at either sampling temperature used, a delicate but visibly obvious variable plaque phenotype was observed for the host expressing R1.LlaJI alone (see Table 2). Perhaps most intriguing was the complete absence of any detectable endonuclease activity in crude extracts of lactococcal strains containing either pJO-J or the pNZ-R1R2 over-expression plasmid (which as stated was always accompanied with plasmid pJO-M1M2).
Purification of the individual LlaJI subunits
As it was established that both LlaJI subunits were required for "complete" activity, it was assumed that expression of the individual subunits separately for purification purposes could be achieved in the absence of the LlaJI MTases. To this end, both llaJIR1 and llaJIR2 were cloned in plasmid pQE60 to produce derivatives pQE-R1 and pQE-R2, respectively which were subsequently established in E.coli (see materials and methods). Over-expression with pQE-R1 allowed purification of R1.LlaJI as a C-terminal 6xhis tagged derivative. Mass spectrometry and N-terminal sequencing confirmed the predicted molecular mass (68.139 kDa inclusive of the mass of the 6x His tag) and identity of purified R1.LlaJI.
Purification of R2.LlaJI by this approach failed, possibly due to instability of the resulting recombinant protein. To overcome this presumed instability, llaJIR2 was inserted downstream of the malE gene of plasmid pMal-c2x under the control of the Ptac promoter (see materials and methods) which, upon induction resulted in the expression of an R2-MBP fusion of the predicted molecular mass (~90.8 kDa). This fusion was purified by affinity chromatography using amylose resin and the R2.LlaJI subunit was subsequently released from the R2-MBP preparation by cleavage with factor Xa protease. Rebinding of MBP to amylose resin resulted in purification of R2.LlaJI to ~60 % homogeneity.
Assignment of the specificity determining LlaJI subunit
Given the domain architecture shared between R1.LlaJI and McrB, and the established DNA binding activity of the latter, we investigated the possibility of R1.LlaJI constituting the DNA recognition element of LlaJI.
A similar analysis using partially purified R2.LlaJI did not reveal any DNA binding to the LlaJI site containing probes. In addition, R2.LlaJI did not appear to have any discernable effect on the R1.LlaJI binding activity.
Investigation of the biological activity of the purified LlaJI subunits
As mentioned above, it was not possible to detect specific endonuclease activity in crude cell extracts of lactococcal strains (over)expressing the LlaJI system (as measured by its phage resistance phenotype). Subsequent attempts to reconstitute active LlaJI by combination of the partially purified subunits in various ratios, and with an assortment of cofactors (see materials and methods) failed to reveal any specific nuclease activity. Furthermore, neither subunit was found to exhibit a DNA nicking activity. Possible explanations for these observations are discussed below.
The presence of the putative GTPase domain in R1.LlaJI prompted us to investigate whether or not this subunit displayed such an activity. Analysis of GTP hydrolysis, conducted by colorimetric assay (see materials and methods) suggested that R1.LlaJI did not exhibit an intrinsic GTPase activity. Surprisingly, R1.LlaJI dependent GTP hydrolysis could not be activated by the inclusion of "specifc" DNA and/or partially purified R2.LlaJI in the reaction. Similar results were obtained with alternative nucleotides.
Phenotypic analysis of heterologous LlaJI expression
LlaJI-mediated resistance to bacteriophage λ. Results reflect values obtained for assays conducted at 30°C.
7 × 10-2
Relationship between EOP and frequency of LlaJI recognition sites on selected phage genomes.
Genome size (kb)
No. LlaJI sites
3.3 × 10-3
1 × 10-4
3 × 10-2
7 × 10-2
Investigation of the temperature induced efficacy of LlaJI in E.coli further revealed that, although the observed resistance (with plasmid pUC-R1R2 at 30°C) was not evident at 37°C, it was not significantly enhanced with growth at low temperatures (data not shown). Again, as with the lactococcal hosts, crude extracts of E.coli over-expressing R1.LlaJI and R2.LlaJI (from plasmid pUC-R1R2) did not exhibit detectable endonuclease activity.
The molecular mechanism by which the lactococcal plasmid-encoded LlaJI R/M system achieves DNA restriction has remained largely elusive. Previous studies have characterized the phage resistance capabilities of this four-gene cluster which is accompanied by an intricate transcriptional regulatory mechanism [1, 20].
In the present study we offer additional phenotypic and genetic data which alludes to an unusual mode of action for LlaJI. Analysis of the deduced amino acid sequence for R1.LlaJI revealed a segmented architecture, with a C-terminally located putative GTP binding and hydrolysis domain, and a largely unique N-terminal portion. This organization was similar to that observed for the McrB subunit of the McrBC methyl-specific endonuclease [25, 37].
A mutational analysis of the three motifs which comprise the GTPase domain of McrB confirmed the assignment for motif I only, which is the most highly conserved of the three with respect to R1.LlaJI and its homologues (see Fig 1B). Conservative substitutions in motifs II and III failed to unequivocally confirm the proposed assignments on McrB . It was suggested that the proposed motif III, although essential for GTP binding and hydrolysis (N333 in particular), rather than define nucleotide specificity (predicted for D336), was more likely to be involved in McrC interaction (a suggestion corroborated by the significant McrC-stimulated GTPase activity of McrB) . The apparent replacement of the conserved lysine of motif III with a threonine (a divergence also preserved on R1.LlaJI, Fig1A) further substantiated the likelihood of this motif contributing an atypical function.
In addition, the relative spacing between motif's II and III on McrB (29 amino acids), appeared too short to be structurally analogous to the classical GTPases (>49 amino acids)[29, 38]. However, in the case of R1.LlaJI, the spacing between these proposed motifs is, although relatively large, more canonical (112 amino acids), and in fact the alignment reveals an alternative candidate for the McrB motif II (although the conserved glycine would be replaced by glutamine, Fig 1B).
It has been shown that it is the N-terminal 160 amino acids of McrB which, like R1.LlaJI did not exhibit significant homology to the database, specified the DNA binding activity of McrBC . This, in addition to the placement of the GTPase domain suggested that an equivalent activity could be attributable to the unique N-terminal portion of R1.LlaJI. We examined these predicted properties for R1.LlaJI and indeed established a specific DNA binding activity for this subunit to LlaJI-site containing DNA fragments.
Interestingly, and in contrast to McrB, R1.LlaJI did not exhibit an intrinsic GTPase activity, which could neither be stimulated by specific DNA nor partially purified R2.LlaJI. This latter finding was again in contrast to McrB, whose intrinsic GTPase activity was stimulated 30-fold by McrC . A putative catalytic centre was identified within the amino acid sequence for R2.LlaJI, which suggested a role for this subunit in substrate strand hydrolysis; however, neither R2.LlaJI nor R1.LlaJI (and combinations of the two) exhibited cleavage or nicking activity.
The absolute requirement for both R1.LlaJI and R2.LlaJI for complete in vivo activity was clearly demonstrated (Table 2). A mild phenotype was associated with over-expression of R1.LlaJI alone (in Lactococcus), although this was most likely due to DNA binding only, which may have served to partially hinder the phage replication machinery.
Although less dramatic, heterologous expression of R1.LlaJI and R2.LlaJI in E.coli provided resistance to phage infection where, again an absolute requirement for both R1.LlaJI and R2.LlaJI was observed (Table 3). The relative poor efficacy of the resistance provided in this host background was attributed to instability of the R2.LlaJI subunit. This subunit, in contrast to R1.LlaJI, could not be observed in over-expressing cell extracts when analyzed by SDS-PAGE.
Purification of R2.LlaJI was only achieved when stabilized by fusion to a carrier protein. It is possible that in this fusion correct folding of the R2.LlaJI subunit had been inhibited, and was not subsequently corrected following the proteolytic removal of the carrier. This could account for our inability to reconstitute an active endonucleolytic complex with R1.LlaJI despite the clearly evident DNA binding activity of the latter.
Furthermore, it is possible that the predicted GTP hydrolyzing activity (which presumably facilitates DNA translocation) of R1.LlaJI would only become apparent following assembly of the "active" heterooligomer on its target DNA. It must be noted however, that an additional lactococcal derived cofactor, proteinacious or otherwise, can not be entirely excluded at this point, despite the absence of endonuclease activity in LlaJI-expressing lactococcal crude extracts.
In summary, LlaJI is a heterooligomeric endonuclease wherein the R1.LlaJI subunit is responsible for determining specificity. It is suggested that R2.LlaJI harbours the catalytic centre for DNA cleavage, whose site, although undetermined, is likely to be defined by the relative "placement" of R2.LlaJI on the target DNA (due to its proposed interaction with R1.LlaJI).
LlaJI was shown to retain behavioral characteristics of a conventional Type IIA system, in that unmethylated DNA is targeted in vivo however; unlike any characterized enzyme of its kind it exhibits striking genetic similarities with the Type IV enzyme McrBC, which is specific for methylated DNA. The presence of an McrBC-like enzyme as a component of an R/M system suggests two mutually exclusive modes of action for these enzyme types whereby they are either in isolation, cleaving methylated DNA (McrBC), or in combination with MTases cleaving unmethylated DNA (LlaJI). Whilst unusual, this is not altogether surprising as such contrasting activities have been described for some conventional endonucleases (DpnI and DpnII) [39, 40]. Similarly organized R/M systems were found in the database and together with LlaJI, are likely to represent a new subgroup of the Type II family, combining features of both Type II and Type IV systems. Further biochemical characterization of LlaJI, and indeed its identified and relatively unknown homologues (Fig 1A) will hopefully disclose further revelations regarding this novel and remarkable endonuclease family.
Bacteria, bacteriophage, plasmids, media and growth conditions
Details of bacterial strains, bacteriophage and plasmids used in this study are summarized in Table 1. All L. lactis strains were grown in M17 broth (Oxoid ltd., Hampshire, UK) containing 0.5 % glucose at 30°C. E. coli was grown at 37°C in Luria-Bertani media . Where appropriate, antibiotics were added as follows: for L. lactis, 5 μg ml-1 tetracycline, and 10 μg ml-1 chloramphenicol. For E. coli, 100 μg ml-1 ampicilin, 25 μg ml-1kanamycin and 10 μg ml-1 chloramphenicol. Induction of the nisin-inducible promoter of pNZ8020 in the NZ9000 host background was achieved by inclusion of nisin (1 ng ml-1) in growth media . Induction of Ptac, T5, and Plac promoters of pMal-c2x, pQE60, and pUC19, respectively, was achieved by inclusion of 1 mM IPTG in growth media of their respective E.coli hosts. Lactococcal bacteriophage were propagated as described previously . Plaque assays were conducted as described elsewhere  and the efficiency of plaquing (EOP) was calculated as the ratio of the number of plaques formed on a tested host to those formed on a sensitive host. Bacteriophage λ.0 was propagated and enumerated using standard techniques .
Restriction enzymes and T4 DNA ligase were purchased from New England Biolabs (Hertfordshire, UK), and used according to the manufacturer's instructions. E.coli plasmid DNA was isolated using the SV wizard plasmid miniprep kit (Promega, Madison, WI, USA) according to the manufacturer's instructions. Lactococcal plasmid DNA was isolated as described previously . Electroporation of plasmid DNA into E.coli was performed using standard techniques  and in L. lactis as described previously . For routine PCR applications Taq DNA polymerase (Qiagen, West Sussex, UK.) was used. For plasmid construction, where high fidelity was required, PCR's were conducted using Phusion DNA polymerase (New England Biolabs) according to the manufacturer's instructions. All PCR products were purified using the Jetquick PCR purification kit (Genomed, Lohne, Germany). The correct orientation and integrity of all constructs was verified by sequencing, performed at MWG Biotech (Ebersberg, Germany). Sequence assembly, verification and analysis were achieved using the seqman program from the DNASTAR package (DNASTAR, Madison, WI.). Database searches were performed using the BLAST suite of programs  and the REBASE database . The nucleotide sequence of the LlaJI system is present in the GenBank database; accession no. AY530537.
For the construction of plasmid pNZ-R1, an llaJIR1 encompassing PCR product (coordinates 2999–4782, using primers R1F 5'gggggatcc catataggaggacaaaataatg'3 and R1R 5'gggctcgag ccatcctctattcttcagtg'3) was inserted into the BamHI – XhoI sites of pNZ8020. Plasmid pNZ-R2 was constructed by insertion of an llaJIR2 encompassing PCR product (coordinates 4769–6031, using primers R2F 5'gggggatcc gaatagaggatggctctatg'3 and R2R 5'gggctcgag ctaggaattcgggtattctat'3) into the BamHI – XhoI sites of pNZ8020. Similarly, for plasmid pNZ-R1R2, an llaJIR1R2 encompassing PCR product (coordinates 2999–6031, using primers R1F and R2R) was inserted into the BamHI – XhoI sites of pNZ8020.
For plasmid pQE-R1, an llaJIR1 encompassing PCR product (coordinates 3017–4772, using primers R1QEF 5'gggccatgg gaaaaatcgtattacca'3 and R1QER 5'gggagatct ttcttcagtgctaacattatc'3) was inserted into the NcoI – BglII sites of pQE60. For plasmid pQE-R2, an llaJIR2 encompassing PCR product (coordinates 4786–6031, using primers R2QEF 5'gggccatgg aaactaacaaaactattaa'3 and R2QER5'gggagatct ggaattcgggtattctattg) was inserted into the NcoI – BglII sites of pQE60. Plasmid pMal-R2 was constructed by insertion of an llaJIR2 encompassing PCR product (coordinates 4786–6031 using primers MalR2F 5'atggaaactaacaaaactattaa'3 and MalR2R 5'gggtctaga ctaggaattcgggtattcta'3) into the XmnI – XbaI sites of pMal-c2x. Plasmid pUC-R1 was constructed by insertion of an llaJIR1 encompassing PCR product (coordinates 3003–4775, using primers UCR1F 5'gggggtcgac ttaaggaggacaaaataatgggaa'3 and UCR1R 5'gggggatcc ctattcttcagtgctaacatta'3) into the SalI – BamHI site of pUC19. For plasmid pUC-R2, an llaJIR2 encompassing PCR product (coordinates 4773–6031, using primers UCR2F 5' gggggtcgac ttaagaggatggctctatggaaa'3 and UCR2R 5'ggggaatcc ctaggaattcgggtattctat'3) was inserted into the SalI – BamHI sites of pUC19. Plasmid pUC-R1R1 was constructed by insertion of an llaJIR1R2 encompassing PCR product (coordinates 3003–6031, using primers UCR1F and UCR2R) into the SalI – BamHI sites of pUC19.
Protein expression and purification
R1.LlaJI was expressed as a C-terminal 6xhis tagged derivative, overproduced from plasmid pQE-R1 in the M15/pRep4 host background and purified essentially to homogeneity using Ni-NTA Agarose (Qiagen) affinity chromatography, according to the manufacturer's instructions. R2.LlaJI was overproduced from plasmid pMal-R2 in ER2683 and purified by affinity chromatography using amlyose resin (New England Biolabs), according to the manufacturer's instructions. Purified proteins were subsequently dialyzed against buffer containing 10 mM Tris (pH 7.5), 50 mM KCl, 1 mM DTT, 0.1 mM EDTA and 50 % glycerol, and stored at -20°C. Initial mass and purity estimations were achieved by SDS-PAGE  under denaturing conditions. Mass spectrometry and N-terminal amino acid sequencing where indicated was performed at the University of Newcastle Upon Tyne using an Applied Biosystems Voyager DESTR mass spectrometer and a Beckman LF 3000 microsequencer.
Gel mobility shift analysis
DNA binding assays were carried out in 20 μl volumes containing 50 mM Tris-HCl (pH 7.5), 150 mM NaCl, 5 mM β-mercaptoethanol, 8 % glycerol and 50 μg ml-1 calf thymus DNA and approximately 2.5–5 nM final concentration of a probe. All probes were generated by PCR amplification of a section of the LlaJI promoter (coordinates 309–539 using primers akjF 5'ctggagactatattgctcatcg'3 and P1R 5'gggattgtgtttgaaatcacacc'3) which contained two LlaJI recognition sites. Probe J represented the "wild type" sequence. Disruption of the 5' proximal and 3'proximal sites individually (by C-A replacement of the internal cytosine of 5'GCGTC'3) to produce probes J1 and J2 respectively, was achieved by oligonucleotide-directed mutagenesis. Both sites were disrupted in this manner to produce probe J3, which was devoid of any LlaJI recognition sites. All probes were subsequently end-labelled with (γ-32P)ATP by use of polynucleotide kinase (New England Biolabs). Binding reactions were allowed to proceed at room temperature for 30 minutes, followed by the addition of 5 μl 50 % glycerol. Samples were then subjected to electrophoresis on a 4 % polyacrylamide gel containing 2.5 % glycerol in 1 × TAE buffer for 3 hours at 120 V. The protein-DNA complexes were detected by autoradiography of the dried gel.
DNA Cleavage assays
Crude cell extracts of L. lactis strains were prepared from 1 l of culture grown in GM17 broth, concentrated by centrifugation (7000 rpm, 15 min) and washed in lysis buffer (10 mM Tris-HCl (pH 7.6), 1 mM EDTA, 100 mM NaCl and 1 mM DTT). The cell pellets were resuspended in 10 ml of lysis buffer and disrupted with a French press (at 2000 lb/in2) after which cell debris and ribosomes were removed by centrifugation (14,000 rpm, 20 min). E. coli crude extracts were generated essentially as above from 1 l of culture grown in LB broth, where lysis was achieved by sonication (5 × 10 sec bursts punctuated by 30 sec on ice). Cleavage activity of LlaJI was surveyed in all 4 NEBuffers (New England Biolabs) and at temperatures ranging from 19 – 30°C.
Typically the reactions contained 1 μg substrate DNA (λ DNA and Sk1 DNA for cleavage assays or pBR322 DNA for nicking assays), 100 μg ml-1 BSA, and 2 μl crude extract in addition to 5 mM GTP in a final volume of 20 μl. For assays with purified proteins, typically reactions contained 0.5 μg R1.LlaJI or R2.LlaJI (or combinations of the two at ratios of 5:1, 2:1, 1:1, 1:2, 1:5, R1.LlaJI:R2.LlaJI).
GTP hydrolysis activity of R1.LlaJI was investigated via quantification of the release of inorganic phosphate by colorimetric assay based on spectrophotometric quantification of a phosphomolybdate-malachite green complex using the Malachite Green Phosphate Assay kit (Bioassay Systems, Hayward, CA 94545, USA) according to the manufacturer's instructions.
This work was funded by Science Foundation Ireland (02/IN1/B198).
We would like to thank Michiel Kleerebezem for supplying plasmid pNZ8020, and L. lactis NZ9000, Todd Klaenhammer for supplying L. lactis NCK690KP1 and Sylvain Moineau for phage Ul36. We also thank Shawn Stickel, Junichi Ryu, Yu Zheng, and Stephen McGrath for helpful discussion.
- O'Driscoll J, Glynn F, Cahalane O, O'Connell-Motherway M, Fitzgerald GF, Van Sinderen D: Lactococcal Plasmid pNP40 Encodes a Novel, Temperature-Sensitive Restriction-Modification System. Appl Environ Microbiol. 2004, 70 (9): 5546-5556. 10.1128/AEM.70.9.5546-5556.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Kunst F, Ogasawara N, Moszer I, Albertini AM, Alloni G, Azevedo V, Bertero MG, Bessieres P, Bolotin A, Borchert S, Borriss R, Boursier L, Brans A, Braun M, Brignell SC, Bron S, Brouillet S, Bruschi CV, Caldwell B, Capuano V, Carter NM, Choi SK, Codani JJ, Connerton IF, Cummings NJ, Daniel RA, Denizot F, Devine KM, Dusterhoft A, Ehrlich SD, Emmerson PT, Entian KD, Errington J, Fabret C, Ferrari E, Foulger D, Fritz C, Fujita M, Fujita Y, Fuma S, Galizzi A, Galleron N, Ghim SY, Glaser P, Goffeau A, Golightly EJ, Grandi G, Guiseppi G, Guy BJ, Haga K, Haiech J, Harwood CR, Henaut A, Hilbert H, Holsappel S, Hosono S, Hullo MF, Itaya M, Jones L, Joris B, Karamata D, Kasahara Y, Klaerr-Blanchard M, Klein C, Kobayashi Y, Koetter P, Koningstein G, Krogh S, Kumano M, Kurita K, Lapidus A, Lardinois S, Lauber J, Lazarevic V, Lee SM, Levine A, Liu H, Masuda S, Mauel C, Medigue C, Medina N, Mellado RP, Mizuno M, Moestl D, Nakai S, Noback M, Noone D, O'Reilly M, Ogawa K, Ogiwara A, Oudega B, Park SH, Parro V, Pohl TM, Portetelle D, Porwollik S, Prescott AM, Presecan E, Pujic P, Purnelle B, Rapoport G, Rey M, Reynolds S, Rieger M, Rivolta C, Rocha E, Roche B, Rose M, Sadaie Y, Sato T, Scanlan E, Schleich S, Schroeter R, Scoffone F, Sekiguchi J, Sekowska A, Seror SJ, Serror P, Shin BS, Soldo B, Sorokin A, Tacconi E, Takagi T, Takahashi H, Takemaru K, Takeuchi M, Tamakoshi A, Tanaka T, Terpstra P, Tognoni A, Tosato V, Uchiyama S, Vandenbol M, Vannier F, Vassarotti A, Viari A, Wambutt R, Wedler E, Wedler H, Weitzenegger T, Winters P, Wipat A, Yamamoto H, Yamane K, Yasumoto K, Yata K, Yoshida K, Yoshikawa HF, Zumstein E, Yoshikawa H, Danchin A: The complete genome sequence of the Gram-positive bacterium Bacillus subtilis. Nature. 1997, 390 (6657): 249-256. 10.1038/36786.View ArticlePubMedGoogle Scholar
- O'Sullivan DJ, Zagula K, Klaenhammer TR: In vivo restriction by LlaI is encoded by three genes, arranged in an operon with llaIM, on the conjugative Lactococcus plasmid pTR2030. J Bacteriol. 1995, 177 (1): 134-143.PubMed CentralPubMedGoogle Scholar
- Ivanova N, Sorokin A, Anderson I, Galleron N, Candelon B, Kapatral V, Bhattacharyya A, Reznik G, Mikhailova N, Lapidus A, Chu L, Mazur M, Goltsman E, Larsen N, D'Souza M, Walunas T, Grechkin Y, Pusch G, Haselkorn R, Fonstein M, Ehrlich SD, Overbeek R, Kyrpides N: Genome sequence of Bacillus cereus and comparative analysis with Bacillus anthracis. Nature. 2003, 423 (6935): 87-91. 10.1038/nature01582.View ArticlePubMedGoogle Scholar
- Alm RA, Ling LS, Moir DT, King BL, Brown ED, Doig PC, Smith DR, Noonan B, Guild BC, deJonge BL, Carmel G, Tummino PJ, Caruso A, Uria-Nickelsen M, Mills DM, Ives C, Gibson R, Merberg D, Mills SD, Jiang Q, Taylor DE, Vovis GF, Trust TJ: Genomic-sequence comparison of two unrelated isolates of the human gastric pathogen Helicobacter pylori. Nature. 1999, 397 (6715): 176-180. 10.1038/16495.View ArticlePubMedGoogle Scholar
- Banks DJ, Porcella SF, Barbian KD, Beres SB, Philips LE, Voyich JM, DeLeo FR, Martin JM, Somerville GA, Musser JM: Progress toward characterization of the group A Streptococcus metagenome: complete genome sequence of a macrolide-resistant serotype M6 strain. J Infect Dis. 2004, 190 (4): 727-738. 10.1086/422697.View ArticlePubMedGoogle Scholar
- Roberts RJ, Vincze T, Posfai J, Macelis D: REBASE – restriction enzymes and DNA methyltransferases. Nucleic Acids Res. 2005, 33 (Database): D230-232. 10.1093/nar/gki029.PubMed CentralPubMedGoogle Scholar
- Roberts RJ, Belfort M, Bestor T, Bhagwat AS, Bickle TA, Bitinaite J, Blumenthal RM, Degtyarev S, Dryden DT, Dybvig K, Firman K, Gromova ES, Gumport RI, Halford SE, Hattman S, Heitman J, Hornby DP, Janulaitis A, Jeltsch A, Josephsen J, Kiss A, Klaenhammer TR, Kobayashi I, Kong H, Kruger DH, Lacks S, Marinus MG, Miyahara M, Morgan RD, Murray NE, Nagaraja V, Piekarowicz A, Pingoud A, Raleigh E, Rao DN, Reich N, Repin VE, Selker EU, Shaw PC, Stein DC, Stoddard BL, Szybalski W, Trautner TA, Van Etten JL, Vitor JM, Wilson GG, Xu SY: 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 CentralView ArticlePubMedGoogle Scholar
- Pingoud A, Jeltsch A: Structure and function of type II restriction endonucleases. Nucleic Acids Res. 2001, 29 (18): 3705-3727. 10.1093/nar/29.18.3705.PubMed CentralView ArticlePubMedGoogle Scholar
- Pieper U, Pingoud A: A mutational analysis of the PD...D/EXK motif suggests that McrC harbors the catalytic center for DNA cleavage by the GTP-dependent restriction enzyme McrBC from Escherichia coli. Biochemistry. 2002, 41 (16): 5236-5244. 10.1021/bi0156862.View ArticlePubMedGoogle Scholar
- Pingoud A, Jeltsch A: Recognition and cleavage of DNA by type-II restriction endonucleases. Eur J Biochem. 1997, 246 (1): 1-22. 10.1111/j.1432-1033.1997.t01-6-00001.x.View ArticlePubMedGoogle Scholar
- Lagunavicius A, Sasnauskas G, Halford SE, Siksnys V: The metal-independent type IIs restriction enzyme BfiI is a dimer that binds two DNA sites but has only one catalytic centre. J Mol Biol. 2003, 326 (4): 1051-1064. 10.1016/S0022-2836(03)00020-2.View ArticlePubMedGoogle Scholar
- Sasnauskas G, Halford SE, Siksnys V: How the BfiI restriction enzyme uses one active site to cut two DNA strands. Proc Natl Acad Sci U S A. 2003, 100 (11): 6410-6415. 10.1073/pnas.1131003100.PubMed CentralView ArticlePubMedGoogle Scholar
- Halford SE, Bilcock DT, Stanford NP, Williams SA, Milsom SE, Gormley NA, Watson MA, Bath AJ, Embleton ML, Gowers DM, Daniels LE, Parry SH, Szczelkun MD: Restriction endonuclease reactions requiring two recognition sites. Biochem Soc Trans. 1999, 27 (4): 696-699.View ArticlePubMedGoogle Scholar
- Heiter DF, Lunnen KD, Wilson GG: Site-specific DNA-nicking mutants of the heterodimeric restriction endonuclease R.BbvCI. J Mol Biol. 2005, 348 (3): 631-640. 10.1016/j.jmb.2005.02.034.View ArticlePubMedGoogle Scholar
- Hsieh PC, Xiao JP, O'Loane D, Xu SY: Cloning, expression, and purification of a thermostable nonhomodimeric restriction enzyme, BslI. J Bacteriol. 2000, 182 (4): 949-955. 10.1128/JB.182.4.949-955.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Stankevicius K, Lubys A, Timinskas A, Vaitkevicius D, Janulaitis A: Cloning and analysis of the four genes coding for Bpu10I restriction-modification enzymes. Nucleic Acids Res. 1998, 26 (4): 1084-1091. 10.1093/nar/26.4.1084.PubMed CentralView ArticlePubMedGoogle Scholar
- Ohshima H, Matsuoka S, Asai K, Sadaie Y: Molecular organization of intrinsic restriction and modification genes BsuM of Bacillus subtilis Marburg. J Bacteriol. 2002, 184 (2): 381-389. 10.1128/JB.184.2.381-389.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Bellamy SR, Milsom SE, Scott DJ, Daniels LE, Wilson GG, Halford SE: Cleavage of individual DNA strands by the different subunits of the heterodimeric restriction endonuclease BbvCI. J Mol Biol. 2005, 348 (3): 641-653. 10.1016/j.jmb.2005.02.035.View ArticlePubMedGoogle Scholar
- O'Driscoll J, Fitzgerald GF, van Sinderen D: A dichotomous epigenetic mechanism governs expression of the LlaJI restriction/modification system. Mol Microbiol. 2005, 57 (6): 1532-1544. 10.1111/j.1365-2958.2005.04769.x.View ArticlePubMedGoogle Scholar
- Raleigh EA: Organization and function of the mcrBC genes of Escherichia coli K-12. Mol Microbiol. 1992, 6 (9): 1079-1086.View ArticlePubMedGoogle Scholar
- Stewart FJ, Panne D, Bickle TA, Raleigh EA: Methyl-specific DNA binding by McrBC, a modification-dependent restriction enzyme. J Mol Biol. 2000, 298 (4): 611-622. 10.1006/jmbi.2000.3697.View ArticlePubMedGoogle Scholar
- Panne D, Raleigh EA, Bickle TA: The McrBC endonuclease translocates DNA in a reaction dependent on GTP hydrolysis. J Mol Biol. 1999, 290 (1): 49-60. 10.1006/jmbi.1999.2894.View ArticlePubMedGoogle Scholar
- Stewart FJ, Raleigh EA: Dependence of McrBC cleavage on distance between recognition elements. Biol Chem. 1998, 379 (4–5): 611-616.PubMedGoogle Scholar
- Sutherland E, Coe L, Raleigh EA: McrBC: a multisubunit GTP-dependent restriction endonuclease. J Mol Biol. 1992, 225 (2): 327-348. 10.1016/0022-2836(92)90925-A.View ArticlePubMedGoogle Scholar
- Dila D, Sutherland E, Moran L, Slatko B, Raleigh EA: Genetic and sequence organization of the mcrBC locus of Escherichia coli K-12. J Bacteriol. 1990, 172 (9): 4888-4900.PubMed CentralPubMedGoogle Scholar
- Kruger T, Grund C, Wild C, Noyer-Weidner M: Characterization of the mcrBC region of Escherichia coli K-12 wild-type and mutant strains. Gene. 1992, 114 (1): 1-12. 10.1016/0378-1119(92)90700-Y.View ArticlePubMedGoogle Scholar
- Pieper U, Brinkmann T, Kruger T, Noyer-Weidner M, Pingoud A: Characterization of the interaction between the restriction endonuclease McrBC from E. coli and its cofactor GTP. J Mol Biol. 1997, 272 (2): 190-199. 10.1006/jmbi.1997.1228.View ArticlePubMedGoogle Scholar
- Pieper U, Schweitzer T, Groll DH, Gast FU, Pingoud A: The GTP-binding domain of McrB: more than just a variation on a common theme?. J Mol Biol. 1999, 292 (3): 547-556. 10.1006/jmbi.1999.3103.View ArticlePubMedGoogle Scholar
- Gast FU, Brinkmann T, Pieper U, Kruger T, Noyer-Weidner M, Pingoud A: The recognition of methylated DNA by the GTP-dependent restriction endonuclease McrBC resides in the N-terminal domain of McrB. Biol Chem. 1997, 378 (9): 975-982.View ArticlePubMedGoogle Scholar
- Panne D, Muller SA, Wirtz S, Engel A, Bickle TA: The McrBC restriction endonuclease assembles into a ring structure in the presence of G nucleotides. Embo J. 2001, 20 (12): 3210-3217. 10.1093/emboj/20.12.3210.PubMed CentralView ArticlePubMedGoogle Scholar
- Raleigh EA, Wilson G: Escherichia coli K-12 restricts DNA containing 5-methylcytosine. Proc Natl Acad Sci U S A. 1986, 83 (23): 9070-9074.PubMed CentralView ArticlePubMedGoogle Scholar
- Panne D, Raleigh EA, Bickle TA: McrBs, a modulator peptide for McrBC activity. Embo J. 1998, 17 (18): 5477-5483. 10.1093/emboj/17.18.5477.PubMed CentralView ArticlePubMedGoogle Scholar
- Kruger T, Wild C, Noyer-Weidner M: McrB: a prokaryotic protein specifically recognizing DNA containing modified cytosine residues. Embo J. 1995, 14 (11): 2661-2669.PubMed CentralPubMedGoogle Scholar
- Bourniquel AA, Bickle TA: Complex restriction enzymes: NTP-driven molecular motors. Biochimie. 2002, 84 (11): 1047-1059. 10.1016/S0300-9084(02)00020-2.View ArticlePubMedGoogle Scholar
- Saraste M, Sibbald PR, Wittinghofer A: The P-loop – a common motif in ATP- and GTP-binding proteins. Trends Biochem Sci. 1990, 15 (11): 430-434. 10.1016/0968-0004(90)90281-F.View ArticlePubMedGoogle Scholar
- Pieper U, Schweitzer T, Groll DH, Pingoud A: Defining the location and function of domains of McrB by deletion mutagenesis. Biol Chem. 1999, 380 (10): 1225-1230. 10.1515/BC.1999.155.View ArticlePubMedGoogle Scholar
- Kjeldgaard M, Nyborg J, Clark BF: The GTP binding motif: variations on a theme. Faseb J. 1996, 10 (12): 1347-1368.PubMedGoogle Scholar
- Lacks S, Greenberg B: Complementary specificity of restriction endonucleases of Diplococcus pneumoniae with respect to DNA methylation. J Mol Biol. 1977, 114 (1): 153-168. 10.1016/0022-2836(77)90289-3.View ArticlePubMedGoogle Scholar
- Lacks S, Greenberg B: A deoxyribonuclease of Diplococcus pneumoniae specific for methylated DNA. J Biol Chem. 1975, 250 (11): 4060-4066.PubMedGoogle Scholar
- Sambrook J, Fritsch EF, Maniatis T: Molecular cloning : a laboratory manual. 1989, Cold Spring Harbor, N.Y.: Cold Spring Harbor Laboratory, 2Google Scholar
- Kuipers OP, de Ruyter PGGA, Kleerebezem M, de Vos WM: Quorum sensing-controlled gene expression in lactic acid bacteria. Journal of Biotechnology. 1998, 64 (1): 15-21. 10.1016/S0168-1656(98)00100-X.View ArticleGoogle Scholar
- Lillehaug D: An improved plaque assay for poor plaque-producing temperate lactococcal bacteriophages. J Appl Microbiol. 1997, 83 (1): 85-90. 10.1046/j.1365-2672.1997.00193.x.View ArticlePubMedGoogle Scholar
- O'Sullivan D, Klaenhammer T: Rapid Mini-Prep Isolation of High-Quality Plasmid DNA from Lactococcus and Lactobacillus spp. Appl Environ Microbiol. 1993, 59 (8): 2730-2733.PubMed CentralPubMedGoogle Scholar
- Wells JM, Wilson PW, Le Page RW: Improved cloning vectors and transformation procedure for Lactococcus lactis. J Appl Bacteriol. 1993, 74 (6): 629-636.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Laemmli UK: Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature. 1970, 227 (5259): 680-685. 10.1038/227680a0.View ArticlePubMedGoogle Scholar
- Law J, Buist G, Haandrikman A, Kok J, Venema G, Leenhouts K: A system to generate chromosomal mutations in Lactococcus lactis which allows fast analysis of targeted genes. J Bacteriol. 1995, 177 (24): 7011-7018.PubMed CentralPubMedGoogle Scholar
- Gasson MJ: Plasmid complements of Streptococcus lactis NCDO 712 and other lactic streptococci after protoplast-induced curing. J Bacteriol. 1983, 154 (1): 1-9.PubMed CentralPubMedGoogle Scholar
- Djordjevic GM, Klaenhammer TR: Bacteriophage-triggered defense systems: phage adaptation and design improvements. Appl Environ Microbiol. 1997, 63 (11): 4370-4376.PubMed CentralPubMedGoogle Scholar
- Chandry PS, Davidson BE, Hillier AJ: Temporal transcription map of the Lactococcus lactis bacteriophage sk1. Microbiology. 1994, 140 (Pt 9): 2251-2261.View ArticlePubMedGoogle Scholar
- Higgins DL, Sanozky-Dawes RB, Klaenhammer TR: Restriction and modification activities from Streptococcus lactis ME2 are encoded by a self-transmissible plasmid, pTN20, that forms cointegrates during mobilization of lactose-fermenting ability. J Bacteriol. 1988, 170 (8): 3435-3442.PubMed CentralPubMedGoogle Scholar
- Labrie S, Moineau S: Complete genomic sequence of bacteriophage ul36: demonstration of phage heterogeneity within the P335 quasi-species of lactococcal phages. Virology. 2002, 296 (2): 308-320. 10.1006/viro.2002.1401.View ArticlePubMedGoogle Scholar
- de Ruyter PG, Kuipers OP, de Vos WM: Controlled gene expression systems for Lactococcus lactis with the food-grade inducer nisin. Appl Environ Microbiol. 1996, 62 (10): 3662-3667.PubMed CentralPubMedGoogle Scholar
- Yanisch-Perron C, Vieira J, Messing J: Improved M13 phage cloning vectors and host strains: nucleotide sequences of the M13mp18 and pUC19 vectors. Gene. 1985, 33 (1): 103-119. 10.1016/0378-1119(85)90120-9.View ArticlePubMedGoogle Scholar
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.