Structure and function of the conserved 690 hairpin in Escherichia coli 16 S ribosomal RNA. II. NMR solution structure

Adjusting the Energy Profile for CH-O Interactions Leads to Improved Stability of RNA Stem-Loop Structures in MD Simulations

The role of ribonucleic acid (RNA) in biology continues to grow, but insight into important aspects of RNA behavior is lacking, such as dynamic structural ensembles in different environments, how flexibility is coupled to function, and how function might be modulated by small molecule binding. In the case of proteins, much progress in these areas has been made by complementing experiments with atomistic simulations, but RNA simulation methods and force fields are less mature. It remains challenging to generate stable RNA simulations, even for small systems where well-defined, thermostable structures have been established by experiments. Many different aspects of RNA energetics have been adjusted in force fields, seeking improvements that are transferable across a variety of RNA structural motifs. In this work, the role of weak CH···O interactions is explored, which are ubiquitous in RNA structure but have received less attention in RNA force field development. By comparing data extracted from high-resolution RNA crystal structures to energy profiles from quantum mechanics and force field calculations, it is shown that CH···O interactions are overly repulsive in the widely used Amber RNA force fields. A simple, targeted adjustment of CH···O repulsion that leaves the remainder of the force field unchanged was developed. Then, the standard and modified force fields were tested using molecular dynamics (MD) simulations with explicit water and salt, amassing over 300 μs of data for multiple RNA systems containing important features such as the presence of loops, base stacking interactions as well as canonical and noncanonical base pairing. In this work and others, standard force fields lead to reproducible unfolding of the NMR-based structures. Including a targeted CH···O adjustment in an otherwise identical protocol dramatically improves the outcome, leading to stable simulations for all RNA systems tested.

Accurately Modeling RNA Stem-Loops in an Implicit Solvent Environment

Ribonucleic acid (RNA) molecules can adopt a variety of secondary and tertiary structures in solution, with stem-loops being one of the more common motifs. Here, we present a systematic analysis of 15 RNA stem-loop sequences simulated with molecular dynamics simulations in an implicit solvent environment. Analysis of RNA cluster ensembles showed that the stem-loop structures can generally adopt the A-form RNA in the stem region. Loop structures are more sensitive, and experimental structures could only be reproduced with modification of CH···O interactions in the force field, combined with an implicit solvent nonpolar correction to better model base stacking interactions. Accurately modeling RNA with current atomistic physics-based models remains challenging, but the RNA systems studied herein may provide a useful benchmark set for testing other RNA modeling methods in the future.

Determining the most accurate 16S rRNA hypervariable region for taxonomic identification from respiratory samples

16S rRNA gene profiling, which contains nine hypervariable regions (V1-V9), is the gold standard for identifying taxonomic units by high-throughput sequencing. Microbiome studies combine two or more region sequences (usually V3-V4) to increase the resolving power for identifying bacterial taxa. We compare the resolving powers of V1-V2, V3-V4, V5-V7, and V7-V9 to improve microbiome analyses in sputum samples from patients with chronic respiratory diseases. DNA were isolated from 33 human sputum samples, and libraries were created using a QIASeq screening panel intended for Illumina platforms (16S/ITS; Qiagen Hilden, Germany). The analysis included a mock community as a microbial standard control (ZymoBIOMICS). We used the Deblur algorithm to identify bacterial amplicon sequence variants (ASVs) at the genus level. Alpha diversity was significantly higher for V1-V2, V3-V4, and V5-V7 compared with V7-V9, and significant compositional dissimilarities in the V1-V2 and V7-V9 analyses versus the V3-V4 and V5-V7 analyses. A cladogram confirmed these compositional differences, with the latter two being very similar in composition. The combined hypervariable regions showed significant differences when discriminating between the relative abundances of bacterial genera. The area under the curve revealed that V1-V2 had the highest resolving power for accurately identifying respiratory bacterial taxa from sputum samples. Our study confirms that 16S rRNA hypervariable regions provide significant differences for taxonomic identification in sputum. Comparing the taxa of microbial community standard control with the taxa samples, V1-V2 combination exhibits the most sensitivity and specificity. Thus, while third generation full-length 16S rRNA sequencing platforms become more available, the V1-V2 hypervariable regions can be used for taxonomic identification in sputum.

Sediment-associated bacterial community and predictive functionalities are influenced by choice of 16S ribosomal RNA hypervariable region(s): An amplicon-based diversity study

Meta-omics approaches such as high-throughput sequencing of 16S hypervariable region(s) [HVR(s)] is extensively applied for profiling microbial community. Several studies have deciphered the influence of HVR(s) on bacterial diversity; most of these were devoted to human body habitats. Extent to which targeted HVR(s) influences the diversity estimates of environmental samples is rather unclear. Here, we evaluated the performance of five widely used universal primer pairs spanning V1-V3, V3-V4, V4, V5-V6 and V7-V9 HVRs to characterize bacterial diversity and predictive functionality of complex marine sediments. Obtained results revealed that the HVR(s) V4 and V5-V6 represented the higher species richness than others while, V1-V3 and V7-V9 were unsuccessful to detect Bacteroidetes and Planctomycetes. Further, PICRUSt analysis showed that the selected HVR(s) also had significant impact on the predictive functional profile. Conclusively, this study proved that HVR selection has a profound effect on overall results and thus should be selected with utmost caution.

The effect of 16S rRNA region choice on bacterial community metabarcoding results

In this work, we compare the resolution of V2-V3 and V3-V4 16S rRNA regions for the purposes of estimating microbial community diversity using paired-end Illumina MiSeq reads, and show that the fragment, including V2 and V3 regions, has higher resolution for lower-rank taxa (genera and species). It allows for a more precise distance-based clustering of reads into species-level OTUs. Statistically convergent estimates of the diversity of major species (defined as those that together are covered by 95% of reads) can be achieved at the sample sizes of 10000 to 15000 reads. The relative error of the Shannon index estimate for this condition is lower than 4%.

Sensitivity and correlation of hypervariable regions in 16S rRNA genes in phylogenetic analysis

Background

Prokaryotic 16S ribosomal RNA (rRNA) sequences are widely used in environmental microbiology and molecular evolution as reliable markers for the taxonomic classification and phylogenetic analysis of microbes. Restricted by current sequencing techniques, the massive sequencing of 16S rRNA gene amplicons encompassing the full length of genes is not yet feasible. Thus, the selection of the most efficient hypervariable regions for phylogenetic analysis and taxonomic classification is still debated. In the present study, several bioinformatics tools were integrated to build an in silico pipeline to evaluate the phylogenetic sensitivity of the hypervariable regions compared with the corresponding full-length sequences.

Results

The correlation of seven sub-regions was inferred from the geodesic distance, a parameter that is applied to quantitatively compare the topology of different phylogenetic trees constructed using the sequences from different sub-regions. The relationship between different sub-regions based on the geodesic distance indicated that V4-V6 were the most reliable regions for representing the full-length 16S rRNA sequences in the phylogenetic analysis of most bacterial phyla, while V2 and V8 were the least reliable regions.

Conclusions

Our results suggest that V4-V6 might be optimal sub-regions for the design of universal primers with superior phylogenetic resolution for bacterial phyla. A potential relationship between function and the evolution of 16S rRNA is also discussed.

Electronic supplementary material

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

RNA-PAIRS: RNA probabilistic assignment of imino resonance shifts

The significant biological role of RNA has further highlighted the need for improving the accuracy, efficiency and the reach of methods for investigating RNA structure and function. Nuclear magnetic resonance (NMR) spectroscopy is vital to furthering the goals of RNA structural biology because of its distinctive capabilities. However, the dispersion pattern in the NMR spectra of RNA makes automated resonance assignment, a key step in NMR investigation of biomolecules, remarkably challenging. Herein we present RNA Probabilistic Assignment of Imino Resonance Shifts (RNA-PAIRS), a method for the automated assignment of RNA imino resonances with synchronized verification and correction of predicted secondary structure. RNA-PAIRS represents an advance in modeling the assignment paradigm because it seeds the probabilistic network for assignment with experimental NMR data, and predicted RNA secondary structure, simultaneously and from the start. Subsequently, RNA-PAIRS sets in motion a dynamic network that reverberates between predictions and experimental evidence in order to reconcile and rectify resonance assignments and secondary structure information. The procedure is halted when assignments and base-parings are deemed to be most consistent with observed crosspeaks. The current implementation of RNA-PAIRS uses an initial peak list derived from proton-nitrogen heteronuclear multiple quantum correlation ((1)H-(15)N 2D HMQC) and proton-proton nuclear Overhauser enhancement spectroscopy ((1)H-(1)H 2D NOESY) experiments. We have evaluated the performance of RNA-PAIRS by using it to analyze NMR datasets from 26 previously studied RNAs, including a 111-nucleotide complex. For moderately sized RNA molecules, and over a range of comparatively complex structural motifs, the average assignment accuracy exceeds 90%, while the average base pair prediction accuracy exceeded 93%. RNA-PAIRS yielded accurate assignments and base pairings consistent with imino resonances for a majority of the NMR resonances, even when the initial predictions are only modestly accurate. RNA-PAIRS is available as a public web-server at http://pine.nmrfam.wisc.edu/RNA/.

NMR structures of (rGCUGAGGCU)2 and (rGCGGAUGCU)2: probing the structural features that shape the thermodynamic stability of GA pairs

The NMR structures of [see text] and [see text] are reported. The internal loop, [see text], is about 2 kcal/mol more stable than [see text] at 37 degrees C. The duplexes assemble into similar global folds characterized by the formation of tandem sheared GA pairs. The different stabilities of the loops are accompanied by differences in the local structure of the closing GU pairs. In the [see text] internal loop, the GU pairs form canonical wobble configurations with two hydrogen bonds, whereas in [see text], the GU pairs form a single hydrogen bond involving the amino group, GH22, and the carbonyl group, UO4. This pairing is similar to the GU closing pair of the 690 hairpin loop found in E. coli 16S rRNA. The [see text] and [see text] structures reveal how the subtle interplay between stacking and hydrogen bonding determines sequence dependent conformation and thermodynamic stability. Thus, this work provides structural and thermodynamic benchmarks for theoreticians in the ongoing effort to understand the sequence dependence of RNA physicochemical properties.

Combinatorial Genetic Technology for the Development of New Anti-infectives

Context.—We previously developed a novel technology known as instant evolution for high-throughput analysis of mutations in Escherichia coli ribosomal RNA.

Objective.—To develop a genetic platform for the isolation of new classes of antiinfectives that are not susceptible to drug resistance based on the instant evolution system.

Design.—Mutation libraries were constructed in the 16S rRNA gene of E coli and analyzed. In addition, the rRNA genes from a number of pathogenic bacteria were cloned and expressed in E coli. The 16S rRNA genes were incorporated into the instant-evolution system in E coli.

Setting.—The Department of Biological Sciences, Wayne State University, Detroit, Mich.

Main Outcome Measures.—Ribosome function was assayed by measuring the amount of green fluorescent protein produced by ribosomes containing mutant or foreign RNA in vivo.

Results.—We have developed a new combinatorial genetic technology (CGT) platform that allows high-throughput in vivo isolation and analysis of rRNA mutations that might lead to drug resistance. This information is being used to develop anti-infectives that recognize the wild type and all viable mutants of the drug target. CGT also provides a novel mechanism for identifying new drug targets.

Conclusions.—Antimicrobials produced using CGT will provide new therapies for the treatment of infections caused by human pathogens that are resistant to current antibiotics. The new therapeutics will be less susceptible to de novo resistance because CGT identifies all mutations of the target that might lead to resistance during the earliest stages of the drug discovery process.

Interactions of yeast ribosomal protein rpS14 with RNA

Yeast ribosomal protein S14 (rpS14) binds to two different RNA molecules: (1). helix 23 of 18S rRNA during its assembly into 40S ribosomal subunits and (2). a stem-loop structure in RPS14B pre-mRNA to repress expression of the RPS14B gene. We used the three-dimensional structure of Thermus thermophilus ribosomal protein S11, a bacterial homologue of rpS14, as a guide to identify conserved, surface-exposed amino acid residues that are likely to contact RNA. Eight residues that met these criteria were mutated to alanine. Most of these mutations affected interaction of rpS14 with either helix 23 or the RPS14B stem-loop RNA or both. Assembly of 40S ribosomal subunits and repression of RPS14B were also affected. S11 contains an extended carboxy-terminal domain rich in basic amino acids, which interacts with rRNA. We systematically evaluated the importance of each of the last ten amino acid residues in the basic, carboxy-terminal tail of yeast rpS14 for binding to RNA, by mutating each to alanine. Mutations in nine of these residues decreased binding of rpS14 to one or both of its RNA ligands. In addition, we examined the importance of four structural motifs in helix 23 of 18S rRNA for binding to rpS14. Mutations that altered either the terminal loop, the G-U base-pair closing the terminal loop, or the internal loop affected binding of rpS14 to helix 23.

Unusual DNA duplex and hairpin motifs

Single-stranded DNA or double-stranded DNA has the potential to adopt a wide variety of unusual duplex and hairpin motifs in the presence (trans) or absence (cis) of ligands. Several principles for the formation of those unusual structures have been established through the observation of a number of recurring structural motifs associated with different sequences. These include: (i) internal loops of consecutive mismatches can occur in a B-DNA duplex when sheared base pairs are adjacent to each other to confer extensive cross- and intra-strand base stacking; (ii) interdigitated (zipper-like) duplex structures form instead when sheared G*A base pairs are separated by one or two pairs of purine*purine mismatches; (iii) stacking is not restricted to base, deoxyribose also exhibits the potential to do so; (iv) canonical G*C or A.T base pairs are flexible enough to exhibit considerable changes from the regular H-bonded conformation. The paired bases become stacked when bracketed by sheared G.A base pairs, or become extruded out and perpendicular to their neighboring bases in the presence of interacting drugs; (v) the purine-rich and pyrimidine-rich loop structures are notably different in nature. The purine-rich loops form compact triloop structures closed by a sheared G*A, A*A, A*C or sheared-like G(anti)*C(syn) base pair that is stacked by a single residue. On the other hand, the pyrimidine-rich loops with a thymidine in the first position exhibit no base pairing but are characterized by the folding of the thymidine residue into the minor groove to form a compact loop structure. Identification of such diverse duplex or hairpin motifs greatly enlarges the repertoire for unusual DNA structural formation.

Modeling a minimal ribosome based on comparative sequence analysis

We have determined the three-dimensional organization of ribosomal RNAs and proteins essential for minimal ribosome function. Comparative sequence analysis identifies regions of the ribosome that have been evolutionarily conserved, and the spatial organization of conserved domains is determined by mapping these onto structures of the 30S and 50S subunits determined by X-ray crystallography. Several functional domains of the ribosome are conserved in their three-dimensional organization in the Archaea, Bacteria, Eucaryotic nuclear, mitochondria and chloroplast ribosomes. In contrast, other regions from both subunits have shifted their position in three-dimensional space during evolution, including the L11 binding domain and the alpha-sarcin-ricin loop (SRL). We examined conserved bridge interactions between the two ribosomal subunits, giving an indication of which contacts are more significant. The tRNA contacts that are conserved were also determined, highlighting functional interactions as the tRNA moves through the ribosome during protein synthesis. To augment these studies of a large collection of comparative structural models sampled from all major branches on the phylogenetic tree, Caenorhabditis elegans mitochondrial rRNA is considered individually because it is among the smallest rRNA sequences known. The C.elegans model supports the large collection of comparative structure models while providing insight into the evolution of mitochondrial ribosomes.

Functional studies of the 900 tetraloop capping helix 27 of 16S ribosomal RNA

The 900 tetraloop (positions 898-901) of Escherichia coli 16S rRNA caps helix 27, which is involved in a conformational switch crucial for the decoding function of the ribosome. This tetraloop forms a GNRA motif involved in intramolecular RNA-RNA interactions with its receptor in helix 24 of 16S rRNA. It is involved also in an intersubunit bridge, via an interaction with helix 67 in domain IV of 23S rRNA. Using a specialized ribosome system and an instant-evolution procedure, the four nucleotides of this loop were randomized and 15 functional mutants were selected in vivo. Positions 899 and 900, responsible for most of the tetraloop/receptor interactions, were found to be the most critical for ribosome activity. Functional studies showed that mutations in the 900 tetraloop impair subunit association and decrease translational fidelity. Computer modeling of the mutations allows correlation of the effect of mutations with perturbations of the tetraloop/receptor interactions.

Structure and function of the conserved 690 hairpin in Escherichia coli 16 S ribosomal RNA. III. Functional analysis of the 690 loop

An instant-evolution experiment was performed on the eight nucleotides comprising the loop region of the 690 hairpin in Escherichia coli 16 S ribosomal RNA. Positions 690 to 697 were randomly mutated and 101 unique functional mutants were isolated, sequenced and analyzed for function in vivo. Non-random nucleotide distributions were observed at each of the mutated positions except 693 and 694. Nucleotide identity significantly affected ribosome function at positions 690, 695, 696 and 697. Pyrimidines were absent at position 696 in the instant-evolution pool as were C at position 691 and G at position 697. A highly significant covariation was observed between nucleotides 690 and 697. No functional double mutants at positions 691 and 696 were obtained from the instant-evolution pool. In our NMR structure of the 690 loop, both the G690.U697 and G691.A696 form sheared hydrogen-bonded mismatches. To further examine the functional constraints between these paired nucleotides, one set of site-directed mutations was constructed at positions 690:697 and another set was constructed at positions 691:696. Functional analysis of the site-directed mutants is consistent with our instant-evolution findings and revealed constraints on the placement of specific functional groups observed in the NMR structure. Ten instant-evolution mutants were isolated that are more functional than the wild-type. Hyperactivity in these mutants correlates with a single mutation at position 693.