Effects of Sodium Ions on DNA Duplex Oligomers: Improved Predictions of Melting Temperatures

Rapid Systematic Screening of Bispecific Antibody Surrogate Geometries for T-Cell Engagement Using DNA Nanotechnology

Bispecific antibodies (bsAbs) are emerging immune-therapeutics, and many formats exist that differ considerably in structure. However, little systematic data exist about how the spatial organization of their components influences activity, requiring innovative approaches combining empirical and quantitative frameworks. This study presents a modular DNA nanotechnology platform to generate numerous bsAbs with surrogate geometries that span the structural features of the BiTE, IgG-like, and IgG-conjugate platforms to screen for T-cell engagement. Results highlight interesting structure-activity relationships regarding bsAb potency and selectivity and raise questions regarding the molecular phenomena underlying activity. To elucidate some effects, the platform was paired with a simple mathematical model. This work is thus one of the first to systematically investigate and reveal the importance of the spatial organization of bsAb components on activity and equally provides an accessible and convenient tool for rapidly mapping out such trends for other combinations of target epitopes.

Efficient enrichment of free target sequences in an integrated microfluidic device for point-of-care detection systems

Nucleic acid biomarker detection has great importance in the diagnosis of disease, the monitoring of disease progression and the classification of patients according to treatment decision making. Nucleic acid biomarkers found in the blood of patients have generated a lot of interest due to the possibility of being detected non-invasively which makes them ideal for monitoring and screening tests and particularly amenable to point-of-care (POC) or self-testing. A major challenge to POC molecular diagnostics is the need to enrich the target to optimise detection. In this work, we describe a microfabricated device for the enrichment of short dsDNA target sequences, which is especially valuable for potential detection methods, as it improves the probability of effectively detecting the target in downstream analyses. The device integrated a heating element and a temperature sensor with a microfluidic chamber to carry out the denaturation of the dsDNA combined with blocking-probes to enrich the target. This procedure was validated by fluorescence resonance energy transfer (FRET) technique, labelling DNA with a fluorophore and a quencher. As proof of concept, a 23-mer long dsDNA sequence corresponding to the L858R mutation of the EGFR gene was used. The qualitative results obtained determined that the most optimal blocking rate was obtained with the incorporation of 11/12-mer blocking-probes at a total concentration of 6 μM. This device is a powerful DNA preparation tool, which is an indispensable initial step for subsequent detection of sequences via nucleic acid hybridisation methods.

Impact of mobile and stationary phases on siRNA duplex stability in liquid chromatography

Nucleic acid duplexes are typically analyzed in non-denaturing conditions. Melting temperature (Tm) is the property used to measure duplex stability; however, it is not known how the chromatographic conditions and mobile phase composition affect the duplex stability. We employed differential scanning calorimetry (DSC) method to measure the melting temperature of chemically modified silencing RNA duplex (21 base pairs, 0.15 mM duplex concentration) in mobile phases commonly used in reversed-phase, ion-pair reversed-phase, size exclusion and hydrophilic interaction chromatography. We investigated mobile phases consisting of ammonium acetate, alkylammonium ion-pairing reagents, alkali-ion chlorides, magnesium chloride, and additives including methanol, ethanol, acetonitrile and hexafluoroisopropanol. Increasing buffer concentration enhanced the duplex stability (Tm was 67.1 - 78.2 °C for 10-100 mM [Na+] concentration). The melting temperature decreases with the increase in cation size (70.2 °C in 10 mM [Li+], 68.1 °C in 10 mM [NH4+], 65.6 °C in 10 mM [Cs+], and 56.6 °C in 10 mM [triethylammonium+] solutions). Inclusion of 20 % of organic solvent in buffer reduced the melting temperature by 1-3 °C, and denaturation power increases in the order MeOH Collapse

Key Words

• Denaturing conditions
• Differential scanning calorimetry
• Liquid chromatography
• Melting temperature
• siRNA duplex

Collapse

MESH Headings

• RNA, Small Interfering/chemistry
• Calorimetry, Differential Scanning
• RNA Stability
• Chromatography, Reverse-Phase/methods
• Acetonitriles/chemistry
• Acetates/chemistry
• Methanol/chemistry
• Hydrophobic and Hydrophilic Interactions
• Solvents/chemistry
• Propanols/chemistry
• Chromatography, Liquid/methods
• Ethanol/chemistry
• Transition Temperature
• Chromatography, Gel/methods
• Magnesium Chloride/chemistry
• Hydrocarbons, Fluorinated

Collapse

Grants Collapse

Affiliation(s)

Martin Gilar Waters Corporation, 34 Maple Street, Milford, MA, 01757, USA.
Samuel Redstone Waters Corporation, 34 Maple Street, Milford, MA, 01757, USA
Alexandre Gomes Waters Corporation, 34 Maple Street, Milford, MA, 01757, USA

Collapse

Controlling the thermally-driven crystallization of DNA-coated nanoparticles with formamide

DNA-coated nanoparticles, also known as programmable atom equivalents (PAEs), facilitate the construction of materials with nanoscopic precision. Thermal annealing plays a pivotal role by controlling DNA hybridization kinetics and thermodynamics, which ensures the formation of intended structures. While various design handles such as particle size, DNA design, and salt concentration influence the stability of the DNA duplexes linking PAEs in a lattice, their influence on the system's melting temperature (Tm) often follows complicated trends that make rational tuning of self-assembly challenging. In this work, the denaturant formamide is used to precisely tune the thermal response of PAEs. Our results reveal a clear and predictable trend in the PAEs' response to formamide, enabling rational control over the Tm of a diverse set of PAE systems. Unlike adjustments made through alterations to PAE design or solution parameters such as ionic strength, formamide achieves its temperature shift without impacting the kinetics of assembly. As a result, PAEs can be rapidly crystallized at ambient temperatures, producing superlattices with similar quality to PAE crystals assembled through standard protocols that use higher temperatures. This study therefore positions formamide as a useful tool for enhancing the synthesis of complex nanostructures under mild conditions.

A pipeline for the development and analysis of extracellular vesicle-based transcriptomic biomarkers in molecular diagnostics

Extracellular vesicles are shed by every cell type and can be found in any biofluid. They contain different molecules that can be utilized as biomarkers, including several RNA species which they protect from degradation. Here, we present a pipeline for the development and analysis of extracellular vesicle-associated transcriptomic biomarkers that our group has successfully applied multiple times. We highlight the key steps of the pipeline and give particular emphasis to the necessary quality control checkpoints, which are linked to numerous available guidelines that should be considered along the workflow. Our pipeline starts with patient recruitment and continues with blood sampling and processing. The purification and characterization of extracellular vesicles is explained in detail, as well as the isolation and quality control of extracellular vesicle-associated RNA. We point out the possible pitfalls during library preparation and RNA sequencing and present multiple bioinformatic tools to pinpoint biomarker signature candidates from the sequencing data. Finally, considerations and pitfalls during the validation of the biomarker signature using RT-qPCR will be elaborated.

Nearest-neighbour parametrization of DNA single, double and triple mismatches at low sodium concentration

DNA mismatches, that is, base pairs different from the canonical AT and CG, are involved in numerous biological processes and can be a problem for technological applications such as PCR amplification. The nearest-neighbour (NN) model is the standard approach for predicting melting temperatures and is used in methods of secondary structure predictions and modelling of hybridization kinetics. However, despite its biological and technological importance, existing NN parameters that include DNA mismatches are incomplete, and those available were obtained from a limited set of melting temperature at high sodium concentration. To our knowledge, there is currently no NN set of parameters for up to three mismatches covering all configurations at low sodium concentrations. Here, we are applying the NN model to a large set of 4096 published melting temperatures, covering all combinations of single, double and triple mismatches. Dealing with such a large set of temperature is challenging in several ways, bringing new methodological problems. Here, optimizing a large number of 252 independent parameters has required the development of a new method where we readjust the seed parameters using the definition of the Gibbs free energy. The new parameters predict the training set within 1.1 °C and the validation set to 2.7 °C.

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.

Unusual similarity of DNA solvation dynamics in high-salinity crowding with divalent cations of varying concentrations

Double-stranded DNA bears the highest linear negative charge density (2e- per base-pair) among all biopolymers, leading to strong interactions with cations and dipolar water, resulting in the formation of a dense 'condensation layer' around DNA. Interactions involving proteins and ligands binding to DNA are primarily governed by strong electrostatic forces. Increased salt concentrations impede such electrostatic interactions - a situation that prevails in oceanic species due to their cytoplasm being enriched with salts. Nevertheless, how these interactions' dynamics are affected in crowded hypersaline environments remains largely unexplored. Here, we employ steady-state and time-resolved fluorescence Stokes shifts (TRFSS) of a DNA-bound ligand (DAPI) to investigate the static and dynamic solvation properties of DNA in the presence of two divalent cations, magnesium (Mg2+), and calcium (Ca2+) at varying high to very-high concentrations of 0.15 M, 1 M and 2 M. We compare the results to those obtained in physiological concentrations (0.15 M) of monovalent Na+ ions. Combining data from fluorescence femtosecond optical gating (FOG) and time-correlated single photon counting (TCSPC) techniques, dynamic fluorescence Stokes shifts in DNA are analysed over a broad range of time-scales, from 100 fs to 10 ns. We find that while divalent cation crowding strongly influences the DNA stability and ligand binding affinity to DNA, the dynamics of DNA solvation remain remarkably similar across a broad range of five decades in time, even in a high-salinity crowded environment with divalent cations, as compared to the physiological concentration of the Na+ ion. Steady-state and time-resolved data of the DNA-groove-bound ligand are seemingly unaffected by ion-crowding in hypersaline solution, possibly due to ions being mostly displaced by the DNA-bound ligand. Furthermore, the dynamic coupling of cations with nearby water may possibly contribute to a net-neutral effect on the overall collective solvation dynamics in DNA, owing to the strong anti-correlation of their electrostatic interaction energy fluctuations. Such dynamic scenarios may persist within the cellular environment of marine life and other biological cells that experience hypersaline conditions.

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.

Mono-valent salt corrections for RNA secondary structures in the ViennaRNA package

BACKGROUND

RNA features a highly negatively charged phosphate backbone that attracts a cloud of counter-ions that reduce the electrostatic repulsion in a concentration dependent manner. Ion concentrations thus have a large influence on folding and stability of RNA structures. Despite their well-documented effects, salt effects are not handled consistently by currently available secondary structure prediction algorithms. Combining Debye-Hückel potentials for line charges and Manning's counter-ion condensation theory, Einert et al. (Biophys J 100: 2745-2753, 2011) modeled the energetic contributions of monovalent cations on loops and helices.

RESULTS

The model of Einert et al. is adapted to match the structure of the dynamic programming recursion of RNA secondary structure prediction algorithms. An empirical term describing the salt dependence of the duplex initiation energy is added to improve co-folding predictions for two or more RNA strands. The slightly modified model is implemented in the ViennaRNA package in such way that only the energy parameters but not the algorithmic structure is affected. A comparison with data from the literature show that predicted free energies and melting temperatures are in reasonable agreement with experiments.

CONCLUSION

The new feature in the ViennaRNA package makes it possible to study effects of salt concentrations on RNA folding in a systematic manner. Strictly speaking, the model pertains only to mono-valent cations, and thus covers the most important parameter, i.e., the NaCl concentration. It remains a question for future research to what extent unspecific effects of bi- and tri-valent cations can be approximated in a similar manner.

AVAILABILITY

Corrections for the concentration of monovalent cations are available in the ViennaRNA package starting from version 2.6.0.

Performance analysis of DNA crossbar arrays for high-density memory storage applications

Deoxyribonucleic acid (DNA) has emerged as a promising building block for next-generation ultra-high density storage devices. Although DNA has high durability and extremely high density in nature, its potential as the basis of storage devices is currently hindered by limitations such as expensive and complex fabrication processes and time-consuming read-write operations. In this article, we propose the use of a DNA crossbar array architecture for an electrically readable read-only memory (DNA-ROM). While information can be 'written' error-free to a DNA-ROM array using appropriate sequence encodings its read accuracy can be affected by several factors such as array size, interconnect resistance, and Fermi energy deviations from HOMO levels of DNA strands employed in the crossbar. We study the impact of array size and interconnect resistance on the bit error rate of a DNA-ROM array through extensive Monte Carlo simulations. We have also analyzed the performance of our proposed DNA crossbar array for an image storage application, as a function of array size and interconnect resistance. While we expect that future advances in bioengineering and materials science will address some of the fabrication challenges associated with DNA crossbar arrays, we believe that the comprehensive body of results we present in this paper establishes the technical viability of DNA crossbar arrays as low power, high-density storage devices. Finally, our analysis of array performance vis-à-vis interconnect resistance should provide valuable insights into aspects of the fabrication process such as proper choice of interconnects necessary for ensuring high read accuracies.

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.

Structure and Dynamics of dsDNA in Cell-like Environments

Deoxyribonucleic acid (DNA) is a fundamental biomolecule for correct cellular functioning and regulation of biological processes. DNA's structure is dynamic and has the ability to adopt a variety of structural conformations in addition to its most widely known double-stranded DNA (dsDNA) helix structure. Stability and structural dynamics of dsDNA play an important role in molecular biology. In vivo, DNA molecules are folded in a tightly confined space, such as a cell chamber or a channel, and are highly dense in solution; their conformational properties are restricted, which affects their thermodynamics and mechanical properties. There are also many technical medical purposes for which DNA is placed in a confined space, such as gene therapy, DNA encapsulation, DNA mapping, etc. Physiological conditions and the nature of confined spaces have a significant influence on the opening or denaturation of DNA base pairs. In this review, we summarize the progress of research on the stability and dynamics of dsDNA in cell-like environments and discuss current challenges and future directions. We include studies on various thermal and mechanical properties of dsDNA in ionic solutions, molecular crowded environments, and confined spaces. By providing a better understanding of melting and unzipping of dsDNA in different environments, this review provides valuable guidelines for predicting DNA thermodynamic quantities and for designing DNA/RNA nanostructures.

Ab initio predictions for 3D structure and stability of single- and double-stranded DNAs in ion solutions

The three-dimensional (3D) structure and stability of DNA are essential to understand/control their biological functions and aid the development of novel materials. In this work, we present a coarse-grained (CG) model for DNA based on the RNA CG model proposed by us, to predict 3D structures and stability for both dsDNA and ssDNA from the sequence. Combined with a Monte Carlo simulated annealing algorithm and CG force fields involving the sequence-dependent base-pairing/stacking interactions and an implicit electrostatic potential, the present model successfully folds 20 dsDNAs (≤52nt) and 20 ssDNAs (≤74nt) into the corresponding native-like structures just from their sequences, with an overall mean RMSD of 3.4Å from the experimental structures. For DNAs with various lengths and sequences, the present model can make reliable predictions on stability, e.g., for 27 dsDNAs with/without bulge/internal loops and 24 ssDNAs including pseudoknot, the mean deviation of predicted melting temperatures from the corresponding experimental data is only ~2.0°C. Furthermore, the model also quantificationally predicts the effects of monovalent or divalent ions on the structure stability of ssDNAs/dsDNAs.

To determine 3D structures and quantify stability of single- (ss) and double-stranded (ds) DNAs is essential to unveil the mechanisms of their functions and to further guide the production and development of novel materials. Although many DNA models have been proposed to reproduce the basic structural, mechanical, or thermodynamic properties of dsDNAs based on the secondary structure information or preset constraints, there are very few models can be used to investigate the ssDNA folding or dsDNA assembly from the sequence. Furthermore, due to the polyanionic nature of DNAs, metal ions (e.g., Na⁺ and Mg²⁺) in solutions can play an essential role in DNA folding and dynamics. Nevertheless, ab initio predictions for DNA folding in ion solutions are still an unresolved problem. In this work, we developed a novel coarse-grained model to predict 3D structures and thermodynamic stabilities for both ssDNAs and dsDNAs in monovalent/divalent ion solutions from their sequences. As compared with the extensive experimental data and available existing models, we showed that the present model can successfully fold simple DNAs into their native-like structures, and can also accurately reproduce the effects of sequence and monovalent/divalent ions on structure stability for ssDNAs including pseudoknot and dsDNAs with/without bulge/internal loops.

Non-linear Hamiltonian models for DNA

Nucleic acids' physical properties have been investigated by theoretical methods based both on fully atomistic representations and on coarse-grained models, e.g., the worm-like-chain, taken from polymer physics. In this review article, I discuss an intermediate (mesoscopic) approach and show how to build a three-dimensional Hamiltonian model which accounts for the main interactions responsible for the stability of the helical molecules. While the 3D mesoscopic model yields a sufficiently detailed description of the helix at the level of the base pair, it also allows one to predict the thermodynamical and structural properties of molecules in solution. Relying on the idea that the base pair fluctuations can be conceived as trajectories, I have built over the past years a computational method based on the time-dependent path integral formalism to derive the partition function. While the main features of the method are presented, I focus here in particular on a newly developed statistical method to set the maximum amplitude of the base pair fluctuations, a key parameter of the theory. Some applications to the calculation of DNA flexibility properties are discussed together with the available experimental data.

Computationally exploring the mechanism of bacteriophage T7 gp4 helicase translocating along ssDNA

Bacteriophage T7 gp4 helicase has served as a model system for understanding mechanisms of hexameric replicative helicase translocation. The mechanistic basis of how nucleoside 5'-triphosphate hydrolysis and translocation of gp4 helicase are coupled is not fully resolved. Here, we used a thermodynamically benchmarked coarse-grained protein force field, Associative memory, Water mediated, Structure and Energy Model (AWSEM), with the single-stranded DNA (ssDNA) force field 3SPN.2C to investigate gp4 translocation. We found that the adenosine 5'-triphosphate (ATP) at the subunit interface stabilizes the subunit-subunit interaction and inhibits subunit translocation. Hydrolysis of ATP to adenosine 5'-diphosphate enables the translocation of one subunit, and new ATP binding at the new subunit interface finalizes the subunit translocation. The LoopD2 and the N-terminal primase domain provide transient protein-protein and protein-DNA interactions that facilitate the large-scale subunit movement. The simulations of gp4 helicase both validate our coarse-grained protein-ssDNA force field and elucidate the molecular basis of replicative helicase translocation.

Enzymatic Synthesis of Chemical Nuclease Triplex-Forming Oligonucleotides with Gene-Silencing Applications

Triplex-forming oligonucleotides (TFOs) are short, single-stranded oligomers that hybridise to a specific sequence of duplex DNA. TFOs can block transcription and thereby inhibit protein production, making them highly appealing in the field of antigene therapeutics. In this work, a primer extension protocol was developed to enzymatically prepare chemical nuclease TFO hybrid constructs, with gene-silencing applications. Click chemistry was employed to generate novel artificial metallo-nuclease (AMN)-dNTPs, which were selectively incorporated into the TFO strand by a DNA polymerase. This purely enzymatic protocol was then extended to facilitate the construction of 5-methylcytosine (5mC) modified TFOs that displayed increased thermal stability. The utility of the enzymatically synthesised di-(2-picolyl)amine (DPA)-TFOs was assessed and compared to a specifically prepared solid-phase synthesis counterpart through gel electrophoresis, quantitative PCR, and Sanger sequencing, which revealed similar recognition and damage properties to target genes. The specificity was then enhanced through coordinated designer intercalators-DPQ and DPPZ-and high-precision DNA cleavage was achieved. To our knowledge, this is the first example of the enzymatic production of an AMN-TFO hybrid and is the largest base modification incorporated using this method. These results indicate how chemical nuclease-TFOs may overcome limitations associated with non-molecularly targeted metallodrugs and open new avenues for artificial gene-editing technology.

DNA-Crosslinked 2D Photonic Crystal Hydrogels for Detection of Adenosine Actuated by an Adenosine-Binding Aptamer

There is a need to develop versatile sensing motifs that can be used to detect a variety of chemical targets in resource-limited settings, for example, at the point of care. While numerous sensing technologies have been developed toward this effort, these technologies can be overly complex and require a skilled technician, extensive sample preparation, or sophisticated instrumentation to use, limiting their generalizability and application in resource-limited settings. Here, we report a novel sensing motif that utilizes DNA-crosslinked two-dimensional photonic crystal (2DPC) hydrogels. These hydrogel sensors contain a DNA aptamer recognition group that binds a target analyte. As proof of concept, we fabricated 2DPC hydrogels using a well-studied adenosine-binding aptamer. This adenosine aptamer is duplexed with a partially complementary strand and forms responsive crosslinks in the hydrogel polymer network. When adenosine is introduced, aptamer-adenosine binding occurs, breaking the DNA crosslinks and causing the hydrogel to swell. This in turn increases the particle spacing of an embedded 2DPC array, shifting the 2DPC Bragg diffraction. Thus, adenosine concentration can be monitored through 2DPC Bragg diffraction measurements. A linear range of 20 μM to 2 mM was observed. The detection limits were calculated to be 13.9 μM in adenosine-binding buffer and 26.7 μM in fetal bovine serum. This reported sensing motif has a readout that is simple and rapid and requires minimal equipment. We hypothesize that this sensing motif is generalizable and that other sensors can be easily fabricated by simply exchanging the aptamer that serves as a molecular recognition group.

Forming a Double-Helix Phase of Single Polymer Chains by the Cooperation between Local Structure and Nonlocal Attraction

Double-helix structures, such as DNA, are formed in nature to realize many unique functions. Inspired by this, researchers are pursuing strategies to design such structures from polymers. A key question is whether the double helix can be formed from the self-folding of a single polymer chain without specific interactions. Here, using Langevin dynamics simulation and theoretical analysis, we find that a stable double-helix phase can be achieved by the self-folding of single semiflexible polymers as a result of the cooperation between local structure and nonlocal attraction. The critical temperature of double-helix formation approximately follows T^{cri}∼ln(k_{θ}) and T^{cri}∼ln(k_{τ}), where k_{θ} and k_{τ} are the polymer bending and torsion stiffness, respectively. Furthermore, the double helix can exhibit major and minor grooves due to symmetric break for better packing. Our results provide a novel guide to the experimental design of the double helix.

High-Force Application by a Nanoscale DNA Force Spectrometer

The ability to apply and measure high forces (>10 pN) on the nanometer scale is critical to the development of nanomedicine, molecular robotics, and the understanding of biological processes such as chromatin condensation, membrane deformation, and viral packaging. Established force spectroscopy techniques including optical traps, magnetic tweezers, and atomic force microscopy rely on micron-sized or larger handles to apply forces, limiting their applications within constrained geometries including cellular environments and nanofluidic devices. A promising alternative to these approaches is DNA-based molecular calipers. However, this approach is currently limited to forces on the scale of a few piconewtons. To study the force application capabilities of DNA devices, we implemented DNA origami nanocalipers with tunable mechanical properties in a geometry that allows application of force to rupture a DNA duplex. We integrated static and dynamic single-molecule characterization methods and statistical mechanical modeling to quantify the device properties including force output and dynamic range. We found that the thermally driven dynamics of the device are capable of applying forces of at least 20 piconewtons with a nanometer-scale dynamic range. These characteristics could eventually be used to study other biomolecular processes such as protein unfolding or to control high-affinity interactions in nanomechanical devices or molecular robots.

Pairing statistics and melting of random DNA oligomers: Finding your partner in superdiverse environments

Understanding of the pairing statistics in solutions populated by a large number of distinct solute species with mutual interactions is a challenging topic, relevant in modeling the complexity of real biological systems. Here we describe, both experimentally and theoretically, the formation of duplexes in a solution of random-sequence DNA (rsDNA) oligomers of length L = 8, 12, 20 nucleotides. rsDNA solutions are formed by 4^L distinct molecular species, leading to a variety of pairing motifs that depend on sequence complementarity and range from strongly bound, fully paired defectless helices to weakly interacting mismatched duplexes. Experiments and theory coherently combine revealing a hybridization statistics characterized by a prevalence of partially defected duplexes, with a distribution of type and number of pairing errors that depends on temperature. We find that despite the enormous multitude of inter-strand interactions, defectless duplexes are formed, involving a fraction up to 15% of the rsDNA chains at the lowest temperatures. Experiments and theory are limited here to equilibrium conditions.

Several biological processes require that specific partner molecules succeed in binding after negotiating their way through a huge number of interactions with other molecules. How such molecular recognition emerges among millions distinct molecular species is an open problem. We have studied, both experimentally and theoretically, such process of “molecular recognition” in pools of highly diverse random DNA oligomers, which binds preferentially, but not exclusively, to its perfect complementary sequence. We find a complex behavior, in which some perfect pairing takes place with a non-trivial temperature dependence that we understand thorough statistical mechanics modelling. The pairing pattern of short random DNA is relevant in the context of the origin of life since the so-called “RNA World” was most probably based on the mutual recognition of random chains.

Needles in Haystacks: Understanding the Success of Selective Pairing of Nucleic Acids

The effectiveness of several biological and biotechnological processes relies on the remarkably selective pairing of nucleic acids in contexts of molecular complexity. Relevant examples are the on-target binding of primers in genomic PCR and the regulatory efficacy of microRNA via binding on the transcriptome. Here, we propose a statistical framework that enables us to describe and understand such selectivity by means of a model that is extremely cheap from a computational point of view. By re-parametrizing the hybridization thermodynamics on three classes of base pairing errors, we find a convenient way to obtain the free energy of pairwise interactions between nucleic acids. We thus evaluate the hybridization statistics of a given oligonucleotide within a large number of competitive sites that we assume to be random, and we compute the probability of on-target binding. We apply our strategy to PCR amplification and microRNA-based gene regulation, shedding new light on their selectivity. In particular, we show the relevance of the defectless pairing of 3′ terminals imposed by the polymerase in PCR selection. We also evaluate the selectivity afforded by the microRNA seed region, thus quantifying the extra contributions given by mechanisms beyond pairing statistics.

Characterization of a Programmable Argonaute Nuclease from the Mesophilic Bacterium Rummeliibacillus suwonensis

Prokaryotic Argonautes (pAgos) from mesophilic bacteria are attracting increasing attention for their genome editing potential. So far, it has been reported that KmAgo from Kurthia massiliensis can utilize DNA and RNA guide of any sequence to effectively cleave DNA and RNA targets. Here we find that three active pAgos, which have about 50% sequence identity with KmAgo, possess typical DNA-guided DNA target cleavage ability. Among them, RsuAgo from Rummeliibacillus suwonensis is mainly explored for which can cleave both DNA and RNA targets. Interestingly, RsuAgo-mediated RNA target cleavage occurs only with short guide DNAs in a narrow length range (16–20 nt), and mismatches between the guide and target sequence greatly affect the efficiency of RNA target cleavage. RsuAgo-mediated target cleavage shows a preference for a guide strand with a 5′-terminal A residue. Furthermore, we have found that RsuAgo can cleave double-stranded DNA in a low-salt buffer at 37 °C. These properties of RsuAgo provide a new tool for DNA and RNA manipulation at moderate temperatures.

Nonalternating purine pyrimidine sequences can form stable left-handed DNA duplex by strong topological constraint

In vivo, left-handed DNA duplex (usually refers to Z-DNA) is mainly formed in the region of DNA with alternating purine pyrimidine (APP) sequence and plays significant biological roles. It is well known that d(CG)_n sequence can form Z-DNA most easily under negative supercoil conditions, but its essence has not been well clarified. The study on sequence dependence of Z-DNA stability is very difficult without modification or inducers. Here, by the strong topological constraint caused by hybridization of two complementary short circular ssDNAs, left-handed duplex part was generated for various sequences, and their characteristics were investigated by using gel-shift after binding to specific proteins, CD and T_m analysis, and restriction enzyme cleavage. Under the strong topological constraint, non-APP sequences can also form left-handed DNA duplex as stable as that of APP sequences. As compared with non-APP sequences, the thermal stability difference for APP sequences between Z-form and B-form is smaller, which may be the reason that Z-DNA forms preferentially for APP ones. This result can help us to understand why nature selected APP sequences to regulate gene expression by transient Z-DNA formation, as well as why polymer with chirality can usually form both duplexes with left- or right-handed helix.

Calculation of Energy for RNA/RNA and DNA/RNA Duplex Formation by Molecular Dynamics Simulation

Abstract

The development of approaches for predictive calculation of hybridization properties of various nucleic acid (NA) derivatives is the basis for the rational design of the NA-based constructs. Modern advances in computer modeling methods provide the feasibility of these calculations. We have analyzed the possibility of calculating the energy of DNA/RNA and RNA/RNA duplex formation using representative sets of complexes (65 and 75 complexes, respectively). We used the classical molecular dynamics (MD) method, the MMPBSA or MMGBSA approaches to calculate the enthalpy (ΔH°) component, and the quasi-harmonic approximation (Q-Harm) or the normal mode analysis (NMA) methods to calculate the entropy (ΔS°) contribution to the Gibbs energy ($$\Delta G_{{37}}^{^\circ }$$ ) of the NA complex formation. We have found that the MMGBSA method in the analysis of the MD trajectory of only the NA duplex and the empirical linear approximation allow calculation of the enthalpy of formation of the DNA, RNA, and hybrid duplexes of various lengths and GC content with an accuracy of 8.6%. Within each type of complex, the combination of rather efficient MMGBSA and Q-Harm approaches being applied to the trajectory of only the bimolecular complex makes it possible to calculate the $$\Delta G_{{37}}^{^\circ }$$ of the duplex formation with an error value of 10%. The high accuracy of predictive calculation for different types of natural complexes (DNA/RNA, DNA/RNA, and RNA/RNA) indicates the possibility of extending the considered approach to analogs and derivatives of nucleic acids, which gives a fundamental opportunity in the future to perform rational design of new types of NA-targeted sequence-specific compounds.

Oligonucleotide conjugated antibody strategies for cyclic immunostaining

A number of highly multiplexed immunostaining and imaging methods have advanced spatial proteomics of cancer for improved treatment strategies. While a variety of methods have been developed, the most widely used methods are limited by harmful signal removal techniques, difficulties with reagent production and antigen sensitivity. Multiplexed immunostaining employing oligonucleotide (oligos)-barcoded antibodies is an alternative approach that is growing in popularity. However, challenges remain in consistent conjugation of oligos to antibodies with maintained antigenicity as well as non-destructive, robust and cost-effective signal removal methods. Herein, a variety of oligo conjugation and signal removal methods were evaluated in the development of a robust oligo conjugated antibody cyclic immunofluorescence (Ab-oligo cyCIF) methodology. Both non- and site-specific conjugation strategies were assessed to label antibodies, where site-specific conjugation resulted in higher retained binding affinity and antigen-specific staining. A variety of fluorescence signal removal methods were also evaluated, where incorporation of a photocleavable link (PCL) resulted in full fluorescence signal removal with minimal tissue disruption. In summary, this work resulted in an optimized Ab-oligo cyCIF platform capable of generating high dimensional images to characterize the spatial proteomics of the hallmarks of cancer.

A Methylation-Directed, Synthetic Pap Switch Based on Self-Complementary Regulatory DNA Reconstituted in an All E. coli Cell-Free Expression System

Pyelonephritis-associated pili (pap) enable migration of the uropathogenic Escherichia coli strain (UPEC) through the urinary tract. UPEC can switch between a stable 'ON phase' where the corresponding pap genes are expressed and a stable 'OFF phase' where their transcription is repressed. Hereditary DNA methylation of either one of two GATC motives within the regulatory region stabilizes the respective phase over many generations. The underlying molecular mechanism is only partly understood. Previous investigations suggest that in vivo phase-variation stability results from cooperative action of the transcriptional regulators Lrp and PapI. Here, we use an E. coli cell-free expression system to study molecular functions of the pap regulatory region based on a specially designed, synthetic construct flanked by two reporter genes encoding fluorescent proteins for simple readout. On the basis of our observations we suggest that besides Lrp, the conformation of the self-complementary regulatory DNA plays a strong role in the regulation of phase-variation. Our work not only contributes to better understand the phase variation mechanism, but it represents a successful start for mimicking stable, hereditary, and strong expression control based on methylation. The conformation of the regulatory DNA corresponds to a Holliday junction. Gene expression must be expected to respond if opposite arms of the junction are drawn outward.

Ionic Atmosphere Effect on the Absorption Spectrum of a Flavoprotein: A Reminder to Consider Solution Ions

This study utilizes the FMN-dependent NADH:quinone oxidoreductase from Pseudomonas aeruginosa PAO1 to investigate the effect of introducing an active site negative charge on the flavin absorption spectrum both in the absence and presence of a long-range electrostatic potential coming from solution ions. There were no observed changes in the flavin UV-visible spectrum when an active site tyrosine (Y277) becomes deprotonated in vitro. These results could only be reproduced computationally using average solvent electrostatic configuration (ASEC) QM/MM simulations that include both positive and negative solution ions. The same calculations performed with minimal ions to neutralize the total protein charge predicted that deprotonating Y277 would significantly alter the flavin absorption spectrum. Analyzing the distribution of solution ions indicated that the ions reorganize around the protein surface upon Y277 deprotonation to cancel the effect of the tyrosinate on the flavin absorption spectrum. Additional biochemical experiments were performed to test this hypothesis.

Complete Mesoscopic Parameterization of Single LNA Modifications in DNA Applied to Oncogene Probe Design

The use of mesoscopic models to describe the thermodynamic properties of locked nucleic acid (LNA)-modified nucleotides can provide useful insights into their properties, such as hydrogen-bonding and stacking interactions. In addition, the mesoscopic parameters can be used to optimize LNA insertion in probes, to achieve accurate melting temperature predictions, and to obtain duplex opening profiles at the base-pair level. Here, we applied this type of model to parameterize a large set of melting temperatures for LNA-modified sequences, from published sources, covering all possible nearest-neighbor configurations. We have found a very large increase in Morse potentials, which indicates very strong hydrogen bonding as the main cause of improved LNA thermodynamic stability. LNA-modified adenine-thymine (AT) was found to have similar hydrogen bonding to unmodified cytosine-guanine (CG) base pairs, while for LNA CG, we found exceptionally large hydrogen bonding. In contrast, stacking interactions, which were thought to be behind the stability of LNA, were similar to unmodified DNA in most cases. We applied the new LNA parameters to the design of BRAF, KRAS, and EGFR oncogene variants by testing all possible LNA modifications. Selected sequences were then synthesized and had their hybridization temperatures measured, achieving a prediction accuracy within 1 °C. We performed a detailed base-pair opening analysis to discuss specific aspects of these probe hybridizations that may be relevant for probe design.

Bottom-Up Fabrication of DNA-Templated Electronic Nanomaterials and Their Characterization

Bottom-up fabrication using DNA is a promising approach for the creation of nanoarchitectures. Accordingly, nanomaterials with specific electronic, photonic, or other functions are precisely and programmably positioned on DNA nanostructures from a disordered collection of smaller parts. These self-assembled structures offer significant potential in many domains such as sensing, drug delivery, and electronic device manufacturing. This review describes recent progress in organizing nanoscale morphologies of metals, semiconductors, and carbon nanotubes using DNA templates. We describe common substrates, DNA templates, seeding, plating, nanomaterial placement, and methods for structural and electrical characterization. Finally, our outlook for DNA-enabled bottom-up nanofabrication of materials is presented.

PhyloPrimer: a taxon-specific oligonucleotide design platform

Many environmental and biomedical biomonitoring and detection studies aim to explore the presence of specific organisms or gene functionalities in microbiome samples. In such cases, when the study hypotheses can be answered with the exploration of a small number of genes, a targeted PCR-approach is appropriate. However, due to the complexity of environmental microbial communities, the design of specific primers is challenging and can lead to non-specific results. We designed PhyloPrimer, the first user-friendly platform to semi-automate the design of taxon-specific oligos (i.e., PCR primers) for a gene of interest. The main strength of PhyloPrimer is the ability to retrieve and align GenBank gene sequences matching the user’s input, and to explore their relationships through an online dynamic tree. PhyloPrimer then designs oligos specific to the gene sequences selected from the tree and uses the tree non-selected sequences to look for and maximize oligo differences between targeted and non-targeted sequences, therefore increasing oligo taxon-specificity (positive/negative consensus approach). Designed oligos are then checked for the presence of secondary structure with the nearest-neighbor (NN) calculation and the presence of off-target matches with in silico PCR tests, also processing oligos with degenerate bases. Whilst the main function of PhyloPrimer is the design of taxon-specific oligos (down to the species level), the software can also be used for designing oligos to target a gene without any taxonomic specificity, for designing oligos from preselected sequences and for checking predesigned oligos. We validated the pipeline on four commercially available microbial mock communities using PhyloPrimer to design genus- and species-specific primers for the detection of Streptococcus species in the mock communities. The software performed well on these mock microbial communities and can be found at https://www.cerealsdb.uk.net/cerealgenomics/phyloprimer.

Dinucleotides as simple models of the base stacking-unstacking component of DNA 'breathing' mechanisms

Regulatory protein access to the DNA duplex ‘interior’ depends on local DNA ‘breathing’ fluctuations, and the most fundamental of these are thermally-driven base stacking-unstacking interactions. The smallest DNA unit that can undergo such transitions is the dinucleotide, whose structural and dynamic properties are dominated by stacking, while the ion condensation, cooperative stacking and inter-base hydrogen-bonding present in duplex DNA are not involved. We use dApdA to study stacking-unstacking at the dinucleotide level because the fluctuations observed are likely to resemble those of larger DNA molecules, but in the absence of constraints introduced by cooperativity are likely to be more pronounced, and thus more accessible to measurement. We study these fluctuations with a combination of Molecular Dynamics simulations on the microsecond timescale and Markov State Model analyses, and validate our results by calculations of circular dichroism (CD) spectra, with results that agree well with the experimental spectra. Our analyses show that the CD spectrum of dApdA is defined by two distinct chiral conformations that correspond, respectively, to a Watson–Crick form and a hybrid form with one base in a Hoogsteen configuration. We find also that ionic structure and water orientation around dApdA play important roles in controlling its breathing fluctuations.

Selective extraction of low-abundance BRAF V600E mutation from plasma, urine, and sputum using ion-tagged oligonucleotides and magnetic ionic liquids

Sequence-specific DNA extractions have the potential to improve the detection of low-abundance mutations from complex matrices, making them ideal for circulating tumor DNA analysis during the early stages of cancer. Ion-tagged oligonucleotides (ITOs) are oligonucleotides modified with an allylimidazolium salt via thiolene click chemistry. The allylimidazolium-based tag allows the ITO-DNA duplex to be selectively captured by a hydrophobic magnetic ionic liquid (MIL). In this study, the selectivity of the ITO-MIL method was examined by extracting low abundance of the BRAF V600E mutation-a common single-nucleotide polymorphism associated with several different cancers-from diluted human plasma, artificial urine, and diluted artificial sputum. Quantitative polymerase chain reaction (qPCR) was not able to distinguish a 9% BRAF V600E standard (50 fg·μL^-1 BRAF V600E, 500 fg·μL^-1 wild-type BRAF) from the 100% wild-type BRAF (50 fg·μL^-1) standard. However, introducing the ITO-MIL extraction prior to qPCR allowed for samples consisting of 0.1% BRAF V600E (50 fg·μL^-1 V600E BRAF, 50,000 fg·μL^-1 wild-type BRAF) to be distinguished from the 100% wild-type BRAF standard. Ion-tagged oligonucleotides (ITOs) are combined with magnetic ionic liquids (MILs) to extract low-abundance BRAF V600E mutation from diluted human plasma, artificial urine, and diluted artificial sputum. The ITO-MIL extraction prior to qPCR allowed for samples consisting of 0.1% BRAF V600E to be distinguished from the 100% wild-type BRAF standard.

Presence of an EML4-ALK gene fusion detected by microfluidic chip DNA hybridization

Non-small cell lung cancer (NSCLC) accounts for ∼80-85% of all lung cancer cases, and the EML4-ALK fusion oncogene is a well-known contributor to NSCLC cases. Expensive methods such as FISH, IHC, and NGS have been used to detect the EML4-ALK fusion oncogene. Here, a cost-effective and facile method of detecting and differentiating an EML4-ALK fusion oncogene from the wild-type gene has been accomplished by DNA hybridization using the microfluidic biochip. First, oligonucleotide probes were confirmed for successful detection of immobilized sense strands. Second, capture of the sense PCR product strands (fusion and WT) and their subsequent detection and differentiation were accomplished. Our proof-of-concept study shows the ability to detect 1% fusion products, among WT ones.

Thermodynamics of DNA Hybridization from Atomistic Simulations

Studying DNA hybridization equilibrium at atomistic length scales, either via molecular dynamics (MD) or through commonly used advanced sampling approaches, is notoriously difficult. In this work, we describe an order-parameter-based advanced sampling technique to calculate the free energy of hybridization, and estimate the melting temperature of DNA oligomers at atomistic resolution. The free energy landscapes are reported as a function of a native-topology-based order parameter for the Drew-Dickerson dodecamer and for a range of DNA decamer sequences of different GC content. Our estimated melting temperatures match the experimental numbers within ±15 °C. As a test of the numerical reliability of the procedures employed, it was verified that the predicted free energy surfaces and melting temperatures of the d- and l-enantiomers of the Drew-Dickerson dodecamer were indistinguishable within numerical accuracy.

Carbon nanotube transmembrane channel formation and single-stranded DNA spontaneous internalization: a dissipative particle dynamics study

Single-walled carbon nanotube (SWCNT) transmembrane channel formation in a pure 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) bilayer, and the spontaneous internalization of single-stranded DNA (ssDNA) into the formed pore were simulated. A combination of computational techniques, Dissipative Particle Dynamics-Monte Carlo hybrid simulations and quantum mechanical calculations at the hybrid-DFT level, was used as a new proposal to perform DPD simulations granting specific chemical identity to the model particles. The simulated transmembrane channels showed that, in the case of pristine SWCNTs and upon increasing the nanotube length, a higher tilt angle with respect to the bilayer normal is observed and more time is needed for the nanotube to stabilize. On the other hand, for SWCNTs with polar rims an almost perpendicular orientation is preferred with less than 15° of tilt with respect to the bilayer normal once the nanotubes have pierced both monolayers. These findings are supported by experimental observations where CNTs of average inner diameters of 1.51 ± 0.21 nm and lengths in the 5-15 nm range were inserted in DOPC membranes [J. Geng, et al., Nature, 2014, 514(7524), 612-615]. Moreover, the narrower the SWCNTs, the slower the spontaneous internalization of ssDNA becomes, and ssDNA ends hydrophobically trapped inside the artificial pore. A dependence on the nucleotide content is found indicating that the higher the presence of adenine and thymine in the ssDNA chains the slower the internalization becomes, in agreement with the experimental [A. M. Ababneh, et al., Biophys. J., 2003, 85(2), 1111-1127] and predicted solvation tendency in water for nucleic acid bases.

Stabilization of DNA by sodium and magnesium ions during the synthesis of DNA-bridged gold nanoparticles

Nanostructures synthesized using DNA-conjugated gold nanoparticles have a wide range of applications in the field of biosensorics. The stability of the DNA duplex plays a critical role as it determines the final geometry of these nanostructures. The main way to control DNA stability is to maintain a high ionic strength of the buffer solution; at the same time, high salt concentrations lead to an aggregation of nanoparticles. In this study, by means of the instrumentality of DNA-bridged seeds using tris(hydroxymethyl)aminomethane as a soft reducing agent the dumbbell-like gold nanoparticles up to 35 nm were synthesized with a high concentration of sodium ions of up to 100 mM and magnesium ions up to 1 mM. We also examined at the atomic level the details of the effect of the gold nanoparticle surface, as well as Na⁺ and Mg²⁺ ions, on the stability of nucleotide pairs located in close proximity to the grafting site.

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.

Liquid Crystal Peptide/DNA Coacervates in the Context of Prebiotic Molecular Evolution

Liquid–liquid phase separation (LLPS) phenomena are ubiquitous in biological systems, as various cellular LLPS structures control important biological processes. Due to their ease of in vitro assembly into membraneless compartments and their presence within modern cells, LLPS systems have been postulated to be one potential form that the first cells on Earth took on. Recently, liquid crystal (LC)-coacervate droplets assembled from aqueous solutions of short double-stranded DNA (s-dsDNA) and poly-L-lysine (PLL) have been reported. Such LC-coacervates conjugate the advantages of an associative LLPS with the relevant long-range ordering and fluidity properties typical of LC, which reflect and propagate the physico-chemical properties of their molecular constituents. Here, we investigate the structure, assembly, and function of DNA LC-coacervates in the context of prebiotic molecular evolution and the emergence of functional protocells on early Earth. We observe through polarization microscopy that LC-coacervate systems can be dynamically assembled and disassembled based on prebiotically available environmental factors including temperature, salinity, and dehydration/rehydration cycles. Based on these observations, we discuss how LC-coacervates can in principle provide selective pressures effecting and sustaining chemical evolution within partially ordered compartments. Finally, we speculate about the potential for LC-coacervates to perform various biologically relevant properties, such as segregation and concentration of biomolecules, catalysis, and scaffolding, potentially providing additional structural complexity, such as linearization of nucleic acids and peptides within the LC ordered matrix, that could have promoted more efficient polymerization. While there are still a number of remaining open questions regarding coacervates, as protocell models, including how modern biologies acquired such membraneless organelles, further elucidation of the structure and function of different LLPS systems in the context of origins of life and prebiotic chemistry could provide new insights for understanding new pathways of molecular evolution possibly leading to the emergence of the first cells on Earth.

Necessity and Challenges of Sample Preconcentration in Analysis of Multiple MicroRNAs by Capillary Electrophoresis

Thousands of putative microRNA (miRNA)-based cancer biomarkers have been reported, but none has been validated for approval by the Food and Drug Administration. One of the reasons for this alarming discrepancy is the lack of a method that is sufficiently robust for carrying out validation studies, which may require analysis of samples from hundreds of patients across multiple institutions and pooling the results together. The capillary electrophoresis (CE)-based hybridization assay proved to be more robust than reversed transcription polymerase chain reaction (the current standard), but its limit of quantification (LOQ) exceeds 10 pM while miRNA concentrations in cell lysates are below 1 pM. Thus, CE-based separation must be preceded by on-column sample preconcentration. Here, we explain the challenges of sample preconcentration for CE-based miRNA analyses and introduce a preconcentration method that can suit CE-based miRNA analysis utilizing peptide nucleic acid (PNA) hybridization probes. The method combines field-amplified sample stacking (FASS) with isotachophoresis (ITP). We proved that FASS-ITP could retain and concentrate both near-neutral PNA with highly negatively charged PNA-miRNA hybrids. We demonstrated that preconcentration by FASS-ITP could be combined with the CE-based separation of the unreacted PNA probes from the PNA-miRNA hybrids and facilitate improvement in LOQ by a factor of 140, down to 0.1 pM. Finally, we applied FASS-ITP-CE for the simultaneous detection of two miRNAs in crude cell lysates and proved that the method was robust when used in complex biological matrices. The 140-fold improvement in LOQ and the robustness to biological matrices will significantly expand the applicability of CE-based miRNA analysis, bringing it closer to becoming a practical tool for validation of miRNA biomarkers.

Duplex-tetraplex equilibria in guanine- and cytosine-rich DNA

Noncanonical four-stranded DNA structures, including G-quadruplexes and i-motifs, have been discovered in the cell and are implicated in a variety of genomic regulatory functions. The tendency of a specific guanine- and cytosine-rich region of genomic DNA to adopt a four-stranded conformation depends on its ability to overcome the constraints of duplex base-pairing by undergoing consecutive duplex-to-coil and coil-to-tetraplex transitions. The latter ability is determined by the balance between the free energies of participating ordered and disordered structures. In this review, we present an overview of the literature on the stability of G-quadruplex and i-motif structures and discuss the extent of duplex-tetraplex competition as a function of the sequence context of the DNA and environmental conditions including temperature, pH, salt, molecular crowding, and the presence of G-quadruplex-binding ligands. We outline how the results of in vitro studies can be expanded to understanding duplex-tetraplex equilibria in vivo.

Heat Capacity Changes (ΔCp) for Interconversions between Differentially-Ordered DNA States within Physiological Temperature Domains: Implications for Biological Regulatory Switches

Knowledge of differences in heat capacity changes (ΔC_p) between biopolymer states provides essential information about the temperature dependence of the thermodynamic properties of these states, while also revealing insights into the nature of the forces that drive the formation of functional and dysfunctional biopolymer "order." In contrast to proteins, for nucleic acids there is a dearth of direct experimental determination of this information-rich parameter, a deficiency that compromises interpretations of the ever-increasing thermodynamic analyses of nucleic acid properties; particularly as they relate to differential nucleic acid (meta)stability states and their potential biological functions. Here we demonstrate that such heat capacity differences, in fact, exist not only between traditionally measured native to fully unfolded (assumed "random coil") DNA states, but also between competing order-to-order transformations. We illustrate the experimental approach by measuring the heat capacity change between "native"/ordered, sequence homologous, "isomeric" DNA states that differ in conformation but not sequence. Importantly, these heat capacity differences occur within biologically relevant temperature ranges. In short, we describe a new and general method to measure the value of such heat capacity differences anywhere in experimentally accessible conformational and temperature space; in this case, between two metastable bulge loop states, implicated in DNA expansion diseases, and their competing, fully paired, thermodynamically more stable duplex states. This measurement reveals a ΔC_p of 61 ± 7 cal mol_bp ^-1 K ^-1. Such heat capacity differences between competing DNA "native" ensemble states must be considered when evaluating equilibria between different DNA "ordered" conformations, including the assessment of the differential stabilizing forces and potential biological functions of competing DNA "structured" motifs.

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.