Mismatches and bubbles in DNA

Effect of Internal and Bulge Loops on the Thermal Stability of Small DNA Duplexes

The thermal stabilities of DNA duplexes analogous to the let-7 microRNA:lin-41 mRNA complex from Caenorhabditis elegans have been measured by free solution capillary electrophoresis. DNA duplexes with the same stems but different types of internal or bulge loops and a control with no loop have also been studied. The melting temperatures of the DNA derivatives increased linearly with the logarithm of the Na+ or K+ ion concentration in the solution. Peaks in the electropherograms corresponding to duplexes with internal or bulge loops exhibited extensive tailing at high temperatures, suggesting that denaturation occurred by slow exchange between the duplexes and their component single strands. The single strands did not separate completely from the duplexes in aqueous solutions; instead, they appeared as small subpeaks on the tails of the duplex peaks. However, complete separation of the duplexes from their component single strands was observed at 20 °C in solutions containing 300 mM tetrapropylammonium ions. In addition, counterion condensation appears to be significantly reduced in DNA duplexes containing internal or bulge loops.

Disruption of energetic and dynamic base pairing cooperativity in DNA duplexes by an abasic site

DNA duplex stability arises from cooperative interactions between multiple adjacent nucleotides that favor base pairing and stacking when formed as a continuous stretch rather than individually. Lesions and nucleobase modifications alter this stability in complex manners that remain challenging to understand despite their centrality to biology. Here, we investigate how an abasic site destabilizes small DNA duplexes and reshapes base pairing dynamics and hybridization pathways using temperature-jump infrared spectroscopy and coarse-grained molecular dynamics simulations. We show how an abasic site splits the cooperativity in a short duplex into two segments, which destabilizes small duplexes as a whole and enables metastable half-dissociated configurations. Dynamically, it introduces an additional barrier to hybridization by constraining the hybridization mechanism to a step-wise process of nucleating and zipping a stretch on one side of the abasic site and then the other.

Bubble Relaxation Dynamics in Homopolymer DNA Sequences

Understanding the inherent timescales of large bubbles in DNA is critical to a thorough comprehension of its physicochemical characteristics, as well as their potential role on helix opening and biological function. In this work, we employ the coarse-grained Peyrard-Bishop-Dauxois model of DNA to study relaxation dynamics of large bubbles in homopolymer DNA, using simulations up to the microsecond time scale. By studying energy autocorrelation functions of relatively large bubbles inserted into thermalised DNA molecules, we extract characteristic relaxation times from the equilibration process for both adenine-thymine (AT) and guanine-cytosine (GC) homopolymers. Bubbles of different amplitudes and widths are investigated through extensive statistics and appropriate fittings of their relaxation. Characteristic relaxation times increase with bubble amplitude and width. We show that, within the model, relaxation times are two orders of magnitude longer in GC sequences than in AT sequences. Overall, our results confirm that large bubbles leave a lasting impact on the molecule's dynamics, for times between 0.5-500 ns depending on the homopolymer type and bubble shape, thus clearly affecting long-time evolutions of the molecule.

Enzymatic Weight Update Algorithm for DNA-Based Molecular Learning

Recent research in DNA nanotechnology has demonstrated that biological substrates can be used for computing at a molecular level. However, in vitro demonstrations of DNA computations use preprogrammed, rule-based methods which lack the adaptability that may be essential in developing molecular systems that function in dynamic environments. Here, we introduce an in vitro molecular algorithm that ‘learns’ molecular models from training data, opening the possibility of ‘machine learning’ in wet molecular systems. Our algorithm enables enzymatic weight update by targeting internal loop structures in DNA and ensemble learning, based on the hypernetwork model. This novel approach allows massively parallel processing of DNA with enzymes for specific structural selection for learning in an iterative manner. We also introduce an intuitive method of DNA data construction to dramatically reduce the number of unique DNA sequences needed to cover the large search space of feature sets. By combining molecular computing and machine learning the proposed algorithm makes a step closer to developing molecular computing technologies for future access to more intelligent molecular systems.

Direct observation of breathing dynamics at the mismatch induced DNA bubble with nanometre accuracy: a smFRET study

The detailed conformational dynamics of the melted region in double-stranded DNA has been studied using a combination of ensemble and single-molecule FRET techniques. We monitored the millisecond time scale fluctuation kinetics of the two strands at the bubble region that varies with the size of the bubble. As the individual strands at the melting bubble behave as single-stranded DNA, and hence fluctuate dynamically to attain energetically favored configurations, the rates of these fluctuations increase with increase in the bubble size. In different short DNAs under investigation, the two strands never cross each other to form a knot, irrespective of the number of base pair mismatches present. Rather, they prefer to stay apart from each other, as the size of the bubble increases and follow exactly an opposite trend for bubbles of smaller size. The range within which the bubble strands fluctuate are monitored with great accuracy in the nanometre resolution from the single-molecule FRET measurements. The shape of the bubble that plays a crucial role in determining the activity of the DNA was speculated. These results shall be useful in quantifying the chemical processes within DNA as well as to develop a deeper understanding of the activity of the DNA due to induced mismatches.

In vitro molecular machine learning algorithm via symmetric internal loops of DNA

Programmable biomolecules, such as DNA strands, deoxyribozymes, and restriction enzymes, have been used to solve computational problems, construct large-scale logic circuits, and program simple molecular games. Although studies have shown the potential of molecular computing, the capability of computational learning with DNA molecules, i.e., molecular machine learning, has yet to be experimentally verified. Here, we present a novel molecular learning in vitro model in which symmetric internal loops of double-stranded DNA are exploited to measure the differences between training instances, thus enabling the molecules to learn from small errors. The model was evaluated on a data set of twenty dialogue sentences obtained from the television shows Friends and Prison Break. The wet DNA-computing experiments confirmed that the molecular learning machine was able to generalize the dialogue patterns of each show and successfully identify the show from which the sentences originated. The molecular machine learning model described here opens the way for solving machine learning problems in computer science and biology using in vitro molecular computing with the data encoded in DNA molecules.

Dynamic Release of Bending Stress in Short dsDNA by Formation of a Kink and Forks

Bending with high curvature is one of the major mechanical properties of double-stranded DNA (dsDNA) that is essential for its biological functions. The emergence of a kink arising from local melting in the middle of dsDNA has been suggested as a mechanism of releasing the energy cost of bending. Herein, we report that strong bending induces two types of short dsDNA deformations, induced by two types of local melting, namely, a kink in the middle and forks at the ends, which we demonstrate using D-shaped DNA nanostructures. The two types of deformed dsDNA structures dynamically interconvert on a millisecond timescale. The transition from a fork to a kink is dominated by entropic contribution (anti-Arrhenius behavior), while the transition from a kink to a fork is dominated by enthalpic contributions. The presence of mismatches in dsDNA accelerates kink formation, and the transition from a kink to a fork is removed when the mismatch size is three base pairs.

Temperature effect on poly(dA).poly(dT): molecular dynamics simulation studies of polymeric and oligomeric constructs

Understanding unwinding and melting of double helical DNA is very important to characterize role of DNA in replication, transcription, translation etc. Sequence dependent melting thermodynamics is used extensively for detecting promoter regions but melting studies are generally done for short oligonucleotides. This study reports several molecular dynamics (MD) simulations of homopolymeric poly(dA).poly(dT) as regular oligonucleotide fragments as well as its corresponding polymeric constructs with water and charge-neutralizing counterions at different temperatures ranging from 300 to 400 K. We have eliminated the end-effect or terminal peeling propensity by employing MD simulation of DNA oligonucleotides in such a manner that gives rise to properties of polymeric DNA of infinite length. The dynamic properties such as basepairing and stacking geometry, groove width, backbone conformational parameters, bending, distribution of counter ions and number of hydrogen bonds of oligomeric and polymeric constructs of poly(dA).poly(dT) have been analyzed. The oligomer shows terminal fraying or peeling effect at temperatures above 340 K. The polymer shows partial melting at elevated temperatures although complete denaturations of basepairs do not take place. The analysis of cross strand hydrogen bonds shows that the number of N-H···O hydrogen bonds increases with increase in temperature while C-H···O hydrogen bond frequencies decrease with temperature. Restructuring of counterions in the minor groove with temperature appear as initiation of melting in duplex structures.

Unfolding kinetics of periodic DNA hairpins

DNA hairpin molecules with periodic base sequences can be expected to exhibit a regular coarse-grained free energy landscape (FEL) as a function of the number of open base pairs and applied mechanical force. Using a commonly employed model, we first analyze for which types of sequences a particularly simple landscape structure is predicted, where forward and backward energy barriers between partly unfolded states are decreasing linearly with force. Stochastic unfolding trajectories for such molecules with simple FEL are subsequently generated by kinetic Monte Carlo simulations. Introducing probabilities that can be sampled from these trajectories, it is shown how the parameters characterizing the FEL can be estimated. Already 300 trajectories, as typically generated in experiments, provide faithful results for the FEL parameters.

Effect of temperature on DNA double helix: An insight from molecular dynamics simulation

The three-dimensional structure of DNA contains various sequence-dependent structural information, which control many cellular processes in life, such as replication, transcription, DNA repair, etc. For the above functions, DNA double helices need to unwind or melt locally, which is different from terminal melting, as often seen in molecular dynamics (MD) simulations or even in many DNA crystal structures. We have carried out detailed MD simulations of DNA double helices of regular oligonucleotide fragments as well as in polymeric constructs with water and charge-neutralizing counter-ions at several different temperatures. We wanted to eliminate the end-effect or terminal melting propensity by employing MD simulation of DNA oligonucleotides in such a manner that gives rise to properties of polymeric DNA of infinite length. The polymeric construct is expected to allow us to see local melting at elevated temperatures. Comparative structural analysis of oligonucleotides and its corresponding virtual polymer at various temperatures ranging from 300 K to 400 K is discussed. The general behaviour, such as volume expansion coefficients of both the simulations show high similarity, indicating polymeric construct, does not give many artificial constraints. Local melting of a polymer, even at elevated temperature, may need a high nucleation energy that was not available in the short (7 ns) simulations. We expected to observe such nucleation followed by cooperative melting of the polymers in longer MD runs. Such simulations of different polymeric sequences would facilitate us to predict probable melting origins in a polymeric DNA.

Tandem GA residues on opposite sides of the loop in molecular beacon-like DNA hairpins compact the loop and increase hairpin stability

The free solution electrophoretic mobilities and thermal stabilities of hairpins formed by two complementary 26-nucleotide oligomers have been measured by capillary electrophoresis. The oligomers are predicted to form molecular beacon-like hairpins with 5 bp stems and 16 nucleotides in the loop. One hairpin, called hairpin2 (hp2), migrates with a relatively fast free solution mobility and exhibits melting temperatures that are reasonably well predicted by the popular structure-prediction program Mfold. Its complement, called hairpin1 (hp1), migrates with a slower free solution mobility and forms a stable hairpin only in solutions containing ≥200 mM Na(+). The melting temperatures observed for hp1 are ~18 °C lower than those observed for hp2 and ~20 °C lower than those predicted by Mfold. The greater thermal stability of hp2 is due to the presence of tandem GA residues on opposite sides of the loop. If the corresponding TC residues in the hp1 loop are replaced by tandem GA residues, the melting temperatures of the modified hairpin are close to those observed for hp2. Eliminating the tandem GA residues in the hp2 loop significantly decreases the thermal stability of hp2. If the loops are replaced by a loop of 16 thymine residues, the free solution mobilities and thermal stabilities of the T-loop hairpin are equal to those observed for hp1. Hence, the loop of hp1 appears to be relatively unstructured, with few base-base stacking interactions. Interactions between tandem GA residues on opposite sides of the hp2 loop appear to compact the loop and increase hairpin stability.

A unified Poland-Scheraga model of oligo- and polynucleotide DNA melting: salt effects and predictive power

Key biological and nano-technological processes require the partial or complete association and dissociation of complementary DNA strands. We present a variant of the Poland-Scheraga model for DNA melting where we introduce a local, sequence-dependent salt correction of the nearest-neighbor parameters. Furthermore, our formulation accounts for capping and interfacial energies of helical and coiled chain sections. We show that the model reproduces experimental data for melting temperatures over the full experimental range of strand length, strand concentration, and ionic strength of the solution. In particular, we reproduce a phenomenological relation by Frank-Kamenetskii for very long chains using a parameterization based on melting curves for short oligomers. However, we also show that the parameters of the Poland-Scheraga model are still not known with sufficient precision to quantitatively predict the fine structure of melting curves. This formulation of the Poland-Scheraga model opens the possibility to overcome this limitation by optimizing parameters with respect to an extended base of experimental data for short-, medium-, and long-chain melting. We argue that the often-discarded melting data for longer oligomers exhibiting non-two-state transitions could play a particularly important role.

Bubbles in DNA melting

This paper presents a review of largely our own work on the DNA melting transition, and some new measurements of the elastic energy of sharp bends in single-stranded DNA and RNA. The purpose is to present the point of view that studying the transition of intermediate size oligomers leads to valuable tests of the models, in particular the ingredients most important for a reduced-degrees-of-freedom description, such as the different role of base pairing and base stacking. We make the case that, with intermediate size oligomers, one can actually measure the bubble length, which exhibits a more interesting behavior than the fraction of dissociated bases alone. Here is where more work seems necessary, both on the experimental and the modeling side, to understand the differences between theory and experiments. We summarize our previous results on the cooperativity parameters, which suggest that the transition is never exactly two-state no matter how short the molecule, or in other words the nucleation size for bubbles opening at the ends of the molecule is essentially 1 base pair. We briefly discuss our own modification of the nearest-neighbor model which treats pairing and stacking separately, as a way to fit the experimental melting profiles in this intermediate length regime. Finally we go on to present some new measurements on the stability of DNA and RNA hairpins with very short loops.

Biophysical studies with parallel stranded oligoduplex

The thermal denaturation profiles of parallel stranded oligoduplexes during UV-spectroscopy and fluorescence spectroscopy, for the first time, were shown to follow the similar patterns as were described for antiparallel stranded oligoduplexes. The B-form structure of these center-bubble-forming oligoduplexes was determined with circular dichroism. Thermodynamic properties were also discussed.

Biophysical characterization of double-stranded oligonucleotides using ETBR and isothermal fluorescence spectroscopy: implication for SNP genotyping

The UV-absorption, fluorescence and CD spectra of aps 23 bp oligoduplexes were performed for potential diagnostic purpose. These oligonucleotide sequences were mimicked from natural mutations (mitochondrial genome) of human population (unpublished). This work was designed on the basis of hybridization of non-self complementary oligoduplexes (aps) containing no mismatch, one-mismatch and two-mismatches. Since melting temperature is dependent on concentration of the oligoduplex, various concentrations were used in this study protocol. The thermal spectra profiles (UV absorbance and fluorescence) of these oligoduplexes (aps) are different for a particular concentration, and can be implicated for mutations. -dF/dT (or dA/dT) vs T, lnK (or RlnK) vs TM, DeltaG vs TM, DeltaS vs TM and DeltaH vs TM are also variable for those sequences. All these thermodynamic data were calculated from absorbance (at 260 nm) data. On the contrary to the 23 bp oligoduplexes (aps), the PCR products of 97 bp and 256 bp length were genotyped with ETBR (excitation 530 nm, emission 600 nm) fluorimetrically. But our attempts to genotype these PCR sequences with isothermal UV absorbance spectroscopy were unsuccessful. Isothermal UV absorbance spectra has a limitation of sequence length. However, the structural conformation (all B-type) of the oligoduplexes (aps) was determined using CD. The minor discrepancy in CD spectra of these oligoduplexes are not significant for mutational analysis. 97 bp nested PCR product was an amplicon having either GcT or AcC mutation of mitochondria of normal human population, whereas 256 bp PCR product was an amplicon of human BRCA2 gene (NCBI Accession No. AY151039) of chromosome 13 having either A or G mutation at position -26.