Conserved G-Quadruplex-Forming Sequences in Mammalian TERT Promoters and Their Effect on Mutation Frequency.

Somatic mutations in the promoter region of the human telomerase reverse transcriptase (hTERT) gene have been identified in many types of cancer. The hTERT promoter is known to be enriched with sequences that enable the formation of G-quadruplex (G4) structures, whose presence is associated with elevated mutagenicity and genome instability. Here, we used a bioinformatics tool (QGRS mapper) to search for G4-forming sequences (G4 motifs) in the 1000 bp TERT promoter regions of 141 mammalian species belonging to 20 orders, 5 of which, including primates and predators, contain more than 10 species. Groups of conserved G4 motifs and single-nucleotide variants within these groups were discovered using a block alignment approach (based on the Nucleotide PanGenome explorer). It has been shown that: (i) G4 motifs are predominantly located in the region proximal to the transcription start site (up to 400 bp) and are over-represented on the non-coding strand of the TERT promoters, (ii) 11 to 22% of the G4 motifs found are evolutionarily conserved across the related organisms, and (iii) a statistically significant higher frequency of nucleotide substitutions in the conserved G4 motifs compared to the surrounding regions was confirmed only for the order Primates. These data support the assumption that G4s can interfere with the DNA repair process and affect the evolutionary adaptation of organisms and species.

Corrigendum to "Comparative Analysis of Context-Dependent Mutagenesis Using Human and Mouse Models".

[This corrects the article DOI: 10.1155/2013/989410.].

Conserved features of complexes of TATA-box binding proteins with DNA.

A comparative analysis of all available structures of complexes of TATA-box binding proteins (TBPs) with DNA is performed. Conserved features of DNA-protein interaction are described, including nine amino acid residues that form conserved hydrogen bonds, 13 residues participating in formation of two conserved hydrophobic clusters at DNA-protein interface, and four conserved water-mediated contacts. Partial symmetry of conserved contacts reflects quasi-symmetry of TBP structure.

An updated version of NPIDB includes new classifications of DNA-protein complexes and their families.

The recent upgrade of nucleic acid-protein interaction database (NPIDB, http://npidb.belozersky.msu.ru/) includes a newly elaborated classification of complexes of protein domains with double-stranded DNA and a classification of families of related complexes. Our classifications are based on contacting structural elements of both DNA: the major groove, the minor groove and the backbone; and protein: helices, beta-strands and unstructured segments. We took into account both hydrogen bonds and hydrophobic interaction. The analyzed material contains 1942 structures of protein domains from 748 PDB entries. We have identified 97 interaction modes of individual protein domain-DNA complexes and 17 DNA-protein interaction classes of protein domain families. We analyzed the sources of diversity of DNA-protein interaction modes in different complexes of one protein domain family. The observed interaction mode is sometimes influenced by artifacts of crystallization or diversity in secondary structure assignment. The interaction classes of domain families are more stable and thus possess more biological sense than a classification of single complexes. Integration of the classification into NPIDB allows the user to browse the database according to the interacting structural elements of DNA and protein molecules. For each family, we present average DNA shape parameters in contact zones with domains of the family.

Two Cytoplasmic Acylation Sites and an Adjacent Hydrophobic Residue, but No Other Conserved Amino Acids in the Cytoplasmic Tail of HA from Influenza A Virus Are Crucial for Virus Replication.

Recruitment of the matrix protein M1 to the assembly site of the influenza virus is thought to be mediated by interactions with the cytoplasmic tail of hemagglutinin (HA). Based on a comprehensive sequence comparison of all sequences present in the database, we analyzed the effect of mutating conserved residues in the cytosol-facing part of the transmembrane region and cytoplasmic tail of HA (A/WSN/33 (H1N1) strain) on virus replication and morphology of virions. Removal of the two cytoplasmic acylation sites and substitution of a neighboring isoleucine by glutamine prevented rescue of infectious virions. In contrast, a conservative exchange of the same isoleucine, non-conservative exchanges of glycine and glutamine, deletion of the acylation site at the end of the transmembrane region and shifting it into the tail did not affect virus morphology and had only subtle effects on virus growth and on the incorporation of M1 and Ribo-Nucleoprotein Particles (RNPs). Thus, assuming that essential amino acids are conserved between HA subtypes we suggest that, besides the two cytoplasmic acylation sites (including adjacent hydrophobic residues), no other amino acids in the cytoplasmic tail of HA are indispensable for virus assembly and budding.

Lifespan of restriction-modification systems critically affects avoidance of their recognition sites in host genomes.

BACKGROUND: Avoidance of palindromic recognition sites of Type II restriction-modification (R-M) systems was shown for many R-M systems in dozens of prokaryotic genomes. However the phenomenon has not been investigated systematically for all presently available genomes and annotated R-M systems. We have studied all known recognition sites in thousands of prokaryotic genomes and found factors that influence their avoidance. RESULTS: Only Type II R-M systems consisting of independently acting endonuclease and methyltransferase (called 'orthodox' here) cause avoidance of their sites, both palindromic and asymmetric, in corresponding prokaryotic genomes; the avoidance takes place for ~ 50 % of 1774 studied cases. It is known that prokaryotes can acquire and lose R-M systems. Thus it is possible to talk about the lifespan of an R-M system in a genome. We have shown that the recognition site avoidance correlates with the lifespan of R-M systems. The sites of orthodox R-M systems that are encoded in host genomes for a long time are avoided more often (up to 100 % in certain cohorts) than the sites of recently acquired ones. We also found cases of site avoidance in absence of the corresponding R-M systems in the genome. An analysis of closely related bacteria shows that such avoidance can be a trace of lost R-M systems. Sites of Type I, IIС/G, IIM, III, and IV R-M systems are not avoided in vast majority of cases. CONCLUSIONS: The avoidance of orthodox Type II R-M system recognition sites in prokaryotic genomes is a widespread phenomenon. Presence of an R-M system without an underrepresentation of its site may indicate that the R-M system was acquired recently. At the same time, a significant underrepresentation of a site may be a sign of presence of the corresponding R-M system in this organism or in its ancestors for a long time. The drastic difference between site avoidance for orthodox Type II R-M systems and R-M systems of other types can be explained by a higher rate of specificity changes or a less self-toxicity of the latter.

Site-specific S-acylation of influenza virus hemagglutinin: the location of the acylation site relative to the membrane border is the decisive factor for attachment of stearate.

S-Acylation of hemagglutinin (HA), the main glycoprotein of influenza viruses, is an essential modification required for virus replication. Using mass spectrometry, we have previously demonstrated specific attachment of acyl chains to individual acylation sites. Whereas the two cysteines in the cytoplasmic tail of HA contain only palmitate, stearate is exclusively attached to a cysteine positioned at the end of the transmembrane region (TMR). Here we analyzed recombinant viruses containing HA with exchange of conserved amino acids adjacent to acylation sites or with a TMR cysteine shifted to a cytoplasmic location to identify the molecular signal that determines preferential attachment of stearate. We first developed a new protocol for sample preparation that requires less material and might thus also be suitable to analyze cellular proteins. We observed cell type-specific differences in the fatty acid pattern of HA: more stearate was attached if human viruses were grown in mammalian compared with avian cells. No underacylated peptides were detected in the mass spectra, and even mutations that prevented generation of infectious virus particles did not abolish acylation of expressed HA as demonstrated by metabolic labeling experiments with [(3)H]palmitate. Exchange of conserved amino acids in the vicinity of an acylation site had a moderate effect on the stearate content. In contrast, shifting the TMR cysteine to a cytoplasmic location virtually eliminated attachment of stearate. Thus, the location of an acylation site relative to the transmembrane span is the main signal for stearate attachment, but the sequence context and the cell type modulate the fatty acid pattern.

Co-evolution analysis to predict protein-protein interactions within influenza virus envelope.

Interactions between integral membrane proteins hemagglutinin (HA), neuraminidase (NA), M2 and membrane-associated matrix protein M1 of influenza A virus are thought to be crucial for assembly of functionally competent virions. We hypothesized that the amino acid residues located at the interface of two different proteins are under physical constraints and thus probably co-evolve. To predict co-evolving residue pairs, the EvFold ( http://evfold.org ) program searching the (nontransitive) Direct Information scores was applied for large samplings of amino acid sequences from Influenza Research Database ( http://www.fludb.org/ ). Having focused on the HA, NA, and M2 cytoplasmic tails as well as C-terminal domain of M1 (being the less conserved among the protein domains) we captured six pairs of correlated positions. Among them, there were one, two, and three position pairs for HA-M2, HA-M1, and M2-M1 protein pairs, respectively. As expected, no co-varying positions were found for NA-HA, NA-M1, and NA-M2 pairs obviously due to high conservation of the NA cytoplasmic tail. The sum of frequencies calculated for two major amino acid patterns observed in pairs of correlated positions was up to 0.99 meaning their high to extreme evolutionary sustainability. Based on the predictions a hypothetical model of pair-wise protein interactions within the viral envelope was proposed.

The broken MLL gene is frequently located outside the inherent chromosome territory in human lymphoid cells treated with DNA topoisomerase II poison etoposide.

The mixed lineage leukaemia (MLL) gene is frequently rearranged in secondary leukaemias, in which it could fuse to a variety of different partners. Breakage in the MLL gene preferentially occurs within a ~8 kb region that possesses a strong DNA topoisomerase II cleavage site. It has been proposed that DNA topoisomerase II-mediated DNA cleavage within this and other regions triggers translocations that occur due to incorrect joining of broken DNA ends. To further clarify a possible mechanism for MLL rearrangements, we analysed the frequency of MLL cleavage in cells exposed to etoposide, a DNA topoisomerase II poison commonly used as an anticancer drug, and positioning of the broken 3'-end of the MLL gene in respect to inherent chromosomal territories. It was demonstrated that exposure of human Jurkat cells to etoposide resulted in frequent cleavage of MLL genes. Using MLL-specific break-apart probes we visualised cleaved MLL genes in ~17% of nuclei. Using confocal microscopy and 3D modelling, we demonstrated that in cells treated with etoposide and cultivated for 1 h under normal conditions, ~9% of the broken MLL alleles were present outside the chromosome 11 territory, whereas in both control cells and cells inspected immediately after etoposide treatment, virtually all MLL alleles were present within the chromosomal territory. The data are discussed in the framework of the "breakage first" model of juxtaposing translocation partners. We propose that in the course of repairing DNA topoisomerase II-mediated DNA lesions (removal of stalled DNA topoisomerase II complexes and non-homologous end joining), DNA ends acquire additional mobility, which allows the meeting and incorrect joining of translocation partners.

NPIDB: Nucleic acid-Protein Interaction DataBase.

The Nucleic acid-Protein Interaction DataBase (http://npidb.belozersky.msu.ru/) contains information derived from structures of DNA-protein and RNA-protein complexes extracted from the Protein Data Bank (3846 complexes in October 2012). It provides a web interface and a set of tools for extracting biologically meaningful characteristics of nucleoprotein complexes. The content of the database is updated weekly. The current version of the Nucleic acid-Protein Interaction DataBase is an upgrade of the version published in 2007. The improvements include a new web interface, new tools for calculation of intermolecular interactions, a classification of SCOP families that contains DNA-binding protein domains and data on conserved water molecules on the DNA-protein interface.

Solitary restriction endonucleases in prokaryotic genomes.

Prokaryotic restriction-modification (R-M) systems defend the host cell from the invasion of a foreign DNA. They comprise two enzymatic activities: specific DNA cleavage activity and DNA methylation activity preventing cleavage. Typically, these activities are provided by two separate enzymes: a DNA methyltransferase (MTase) and a restriction endonuclease (RE). In the absence of a corresponding MTase, an RE of Type II R-M system is highly toxic for the cell. Genes of the R-M system are linked in the genome in the vast majority of annotated cases. There are only a few reported cases in which the genes of MTase and RE from one R-M system are not linked. Nevertheless, a few hundreds solitary RE genes are present in the Restriction Enzyme Database (http://rebase.neb.com) annotations. Using the comparative genomic approach, we analysed 272 solitary RE genes. For 57 solitary RE genes we predicted corresponding MTase genes located distantly in a genome. Of the 272 solitary RE genes, 99 are likely to be fragments of RE genes. Various explanations for the existence of the remaining 116 solitary RE genes are also discussed.

SheeP: A tool for description of ß-sheets in protein 3D structures.

The description of a protein fold is a hard problem due to significant variability of main structural units, ß-sheets and α-helixes, and their mutual arrangements. An adequate description of the structural units is an important step in objective protein structure classification, which to date is based on expert judgment in a number of cases. Explicit determination and description of structural units is more complicated for ß-sheets than for α-helixes due to ß-sheets variability both in composition and geometry. We have developed an algorithm that can significantly modify ß-sheets detected by commonly used DSSP and Stride algorithms and represent the result as a "ß-sheet map," a table describing certain ß-sheet features. In our approach, ß-sheets (rather than ß-strands) are considered as holistic objects. Both hydrogen bonds and geometrical restrains are explored for the determination of ß-sheets. The algorithm is implemented in SheeP program. It was tested for prediction architectures of domains from 93 well-defined all-ß and α/ß SCOP protein domain families, and showed 93% of correct results. The Web-service allows to detect ß-sheets in a given protein structure, visualize ß-sheet maps, as well as input three-dimensional structures with highlighted ß-sheets and their structural features.

New words in human mutagenesis.

BACKGROUND: The substitution rates within different nucleotide contexts are subject to varying levels of bias. The most well known example of such bias is the excess of C to T (C > T) mutations in CpG (CG) dinucleotides. The molecular mechanisms underlying this bias are important factors in human genome evolution and cancer development. The discovery of other nucleotide contexts that have profound effects on substitution rates can improve our understanding of how mutations are acquired, and why mutation hotspots exist. RESULTS: We compared rates of inherited mutations in 1-4 bp nucleotide contexts using reconstructed ancestral states of human single nucleotide polymorphisms (SNPs) from intergenic regions. Chimp and orangutan genomic sequences were used as outgroups. We uncovered 3.5 and 3.3-fold excesses of T > C mutations in the second position of ATTG and ATAG words, respectively, and a 3.4-fold excess of A > C mutations in the first position of the ACAA word. CONCLUSIONS: Although all the observed biases are less pronounced than the 5.1-fold excess of C > T mutations in CG dinucleotides, the three 4 bp mutation contexts mentioned above (and their complementary contexts) are well distinguished from all other mutation contexts. This provides a challenge to discover the underlying mechanisms responsible for the observed excesses of mutations.

Linker and/or transmembrane regions of influenza A/Group-1, A/Group-2, and type B virus hemagglutinins are packed differently within trimers.

Influenza virus hemagglutinin is a homotrimeric spike glycoprotein crucial for virions' attachment, membrane fusion, and assembly reactions. X-ray crystallography data are available for hemagglutinin ectodomains of various types/subtypes but not for anchoring segments. To get structural information for the linker and transmembrane regions of hemagglutinin, influenza A (H1-H16 subtypes except H8 and H15) and B viruses were digested with bromelain or subtilisin Carlsberg, either within virions or in non-ionic detergent micelles. Proteolytical fragments were analyzed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Within virions, hemagglutinins of most influenza A/Group-1 and type B virus strains were more susceptible to digestion with bromelain and/or subtilisin compared to A/Group-2 hemagglutinins. The cleavage sites were always located in the hemagglutinin linker sequence. In detergent, 1) bromelain cleaved hemagglutinin of every influenza A subtype in the linker region; 2) subtilisin cleaved Group-2 hemagglutinins in the linker region; 3) subtilisin cleaved Group-1 hemagglutinins in the transmembrane region; 4) both enzymes cleaved influenza B virus hemagglutinin in the transmembrane region. We propose that the A/Group-2 hemagglutinin linker and/or transmembrane regions are more tightly associated within trimers than type A/Group-1 and particularly type B ones. This hypothesis is underpinned by spatial trimeric structure modeling performed for transmembrane regions of both Group-1 and Group-2 hemagglutinin representatives. Differential S-acylation of the hemagglutinin C-terminal anchoring segment with palmitate/stearate residues possibly contributes to fine tuning of transmembrane trimer packing and stabilization since decreased stearate amount correlated with deeper digestion of influenza B and some A/Group-1 hemagglutinins.

Identification of conserved features of LAGLIDADG homing endonucleases.

LAGLIDADG family of homing endonucleases are rare-cutting enzymes which recognize long target sequences and are of great interest in genome engineering. Despite advances in homing endonuclease engineering, effective methods of broadening the range of cleaved sequences are still lacking. Here, we present a study of conserved structural features of LAGLIDADG homing endonucleases that might aid further development of such methods. The protein-DNA interface of LAGLIDADG homing endonucleases differs considerably with the particular nuclease, and the analysis of conserved protein-DNA interactions could not identify any residues crucial for DNA binding and common to most nucleases of the family. For the homing endonuclease PI-SceI, a comparison of structural and experimental data derived from literature helped to identify 23 residues that are likely to be important for DNA binding. Analysis of the LAGLIDADG domain dimerization interface allowed the choosing of six positions that contribute to dimerization specificity most, while comparison of 446 sequences of LAGLIDADG endonucleases revealed groups of residues in these positions that appear to be most favorable for dimerization.

Malakite: an automatic tool for characterisation of structure of reliable blocks in multiple alignments of protein sequences.

It makes sense to speak of alignment of protein sequences only within the regions, where the sequences are related to each other. This simple consideration is often disregarded by programs of multiple alignment construction. A package for alignment analysis MAlAKiTE (Multiple Alignment Automatic Kinship Tiling Engine) is introduced. It aims to find the blocks of reliable alignment, which contain related regions only, within the whole alignment and allows for dealing with them. The validity of the detection of reliable blocks' was verified by comparison with structural data.

Exclusive sequences of different genomes.

We studied the distribution of 1-7 bp words in a dataset that includes 139 complete eukaryotic genomes, 33 masked eukaryotic genomes and coding regions from 35 genomes. We tested different statistical models to determine over- and under-represented words. The method described by Karlin et al. has the strongest predictive power compared to other methods. Using this method we identified over- and under-represented words consistent within a large array of taxonomic groups. Some of those words have not yet been described as exclusive. For example, CGCG is over-represented in CG-deficient organisms. We also describe exceptions for widely known exclusive words, such as CG and TA.

Conserved water molecules in X-ray structures highlight the role of water in intramolecular and intermolecular interactions.

Water molecules immobilized on a protein or DNA surface are known to play an important role in intramolecular and intermolecular interactions. Comparative analysis of related three-dimensional (3D) structures allows to predict the locations of such water molecules on the protein surface. We have developed and implemented the algorithm WLAKE detecting "conserved" water molecules, i.e. those located in almost the same positions in a set of superimposed structures of related proteins or macromolecular complexes. The problem is reduced to finding maximal cliques in a certain graph. Despite exponential algorithm complexity, the program works appropriately fast for dozens of superimposed structures. WLAKE was used to predict functionally significant water molecules in enzyme active sites (transketolases) as well as in intermolecular (ETS-DNA complexes) and intramolecular (thiol-disulfide interchange protein) interactions. The program is available online at http://monkey.belozersky.msu.ru/~evgeniy/wLake/wLake.html.

NPIDB: a database of nucleic acids-protein interactions.

UNLABELLED: The database NPIDB (Nucleic Acids-Protein Interaction DataBase) contains information derived from structures of DNA-protein and RNA-protein complexes extracted from PDB (1834 complexes in July 2007). It is organized as a collection of files in PDB format and is equipped with a web-interface and a set of tools for extracting biologically meaningful characteristics of complexes. The content of the database is weekly updated. AVAILABILITY: http://monkey.belozersky.msu.ru/NPIDB/

Changes in chromosome positioning may contribute to the development of diseases related to X-chromosome aneuploidy.

The radial positions of the centromeric regions of chromosomes 1 and X were determined in normal male fibroblasts (XY) and in fibroblasts from a patient with a rare case of XXXXY polysomy. The centromeric regions and presumably the whole territories of active X chromosomes were demonstrated to occupy similar, although not identical, positions in XY and XXXXY cells. The centromeres of inactive X chromosomes (Barr bodies) were located closer to the nuclear periphery as compared with the centromeres of active X chromosomes. In addition, it was established that the nuclear radial position of gene-rich chromosome 1 was changed in XXXXY cells as compared to normal XY cells. The data are discussed in the context of the hypothesis postulating that changes in nuclear positioning of chromosomal territories induced by the presence of extra copies of individual chromosomes may contribute to the development of diseases related to different polysomies.