Thermodynamic characterization of the complete set of sequence symmetric tandem mismatches in RNA and an improved model for predicting the free energy contribution of sequence asymmetric tandem mismatches

A Test and Refinement of Folding Free Energy Nearest Neighbor Parameters for RNA Including N6-Methyladenosine

RNA folding free energy change parameters are widely used to predict RNA secondary structure and to design RNA sequences. These parameters include terms for the folding free energies of helices and loops. Although the full set of parameters has only been traditionally available for the four common bases and backbone, it is well known that covalent modifications of nucleotides are widespread in natural RNAs. Covalent modifications are also widely used in engineered sequences. We recently derived a full set of nearest neighbor terms for RNA that includes N6-methyladenosine (m6A). In this work, we test the model using 98 optical melting experiments, matching duplexes with or without N6-methylation of A. Most experiments place RRACH, the consensus site of N6-methylation, in a variety of contexts, including helices, bulge loops, internal loops, dangling ends, and terminal mismatches. For matched sets of experiments that include either A or m6A in the same context, we find that the parameters for m6A are as accurate as those for A. Across all experiments, the root mean squared deviation between estimated and experimental free energy changes is 0.67 kcal/mol. We used the new experimental data to refine the set of nearest neighbor parameter terms for m6A. These parameters enable prediction of RNA secondary structures including m6A, which can be used to model how N6-methylation of A affects RNA structure.

Thermodynamic characterization of naturally occurring RNA pentaloops

RNA folding is hierarchical; therefore, predicting RNA secondary structure from sequence is an intermediate step in predicting tertiary structure. Secondary structure prediction is based on a nearest neighbor model using free energy minimization. To improve secondary structure prediction, all types of naturally occurring secondary structure motifs need to be thermodynamically characterized. However, not all secondary structure motifs are well characterized. Pentaloops, the second most abundant hairpin size, is one such uncharacterized motif. In fact, the current thermodynamic model used to predict the stability of pentaloops was derived from a small data set of pentaloops and from data for other hairpins of different sizes. Here, the most commonly occurring pentaloops were identified and optically melted. New experimental data for 22 pentaloop sequences were combined with previously published data for nine pentaloop sequences. Using linear regression, a pentaloop-specific model was derived. This new model is simpler and more accurate than the current model. The new experimental data and improved model can be incorporated into software that is used to predict RNA secondary structure from sequence.

Efficient Biosynthesis of R-(-)-Linalool through Adjusting the Expression Strategy and Increasing GPP Supply in Escherichia coli

R-(-)-linalool is widely used in pharmaceutical, agrochemical, and fragrance industries. However, plant extraction furnishes only limited and unstable R-(-)-linalool yields that do not satisfy market demand. Therefore, a sustainable yet efficient and productive method is urgently needed. To induce the R-(-)-linalool biosynthesis pathway in Escherichia coli, we expressed several heterologous (3R)-linalool synthases (LISs) and then chose a suitable LIS from Streptomyces clavuligerus (bLIS) for further study. The bLIS expression was markedly elevated by using optimized ribosomal binding sites and protein fusion tags. To increase the geranyl diphosphate content, we tested various alterations in prenyltransferases and their mutants. The final strain accumulated 100.1 and 1027.3 mg L^-1 R-(-)-linalool under shake flask and fed-batch fermentation conditions, respectively. The latter is the highest reported R-(-)-linalool yield to date. This work could lay theoretical and empirical foundations for engineering terpenoid pathways and optimizing other metabolic pathways.

Nuclear Magnetic Resonance Reveals That GU Base Pairs Flanking Internal Loops Can Adopt Diverse Structures

RNA thermodynamics play an important role in determining the two- and three-dimensional structures of RNA. Internal loops of the sequence 5'-GMNU/3'-UNMG are relatively unstable thermodynamically. Here, five duplexes with GU-flanked 2 × 2 nucleotide internal loops were structurally investigated to reveal determinants of their instability. The following internal loops were investigated: 5'-GCAU/3'-UACG, 5'-UUCG/3'-GCUU, 5'-GCUU/3'-UUCG, 5'-GUCU/3'-UCUG, and 5'-GCCU/3'-UCCG. Two-dimensional nuclear magnetic resonance spectra indicate the absence of GU wobble base pairing in 5'-GCUU/3'-UUCG, 5'-GUCU/3'-UCUG, and 5'-GCCU/3'-UCCG. The 5'-GCUU/3'-UUCG loop has an unusual conformation of the GU base pairs, in which U's O2 carbonyl forms a bifurcated hydrogen bond with G's amino and imino protons. The internal loop of 5'-GUCU/3'-UCUG displays a shifted configuration in which GC pairs flank a U-U pair and several U's are in fast exchange between positions inside and outside the helix. In contrast, 5'-GCAU/3'-UACG and 5'-UUCG/3'-GCUU both have the expected GU wobble base pairs flanking the internal loop. Evidently, GU base pairs flanking internal loops are more likely to display atypical structures relative to Watson-Crick base pairs flanking internal loops. This appears to be more likely when the G of the GU pair is 5' to the loop. Such unusual structures could serve as recognition elements for biological function and as benchmarks for structure prediction methods.

Stability of RNA duplexes containing inosine·cytosine pairs

Inosine is found naturally in the anticodon loop of tRNA, is a product of adenosine deaminases that act on RNA, and can be used in oligonucleotide probes or to investigate the role of the exocyclic amino group of guanosine. Although the thermodynamics of I·U pairs in RNA have been systematically studied [Wright, D. J., Rice, J. L., Yanker, D. M., and Znosko, B. M. (2007) Biochemistry 46, 4625-4634], the thermodynamics of I·C pairs in RNA have not. Here, we have performed optical melting experiments on a series of RNA duplexes containing I·C pairs and compared their thermodynamics to the same duplexes containing A·C and G-C pairs. Nearest neighbor parameters for single I·C pairs adjacent to Watson-Crick pairs were derived. The derived nearest neighbor parameters are compared to those previously predicted blindly through a reweighting of energy-function collection with conformational ensemble sampling in Rosetta [Chou, F.-C., Kladwang, W., Kappel, K., and Das, R. (2016) Proc. Natl. Acad. Sci. U.S.A. 113, 8430-8435]. Scientists can use these nearest neighbor parameters to calculate the stability of ADAR products and to calculate the stability of an RNA duplex in which G-to-I substitution was used to determine the role of the exocyclic amino group of G.

Analysis of RNA nearest neighbor parameters reveals interdependencies and quantifies the uncertainty in RNA secondary structure prediction

RNA secondary structure prediction is often used to develop hypotheses about structure-function relationships for newly discovered RNA sequences, to identify unknown functional RNAs, and to design sequences. Secondary structure prediction methods typically use a thermodynamic model that estimates the free energy change of possible structures based on a set of nearest neighbor parameters. These parameters were derived from optical melting experiments of small model oligonucleotides. This work aims to better understand the precision of structure prediction. Here, the experimental errors in optical melting experiments were propagated to errors in the derived nearest neighbor parameter values and then to errors in RNA secondary structure prediction. To perform this analysis, the optical melting experimental values were systematically perturbed within the estimates of experimental error and alternative sets of nearest neighbor parameters were then derived from these error-bounded values. Secondary structure predictions using either the perturbed or reference parameter sets were then compared. This work demonstrated that the precision of RNA secondary structure prediction is more robust than suggested by previous work based on perturbation of the nearest neighbor parameters. This robustness is due to correlations between parameters. Additionally, this work identified weaknesses in the parameter derivation that makes accurate assessment of parameter uncertainty difficult. Considerations for experimental design are provided to mitigate these weaknesses are provided.

Molecular dynamics correctly models the unusual major conformation of the GAGU RNA internal loop and with NMR reveals an unusual minor conformation

The RNA "GAGU" duplex, (5'GACGAGUGUCA)2, contains the internal loop (5'-GAGU-3')2 , which has two conformations in solution as determined by NMR spectroscopy. The major conformation has a loop structure consisting of trans-Watson-Crick/Hoogsteen GG pairs, A residues stacked on each other, U residues bulged outside the helix, and all sugars with a C2'-endo conformation. This differs markedly from the internal loops, (5'-GAGC-3')2, (5'-AAGU-3')2, and (5'-UAGG-3')2, which all have cis-Watson-Crick/Watson-Crick AG "imino" pairs flanked by cis-Watson-Crick/Watson-Crick canonical pairs resulting in maximal hydrogen bonding. Here, molecular dynamics was used to test whether the Amber force field (ff99 + bsc0 + OL3) approximates molecular interactions well enough to keep stable the unexpected conformation of the GAGU major duplex structure and the NMR structures of the duplexes containing (5'-GAGC-3')2, (5'-AAGU-3')2, and (5'-UAGG-3')2 internal loops. One-microsecond simulations were repeated four times for each of the duplexes starting in their NMR conformations. With the exception of (5'-UAGG-3')2, equivalent simulations were also run starting with alternative conformations. Results indicate that the Amber force field keeps the NMR conformations of the duplexes stable for at least 1 µsec. They also demonstrate an unexpected minor conformation for the (5'-GAGU-3')2 loop that is consistent with newly measured NMR spectra of duplexes with natural and modified nucleotides. Thus, unrestrained simulations led to the determination of the previously unknown minor conformation. The stability of the native (5'-GAGU-3')2 internal loop as compared to other loops can be explained by changes in hydrogen bonding and stacking as the flanking bases are changed.

Surprising Sequence Effects on GU Closure of Symmetric 2 × 2 Nucleotide RNA Internal Loops

GU base pairs are important RNA structural motifs and often close loops. Accurate prediction of RNA structures relies upon understanding the interactions determining structure. The thermodynamics of some 2 × 2 nucleotide internal loops closed by GU pairs are not well understood. Here, several self-complementary oligonucleotide sequences expected to form duplexes with 2 × 2 nucleotide internal loops closed by GU pairs were investigated. Surprisingly, nuclear magnetic resonance revealed that many of the sequences exist in equilibrium between hairpin and duplex conformations. This equilibrium is not observed with loops closed by Watson-Crick pairs. To measure the thermodynamics of some 2 × 2 nucleotide internal loops closed by GU pairs, non-self-complementary sequences that preclude formation of hairpins were designed. The measured thermodynamics indicate that some internal loops closed by GU pairs are unusually unstable. This instability accounts for the observed equilibria between duplex and hairpin conformations. Moreover, it suggests that future three-dimensional structures of loops closed by GU pairs may reveal interactions that unexpectedly destabilize folding.

Advancing viral RNA structure prediction: measuring the thermodynamics of pyrimidine-rich internal loops

Accurate thermodynamic parameters improve RNA structure predictions and thus accelerate understanding of RNA function and the identification of RNA drug binding sites. Many viral RNA structures, such as internal ribosome entry sites, have internal loops and bulges that are potential drug target sites. Current models used to predict internal loops are biased toward small, symmetric purine loops, and thus poorly predict asymmetric, pyrimidine-rich loops with >6 nucleotides (nt) that occur frequently in viral RNA. This article presents new thermodynamic data for 40 pyrimidine loops, many of which can form UU or protonated CC base pairs. Uracil and protonated cytosine base pairs stabilize asymmetric internal loops. Accurate prediction rules are presented that account for all thermodynamic measurements of RNA asymmetric internal loops. New loop initiation terms for loops with >6 nt are presented that do not follow previous assumptions that increasing asymmetry destabilizes loops. Since the last 2004 update, 126 new loops with asymmetry or sizes greater than 2 × 2 have been measured. These new measurements significantly deepen and diversify the thermodynamic database for RNA. These results will help better predict internal loops that are larger, pyrimidine-rich, and occur within viral structures such as internal ribosome entry sites.

Improved Model for Predicting the Free Energy Contribution of Dinucleotide Bulges to RNA Duplex Stability

Predicting the secondary structure of RNA is an intermediate in predicting RNA three-dimensional structure. Commonly, determining RNA secondary structure from sequence uses free energy minimization and nearest neighbor parameters. Current algorithms utilize a sequence-independent model to predict free energy contributions of dinucleotide bulges. To determine if a sequence-dependent model would be more accurate, short RNA duplexes containing dinucleotide bulges with different sequences and nearest neighbor combinations were optically melted to derive thermodynamic parameters. These data suggested energy contributions of dinucleotide bulges were sequence-dependent, and a sequence-dependent model was derived. This model assigns free energy penalties based on the identity of nucleotides in the bulge (3.06 kcal/mol for two purines, 2.93 kcal/mol for two pyrimidines, 2.71 kcal/mol for 5'-purine-pyrimidine-3', and 2.41 kcal/mol for 5'-pyrimidine-purine-3'). The predictive model also includes a 0.45 kcal/mol penalty for an A-U pair adjacent to the bulge and a -0.28 kcal/mol bonus for a G-U pair adjacent to the bulge. The new sequence-dependent model results in predicted values within, on average, 0.17 kcal/mol of experimental values, a significant improvement over the sequence-independent model. This model and new experimental values can be incorporated into algorithms that predict RNA stability and secondary structure from sequence.

A U⋅U Pair-to-U⋅C Pair Mutation-Induced RNA Native Structure Destabilisation and Stretching-Force-Induced RNA Misfolding

Little is known about how a non-Watson-Crick pair affects the RNA folding dynamics. We studied the effects of a U⋅U-to-U⋅C pair mutation on the folding of a hairpin in human telomerase RNA. The ensemble thermal melting of the hairpins shows an on-pathway intermediate with the disruption of the internal loop structure containing the U⋅U/U⋅C pairs. By using optical tweezers, we applied a stretching force on the terminal ends of the hairpins to probe directly the non-nearest-neighbour effects upon the mutations. The single U⋅U to U⋅C mutations are observed to 1) lower the mechanical unfolding force by approximately 1 picoNewton (pN) per mutation without affecting the unfolding reaction transition-state position (thus suggesting that removing a single hydrogen bond affects the structural dynamics at least two base pairs away), 2) result in more frequent misfolding into a small hairpin at approximately 10 pN and 3) shift the folding reaction transition-state position towards the native hairpin structure and slightly increase the mechanical folding kinetics (thus suggesting that untrapping from the misfolded state is not the rate-limiting step).

Geranyl diphosphate synthase: an important regulation point in balancing a recombinant monoterpene pathway in Escherichia coli

The expression level of geranyl diphosphate synthase (GPPS) was suspected to play a key role for geraniol production in recombinant Escherichia coli harboring an entire mevalonate pathway operon and a geraniol synthesis operon. The expression of GPPS was optimized by using ribosomal binding sites (RBSs) designed to have different translation initiation rates (TIRs). The RBS strength in TIR window of 500 arbitrary unit (au)-1400 au for GPPS appears to be suitable for balancing the geraniol biosynthesis pathway in this study. With the TIR of 500 au, the highest production titer of geraniol was obtained at a level of 1119mg/L, which represented a 6-fold increase in comparison with the previous titer of 183mg/L. The TIRs of GPPS locating out of range of the optimal window (500-1400 au) caused significant decreases of cell growth and geraniol production. It was suspected to result from metabolic imbalance and plasmid instability in geraniol production by inappropriate expression level of GPP synthase. Our results collectively indicated GPPS as an important regulation point in balancing a recombinant geraniol synthesis pathway. The GPPS-based regulation approach could be applicable for optimizing microbial production of other monoterpenes.

Improved model to predict the free energy contribution of trinucleotide bulges to RNA duplex stability

Trinucleotide bulges in RNA commonly occur in nature. Yet, little data exists concerning the thermodynamic parameters of this motif. Algorithms that predict RNA secondary structure from sequence currently attribute a constant free energy value of 3.2 kcal/mol to all trinucleotide bulges, regardless of bulge sequence. To test the accuracy of this model, RNA duplexes that contain frequent naturally occurring trinucleotide bulges were optically melted, and their thermodynamic parameters—enthalpy, entropy, free energy, and melting temperature—were determined. The thermodynamic data were used to derive a new model to predict the free energy contribution of trinucleotide bulges to RNA duplex stability: ΔG°_{37, trint bulge} = ΔG°_37, bulge + ΔG°_37, AU + ΔG°_37, GU. The parameter ΔG°_37, bulge is variable depending upon the purine and pyrimidine composition of the bulge, ΔG°_37, AU is a 0.49 kcal/mol penalty for an A-U closing pair, and ΔG° _37, GU is a −0.56 kcal/mol bonus for a G-U closing pair. With both closing pair and bulge sequence taken into account, this new model predicts free energy values within 0.30 kcal/mol of the experimental value. The new model can be used by algorithms that predict RNA free energies as well as algorithms that use free energy minimization to predict RNA secondary structure from sequence.

Ensemble analysis of primary microRNA structure reveals an extensive capacity to deform near the Drosha cleavage site

Most noncoding RNAs function properly only when folded into complex three-dimensional (3D) structures, but the experimental determination of these structures remains challenging. Understanding of primary microRNA (miRNA) maturation is currently limited by a lack of determined structures for nonprocessed forms of the RNA. SHAPE chemistry efficiently determines RNA secondary structural information with single-nucleotide resolution, providing constraints suitable for input into MC-Pipeline for refinement of 3D structure models. Here we combine these approaches to analyze three structurally diverse primary microRNAs, revealing deviations from canonical double-stranded RNA structure in the stem adjacent to the Drosha cut site for all three. The necessity of these deformable sites for efficient processing is demonstrated through Drosha processing assays. The structure models generated herein support the hypothesis that deformable sequences spaced roughly once per turn of A-form helix, created by noncanonical structure elements, combine with the necessary single-stranded RNA-double-stranded RNA junction to define the correct Drosha cleavage site.

The ribosome binding site calculator

The Ribosome Binding Site (RBS) Calculator is a design method for predicting and controlling translation initiation and protein expression in bacteria. The method can predict the rate of translation initiation for every start codon in an mRNA transcript. The method may also optimize a synthetic RBS sequence to achieve a targeted translation initiation rate. Using the RBS Calculator, a protein coding sequence's translation rate may be rationally controlled across a 100,000+ fold range. We begin by providing an overview of the potential biotechnology applications of the RBS Calculator, including the optimization of synthetic metabolic pathways and genetic circuits. We then detail the definitions, methodologies, and algorithms behind the RBS Calculator's thermodynamic model and optimization method. Finally, we outline a protocol for precisely measuring steady-state fluorescent protein expression levels. These methods and protocols provide a clear explanation of the RBS Calculator and its uses.

Triplex-forming MicroRNAs form stable complexes with HIV-1 provirus and inhibit its replication

BACKGROUND

One of the most fascinating discoveries in biology in recent years is unquestionably the identification of the family of small, noncoding RNAs known as microRNAs (miRNAs). Each miRNA targets multiple mRNA species through recognition of complementary sequences, typically located at multiple sites within the 3 untranslated region. In animals, single-stranded miRNA binds specific messenger RNA (mRNA) by a mechanism that is yet to be fully characterized. The bound mRNA remains untranslated resulting in reduced levels of the corresponding protein; however, if the sequence match between the miRNA and its target is precise, the bound mRNA can be degraded resulting in reduced levels of the corresponding transcript. Eukaryotic genes are also regulated by triplex formation between double helix and a third small RNA or DNA molecule. Thousands of triplex-forming (TF) islands in human genomes are mapped. However, the role of TF miRNAs within the hairpin structures of miRNA and the target mRNA has not been reported. We have explored TF complexes between human miRNAs (hsa-miR) that are complementary to human immunodeficiency virus (HIV)-1 and their antiviral potential as therapeutic agents.

METHODS

We downloaded mature miRNA sequences from the human miRBase Sequence Database (http://microrna.sanger.ac.uk/sequences/), and computationally analyzed miRNAs that have significant homologies to HIV-1 genome (pNL 4-3 Accession #AF324493). We developed an algorithm to look for triplex-binding motifs (C+CG and T AT) and selected 4 miRNAs with 3 negative controls. TF stability was tested by using fluorophore-labeled duplexes connected by a single hexaethylene glycol moiety, representing HIV-1 proviral motifs, and black-hole quencher-1 labeled oligonucleotides, representing miRNA.

RESULTS

Fifty miRNAs were discovered that showed greater than 80% homology to HIV-1, of which 4 hsa-miR that exhibited an ability to form stable triplex with double stranded-HIV-1 sequences were selected. Three negative controls were used. The TF stability of the 4 hsa-miRs and the negative controls were confirmed and measured. Stably transfected Hela-CD4+ cell lines expressing each of the hsa-miR were developed. All 4 miRNAs exhibited a significant inhibition of HIV-1 as measured by HIV-1 p24 enzyme-linked immunosorbent assay (>90%; P>0.001) when compared with the 3 negative controls. By using immunohistochemical staining with triplex binding monoclonal antibodies, significant expression of TF miRNAs was detected in the cell lines, but not in the negative controls (P<0.001).

CONCLUSIONS

In this study, we demonstrated for the first time that besides the well-established post-transcriptional silencing based on mRNA degradation, miRNAs may be responsible for long-term latency of HIV-1 by TF, a different mechanism. We provide a possible molecular mechanism by which HIV-1 homologous miRNAs may impart resistance to HIV-1 and suggest a new miRNA-based therapeutic strategy for HIV-1.

Computational approaches for RNA energy parameter estimation

Methods for efficient and accurate prediction of RNA structure are increasingly valuable, given the current rapid advances in understanding the diverse functions of RNA molecules in the cell. To enhance the accuracy of secondary structure predictions, we developed and refined optimization techniques for the estimation of energy parameters. We build on two previous approaches to RNA free-energy parameter estimation: (1) the Constraint Generation (CG) method, which iteratively generates constraints that enforce known structures to have energies lower than other structures for the same molecule; and (2) the Boltzmann Likelihood (BL) method, which infers a set of RNA free-energy parameters that maximize the conditional likelihood of a set of reference RNA structures. Here, we extend these approaches in two main ways: We propose (1) a max-margin extension of CG, and (2) a novel linear Gaussian Bayesian network that models feature relationships, which effectively makes use of sparse data by sharing statistical strength between parameters. We obtain significant improvements in the accuracy of RNA minimum free-energy pseudoknot-free secondary structure prediction when measured on a comprehensive set of 2518 RNA molecules with reference structures. Our parameters can be used in conjunction with software that predicts RNA secondary structures, RNA hybridization, or ensembles of structures. Our data, software, results, and parameter sets in various formats are freely available at http://www.cs.ubc.ca/labs/beta/Projects/RNA-Params.

Positional and neighboring base pair effects on the thermodynamic stability of RNA single mismatches

Many naturally occurring RNA structures contain single mismatches, many of which occur near the ends of helices. However, previous thermodynamic studies have focused their efforts on thermodynamically characterizing centrally placed single mismatches. Additionally, algorithms currently used to predict secondary structure from sequence are based on two assumptions for predicting the stability of RNA duplexes containing this motif. It has been assumed that the thermodynamic contribution of small RNA motifs is independent of both its position in the duplex and the identity of the non-nearest neighbors. Thermodynamically characterizing single mismatches three nucleotides from both the 3' and 5' ends (i.e., off-center) of an RNA duplex and comparing these results to those of the same single mismatch-nearest neighbor combination centrally located have allowed for the investigation of these effects. The thermodynamic contributions of 13 single mismatch-nearest neighbor combinations are reported, but only nine combinations are studied at all three duplex positions and are used to determine trends and patterns. In general, the 5'- and 3'-shifted single mismatches are relatively similar, on average, and more favorable in free energy than centrally placed single mismatches. However, close examination and comparison shows there are several associated idiosyncrasies with these identified general trends. These peculiarities may be due, in part, to the identities of the single mismatch, the nearest neighbors, and the non-nearest neighbors, along with the effects of the single mismatch position in the duplex. The prediction algorithm recently proposed by Davis and Znosko [Davis, A. R., and Znosko, B. M. (2008) Biochemistry 47, 10178-10187] is used to predict the thermodynamic parameters of single mismatch contribution, and those values are compared to the measured values presented here. This comparison suggests the proposed model is a good approximation but could be improved by the addition of parameters that account for positional and/or non-nearest neighbor effects. However, more data are required to improve our understanding of these effects and to accurately account for them.

RNA internal loops with tandem AG pairs: the structure of the 5'GAGU/3'UGAG loop can be dramatically different from others, including 5'AAGU/3'UGAA

Thermodynamic stabilities of 2 × 2 nucleotide tandem AG internal loops in RNA range from −1.3 to +3.4 kcal/mol at 37 °C and are not predicted well with a hydrogen-bonding model. To provide structural information to facilitate development of more sophisticated models for the sequence dependence of stability, we report the NMR solution structures of five RNA duplexes: (rGACGAGCGUCA)₂, (rGACUAGAGUCA)₂, (rGACAAGUGUCA)₂, (rGGUAGGCCA)₂, and (rGACGAGUGUCA)₂. The structures of these duplexes are compared to that of the previously solved (rGGCAGGCC)₂ (Wu, M., SantaLucia, J., Jr., and Turner, D. H. (1997) Biochemistry 36, 4449−4460). For loops bounded by Watson−Crick pairs, the AG and Watson−Crick pairs are all head-to-head imino-paired (cis Watson−Crick/Watson−Crick). The structures suggest that the sequence-dependent stability may reflect non-hydrogen-bonding interactions. Of the two loops bounded by G-U pairs, only the 5′UAGG/3′GGAU loop adopts canonical UG wobble pairing (cis Watson−Crick/Watson−Crick), with AG pairs that are only weakly imino-paired. Strikingly, the 5′GAGU/3′UGAG loop has two distinct duplex conformations, the major of which has both guanosine residues (G4 and G6 in (rGACGAGUGUCA)₂) in a syn glycosidic bond conformation and forming a sheared GG pair (G4-G6*, GG trans Watson−Crick/Hoogsteen), both uracils (U7 and U7*) flipped out of the helix, and an AA pair (A5-A5*) in a dynamic or stacked conformation. These structures provide benchmarks for computational investigations into interactions responsible for the unexpected differences in loop free energies and structure.

Thermodynamic characterization of naturally occurring RNA tetraloops

Although tetraloops are one of the most frequently occurring secondary structure motifs in RNA, less than one-third of the 30 most frequently occurring RNA tetraloops have been thermodynamically characterized. Therefore, 24 stem-loop sequences containing common tetraloops were optically melted, and the thermodynamic parameters DeltaH degrees , DeltaS degrees , DeltaG degrees (37,) and T(M) for each stem-loop were determined. These new experimental values, on average, are 0.7 kcal/mol different from the values predicted for these tetraloops using the model proposed by Vecenie CJ, Morrow CV, Zyra A, Serra MJ. 2006. Biochemistry 45: 1400-1407. The data for the 24 tetraloops reported here were then combined with the data for 28 tetraloops that were published previously. A new model, independent of terminal mismatch data, was derived to predict the free energy contribution of previously unmeasured tetraloops. The average absolute difference between the measured values and the values predicted using this proposed model is 0.4 kcal/mol. This new experimental data and updated predictive model allow for more accurate calculations of the free energy of RNA stem-loops containing tetraloops and, furthermore, should allow for improved prediction of secondary structure from sequence. It was also shown that tetraloops within the sequence 5'-GCCNNNNGGC-3' are, on average, 0.6 kcal/mol more stable than the same tetraloop within the sequence 5'-GGCNNNNGCC-3'. More systemic studies are required to determine the full extent of non-nearest-neighbor effects on tetraloop stability.

Automated design of synthetic ribosome binding sites to control protein expression

Microbial engineering often requires fine control over protein expression; for example, to connect genetic circuits 1-7 or control flux through a metabolic pathway 8-13. We have developed a predictive design method for synthetic ribosome binding sites that enables the rational control of a protein's production rate on a proportional scale. Experimental validation of over 100 predictions in Escherichia coli shows that the method is accurate to within a factor of 2.3 over a range of 100,000-fold. The design method also correctly predicts that reusing a ribosome binding site sequence in different genetic contexts can result in different protein expression levels. We demonstrate the method's utility by rationally optimizing a protein's expression level to connect a genetic sensor to a synthetic circuit. The proposed forward engineering approach will accelerate the construction and systematic optimization of large genetic systems.

Thermodynamic characterization of tandem mismatches found in naturally occurring RNA

Although all sequence symmetric tandem mismatches and some sequence asymmetric tandem mismatches have been thermodynamically characterized and a model has been proposed to predict the stability of previously unmeasured sequence asymmetric tandem mismatches [Christiansen,M.E. and Znosko,B.M. (2008) Biochemistry, 47, 4329–4336], experimental thermodynamic data for frequently occurring tandem mismatches is lacking. Since experimental data is preferred over a predictive model, the thermodynamic parameters for 25 frequently occurring tandem mismatches were determined. These new experimental values, on average, are 1.0 kcal/mol different from the values predicted for these mismatches using the previous model. The data for the sequence asymmetric tandem mismatches reported here were then combined with the data for 72 sequence asymmetric tandem mismatches that were published previously, and the parameters used to predict the thermodynamics of previously unmeasured sequence asymmetric tandem mismatches were updated. The average absolute difference between the measured values and the values predicted using these updated parameters is 0.5 kcal/mol. This updated model improves the prediction for tandem mismatches that were predicted rather poorly by the previous model. This new experimental data and updated predictive model allow for more accurate calculations of the free energy of RNA duplexes containing tandem mismatches, and, furthermore, should allow for improved prediction of secondary structure from sequence.

Thermodynamic characterization of naturally occurring RNA single mismatches with G-U nearest neighbors

Due to their prevalence and roles in biological systems, single mismatches adjacent to G-U pairs are important RNA structural elements. Since there are only limited experimental values for the stability of single mismatches adjacent to G-U pairs, current algorithms using free energy minimization to predict RNA secondary structure from sequence assign predicted thermodynamic values to these types of single mismatches. Here, thermodynamic data are reported for frequently occurring single mismatches adjacent to at least one G-U pair. This experimental data can be used in place of predicted thermodynamic values in algorithms that predict secondary structure from sequence using free energy minimization. When predicting the thermodynamic contributions of previously unmeasured single mismatches, most algorithms apply the same thermodynamic penalty for an A-U pair adjacent to a single mismatch and a G-U pair adjacent to a single mismatch. A recent study, however, suggests that the penalty for a G-U pair adjacent to a tandem mismatch should be 1.2 +/- 0.1 kcal/mol, and the penalty for an A-U pair adjacent to a tandem mismatch should be 0.5 +/- 0.2 kcal/mol [Christiansen, M. E. and Znosko, B. M. (2008) Biochemistry 47, 4329-4336]. Therefore, the data reported here are combined with the existing thermodynamic dataset of single mismatches, and nearest neighbor parameters are derived for an A-U pair adjacent to a single mismatch (1.1 +/- 0.1 kcal/mol) and a G-U pair adjacent to a single mismatch (1.4 +/- 0.1 kcal/mol).