The impact of various seed, accessibility and interaction constraints on sRNA target prediction- a systematic assessment.

BACKGROUND: Seed and accessibility constraints are core features to enable highly accurate sRNA target screens based on RNA-RNA interaction prediction. Currently, available tools provide different (sets of) constraints and default parameter sets. Thus, it is hard to impossible for users to estimate the influence of individual restrictions on the prediction results. RESULTS: Here, we present a systematic assessment of the impact of established and new constraints on sRNA target prediction both on a qualitative as well as computational level. This is done exemplarily based on the performance of IntaRNA, one of the most exact sRNA target prediction tools. IntaRNA provides various ways to constrain considered seed interactions, e.g. based on seed length, its accessibility, minimal unpaired probabilities, or energy thresholds, beside analogous constraints for the overall interaction. Thus, our results reveal the impact of individual constraints and their combinations. CONCLUSIONS: This provides both a guide for users what is important and recommendations for existing and upcoming sRNA target prediction approaches.We show on a large sRNA target screen benchmark data set that only by altering the parameter set, IntaRNA recovers 30% more verified interactions while becoming 5-times faster. This exemplifies the potential of seed, accessibility and interaction constraints for sRNA target prediction.

Probing RNA structure and interaction dynamics at the single molecule level.

RNA structures and their dynamic fluctuations lie at the heart of understanding key biological process such as transcription, splicing, translation and RNA decay. While conventional bulk assays have proven to identify and characterize key pathway intermediates, the generally dynamic nature of RNA structures renders the information obtained from time and ensemble averaging techniques necessarily lacking in critical details. Here we detail Single-Molecule Kinetic Analysis of RNA Transient Structure (SiM-KARTS), a method that readily monitors structural fluctuations of single RNA molecules through the repetitive interaction of fluorescent probes with an unlabeled, surface-immobilized RNA target of virtually any length and in any biological context. In addition, we demonstrate the broad applicability of SiM-KARTS by kinetically fingerprinting the binding of cognate tRNA ligand to single immobilized T-box riboswitch molecules. SiM-KARTS represents a valuable tool for probing biologically relevant structure and interaction features of potentially many diverse RNA metabolic pathways.

Conformational Dynamics of thiM Riboswitch To Understand the Gene Regulation Mechanism Using Markov State Modeling and the Residual Fluctuation Network Approach.

Thiamine pyrophosphate (TPP) riboswitch is a cis-regulatory element in the noncoding region of mRNA. The aptamer domain of TPP riboswitch detects the high abundance of coenzyme thiamine pyrophosphate (TPP) and modulates the gene expression for thiamine synthetic gene. The mechanistic understanding in recognition of TPP in aptamer domain and ligand-induced compactness for folding of expression platform are most important to designing novel modulators. To understand the dynamic behavior of TPP riboswitch upon TPP binding, molecular dynamics simulations were performed for 400 ns in both apo and TPP bound forms of thiM riboswitch from E. coli and analyzed in terms of eRMSD-based Markov state modeling and residual fluctuation network. Markov state models show good correlations in transition probability among metastable states from simulated trajectory and generated models. Structural compactness in TPP bound form is observed which is correlated with SAXS experiment. The importance of junction of P4 and P5 is evident during dynamics, which correlates with FRET analysis. The dynamic nature of two sensor forearms is due to the flexible P1 helix, which is its intrinsic property. The transient state in TPP-bound form was observed in the Markov state model, along with stable states. We believe that this transient state is responsible to assist the influx and outflux of ligand molecule by creating a solvent channel around the junction region of P4 and P5 and such a structure was anticipated in FRET analysis. The dynamic nature of riboswitch is dependent on the interaction between residues on distal loops L3 and L5/P3 and junction P4 and P5, J3/2 which stabilize the J2/4. It helps in the transfer of allosteric information between J2/4 and P3/L5 tertiary docking region through the active site residues. Understanding such information flow will benefit in highlighting crucial residues in highly dynamic and kinetic systems. Here, we report the residues and segments in riboswitch that play vital roles in providing stability and this can be exploited in designing inhibitors to regulate the functioning of riboswitches.

Assessment of External Guide Sequences' (EGS) Efficiency as Inducers of RNase P-Mediated Cleavage of mRNA Target Molecules.

RNase P is a ribozyme consisting of a catalytic RNA molecule and, depending on the organism, one or more cofactor proteins. It was initially identified as the enzyme that mediates cleavage of precursor tRNAs at the 5'-end termini to generate the mature tRNAs. An important characteristic of RNase P is that its specificity depends on the structure rather than the sequence of the RNA substrate. Any RNA species that interacts with an antisense molecule (called external guide sequence, EGS) and forms the appropriate structure can be cleaved by RNase P. This property is the basis for EGS technology, an antisense methodology for inhibiting gene expression by eliciting RNase P-mediated cleavage of a target mRNA molecule. EGS technology is being developed to design therapies against a large variety of diseases. An essential milestone in developing EGSs as therapies is the assessment of the efficiency of antisense molecules to induce cleavage of the target mRNA and evaluate their effect in vivo. Here, we describe simple protocols to test the ability of EGSs to induce cleavage of a target mRNA in vitro and to induce a phenotypic change in growing cells.

Interactions between horizontally acquired genes create a fitness cost in Pseudomonas aeruginosa.

Horizontal gene transfer (HGT) plays a key role in bacterial evolution, especially with respect to antibiotic resistance. Fitness costs associated with mobile genetic elements (MGEs) are thought to constrain HGT, but our understanding of these costs remains fragmentary, making it difficult to predict the success of HGT events. Here we use the interaction between P. aeruginosa and a costly plasmid (pNUK73) to investigate the molecular basis of the cost of HGT. Using RNA-Seq, we show that the acquisition of pNUK73 results in a profound alteration of the transcriptional profile of chromosomal genes. Mutations that inactivate two genes encoded on chromosomally integrated MGEs recover these fitness costs and transcriptional changes by decreasing the expression of the pNUK73 replication gene. Our study demonstrates that interactions between MGEs can compromise bacterial fitness via altered gene expression, and we argue that conflicts between mobile elements impose a general constraint on evolution by HGT.

Molecular cartography of the human skin surface in 3D.

The human skin is an organ with a surface area of 1.5-2 m(2) that provides our interface with the environment. The molecular composition of this organ is derived from host cells, microbiota, and external molecules. The chemical makeup of the skin surface is largely undefined. Here we advance the technologies needed to explore the topographical distribution of skin molecules, using 3D mapping of mass spectrometry data and microbial 16S rRNA amplicon sequences. Our 3D maps reveal that the molecular composition of skin has diverse distributions and that the composition is defined not only by skin cells and microbes but also by our daily routines, including the application of hygiene products. The technological development of these maps lays a foundation for studying the spatial relationships of human skin with hygiene, the microbiota, and environment, with potential for developing predictive models of skin phenotypes tailored to individual health.

Estimation of GFP-tagged RNA numbers from temporal fluorescence intensity data.

MOTIVATION: MS2-GFP-tagging of RNA is currently the only method to measure intervals between consecutive transcription events in live cells. For this, new transcripts must be accurately detected from intensity time traces. RESULTS: We present a novel method for automatically estimating RNA numbers and production intervals from temporal data of cell fluorescence intensities that reduces uncertainty by exploiting temporal information. We also derive a robust variant, more resistant to outliers caused e.g. by RNAs moving out of focus. Using Monte Carlo simulations, we show that the quantification of RNA numbers and production intervals is generally improved compared with previous methods. Finally, we analyze data from live Escherichia coli and show statistically significant differences to previous methods. The new methods can be used to quantify numbers and production intervals of any fluorescent probes, which are present in low copy numbers, are brighter than the cell background and degrade slowly. AVAILABILITY: Source code is available under Mozilla Public License at http://www.cs.tut.fi/%7ehakkin22/jumpdet/.

Comparison of extracellular DNase- and protease-producing spoilage bacteria isolated from Delaware pond-sourced and retail channel catfish (Ictalurus punctatus).

BACKGROUND: Spoilage of fishery products begins immediately following filleting due to microbial growth that degrades fish tissue quality prior to consumption. Extensive research has been conducted to identify such bacterial populations. A better understanding of the mechanisms involved in fish spoilage is necessary as a novel remedy for microbial spoilage inhibition has yet to be established for fish tissue. The present study identified, for the first time, bacterial populations that produce extracellular DNase and protease from Delaware and local retail distributed channel catfish (Ictalurus punctatus) fillets. RESULTS: A clear trend was identified between bacteria derived from catfish filleted under aseptic conditions where Pseudomonas was the dominant genus. Bacteria isolated from retail catfish contained high quantities of DNase-producing isolates, in contrast to aseptic-filleted catfish tissue which had none. Both types of catfish sample maintained high populations of protease-producing bacterial colonies throughout the duration of the study. Most bacteria isolated from catfish intestines exhibited DNase production with no protease production. CONCLUSION: Specific spoilage organism populations were significantly higher on retail-derived catfish in comparison to lab-filleted Delaware cultured catfish tissue. It is suggested that DNase production and protease production contribute to the spoilage of fish tissue as a result of mishandling and septic filleting being the major cause of rapid catfish tissue spoilage.

Structural model of an mRNA in complex with the bacterial chaperone Hfq.

The Sm-like protein Hfq (host factor Q-beta phage) facilitates regulation by bacterial small noncoding RNAs (sRNAs) in response to stress and other environmental signals. Here, we present a low-resolution model of Escherichia coli Hfq bound to the rpoS mRNA, a bacterial stress response gene that is targeted by three different sRNAs. Selective 2'-hydroxyl acylation and primer extension, small-angle X-ray scattering, and Monte Carlo molecular dynamics simulations show that the distal face and lateral rim of Hfq interact with three sites in the rpoS leader, folding the RNA into a compact tertiary structure. These interactions are needed for sRNA regulation of rpoS translation and position the sRNA target adjacent to an sRNA binding region on the proximal face of Hfq. Our results show how Hfq specifically distorts the structure of the rpoS mRNA to enable sRNA base pairing and translational control.

A universal TaqMan-based RT-PCR protocol for cost-efficient detection of small noncoding RNA.

Several methods for the detection of RNA have been developed over time. For small RNA detection, a stem-loop reverse primer-based protocol relying on TaqMan RT-PCR has been described. This protocol requires an individual specific TaqMan probe for each target RNA and, hence, is highly cost-intensive for experiments with small sample sizes or large numbers of different samples. We describe a universal TaqMan-based probe protocol which can be used to detect any target sequence and demonstrate its applicability for the detection of endogenous as well as artificial eukaryotic and bacterial small RNAs. While the specific and the universal probe-based protocol showed the same sensitivity, the absolute sensitivity of detection was found to be more than 100-fold lower for both than previously reported. In subsequent experiments, we found previously unknown limitations intrinsic to the method affecting its feasibility in determination of mature template RISC incorporation as well as in multiplexing. Both protocols were equally specific in discriminating between correct and incorrect small RNA targets or between mature miRNA and its unprocessed RNA precursor, indicating the stem-loop RT-primer, but not the TaqMan probe, triggers target specificity. The presented universal TaqMan-based RT-PCR protocol represents a cost-efficient method for the detection of small RNAs.

Ion Torrent-based transcriptional assessment of a Corynebacterium pseudotuberculosis equi strain reveals denaturing high-performance liquid chromatography a promising rRNA depletion method.

Corynebacterium pseudotuberculosis equi is a Gram-positive pathogenic bacterium which affects a variety of hosts. Besides the great economic losses it causes to horse-breeders, this organism is also known to be an important infectious agent to cattle and buffaloes. As an outcome of the efforts in characterizing the molecular basis of its virulence, several complete genome sequences were made available in recent years, enabling the large-scale assessment of genes throughout distinct isolates. Meanwhile, the RNA-seq stood out as the technology of choice for comprehensive transcriptome studies, which may bring valuable information regarding active genomic regions, despite of the still impeditive associated costs. In an attempt to increase the use of generated reads per instrument run, by effectively eliminating unwanted rRNAs from total RNA samples without relying on any commercially available kits, we applied denaturing high-performance liquid chromatography (DHPLC) as an alternative method to assess the transcriptional profile of C. pseudotuberculosis. We have found that the DHPLC depletion method, allied to Ion Torrent sequencing, allows mapping of transcripts in a comprehensive way and identifying novel transcripts when a de novo approach is used. These data encourage us to use DHPLC in future transcriptional evaluations in C. pseudotuberculosis.

Genotypic characterization of Brochothrix thermosphacta isolated during storage of minced pork under aerobic or modified atmosphere packaging conditions.

A total of 306 colonies were isolated from the selective medium for Brochothrix spp., during the spoilage of minced pork stored at 0, 5, 10 and 15 °C and packed aerobically and under modified atmosphere packaging conditions (MAP). Brochothrix biodiversity was assessed by Pulsed Field Gel Electrophoresis (PFGE), and representative strains were further analysed by Rep-PCR using primer (GTG)5 and Sau-PCR with primers SAG1 and SAG2. Although, different results were obtained from the different methods, a significant diversity among isolates recovered from aerobic conditions was observed. On the contrary, isolates from MAP showed a lower degree of heterogeneity. The storage conditions affected the Brochothrix diversity, the strains isolated in the initial stage being different from the ones present at the final stage of storage at chill temperatures. A representative number of isolates, based on the results of the clustering by molecular methods, were subjected to 16S rRNA gene sequencing, revealing that all belonged to Brochothrix thermosphacta.

Global parameter estimation methods for stochastic biochemical systems.

BACKGROUND: The importance of stochasticity in cellular processes having low number of molecules has resulted in the development of stochastic models such as chemical master equation. As in other modelling frameworks, the accompanying rate constants are important for the end-applications like analyzing system properties (e.g. robustness) or predicting the effects of genetic perturbations. Prior knowledge of kinetic constants is usually limited and the model identification routine typically includes parameter estimation from experimental data. Although the subject of parameter estimation is well-established for deterministic models, it is not yet routine for the chemical master equation. In addition, recent advances in measurement technology have made the quantification of genetic substrates possible to single molecular levels. Thus, the purpose of this work is to develop practical and effective methods for estimating kinetic model parameters in the chemical master equation and other stochastic models from single cell and cell population experimental data. RESULTS: Three parameter estimation methods are proposed based on the maximum likelihood and density function distance, including probability and cumulative density functions. Since stochastic models such as chemical master equations are typically solved using a Monte Carlo approach in which only a finite number of Monte Carlo realizations are computationally practical, specific considerations are given to account for the effect of finite sampling in the histogram binning of the state density functions. Applications to three practical case studies showed that while maximum likelihood method can effectively handle low replicate measurements, the density function distance methods, particularly the cumulative density function distance estimation, are more robust in estimating the parameters with consistently higher accuracy, even for systems showing multimodality. CONCLUSIONS: The parameter estimation methodologies described in this work have provided an effective and practical approach in the estimation of kinetic parameters of stochastic systems from either sparse or dense cell population data. Nevertheless, similar to kinetic parameter estimation in other modelling frameworks, not all parameters can be estimated accurately, which is a common problem arising from the lack of complete parameter identifiability from the available data.

Loop-mediated isothermal amplification (LAMP) for the rapid and sensitive detection of Salmonella Typhimurium from pork.

UNLABELLED: Loop-mediated isothermal amplification (LAMP) is a novel molecular detection method that is more rapid and simpler than PCR. Products can be detected by turbidity using one temperature without the need for expensive PCR equipment. Our objective was to sensitively detect Salmonella Typhimurium from pork products within 1 d using the LAMP assay. Pork chop and pork sausage samples (25 g) were inoculated with high (10(8) to 10(6) CFU) and low (10(5) to 10(0) CFU) inocula of S. Typhimurium. Serial dilutions in phosphate buffered saline were plated on XLT4 agar either immediately or after selective preenrichment in tetrathionate broth (225 mL) for 10-h at 37 degrees C. Nucleic acid was extracted using the TRIzol method from 1-mL samples. The LAMP assay using 6 specific invA gene primers and Bst DNA polymerase reaction mix was carried out at 62 degrees C for 90 min in a waterbath. Turbid products were detected visually and by agarose gel electrophoresis. Improved Salmonella detection at 10(2) CFU/25 g for both pork chop and sausage was obtained after 10-h enrichment and 10(6) CFU/25 g without enrichment for both products. This assay can detect Salmonella from pork within 1 d, significantly faster than traditional methods that take >5 d. This method shows tremendous potential for routine diagnostics and monitoring of Salmonella by the pork industry. PRACTICAL APPLICATION: The novel loop-mediated isothermal amplification (LAMP) assay is a rapid, specific, and sensitive method that has potential application for routine diagnostics of Salmonella from pork products. The isothermal method does not require expensive equipment such as a PCR thermocycler but only a simple waterbath for amplification within 90 min. Detection is even simpler by visual eye or turbidimeters that are less expensive than fluorescent spectrophotometers or real-time PCR machines. All these advantages make it a practical approach for routine use by processing industries to rapidly detect Salmonella in their environment and to implement appropriate control strategies. To improve detection sensitivities, preenrichment followed by selective enrichment may be necessary. Even so, the entire assay can be completed at the most within two 8-h working shifts.

Disparity between changes in mRNA abundance and enzyme activity in Corynebacterium glutamicum: implications for DNA microarray analysis.

The relationship between changes in mRNA abundance and enzyme activity was determined for three genes over a span of nearly 3 h during amino acid production in Corynebacterium glutamicum. Gene expression changes during C. glutamicum fermentations were examined by complementary DNA (cDNA) microarrays and by a second method for quantitating RNA levels, competitive reverse transcriptase-PCR (RT-PCR). The results obtained independently by both methods were compared and found to be in agreement, thus validating the quantitative potential of DNA microarrays for gene expression profiling. Evidence of a disparity between mRNA abundance and enzyme activity is presented and supports our belief that it is difficult to generally predict protein activity from quantitative transcriptome data. Homoserine dehydrogenase, threonine dehydratase, and homoserine kinase are enzymes involved in the biosynthesis of l-isoleucine and other aspartate-derived amino acids in C. glutamicum. Our data suggest that different underlying regulatory mechanisms may be connected with the expression of the genes encoding each of these three enzymes. Indeed, whereas in one case the increases in enzyme activity exceeded those in the corresponding mRNA abundance, in another case large increases in the levels of gene expression were not congruent with changes in enzyme activity.

Conserved economics of transcription termination in eubacteria.

A secondary structure in the nascent RNA followed by a trail of U residues is believed to be necessary and sufficient to terminate transcription. Such structures represent an extremely economical mechanism of transcription termination since they function in the absence of any additional protein factors. We have developed a new algorithm, GeSTer, to identify putative terminators and analysed all available complete bacterial genomes. The algorithm classifies the structures into five classes. We find that potential secondary structure sequences are concentrated downstream of coding regions in most bacterial genomes. Interestingly, many of these structures are not followed by a discernible U-trail. However, irrespective of the nature of the trail sequence, the structures show a similar distribution, indicating that they serve the same purpose. In contrast, such a distribution is absent in archaeal genomes, indicating that they employ a distinct mechanism for transcription termination. The present algorithm represents the fastest and most accurate algorithm for identifying terminators in eubacterial genomes without being restricted by the classical Escherichia coli paradigm.

Acquisition of novel catalytic activity by the M1 RNA ribozyme: the cost of molecular adaptation.

The ribonucleoprotein RNase P is a critical component of metabolism in all known organisms. In Escherichia coli, RNase P processes a vast array of substrates, including precursor-tRNAs and precursor 4. 5S RNA. In order to understand how such catalytic versatility is achieved and how novel catalytic activity can be acquired, we evolve the M1 RNA ribozyme (the catalytic component of E. coli RNase P) in vitro for cleavage of a DNA substrate. In so doing, we probe the consequences of enhancing catalytic activity on a novel substrate and investigate the cost this versatile enzyme pays for molecular adaptation. A total of 25 generations of in vitro evolution yield a population showing more than a 1000-fold increase in DNA substrate cleavage efficiency (kcat/KM) relative to wild-type M1 RNA. This enhancement is accompanied by a significant reduction in the ability of evolved ribozymes to process the ptRNA class of substrates but also a contrasting increase in activity on the p4.5S RNA class of substrates. This change in the catalytic versatility of the evolved ribozymes suggests that the acquired activity comes at the cost of substrate versatility, and indicates that E. coli RNase P catalytic flexibility is maintained in vivo by selection for the processing of multiple substrates. M1 RNA derivatives enhance cleavage of the DNA substrate by accelerating the catalytic step (kcat) of DNA cleavage, although overall processing efficiency is offset by reduced substrate binding. The enhanced ability to cleave a DNA substrate cannot be readily traced to any of the predominant mutations found in the evolved population, and must instead be due to multiple sequence changes dispersed throughout the molecule. This conclusion underscores the difficulty of correlating observed mutations with changes in catalytic behavior, even in simple biological catalysts for which three-dimensional models are available.

The proofreading domain of Escherichia coli DNA polymerase I and other DNA and/or RNA exonuclease domains.

Prior sequence analysis studies have suggested that bacterial ribonuclease (RNase) Ds comprise a complete domain that is found also in Homo sapiens polymyositis-scleroderma overlap syndrome 100 kDa autoantigen and Werner syndrome protein. This RNase D 3'-->5' exoribonuclease domain was predicted to have a structure and mechanism of action similar to the 3'-->5' exodeoxyibonuclease (proofreading) domain of DNA polymerases. Here, hidden Markov model (HMM) and phylogenetic studies have been used to identify and characterise other sequences that may possess this exonuclease domain. Results indicate that it is also present in the RNase T family; Borrelia burgdorferi P93 protein, an immunodominant antigen in Lyme disease; bacteriophage T4 dexA and Escherichia coli exonuclease I, processive 3'-->5' exodeoxyribonucleases that degrade single-stranded DNA; Bacillus subtilis dinG, a probable helicase involved in DNA repair and possibly replication, and peptide synthase 1; Saccharomyces cerevisiae Pab1p-dependent poly(A) nuclease PAN2 subunit, required for shortening mRNA poly(A) tails; Caenorhabditis elegans and Mus musculus CAF1, transcription factor CCR4-associated factor 1; Xenopus laevis XPMC2, prevention of mitotic catastrophe in fission yeast; Drosophila melanogaster egalitarian, oocyte specification and axis determination, and exuperantia, establishment of oocyte polarity; H.sapiens HEM45, expressed in tumour cell lines and uterus and regulated by oestrogen; and 31 open reading frames including one in Methanococcus jannaschii . Examination of a multiple sequence alignment and two three-dimensional structures of proofreading domains has allowed definition of the core sequence, structural and functional elements of this exonuclease domain.

Assessment of quantitative models for plasmid ColE1 copy number control.

Two quantitative models of plasmid ColE1 copy number control are compared with respect to mathematical logic of derivation and application to experimental observations. Explanatory background material and clarifications are supplied for selected aspects of each model. Contrasting features are emphasized and experiments are suggested to distinguish between predictions of the models.