PMCC PMCC

Search tips
Search criteria

Advanced
Results 1-25 (122)
 

Clipboard (0)
None

Select a Filter Below

Year of Publication
more »
1.  RNA methyltransferases involved in 5′ cap biosynthesis 
RNA Biology  2015;11(12):1597-1607.
In eukaryotes and viruses that infect them, the 5′ end of mRNA molecules, and also many other functionally important RNAs, are modified to form a so-called cap structure that is important for interactions of these RNAs with many nuclear and cytoplasmic proteins. The RNA cap has multiple roles in gene expression, including enhancement of RNA stability, splicing, nucleocytoplasmic transport, and translation initiation. Apart from guanosine addition to the 5′ end in the most typical cap structure common to transcripts produced by RNA polymerase II (in particular mRNA), essentially all cap modifications are due to methylation. The complexity of the cap structure and its formation can range from just a single methylation of the unprocessed 5′ end of the primary transcript, as in mammalian U6 and 7SK, mouse B2, and plant U3 RNAs, to an elaborate m7Gpppm6,6AmpAmpCmpm3Um structure at the 5′ end of processed RNA in trypanosomes, which are formed by as many as 8 methylation reactions. While all enzymes responsible for methylation of the cap structure characterized to date were found to belong to the same evolutionarily related and structurally similar Rossmann Fold Methyltransferase superfamily, that uses the same methyl group donor, S-adenosylmethionine; the enzymes also exhibit interesting differences that are responsible for their distinct functions. This review focuses on the evolutionary classification of enzymes responsible for cap methylation in RNA, with a focus on the sequence relationships and structural similarities and dissimilarities that provide the basis for understanding the mechanism of biosynthesis of different caps in cellular and viral RNAs. Particular attention is paid to the similarities and differences between methyltransferases from human cells and from human pathogens that may be helpful in the development of antiviral and antiparasitic drugs.
doi:10.1080/15476286.2015.1004955
PMCID: PMC4615557  PMID: 25626080
antiviral drugs; cap; crystallography; methylation; modified nucleotides; mRNA; post-transcriptional modification; RNA maturation; RNA modification; trypanosomes
2.  A composite double-/single-stranded RNA-binding region in protein Prp3 supports tri-snRNP stability and splicing 
eLife  null;4:e07320.
Prp3 is an essential U4/U6 di-snRNP-associated protein whose functions and molecular mechanisms in pre-mRNA splicing are presently poorly understood. We show by structural and biochemical analyses that Prp3 contains a bipartite U4/U6 di-snRNA-binding region comprising an expanded ferredoxin-like fold, which recognizes a 3′-overhang of U6 snRNA, and a preceding peptide, which binds U4/U6 stem II. Phylogenetic analyses revealed that the single-stranded RNA-binding domain is exclusively found in Prp3 orthologs, thus qualifying as a spliceosome-specific RNA interaction module. The composite double-stranded/single-stranded RNA-binding region assembles cooperatively with Snu13 and Prp31 on U4/U6 di-snRNAs and inhibits Brr2-mediated U4/U6 di-snRNA unwinding in vitro. RNP-disrupting mutations in Prp3 lead to U4/U6•U5 tri-snRNP assembly and splicing defects in vivo. Our results reveal how Prp3 acts as an important bridge between U4/U6 and U5 in the tri-snRNP and comparison with a Prp24-U6 snRNA recycling complex suggests how Prp3 may be involved in U4/U6 reassembly after splicing.
DOI: http://dx.doi.org/10.7554/eLife.07320.001
eLife digest
Proteins are built following instructions contained within the DNA of gene sequences. This genetic information is copied into short-lived molecules, called messenger RNAs (or mRNAs), which move away from the DNA and are then decoded by the molecular machines that build proteins. However, mRNA sequences often have to be edited before they are used. Another molecular machine, called a spliceosome, carries out some of this editing.
A spliceosome is formed from a number of smaller subunits, including three RNA-protein particles that each contain one RNA molecule (called U1, U2 and U5), and one particle that contains two RNA molecules (called U4 and U6). These subunits must assemble around an unedited mRNA in a particular order so that the spliceosome can work correctly. Once the mRNA has been edited, and the spliceosome has performed its job, these complexes need to disassemble so that they are ready to be reassembled around a new mRNA molecule. A protein called Prp3 is known to be involved in these assembly, disassembly and reassembly steps. However, it is unclear how this protein performs these activities.
Liu et al. have now used structural biology and biochemical techniques to determine the three-dimensional structure of Prp3, and have shown that this protein has a “two-part” binding site that binds to the RNA molecules in the U4/U6 subunit of the spliceosome. Further analyses revealed that one of these features is only found in Prp3 and not in other types of RNA-binding proteins.
Together with previous work, Liu et al. also reveal that Prp3 can serve as a ‘bridge’ between the U4/U6 and U5 subunits of the spliceosome, and suggest how these features allow the two subunits to group together before they are incorporated into a spliceosome.
Notably, certain mutations in the gene for the Prp3 protein lead to a human eye disease called retinitis pigmentosa. In the future it will be important to investigate if the above activities are affected in the mutant variants of the Prp3 protein.
DOI: http://dx.doi.org/10.7554/eLife.07320.002
doi:10.7554/eLife.07320
PMCID: PMC4520091  PMID: 26161500
small nuclear ribonucleoprotein particle; pre-mRNA splicing; protein–RNA interaction; X-ray crystallography; human; S. cerevisiae
3.  Distribution and frequencies of post-transcriptional modifications in tRNAs 
RNA Biology  2015;11(12):1619-1629.
Functional tRNA molecules always contain a wide variety of post-transcriptionally modified nucleosides. These modifications stabilize tRNA structure, allow for proper interaction with other macromolecules and fine-tune the decoding of mRNAs during translation. Their presence in functionally important regions of tRNA is conserved in all domains of life. However, the identities of many of these modified residues depend much on the phylogeny of organisms the tRNAs are found in, attesting for domain-specific strategies of tRNA maturation. In this work we present a new tool, tRNAmodviz web server (http://genesilico.pl/trnamodviz) for easy comparative analysis and visualization of modification patterns in individual tRNAs, as well as in groups of selected tRNA sequences. We also present results of comparative analysis of tRNA sequences derived from 7 phylogenetically distinct groups of organisms: Gram-negative bacteria, Gram-positive bacteria, cytosol of eukaryotic single cell organisms, Fungi and Metazoa, cytosol of Viridiplantae, mitochondria, plastids and Euryarchaeota. These data update the study conducted 20 y ago with the tRNA sequences available at that time.
doi:10.4161/15476286.2014.992273
PMCID: PMC4615829  PMID: 25611331
comparative analysis; evolution; modified nucleotides; post-transcriptional modification; RNA maturation; tRNA; tRNA modifications; tRNA sequence; web server
4.  Computational modeling of RNA 3D structures, with the aid of experimental restraints 
RNA Biology  2014;11(5):522-536.
In addition to mRNAs whose primary function is transmission of genetic information from DNA to proteins, numerous other classes of RNA molecules exist, which are involved in a variety of functions, such as catalyzing biochemical reactions or performing regulatory roles. In analogy to proteins, the function of RNAs depends on their structure and dynamics, which are largely determined by the ribonucleotide sequence. Experimental determination of high-resolution RNA structures is both laborious and difficult, and therefore, the majority of known RNAs remain structurally uncharacterized. To address this problem, computational structure prediction methods were developed that simulate either the physical process of RNA structure formation (“Greek science” approach) or utilize information derived from known structures of other RNA molecules (“Babylonian science” approach). All computational methods suffer from various limitations that make them generally unreliable for structure prediction of long RNA sequences. However, in many cases, the limitations of computational and experimental methods can be overcome by combining these two complementary approaches with each other. In this work, we review computational approaches for RNA structure prediction, with emphasis on implementations (particular programs) that can utilize restraints derived from experimental analyses. We also list experimental approaches, whose results can be relatively easily used by computational methods. Finally, we describe case studies where computational and experimental analyses were successfully combined to determine RNA structures that would remain out of reach for each of these approaches applied separately.
doi:10.4161/rna.28826
PMCID: PMC4152360  PMID: 24785264
RNA structure; RNA structure prediction; macromolecular modeling; bioinformatics; chemical probing
5.  RNA Bricks—a database of RNA 3D motifs and their interactions 
Nucleic Acids Research  2013;42(Database issue):D123-D131.
The RNA Bricks database (http://iimcb.genesilico.pl/rnabricks), stores information about recurrent RNA 3D motifs and their interactions, found in experimentally determined RNA structures and in RNA–protein complexes. In contrast to other similar tools (RNA 3D Motif Atlas, RNA Frabase, Rloom) RNA motifs, i.e. ‘RNA bricks’ are presented in the molecular environment, in which they were determined, including RNA, protein, metal ions, water molecules and ligands. All nucleotide residues in RNA bricks are annotated with structural quality scores that describe real-space correlation coefficients with the electron density data (if available), backbone geometry and possible steric conflicts, which can be used to identify poorly modeled residues. The database is also equipped with an algorithm for 3D motif search and comparison. The algorithm compares spatial positions of backbone atoms of the user-provided query structure and of stored RNA motifs, without relying on sequence or secondary structure information. This enables the identification of local structural similarities among evolutionarily related and unrelated RNA molecules. Besides, the search utility enables searching ‘RNA bricks’ according to sequence similarity, and makes it possible to identify motifs with modified ribonucleotide residues at specific positions.
doi:10.1093/nar/gkt1084
PMCID: PMC3965019  PMID: 24220091
6.  A toolbox for developing bioinformatics software 
Briefings in Bioinformatics  2011;13(2):244-257.
Creating useful software is a major activity of many scientists, including bioinformaticians. Nevertheless, software development in an academic setting is often unsystematic, which can lead to problems associated with maintenance and long-term availibility. Unfortunately, well-documented software development methodology is difficult to adopt, and technical measures that directly improve bioinformatic programming have not been described comprehensively. We have examined 22 software projects and have identified a set of practices for software development in an academic environment. We found them useful to plan a project, support the involvement of experts (e.g. experimentalists), and to promote higher quality and maintainability of the resulting programs. This article describes 12 techniques that facilitate a quick start into software engineering. We describe 3 of the 22 projects in detail and give many examples to illustrate the usage of particular techniques. We expect this toolbox to be useful for many bioinformatics programming projects and to the training of scientific programmers.
doi:10.1093/bib/bbr035
PMCID: PMC3294241  PMID: 21803787
software development; programming; project management; software quality
7.  Rational engineering of sequence specificity in R.MwoI restriction endonuclease 
Nucleic Acids Research  2012;40(17):8579-8592.
R.MwoI is a Type II restriction endonucleases enzyme (REase), which specifically recognizes a palindromic interrupted DNA sequence 5′-GCNNNNNNNGC-3′ (where N indicates any nucleotide), and hydrolyzes the phosphodiester bond in the DNA between the 7th and 8th base in both strands. R.MwoI exhibits remote sequence similarity to R.BglI, a REase with known structure, which recognizes an interrupted palindromic target 5′-GCCNNNNNGGC-3′. A homology model of R.MwoI in complex with DNA was constructed and used to predict functionally important amino acid residues that were subsequently targeted by mutagenesis. The model, together with the supporting experimental data, revealed regions important for recognition of the common bases in DNA sequences recognized by R.BglI and R.MwoI. Based on the bioinformatics analysis, we designed substitutions of the S310 residue in R.MwoI to arginine or glutamic acid, which led to enzyme variants with altered sequence selectivity compared with the wild-type enzyme. The S310R variant of R.MwoI preferred the 5′-GCCNNNNNGGC-3′ sequence as a target, similarly to R.BglI, whereas the S310E variant preferentially cleaved a subset of the MwoI sites, depending on the identity of the 3rd and 9th nucleotide residues. Our results represent a case study of a REase sequence specificity alteration by a single amino acid substitution, based on a theoretical model in the absence of a crystal structure.
doi:10.1093/nar/gks570
PMCID: PMC3458533  PMID: 22735699
8.  Databases and Bioinformatics Tools for the Study of DNA Repair 
DNA is continuously exposed to many different damaging agents such as environmental chemicals, UV light, ionizing radiation, and reactive cellular metabolites. DNA lesions can result in different phenotypical consequences ranging from a number of diseases, including cancer, to cellular malfunction, cell death, or aging. To counteract the deleterious effects of DNA damage, cells have developed various repair systems, including biochemical pathways responsible for the removal of single-strand lesions such as base excision repair (BER) and nucleotide excision repair (NER) or specialized polymerases temporarily taking over lesion-arrested DNA polymerases during the S phase in translesion synthesis (TLS). There are also other mechanisms of DNA repair such as homologous recombination repair (HRR), nonhomologous end-joining repair (NHEJ), or DNA damage response system (DDR). This paper reviews bioinformatics resources specialized in disseminating information about DNA repair pathways, proteins involved in repair mechanisms, damaging agents, and DNA lesions.
doi:10.4061/2011/475718
PMCID: PMC3200286  PMID: 22091405
9.  Structural basis for the methylation of A1408 in 16S rRNA by a panaminoglycoside resistance methyltransferase NpmA from a clinical isolate and analysis of the NpmA interactions with the 30S ribosomal subunit 
Nucleic Acids Research  2010;39(5):1903-1918.
NpmA, a methyltransferase that confers resistance to aminoglycosides was identified in an Escherichia coli clinical isolate. It belongs to the kanamycin–apramycin methyltransferase (Kam) family and specifically methylates the 16S rRNA at the N1 position of A1408. We determined the structures of apo-NpmA and its complexes with S-adenosylmethionine (AdoMet) and S-adenosylhomocysteine (AdoHcy) at 2.4, 2.7 and 1.68 Å, respectively. We generated a number of NpmA variants with alanine substitutions and studied their ability to bind the cofactor, to methylate A1408 in the 30S subunit, and to confer resistance to kanamycin in vivo. Residues D30, W107 and W197 were found to be essential. We have also analyzed the interactions between NpmA and the 30S subunit by footprinting experiments and computational docking. Helices 24, 42 and 44 were found to be the main NpmA-binding site. Both experimental and theoretical analyses suggest that NpmA flips out the target nucleotide A1408 to carry out the methylation. NpmA is plasmid-encoded and can be transferred between pathogenic bacteria; therefore it poses a threat to the successful use of aminoglycosides in clinical practice. The results presented here will assist in the development of specific NpmA inhibitors that could restore the potential of aminoglycoside antibiotics.
doi:10.1093/nar/gkq1033
PMCID: PMC3061052  PMID: 21062819
10.  Structural basis for the methylation of G1405 in 16S rRNA by aminoglycoside resistance methyltransferase Sgm from an antibiotic producer: a diversity of active sites in m7G methyltransferases 
Nucleic Acids Research  2010;38(12):4120-4132.
Sgm (Sisomicin-gentamicin methyltransferase) from antibiotic-producing bacterium Micromonospora zionensis is an enzyme that confers resistance to aminoglycosides like gentamicin and sisomicin by specifically methylating G1405 in bacterial 16S rRNA. Sgm belongs to the aminoglycoside resistance methyltransferase (Arm) family of enzymes that have been recently found to spread by horizontal gene transfer among disease-causing bacteria. Structural characterization of Arm enzymes is the key to understand their mechanism of action and to develop inhibitors that would block their activity. Here we report the structure of Sgm in complex with cofactors S-adenosylmethionine (AdoMet) and S-adenosylhomocysteine (AdoHcy) at 2.0 and 2.1 Å resolution, respectively, and results of mutagenesis and rRNA footprinting, and protein-substrate docking. We propose the mechanism of methylation of G1405 by Sgm and compare it with other m7G methyltransferases, revealing a surprising diversity of active sites and binding modes for the same basic reaction of RNA modification. This analysis can serve as a stepping stone towards developing drugs that would specifically block the activity of Arm methyltransferases and thereby re-sensitize pathogenic bacteria to aminoglycoside antibiotics.
doi:10.1093/nar/gkq122
PMCID: PMC2896518  PMID: 20194115
11.  MODOMICS: a database of RNA modification pathways. 2008 update 
Nucleic Acids Research  2008;37(Database issue):D118-D121.
MODOMICS, a database devoted to the systems biology of RNA modification, has been subjected to substantial improvements. It provides comprehensive information on the chemical structure of modified nucleosides, pathways of their biosynthesis, sequences of RNAs containing these modifications and RNA-modifying enzymes. MODOMICS also provides cross-references to other databases and to literature. In addition to the previously available manually curated tRNA sequences from a few model organisms, we have now included additional tRNAs and rRNAs, and all RNAs with 3D structures in the Nucleic Acid Database, in which modified nucleosides are present. In total, 3460 modified bases in RNA sequences of different organisms have been annotated. New RNA-modifying enzymes have been also added. The current collection of enzymes includes mainly proteins for the model organisms Escherichia coli and Saccharomyces cerevisiae, and is currently being expanded to include proteins from other organisms, in particular Archaea and Homo sapiens. For enzymes with known structures, links are provided to the corresponding Protein Data Bank entries, while for many others homology models have been created. Many new options for database searching and querying have been included. MODOMICS can be accessed at http://genesilico.pl/modomics.
doi:10.1093/nar/gkn710
PMCID: PMC2686465  PMID: 18854352
12.  The architecture of the Schizosaccharomyces pombe CCR4-NOT complex 
Nature Communications  2016;7:10433.
CCR4-NOT is a large protein complex present both in cytoplasm and the nucleus of eukaryotic cells. Although it is involved in a variety of distinct processes related to expression of genetic information such as poly(A) tail shortening, transcription regulation, nuclear export and protein degradation, there is only fragmentary information available on some of its nine subunits. Here we show a comprehensive structural characterization of the native CCR4-NOT complex from Schizosaccharomyces pombe. Our cryo-EM 3D reconstruction of the complex, combined with techniques such as immunomicroscopy, RNA-nanogold labelling, docking of the available high-resolution structures and models of different subunits and domains, allow us to propose its full molecular architecture. We locate all functionally defined domains endowed with deadenylating and ubiquitinating activities, the nucleus-specific RNA-interacting subunit Mmi1, as well as surfaces responsible for protein–protein interactions. This information provides insight into cooperation of the different CCR4-NOT complex functions.
CCR4-NOT is a protein complex involved in a variety of important genetic processes. Here, the authors report the mid-resolution structure of this complex, and model the positions and contacts between the subunits, providing structural support for the previously reported functions of the complex.
doi:10.1038/ncomms10433
PMCID: PMC4737751  PMID: 26804377
13.  Loss of Conserved Noncoding RNAs in Genomes of Bacterial Endosymbionts 
Genome Biology and Evolution  2016;8(2):426-438.
The genomes of intracellular symbiotic or pathogenic bacteria, such as of Buchnera, Mycoplasma, and Rickettsia, are typically smaller compared with their free-living counterparts. Here we showed that noncoding RNA (ncRNA) families, which are conserved in free-living bacteria, frequently could not be detected by computational methods in the small genomes. Statistical tests demonstrated that their absence is not an artifact of low GC content or small deletions in these small genomes, and thus it was indicative of an independent loss of ncRNAs in different endosymbiotic lineages. By analyzing the synteny (conservation of gene order) between the reduced and nonreduced genomes, we revealed instances of protein-coding genes that were preserved in the reduced genomes but lost cis-regulatory elements. We found that the loss of cis-regulatory ncRNA sequences, which regulate the expression of cognate protein-coding genes, is characterized by the reduction of secondary structure formation propensity, GC content, and length of the corresponding genomic regions.
doi:10.1093/gbe/evw007
PMCID: PMC4779614  PMID: 26782934
endosymbionts; noncoding RNA loss; covariance models; Rfam
14.  Crystal structure of Bacillus subtilis TrmB, the tRNA (m7G46) methyltransferase 
Nucleic Acids Research  2006;34(6):1925-1934.
The structure of Bacillus subtilis TrmB (BsTrmB), the tRNA (m7G46) methyltransferase, was determined at a resolution of 2.1 Å. This is the first structure of a member of the TrmB family to be determined by X-ray crystallography. It reveals a unique variant of the Rossmann-fold methyltransferase (RFM) structure, with the N-terminal helix folded on the opposite site of the catalytic domain. The architecture of the active site and a computational docking model of BsTrmB in complex with the methyl group donor S-adenosyl-l-methionine and the tRNA substrate provide an explanation for results from mutagenesis studies of an orthologous enzyme from Escherichia coli (EcTrmB). However, unlike EcTrmB, BsTrmB is shown here to be dimeric both in the crystal and in solution. The dimer interface has a hydrophobic core and buries a potassium ion and five water molecules. The evolutionary analysis of the putative interface residues in the TrmB family suggests that homodimerization may be a specific feature of TrmBs from Bacilli, which may represent an early stage of evolution to an obligatory dimer.
doi:10.1093/nar/gkl116
PMCID: PMC1447647  PMID: 16600901
15.  Structural and functional insights into tRNA binding and adenosine N1-methylation by an archaeal Trm10 homologue 
Nucleic Acids Research  2015;44(2):940-953.
Purine nucleosides on position 9 of eukaryal and archaeal tRNAs are frequently modified in vivo by the post-transcriptional addition of a methyl group on their N1 atom. The methyltransferase Trm10 is responsible for this modification in both these domains of life. While certain Trm10 orthologues specifically methylate either guanosine or adenosine at position 9 of tRNA, others have a dual specificity. Until now structural information about this enzyme family was only available for the catalytic SPOUT domain of Trm10 proteins that show specificity toward guanosine. Here, we present the first crystal structure of a full length Trm10 orthologue specific for adenosine, revealing next to the catalytic SPOUT domain also N- and C-terminal domains. This structure hence provides crucial insights in the tRNA binding mechanism of this unique monomeric family of SPOUT methyltransferases. Moreover, structural comparison of this adenosine-specific Trm10 orthologue with guanosine-specific Trm10 orthologues suggests that the N1 methylation of adenosine relies on additional catalytic residues.
doi:10.1093/nar/gkv1369
PMCID: PMC4737155  PMID: 26673726
16.  Reassignment of specificities of two cap methyltransferase domains in the reovirus lambda2 protein 
Genome Biology  2001;2(9):research0038.1-research0038.6.
Background
The reovirus λ2 protein catalyzes mRNA capping, that is, addition of a guanosine to the 5' end of each transcript in a 5'-to-5' orientation, as well as transfer of a methyl group from S-adenosyl-L-methionine (AdoMet) to the N7 atom of the added guanosyl moiety and subsequently to the ribose 2'-O atom of the first template-encoded nucleotide. The structure of the human reovirus core has been solved at 3.6 Å resolution, revealing a series of domains that include a putative guanylyltransferase domain and two putative methyltransferase (MTase) domains. It has been suggested that the order of domains in the λ2 protein corresponds to the order of reactions in the pathway and that the m7G (cap 0) and the 2'-O-ribose (cap 1) MTase activities may be exerted by the MTase 1 and the MTase 2 domains, respectively.
Results
We show that the reovirus MTase 1 domain shares a putative active site with the structurally characterized 2'-O-ribose MTases, including vaccinia virus cap 1 MTase, whereas the MTase 2 domain is structurally similar to glycine N-MTase.
Conclusions
On the basis of our analysis of the structural details we propose that the previously suggested functional assignments of the MTase 1 and MTase 2 domains should be swapped.
PMCID: PMC56899  PMID: 11574057
17.  Phylogenomics and sequence-structure-function relationships in the GmrSD family of Type IV restriction enzymes 
BMC Bioinformatics  2015;16:336.
Background
GmrSD is a modification-dependent restriction endonuclease that specifically targets and cleaves glucosylated hydroxymethylcytosine (glc-HMC) modified DNA. It is encoded either as two separate single-domain GmrS and GmrD proteins or as a single protein carrying both domains. Previous studies suggested that GmrS acts as endonuclease and NTPase whereas GmrD binds DNA.
Methods
In this work we applied homology detection, sequence conservation analysis, fold recognition and homology modeling methods to study sequence-structure-function relationships in the GmrSD restriction endonucleases family. We also analyzed the phylogeny and genomic context of the family members.
Results
Results of our comparative genomics study show that GmrS exhibits similarity to proteins from the ParB/Srx fold which can have both NTPase and nuclease activity. In contrast to the previous studies though, we attribute the nuclease activity also to GmrD as we found it to contain the HNH endonuclease motif. We revealed residues potentially important for structure and function in both domains. Moreover, we found that GmrSD systems exist predominantly as a fused, double-domain form rather than as a heterodimer and that their homologs are often encoded in regions enriched in defense and gene mobility-related elements. Finally, phylogenetic reconstructions of GmrS and GmrD domains revealed that they coevolved and only few GmrSD systems appear to be assembled from distantly related GmrS and GmrD components.
Conclusions
Our study provides insight into sequence-structure-function relationships in the yet poorly characterized family of Type IV restriction enzymes. Comparative genomics allowed to propose possible role of GmrD domain in the function of the GmrSD enzyme and possible active sites of both GmrS and GmrD domains. Presented results can guide further experimental characterization of these enzymes.
Electronic supplementary material
The online version of this article (doi:10.1186/s12859-015-0773-z) contains supplementary material, which is available to authorized users.
doi:10.1186/s12859-015-0773-z
PMCID: PMC4619093  PMID: 26493560
Restriction-Modification systems; Modification-Dependent systems; Type IV; Comparative genomics; Defense islands; Fold recognition; HNH endonuclease; ParB/Srx fold
19.  Biochemical Characterization and Validation of a Catalytic Site of a Highly Thermostable Ts2631 Endolysin from the Thermus scotoductus Phage vB_Tsc2631 
PLoS ONE  2015;10(9):e0137374.
Phage vB_Tsc2631 infects the extremophilic bacterium Thermus scotoductus MAT2631 and uses the Ts2631 endolysin for the release of its progeny. The Ts2631 endolysin is the first endolysin from thermophilic bacteriophage with an experimentally validated catalytic site. In silico analysis and computational modelling of the Ts2631 endolysin structure revealed a conserved Zn2+ binding site (His30, Tyr58, His131 and Cys139) similar to Zn2+ binding site of eukaryotic peptidoglycan recognition proteins (PGRPs). We have shown that the Ts2631 endolysin lytic activity is dependent on divalent metal ions (Zn2+ and Ca2+). The Ts2631 endolysin substitution variants H30N, Y58F, H131N and C139S dramatically lost their antimicrobial activity, providing evidence for the role of the aforementioned residues in the lytic activity of the enzyme. The enzyme has proven to be not only thermoresistant, retaining 64.8% of its initial activity after 2 h at 95°C, but also highly thermodynamically stable (Tm = 99.82°C, ΔHcal = 4.58 × 104 cal mol-1). Substitutions of histidine residues (H30N and H131N) and a cysteine residue (C139S) resulted in variants aggregating at temperatures ≥75°C, indicating a significant role of these residues in enzyme thermostability. The substrate spectrum of the Ts2631 endolysin included extremophiles of the genus Thermus but also Gram-negative mesophiles, such as Escherichia coli, Salmonella panama, Pseudomonas fluorescens and Serratia marcescens. The broad substrate spectrum and high thermostability of this endolysin makes it a good candidate for use as an antimicrobial agent to combat Gram-negative pathogens.
doi:10.1371/journal.pone.0137374
PMCID: PMC4573324  PMID: 26375388
21.  An atlas of RNA base pairs involving modified nucleobases with optimal geometries and accurate energies 
Nucleic Acids Research  2015;43(14):6714-6729.
Posttranscriptional modifications greatly enhance the chemical information of RNA molecules, contributing to explain the diversity of their structures and functions. A significant fraction of RNA experimental structures available to date present modified nucleobases, with half of them being involved in H-bonding interactions with other bases, i.e. ‘modified base pairs’. Herein we present a systematic investigation of modified base pairs, in the context of experimental RNA structures. To this end, we first compiled an atlas of experimentally observed modified base pairs, for which we recorded occurrences and structural context. Then, for each base pair, we selected a representative for subsequent quantum mechanics calculations, to find out its optimal geometry and interaction energy. Our structural analyses show that most of the modified base pairs are non Watson–Crick like and are involved in RNA tertiary structure motifs. In addition, quantum mechanics calculations quantify and provide a rationale for the impact of the different modifications on the geometry and stability of the base pairs they participate in.
doi:10.1093/nar/gkv606
PMCID: PMC4538814  PMID: 26117545
22.  NPDock: a web server for protein–nucleic acid docking 
Nucleic Acids Research  2015;43(Web Server issue):W425-W430.
Protein–RNA and protein–DNA interactions play fundamental roles in many biological processes. A detailed understanding of these interactions requires knowledge about protein–nucleic acid complex structures. Because the experimental determination of these complexes is time-consuming and perhaps futile in some instances, we have focused on computational docking methods starting from the separate structures. Docking methods are widely employed to study protein–protein interactions; however, only a few methods have been made available to model protein–nucleic acid complexes. Here, we describe NPDock (Nucleic acid–Protein Docking); a novel web server for predicting complexes of protein–nucleic acid structures which implements a computational workflow that includes docking, scoring of poses, clustering of the best-scored models and refinement of the most promising solutions. The NPDock server provides a user-friendly interface and 3D visualization of the results. The smallest set of input data consists of a protein structure and a DNA or RNA structure in PDB format. Advanced options are available to control specific details of the docking process and obtain intermediate results. The web server is available at http://genesilico.pl/NPDock.
doi:10.1093/nar/gkv493
PMCID: PMC4489298  PMID: 25977296
23.  Magnesium-binding architectures in RNA crystal structures: validation, binding preferences, classification and motif detection 
Nucleic Acids Research  2015;43(7):3789-3801.
The ubiquitous presence of magnesium ions in RNA has long been recognized as a key factor governing RNA folding, and is crucial for many diverse functions of RNA molecules. In this work, Mg2+-binding architectures in RNA were systematically studied using a database of RNA crystal structures from the Protein Data Bank (PDB). Due to the abundance of poorly modeled or incorrectly identified Mg2+ ions, the set of all sites was comprehensively validated and filtered to identify a benchmark dataset of 15 334 ‘reliable’ RNA-bound Mg2+ sites. The normalized frequencies by which specific RNA atoms coordinate Mg2+ were derived for both the inner and outer coordination spheres. A hierarchical classification system of Mg2+ sites in RNA structures was designed and applied to the benchmark dataset, yielding a set of 41 types of inner-sphere and 95 types of outer-sphere coordinating patterns. This classification system has also been applied to describe six previously reported Mg2+-binding motifs and detect them in new RNA structures. Investigation of the most populous site types resulted in the identification of seven novel Mg2+-binding motifs, and all RNA structures in the PDB were screened for the presence of these motifs.
doi:10.1093/nar/gkv225
PMCID: PMC4402538  PMID: 25800744
24.  Brickworx builds recurrent RNA and DNA structural motifs into medium- and low-resolution electron-density maps 
A computer program that builds crystal structure models of nucleic acid molecules is presented. It can be accessed at http://iimcb.genesilico.pl/brickworx.
Brickworx is a computer program that builds crystal structure models of nucleic acid molecules using recurrent motifs including double-stranded helices. In a first step, the program searches for electron-density peaks that may correspond to phosphate groups; it may also take into account phosphate-group positions provided by the user. Subsequently, comparing the three-dimensional patterns of the P atoms with a database of nucleic acid fragments, it finds the matching positions of the double-stranded helical motifs (A-RNA or B-DNA) in the unit cell. If the target structure is RNA, the helical fragments are further extended with recurrent RNA motifs from a fragment library that contains single-stranded segments. Finally, the matched motifs are merged and refined in real space to find the most likely conformations, including a fit of the sequence to the electron-density map. The Brickworx program is available for download and as a web server at http://iimcb.genesilico.pl/brickworx.
doi:10.1107/S1399004715000383
PMCID: PMC4356372  PMID: 25760616
Brickworx; model building; nucleic acids
25.  Sequence-specific cleavage of dsRNA by Mini-III RNase 
Nucleic Acids Research  2015;43(5):2864-2873.
Ribonucleases (RNases) play a critical role in RNA processing and degradation by hydrolyzing phosphodiester bonds (exo- or endonucleolytically). Many RNases that cut RNA internally exhibit substrate specificity, but their target sites are usually limited to one or a few specific nucleotides in single-stranded RNA and often in a context of a particular three-dimensional structure of the substrate. Thus far, no RNase counterparts of restriction enzymes have been identified which could cleave double-stranded RNA (dsRNA) in a sequence-specific manner. Here, we present evidence for a sequence-dependent cleavage of long dsRNA by RNase Mini-III from Bacillus subtilis (BsMiniIII). Analysis of the sites cleaved by this enzyme in limited digest of bacteriophage Φ6 dsRNA led to the identification of a consensus target sequence. We defined nucleotide residues within the preferred cleavage site that affected the efficiency of the cleavage and were essential for the discrimination of cleavable versus non-cleavable dsRNA sequences. We have also determined that the loop α5b-α6, a distinctive structural element in Mini-III RNases, is crucial for the specific cleavage, but not for dsRNA binding. Our results suggest that BsMiniIII may serve as a prototype of a sequence-specific dsRNase that could possibly be used for targeted cleavage of dsRNA.
doi:10.1093/nar/gkv009
PMCID: PMC4357697  PMID: 25634891

Results 1-25 (122)