Atomistic investigation of the effect of incremental modification of deoxyribose sugars by locked nucleic acid (β-D-LNA and α-L-LNA) moieties on the structures and thermodynamics of DNA-RNA hybrid duplexes

SN- and NS-puckered sugar conformers are precursors of the (6-4) photoproduct in thymine dinucleotide

Some amount of furanose in a southern conformation, possibly in both, but certainly in one of the two adjacent nucleotides of a dipyrimidine site, is necessary for (6-4) photoproduct formation in oligonucleotides. To explore the necessity, role, and most favorable location of each South sugar conformer in the formation of the (6-4) adduct in the thymine dinucleotide TpT, the photochemical behavior of two synthetic analogues, in which the South sugar conformation is prohibited for one of their two sugars, has been examined. Herein, we experimentally demonstrate that the presence of one sugar presenting some amount of South puckering, at any of the extremities, is sufficient to trigger (6-4) adduct formation. Nonetheless, the photochemical behavior of the dinucleotide with a South-puckered conformation at the 5'-end, mimics more closely that of TpT. In addition, using the 5' North 3' South-dilocked dinucleotide, we demonstrate that the flexibility of the South pucker at the 3'-end has little influence on the (6-4) adduct formation.

Molecular Dynamics Study of the Hybridization between RNA and Modified Oligonucleotides

MicroRNAs (miRNAs) are attractive drug candidates for many diseases as they can modulate the expression of gene networks. Recently, we discovered that DNAs targeting microRNA-22-3p (miR-22-3p) hold the potential for treating obesity and related metabolic disorders (type 2 diabetes mellitus, hyperlipidemia, and nonalcoholic fatty liver disease (NAFLD)) by turning fat-storing white adipocytes into fat-burning adipocytes. In this work, we explored the effects of chemical modifications, including phosphorothioate (PS), locked nucleic acid (LNA), and peptide nucleic acid (PNA), on the structure and energy of DNA analogs by using molecular dynamics (MD) simulations. To achieve a reliable prediction of the hybridization free energy, the AMOEBA polarizable force field and the free energy perturbation technique were employed. The calculated hybridization free energies are generally compatible with previous experiments. For LNA and PNA, the enhanced duplex stability can be explained by the preorganization mechanism, i.e., the single strands adopt stable helical structures similar to those in the duplex. For PS, the S and R isomers (Sp and Rp) have preferences for C2'-endo and C3'-endo sugar puckering conformations, respectively, and therefore Sp is less stable than Rp in DNA/RNA hybrids. In addition, the solvation penalty of Rp accounts for its destabilization effect. PS-LNA is similar to LNA as the sugar puckering is dominated by the locked sugar ring. This work demonstrated that MD simulations with polarizable force fields are useful for the understanding and design of modified nucleic acids.

The ability of locked nucleic acid oligonucleotides to pre-structure the double helix: A molecular simulation and binding study

Locked nucleic acid (LNA) oligonucleotides bind DNA target sequences forming Watson-Crick and Hoogsteen base pairs, and are therefore of interest for medical applications. To be biologically active, such an oligonucleotide has to efficiently bind the target sequence. Here we used molecular dynamics simulations and electrophoresis mobility shift assays to elucidate the relation between helical structure and affinity for LNA-containing oligonucleotides. In particular, we have studied how LNA substitutions in the polypyrimidine strand of a duplex (thus forming a hetero duplex, i.e. a duplex with a DNA polypurine strand and an LNA/DNA polypyrimidine strand) enhance triplex formation. Based on seven polypyrimidine single strand oligonucleotides, having LNAs in different positions and quantities, we show that alternating LNA with one or more non-modified DNA nucleotides pre-organizes the hetero duplex toward a triple-helical-like conformation. This in turn promotes triplex formation, while consecutive LNAs distort the duplex structure disfavoring triplex formation. The results support the hypothesis that a pre-organization in the hetero duplex structure enhances the binding of triplex forming oligonucleotides. Our findings may serve as a criterion in the design of new tools for efficient oligonucleotide hybridization.

Modeling and Simulation of Oligonucleotide Hybrids: Outlining a Strategy

Molecular dynamics simulations with a state-of-the-art force field provide an atomistic detailed description of the structural and thermodynamic features of biomolecules. Effects of chemical modifications and of the environment such as sequence, solvent, and ionic strength can explicitly be taken into account. Molecular simulation techniques can also provide insight in change in binding affinity, in protonation (pK_a shift) and tautomeric propensity due to changes in the environment or in the molecular system. The quality and reliability of a simulation depend on the quality of the force field and on the reproducibility of the data, and validation depends on the availability of suitable experimental data. Here, we describe the workflow to investigate oligonucleotide hybrids using molecular simulation including hardware and software information.

Systematics for types and effects of DNA variations

BACKGROUND

Numerous different types of variations can occur in DNA and have diverse effects and consequences. The Variation Ontology (VariO) was developed for systematic descriptions of variations and their effects at DNA, RNA and protein levels.

RESULTS

VariO use and terms for DNA variations are described in depth. VariO provides systematic names for variation types and detailed descriptions for changes in DNA function, structure and properties. The principles of VariO are presented along with examples from published articles or databases, most often in relation to human diseases. VariO terms describe local DNA changes, chromosome number and structure variants, chromatin alterations, as well as genomic changes, whether of genetic or non-genetic origin.

CONCLUSIONS

DNA variation systematics facilitates unambiguous descriptions of variations and their effects and further reuse and integration of data from different sources by both human and computers.

Modeling p K Shift in DNA Triplexes Containing Locked Nucleic Acids

The protonation states for nucleic acid bases are difficult to assess experimentally. In the context of DNA triplex, the protonation state of cytidine in the third strand is particularly important, because it needs to be protonated in order to form Hoogsteen hydrogen bonds. A sugar modification, locked nucleic acid (LNA), is widely used in triplex forming oligonucleotides to target sites in the human genome. In this study, the parameters for LNA are developed in line with the CHARMM nucleic acid force field and validated toward the available structural experimental data. In conjunction, two computational methods were used to calculate the protonation state of the third strand cytidine in various DNA triplex environments: λ-dynamics and multiple pH regime. Both approaches predict p K of this cytidine shifted above physiological pH when cytidine is in the third strand in a triplex environment. Both methods show an upshift due to cytidine methylation, and a small downshift when the sugar configuration is locked. The predicted p K values for cytidine in DNA triplex environment can inform the design of better-binding oligonucleotides.

The free energy of locking a ring: Changing a deoxyribonucleoside to a locked nucleic acid

Locked nucleic acid (LNA), a modified nucleoside which contains a bridging group across the ribose ring, improves the stability of DNA/RNA duplexes significantly, and therefore is of interest in biotechnology and gene therapy applications. In this study, we investigate the free energy change between LNA and DNA nucleosides. The transformation requires the breaking of the bridging group across the ribose ring, a problematic transformation in free energy calculations. To address this, we have developed a 3-step (easy to implement) and a 1-step protocol (more efficient, but more complicated to setup), for single and dual topologies in classical molecular dynamics simulations, using the Bennett Acceptance Ratio method to calculate the free energy. We validate the approach on the solvation free energy difference for the nucleosides thymidine, cytosine, and 5-methyl-cytosine. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Ligand-Induced Stabilization of a Duplex-like Architecture Is Crucial for the Switching Mechanism of the SAM-III Riboswitch

Riboswitches are structured RNA motifs that control gene expression by sensing the concentrations of specific metabolites and make up a promising new class of antibiotic targets. S-Adenosylmethionine (SAM)-III riboswitch, mainly found in lactic acid bacteria, is involved in regulating methionine and SAM biosynthetic pathways. SAM-III riboswitch regulates the gene expression by switching the translation process on and off with respect to the absence and presence of the SAM ligand, respectively. In this study, an attempt is made to understand the key conformational transitions involved in ligand binding using atomistic molecular dynamics (MD) simulations performed in an explicit solvent environment. G26 is found to recognize the SAM ligand by forming hydrogen bonds, whereas the absence of the ligand leads to opening of the binding pocket. Consistent with experimental results, the absence of the SAM ligand weakens the base pairing interactions between the nucleobases that are part of the Shine-Dalgarno (SD) and anti-Shine-Dalgarno (aSD) sequences, which in turn facilitates recognition of the SD sequence by ribosomes. Detailed analysis reveals that a duplex-like structure formed by nucleotides from different parts of the RNA and the adenine base of the ligand is crucial for the stability of the completely folded state in the presence of the ligand. Previous experimental studies have shown that the SAM-III riboswitch exists in equilibrium between the unfolded and partially folded states in the absence of the ligand, which completely folds upon binding of the ligand. Comparison of the results presented here to the available experimental data indicates the structures obtained using the MD simulations resemble the partially folded state. Thus, this study provides a detailed understanding of the fully and partially folded structures of the SAM-III riboswitch in the presence and absence of the ligand, respectively. This study hypothesizes a dual role for the SAM ligand, which facilitates conformational switching between partially and fully folded states by forming a stable duplex-like structure and strengthening the interactions between SD and aSD nucleotides.

Therapeutic Oligonucleotides Targeting Liver Disease: TTR Amyloidosis

The liver has become an increasingly interesting target for oligonucleotide therapy. Mutations of the gene encoding transthyretin (TTR), expressed in vast amounts by the liver, result in a complex degenerative disease, termed familial amyloid polyneuropathy (FAP). Misfolded variants of TTR are linked to the establishment of extracellular protein deposition in various tissues, including the heart and the peripheral nervous system. Recent progress in the chemistry and formulation of antisense (ASO) and small interfering RNA (siRNA) designed for a knockdown of TTR mRNA in the liver has allowed to address the issue of gene-specific molecular therapy in a clinical setting of FAP. The two therapeutic oligonucleotides bind to RNA in a sequence specific manner but exploit different mechanisms. Here we describe major developments that have led to the advent of therapeutic oligonucleotides for treatment of TTR-related disease.

DNA–RNA hybrid duplexes with decreasing pyrimidine content in the DNA strand provide structural snapshots for the A- to B-form conformational transition of nucleic acids

A gradual increase in the deoxypyrimidine content in DNA–RNA hybrids leads to B- to A-form nucleic acid transition. Possible factors that govern nuclease activity on hybrid duplexes are presented.