Molecular Dynamics Simulations Highlight the Structural Differences among DNA:DNA, RNA:RNA, and DNA:RNA Hybrid Duplexes

Nanopore Translocation Reveals Electrophoretic Force on Noncanonical RNA:DNA Double Helix

Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on the translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of noncanonical RNA:DNA hybrids in electrophoretic transport to the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a noncanonical helix, with distinct transport properties from B-form DD molecules. We find that RD and DD molecules, with the same contour length, move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and the molecular understanding of electrophoretic transport.

Nanopore translocation reveals electrophoretic force on non-canonical RNA:DNA double helix

Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of non-canonical RNA:DNA hybrids in electrophoretic transport with the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a non-canonical helix with distinct transport properties from B-form DD molecules. We find RD and DD molecules with the same contour length move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same within experimental errors. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and molecular understanding of electrophoretic transport.

Coarse-grained modeling of DNA-RNA hybrids

We introduce oxNA, a new model for the simulation of DNA-RNA hybrids that is based on two previously developed coarse-grained models-oxDNA and oxRNA. The model naturally reproduces the physical properties of hybrid duplexes, including their structure, persistence length, and force-extension characteristics. By parameterizing the DNA-RNA hydrogen bonding interaction, we fit the model's thermodynamic properties to experimental data using both average-sequence and sequence-dependent parameters. To demonstrate the model's applicability, we provide three examples of its use-calculating the free energy profiles of hybrid strand displacement reactions, studying the resolution of a short R-loop, and simulating RNA-scaffolded wireframe origami.

Assessment of A- to B- DNA Transitions Utilizing the Drude Polarizable Force Field

In addition to the well-characterized B-form of DNA, duplex DNA can adopt various conformations, such as A or Z-DNA. Though less common, these structures can be induced biologically through protein or ligand interactions or experimentally with niche environmental conditions, such as high salt concentrations or in mixed water-ethanol. Reproducing these alternate structures through molecular dynamics simulations in recent years has been quite challenging with the currently available force fields, simulation techniques, and time scales. In this study, the Drude polarizable force field is tested for its ability to facilitate transitions between A-DNA and B-DNA or maintain A-DNA. Though transitions away from B-DNA were observed in high concentrations of ethanol, the resulting structures had hybrid properties taken from both B-DNA and A-DNA structures. This was also true for A-DNA in ethanol, which lost some of the A-DNA properties that it was expected to maintain. When B-DNA was tested in high salt environments, the resulting B-DNA structures showed no distinguishable differences with the increasing salt concentrations tested. These results with the Drude FF and recent results with additive force fields suggest that at present the current additive and polarizable force fields do not facilitate a complete transition between B- to A-DNA conformations under the conditions simulated. At present, the Drude FF favors A-B DNA hybrid structures when simulated in nonphysiological conditions.

Highly nonlinear optic nucleic acid thin-solid film to generate short pulse laser

Using aqueous precursors, we report successfully fabricating thin-solid films of two nucleic acids, ribonucleic acid (RNA) and deoxyribonucleic acid (DNA). We investigated the potential of these films deposited on a fiber optic platform as all-fiber integrated saturable absorbers (SAs) for ultrafast nonlinear optics. RNA-SA performances were comparable to those of DNA-SA in terms of its nonlinear transmission, modulation depth, and saturation intensity. Upon insertion of these devices into an Erbium-doped fiber ring-laser cavity, both RNA and DNA SAs enabled efficient passive Q-switching operation. RNA-SA application further facilitated robust mode-locking and generated a transform-limited soliton pulse, exhibiting a pulse duration of 633 femtoseconds. A detailed analysis of these pulsed laser characteristics compared RNA and DNA fiber optic SAs with other nonlinear optic materials. The findings of this research establish the feasibility of utilizing RNA as a saturable absorber in ultrafast laser systems with an equal or higher potential as DNA, which presents novel possibilities for the nonlinear photonic applications of nucleic acid thin solid films.

Probing the structure of human tRNA3Lys in the presence of ligands using docking, MD simulations and MSM analysis

The tRNA3Lys, which acts as a primer for human immunodeficiency virus type 1 (HIV-1) reverse transcription, undergoes structural changes required for the formation of a primer-template complex. Small molecules have been targeted against tRNA3Lys to inhibit the primer-template complex formation. The present study aims to understand the kinetics of the conformational landscape spanned by tRNA3Lys in apo form using molecular dynamics simulations and Markov state modeling. The study is taken further to investigate the effect of small molecules like 1,4T and 1,5T on structural conformations and kinetics of tRNA3Lys, and comparative analysis is presented. Markov state modeling of tRNA3Lys apo resulted in three metastable states where the conformations have shown the non-canonical structures of the anticodon loop. Based on analyses of ligand-tRNA3Lys interactions, crucial ion and water mediated H-bonds and free energy calculations, it was observed that the 1,4-triazole more strongly binds to the tRNA3Lys compared to 1,5-triazole. However, the MSM analysis suggest that the 1,5-triazole binding to tRNA3Lys has brought rigidity not only in the binding pocket (TΨC arm, D-TΨC loop) but also in the whole structure of tRNA3Lys. This may affect the easy opening of primer tRNA3Lys required for HIV-1 reverse transcription.

Sequence-Dependent Properties of the RNA Duplex

Sequence-dependent properties of the DNA duplex have been accurately described using extensive molecular dynamics simulations. The RNA duplex meanwhile─which is typically represented as a sequence-averaged rigid rod─does not benefit from having equivalent molecular dynamics simulations. In this paper, we present a massive simulation effort using a set of ABC-optimized duplexes from which we derived tetramer-resolution properties of the RNA duplex and a simple mesoscopic model that can represent elastic properties of long RNA duplexes. Despite the extreme chemical similarity between DNA and RNA, the local and global elastic properties of the duplexes are very different. DNA duplexes show a complex and nonelastic pattern of flexibility, for instance, while RNA duplexes behave as an elastic system whose deformations can be represented by simple harmonic potentials. In RNA duplexes (RNA2), not only are intra- and interbase pair parameters (equilibrium and mechanical) different from those in the equivalent DNA duplex sequences (DNA2) but the correlations between movements also differ. Simple statements on the relative flexibility or stability of both polymers are meaningless and should be substituted by a more detailed description depending on the sequence and the type of deformation considered.

Self-assembly of a multimeric genomic hydrogel via multi-primed chain reaction of dual single-stranded circular plasmids for cell-free protein production

Recent technical advances in cell-free protein synthesis (CFPS) offer several advantages over cell-based expression systems, including the application of cellular machinery, such as transcription and translation, in the test tube. Inspired by the advantages of CFPS, we have fabricated a multimeric genomic DNA hydrogel (mGD-gel) via rolling circle chain amplification (RCCA) using dual single-stranded circular plasmids with multiple primers. The mGD-gel exhibited significantly enhanced protein yield. In addition, mGD-gel can be reused at least five times, and the shape of the mGD-gel can be easily manipulated without losing the feasibility of protein expression. The mGD-gel platform based on the self-assembly of multimeric genomic DNA strands (mGD strands) has the potential to be used in CFPS systems for a variety of biotechnological applications.

Sensitivity of the RNA Structure to Ion Conditions as Probed by Molecular Dynamics Simulations of Common Canonical RNA Duplexes

RNA molecules play a key role in countless biochemical processes. RNA interactions, which are of highly diverse nature, are determined by the fact that RNA is a highly negatively charged polyelectrolyte, which leads to intimate interactions with an ion atmosphere. Although RNA molecules are formally single-stranded, canonical (Watson-Crick) duplexes are key components of folded RNAs. A double-stranded (ds) RNA is also important for the design of RNA-based nanostructures and assemblies. Despite the fact that the description of canonical dsRNA is considered the least problematic part of RNA modeling, the imperfect shape and flexibility of dsRNA can lead to imbalances in the simulations of larger RNAs and RNA-containing assemblies. We present a comprehensive set of molecular dynamics (MD) simulations of four canonical A-RNA duplexes. Our focus was directed toward the characterization of the influence of varying ion concentrations and of the size of the solvation box. We compared several water models and four RNA force fields. The simulations showed that the A-RNA shape was most sensitive to the RNA force field, with some force fields leading to a reduced inclination of the A-RNA duplexes. The ions and water models played a minor role. The effect of the box size was negligible, and even boxes with a small fraction of the bulk solvent outside the RNA hydration sphere were sufficient for the simulation of the dsRNA.

A Machine Learning Approach to Identify the Importance of Novel Features for CRISPR/Cas9 Activity Prediction

The reprogrammable CRISPR/Cas9 genome editing tool's growing popularity is hindered by unwanted off-target effects. Efforts have been directed toward designing efficient guide RNAs as well as identifying potential off-target threats, yet factors that determine efficiency and off-target activity remain obscure. Based on sequence features, previous machine learning models performed poorly on new datasets, thus there is a need for the incorporation of novel features. The binding energy estimation of the gRNA-DNA hybrid as well as the Cas9-gRNA-DNA hybrid allowed generating better performing machine learning models for the prediction of Cas9 activity. The analysis of feature contribution towards the model output on a limited dataset indicated that energy features played a determining role along with the sequence features. The binding energy features proved essential for the prediction of on-target activity and off-target sites. The plateau, in the performance on unseen datasets, of current machine learning models could be overcome by incorporating novel features, such as binding energy, among others. The models are provided on GitHub (GitHub Inc., San Francisco, CA, USA).

Biophysical interaction between lanthanum chloride and (CG)n or (GC)n repeats: A reversible B-to-Z DNA transition

The transition from right-handed to left-handed DNA is not only acts as the controlling factor for switching gene expression but also has equal importance in designing nanomechanical devices. The (CG)_n and (GC)_n repeat sequences are well known model molecules to study B-Z transition in the presence of higher concentration of monovalent cations. In this communication, we report a cyclic transition in (CG)₆ DNA using millimolar concentration of trivalent lanthanide salt LaCl₃. The controlled and reversible transition was seen in (CG)₁₂, and (GC)₁₂ DNA employing CD spectroscopy. While LaCl₃ failed to induce B-Z transition in shorter oligonucleotides such as (CG)₃ and (GC)₃, a smooth B-Z transition was recorded for (CG)₆, (CG)₁₂ and (GC)₁₂ sequences. Interestingly, the phenomenon was reversible (Z-B transition) with addition of EDTA. Particularly, two rounds of cyclic transition (B-Z-B-Z-B) have been noticed in (CG)₆ DNA in presence of LaCl₃ and EDTA which strongly suggest that B-Z transition is reversible in short repeat sequences. Thermal melting and annealing behaviour of B-DNA are reversible while the thermal melting of LaCl₃-induced Z-DNA is irreversible which suggest a stronger binding of LaCl₃ to the phosphate backbone of Z-DNA. This was further supported by isothermal titration calorimetric study. Molecular dynamics (MD) simulation indicates that the mode of binding of La³⁺ (of LaCl₃) with d(CG)₈.d(CG)₈ is through the minor groove, wherein, 3 out of 11 La³⁺ bridge the anionic oxygens of the complementary strands. Such a tight coordination of La³⁺ with the anionic oxygens at the minor groove surface may be the reason for the experimentally observed irreversibility of LaCl₃-induced Z-DNA seen in longer DNA fragments. Thus, these results indicate LaCl₃ can easily be adopted as an inducer of left-handed DNA in other short oligonucleotides sequences to facilitate the understanding of the molecular mechanism of B-Z transition.

Iso-FRET: an isothermal competition assay to analyze quadruplex formation in vitro

Algorithms have been widely used to predict G-quadruplexes (G4s)-prone sequences. However, an experimental validation of these predictions is generally required. We previously reported a high-throughput technique to evidence G4 formation in vitro called FRET-MC. This method, while convenient and reproducible, has one known weakness: its inability to pin point G4 motifs of low thermal stability. As such quadruplexes may still be biologically relevant if formed at physiological temperature, we wanted to develop an independent assay to overcome this limitation. To this aim, we introduced an isothermal version of the competition assay, called iso-FRET, based on a duplex-quadruplex competition and a well-characterized bis-quinolinium G4 ligand, PhenDC3. G4-forming competitors act as decoys for PhenDC3, lowering its ability to stabilize the G4-forming motif reporter oligonucleotide conjugated to a fluorescence quencher (37Q). The decrease in available G4 ligand concentration restores the ability of 37Q to hybridize to its FAM-labeled short complementary C-rich strand (F22), leading to a decrease in fluorescence signal. In contrast, when no G4-forming competitor is present, PhenDC3 remains available to stabilize the 37Q quadruplex, preventing the formation of the F22 + 37Q complex. Iso-FRET was first applied to a reference panel of 70 sequences, and then used to investigate 23 different viral sequences.

Smart Metal-Organic Frameworks for Biotechnological Applications: A Mini-Review

In this pandemic situation it is evident that viruses and bacteria, more specifically, multiple drug resistant (MDR) bacteria, endanger human civilization severely. It is high time to design smart weapons to combat these pathogens for the prevention and cure of allied ailments. Metal-organic frameworks (MOFs) are porous materials designed from metal ions or inorganic clusters and multidentate organic ligands. Due to some unique features like high porosity, tunable pore shape and size, numerous possible metal-ligand combinations, etc., MOFs are ideal candidates to design "smart biotechnological tools". MOFs construct promising fluorescence based biosensing platforms for detection of viruses. MOFs also exhibit excellent antibacterial activity due to their ability for sustained release of active biocidal agents. There are several reviews that summarize the antibacterial applications of MOFs, but the biosensing platforms based on MOFs for detection of viruses have scarcely been summarized. This review carefully covers both the aspects including virus detection (nucleic acid recognition and immunological detection) with underlying mechanisms as well as antibacterial application of MOFs and doped MOFs or composites. This review will deliver valuable information and references for designing new, smarter antimicrobial agents based on MOFs.

Inhibition of off-target cleavage by RNase H using an artificial cationic oligosaccharide

Sequence-dependent off-target effects are a serious problem of antisense oligonucleotide-based drugs. Some of these side effects are induced by ribonuclease H (RNase H)-mediated cleavage of non-target RNAs with base sequences similar to that of the target RNA. We found that an artificial cationic oligosaccharide, ODAGal4, improved single-base discrimination for RNase H cleavage.

Transferable Ion Force Fields in Water from a Simultaneous Optimization of Ion Solvation and Ion-Ion Interaction

The poor performance of many existing nonpolarizable ion force fields is typically blamed on either the lack of explicit polarizability, the absence of charge transfer, or the use of unreduced Coulomb interactions. However, this analysis disregards the large and mostly unexplored parameter range offered by the Lennard-Jones potential. We use a global optimization procedure to develop water-model-transferable force fields for the ions K⁺, Na⁺, Cl^–, and Br^– in the complete parameter space of all Lennard-Jones interactions using standard mixing rules. No extra-thermodynamic assumption is necessary for the simultaneous optimization of the four ion pairs. After an optimization with respect to the experimental solvation free energy and activity, the force fields reproduce the concentration-dependent density, ionic conductivity, and dielectric constant with high accuracy. The force field is fully transferable between simple point charge/extended and transferable intermolecular potential water models. Our results show that a thermodynamically consistent force field for these ions needs only Lennard-Jones and standard Coulomb interactions.

Structural insight into locked nucleic acid based novel antisense modifications: A DFT calculations at monomer and MD simulations at oligomer level

In the present study, five novel LNA based antisense modifications have been proposed. A conformational search was carried out using TANGO, followed by geometry optimization using MOPAC. Based on their electronic energies the most stable conformation for each modification was identified. Further, DFT based full geometry optimization on the most stable conformations at the gas phase B3LYP/6-31G(d,p) using a Gaussian03 and single point energy calculations on the optimized structures at the solvent phase B3LYP/6-311G(d,p) level of theory were done to derive their quantum chemical descriptors using the Gaussian09. A comparison of global reactivity descriptors confirmed that the LNA based modifications were the most reactive. Base-pair stability was recorded by observing the binding energies and base-pairing conformations of modified GC base pairs at the B3LYP/6-311G(d,p) level of theory. Molecular dynamics simulations have been performed at the oligomer duplex level by incorporating individual modifications on 20-mer RNA-RNA duplexes using AMBER16. Free energy calculations of duplex structures suggested that incorporation of A2 modification into the RNA-RNA duplex increased the duplex binding affinity similar to LNA. Whereas, the A3 modification showed less binding compared to LNA but improved binding compared to MOE. This computational approach using quantum chemical methods may be very useful to propose better modifications than the existing ones before performing the experiments in the area of antisense technology.

Towards a Better Understanding of Computational Models for Predicting DNA Methylation Effects at the Molecular Level

Human DNA is a very sensitive macromolecule and slight changes in the structure of DNA can have disastrous effects on the organism. When nucleotides are modified, or changed, the resulting DNA sequence can lose its information, if it is part of a gene, or it can become a problem for replication and repair. Human cells can regulate themselves by using a process known as DNA methylation. This methylation is vitally important in cell differentiation and expression of genes. When the methylation is uncontrolled, however, or does not occur in the right place, serious pathophysiological consequences may result. Excess methylation causes changes in the conformation of the DNA double helix. The secondary structure of DNA is highly dependent upon the sequence. Therefore, if the sequence changes slightly the secondary structure can change as well. These slight changes will then cause the doublestranded DNA to be more open and available in some places where large adductions can come in and react with the DNA base pairs. Computer models have been used to simulate a variety of biological processes including protein function and binding, and there is a growing body of evidence that in silico methods can shed light on DNA methylation. Understanding the anomeric effect that contributes to the structural and conformational flexibility of furanose rings through a combination of quantum mechanical and experimental studies is critical for successful molecular dynamic simulations.

Improved nearest-neighbor parameters for the stability of RNA/DNA hybrids under a physiological condition

The stability of Watson–Crick paired RNA/DNA hybrids is important for designing optimal oligonucleotides for ASO (Antisense Oligonucleotide) and CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)–Cas9 techniques. Previous nearest-neighbour (NN) parameters for predicting hybrid stability in a 1 M NaCl solution, however, may not be applicable for predicting stability at salt concentrations closer to physiological condition (e.g. ∼100 mM Na⁺ or K⁺ in the presence or absence of Mg²⁺). Herein, we report measured thermodynamic parameters of 38 RNA/DNA hybrids at 100 mM NaCl and derive new NN parameters to predict duplex stability. Predicted ΔG°₃₇ and T_m values based on the established NN parameters agreed well with the measured values with 2.9% and 1.1°C deviations, respectively. The new results can also be used to make precise predictions for duplexes formed in 100 mM KCl or 100 mM NaCl in the presence of 1 mM Mg²⁺, which can mimic an intracellular and extracellular salt condition, respectively. Comparisons of the predicted thermodynamic parameters with published data using ASO and CRISPR–Cas9 may allow designing shorter oligonucleotides for these techniques that will diminish the probability of non-specific binding and also improve the efficiency of target gene regulation.

Metal-organic frameworks for virus detection

Virus severely endangers human life and health, and the detection of viruses is essential for the prevention and treatment of associated diseases. Metal-organic framework (MOF), a novel hybrid porous material which is bridged by the metal clusters and organic linkers, has become a promising biosensor platform for virus detection due to its outstanding properties including high surface area, adjustable pore size, easy modification, etc. However, the MOF-based sensing platforms for virus detection are rarely summarized. This review systematically divided the detection platforms into nucleic acid and immunological (antigen and antibody) detection, and the underlying sensing mechanisms were interpreted. The nucleic acid sensing was discussed based on the properties of MOF (such as metal ion, functional group, geometry structure, size, porosity, stability, etc.), revealing the relationship between the sensing performance and properties of MOF. Moreover, antibodies sensing based on the fluorescence detection and antigens sensing based on molecular imprinting or electrochemical immunoassay were highlighted. Furthermore, the remaining challenges and future development of MOF for virus detection were further discussed and proposed. This review will provide valuable references for the construction of sophisticated sensing platform for the detection of viruses, especially the 2019 coronavirus.

Structure-mechanics statistical learning unravels the linkage between local rigidity and global flexibility in nucleic acids

The mechanical properties of nucleic acids underlie biological processes ranging from genome packaging to gene expression, but tracing their molecular origin has been difficult due to the structural and chemical complexity. We posit that concepts from machine learning can help to tackle this long-standing challenge. Here, we demonstrate the feasibility and advantage of this strategy through developing a structure-mechanics statistical learning scheme to elucidate how local rigidity in double-stranded (ds)DNA and dsRNA may lead to their global flexibility in bend, stretch, and twist. Specifically, the mechanical parameters in a heavy-atom elastic network model are computed from the trajectory data of all-atom molecular dynamics simulation. The results show that the inter-atomic springs for backbone and ribose puckering in dsRNA are stronger than those in dsDNA, but are similar in strengths for base-stacking and base-pairing. Our analysis shows that the experimental observation of dsDNA being easier to bend but harder to stretch than dsRNA comes mostly from the respective B- and A-form topologies. The computationally resolved composition of local rigidity indicates that the flexibility of both nucleic acids is mostly due to base-stacking. But for properties like twist-stretch coupling, backbone springs are shown to play a major role instead. The quantitative connection between local rigidity and global flexibility sets foundation for understanding how local binding and chemical modification of genetic materials effectuate longer-ranged regulatory signals.

The mechanical properties of nucleic acids underlie biological processes ranging from genome packaging to gene expression. We devise structural mechanics statistical learning method to reveal their molecular origin in terms of chemical interactions.

Synthesis and properties of DNA oligomers containing stereopure phosphorothioate linkages and C-5 modified deoxyuridine derivatives

Phosphorothioate (PS) modification, where a non-bridging oxygen atom in a phosphodiester linkage is replaced by a sulfur atom, is widely used to improve the properties of nucleic acid drugs. Each PS linkage can be found in two stereoisomers, Rp and Sp. Since one non-bridging oxygen or sulfur atom in Sp-PS or Rp-PS, respectively, is located close to the C-5 substituent of uracil in a DNA/RNA hybrid duplex, the combination of the stereochemistry of the PS linkages and the type of the C-5 modification of uracil bases is expected to affect the properties of the hybrid duplexes. Herein, DNA oligomers containing both stereopure phosphorohioate linkages and C-5 modified deoxyuridine derivatives were synthesized. The thermodynamic stability of the DNA/RNA and DNA/DNA duplexes and RNase H activity of the DNA/RNA duplexes were evaluated. The combination of 5-propynyluracil and (Rp)-PS linkages in a DNA strand could significantly increase the thermal stability of a DNA/RNA hybrid duplex without reducing its RNase H activity.

The combination of 5-propynyluracil and (Rp)-PS linkages in a DNA strand could significantly increase the thermal stability of a DNA/RNA hybrid duplex.

The Mechanical Properties of RNA-DNA Hybrid Duplex Stretched by Magnetic Tweezers

RNA can anneal to its DNA template to generate an RNA-DNA hybrid (RDH) duplex and a displaced DNA strand, termed R-loop. RDH duplex occupies up to 5% of the mammalian genome and plays important roles in many biological processes. The functions of RDH duplex are affected by its mechanical properties, including the elasticity and the conformation transitions. The mechanical properties of RDH duplex, however, are still unclear. In this work, we studied the mechanical properties of RDH duplex using magnetic tweezers in comparison with those of DNA and RNA duplexes with the same sequences. We report that the contour length of RDH duplex is ∼0.30 nm/bp, and the stretching modulus of RDH duplex is ∼660 pN, neither of which is sensitive to NaCl concentration. The persistence length of RDH duplex depends on NaCl concentration, decreasing from ∼63 nm at 1 mM NaCl to ∼49 nm at 500 mM NaCl. Under high tension of ∼60 pN, the end-opened RDH duplex undergoes two distinct overstretching transitions; at high salt in which the basepairs are stable, it undergoes the nonhysteretic transition, leading to a basepaired elongated structure, whereas at low salt, it undergoes a hysteretic peeling transition, leading to the single-stranded DNA strand under force and the single-stranded RNA strand coils. The peeled RDH is difficult to reanneal back to the duplex conformation, which may be due to the secondary structures formed in the coiled single-stranded RNA strand. These results help us understand the full picture of the structures and mechanical properties of nucleic acid duplexes in solution and provide a baseline for studying the interaction of RDH with proteins at the single-molecule level.

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.

Molecular dynamics simulations of theoretical cellulose nanotube models

Nanotubes are remarkable nanoscale architectures for a wide range of potential applications. In the present paper, we report a molecular dynamics (MD) study of the theoretical cellulose nanotube (CelNT) models to evaluate their dynamic behavior in solution (either chloroform or benzene). Based on the one-quarter chain staggering relationship, we constructed six CelNT models by combining the two chain polarities (parallel (P) and antiparallel (AP)) and three symmetry operations (helical right (H_R), helical left (H_L), and rotation (R)) to generate a circular arrangement of molecular chains. Among the four models that retained the tubular form (P-H_R, P-H_L, P-R, and AP-R), the P-R and AP-R models have the lowest steric energies in benzene and chloroform, respectively. The structural features of the CelNT models were characterized in terms of the hydroxymethyl group conformation and intermolecular hydrogen bonds. Solvent structuring more clearly occurred with benzene than chloroform, suggesting that the CelNT models may disperse in benzene.

Dynamics changes of CRISPR-Cas9 systems induced by high fidelity mutations

CRISPR-Cas9, a powerful genome editing tool, has widely been applied in biological fields.

Biophysical properties, thermal stability and functional impact of 8-oxo-7,8-dihydroguanine on oligonucleotides of RNA-a study of duplex, hairpins and the aptamer for preQ1 as models

A better understanding of the effects that oxidative lesions have on RNA is of importance to understand their role in the development/progression of disease. 8-oxo-7,8-dihydroguanine was incorporated into RNA to understand its structural and functional impact on RNA:RNA and RNA:DNA duplexes, hairpins and pseudoknots. One to three modifications were incorporated into dodecamers of RNA [AAGAGGGAUGAC] resulting in thermal destabilization (ΔT_m – 10°C per lesion). Hairpins with tetraloops c-UUCG*-g* (8-10), a-ACCG-g* (11-12), c-UUG*G*-g* (13-16) and c-ACG*G*-g* (17-20) were modified and used to determine thermal stabilities, concluding that: (i) modifying the stem leads to destabilization unless adenosine is the opposing basepair of 8-oxoGua; (ii) modification at the loop is position- and sequence-dependent and varies from slight stabilization to large destabilization, in some cases leading to formation of other secondary structures (hairpin→duplex). Functional effects were established using the aptamer for preQ₁ as model. Modification at G5 disrupted the stem P1 and inhibited recognition of the target molecule 7-methylamino-7-deazaguanine (preQ₁). Modifying G11 results in increased thermal stability, albeit with a K_d 4-fold larger than its canonical analog. These studies show the capability of 8-oxoG to affect structure and function of RNA, resulting in distinct outcomes as a function of number and position of the lesion.

DNA recognition by an RNA-guided bacterial Argonaute

Argonaute (Ago) proteins are widespread in prokaryotes and eukaryotes and share a four-domain architecture capable of RNA- or DNA-guided nucleic acid recognition. Previous studies identified a prokaryotic Argonaute protein from the eubacterium Marinitoga piezophila (MpAgo), which binds preferentially to 5′-hydroxylated guide RNAs and cleaves single-stranded RNA (ssRNA) and DNA (ssDNA) targets. Here we present a 3.2 Å resolution crystal structure of MpAgo bound to a 21-nucleotide RNA guide and a complementary 21-nucleotide ssDNA substrate. Comparison of this ternary complex to other target-bound Argonaute structures reveals a unique orientation of the N-terminal domain, resulting in a straight helical axis of the entire RNA-DNA heteroduplex through the central cleft of the protein. Additionally, mismatches introduced into the heteroduplex reduce MpAgo cleavage efficiency with a symmetric profile centered around the middle of the helix. This pattern differs from the canonical mismatch tolerance of other Argonautes, which display decreased cleavage efficiency for substrates bearing sequence mismatches to the 5′ region of the guide strand. This structural analysis of MpAgo bound to a hybrid helix advances our understanding of the diversity of target recognition mechanisms by Argonaute proteins.

Focus on PNA Flexibility and RNA Binding using Molecular Dynamics and Metadynamics

Peptide Nucleic Acids (PNAs) can efficiently target DNA or RNA acting as chemical tools for gene regulation. Their backbone modification and functionalization is often used to increase the affinity for a particular sequence improving selectivity. The understanding of the trading forces that lead the single strand PNA to bind the DNA or RNA sequence is preparatory for any further rational design, but a clear and unique description of this process is still not complete. In this paper we report further insights into this subject, by a computational investigation aiming at the characterization of the conformations of a single strand PNA and how these can be correlated to its capability in binding DNA/RNA. Employing Metadynamics we were able to better define conformational pre-organizations of the single strand PNA and γ-modified PNA otherwise unrevealed through classical molecular dynamics. Our simulations driven on backbone modified PNAs lead to the conclusion that this γ-functionalization affects the single strand preorganization and targeting properties to the DNA/RNA, in agreement with circular dichroism (CD) spectra obtained for this class of compounds. MD simulations on PNA:RNA dissociation and association mechanisms allowed to reveal the critical role of central bases and preorganization in the binding process.

Synthesis, Thermal Stability, Biophysical Properties, and Molecular Modeling of Oligonucleotides of RNA Containing 2'-O-2-Thiophenylmethyl Groups

Dodecamers of RNA [CUACGGAAUCAU] were functionalized with C2'-O-2-thiophenylmethyl groups to obtain oligonucleotides 10-14 and 17. The modified nucleotides were incorporated into RNA strands via solid-phase synthesis. The biophysical properties of these ONs were used to quantify the effects of this modification on RNA:RNA and RNA:DNA duplexes. A combination of UV-vis and circular dichroism were used to determine thermal stabilities of all strands, which hybridized into A-form geometries. Destabilization of the double stranded RNA was measured as a function of number of consecutive modifications, reflected in decreased thermal denaturation values (ΔT_m, ca. 2.5-11.5 °C). Van't Hoff plots on a duplex containing one modification (10:15) displayed a ca. ΔΔG° of +4 kcal/mol with respect to its canonical analogue. Interestingly, hybridization of two modified strands (13:17, containing a total of eight modifications) resulted in increased stability and a distinct secondary structure, reflected in its CD spectrum. Molecular modeling based on DFT calculations shed light on the nature of this stability, with induced changes in the torsional angle δ (C5'-C4'-C3'-O3) and phosphate-phosphate distances that are in agreement with a compacted structure. The described synthetic methodology and structural information will be useful in the design of thermodynamically stable structures containing chemically reactive modifications.

Long-range electrostatic corrections in multipolar/polarizable QM/MM simulations

Taking long-range electrostatic effects into account in classical and hybrid quantum mechanics-molecular mechanics (QM/MM) simulations is necessary for an accurate description of the system under study. We have recently developed a method, termed long-range electrostatic corrections (LREC), for monopolar QM/MM calculations. Here, we present an extension of LREC for multipolar/polarizable QM/MM simulations within the LICHEM software package. Reaction barriers and QM-MM interaction energies converge with a LREC cutoff between 20 and 25 Å, in agreement with our previous results. Additionally, the LREC approach for the QM-MM interactions can be smoothly combined with standard shifting or Ewald summation methods in the MM calculations. We recommend the use of QM(LREC)/MM(PME), where the QM region is treated with LREC and the MM region is treated with particle mesh Ewald (PME) summation. This combination is an excellent compromise between simplicity, speed, and accuracy for large QM/MM simulations.

Biophysical properties and thermal stability of oligonucleotides of RNA containing 7,8-dihydro-8-hydroxyadenosine

Circular dichroism (CD) was used to assess the stabilization/destabilization imposed by oxidative lesion 7,8-dihydro-8-hydroxyadenosine (8-oxoA) on strands of RNA with different structural motifs. RNA:RNA homoduplex destabilization was observed in a position dependent manner using 10-mers as models that displayed differences between 12.7 and 15.1°C. We found that increasing the number of modifications resulted in depressed Tm values of about 12-15°C per lesion. The same effect was observed on RNA:DNA heteroduplex samples. We also tested the effects of this lesion in short hairpins containing the tetraloop UUCX (X = A, 8-oxoA). We found that the stem was hypersensitive to substitution of A by 8-oxoA and that it destabilized the structure by >23°C. Concomitant substitution at the stem and loop prevented formation of this secondary structure or lead to other less-stable hairpins. Incorporation of this lesion at the first base of the loop had no effect on either structure. Overall, we found that the effects of 8-oxoA on RNA structure are position dependent and that its stabilization may vary from sharp decreases to small increments, in some cases, leading to the formation of other more/less stable structures. These structural changes may have larger biological implications, particularly if the oxidatively modified RNA persists, thus leading to changes in RNA reactivity and function.

Discovery and Structural Characterization of G-quadruplex DNA in Human Acetyl-CoA Carboxylase Gene Promoters: Its Role in Transcriptional Regulation and as a Therapeutic Target for Human Disease

Accumulating evidence suggests that G-quadruplexes play vital roles in gene expression, DNA replication, and recombination. Three distinct promoters (PI, PII, and PIII) regulate human acetyl-CoA carboxylase 1 (ACC1) gene expression. In this study, we asked whether the G-rich sequences within the human ACC1 (PI and PII) promoters can form G-quadruplex structures and regulate normal DNA transactions. Using multiple complementary methods, we show that G-rich sequences of PI and PII promoters form intramolecular G-quadruplex structures and then establish unambiguously the topologies of these structures. Importantly, G-quadruplex formation in ACC1 gene promoter region blocks DNA replication and suppresses transcription, and this effect was further augmented by G-quadruplex stabilizing ligands. Altogether, these results are consistent with the notion that G-quadruplex structures exist within the human ACC1 gene promoter region, whose activity can be suppressed by G-quadruplex stabilizing ligands, thereby revealing a novel regulatory mechanism of ACC1 gene expression and as a possible therapeutic target.

Platforms Formed from a Three-Dimensional Cu-Based Zwitterionic Metal-Organic Framework and Probe ss-DNA: Selective Fluorescent Biosensors for Human Immunodeficiency Virus 1 ds-DNA and Sudan Virus RNA Sequences

We herein report a water-stable three-dimensional Cu-based metal-organic framework (MOF) 1 supported by a tritopic quaternized carboxylate and 4,4'-dipyridyl sulfide as an ancillary ligand. This MOF exhibits unique pore shapes with aromatic rings, positively charged pyridinium and unsaturated Cu(II) cation centers, free carboxylates, tessellating H2O, and coordinating SO4(2-) on the pore surface. Compound 1 can interact with two carboxyfluorescein (FAM)-labeled single-stranded DNA sequences (probe ss-DNA, delineated as P-DNA) through electrostatic, π-stacking, and/or hydrogen-bonding interactions to form two P-DNA@1 systems, and thus quench the fluorescence of FAM via a photoinduced electron-transfer process. These P-DNA@1 systems can be used as effective fluorescent sensors for human immunodeficiency virus 1 double-stranded DNA and Sudan virus RNA sequences, respectively, with detection limits of 196 and 73 pM, respectively.

Refinement of Generalized Born Implicit Solvation Parameters for Nucleic Acids and Their Complexes with Proteins

The Generalized Born (GB) implicit solvent model has undergone significant improvements in accuracy for modeling of proteins and small molecules. However, GB still remains a less widely explored option for nucleic acid simulations, in part because fast GB models are often unable to maintain stable nucleic acid structures or they introduce structural bias in proteins, leading to difficulty in application of GB models in simulations of protein-nucleic acid complexes. Recently, GB-neck2 was developed to improve the behavior of protein simulations. In an effort to create a more accurate model for nucleic acids, a similar procedure to the development of GB-neck2 is described here for nucleic acids. The resulting parameter set significantly reduces absolute and relative energy error relative to Poisson-Boltzmann for both nucleic acids and nucleic acid-protein complexes, when compared to its predecessor GB-neck model. This improvement in solvation energy calculation translates to increased structural stability for simulations of DNA and RNA duplexes, quadruplexes, and protein-nucleic acid complexes. The GB-neck2 model also enables successful folding of small DNA and RNA hairpins to near native structures as determined from comparison with experiment. The functional form and all required parameters are provided here and also implemented in the AMBER software.

A new smoothing function to introduce long-range electrostatic effects in QM/MM calculations

A new method to account for long range electrostatic contributions is proposed and implemented for quantum mechanics/molecular mechanics long range electrostatic correction (QM/MM-LREC) calculations. This method involves the use of the minimum image convention under periodic boundary conditions and a new smoothing function for energies and forces at the cutoff boundary for the Coulomb interactions. Compared to conventional QM/MM calculations without long-range electrostatic corrections, the new method effectively includes effects on the MM environment in the primary image from its replicas in the neighborhood. QM/MM-LREC offers three useful features including the avoidance of calculations in reciprocal space (k-space), with the concomitant avoidance of having to reproduce (analytically or approximately) the QM charge density in k-space, and the straightforward availability of analytical Hessians. The new method is tested and compared with results from smooth particle mesh Ewald (PME) for three systems including a box of neat water, a double proton transfer reaction, and the geometry optimization of the critical point structures for the rate limiting step of the DNA dealkylase AlkB. As with other smoothing or shifting functions, relatively large cutoffs are necessary to achieve comparable accuracy with PME. For the double-proton transfer reaction, the use of a 22 Å cutoff shows a close reaction energy profile and geometries of stationary structures with QM/MM-LREC compared to conventional QM/MM with no truncation. Geometry optimization of stationary structures for the hydrogen abstraction step by AlkB shows some differences between QM/MM-LREC and the conventional QM/MM. These differences underscore the necessity of the inclusion of the long-range electrostatic contribution.

Understanding the structural and dynamic consequences of DNA epigenetic modifications: computational insights into cytosine methylation and hydroxymethylation

We report a series of molecular dynamics (MD) simulations of up to a microsecond combined simulation time designed to probe epigenetically modified DNA sequences. More specifically, by monitoring the effects of methylation and hydroxymethylation of cytosine in different DNA sequences, we show, for the first time, that DNA epigenetic modifications change the molecule's dynamical landscape, increasing the propensity of DNA toward different values of twist and/or roll/tilt angles (in relation to the unmodified DNA) at the modification sites. Moreover, both the extent and position of different modifications have significant effects on the amount of structural variation observed. We propose that these conformational differences, which are dependent on the sequence environment, can provide specificity for protein binding.

Roles of the amino group of purine bases in the thermodynamic stability of DNA base pairing

The energetic aspects of hydrogen-bonded base-pair interactions are important for the design of functional nucleotide analogs and for practical applications of oligonucleotides. The present study investigated the contribution of the 2-amino group of DNA purine bases to the thermodynamic stability of oligonucleotide duplexes under different salt and solvent conditions, using 2'-deoxyriboinosine (I) and 2'-deoxyribo-2,6-diaminopurine (D) as non-canonical nucleotides. The stability of DNA duplexes was changed by substitution of a single base pair in the following order: G•C > D•T ≈ I•C > A•T > G•T > I•T. The apparent stabilization energy due to the presence of the 2-amino group of G and D varied depending on the salt concentration, and decreased in the water-ethanol mixed solvent. The effects of salt concentration on the thermodynamics of DNA duplexes were found to be partially sequence-dependent, and the 2-amino group of the purine bases might have an influence on the binding of ions to DNA through the formation of a stable base-paired structure. Our results also showed that physiological salt conditions were energetically favorable for complementary base recognition, and conversely, low salt concentration media and ethanol-containing solvents were effective for low stringency oligonucleotide hybridization, in the context of conditions employed in this study.

Structural and energetic analysis of 2-aminobenzimidazole inhibitors in complex with the hepatitis C virus IRES RNA using molecular dynamics simulations

Despite the many biological functions of RNA, very few drugs have been designed or found to target RNA. Here we report the results of molecular dynamics (MD) simulations and binding energy analyses on hepatitis C virus internal ribosome entry site (IRES) RNA in complex with highly charged 2-aminobenzimidazole inhibitors. Initial coordinates were taken from NMR and crystallography studies that had yielded different binding modes. During MD simulations, the RNA-inhibitor complex is stable in the crystal conformation but not in the NMR conformation. Additionally, we found that existing and standard MD trajectory postprocessing free energy methods, such as the MM-GBSA and MM-PBSA approaches available in AMBER, seem unsuitable to properly rank the binding energies of complexes between highly charged molecules. A better correlation with the experimental data was found using a rather simple binding enthalpy calculation based on the explicitly solvated potential energies. In anticipation of further growth in the use of small molecules to target RNA, we include results addressing the impact of charge assignment on docking, the structural role of magnesium in the IRES-inhibitor complex, the entropic contribution to binding energy, and simulations of a plausible scaffold design for new inhibitors.

An investigation of G-quadruplex structural polymorphism in the human telomere using a combined approach of hydrodynamic bead modeling and molecular dynamics simulation

Guanine-rich oligonucleotides can adopt noncanonical tertiary structures known as G-quadruplexes, which can exist in different forms depending on experimental conditions. High-resolution structural methods, such as X-ray crystallography and NMR spectroscopy, have been of limited usefulness in resolving the inherent structural polymorphism associated with G-quadruplex formation. The lack of, or the ambiguous nature of, currently available high-resolution structural data, in turn, has severely hindered investigations into the nature of these structures and their interactions with small-molecule inhibitors. We have used molecular dynamics in conjunction with hydrodynamic bead modeling to study the structures of the human telomeric G-quadruplex-forming sequences at the atomic level. We demonstrated that molecular dynamics can reproduce experimental hydrodynamic measurements and thus can be a powerful tool in the structural study of existing G-quadruplex sequences or in the prediction of new G-quadruplex structures.

Substrate recognition and specificity of double-stranded RNA binding proteins

Recognition of double-stranded (ds) RNA is an important part of many cellular pathways, including RNA silencing, viral recognition, RNA editing, processing, and transport. dsRNA recognition is often achieved by dsRNA binding domains (dsRBDs). We use atomistic molecular dynamics simulations to examine the binding interface of the transactivation response RNA binding protein (TRBP) dsRBDs to dsRNA substrates. Our results explain the exclusive selectivity of dsRBDs toward dsRNA and against DNA–RNA hybrid and dsDNA duplexes. We also provide corresponding experimental evidence. The dsRNA duplex is recognized by dsRBDs through the A-form of three duplex grooves and by the chemical properties of RNA bases, which have 2′-hydroxyl groups on their sugar rings. Our simulations show that TRBP dsRBD discriminates dsRNA- from DNA-containing duplexes primarily through interactions at two duplex grooves. The simulations also reveal that the conformation of the DNA–RNA duplex can be altered by dsRBD proteins, resulting in a weak binding of dsRBDs to DNA–RNA hybrids. Our study reveals the structural and molecular basis of protein–RNA interaction that gives rise to the observed substrate specificity of dsRNA binding proteins.

Molecular dynamics analysis of the interaction between the human BCL6 BTB domain and its SMRT, NcoR and BCOR corepressors: the quest for a consensus dynamic pharmacophore

Targeting the BCL6 protein is a promising therapeutic strategy for the treatment of B cell lymphomas. One approach to treat these diseases consists of finding drug candidates able to disrupt the interactions established between BCL6 and its corepressors. Thus, this work presents a thorough comparative analysis of the interactions between the BCL6 BTB (bric-a-brac tramtrack broad complex) protein domain and its SMRT, NcoR and BCOR corepressor BBDs (BCL6 binding domain) through molecular dynamics. Moreover, a theoretical structure is presented and checked for the BCL6(BTB)-NcoR(BBD) complex. Considering the BBDs to be composed of 17 amino acids, our analyses show the region involving residues 4-15 of these 17 to play a main role in the protein-corepressor interactions. Particularly SER(11) seems to have a high relevance as it establishes specific bonds with BCL6(BTB) and is one of the only two residues sequence equivalent for the three studied corepressors. From this study, 14 pharmacophoric points have been proposed divided in two groups which coincide with residues 4-11 and 11-15, being SER(11) a hinge point. This finding suggests the possibility of searching for 2 small molecule inhibitors, mimicking 8 and 7 pharmacophoric points, respectively, which could incorporate a hydrogen donor pharmacophoric point mimicking SER(11) in any or both molecules. In short, the present work aims to contribute further knowledge in the modeling of drugs mimicking BCL6(BTB)-corepressor complexes.

Detection of cyclin D1 mRNA by hybridization sensitive NIC-oligonucleotide probe

A large group of fluorescent hybridization probes, includes intercalating dyes for example thiazole orange (TO). Usually TO is coupled to nucleic acids post-synthetically which severely limits its use. Here, we have developed a phosphoramidite monomer, 10, and prepared a 2'-OMe-RNA probe, labeled with 5-(trans-N-hexen-1-yl-)-TO-2'-deoxy-uridine nucleoside, dU(TO), (Nucleoside bearing an Inter-Calating moiety, NIC), for selective mRNA detection. We investigated a series of 15-mer 2'-OMe-RNA probes, targeting the cyclin D1 mRNA, containing one or several dU(TO) at various positions. dU(TO)-2'-OMe-RNA exhibited up to 7-fold enhancement of TO emission intensity upon hybridization with the complementary RNA versus that of the oligomer alone. This NIC-probe was applied for the specific detection of a very small amount of a breast cancer marker, cyclin D1 mRNA, in total RNA extract from cancerous cells (250 ng/μl). Furthermore, this NIC-probe was found to be superior to our related NIF (Nucleoside with Intrinsic Fluorescence)-probe which could detect cyclin D1 mRNA target only at high concentrations (1840 ng/μl). Additionally, dU(T) can be used as a monomer in solid-phase oligonucleotide synthesis, thus avoiding the need for post-synthetic modification of oligonucleotide probes. Hence, we propose dU(TO) oligonucleotides, as hybridization probes for the detection of specific RNA in homogeneous solutions and for the diagnosis of breast cancer.