Structures, dynamics, and stabilities of fully modified locked nucleic acid (β-D-LNA and α-L-LNA) duplexes in comparison to pure DNA and RNA duplexes

Inhibition of Chikungunya virus genome replication by targeting essential RNA structures within the virus genome

Chikungunya virus (CHIKV) is a pathogenic arbovirus spread by Aedes spp. mosquitos. CHIKV has a wide global prevalence and represents a significant health burden in affected populations. Symptoms of CHIKV infection include fever, rashes and debilitating joint and muscle pain, which can persist for several months to years in some patients. To date, there remains no vaccine or specific antiviral therapy against this important human pathogen. Based on our previously published structural and phenotypic analysis of the 5' region of the CHIKV genome, we designed a panel of locked nucleic acid oligonucleotides to bind structured RNA replication elements within the virus genome, which are essential for efficient CHIKV replication. Using electromobility shift assays, we confirmed the relative binding efficiencies of each LNA to target CHIKV genomic RNA. We then went on to demonstrate, using both sub-genomic replicon and infectious virus systems, that targeting individual RNA replication elements inhibits CHIKV genome replication and production of infectious virus. Time course assays demonstrated that LNAs can access the CHIKV replication complex and virus genome, during active virus replication. For the first time, these findings show that functional RNA elements can be specifically targeted during the CHIKV lifecycle and consequently represent potential novel antiviral targets.

Locked nucleic acid (LNA): A modern approach to cancer diagnosis and treatment

Cancer is responsible for about one in six deaths in the world. Conventional cancer treatments like chemotherapy, radiotherapy, and surgery are associated with drug poisoning and poor prognosis. Thanks to advances in RNA delivery and target selection, new cancer medicines are now conceivable to improve the quality of life and extend the lives of cancer patients. Antisense oligonucleotides (ASOs) and siRNAs are the most important tools in RNA therapies. Locked Nucleic Acids (LNAs) are one of the newest RNA analogs, exhibiting more affinity to binding, sequence specificity, thermal stability, and nuclease resistance due to their unique properties. Assays using LNA are also used in molecular diagnostic methods and provide accurate and rapid mutation detection that improves specificity and sensitivity. This study aims to review the special properties of LNA oligonucleotides that make them safe and effective antisense drugs for cancer treatment by controlling gene expression. Following that, we go over all of the molecular detection methods and cancer treatment antisense tactics that are possible with LNA technology.

Probe-Based Real-Time qPCR Assays for a Reliable Differentiation of Capripox Virus Species

Outbreaks of the three capripox virus species, namely lumpy skin disease virus, sheeppox virus, and goatpox virus, severely affect animal health and both national and international economies. Therefore, the World Organization for Animal Health (OIE) classified them as notifiable diseases. Until now, discrimination of capripox virus species was possible by using different conventional PCR protocols. However, more sophisticated probe-based real-time qPCR systems addressing this issue are, to our knowledge, still missing. In the present study, we developed several duplex qPCR assays consisting of different types of fluorescence-labelled probes that are highly sensitive and show a high analytical specificity. Finally, our assays were combined with already published diagnostic methods to a diagnostic workflow that enables time-saving, reliable, and robust detection, differentiation, and characterization of capripox virus isolates.

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.

Structural insight into antisense gapmer-RNA oligomer duplexes through molecular dynamics simulations

There is an extensive research carrying out on antisense technology and the molecules entering into clinical trials are increasing rapidly. Phosphorothioate (PS) is a chemical modification in which nonbridged oxygen is replaced with a sulfur, consequently providing resistance against nuclease activity. The 2'-4' conformationally restricted nucleoside has the structural features of both 2'-O-methoxy ethyl RNA (MOE), which shows good toxicity profile, and locked nucleic acid (LNA), which shows good binding affinity towards the target RNA. These modifications have been studied and suggested that they can be a potential therapeutic agents in antisense therapy. Mipomersen (ISIS 301012), which contains the novel nucleoside modification has been used to target to apolipoprotein (Apo B), which reduces LDL cholesterol by 6-41%. In this study, classical molecular dynamics (MD) simulations were performed on six different antisense gapmer/target-RNA oligomer duplexes (LNA-PS-LNA/RNA, RcMOE-PS-RcMOE/RNA, ScMOE-PS-ScMOE/RNA, MOE-PS-MOE/RNA, PS-DNA/RNA and DNA/RNA) to investigate the structural dynamics, stability and solvation properties. The LNA, MOE nucleotides present in respective duplexes are showing the structure of A-form and the PS-DNA nucleotides resemble the structure of B-form helix with respect to some of the helical parameters. Free energy calculations suggest that the oligomer, which contains LNA binds to the RNA strongly than other modifications as shown in experimental results. The MOE modified nucleotide, which although had a lower binding affinity but higher solvent accessible surface area (SASA) compared to the other modifications, may be influencing the toxicity and hence may be used it in Mipomersen, the second antisense molecule which is approved by FDA. Communicated by Ramaswamy H. Sarma.

A highly sensitive 1-tube nested real-time RT-PCR assay using LNA-modified primers for detection of respiratory syncytial virus

Respiratory syncytial virus (RSV) causes serious respiratory tract infection worldwide. The relatively low RSV load makes it difficult to detect in frail, elderly, and severely immune-compromised patients. In the present study, we developed a locked nucleic acid–-based 1-tube nested real-time RT-PCR (OTNRT-PCR) assay with the advantages of extremely high sensitivity, facile operability, and less likelihood of cross-contamination. The sensitivity, specificity, and clinical performance of the OTNRT-PCR assay were compared in parallel with a conventional TaqMan probe-based real-time PCR (qRT-PCR) assay and a traditional 2-step nested RT-PCR assay. The limit of detection of the OTNRT-PCR assay was 1.02 × 10⁻¹ TCID50/mL, equivalent to the traditional 2-step nested RT-PCR assay and 25-fold lower than the qRT-PCR assay. Of 616 nasopharyngeal aspirates tested, 143 RSV-negative samples by qRT-PCR were confirmed as positive by sequencing the OTNRT-PCR products. We therefore conclude that OTNRT-PCR is more sensitive than qRT-PCR for detection of RSV in clinical samples.

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.

Inclusion of methoxy groups inverts the thermodynamic stabilities of DNA-RNA hybrid duplexes: A molecular dynamics simulation study

Modified nucleic acids have found profound applications in nucleic acid based technologies such as antisense and antiviral therapies. Previous studies on chemically modified nucleic acids have suggested that modifications incorporated in furanose sugar especially at 2'-position attribute special properties to nucleic acids when compared to other modifications. 2'-O-methyl modification to deoxyribose sugars of DNA-RNA hybrids is one such modification that increases nucleic acid stability and has become an attractive class of compounds for potential antisense applications. It has been reported that modification of DNA strands with 2'-O-methyl group reverses the thermodynamic stability of DNA-RNA hybrid duplexes. Molecular dynamics simulations have been performed on two hybrid duplexes (DR and RD) which differ from each other and 2'-O-methyl modified counterparts to investigate the effect of 2'-O-methyl modification on their duplex stability. The results obtained suggest that the modification drives the conformations of both the hybrid duplexes towards A-RNA like conformation. The modified hybrid duplexes exhibit significantly contrasting dynamics and hydration patterns compared to respective parent duplexes. In line with the experimental results, the relative binding free energies suggest that the introduced modifications stabilize the less stable DR hybrid, but destabilize the more stable RD duplex. Binding free energy calculations suggest that the increased hydrophobicity is primarily responsible for the reversal of thermodynamic stability of hybrid duplexes. Free energy component analysis further provides insights into the stability of modified duplexes.

Interplay of LNA and 2'-O-methyl RNA in the structure and thermodynamics of RNA hybrid systems: a molecular dynamics study using the revised AMBER force field and comparison with experimental results

When used in nucleic acid duplexes, locked nucleic acid (LNA) and 2'-O-methyl RNA residues enhance the duplex stabilities, and this makes it possible to create much better RNA aptamers to target specific molecules in cells. Thus, LNA and 2'-O-methyl RNA residues are finding increasingly widespread use in RNA-based therapeutics. Herein, we utilize molecular dynamics (MD) simulations and UV melting experiments to investigate the structural and thermodynamic properties of 13 nucleic acid duplexes, including full DNA, RNA, LNA, and 2'-O-methyl RNA duplexes as well as hybrid systems such as LNA:RNA, 2'-O-methyl RNA:RNA, LNA/2'-O-methyl RNA:RNA, and RNA/2'-O-methyl RNA:RNA duplexes. The MD simulations are based on a version of the Amber force field revised specifically for RNA and LNA residues. Our results indicate that LNA and 2'-O-methyl RNA residues have two different hybridization mechanisms when included in hybrid duplexes with RNA wherein the former underwinds while the latter overwinds the duplexes. These computational predictions are supported by X-ray structures of LNA and 2'-O-methyl RNA duplexes that were recently presented by different groups, and there is also good agreement with the measured thermal stabilities of the duplexes. We find out that the "underwinding" phenomenon seen in LNA and LNA:RNA hybrid duplexes happens due to expansion of the major groove widths (Mgw) of the duplexes that is associated with decrease in the slide and twist values in base-pair steps. In contrast, 2'-O-methyl RNA residues in RNA duplexes slightly overwind the duplexes while the backbone is forced to stay in C3'-endo. Moreover, base-pair stacking in the LNA and LNA:RNA hybrid systems is gradually reduced with the inclusion of LNA residues in the duplexes while no such effect is seen in the 2'-O-methyl RNA systems. Our results show how competition between base stacking and structural rigidity in these RNA hybrid systems influences structures and stabilities. Even though both LNA and 2'-O-methyl RNA residues have C3'-endo sugar puckering, structurally LNA residues have a frozen sugar backbone which provides entropic enhancement of stabilities while the 2'-O-methyl RNA residues are more flexible and maintain base stacking that is almost untouched compared to RNA. Thus, enhancement of the structural stabilities of RNA duplexes by 2'-O-methyl RNA modifications is smaller than for the corresponding LNA modifications. Indeed, our experimental measurements show that on average each 2'-O-methyl RNA and LNA substitution in a RNA duplex enhances duplex stability by 0.2 and 1.4 kcal/mol, respectively. Our computational binding free energy predictions are qualitatively in line with these results. The only exception is for the full 2'-O-methyl RNA duplex, which is overstabilized, implying that further force field revisions are needed. Collectively, the results presented in this paper explain the atomistic details of the structural and thermodynamic roles of LNA and 2'-O-methyl RNA residues in RNA hybrid duplexes, shedding light on the mechanism behind targeting endogenous micro RNA (miRNA) in order to regulate mRNA activity and inhibit gene expression in the cell.

Structural and energetic characterization of the major DNA adduct formed from the food mutagen ochratoxin A in the NarI hotspot sequence: influence of adduct ionization on the conformational preferences and implications for the NER propensity

The nephrotoxic food mutagen ochratoxin A (OTA) produces DNA adducts in rat kidneys, the major lesion being the C8-linked-2′-deoxyguanosine adduct (OTB-dG). Although research on other adducts stresses the importance of understanding the structure of the associated adducted DNA, site-specific incorporation of OTB-dG into DNA has yet to be attempted. The present work uses a robust computational approach to determine the conformational preferences of OTB-dG in three ionization states at three guanine positions in the NarI recognition sequence opposite cytosine. Representative adducted DNA helices were derived from over 2160 ns of simulation and ranked via free energies. For the first time, a close energetic separation between three distinct conformations is highlighted, which indicates OTA-adducted DNA likely adopts a mixture of conformations regardless of the sequence context. Nevertheless, the preferred conformation depends on the flanking bases and ionization state due to deviations in discrete local interactions at the lesion site. The structural characteristics of the lesion thus discerned have profound implications regarding its repair propensity and mutagenic outcomes, and support recent experiments suggesting the induction of double-strand breaks and deletion mutations upon OTA exposure. This combined structural and energetic characterization of the OTB-dG lesion in DNA will encourage future biochemical experiments on this potentially genotoxic lesion.

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

Chemically modified oligonucleotides offer many possibilities in utilizing their special features for a vast number of applications in nucleic acid based therapies and synthetic molecular biology. Locked nucleic acid analogues (α-/β-LNA) are modifications having an extra ring of 2'-O,4'-C-methylene group in the furanose sugar. LNA strands have been shown to exhibit high binding affinity toward RNA and DNA strands, and the resultant duplexes show significantly high melting temperatures. In the present study, molecular dynamics (MD) simulations were performed on DNA-RNA hybrid duplexes by systematically modifying their deoxyribose sugars with locked nucleic acid analogues. Several geometrical and energetic analyses were performed using principal component (PCA) analysis and binding free energy methods to understand the consequence of incorporated isomeric LNA modifications on the structure, dynamics, and stability of DNA-RNA hybrid duplex. The β-modification systematically changes the conformation of the DNA-RNA hybrid duplex whereas drastic changes are observed for α-modification. The fully modified duplexes have distinct properties compared to partial and unmodified duplexes, and the partly modified duplexes have properties intermediate to full strand and unmodified duplexes. The distribution of BI versus BII populations suggests that backbone rearrangement is minimal for β-LNA modification in order to accommodate it in duplexes whereas extensive backbone rearrangement is necessary in order to incorporate α-LNA modification which subsequently alters the energetic and structural properties of the duplexes. The simulation results also suggest that the alteration of DNA-RNA hybrid properties depends on the position of modification and the gap between the modifications.

Recognition of 2',5'-linked oligoadenylates by human ribonuclease L: molecular dynamics study

The capability of current MD simulations to be used as a tool in rational design of agonists of medically interesting enzyme RNase L was tested. Dimerization and enzymatic activity of RNase L is stimulated by 2',5'-linked oligoadenylates (pA₂₅A₂₅A; 2-5A). First, it was necessary to ensure that a complex of monomeric human RNase L and 25A was stable in MD simulations. It turned out that Glu131 had to be protonated. The non-protonated Glu131 caused dissociation of 2-5A from RNase L. Because of the atypical 2'-5' internucleotide linkages and a specific spatial arrangement of the 25A trimer, when a single molecule carries all possible conformers of the glycosidic torsion angle, several versions of the AMBER force field were tested. One that best maintained functionally important interactions of 25A and RNase L was selected for subsequent MD simulations. Furthermore, we wonder whether powerful GPUs are able to produce MD trajectories long enough to convincingly demonstrate effects of subtle perturbations of interactions between 25A and RNase L. Detrimental impacts of various point mutations of RNase L (R155A, F126A, W60A, K89A) on 2-5A binding were observed on a time scale of 200 ns. Finally, 2-5A analogues with a bridged 3'--O,4'--C-alkylene linkage (B) introduced into the adenosine units (A) were used to assess ability of MD simulations to distinguish on the time scale of hundreds of nanoseconds between agonists of RNase L (pA₂₅A₂₅B, pB₂₅A₂₅A, pB₂₅A₂₅B) and inactive analogs (pA₂₅B₂₅A, pA₂₅B₂₅B, pB₂₅B₂₅A, pB₂₅B₂₅B). Agonists were potently bound to RNase L during 200 ns MD runs. For inactive 2-5A analogs, by contrast, significant disruptions of their interactions with RNase L already within 100 ns MD runs were found.

Optimization of an AMBER force field for the artificial nucleic acid, LNA, and benchmarking with NMR of L(CAAU)

Locked Nucleic Acids (LNAs) are RNA analogues with an O2'-C4' methylene bridge which locks the sugar into a C3'-endo conformation. This enhances hybridization to DNA and RNA, making LNAs useful in microarrays and potential therapeutics. Here, the LNA, L(CAAU), provides a simplified benchmark for testing the ability of molecular dynamics (MD) to approximate nucleic acid properties. LNA χ torsions and partial charges were parametrized to create AMBER parm99_LNA. The revisions were tested by comparing MD predictions with AMBER parm99 and parm99_LNA against a 200 ms NOESY NMR spectrum of L(CAAU). NMR indicates an A-Form equilibrium ensemble. In 3000 ns simulations starting with an A-form structure, parm99_LNA and parm99 provide 66% and 35% agreement, respectively, with NMR NOE volumes and (3)J-couplings. In simulations of L(CAAU) starting with all χ torsions in a syn conformation, only parm99_LNA is able to repair the structure. This implies methods for parametrizing force fields for nucleic acid mimics can reasonably approximate key interactions and that parm99_LNA will improve reliability of MD studies for systems with LNA. A method for approximating χ population distribution on the basis of base to sugar NOEs is also introduced.

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.