Predicting sequence-dependent melting stability of short duplex DNA oligomers

Unleashing inkjet-printed nanostructured electrodes and battery-free potentiostat for the DNA-based multiplexed detection of SARS-CoV-2 genes

Following the global COVID-19 pandemic triggered by SARS-CoV-2, the need for rapid, specific and cost-effective point-of-care diagnostic solutions remains paramount. Even though COVID-19 is no longer a public health emergency, the disease still poses a global threat leading to deaths, and it continues to change with the risk of new variants emerging causing a new surge in cases and deaths. Here, we address the urgent need for rapid, cost-effective and point-of-care diagnostic solutions for SARS-CoV-2. We propose a multiplexed DNA-based sensing platform that utilizes inkjet-printed nanostructured gold electrodes and an inkjet-printed battery-free near-field communication (NFC) potentiostat for the simultaneous quantitative detection of two SARS-CoV-2 genes, the ORF1ab and the N gene. The detection strategy based on the formation of an RNA-DNA sandwich structure leads to a highly specific electrochemical output. The inkjet-printed nanostructured gold electrodes providing a large surface area enable efficient binding and increase the sensitivity. The inkjet-printed battery-free NFC potentiostat enables rapid measurements and real-time data analysis via a smartphone application, making the platform accessible and portable. With the advantages of speed (5 min), simplicity, sensitivity (low pM range, ∼450% signal gain) and cost-effectiveness, the proposed platform is a promising alternative for point-of-care diagnostics and high-throughput analysis that complements the COVID-19 diagnostic toolkit.

The impact of the sequence-dependent physical properties of DNA on chromatin dynamics

The local mechanical properties of DNA depend on local sequence. Here we review recent genomic, structural, and computational efforts at deciphering the "mechanical code", i.e., the mapping between sequence and mechanics. We then discuss works that suggest how evolution has exploited the mechanical code to control the energetics of DNA-deforming biological processes such as nucleosome organization, transcription factor binding, DNA supercoiling, gene regulation, and 3D chromatin organization. As a whole, these recent works suggest that DNA sequence in diverse organisms can encode regulatory information governing diverse processes via the mechanical code.

Generation of DNA oligomers with similar chemical kinetics via in-silico optimization

Networks of interacting DNA oligomers are useful for applications such as biomarker detection, targeted drug delivery, information storage, and photonic information processing. However, differences in the chemical kinetics of hybridization reactions, referred to as kinetic dispersion, can be problematic for some applications. Here, it is found that limiting unnecessary stretches of Watson-Crick base pairing, referred to as unnecessary duplexes, can yield exceptionally low kinetic dispersions. Hybridization kinetics can be affected by unnecessary intra-oligomer duplexes containing only 2 base-pairs, and such duplexes explain up to 94% of previously reported kinetic dispersion. As a general design rule, it is recommended that unnecessary intra-oligomer duplexes larger than 2 base-pairs and unnecessary inter-oligomer duplexes larger than 7 base-pairs be avoided. Unnecessary duplexes typically scale exponentially with network size, and nearly all networks contain unnecessary duplexes substantial enough to affect hybridization kinetics. A new method for generating networks which utilizes in-silico optimization to mitigate unnecessary duplexes is proposed and demonstrated to reduce in-vitro kinetic dispersions as much as 96%. The limitations of the new design rule and generation method are evaluated in-silico by creating new oligomers for several designs, including three previously programmed reactions and one previously engineered structure.

Calorimetric analysis using DNA thermal stability to determine protein concentration

It was recently reported for two globular proteins and a short DNA hairpin in NaCl buffer that values of the transition heat capacities, Cp,DNA and Cp,PRO, for equal concentrations (mg/mL) of DNA and proteins, are essentially equivalent (differ by less than 1%). Additional evidence for this equivalence is presented that reveals this phenomenon does not depend on DNA sequence, buffer salt, or Tm. Sequences of two DNA hairpins were designed to confer a near 20°C difference in their Tm's. For the molecules, in NaCl and CsCl buffer the evaluated Cp,PRO and Cp,DNA were equivalent. Based on the equivalence of transition heat capacities, a calorimetric method was devised to determine protein concentrations in pure and complex solutions. The scheme uses direct comparisons between the thermodynamic stability of a short DNA hairpin standard of known concentration, and thermodynamic stability of protein solutions of unknown concentrations. In all cases, evaluated protein concentrations determined from the DNA standard curve agreed with the UV-Vis concentration for monomeric proteins. For samples of multimeric proteins, streptavidin (tetramer), Herpes Simplex Virus glycoprotein D (trimer/dimer), and a 16 base pair DNA duplex (dimer), evaluated concentrations were greater than determined by UV-Vis by factors of 3.94, 2.65, and 2.15, respectively.

The stability and number of nucleating interactions determine DNA hybridization rates in the absence of secondary structure

The kinetics of DNA hybridization are fundamental to biological processes and DNA-based technologies. However, the precise physical mechanisms that determine why different DNA sequences hybridize at different rates are not well understood. Secondary structure is one predictable factor that influences hybridization rates but is not sufficient on its own to fully explain the observed sequence-dependent variance. In this context, we measured hybridization rates of 43 different DNA sequences that are not predicted to form secondary structure and present a parsimonious physically justified model to quantify our observations. Accounting only for the combinatorics of complementary nucleating interactions and their sequence-dependent stability, the model achieves good correlation with experiment with only two free parameters. Our results indicate that greater repetition of Watson-Crick pairs increases the number of initial states able to proceed to full hybridization, with the stability of those pairings dictating the likelihood of such progression, thus providing new insight into the physical factors underpinning DNA hybridization rates.

Developing predictive hybridization models for phosphorothioate oligonucleotides using high-resolution melting

The ability to predict nucleic acid hybridization energies has been greatly enabling for many applications, but predictive models require painstaking experimentation, which may limit expansion to non-natural nucleic acid analogues and chemistries. We have assessed the utility of dye-based, high-resolution melting (HRM) as an alternative to UV-Vis determinations of hyperchromicity in order to more quickly acquire parameters for duplex stability prediction. The HRM-derived model for phosphodiester (PO) DNA can make comparable predictions to previously established models. Using HRM, it proved possible to develop predictive models for DNA duplexes containing phosphorothioate (PS) linkages, and we found that hybridization stability could be predicted as a function of sequence and backbone composition for a variety of duplexes, including PS:PS, PS:PO, and partially modified backbones. Individual phosphorothioate modifications destabilize helices by around 0.12 kcal/mol on average. Finally, we applied these models to the design of a catalytic hairpin assembly circuit, an enzyme-free amplification method used for nucleic acid-based molecular detection. Changes in PS circuit behavior were consistent with model predictions, further supporting the addition of HRM modeling and parameters for PS oligonucleotides to the rational design of nucleic acid hybridization.

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.

Design and Experimental Evaluation of a New RNA-FISH Probe to Detect and Identify Paenibacillus sp

Paenibacillus, rod-saped gram-positive endospores forming aerobic or facultative anaerobic bacteria, colonize diverse ecosystems and are involved in the biodegradation of cultural heritage assets. Biodeteriogenic microorganisms can be easily detected/identified by ribonucleic acid- fluorescent in situ hybridization RNA-FISH with specific probes. In this work, probes designed in silico were analyzed to calculate hybridization efficiency and specificity by varying the formamide concentration in the hybridization. The Pab489 probe showed excellent in silico performance with high theoretical maximum efficiency hybridization (99.99%) and specificity and was selected for experimental assays with target Paenibacillus sp. and non-target biodeteriogenic microorganisms. Results assessed by epifluorescence microscopy and flow cytometry revealed that, regardless of the formamide concentration, it was possible to observe that the Pab489-Cy3 probe had a similar signal intensity to the EUB338-Cy3 probe (positive control), so the presence of formamide, a highly toxic and carcinogenic compound used to aid the hybridization process, is not necessary. The designed probe used in FISH assays allows specific in situ identification of Paenibacillus spp. in microbial communities in a culture-independent way. This approach can be employed for screening Paenibacillus spp., showing great potential for future application in biodeterioration of heritage assets, in the search for Paenibacillus strains that produce compounds with biotechnological or medical potential.

CROPSR: an automated platform for complex genome-wide CRISPR gRNA design and validation

BACKGROUND

CRISPR/Cas9 technology has become an important tool to generate targeted, highly specific genome mutations. The technology has great potential for crop improvement, as crop genomes are tailored to optimize specific traits over generations of breeding. Many crops have highly complex and polyploid genomes, particularly those used for bioenergy or bioproducts. The majority of tools currently available for designing and evaluating gRNAs for CRISPR experiments were developed based on mammalian genomes that do not share the characteristics or design criteria for crop genomes.

RESULTS

We have developed an open source tool for genome-wide design and evaluation of gRNA sequences for CRISPR experiments, CROPSR. The genome-wide approach provides a significant decrease in the time required to design a CRISPR experiment, including validation through PCR, at the expense of an overhead compute time required once per genome, at the first run. To better cater to the needs of crop geneticists, restrictions imposed by other packages on design and evaluation of gRNA sequences were lifted. A new machine learning model was developed to provide scores while avoiding situations in which the currently available tools sometimes failed to provide guides for repetitive, A/T-rich genomic regions. We show that our gRNA scoring model provides a significant increase in prediction accuracy over existing tools, even in non-crop genomes.

CONCLUSIONS

CROPSR provides the scientific community with new methods and a new workflow for performing CRISPR/Cas9 knockout experiments. CROPSR reduces the challenges of working in crops, and helps speed gRNA sequence design, evaluation and validation. We hope that the new software will accelerate discovery and reduce the number of failed experiments.

Determining Sequence-Dependent DNA Oligonucleotide Hybridization and Dehybridization Mechanisms Using Coarse-Grained Molecular Simulation, Markov State Models, and Infrared Spectroscopy

A robust understanding of the sequence-dependent thermodynamics of DNA hybridization has enabled rapid advances in DNA nanotechnology. A fundamental understanding of the sequence-dependent kinetics and mechanisms of hybridization and dehybridization remains comparatively underdeveloped. In this work, we establish new understanding of the sequence-dependent hybridization/dehybridization kinetics and mechanism within a family of self-complementary pairs of 10-mer DNA oligomers by integrating coarse-grained molecular simulation, machine learning of the slow dynamical modes, data-driven inference of long-time kinetic models, and experimental temperature-jump infrared spectroscopy. For a repetitive ATATATATAT sequence, we resolve a rugged dynamical landscape comprising multiple metastable states, numerous competing hybridization/dehybridization pathways, and a spectrum of dynamical relaxations. Introduction of a G:C pair at the terminus (GATATATATC) or center (ATATGCATAT) of the sequence reduces the ruggedness of the dynamics landscape by eliminating a number of metastable states and reducing the number of competing dynamical pathways. Only by introducing a G:C pair midway between the terminus and the center to maximally disrupt the repetitive nature of the sequence (ATGATATCAT) do we recover a canonical "all-or-nothing" two-state model of hybridization/dehybridization with no intermediate metastable states. Our results establish new understanding of the dynamical richness of sequence-dependent kinetics and mechanisms of DNA hybridization/dehybridization by furnishing quantitative and predictive kinetic models of the dynamical transition network between metastable states, present a molecular basis with which to understand experimental temperature jump data, and furnish foundational design rules by which to rationally engineer the kinetics and pathways of DNA association and dissociation for DNA nanotechnology applications.

The shaping of a molecular linguist: How a career studying DNA energetics revealed the language of molecular communication

My personal and professional journeys have been far from predictable based on my early childhood. Owing to a range of serendipitous influences, I miraculously transitioned from a rebellious, apathetic teenage street urchin who did poorly in school to a highly motivated, disciplined, and ambitious academic honors student. I was the proverbial “late bloomer.” Ultimately, I earned my PhD in biophysical chemistry at Yale, followed by a postdoc fellowship at Berkeley. These two meccas of thermodynamics, coupled with my deep fascination with biology, instilled in me a passion to pursue an academic career focused on mapping the energy landscapes of biological systems. I viewed differential energetics as the language of molecular communication that would dictate and control biological structures, as well as modulate the modes of action associated with biological functions. I wanted to be a “molecular linguist.” For the next 50 years, my group and I used a combination of spectroscopic and calorimetric techniques to characterize the energy profiles of the polymorphic conformational space of DNA molecules, their differential ligand-binding properties, and the energy landscapes associated with mutagenic DNA damage recognition, repair, and replication. As elaborated below, the resultant energy databases have enabled the development of quantitative molecular biology through the rational design of primers, probes, and arrays for diagnostic, therapeutic, and molecular-profiling protocols, which collectively have contributed to a myriad of biomedical assays. Such profiling is further justified by yielding unique energy-based insights that complement and expand elegant, structure-based understandings of biological processes.

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.

High-throughput methods for measuring DNA thermodynamics

Understanding the thermodynamics of DNA motifs is important for prediction and design of probes and primers, but melt curve analyses are low-throughput and produce inaccurate results for motifs such as bulges and mismatches. Here, we developed a new, accurate and high-throughput method for measuring DNA motif thermodynamics called TEEM (Toehold Exchange Energy Measurement). It is a refined framework of comparing two toehold exchange reactions, which are competitive strand displacement between oligonucleotides. In a single experiment, TEEM can measure over 1000 ΔG° values with standard error of roughly 0.05 kcal/mol.

Design of Split Proximity Circuit as a Plug-and-Play Translator for Point Mutation Discrimination

Point mutations are a common form of genetic variation and have been identified as important disease biomarkers. Conventional methods for analyzing point mutations, e.g., polymerase chain reaction (PCR), are based on differences in thermal stability of the DNA duplex, which require extensive optimization of the reaction condition and nontrivial design of sequence-selective primers. This motivated the design of molecular translators to convert molecular inputs into generic output sequences, which allows for the target recognition and signal generation regions to be designed independently. In this work, we propose a translator design based on the concept of split proximity circuit (SPC) to achieve both high sequence selectivity and assay robustness using a universal reaction condition, i.e., room temperature and constant ionic concentration. We discussed the design aspects of the SPC recognition regions and demonstrated its plug-and-play capability to discriminate different point mutations for both DNA (seven G6PD mutations) and RNA (let-7 microRNA family members) targets while retaining the same signal generation region. Despite its simple design and nonstringent assay condition requirements, the SPC retained good analytical performance to detect subnanomolar target concentration within a reasonable time of an hour.

Interlocked DNA Nanojoints for Reversible Thermal Sensing

The ability to precisely measure and monitor temperature at high resolution at the nanoscale is an important task for better understanding the thermodynamic properties of functional entities at the nanoscale in complex systems, or at the level of a single cell. However, the development of high-resolution and robust thermal nanosensors is challenging. The design, assembly, and characterization of a group of thermal-responsive deoxyribonucleic acid (DNA) joints, consisting of two interlocked double-stranded DNA (dsDNA) rings, is described. The DNA nanojoints reversibly switch between the static and mobile state at different temperatures without a special annealing process. The temperature response range of the DNA nanojoint can be easily tuned by changing the length or the sequence of the hybridized region in its structure, and because of its interlocked structure the temperature response range of the DNA nanojoint is largely unaffected by its own concentration; this contrasts with systems that consist of separated components.

Hydrogen bonding of ionic liquids in the groove region of DNA controls the extent of its stabilization: synthesis, spectroscopic and simulation studies

Ionic liquids (ILs) have been extensively used for stabilization and long-term DNA storage. However, molecular level understanding of the role of the hydrogen bond of DNA with ILs in its stabilization is still inadequate. Two ILs, namely, 1,1,3,3-tetramethylguanidinium acetate (TMG) and 2,2-diethyl-1,1,3,3-tetramethylguanidinium acetate (DETMG), have been synthesized, of which TMG has a hydrogen bonding N-H group whereas DETMG does not contain any hydrogen bonding site. It has been found that both TMG and DETMG cations interact in the groove region of DNA; however, their mode of interaction is distinctly different, which causes the stabilization of DNA in the presence of TMG, whereas the effect is opposite in the case of DETMG. It is apparent from the data that only the accommodation of ILs in the groove region is not enough for the stabilization of DNA. MD simulation and spectroscopic studies combinedly indicate that the hydrogen bonding capability of the TMG cation enhances the hydrogen bonding between the Watson-Crick base pairs of DNA, resulting in its stabilization. In contrast, the bigger size as well as the absence of the hydrogen bonding site of the DETMG cation perturbs the minor groove width and base pair step parameters of DNA during its intrusion into the minor groove, which decreases the hydrogen bond between the Watson-Crick base pairs of DNA, leading to destabilization.

A convenient renewable surface plasmon resonance chip for relative quantification of genetically modified soybean in food and feed

The cultivation of genetically modified organisms (GMO) continues to expand worldwide. Still, many consumers express concerns about the use of GMO in food or feed, and many countries have legislated on labelling systems to indicate the presence of GMO in commercial products. To deal with the increased number of GMO events and to address related regulations, alternative detection methods for GMO inspection are required. In this work, a genosensor based on Surface Plasmon Resonance under continuous flow was developed for the detection and quantification of a genetically modified soybean (event GTS 40-3-2). In a single chip, the simultaneous detection of the event-specific and the taxon-specific samples were achieved, whose detection limits were 20 pM and 16 pM, respectively. The reproducibility was 1.4%, which supports the use of the chip as a reliable and cost-effective alternative to other DNA-based techniques. The results indicate that the proposed method is a versatile tool for GMO quantification in food and feed samples.

A simple procedure for detecting Dekkera bruxellensis in wine environment by RNA-FISH using a novel probe

Dekkera bruxellensis, considered the major microbial contaminant in wine production, produces 4-ethylphenol, a cause of unpleasant odors. Thus, identification of this yeast before wine spoilage is crucial. Although challenging, it could be achieved using a simple technique: RNA-FISH. To reach it is necessary to design probes that allow specific detection/identification of D. bruxellensis among the wine microorganisms and in the wine environment and, if possible, using low formamide concentrations. Therefore, this study was focused on: a) designing a DNA-FISH probe to identify D. bruxellensis that matches these requirements and b) determining the applicability of the RNA-FISH procedure after the end of the alcoholic fermentation and in wine. A novel DNA-FISH D. bruxellensis probe with good performance and specificity was designed. The application of this probe using an in-suspension RNA-FISH protocol (applying only 5% of formamide) allowed the early detection/identification of D. bruxellensis at low cell densities (5 × 10² cell/mL). This was possible by flow cytometry independently of the growth stage of the target cells, both at the end of the alcoholic fermentation and in wine even in the presence of high S. cerevisiae cell densities. Thus, this study aims to contribute to facilitate the identification of D. bruxellensis before wine spoilage occurs, preventing economic losses to the wine industry.

Gold-DNA nanosunflowers for efficient gene silencing with controllable transformation

The development of an efficient delivery system for enhanced and controlled gene interference-based therapeutics is still facing great challenges. Fortunately, the flourishing field of nanotechnology provides more effective strategies for nucleic acid delivery. Here, the triplex-forming oligonucleotide sequence and its complementary strand were used to mediate self-assembly of ultrasmall gold nanoparticles. The obtained sunflower-like nanostructures exhibited strong near-infrared (NIR) absorption and photothermal conversion ability. Upon NIR irradiation, the large-sized nanostructure could disassemble and generate ultrasmall nanoparticles modified with c-myc oncogene silencing sequence, which could directly target the cell nucleus. Moreover, the controlled gene silencing effect could be realized by synergistically controlling the preincubation time with the self-assembled nanostructure (in vitro and in vivo) and NIR irradiation time point. This study provides a new approach for constructing more efficient and tailorable nanocarriers for gene interference applications.

In-silico Design of DNA Oligonucleotides: Challenges and Approaches

DNA oligonucleotides are essential components of a high number of technologies in molecular biology. The key event of each oligonucleotide-based assay is the specific binding between oligonucleotides and their target DNA. However, single-stranded DNA molecules also tend to bind to unintended targets or themselves. The probability of such unspecific binding increases with the complexity of an assay. Therefore, accurate data management and design workflows are necessary to optimize the in-silico design of primers and probes. Important considerations concerning computational infrastructure and run time need to be made for both data management and the design process. Data retrieval, data updates, storage, filtering and analysis are the main parts of a sequence data management system. Each part needs to be well-implemented as the resulting sequences form the basis for the oligonucleotide design. Important key features, such as the oligonucleotide length, melting temperature, secondary structures and primer dimer formation, as well as the specificity, should be considered for the in-silico selection of oligonucleotides. The development of an efficient oligonucleotide design workflow demands the right balance between the precision of the applied computer models, the general expenditure of time, and computational workload. This paper gives an overview of important parameters during the design process, starting from the data retrieval, up to the design parameters for optimized oligonucleotide design.

Revealing the Presence of a Symbolic Sequence Representing Multiple Nucleotides Based on K-Means Clustering of Oligonucleotides

In biological systems, a few sequence differences diversify the hybridization profile of nucleotides and enable the quantitative control of cellular metabolism in a cooperative manner. In this respect, the information required for a better understanding may not be in each nucleotide sequence, but representative information contained among them. Existing methodologies for nucleotide sequence design have been optimized to track the function of the genetic molecule and predict interaction with others. However, there has been no attempt to extract new sequence information to represent their inheritance function. Here, we tried to conceptually reveal the presence of a representative sequence from groups of nucleotides. The combined application of the K-means clustering algorithm and the social network analysis theorem enabled the effective calculation of the representative sequence. First, a “common sequence” is made that has the highest hybridization property to analog sequences. Next, the sequence complementary to the common sequence is designated as a ‘representative sequence’. Based on this, we obtained a representative sequence from multiple analog sequences that are 8–10-bases long. Their hybridization was empirically tested, which confirmed that the common sequence had the highest hybridization tendency, and the representative sequence better alignment with the analogs compared to a mere complementary.

Classical Triplex Molecular Beacons for MicroRNA-21 and Vascular Endothelial Growth Factor Detection

Triplex molecular beacons (tMBs) possess great potential in biological sensing because of the pH responsiveness and controllability of binding strength. Here, we systematically investigate and rationally design a classical tMB for convenient detection of microRNA-21, a well-known biomarker of cardio-cerebrovascular diseases. In the tMB, we employ the complementary sequence of miR-21 as the loop and the sequences of protonated cytosine-guanine-cytosine (C-G•C+) and thymine-adenine-thymine (T-A•T) as the triplex stem, in which both the Watson-Crick and Hoogsteen base-pairing control the binding strength in cooperation. It is demonstrated for the first time that the presence of miR-21 would only break the Hoogsteen base-pairing in the stem and hybridize with the tMB to form the rigid heterozygous hybrid duplex structure. These would hinder the fluorescence resonance energy transfer (FRET) between the fluorophore (FAM) and quencher (BHQ1) labeled at the ends of the oligonucleotide, and the fluorescence recovery degree of FAM can be used as the standard to quantitate the miR-21. More significantly, the excellent adjustability and sensitivity of our tMBs have been confirmed by constructing the corresponding duplex molecular beacon (dMB) for comparison. The fluorophore FAM in the tMB could be replaced by the fluorescent DNA/silver nanoclusters, which exhibits the universal applicability of energy donor and receptor selection for tMB. Furthermore, our proposed tMB could also be developed as an aptasensor for the detection of vascular endothelial growth factor (VEGF) by only introducing the complementary sequence of its aptamer into the tMB. This work is of great significance for the systematic study of tMBs for the detection of biomarkers such as nucleic acids and proteins.

Derivation of nearest-neighbor DNA parameters in magnesium from single molecule experiments

DNA hybridization is an essential molecular reaction in biology with many applications. The nearest-neighbor (NN) model for nucleic acids predicts DNA thermodynamics using energy values for the different base pair motifs. These values have been derived from melting experiments in monovalent and divalent salt and applied to predict melting temperatures of oligos within a few degrees. However, an improved determination of the NN energy values and their salt dependencies in magnesium is still needed for current biotechnological applications seeking high selectivity in the hybridization of synthetic DNAs. We developed a methodology based on single molecule unzipping experiments to derive accurate NN energy values and initiation factors for DNA. A new set of values in magnesium is derived, which reproduces unzipping data and improves melting temperature predictions for all available oligo lengths, in a range of temperature and salt conditions where correlation effects between the magnesium bound ions are weak. The NN salt correction parameters are shown to correlate to the GC content of the NN motifs. Our study shows the power of single-molecule force spectroscopy assays to unravel novel features of nucleic acids such as sequence-dependent salt corrections.

DNA melting and energetics of the double helix

Studying melting and energetics of the DNA double helix has been one of the major topics of molecular biophysics over the past six decades. The main objective of this article is to overview the current state of the field and to emphasize that there are still serious gaps in our understanding of the issue. We start with a concise description of the commonly accepted theoretical model of the DNA melting. We then concentrate on studies devoted to the comparison with experiment of theoretically predicted melting profiles of DNAs with known sequences. For long DNA molecules, such comparison is significant from the basic-science viewpoint while an accurate theoretical description of melting of short duplexes is necessary for various very important applications in biotechnology. Several sets of DNA melting parameters, proposed within the framework of the nearest neighbor model, are compared and analyzed. The analysis leads to a conclusion that in case of long DNA molecules the consensus set of nearest neighbor parameters describes well the experimental melting profiles. Unexpectedly, for short DNA duplexes the same set of parameters hardly yields any improvement as compared to the simplest model, which completely ignores the effect of heterogeneous stacking. Possible causes of this striking observation are discussed. We then overview the issue of separation of base-pairing and base-stacking contributions into the double helix stability. The recent experimental attempts to solve the problem are extensively analyzed. It is concluded that the double helix is essentially stabilized by stacking interaction between adjacent base pairs. Base pairing between complementary pairs does not appreciably contribute into the duplex stability. In the final section of the article, kinetic aspects of the DNA melting phenomenon are discussed. The main emphasis is made on the hysteresis effects often observed in melting of long DNA molecules. It is argued that the phenomenon can be well described via an accurate theoretical treatment of the random-walk model of melting kinetics of an isolated helical segment in DNA.

Unique Thermal Stability of Unnatural Hydrophobic Ds Bases in Double-Stranded DNAs

Genetic alphabet expansion technology, the introduction of unnatural bases or base pairs into replicable DNA, has rapidly advanced as a new synthetic biology area. A hydrophobic unnatural base pair between 7-(2-thienyl)imidazo[4,5-b]pyridine (Ds) and 2-nitro-4-propynylpyrrole (Px) exhibited high fidelity as a third base pair in PCR. SELEX methods using the Ds-Px pair enabled high-affinity DNA aptamer generation, and introducing a few Ds bases into DNA aptamers extremely augmented their affinities and selectivities to target proteins. Here, to further scrutinize the functions of this highly hydrophobic Ds base, the thermal stabilities of double-stranded DNAs (dsDNA) containing a noncognate Ds-Ds or G-Ds pair were examined. The thermal stability of the Ds-Ds self-pair was as high as that of the natural G-C pair, and apart from the generally higher stability of the G-C pair than that of the A-T pair, most of the 5'-pyrimidine-Ds-purine-3' sequences, such as CDsA and TDsA, exhibited higher stability than the 5'-purine-Ds-pyrimidine-3' sequences, such as GDsC and ADsC, in dsDNAs. This trait enabled the GC-content-independent control of the thermal stability of the designed dsDNA fragments. The melting temperatures of dsDNA fragments containing the Ds-Ds pair can be predicted from the nearest-neighbor parameters including the Ds base. In addition, the noncognate G-Ds pair can efficiently distinguish its neighboring cognate natural base pairs from noncognate pairs. We demonstrated that real-time PCR using primers containing Ds accurately detected a single-nucleotide mismatch in target DNAs. These unique properties of the Ds base that affect the stabilities of the neighboring base pairs could impart new functions to DNA molecules and technologies.

Ionic effects on the temperature-force phase diagram of DNA

Double-stranded DNA (dsDNA) undergoes a structural transition to single-stranded DNA (ssDNA) in many biologically important processes such as replication and transcription. This strand separation arises in response either to thermal fluctuations or to external forces. The roles of ions are twofold, shortening the range of the interstrand potential and renormalizing the DNA elastic modulus. The dsDNA-to-ssDNA transition is studied on the basis that dsDNA is regarded as a bound state while ssDNA is regarded as an unbound state. The ground state energy of DNA is obtained by mapping the statistical mechanics problem to the imaginary time quantum mechanics problem. In the temperature-force phase diagram the critical force F _c (T) increases logarithmically with the Na⁺ concentration in the range from 32 to 110 mM. Discussing this logarithmic dependence of F _c (T) within the framework of polyelectrolyte theory, it inevitably suggests a constraint on the difference between the interstrand separation and the length per unit charge during the dsDNA-to-ssDNA transition.

Sufficient minimal model for DNA denaturation: Integration of harmonic scalar elasticity and bond energies

We study DNA denaturation by integrating elasticity - as described by the Gaussian network model - with bond binding energies, distinguishing between different base pairs and stacking energies. We use exact calculation, within the model, of the Helmholtz free-energy of any partial denaturation state, which implies that the entropy of all formed "bubbles" ("loops") is accounted for. Considering base pair bond removal single events, the bond designated for opening is chosen by minimizing the free-energy difference for the process, over all remaining base pair bonds. Despite of its great simplicity, for several known DNA sequences our results are in accord with available theoretical and experimental studies. Moreover, we report free-energy profiles along the denaturation pathway, which allow to detect stable or meta-stable partial denaturation states, composed of bubble, as local free-energy minima separated by barriers. Our approach allows to study very long DNA strands with commonly available computational power, as we demonstrate for a few random sequences in the range 200-800 base-pairs. For the latter, we also elucidate the self-averaging property of the system. Implications for the well known breathing dynamics of DNA are elucidated.

Competitive Assays of Label-Free DNA Hybridization with Single-Molecule Fluorescence Imaging Detection

Single-molecule imaging of fluorescently labeled biomolecules is a powerful technique for measuring association interactions; however, care must be taken to ensure that the fluorescent labels do not influence the system being probed. Label-free techniques are needed to understand biomolecule interactions free from the influence of an attached label, but these techniques often lack sensitivity and specificity. To solve these challenges, we have developed a competitive assay that uses single-molecule detection to track the population of unlabeled target single-stranded DNA (ssDNA) hybridized with probe DNA immobilized at a glass interface by detecting individual duplexes with a fluorescently labeled "tracer" ssDNA. By labeling a small fraction (<0.2%) of target molecules, the "tracer" DNA tracks the available probe DNA sites without significant competition with the unlabeled target population. Single-molecule fluorescence imaging is a good read-out scheme for competitive assays, as it is sufficiently sensitive to detect tracer DNA on substrates with relatively low densities of probe DNA, ∼10(-3) of a monolayer, so that steric interactions do not hinder DNA hybridization. Competitive assays are used to measure the association constant of complementary strand DNA hybridization of 9- and 10-base pair targets, where the tracer assay predicts the same association constant as a traditional displacement competitive assay. This methodology was used to compare the Ka of hybridization for identical DNA strands differing only by the presence of a fluorescent label tethered to the 5' end of the solution-phase target. The addition of the fluorescent label significantly stabilizes the DNA duplex by 3.6 kJmol(-1), adding more stability than an additional adenine-thymine base-pairing interaction, 2.7 kJmol(-1). This competitive tracer assay could be used to screen a number of labeled and unlabeled target DNA strands to measure the impact of fluorescent labeling on duplex stability. This single-molecule competitive hybridization scheme could be easily adapted into a sensitive assay, where competition between tracer and target oligonucleotides for probe sites could be used to measure concentrations of unlabeled DNA or RNA.

Designing the structure and folding pathway of modular topological bionanostructures

Polypeptides and polynucleotides are programmable natural polymers whose linear sequence can be easily designed and synthesized by the cellular transcription/translation machinery. Nature primarily uses proteins as the molecular machines and nucleic acids as the medium for the manipulation of heritable information. A protein's tertiary structure and function is defined by multiple cooperative weak long-range interactions that have been optimized through evolution. DNA nanotechnology uses orthogonal pairwise interacting modules of complementary nucleic acids as a strategy to construct defined complex 3D structures. A similar approach has recently been applied to protein design, using orthogonal dimerizing coiled-coil segments as interacting modules. When concatenated into a single polypeptide chain, they self-assemble into the 3D structure defined by the topology of interacting modules within the chain. This approach allows the construction of geometric polypeptide scaffolds, bypassing the folding problem of compact proteins by relying on decoupled pairwise interactions. However, the folding pathway still needs to be optimized in order to allow rapid self-assembly under physiological conditions. Again the modularity of designed topological structures can be used to define the rules that guide the folding pathway of long polymers, such as DNA, based on the stability and topology of connected building modules. This approach opens the way towards incorporation of designed foldamers in biological systems and their functionalization.

PRIMEGENSw3: a web-based tool for high-throughput primer and probe design

Highly specific and efficient primer and probe design has been a major hurdle in many high-throughput techniques. Successful implementation of any PCR or probe hybridization technique depends on the quality of primers and probes used in terms of their specificity and cross-hybridization. Here we describe PRIMEGENSw3, a set of web-based utilities for high-throughput primer and probe design. These utilities allow users to select genomic regions and to design primer/probe for selected regions in an interactive, user-friendly, and automatic fashion. The system runs the PRIMEGENS algorithm in the back-end on the high-performance server with the stored genomic database or user-provided custom database for cross-hybridization check. Cross-hybridization is checked not only using BLAST but also by checking mismatch positions and energy calculation of potential hybridization hits. The results can be visualized online and also can be downloaded. The average success rate of primer design using PRIMEGENSw3 is ~90 %. The web server also supports primer design for methylated sequences, which is used in epigenetic studies. Stand-alone version of the software is also available for download at the website.

Low-concentration initiator primers improve the amplification of gene targets with high sequence variability

The amplification and detection of diverse strains of an infectious virus or bacteria, or variants within a gene family is important for both clinical and basic research but can be difficult using conventional PCR. This report describes and illustrates a novel closed-tube method for amplifying and characterizing heterogeneous target sequences using members of the CTX-M beta-lactamase gene family. Different subgroups of CTX-M genes exhibit low sequence identity, but accurate and efficient detection of these variants is critical because they all confer resistance to penicillin, cefotaxime, and other antibiotics of the beta-lactam class. The method combines a single pair of "thermodynamic consensus primers" (tcPrimers) with one or more "initiator primers" (iPrimers), added at low concentration (5-10 nM). Each iPrimer improves the initial amplification of one or more variants because it has fewer mismatches to its intended target than the more abundant tcPrimers. As a result of initial amplification, each heterogeneous sequence is shifted stepwise toward a better match with the tcPrimers. As soon as the tcPrimer hybridization takes place, amplification proceeds with high efficiency. The tcPrimer pairs can be designed for symmetric PCR or for Linear-After-The-Exponential (LATE)-PCR. LATE-PCR offers the advantage of generating single-stranded DNA that can be characterized for different gene variants in the same closed tube, using low-temperature mismatch-tolerant fluorescent probes.

Rationally designed molecular beacons for bioanalytical and biomedical applications

Nucleic acids hold promise as biomolecules for future applications in biomedicine and biotechnology. Their well-defined structures and compositions afford unique chemical properties and biological functions. Moreover, the specificity of hydrogen-bonded Watson-Crick interactions allows the construction of nucleic acid sequences with multiple functions. In particular, the development of nucleic acid probes as essential molecular engineering tools will make a significant contribution to advancements in biosensing, bioimaging and therapy. The molecular beacon (MB), first conceptualized by Tyagi and Kramer in 1996, is an excellent example of a double-stranded nucleic acid (dsDNA) probe. Although inactive in the absence of a target, dsDNA probes can report the presence of a specific target through hybridization or a specific recognition-triggered change in conformation. MB probes are typically fluorescently labeled oligonucleotides that range from 25 to 35 nucleotides (nt) in length, and their structure can be divided into three components: stem, loop and reporter. The intrinsic merit of MBs depends on predictable design, reproducibility of synthesis, simplicity of modification, and built-in signal transduction. Using resonance energy transfer (RET) for signal transduction, MBs are further endowed with increased sensitivity, rapid response and universality, making them ideal for chemical sensing, environmental monitoring and biological imaging, in contrast to other nucleic acid probes. Furthermore, integrating MBs with targeting ligands or molecular drugs can substantially support their in vivo applications in theranositics. In this review, we survey advances in bioanalytical and biomedical applications of rationally designed MBs, as they have evolved through the collaborative efforts of many researchers. We first discuss improvements to the three components of MBs: stem, loop and reporter. The current applications of MBs in biosensing, bioimaging and therapy will then be described. In particular, we emphasize recent progress in constructing MB-based biosensors in homogeneous solution or on solid surfaces. We expect that such rationally designed and functionalized MBs will open up new and exciting avenues for biological and medical research and applications.

A statistical thermodynamic model for investigating the stability of DNA sequences from oligonucleotides to genomes

We describe the development and testing of a simple statistical mechanics methodology for duplex DNA applicable to sequences of any composition and extensible to genomes. The microstates of a DNA sequence are modeled in terms of blocks of basepairs that are assumed to be fully closed (paired) or open. This approach generates an ensemble of bubblelike microstates that are used to calculate the corresponding partition function. The energies of the microstates are calculated as additive contributions from hydrogen bonding, basepair stacking, and solvation terms parameterized from a comprehensive series of molecular dynamics simulations including solvent and ions. Thermodynamic properties and nucleotide stability constants for DNA sequences follow directly from the partition function. The methodology was tested by comparing computed free energies per basepair with the experimental melting temperatures of 60 oligonucleotides, yielding a correlation coefficient of -0.96. The thermodynamic stability of genic/nongenic regions was tested in terms of nucleotide stability constants versus sequence for the Escherichia coli K-12 genome. It showed clear differentiation of the genes from promoters and captures genic regions with a sensitivity of 0.94. The statistical thermodynamic model presented here provides a seemingly new handle on the challenging problem of interpreting genomic sequences.

Prediction of melting temperatures in fluorescence in situ hybridization (FISH) procedures using thermodynamic models

The thermodynamics and kinetics of DNA hybridization, i.e. the process of self-assembly of one, two or more complementary nucleic acid strands, has been studied for many years. The appearance of the nearest-neighbor model led to several theoretical and experimental papers on DNA thermodynamics that provide reasonably accurate thermodynamic information on nucleic acid duplexes and allow estimation of the melting temperature. Because there are no thermodynamic models specifically developed to predict the hybridization temperature of a probe used in a fluorescence in situ hybridization (FISH) procedure, the melting temperature is used as a reference, together with corrections for certain compounds that are used during FISH. However, the quantitative relation between melting and experimental FISH temperatures is poorly described. In this review, various models used to predict the melting temperature for rRNA targets, for DNA oligonucleotides and for nucleic acid mimics (chemically modified oligonucleotides), will be addressed in detail, together with a critical assessment of how this information should be used in FISH.

Cy3 and Cy5 dyes attached to oligonucleotide terminus stabilize DNA duplexes: predictive thermodynamic model

Cyanine dyes are important chemical modifications of oligonucleotides exhibiting intensive and stable fluorescence at visible light wavelengths. When Cy3 or Cy5 dye is attached to 5' end of a DNA duplex, the dye stacks on the terminal base pair and stabilizes the duplex. Using optical melting experiments, we have determined thermodynamic parameters that can predict the effects of the dyes on duplex stability quantitatively (ΔG°, Tm). Both Cy dyes enhance duplex formation by 1.2 kcal/mol on average, however, this Gibbs energy contribution is sequence-dependent. If the Cy5 is attached to a pyrimidine nucleotide of pyrimidine-purine base pair, the stabilization is larger compared to the attachment to a purine nucleotide. This is likely due to increased stacking interactions of the dye to the purine of the complementary strand. Dangling (unpaired) nucleotides at duplex terminus are also known to enhance duplex stability. Stabilization originated from the Cy dyes is significantly larger than the stabilization due to the presence of dangling nucleotides. If both the dangling base and Cy3 are present, their thermodynamic contributions are approximately additive. New thermodynamic parameters improve predictions of duplex folding, which will help design oligonucleotide sequences for biophysical, biological, engineering, and nanotechnology applications.

Multiplexed programmable release of captured DNA

Nucleic-acid hybridization is widely used for the specific capture of complementary sequences from complex samples. It is useful for both analytical methodologies, such as array hybridization (e.g. transcriptome analysis, genetic-variation analysis), and preparative strategies such as exome sequencing and sequence-specific proteome capture and analysis (PICh, HyCCAPP). It has not generally been possible to selectively elute particular captured subsequences, however, as the conditions employed for disruption of a duplex can lack the specificity needed to discriminate between different sequences. We show here that it is possible to bind and selectively release multiple sets of sequences by using toehold-mediated DNA branch migration. The strategy is illustrated for simple mixtures of oligonucleotides, for the sequence-specific capture and specific release of crosslinked yeast chromatin, and for the specific release of oligonucleotides hybridized to DNA microarrays.

Emerging technologies for hybridization based single nucleotide polymorphism detection

Detection of single nucleotide polymorphisms (SNPs) is a crucial challenge in the development of a novel generation of diagnostic tools. Accurate detection of SNPs can prove elusive, as the impact of a single variable nucleotide on the properties of a target sequence is limited, even if this sequence consists of only a few nucleotides. New, accurate and facile strategies for the detection of point mutations are therefore absolutely necessary for the increased adoption of point-of-care molecular diagnostics. Currently, PCR and sequencing are mostly applied for diagnosing SNPs. However these methods have serious drawbacks as routine diagnostic tools because of their labour intensity and cost. Several new, more suitable methods can be applied to enable sensitive detection of mutations based on specially designed hybridization probes, mutation recognizing enzymes and thermal denaturation. Here, an overview is presented of the most recent advances in the field of fast and sensitive SNP detection assays with strong potential for integration in point-of-care tests.

Application of nucleic acid-lipid conjugates for the programmable organisation of liposomal modules

We present a critical review of recent work related to the assembly of multicompartment liposome clusters using nucleic acids as a specific recognition unit to link liposomal modules. The asymmetry in nucleic acid binding to its non-self complementary strand allows the controlled association of different compartmental modules into composite systems. These biomimetic multicompartment architectures could have future applications in chemical process control, drug delivery and synthetic biology. We assess the different methods of anchoring DNA to lipid membrane surfaces and discuss how lipid and DNA properties can be tuned to control the morphology and properties of liposome superstructures. We consider different methods for chemical communication between the contents of liposomal compartments within these clusters and assess the progress towards making this chemical mixing efficient, switchable and chemically specific. Finally, given the current state of the art, we assess the outlook for future developments towards functional modular networks of liposomes.

DNA-functionalized polystyrene particles and their controlled self-assembly

DNA-grafted particles can assemble together, and the interaction strength is related to the hydrogen bonding from DNA linker between DNA-grafted particles. Two types of DNA linker are proposed.

Comparative studies on species identification of Noctuoidea moths in two nature reserve conservation zones (Beijing, China) using DNA barcodes and thin-film biosensor chips

Rapid and accurate identification of species is required for the biological control of pest Noctuoidea moths. DNA barcodes and thin-film biosensor chips are two molecular approaches that have gained wide attention. Here, we compare these two methods for the identification of a limited number of Noctuoidea moth species. Based on the commonly used mitochondrial gene cytochrome c oxidase I (the standard DNA barcode for animal species), 14 probes were designed and synthesized for 14 species shared by two national nature reserves in Beijing and Hebei, China. Probes ranged in length from 18 to 27 bp and were designed as mismatch probes to guarantee that there were at least three base differences between the probe and nontarget sequences. The results on the chip could be detected by the naked eye without needing special equipment. No cross-hybridizations were detected although we tested all probes on the 14 target and 24 nontarget Noctuoidea species. The neighbour-joining tree of the 38 species based on COI sequences gave 38 highly supported independent groups. Both DNA barcoding and thin-film biosensor chips, based on the COI gene, are able to accurately identify and discriminate the 14 targeted moth species in this study. Because of its speed, high accuracy and low cost, the thin-film biosensor chip is a very practical means of species identification. Now, a more comprehensive chip will be developed for the identification of additional Noctuoidea moths for pest control and ecological protection.

Isothermal hybridization kinetics of DNA assembly of two-dimensional DNA origami

The Watson-Crick base-pairing with specificity and predictability makes DNA molecules suitable for building versatile nanoscale structures and devices, and the DNA origami method enables researchers to incorporate more complexities into DNA-based devices. Thermally controlled atomic force microscopy in combination with nanomechanical spectroscopy with forces controlled in the pico Newton (pN) range as a novel technique is introduced to directly investigate the kinetics of multistrand DNA hybridization events on DNA origami nanopores under defined isothermal conditions. For the synthesis of DNA nanostructures under isothermal conditions at 60 °C, a higher hybridization rate, fewer defects, and a higher stability are achieved compared to room-temperature studies. By quantifying the assembly times for filling pores in origami structures at several constant temperatures, the fill factors show a consistent exponential increase over time. Furthermore, the local hybridization rate can be accelerated by adding a higher concentration of the staples. The new insight gained on the kinetics of staple-scaffold hybridization on the synthesis of two dimensional DNA origami structures may open up new routes and ideas for designing DNA assembly systems with increased potential for their application.

Melting of DNA oligomers: dynamical models and comparison with experimental results

We compare experimental melting curves of short heterogeneous DNA oligomers with theoretical curves derived from statistical mechanics. Partition functions are computed with the one-dimensional Peyrard-Bishop (PB) Hamiltonian, already used in the study of the melting of long DNA chains. Working with short chains we take into account, in the computations, not only the breaking of the interstrand hydrogen bonds, but also the complete dissociation of the double helix into separate single strands. Since this dissociation equilibrium is of general relevance, independent of the particular microscopic model, we give some details of its treatment. We discuss how the non bonded three-dimensional interactions, not explicitly considered in the one-dimensional PB model, are taken into account through the treatment of the dissociation equilibrium. We also evaluate the relevance of the dissociation as a function of the chain length.