Purification and characterization of C.BamHI, a regulator of the BamHI restriction-modification system

Controller protein of restriction-modification system Kpn2I affects transcription of its gene by acting as a transcription elongation roadblock

C-proteins control restriction-modification (R-M) systems' genes transcription to ensure sufficient levels of restriction endonuclease to allow protection from foreign DNA while avoiding its modification by excess methyltransferase. Here, we characterize transcription regulation in C-protein dependent R-M system Kpn2I. The Kpn2I restriction endonuclease gene is transcribed from a constitutive, weak promoter, which, atypically, is C-protein independent. Kpn2I C-protein (C.Kpn2I) binds upstream of the strong methyltransferase gene promoter and inhibits it, likely by preventing the interaction of the RNA polymerase sigma subunit with the -35 consensus element. Diminished transcription from the methyltransferase promoter increases transcription from overlapping divergent C-protein gene promoters. All known C-proteins affect transcription initiation from R-M genes promoters. Uniquely, the C.Kpn2I binding site is located within the coding region of its gene. C.Kpn2I acts as a roadblock stalling elongating RNA polymerase and decreasing production of full-length C.Kpn2I mRNA. Mathematical modeling shows that this unusual mode of regulation leads to the same dynamics of accumulation of R-M gene transcripts as observed in systems where C-proteins act at transcription initiation stage only. Bioinformatics analyses suggest that transcription regulation through binding of C.Kpn2I-like proteins within the coding regions of their genes may be widespread.

Understanding key features of bacterial restriction-modification systems through quantitative modeling

Background

Restriction-modification (R-M) systems are rudimentary bacterial immune systems. The main components include restriction enzyme (R), which cuts specific unmethylated DNA sequences, and the methyltransferase (M), which protects the same DNA sequences. The expression of R-M system components is considered to be tightly regulated, to ensure successful establishment in a naïve bacterial host. R-M systems are organized in different architectures (convergent or divergent) and are characterized by different features, i.e. binding cooperativities, dissociation constants of dimerization, translation rates, which ensure this tight regulation. It has been proposed that R-M systems should exhibit certain dynamical properties during the system establishment, such as: i) a delayed expression of R with respect to M, ii) fast transition of R from “OFF” to “ON” state, iii) increased stability of the toxic molecule (R) steady-state levels. It is however unclear how different R-M system features and architectures ensure these dynamical properties, particularly since it is hard to address this question experimentally.

Results

To understand design of different R-M systems, we computationally analyze two R-M systems, representative of the subset controlled by small regulators called ‘C proteins’, and differing in having convergent or divergent promoter architecture. We show that, in the convergent system, abolishing any of the characteristic system features adversely affects the dynamical properties outlined above. Moreover, an extreme binding cooperativity, accompanied by a very high dissociation constant of dimerization, observed in the convergent system, but absent from other R-M systems, can be explained in terms of the same properties. Furthermore, we develop the first theoretical model for dynamics of a divergent R-M system, which does not share any of the convergent system features, but has overlapping promoters. We show that i) the system dynamics exhibits the same three dynamical properties, ii) introducing any of the convergent system features to the divergent system actually diminishes these properties.

Conclusions

Our results suggest that different R-M architectures and features may be understood in terms of constraints imposed by few simple dynamical properties of the system, providing a unifying framework for understanding these seemingly diverse systems. We also provided predictions for the perturbed R-M systems dynamics, which may in future be tested through increasingly available experimental techniques, such as re-engineering R-M systems and single-cell experiments.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0377-x) contains supplementary material, which is available to authorized users.

Highlights of the DNA cutters: a short history of the restriction enzymes

In the early 1950's, 'host-controlled variation in bacterial viruses' was reported as a non-hereditary phenomenon: one cycle of viral growth on certain bacterial hosts affected the ability of progeny virus to grow on other hosts by either restricting or enlarging their host range. Unlike mutation, this change was reversible, and one cycle of growth in the previous host returned the virus to its original form. These simple observations heralded the discovery of the endonuclease and methyltransferase activities of what are now termed Type I, II, III and IV DNA restriction-modification systems. The Type II restriction enzymes (e.g. EcoRI) gave rise to recombinant DNA technology that has transformed molecular biology and medicine. This review traces the discovery of restriction enzymes and their continuing impact on molecular biology and medicine.

Tuning the relative affinities for activating and repressing operators of a temporally regulated restriction-modification system

Most type II restriction-modification (R-M) systems produce separate endonuclease (REase) and methyltransferase (MTase) proteins. After R-M genes enter a new cell, MTase activity must appear before REase or the host chromosome will be cleaved. Temporal control of these genes thus has life-or-death consequences. PvuII and some other R-M systems delay endonuclease expression by cotranscribing the REase gene with the upstream gene for an autogenous activator/repressor (C protein). C.PvuII was previously shown to have low levels early, but positive feedback later boosts transcription of the C and REase genes. The MTase is expressed without delay, and protects the host DNA. C.PvuII binds to two sites upstream of its gene: O_L, associated with activation, and O_R, associated with repression. Even when symmetry elements of each operator are made identical, C.PvuII binds preferentially to O_L. In this study, the intra-operator spacers are shown to modulate relative C.PvuII affinity. In light of a recently reported C.Esp1396I-DNA co-crystal structure, in vitro and in vivo effects of altering O_L and O_R spacers were determined. The results suggest that the GACTnnnAGTC consensus is the primary determinant of C.PvuII binding affinity, with intra-operator spacers playing a fine-tuning role that affects mobility of this R-M system.

Systematic prediction of control proteins and their DNA binding sites

We present here the results of a systematic bioinformatics analysis of control (C) proteins, a class of DNA-binding regulators that control time-delayed transcription of their own genes as well as restriction endonuclease genes in many type II restriction-modification systems. More than 290 C protein homologs were identified and DNA-binding sites for ∼70% of new and previously known C proteins were predicted by a combination of phylogenetic footprinting and motif searches in DNA upstream of C protein genes. Additional analysis revealed that a large proportion of C protein genes are translated from leaderless RNA, which may contribute to time-delayed nature of genetic switches operated by these proteins. Analysis of genetic contexts of newly identified C protein genes revealed that they are not exclusively associated with restriction-modification genes; numerous instances of associations with genes originating from mobile genetic elements were observed. These instances might be vestiges of ancient horizontal transfers and indicate that during evolution ancestral restriction-modification system genes were the sites of mobile elements insertions.

Bioinformatic and partial functional analysis of pEspA and pEspB, two plasmids from Exiguobacterium arabatum sp. nov. RFL1109

The complete nucleotide sequences of two plasmids from Exiguobacterium arabatum sp. nov. RFL1109, pEspA (4563bp) and pEspB (38,945bp), have been determined. Five ORFs were identified in the pEspA plasmid, and putative functions were assigned to two of them. Using deletion mapping approach, the Rep-independent replication region of pEspA, which functions in Bacillus subtilis, was localized within a 0.6kb DNA region. Analysis of the pEspB sequence revealed 42 ORFs. From these, function of two genes encoding enzymes of the Lsp1109I restriction-modification system was confirmed experimentally, while putative functions of another 18 ORFs were suggested based on comparative analysis. Three functional regions have been proposed for the pEspB plasmid: the putative conjugative transfer region, the region involved in plasmid replication and maintenance, and the region responsible for transposition of the IS21 family-like transposable elements.

Regulatory circuit based on autogenous activation-repression: roles of C-boxes and spacer sequences in control of the PvuII restriction-modification system

Type II restriction-modification (R-M) systems comprise a restriction endonuclease (REase) and a protective methyltransferase (MTase). After R-M genes enter a new cell, MTase must appear before REase or the chromosome will be cleaved. PvuII and some other R-M systems achieve this delay by cotranscribing the REase gene with the gene for an autogenous transcription activator (the controlling or 'C' protein C.PvuII). This study reveals, through in vivo titration, that C.PvuII is not only an activator but also a repressor for its own gene. In other systems, this type of circuit can result in oscillatory behavior. Despite the use of identical, symmetrical C protein-binding sequences (C-boxes) in the left and right operators, C.PvuII showed higher in vitro affinity for O(L) than for O(R), implicating the spacer sequences in this difference. Mutational analysis associated the repression with O(R), which overlaps the promoter -35 hexamer but is otherwise dispensable for activation. A nonrepressing mutant exhibited poor establishment in new cells. Comparing promoter-operator regions from PvuII and 29 R-M systems controlled by C proteins revealed that the most-highly conserved sequence is the tetranucleotide spacer separating O(L) from O(R). Any changes in that spacer reduced the stability of C.PvuII-operator complexes and abolished activation.

Nature of the promoter activated by C.PvuII, an unusual regulatory protein conserved among restriction-modification systems

A widely distributed family of small regulators, called C proteins, controls a subset of restriction-modification systems. The C proteins studied to date activate transcription of their own genes and that of downstream endonuclease genes; this arrangement appears to delay endonuclease expression relative to that of the protective methyltransferase when the genes enter a new cell. C proteins bind to conserved sequences called C boxes. In the PvuII system, the C boxes have been reported to extend from -23 to +3 relative to the transcription start for the gene for the C protein, an unexpected starting position relative to a bound activator. This study suggests that transcript initiation within the C boxes represents initial, C-independent transcription of pvuIICR. The major C protein-dependent transcript appears to be a leaderless mRNA starting farther downstream, at the initiation codon for the pvuIIC gene. This conclusion is based on nuclease S1 transcript mapping and the effects of a series of nested deletions in the promoter region. Furthermore, replacing the region upstream of the pvuIIC initiation codon with a library of random oligonucleotides, followed by selection for C-dependent transcription, yielded clones having sequences that resemble -10 promoter hexamers. The -35 hexamer of this promoter would lie within the C boxes. However, the spacing between C boxes/-35 and the apparent -10 hexamer can be varied by +/-4 bp with little effect. This suggests that, like some other activator-dependent promoters, PpvuIICR may not require a -35 hexamer. Features of this transcription activation system suggest explanations for its broad host range.

A nomenclature for restriction enzymes, DNA methyltransferases, homing endonucleases and their genes

A nomenclature is described for restriction endonucleases, DNA methyltransferases, homing endonucleases and related genes and gene products. It provides explicit categories for the many different Type II enzymes now identified and provides a system for naming the putative genes found by sequence analysis of microbial genomes.

Behavior of restriction-modification systems as selfish mobile elements and their impact on genome evolution

Restriction-modification (RM) systems are composed of genes that encode a restriction enzyme and a modification methylase. RM systems sometimes behave as discrete units of life, like viruses and transposons. RM complexes attack invading DNA that has not been properly modified and thus may serve as a tool of defense for bacterial cells. However, any threat to their maintenance, such as a challenge by a competing genetic element (an incompatible plasmid or an allelic homologous stretch of DNA, for example) can lead to cell death through restriction breakage in the genome. This post-segregational or post-disturbance cell killing may provide the RM complexes (and any DNA linked with them) with a competitive advantage. There is evidence that they have undergone extensive horizontal transfer between genomes, as inferred from their sequence homology, codon usage bias and GC content difference. They are often linked with mobile genetic elements such as plasmids, viruses, transposons and integrons. The comparison of closely related bacterial genomes also suggests that, at times, RM genes themselves behave as mobile elements and cause genome rearrangements. Indeed some bacterial genomes that survived post-disturbance attack by an RM gene complex in the laboratory have experienced genome rearrangements. The avoidance of some restriction sites by bacterial genomes may result from selection by past restriction attacks. Both bacteriophages and bacteria also appear to use homologous recombination to cope with the selfish behavior of RM systems. RM systems compete with each other in several ways. One is competition for recognition sequences in post-segregational killing. Another is super-infection exclusion, that is, the killing of the cell carrying an RM system when it is infected with another RM system of the same regulatory specificity but of a different sequence specificity. The capacity of RM systems to act as selfish, mobile genetic elements may underlie the structure and function of RM enzymes.

Transcriptional analysis and regulation of expression of the ScrFI restriction-modification system of Lactococcus lactis subsp. cremoris UC503

ScrFI is a type II restriction-modification system from Lactococcus lactis which recognizes the nucleotide sequence 5'-CC downward arrow NGG-3', cleaving at the point indicated by the arrow, and it comprises an endonuclease gene that is flanked on either side by genes encoding two 5-methylcytosine methylases. An open reading frame (orfX) of unknown function is located immediately upstream of these genes. In this study Northern analysis was performed, and it revealed that orfX, scrFIBM, and scrFIR are cotranscribed as a single polygenic mRNA molecule, while scrFIAM is transcribed independently. 5' extension analysis indicated that the start site for the scrFIAM promoter was a thymine located 4 bp downstream of the -10 motif. The transcriptional start site for the orfX promoter was also found to be a thymine which is more atypically located 24 bp downstream of the -10 motif proximal to the start codon. A helix-turn-helix motif was identified at the N-terminal end of one of the methylases (M.ScrFIA). In order to determine if this motif played a role in regulation of the ScrFI locus, M.ScrFIA was purified. It was then employed in gel retardation assays using fragments containing the two promoters found on the ScrFI operon, one located upstream of orfX and the other located just upstream of scrFIAM. M.ScrFIA was found to bind to the promoter region upstream of the gene encoding it, indicating that it may have a regulatory role. In further studies the two putative promoters were introduced into a vector (pAK80) upstream of a promoterless lacZ gene, and cloned fragments of the ScrFI locus were introduced in trans with each of these promoter constructs to investigate the effect on promoter activity. These results implicated M.ScrFIA in regulation of both promoters on the ScrFI locus.

DNA methylation at the CfrBI site is involved in expression control in the CfrBI restriction-modification system

We have previously found that genes of the CFR:BI restriction-modification (R-M) system from Citrobacter freundii are oriented divergently and that their promoter regions overlap. The overlapping promoters suggest regulation of gene expression at the transcriptional level. In this study the transcription regulation of CFR:BI R-M genes was analyzed in vivo and in vitro in Escherichia coli. It was shown that in the presence of CFR:BI methyltransferase (M.CFR:BI), cell galactokinase activity decreases 10-fold when the galactokinase gene (galK) is under the control of the cfrBIM promoter and increases 20-fold when galK is under the control of the cfrBIR promoter. The CFR:BI site, proven to be unique for the entire CFR:BI R-M gene sequence, is located in the -35 cfrBIM promoter region and is in close vicinity of the -10 cfrBIR promoter region. A comparison of the cfrBIM and the cfrBIR promoter activities in the in vitro transcription system using methylated and unmethylated DNA fragments as templates demonstrated that the efficiency of CFR:BI R-M gene transcription is regulated by enzymatic modification at the N-4-position of cytosine bases of the CFR:BI site by M.CFR:BI. From the results of the in vivo and in vitro experiments we suggest a new model of gene expression regulation in type II R-M systems.

Recognition of native DNA methylation by the PvuII restriction endonuclease

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