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

Molecular insight into how the position of an abasic site modifies DNA duplex stability and dynamics

Local perturbations to DNA base-pairing stability from lesions and chemical modifications can alter the stability and dynamics of an entire oligonucleotide. End effects may cause the position of a disruption within a short duplex to influence duplex stability and structural dynamics, yet this aspect of nucleic acid modifications is often overlooked. We investigate how the position of an abasic site (AP site) impacts the stability and dynamics of short DNA duplexes. Using a combination of steady-state and time-resolved spectroscopy and molecular dynamics simulations, we unravel an interplay between AP-site position and nucleobase sequence that controls energetic and dynamic disruption to the duplex. The duplex is disrupted into two segments by an entropic barrier for base-pairing on each side of the AP site. The barrier induces fraying of the short segment when an AP site is near the termini. Shifting the AP site inward promotes a transition from short-segment fraying to fully encompassing the barrier into the thermodynamics of hybridization, leading to further destabilization of the duplex. Nucleobase sequence determines the length scale for this transition by tuning the barrier height and base-pair stability of the short segment, and certain sequences enable out-of-register base-pairing to minimize the barrier height.

Tailor made: the art of therapeutic mRNA design

mRNA medicine is a new and rapidly developing field in which the delivery of genetic information in the form of mRNA is used to direct therapeutic protein production in humans. This approach, which allows for the quick and efficient identification and optimization of drug candidates for both large populations and individual patients, has the potential to revolutionize the way we prevent and treat disease. A key feature of mRNA medicines is their high degree of designability, although the design choices involved are complex. Maximizing the production of therapeutic proteins from mRNA medicines requires a thorough understanding of how nucleotide sequence, nucleotide modification and RNA structure interplay to affect translational efficiency and mRNA stability. In this Review, we describe the principles that underlie the physical stability and biological activity of mRNA and emphasize their relevance to the myriad considerations that factor into therapeutic mRNA design.

DNA Structure Design Is Improved Using an Artificially Expanded Alphabet of Base Pairs Including Loop and Mismatch Thermodynamic Parameters

We show that in silico design of DNA secondary structures is improved by extending the base pairing alphabet beyond A-T and G-C to include the pair between 2-amino-8-(1'-β-d-2'-deoxyribofuranosyl)-imidazo-[1,2-a]-1,3,5-triazin-(8H)-4-one and 6-amino-3-(1'-β-d-2'-deoxyribofuranosyl)-5-nitro-(1H)-pyridin-2-one, abbreviated as P and Z. To obtain the thermodynamic parameters needed to include P-Z pairs in the designs, we performed 47 optical melting experiments and combined the results with previous work to fit free energy and enthalpy nearest neighbor folding parameters for P-Z pairs and G-Z wobble pairs. We find G-Z pairs have stability comparable to that of A-T pairs and should therefore be included as base pairs in structure prediction and design algorithms. Additionally, we extrapolated the set of loop, terminal mismatch, and dangling end parameters to include the P and Z nucleotides. These parameters were incorporated into the RNAstructure software package for secondary structure prediction and analysis. Using the RNAstructure Design program, we solved 99 of the 100 design problems posed by Eterna using the ACGT alphabet or supplementing it with P-Z pairs. Extending the alphabet reduced the propensity of sequences to fold into off-target structures, as evaluated by the normalized ensemble defect (NED). The NED values were improved relative to those from the Eterna example solutions in 91 of 99 cases in which Eterna-player solutions were provided. P-Z-containing designs had average NED values of 0.040, significantly below the 0.074 of standard-DNA-only designs, and inclusion of the P-Z pairs decreased the time needed to converge on a design. This work provides a sample pipeline for inclusion of any expanded alphabet nucleotides into prediction and design workflows.

Molecular insight into how the position of an abasic site and its sequence environment influence DNA duplex stability and dynamics

Local perturbations to DNA base-pairing stability from lesions and chemical modifications can alter the stability and dynamics of an entire oligonucleotide. End effects may cause the position of a disruption within a short duplex to influence duplex stability and structural dynamics, yet this aspect of nucleic acid modifications is often overlooked. We investigate how the position of an abasic site (AP site) impacts the stability and dynamics of short DNA duplexes. Using a combination of steady-state and time-resolved spectroscopy and molecular dynamics simulations, we unravel an interplay between AP-site position and nucleobase sequence that controls energetic and dynamic disruption to the duplex. The duplex is disrupted into two segments by an entropic barrier for base pairing on each side of the AP site. The barrier induces fraying of the short segment when an AP site is near the termini. Shifting the AP site inward promotes a transition from short-segment fraying to fully encompassing the barrier into the thermodynamics of hybridization, leading to further destabilization the duplex. Nucleobase sequence determines the length scale for this transition by tuning the barrier height and base-pair stability of the short segment, and certain sequences enable out-of-register base pairing to minimize the barrier height.

RNA Secondary Structure Analysis Using RNAstructure

RNAstructure is a user-friendly program for the prediction and analysis of RNA secondary structure. It is available as a web server, a program with a graphical user interface, or a set of command line tools. The programs are available for Microsoft Windows, macOS, or Linux. This article provides protocols for prediction of RNA secondary structure (using the web server, the graphical user interface, or the command line) and high-affinity oligonucleotide binding sites to a structured RNA target (using the graphical user interface). © 2023 Wiley Periodicals LLC. Basic Protocol 1: Predicting RNA secondary structure using the RNAstructure web server Alternate Protocol 1: Predicting secondary structure and base pair probabilities using the RNAstructure graphical user interface Alternate Protocol 2: Predicting secondary structure and base pair probabilities using the RNAstructure command line interface Basic Protocol 2: Predicting binding affinities of oligonucleotides complementary to an RNA target using OligoWalk.

DNA Structure Design Is Improved Using an Artificially Expanded Alphabet of Base Pairs Including Loop and Mismatch Thermodynamic Parameters

We show that in silico design of DNA secondary structures is improved by extending the base pairing alphabet beyond A-T and G-C to include the pair between 2-amino-8-(1'-β-D-2'-deoxyribofuranosyl)-imidazo-[1,2- a ]-1,3,5-triazin-(8 H )-4-one and 6-amino-3-(1'-β-D-2'-deoxyribofuranosyl)-5-nitro-(1 H )-pyridin-2-one, simply P and Z. To obtain the thermodynamic parameters needed to include P-Z pairs in the designs, we performed 47 optical melting experiments and combined the results with previous work to fit a new set of free energy and enthalpy nearest neighbor folding parameters for P-Z pairs and G-Z wobble pairs. We find that G-Z pairs have stability comparable to A-T pairs and therefore should be considered quantitatively by structure prediction and design algorithms. Additionally, we extrapolated the set of loop, terminal mismatch, and dangling end parameters to include P and Z nucleotides. These parameters were incorporated into the RNAstructure software package for secondary structure prediction and analysis. Using the RNAstructure Design program, we solved 99 of the 100 design problems posed by Eterna using the ACGT alphabet or supplementing with P-Z pairs. Extending the alphabet reduced the propensity of sequences to fold into off-target structures, as evaluated by the normalized ensemble defect (NED). The NED values were improved relative to those from the Eterna example solutions in 91 of 99 cases where Eterna-player solutions were provided. P-Z-containing designs had average NED values of 0.040, significantly below the 0.074 of standard-DNA-only designs, and inclusion of the P-Z pairs decreased the time needed to converge on a design. This work provides a sample pipeline for inclusion of any expanded alphabet nucleotides into prediction and design workflows.