Avoidance of recognition sites of restriction-modification systems is a widespread but not universal anti-restriction strategy of prokaryotic viruses.

BACKGROUND: Restriction-modification (R-M) systems protect bacteria and archaea from attacks by bacteriophages and archaeal viruses. An R-M system specifically recognizes short sites in foreign DNA and cleaves it, while such sites in the host DNA are protected by methylation. Prokaryotic viruses have developed a number of strategies to overcome this host defense. The simplest anti-restriction strategy is the elimination of recognition sites in the viral genome: no sites, no DNA cleavage. Even a decrease of the number of recognition sites can help a virus to overcome this type of host defense. Recognition site avoidance has been a known anti-restriction strategy of prokaryotic viruses for decades. However, recognition site avoidance has not been systematically studied with the currently available sequence data. We analyzed the complete genomes of almost 4000 prokaryotic viruses with known host species and more than 17,000 restriction endonucleases with known specificities in terms of recognition site avoidance. RESULTS: We observed considerable limitations of recognition site avoidance as an anti-restriction strategy. Namely, the avoidance of recognition sites is specific for dsDNA and ssDNA prokaryotic viruses. Avoidance is much more pronounced in the genomes of non-temperate bacteriophages than in the genomes of temperate ones. Avoidance is not observed for the sites of Type I and Type IIG systems and is very rarely observed for the sites of Type III systems. The vast majority of avoidance cases concern recognition sites of orthodox Type II restriction-modification systems. Even under these constraints, complete or almost complete elimination of sites is observed for approximately one-tenth of viral genomes and a significant under-representation for approximately one-fourth of them. CONCLUSIONS: Avoidance of recognition sites of restriction-modification systems is a widespread but not universal anti-restriction strategy of prokaryotic viruses.

Comparison of Methods of Detection of Exceptional Sequences in Prokaryotic Genomes.

Many proteins need recognition of specific DNA sequences for functioning. The number of recognition sites and their distribution along the DNA might be of biological importance. For example, the number of restriction sites is often reduced in prokaryotic and phage genomes to decrease the probability of DNA cleavage by restriction endonucleases. We call a sequence an exceptional one if its frequency in a genome significantly differs from one predicted by some mathematical model. An exceptional sequence could be either under- or over-represented, depending on its frequency in comparison with the predicted one. Exceptional sequences could be considered biologically meaningful, for example, as targets of DNA-binding proteins or as parts of abundant repetitive elements. Several methods to predict frequency of a short sequence in a genome, based on actual frequencies of certain its subsequences, are used. The most popular are methods based on Markov chain models. But any rigorous comparison of the methods has not previously been performed. We compared three methods for the prediction of short sequence frequencies: the maximum-order Markov chain model-based method, the method that uses geometric mean of extended Markovian estimates, and the method that utilizes frequencies of all subsequences including discontiguous ones. We applied them to restriction sites in complete genomes of 2500 prokaryotic species and demonstrated that the results depend greatly on the method used: lists of 5% of the most under-represented sites differed by up to 50%. The method designed by Burge and coauthors in 1992, which utilizes all subsequences of the sequence, showed a higher precision than the other two methods both on prokaryotic genomes and randomly generated sequences after computational imitation of selective pressure. We propose this method as the first choice for detection of exceptional sequences in prokaryotic genomes.

Endonuclease Activity of MutL Protein of the Rhodobacter sphaeroides Mismatch Repair System.

We have purified the MutL protein from Rhodobacter sphaeroides mismatch repair system (rsMutL) for the first time. rsMutL demonstrated endonuclease activity in vitro, as predicted by bioinformatics analysis. Based on the alignment of 1483 sequences of bacterial MutL homologs with presumed endonuclease activity, conserved functional motifs and amino acid residues in the rsMutL sequence were identified: five motifs comprising the catalytic site responsible for DNA cleavage were found in the C-terminal domain; seven conserved motifs involved in ATP binding and hydrolysis and specific to the GHKL family of ATPases were found in the N-terminal domain. rsMutL demonstrated the highest activity in the presence of Mn2+. The extent of plasmid DNA hydrolysis declined in the row Mn2+ > Co2+ > Mg2+ > Cd2+; Ni2+ and Ca2+ did not activate rsMutL. Divalent zinc ions inhibited rsMutL endonuclease activity in the presence of Mn2+ excess. ATP also suppressed plasmid DNA hydrolysis by rsMutL. Analysis of amino acid sequences and biochemical properties of five studied bacterial MutL homologs with endonuclease activity revealed that rsMutL resembles the MutL proteins from Neisseria gonorrhoeae and Pseudomonas aeruginosa.

A New Program for Detecting the Geometrical Core of a Set of Structures of Macromolecular Complexes.

Comparison of structures of homological proteins often helps to understand functionally significant features of these structures. This concerns not only structures of separate protein chains, but also structures of macromolecular complexes. In particular, a comparison of complexes of homologous proteins with DNA is significant for analysis of the recognition of DNA by proteins. We present program LCore for detecting geometrical cores of a family of structures; a geometrical core is a set of amino acid residues and nucleotides that disposed similarly in all structures of the family. We describe the algorithm of the program, its web interface, and an example of its application to analysis of complexes of homeodomains with DNA.

StructAlign, a Program for Alignment of Structures of DNA-Protein Complexes.

Comparative analysis of structures of complexes of homologous proteins with DNA is important in the analysis of DNA-protein recognition. Alignment is a necessary stage of the analysis. An alignment is a matching of amino acid residues and nucleotides of one complex to residues and nucleotides of the other. Currently, there are no programs available for aligning structures of DNA-protein complexes. We present the program StructAlign, which should fill this gap. The program inputs a pair of complexes of DNA double helix with proteins and outputs an alignment of DNA chains corresponding to the best spatial fit of the protein chains.

Moss Phylogeny Reconstruction Using Nucleotide Pangenome of Complete Mitogenome Sequences.

Stability of composition and sequence of genes was shown earlier in 13 mitochondrial genomes of mosses (Rensing, S. A., et al. (2008) Science, 319, 64-69). It is of interest to study the evolution of mitochondrial genomes not only at the gene level, but also on the level of nucleotide sequences. To do this, we have constructed a "nucleotide pangenome" for mitochondrial genomes of 24 moss species. The nucleotide pangenome is a set of aligned nucleotide sequences of orthologous genome fragments covering the totality of all genomes. The nucleotide pangenome was constructed using specially developed new software, NPG-explorer (NPGe). The stable part of the mitochondrial genome (232 stable blocks) is shown to be, on average, 45% of its length. In the joint alignment of stable blocks, 82% of positions are conserved. The phylogenetic tree constructed with the NPGe program is in good correlation with other phylogenetic reconstructions. With the NPGe program, 30 blocks have been identified with repeats no shorter than 50 bp. The maximal length of a block with repeats is 140 bp. Duplications in the mitochondrial genomes of mosses are rare. On average, the genome contains about 500 bp in large duplications. The total length of insertions and deletions was determined in each genome. The losses and gains of DNA regions are rather active in mitochondrial genomes of mosses, and such rearrangements presumably can be used as additional markers in the reconstruction of phylogeny.

Role of Restriction-Modification Systems in Prokaryotic Evolution and Ecology.

Restriction-modification (R-M) systems are able to methylate or cleave DNA depending on methylation status of their recognition site. It allows them to protect bacterial cells from invasion by foreign DNA. Comparative analysis of a large number of available bacterial genomes and methylomes clearly demonstrates that the role of R-M systems in bacteria is wider than only defense. R-M systems maintain heterogeneity of a bacterial population and are involved in adaptation of bacteria to change in their environmental conditions. R-M systems can be essential for host colonization by pathogenic bacteria. Phase variation and intragenomic recombinations are sources of the fast evolution of the specificity of R-M systems. This review focuses on the influence of R-M systems on evolution and ecology of prokaryotes.

Hydrogen peroxide produced inside mitochondria takes part in cell-to-cell transmission of apoptotic signal.

In monolayer of HeLa cells treated with tumor necrosis factor (TNF), apoptotic cells formed clusters indicating possible transmission of apoptotic signal via the culture media. To investigate this phenomenon, a simple method of enabling two cell cultures to interact has been employed. Two coverslips were placed side by side in a Petri dish, one coverslip covered with apoptogen-treated cells (the inducer) and another with non-treated cells (the recipient). TNF, staurosporine, or H2O2 treatment of the inducer cells is shown to initiate apoptosis on the recipient coverslip. This effect is increased by a catalase inhibitor aminotriazole and is arrested by addition of catalase or by pre-treatment of either the inducer or the recipient cells with nanomolar concentrations of mitochondria-targeted cationic antioxidant MitoQ (10-(6 -ubiquinolyl)decyltriphenylphosphonium), which specifically arrests H2O2-induced apoptosis. The action of MitoQ is abolished by an uncoupler preventing accumulation of MitoQ in mitochondria. It is concluded that reactive oxygen species (ROS) produced by mitochondria in the apoptotic cells initiate the release of H2O2 from these cells. The H2O2 released is employed as a long-distance cell suicide messenger. In processing of such a signal by the recipient cells, mitochondrial ROS production is also involved. It is suggested that the described phenomenon may be involved in expansion of the apoptotic region around a damaged part of the tissue during heart attack or stroke as well as in "organoptosis", i.e. disappearance of organs during ontogenesis.