Nucleic acid duplex stability: influence of base composition on cation effects

Quantitative Effects of the Loop Region on Topology, Thermodynamics, and Cation Binding of DNA G-quadruplexes

The thermal stability of G-quadruplexes is important for their biological roles. G-quadruplexes are stable in the presence of cations such as K+ and Na+ because these cations coordinate in the G-quartet of four guanine bases. It is well known that the number of G-quartets and the configuration of the guanine bases affect the binding affinity of the cation. Recently, structures formed in the loop regions connecting the guanine stretches have attracted significant attention, because the loop region affects G-quadruplex properties, such as topology, thermal stability, and interactions with proteins and small molecules. Considering these effects, the loop region can also affect the binding affinity of the cations. Here, we designed a series of G-quadruplex-forming DNA sequences that contain a hairpin in a loop region and investigated the effects of the sequence and structure of the loop region on the cation binding affinity as well as the thermal stability of the G-quadruplex as a whole. First, structural analysis of the DNA sequences showed that the hairpin at the loop plays a key role in determining G4 topology (strand orientation). Second, in the case of the G-quadruplexes with the hairpin-forming loop region, it was found that a longer loop length led to a higher thermodynamic stability of the G-quadruplex as well as higher cation binding affinity. In contrast, an unstructured loop region did not lead to such effects. Interestingly, the cation binding affinity was correlated to the thermodynamic stability of the hairpin structure at the loop region. It was quantitatively demonstrated that the stable loop region stabilized the whole G-quadruplex structure, which induced higher cation binding affinity. These systematic and quantitative results showed that the loop region is one of the determinants of cation binding and expanded the possibilities of drug development targeting G4s by stabilizing the loop region.

Cooperative control of a DNA origami force sensor

Biomolecular systems are dependent on a complex interplay of forces. Modern force spectroscopy techniques provide means of interrogating these forces, but they are not optimized for studies in constrained environments as they require attachment to micron-scale probes such as beads or cantilevers. Nanomechanical devices are a promising alternative, but this requires versatile designs that can be tuned to respond to a wide range of forces. We investigate the properties of a nanoscale force sensitive DNA origami device which is highly customizable in geometry, functionalization, and mechanical properties. The device, referred to as the NanoDyn, has a binary (open or closed) response to an applied force by undergoing a reversible structural transition. The transition force is tuned with minor alterations of 1 to 3 DNA oligonucleotides and spans tens of picoNewtons (pN). The DNA oligonucleotide design parameters also strongly influence the efficiency of resetting the initial state, with higher stability devices (≳10 pN) resetting more reliably during repeated force-loading cycles. Finally, we show the opening force is tunable in real time by adding a single DNA oligonucleotide. These results establish the potential of the NanoDyn as a versatile force sensor and provide fundamental insights into how design parameters modulate mechanical and dynamic properties.

Competitive Influence of Alkali Metals in the Ion Atmosphere on Nucleic Acid Duplex Stability

The nonspecific atmosphere around nucleic acids, often termed the ion atmosphere, encompasses a collection of weak ion-nucleic acid interactions. Although nonspecific, the ion atmosphere has been shown to influence nucleic acid folding and structural stability. Studies investigating the composition of the ion atmosphere have shown competitive occupancy of the atmosphere between metal ions in the same solution. Many studies have investigated single ion effects on nucleic acid secondary structure stability; however, no comprehensive studies have investigated how the competitive occupancy of mixed ions in the ion atmosphere influences nucleic acid secondary structure stability. Here, six oligonucleotides were optically melted in buffers containing molar quantities, or mixtures, of either XCl (X = Li, K, Rb, or Cs) or NaCl. A correction factor was developed to better predict RNA duplex stability in solutions containing mixed XCl/NaCl. For solutions containing a 1:1 mixture of XCl/NaCl, one alkali metal chloride contributed more to duplex stability than the other. Overall, there was a 54% improvement in predictive capabilities with the correction factor compared with the standard 1.0 M NaCl nearest-neighbor models. This correction factor can be used in models to better predict RNA secondary structure in solutions containing mixed XCl/NaCl.

Detection of a 7SL RNA-derived small non-coding RNA using Molecular Beacons in vitro and in cells

Small non-coding RNAs (sncRNA) are involved in many steps of the gene expression cascade and regulate processing and expression of mRNAs by the formation of ribonucleoprotein complexes (RNP) such as the RNA-induced silencing complex (RISC). By analyzing small RNA Seq data sets, we identified a sncRNA annotated as piR-hsa-1254, which is likely derived from the 3'-end of 7SL RNA2 (RN7SL2), herein referred to as snc7SL RNA. The 7SL RNA is an abundant long non-coding RNA polymerase III transcript and serves as structural component of the cytoplasmic signal recognition particle (SRP). To evaluate a potential functional role of snc7SL RNA, we aimed to define its cellular localization by live cell imaging. Therefore, a Molecular Beacon (MB)-based method was established to compare the subcellular localization of snc7SL RNA with its precursor 7SL RNA. We designed and characterized several MBs in vitro and tested those by live cell fluorescence microscopy. Using a multiplex approach, we show that 7SL RNA localizes mainly to the endoplasmic reticulum (ER), as expected for the SRP, whereas snc7SL RNA predominately localizes to the nucleus. This finding suggests a fundamentally different function of 7SL RNA and its derivate snc7SL RNA.

Stability of DNA Origami Nanostructures in Physiological Media: The Role of Molecular Interactions

Programmable, custom-shaped, and nanometer-precise DNA origami nanostructures have rapidly emerged as prospective and versatile tools in bionanotechnology and biomedicine. Despite tremendous progress in their utilization in these fields, essential questions related to their structural stability under physiological conditions remain unanswered. Here, DNA origami stability is explored by strictly focusing on distinct molecular-level interactions. In this regard, the fundamental stabilizing and destabilizing ionic interactions as well as interactions involving various enzymes and other proteins are discussed, and their role in maintaining, modulating, or decreasing the structural integrity and colloidal stability of DNA origami nanostructures is summarized. Additionally, specific issues demanding further investigation are identified. This review - through its specific viewpoint - may serve as a primer for designing new, stable DNA objects and for adapting their use in applications dealing with physiological media.

Cooperative control of a DNA origami force sensor

Most biomolecular systems are dependent on a complex interplay of forces. Modern force spectroscopy techniques provide means of interrogating these forces. These techniques, however, are not optimized for studies in constrained or crowded environments as they typically require micron-scale beads in the case of magnetic or optical tweezers, or direct attachment to a cantilever in the case of atomic force microscopy. We implement a nanoscale force-sensing device using a DNA origami which is highly customizable in geometry, functionalization, and mechanical properties. The device, referred to as the NanoDyn, functions as a binary (open or closed) force sensor that undergoes a structural transition under an external force. The transition force is tuned with minor alterations of 1 to 3 DNA oligonucleotides and spans tens of picoNewtons (pN). This actuation of the NanoDyn is reversible and the design parameters strongly influence the efficiency of resetting the initial state, with higher stability devices (≳10 pN) resetting more reliably during repeated force-loading cycles. Finally, we show that the opening force can be adjusted in real time by the addition of a single DNA oligonucleotide. These results establish the NanoDyn as a versatile force sensor and provide fundamental insights into how design parameters modulate mechanical and dynamic properties.

Molecular Recognition in Confined Space Elucidated with DNA Nanopores and Single-Molecule Force Microscopy

The binding of ligands to receptors within a nanoscale small space is relevant in biology, biosensing, and affinity filtration. Binding in confinement can be studied with biological systems but under the limitation that essential parameters cannot be easily controlled including receptor type and position within the confinement and its dimensions. Here we study molecular recognition with a synthetic confined nanopore with controllable pore dimension and molecular DNA receptors at different depth positions within the channel. Binding of a complementary DNA strand is studied at the single-molecule level with atomic force microscopy. Following the analysis, kinetic association rates are lower for receptors positioned deeper inside the pore lumen while dissociation is faster and requires less force. The phenomena are explained by the steric constraints on molecular interactions in confinement. Our study is the first to explore recognition in DNA nanostructures with atomic force microscopy and lays out new tools to further quantify the effect of nanoconfinement on molecular interactions.

Population Genetic Study on the European Flounder (Platichthys flesus) from the Southern Baltic Sea Using SNPs and Microsatellite Markers

The European flounder (Platichthys flesus), which is closely related to the recently discovered Baltic flounder (Platichthys solemdali), is currently the third most commercially fished species in the Baltic Sea. According to the available data from the Polish Fisheries Monitoring Center and fishermen's observations, the body condition indices of the species in the Baltic Sea have declined in recent years. The aim of the present study was to obtain information on the current patterns of genetic variability and the population structure of the European flounder and to verify whether the Baltic flounder is present in the southern Baltic Sea. Moreover, we aimed to verify whether the observed decline in the body condition indices of the species in the Baltic Sea might be associated with adaptive alterations in its gene pool due to increased fishing pressure. For this purpose, 190 fish were collected from four locations along the central coastline of Poland, i.e., Mechelinki, Władysławowo, the Vistula Lagoon in 2018, and the Słupsk Bank in 2020. The fish were morphologically analyzed and then genetically screened by the application of nineteen microsatellite DNA and two diagnostic SNP markers. The examined European flounder specimens displayed a high level of genetic diversity (PIC = 0.832-0.903, I = 2.579-2.768). A lack of significant genetic differentiation (Fst = 0.004, p > 0.05) was observed in all the examined fish, indicating that the European flounder in the sampled area constitutes a single genetic cluster. A significant deficiency in heterozygotes (Fis = 0.093, p < 0.05) and overall deviations from Hardy-Weinberg expectations (H-WE) were only detected in fish sampled from the Słupsk Bank. The estimated effective population size (Ne) among the sampled fish groups varied from 712 (Słupsk Bank) to 10,115 (Władysławowo and Mechelinki). However, the recorded values of the Garza-Williamson indicator (M = 0.574-0.600) and the lack of significant (p > 0.05) differences in Heq > He under the SMM model did not support the species' population size changes in the past. The applied SNP markers did not detect the presence of the Baltic flounder among the fish sampled from the studied area. The analysis of an association between biological traits and patterns of genetic diversity did not detect any signs of directional selection or density-dependent adaptive changes in the gene pool of the examined fish that might be caused by increased fishing pressure.

Structural properties and binding mechanism of DNA aptamers sensing saliva melatonin for diagnosis and monitoring of circadian clock and sleep disorders

Circadian desynchrony with the external light-dark cycle influences the rhythmic secretion of melatonin which is among the first signs of circadian rhythm sleep disorders. An accurate dim light melatonin onset (established indicator of circadian rhythm sleep disorders) measurement requires lengthy assays, and antibody affinities alterations, especially in patients with circadian rhythm disorders whose melatonin salivary levels vary significantly, making antibodies detection mostly inadequate. In contrast, aptamers with their numerous advantages (e.g., target selectivity, structural flexibility in tuning binding affinities, small size, etc.) can become preferable biorecognition molecules for salivary melatonin detection with high sensitivity and specificity. This study thoroughly characterizes the structural property and binding mechanism of a single-stranded DNA aptamer full sequence (MLT-C-1) and its truncated versions (MLT-A-2, MLT-A-4) to decipher its optimal characteristics for saliva melatonin detection. We use circular dichroism spectroscopy to determine aptamers' conformational changes under different ionic strengths and showed that aptamers display a hairpin loop structure where few base pairs in the stem play a significant role in melatonin binding and formation of aptamer stabilized structure. Through microscale thermophoresis, aptamers demonstrated a high binding affinity in saliva samples (MLT-C-1F K_d = 12.5 ± 1.7 nM; MLT-A-4F K_d = 11.2 ± 1.6 nM; MLT-A-2F K_d = 2.4 ± 2.8 nM; limit-of-detection achieved in pM, highest sensitivity attained for MLT-A-2F aptamer with the lowest detection limit of 1.35 pM). Our data suggest that aptamers are promising as biorecognition molecules and provide the baseline parameters for the development of an aptamer-based point-of-care diagnostic system for melatonin detection and accurate profiling of its fluctuations in saliva.

Thermodynamic determination of RNA duplex stability in magnesium solutions

The prediction of RNA secondary structure and thermodynamics from sequence relies on free energy minimization and nearest neighbor parameters. Currently, algorithms used to make these predictions are based on parameters from optical melting studies performed in 1 M NaCl. However, many physiological and biochemical buffers containing RNA include much lower concentrations of monovalent cations and the presence of divalent cations. In order to improve these algorithms, thermodynamic data was previously collected for RNA duplexes in solutions containing 71, 121, 221, and 621 mM Na+. From this data, correction factors for free energy (ΔG°37) and melting temperature (Tm) were derived. Despite these newly derived correction factors for sodium, the stabilizing effects of magnesium have been ignored. Here, the same RNA duplexes were melted in solutions containing 0.5, 1.5, 3.0, and 10.0 mM Mg2+ in the absence of monovalent cations. Correction factors for Tm and ΔG°37 were derived to scale the current parameters to a range of magnesium concentrations. The Tm correction factor predicts the melting temperature within 1.2°C, and the ΔG°37 correction factor predicts the free energy within 0.30 kcalmol. These newly derived magnesium correction factors can be incorporated into algorithms that predict RNA secondary structure and stability from sequence.

Nucleic Acid Structure Prediction Including Pseudoknots Through Direct Enumeration of States: A User's Guide to the LandscapeFold Algorithm

Here we detail the LandscapeFold secondary structure prediction algorithm and how it is used. The algorithm was previously described and tested in (Kimchi O et al., Biophys J 117(3):520-532, 2019), though it was not named there. The algorithm directly enumerates all possible secondary structures into which up to two RNA or single-stranded DNA sequences can fold. It uses a polymer physics model to estimate the configurational entropy of structures including complex pseudoknots. We detail each of these steps and ways in which the user can adjust the algorithm as desired. The code is available on the GitHub repository https://github.com/ofer-kimchi/LandscapeFold .

Improving Stability and Specificity of CRISPR/Cas9 System by Selective Modification of Guide RNAs with 2'-fluoro and Locked Nucleic Acid Nucleotides

The genome editing approach using the components of the CRISPR/Cas system has found wide application in molecular biology, fundamental medicine and genetic engineering. A promising method is to increase the efficacy and specificity of CRISPR/Cas-based genome editing systems by modifying their components. Here, we designed and chemically synthesized guide RNAs (crRNA, tracrRNA and sgRNA) containing modified nucleotides (2'-O-methyl, 2'-fluoro, LNA-locked nucleic acid) or deoxyribonucleotides in certain positions. We compared their resistance to nuclease digestion and examined the DNA cleavage efficacy of the CRISPR/Cas9 system guided by these modified guide RNAs. The replacement of ribonucleotides with 2'-fluoro modified or LNA nucleotides increased the lifetime of the crRNAs, while other types of modification did not change their nuclease resistance. Modification of crRNA or tracrRNA preserved the efficacy of the CRISPR/Cas9 system. Otherwise, the CRISPR/Cas9 systems with modified sgRNA showed a remarkable loss of DNA cleavage efficacy. The kinetic constant of DNA cleavage was higher for the system with 2'-fluoro modified crRNA. The 2'-modification of crRNA also decreased the off-target effect upon in vitro dsDNA cleavage.

DeLTa-Seq: direct-lysate targeted RNA-Seq from crude tissue lysate

BACKGROUND

Quantification of gene expression such as RNA-Seq is a popular approach to study various biological phenomena. Despite the development of RNA-Seq library preparation methods and sequencing platforms in the last decade, RNA extraction remains the most laborious and costly step in RNA-Seq of tissue samples of various organisms. Thus, it is still difficult to examine gene expression in thousands of samples.

RESULTS

Here, we developed Direct-RT buffer in which homogenization of tissue samples and direct-lysate reverse transcription can be conducted without RNA purification. The DTT concentration in Direct-RT buffer prevented RNA degradation but not RT in the lysates of several plant tissues, yeast, and zebrafish larvae. Direct reverse transcription on these lysates in Direct-RT buffer produced comparable amounts of cDNA to those synthesized from purified RNA. To maximize the advantage of the Direct-RT buffer, we integrated Direct-RT and targeted RNA-Seq to develop a cost-effective, high-throughput quantification method for the expressions of hundreds of genes: DeLTa-Seq (Direct-Lysate reverse transcription and Targeted RNA-Seq). The DeLTa-Seq method could drastically improve the efficiency and accuracy of gene expression analysis. DeLTa-Seq analysis of 1056 samples revealed the temperature-dependent effects of jasmonic acid and salicylic acid in Arabidopsis thaliana.

CONCLUSIONS

The DeLTa-Seq method can realize large-scale studies using thousands of animal, plant, and microorganism samples, such as chemical screening, field experiments, and studies focusing on individual variability. In addition, Direct-RT is also beneficial for gene expression analysis in small tissues from which it is difficult to purify enough RNA for the experiments.

Self-Assembly of DNA Tiles with G-Quadruplex DNAzyme Catalytic Activity for Sensing Applications

DNA tiles form through self-assembly of a small number of DNA strands that interact through basic repeated interactions, allowing the growth of nanoscale structures seeded by molecular inputs. If an approach for catalytic signal amplification can be integrated into the resultant nanostructure, then one can anticipate biosensing or diagnostic applications mediated by DNA tile self-assembly. Here, two-dimensional DNA tiles with split quadruplexes were designed as diagnostic tools for nucleic acid sensing without the use of protein enzymes. The presence of a target sequence leads to formation of extended microscale structures with arrayed multiple G-quadruplexes across the tile plane, with catalytic activity coupled to a colorimetric reporter. Such a mechanism has potential for low-cost signal amplification using unmodified DNA without the use of protein enzymes for biosensing.

Applicability of the nearest-neighbour model for pseudoknot RNAs

The validity of the nearest-neighbour (NN) model was verified in an RNA pseudoknot (PK) structure. The thermodynamic parameters of the second hairpin stem (S2) region, which separates the PK from a hairpin structure, were monitored using CD and UV melting. Different PKs with identical NN base pairs in the S2 region exhibited similar thermodynamic parameters, highlighting the validity of the NN model in this RNA tertiary structure motif.

Environment-Dependent Stability and Mechanical Properties of DNA Origami Six-Helix Bundles with Different Crossover Spacings

The internal design of DNA nanostructures defines how they behave in different environmental conditions, such as endonuclease-rich or low-Mg²⁺ solutions. Notably, the inter-helical crossovers that form the core of such DNA objects have a major impact on their mechanical properties and stability. Importantly, crossover design can be used to optimize DNA nanostructures for target applications, especially when developing them for biomedical environments. To elucidate this, two otherwise identical DNA origami designs are presented that have a different number of staple crossovers between neighboring helices, spaced at 42- and 21- basepair (bp) intervals, respectively. The behavior of these structures is then compared in various buffer conditions, as well as when they are exposed to enzymatic digestion by DNase I. The results show that an increased number of crossovers significantly improves the nuclease resistance of the DNA origami by making it less accessible to digestion enzymes but simultaneously lowers its stability under Mg²⁺ -free conditions by reducing the malleability of the structures. Therefore, these results represent an important step toward rational, application-specific DNA nanostructure design.

Comparison of rRNA depletion methods for efficient bacterial mRNA sequencing

Current methods of high-throughput RNA sequencing of prokaryotes, including transcriptome analysis or ribosomal profiling, need deep sequencing to achieve sufficient numbers of effective reads (e.g., mapping to mRNA) in order to also find weakly expressed genetic elements. The fraction of high-quality reads mapping to coding RNAs (i.e., mRNA) is mainly influenced by the large content of rRNA and, to a lesser extent, tRNA in total RNA. Thus, depletion of rRNA increases coverage and thus sequencing costs. RiboZero, a depletion kit based on probe hybridisation and rRNA-removal was found to be most efficient in the past, but it was discontinued in 2018. To facilitate comparability with previous experiments and to help choose adequate replacements, we compare three commercially available rRNA depletion kits also based on hybridization and magnetic beads, i.e., riboPOOLs, RiboMinus and MICROBExpress, with the former RiboZero. Additionally, we constructed biotinylated probes for magnetic bead capture and rRNA depletion in this study. Based on E. coli, we found similar efficiencies in rRNA depletion for riboPOOLs and the self-made depletion method; both comparable to the former RiboZero, followed by RiboMinus, succeeded by MICROBExpress. Further, our in-house protocol allows customized species-specific rRNA or even tRNA depletion or depletion of other RNA targets. Both, the self-made biotinylated probes and riboPOOLs, were most successful in reducing the rRNA content and thereby increasing sequencing depth concerning mRNA reads. Additionally, the number of reads matching to weakly expressed genes are increased. In conclusion, the self-made specific biotinylated probes and riboPOOLs are an adequate replacement for the former RiboZero. Both are very efficient in depleting rRNAs, increasing mRNA reads and thus sequencing efficiency.

Star-Polymer-DNA Gels Showing Highly Predictable and Tunable Mechanical Responses

Dynamically crosslinked gels are appealing materials for applications that require time-dependent mechanical responses. DNA duplexes are ideal crosslinkers for building such gels because of their excellent sequence addressability and flexible tunability in bond energy. However, the mechanical responses of most DNA gels are complicated and unpredictable. Here, a DNA gel with a highly homogeneous gel network and well predictable mechanical behaviors is demonstrated by using a pair of star-polymer-DNA precursors with presimulated DNA sequences showing the two-state transition. The melting curve analysis of the DNA gels reveals the good correspondence between the thermodynamic potentials of the DNA crosslinkers and the presimulated values by DNA calculators. Stress-relaxation tests and dissociation kinetics measurements show that the macroscopic relaxation time of the DNA gels is approximately equal to the lifetime of the DNA crosslinkers over 4 orders of magnitude from 0.1-2000 s. Furthermore, a series of durability tests find the DNA gels are hysteresis-less and self-healable after the applications of repeated temperature and mechanical stimuli. These results demonstrate the great potential of star-polymer-DNA precursors for building gels with predictable and tunable viscoelastic properties, suitable for applications such as stress-response extracellular matrices, injectable solids, and soft robotics.

Ion Complexation Explains Orders of Magnitude Changes in the Equilibrium Constant of Biochemical Reactions in Buffers Crowded by Nonionic Compounds

The equilibrium constant (K) of biochemical complex formation in aqueous buffers with high concentration (>20 wt %) of nonionic compounds can vary by orders of magnitude in comparison with the K in a pure buffer. The precise molecular mechanisms of these profound changes are not known. Herein, we show up to a 1000-fold decrease of the K value of DNA hybridization (at nM concentration) in standard molecular crowder systems such as PEG, dextrans, Ficoll, and glycerol. The effect responsible for the decrease of K is the complexation of positively charged ions from a buffer by nonionic polymers/small molecules. We determined the average equilibrium constant for the complexation of ions per monomer (∼0.49 M^-1). We retrieve K's original value for a pure buffer if we properly increase the ionic strength of the buffer crowded by the polymers, compensating for the loss of complexed ions.

Binding of DNA origami to lipids: maximizing yield and switching via strand displacement

Liposomes are widely used as synthetic analogues of cell membranes and for drug delivery. Lipid-binding DNA nanostructures can modify the shape, porosity and reactivity of liposomes, mediated by cholesterol modifications. DNA nanostructures can also be designed to switch conformations by DNA strand displacement. However, the optimal conditions to facilitate stable, high-yield DNA–lipid binding while allowing controlled switching by strand displacement are not known. Here, we characterized the effect of cholesterol arrangement, DNA structure, buffer and lipid composition on DNA–lipid binding and strand displacement. We observed that binding was inhibited below pH 4, and above 200 mM NaCl or 40 mM MgCl₂, was independent of lipid type, and increased with membrane cholesterol content. For simple motifs, binding yield was slightly higher for double-stranded DNA than single-stranded DNA. For larger DNA origami tiles, four to eight cholesterol modifications were optimal, while edge positions and longer spacers increased yield of lipid binding. Strand displacement achieved controlled removal of DNA tiles from membranes, but was inhibited by overhang domains, which are used to prevent cholesterol aggregation. These findings provide design guidelines for integrating strand displacement switching with lipid-binding DNA nanostructures. This paves the way for achieving dynamic control of membrane morphology, enabling broader applications in nanomedicine and biophysics.

Enhancement of the Catalytic Activity of Hammerhead Ribozymes by Organic Cations

Catalytic turnover is important for the application of ribozymes to biotechnology. However, the turnover is often impaired because of the intrinsic high stability of base pairs with cleaved RNA products. Here, organic cations were used as additives to improve the catalytic performance of hammerhead ribozyme constructs that exhibit different kinetic behaviors. Kinetic analysis of substrate cleavage demonstrated that bulky cations, specifically tetra-substituted ammonium ions containing pentyl groups or a benzyl group, have the ability to greatly increase the turnover rate of the ribozymes. Thermal stability analysis of RNA structures revealed that the bulky cations promote the dissociation of cleaved products and refolding of incorrectly folded structures with small disruption of the catalytic structure. The use of bulky cations is a convenient method for enhancing the catalytic activity of hammerhead ribozymes, and the approach may be useful for advancing ribozyme technologies.

Watson-Crick versus Hoogsteen Base Pairs: Chemical Strategy to Encode and Express Genetic Information in Life

Nucleic acids typically form a double helix structure through Watson-Crick base-pairing. In contrast, non-Watson-Crick base pairs can form other three-dimensional structures. Although it is well-known that Watson-Crick base pairs may be more unstable than non-Watson-Crick base pairs under some conditions, the importance of non-Watson-Crick base pairs has not been widely examined. Hoogsteen base pairs, the non-Watson-Crick base pairs, contain important hydrogen-bond patterns that form the helices of nucleic acids, such as in Watson-Crick base pairs, and can form non-double helix structures such as triplexes and quadruplexes. In recent years, non-double helix structures have been discovered in cells and were reported to considerably influence gene expression. The complex behavior of these nucleic acids in cells is gradually being revealed, but the underlying mechanisms remain almost unknown.Quantitatively analyzing the structural stability of nucleic acids is important for understanding their behavior. A nucleic acid is an anionic biopolymer composed of a sugar, base, and phosphoric acid. The physicochemical factors that determine the stability of nucleic acid structures include those derived from the interactions of nucleic acid structures and those derived from the environments surrounding nucleic acids. The Gibbs free energy change (ΔG) of structure formation is the most commonly used physicochemical parameter for analyzing quantitative stability. Quantitatively understanding the intracellular behavior of nucleic acids involves describing the formation of nucleic acid structures and related reactions as ΔG. Based on this concept, we quantitatively analyzed the stability of double helix and non-double helix structures and found that decreased water activity, an important factor in crowded cellular conditions, significantly destabilize the formation of Watson-Crick base pairs but stabilizes Hoogsteen base pairs.Here, we describe a physicochemical approach to understand the regulation of gene expressions based on the stability of nucleic acid structures. We developed new methods for predicting the stability of double and non-double helices in various molecular environments by mimicking intracellular environments. Furthermore, the physicochemical approach used for analyzing gene expression regulated by non-double helix structures is useful for not only determining how gene expression is controlled by cellular environments but also for developing new technologies to chemically regulate gene expression by targeting non-double helix structures. We discuss the roles of Watson-Crick and Hoogsteen base pairs in cells based on our results and why both types of base pairing are required for life. Finally, a new concept in nucleic acid science beyond that of Watson and Crick base pairing is introduced.

DNA hybridisation kinetics using single-molecule fluorescence imaging

Deoxyribonucleic acid (DNA) hybridisation plays a key role in many biological processes and nucleic acid biotechnologies, yet surprisingly there are many aspects about the process which are still unknown. Prior to the invention of single-molecule microscopy, DNA hybridisation experiments were conducted at the ensemble level, and thus it was impossible to directly observe individual hybridisation events and understand fully the kinetics of DNA hybridisation. In this mini-review, recent single-molecule fluorescence-based studies of DNA hybridisation are discussed, particularly for short nucleic acids, to gain more insight into the kinetics of DNA hybridisation. As well as looking at single-molecule studies of intrinsic and extrinsic factors affecting DNA hybridisation kinetics, the influence of the methods used to detect hybridisation of single DNAs is considered. Understanding the kinetics of DNA hybridisation not only gives insight into an important biological process but also allows for further advancements in the growing field of nucleic acid biotechnology.

Salt dependent mesoscopic model for RNA at multiple strand concentrations

Mesoscopic models can be used for the description of the thermodynamic properties of RNA duplexes. With the use of experimental melting temperatures, its parametrization can provide important insights into its hydrogen bonds and stacking interactions as has been done for high sodium concentrations. However, the RNA parametrization for lower salt concentrations is still missing due to the limited amount of published melting temperature data. While the Peyrard-Bishop (PB) parametrization was found to be largely independent of strand concentrations, it requires that all temperatures are provided at the same strand concentrations. Here we adapted the PB model to handle multiple strand concentrations and in this way we were able to make use of an experimental set of temperatures to model the hydrogen bond and stacking interactions at low and intermediate sodium concentrations. For the parametrizations we make a distinction between terminal and internal base pairs, and the resulting potentials were qualitatively similar as we obtained previously for DNA. The main difference from DNA parameters, was the Morse potentials at low sodium concentrations for terminal r(AU) which is stronger than d(AT), suggesting higher hydrogen bond strength.

Combined Effects of Methylated Cytosine and Molecular Crowding on the Thermodynamic Stability of DNA Duplexes

Methylated cytosine within CpG dinucleotides is a key factor for epigenetic gene regulation. It has been revealed that methylated cytosine decreases DNA backbone flexibility and increases the thermal stability of DNA. Although the molecular environment is an important factor for the structure, thermodynamics, and function of biomolecules, there are few reports on the effects of methylated cytosine under a cell-mimicking molecular environment. Here, we systematically investigated the effects of methylated cytosine on the thermodynamics of DNA duplexes under molecular crowding conditions, which is a critical difference between the molecular environment in cells and test tubes. Thermodynamic parameters quantitatively demonstrated that the methylation effect and molecular crowding effect on DNA duplexes are independent and additive, in which the degree of the stabilization is the sum of the methylation effect and molecular crowding effect. Furthermore, the effects of methylation and molecular crowding correlate with the hydration states of DNA duplexes. The stabilization effect of methylation was due to the favorable enthalpic contribution, suggesting that direct interactions of the methyl group with adjacent bases and adjacent methyl groups play a role in determining the flexibility and thermodynamics of DNA duplexes. These results are useful to predict the properties of DNA duplexes with methylation in cell-mimicking conditions.

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.

Prebiotically-relevant low polyion multivalency can improve functionality of membraneless compartments

Multivalent polyions can undergo complex coacervation, producing membraneless compartments that accumulate ribozymes and enhance catalysis, and offering a mechanism for functional prebiotic compartmentalization in the origins of life. Here, we evaluate the impact of lower, more prebiotically-relevant, polyion multivalency on the functional performance of coacervates as compartments. Positively and negatively charged homopeptides with 1–100 residues and adenosine mono-, di-, and triphosphate nucleotides are used as model polyions. Polycation/polyanion pairs are tested for coacervation, and resulting membraneless compartments are analyzed for salt resistance, ability to provide a distinct internal microenvironment (apparent local pH, RNA partitioning), and effect on RNA structure formation. We find that coacervates formed by phase separation of the shorter polyions more effectively generated distinct pH microenvironments, accumulated RNA, and preserved duplexes than those formed by longer polyions. Hence, coacervates formed by reduced multivalency polyions are not only viable as functional compartments for prebiotic chemistries, they can outperform higher molecular weight analogues.

Short cationic peptides and nucleotides can form complex coacervates, but the influence of reduced multivalency on coacervate functionality was not investigated. Here, the authors report that coacervates formed from short polyions generate distinct pH microenvironments, accumulate RNA and preserve nucleic acid duplexes more efficiently than their longer counterparts.

Nearest-neighbor parameters for predicting DNA duplex stability in diverse molecular crowding conditions

The intracellular environment is crowded and heterogeneous. Although the thermodynamic stability of nucleic acid duplexes is predictable in dilute solutions, methods of predicting such stability under specific intracellular conditions are not yet available. We recently showed that the nearest-neighbor model for self-complementary DNA is valid under molecular crowding condition of 40% polyethylene glycol with an average molecular weight of 200 (PEG 200) in 100 mM NaCl. Here, we determined nearest-neighbor parameters for DNA duplex formation under the same crowding condition to predict the thermodynamics of DNA duplexes in the intracellular environment. Preferential hydration of the nucleotides was found to be the key factor for nearest-neighbor parameters in the crowding condition. The determined parameters were shown to predict the thermodynamic parameters (∆H°, ∆S°, and ∆G°₃₇) and melting temperatures (T _m) of the DNA duplexes in the crowding condition with significant accuracy. Moreover, we proposed a general method for predicting the stability of short DNA duplexes in different cosolutes based on the relationship between duplex stability and the water activity of the cosolute solution. The method described herein would be valuable for investigating biological processes that occur under specific intracellular crowded conditions and for the application of DNA-based biotechnologies in crowded environments.

In silico modelling of DNA nanostructures

The rise of material science and nanotechnology created a demand for a next generation of materials and procedures that can transcend the shaping of simple geometrical nano-objects. As a legacy of the technological progress made in the Human Genome Project, DNA was identified as a possible candidate. The low production costs of custom-made DNA molecules and the possibilities concerning the structural manipulation triggered significant advances in the field of DNA nanotechnology in the last decade. To facilitate the development of new DNA nanostructures and provide users an insight in less intuitive complexities and physical properties of the DNA folding, several in silico modelling tools were published. Here, we summarize the main characteristics of these specialised tools, describe the most common design principles and discuss tools and strategies used to predict the properties of DNA nanostructures.

Pioneer factors and their in vitro identification methods

Pioneer transcription factors are a special group of transcription factors that can interact with nucleosomal DNA and initiate regulatory events. Their binding to regulatory regions is the first event in gene activation and can occur in silent or heterochromatin regions. Several research groups have endeavored to define pioneer factors and study their binding characteristics using various techniques. In this review, we describe the in vitro methods used to define and characterize pioneer factors, paying particular attention to differences in methodologies and how these differences can affect results.

Cationic copolymer enhances 8–17 DNAzyme and MNAzyme activities

Cationic copolymer acts as a chaperone to facilitate multiple strand assembly and enhance nucleic acid enzyme activities.

Stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells

This review provides the biophysicochemical background and recent advances in stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells.

Thermodynamic examination of pH and magnesium effect on U6 RNA internal loop

U6 RNA contains a 1 × 2-nt internal loop that folds and unfold during spliceosomal assembly and activation. The 1 × 2 loop consists of a C₆₇•A₇₉ base pair that forms an additional hydrogen bond upon protonation, C₆₇•A⁺₇₉, and uracil (U80) that coordinates the catalytically essential magnesium ions. We designed a series of RNA and DNA constructs with a 1 × 2 loop sequence contained in the ISL, and its modifications, to measure the thermodynamic effects of protonation and magnesium binding using UV-visible thermal denaturation experiments. We show that the wild-type RNA construct gains 0.43 kcal/mol in 1 M KCl upon lowering the pH from 7.5 to 5.5; the presence of magnesium ions increases its stability by 2.17 kcal/mol at pH 7.5 over 1 M KCl. Modifications of the helix closing base pairs from C-G to U•G causes a loss in protonation-dependent stability and a decrease in stability in the presence of magnesium ions, especially in the C68U construct. A79G single-nucleotide bulge loop construct showed the largest gain in stability in the presence of magnesium ions. The DNA wild-type construct shows a smaller effect on stability upon lowering the pH and in the presence of magnesium ions, highlighting differences in RNA and DNA structures. A U6 RNA 1 × 2 loop sequence is rare in the databases examined.

Validation of the nearest-neighbor model for Watson-Crick self-complementary DNA duplexes in molecular crowding condition

Recent advancement in nucleic acid techniques inside cells demands the knowledge of the stability of nucleic acid structures in molecular crowding. The nearest-neighbor model has been successfully used to predict thermodynamic parameters for the formation of nucleic acid duplexes, with significant accuracy in a dilute solution. However, knowledge about the applicability of the model in molecular crowding is still limited. To determine and predict the stabilities of DNA duplexes in a cell-like crowded environment, we systematically investigated the validity of the nearest-neighbor model for Watson–Crick self-complementary DNA duplexes in molecular crowding. The thermodynamic parameters for the duplex formation were measured in the presence of 40 wt% poly(ethylene glycol)200 for different self-complementary DNA oligonucleotides consisting of identical nearest-neighbors in a physiological buffer containing 0.1 M NaCl. The thermodynamic parameters as well as the melting temperatures (T_m) obtained from the UV melting studies revealed similar values for the oligonucleotides having identical nearest-neighbors, suggesting the validity of the nearest-neighbor model in the crowding condition. Linear relationships between the measured ΔG°₃₇ and T_m in crowding condition and those predicted in dilute solutions allowed us to predict ΔG°₃₇, T_m and nearest-neighbor parameters in molecular crowding using existing parameters in the dilute condition, which provides useful information about the thermostability of the self-complementary DNA duplexes in molecular crowding.

Stabilization of DNA Loop Structures by Large Cations

The DNA-binding properties of large cations differ from those of metal ions due to steric exclusion from base-paired regions. In this study, the thermal stability of DNA secondary structures, including duplexes, internal loops, bulge loops, hairpin loops, dangling ends, and G-quadruplexes, was investigated in the presence of cations of different sizes. Large cations, such as tetrabutylammonium and tetrapentylammonium ions, reduced the stability of fully matched duplexes but increased the stability of duplexes with a long loop. The cations also increased the stability of G-quadruplexes with a long loop, and the degree of stabilization was greater for low-stability G-quadruplexes. Analysis of the salt concentration dependence indicates that large cations bind to the loop nucleotides, leading to counteracting the destabilization effect on base pairing. It is likely that binding occurs when loop nucleotides are sufficiently flexible to allow for greater accessibility for large cations. These results provide insight into nucleic acid interactions with large cationic molecules and suggest a potential method for stabilizing noncanonical DNA structures under intracellular conditions.

Nearest-neighbour parameters optimized for melting temperature prediction of DNA/RNA hybrids at high and low salt concentrations

Gene editing technologies sparked a renewed interest in the hybridization of DNA/RNA duplexes, yet little improvement on nearest-neighbour parameters was made over the past two decades. For low sodium concentration no parameter set was yet calculated. Here, we revised the existing experimental datasets and used an expanded set of sequences from which we recalculated the nearest-neighbour parameters, reducing the average temperature prediction uncertainty to 1.6 °C. Two experimental sets using temperatures extracted via different methods were used with similar results, with the curve-fitting method achieving a slight advantage in prediction quality over other methods. Additionally, we obtained new parameters for low salt with an average uncertainty of 0.98 °C. We also tested several types of salt correction factors and concluded that it is advisable to use those originally developed for RNA/RNA rather than for DNA/DNA.

Kinetics and equilibrium constants of oligonucleotides at low concentrations. Hybridization and melting study

Although DNA hybridization/melting is one of the most important biochemical reactions, the non-trivial kinetics of the process is not yet fully understood. In this work, we use Förster resonance energy transfer (FRET) to investigate the influence of temperature, ionic strength, and oligonucleotide length on the kinetic and equilibrium constants of DNA oligonucleotide binding and dissociation. We show that at low reagent concentrations and ionic strength, the time needed to establish equilibrium between single and double strand forms may be of the order of days, even for simple oligonucleotides of a length of 20 base pairs or less. We also identify and discuss the possible artifacts related to fluorescence-based experiments conducted in extremely dilute solutions. The results should prove useful for the judicious design of technologies based on DNA-matching, including sensors, DNA multiplication, sequencing, and gene manipulation.

OligoCOOL: A mobile application for nucleotide sequence analysis

Today smartphones are inseparable parts of modern life and are capable of performing many desktop computers' tasks such as scientific analysis with greater convenience. Here, we present OligoCOOL, which is an Android application for analyzing nucleic sequences. This application enables users to perform several common biomedical analyses for a given nucleotide sequence. OligoCOOL is a freely accessible Android app at http://bioinf.modares.ac.ir/software/OligoCOOL, which can be a suitable tool for the experimental design in the laboratories. This application also can be used to learn the basics of nucleotide sequence analysis. © 2019 International Union of Biochemistry and Molecular Biology, 47(2): 201-206, 2019.

Colloid supported lipid bilayers for self-assembly

The use of colloid supported lipid bilayers (CSLBs) has recently been extended to create colloidal joints, that enable the assembly of structures with internal degrees of flexibility, and to study lipid membranes on curved and closed geometries. These novel applications of CSLBs rely on previously unappreciated properties: the simultaneous fluidity of the bilayer, lateral mobility of inserted (linker) molecules and colloidal stability. Here we characterize every step in the manufacturing of CSLBs in view of these requirements using confocal microscopy and fluorescence recovery after photobleaching (FRAP). Specifically, we have studied the influence of different particle properties (roughness, surface charge, chemical composition, polymer coating) on the quality and mobility of the supported bilayer. We find that the insertion of lipopolymers in the bilayer can affect its homogeneity and fluidity. We improve the colloidal stability by inserting lipopolymers or double-stranded inert DNA into the bilayer. We include surface-mobile DNA linkers and use FRAP to characterize their lateral mobility both in their freely diffusive and bonded state. Finally, we demonstrate the self-assembly of flexibly linked structures from the CSLBs modified with surface-mobile DNA linkers. Our work offers a collection of experimental tools for working with CSLBs in applications ranging from controlled bottom-up self-assembly to model membrane studies.

Replacing salt correction factors with optimized RNA nearest-neighbour enthalpy and entropy parameters

We calculate the nearest-neighbour enthalpies and entropies at 5 salt concentrations of 18 RNA sequences, each for at least 9 different species concentrations, totalling 757 melting temperatures, using a melting temperature optimization method. These new parameters do not need to be salt-corrected and are shown to provide overall improved melting temperature predictions. They show a marked quadratic dependence with salt concentrations which are compensated to form linear Gibbs free energies. Two different parameter schemes were tested, with fixed or variable initial parameters. We have found that using variable initial parameters provides better predictive results than using salt correction factors and that the prediction uncertainty is considerably reduced for a validation set of independent sequences. An interpolation scheme is introduced to generate model parameters for arbitrary salt concentrations which performs better against a validation set than predictions using salt corrections.

High-performance interactive analysis of split aptamer and HIV-1 Tat on multiwall carbon nanotube-modified field-effect transistor

Interaction between split RNA aptamer and the clinically important target, HIV-1 Tat was investigated on a biosensing surface transduced by functionally choreographed multiwall carbon nanotubes (MWCNTs). Acid oxidation was performed to functionalize MWCNTs with carboxyl functional groups. X-ray photoelectron spectroscopy analysis had profound ~2.91% increment in overall oxygen group and ~1% increment was noticed with a specific carboxyl content owing to CO and OCO bonding. The interaction between split RNA aptamer and HIV-1 Tat protein was quantified by electrical measurements with the current signal (I_ds) over a gate voltage (V_gs). Initially, 34.4 mV gate voltage shift was observed by the immobilization of aptamer on MWCNT. With aptamer and HIV-1 Tat interaction, the current flow was decreased with the concomitant gate voltage shift of 23.5 mV. The attainment of sensitivity with split aptamer and HIV-1 Tat interaction on the fabricated device was 600 pM. To ensure the genuine interaction of aptamer with HIV-1 Tat, other HIV-1 proteins, Nef and p24 were interacted with aptamer and they displayed the negligible interferences with gate voltage shift of 3.5 mV and 5.7 mV, which shows 4 and 2.5 folds lesser than HIV-1 Tat interaction, respectively.