Hachimoji DNA and RNA: A genetic system with eight building blocks

Reported here are DNA and RNA-like systems built from eight (hachi-) nucleotide letters (-moji) that form four orthogonal pairs. This synthetic genetic biopolymer meets the structural requirements needed to support Darwinism, including a polyelectrolyte backbone, predictable thermodynamic stability, and stereoregular building blocks that fit a Schrödinger aperiodic crystal. Measured thermodynamic parameters predict the stability of hachimoji duplexes, allowing hachimoji DNA to double the information density of natural terran DNA. Three crystal structures show that the synthetic building blocks do not perturb the aperiodic crystal seen in the DNA double helix. Hachimoji DNA was then transcribed to give hachimoji RNA in the form of a functioning fluorescent hachimoji aptamer. These results expand the scope of molecular structures that might support life, including life throughout the cosmos

Articles Nearest-Neighbor Thermodynamics and NMR of DNA Sequences with Internal A‚A, C‚C, G‚G, and T‚T Mismatches †

ABSTRACT: Thermodynamic measurements are reported for 51 DNA duplexes with A‚A, C‚C, G‚G, and T‚T single mismatches in all possible Watson-Crick contexts. These measurements were used to test the applicability of the nearest-neighbor model and to calculate the 16 unique nearest-neighbor parameters for the 4 single like with like base mismatches next to a Watson-Crick pair. The observed trend in stabilities of mismatches at 37°C is G‚G > T‚T ≈ A‚A > C‚C. The observed stability trend for the closing Watson-Crick pair on the 5′ side of the mismatch is G‚C g C‚G g A‚T g T‚A. The mismatch contribution to duplex stability ranges from -2.22 kcal/mol for GGC‚GGC to +2.66 kcal/mol for ACT‚ ACT. The mismatch nearest-neighbor parameters predict the measured thermodynamics with average deviations of ∆G°3 7 ) 3.3%, ∆H°) 7.4%, ∆S°) 8.1%, and T M ) 1.1°C. The imino proton region of 1-D NMR spectra shows that G‚G and T‚T mismatches form hydrogen-bonded structures that vary depending on the Watson-Crick context. The data reported here combined with our previous work provide for the first time a complete set of thermodynamic parameters for molecular recognition of DNA by DNA with or without single internal mismatches. The results are useful for primer design and understanding the mechanism of triplet repeat diseases. DNA mismatches occur in vivo due to misincorporation of bases during replication (1), heteroduplex formation during homologous recombination (2), mutagenic chemicals (3, 4), ionizing radiation (5), and spontaneous deamination (6). Knowledge of the thermodynamics of DNA mismatches will be useful for elucidating the mechanisms of polymerase fidelity and mismatch repair efficiency. Moreover, thermodynamic parameters for mismatch formation are important for DNA secondary structure prediction (see http://sun2.science.wayne.edu/∼jslsun2 and http://mfold1.wustl.edu/∼mfold/dna/form1.cgi). Recent work has shown that triplet repeat sequences form transiently stable hairpins that contain like with like base mismatche

Bubbles and denaturation in DNA

The local opening of DNA is an intriguing phenomenon from a statistical physics point of view, but is also essential for its biological function. For instance, the transcription and replication of our genetic code can not take place without the unwinding of the DNA double helix. Although these biological processes are driven by proteins, there might well be a relation between these biological openings and the spontaneous bubble formation due to thermal fluctuations. Mesoscopic models, like the Peyrard-Bishop-Dauxois model, have fairly accurately reproduced some experimental denaturation curves and the sharp phase transition in the thermodynamic limit. It is, hence, tempting to see whether these models could be used to predict the biological activity of DNA. In a previous study, we introduced a method that allows to obtain very accurate results on this subject, which showed that some previous claims in this direction, based on molecular dynamics studies, were premature. This could either imply that the present PBD should be improved or that biological activity can only be predicted in a more complex frame work that involves interactions with proteins and super helical stresses. In this article, we give detailed description of the statistical method introduced before. Moreover, for several DNA sequences, we give a thorough analysis of the bubble-statistics as function of position and bubble size and the so-called $l$-denaturation curves that can be measured experimentally. These show that some important experimental observations are missing in the present model. We discuss how the present model could be improved.Comment: 15 pages, 5 figures, published as Eur. Phys. J. E 20 : 421-434 AUG 200

Salerno's model of DNA reanalysed: could solitons have biological significance?

We investigate the sequence-dependent behaviour of localised excitations in a toy, nonlinear model of DNA base-pair opening originally proposed by Salerno. Specifically we ask whether ``breather'' solitons could play a role in the facilitated location of promoters by RNA polymerase. In an effective potential formalism, we find excellent correlation between potential minima and {\em Escherichia coli} promoter recognition sites in the T7 bacteriophage genome. Evidence for a similar relationship between phage promoters and downstream coding regions is found and alternative reasons for links between AT richness and transcriptionally-significant sites are discussed. Consideration of the soliton energy of translocation provides a novel dynamical picture of sliding: steep potential gradients correspond to deterministic motion, while ``flat'' regions, corresponding to homogeneous AT or GC content, are governed by random, thermal motion. Finally we demonstrate an interesting equivalence between planar, breather solitons and the helical motion of a sliding protein ``particle'' about a bent DNA axis.Comment: Latex file 20 pages, 5 figures. Manuscript of paper to appear in J. Biol. Phys., accepted 02/09/0

Accurate Estimates of Microarray Target Concentration from a Simple Sequence-Independent Langmuir Model

Background: Microarray technology is a commonly used tool for assessing global gene expression. Many models for estimation of target concentration based on observed microarray signal have been proposed, but, in general, these models have been complex and platform-dependent. Principal Findings: We introduce a universal Langmuir model for estimation of absolute target concentration from microarray experiments. We find that this sequence-independent model, characterized by only three free parameters, yields excellent predictions for four microarray platforms, including Affymetrix, Agilent, Illumina and a custom-printed microarray. The model also accurately predicts concentration for the MAQC data sets. This approach significantly reduces the computational complexity of quantitative target concentration estimates. Conclusions: Using a simple form of the Langmuir isotherm model, with a minimum of parameters and assumptions, and without explicit modeling of individual probe properties, we were able to recover absolute transcript concentrations with high R 2 on four different array platforms. The results obtained here suggest that with a ‘‘spiked-in’ ’ concentration serie

Application of Equilibrium Models of Solution Hybridization to Microarray Design and Analysis

Background: The probe percent bound value, calculated using multi-state equilibrium models of solution hybridization, is shown to be useful in understanding the hybridization behavior of microarray probes having 50 nucleotides, with and without mismatches. These longer oligonucleotides are in widespread use on microarrays, but there are few controlled studies of their interactions with mismatched targets compared to 25-mer based platforms. Principal Findings: 50-mer oligonucleotides with centrally placed single, double and triple mismatches were spotted on an array. Over a range of target concentrations it was possible to discriminate binding to perfect matches and mismatches, and the type of mismatch could be predicted accurately in the concentration midrange (100 pM to 200 pM) using solution hybridization modeling methods. These results have implications for microarray design, optimization and analysis methods. Conclusions: Our results highlight the importance of incorporating biophysical factors in both the design and the analysis of microarrays. Use of the probe ‘‘percent bound’ ’ value predicted by equilibrium models of hybridization is confirmed to be important for predicting and interpreting the behavior of long oligonucleotide arrays, as has been shown for shor

Development of a Flow-Trough Microarray based Reverse Transcriptase Multiplex Ligation-Dependent Probe Amplification Assay for the Detection of European Bunyaviruses

It is suspected that apart from tick-borne encephalitis virus several additional European Arboviruses such as the sandfly borne Toscana virus, sandfly fever Sicilian virus and sandfly fever Naples virus, mosquito-borne Tahyna virus, Inkoo virus, Batai virus and tick-borne Uukuniemi virus cause aseptic meningo-encephalitis or febrile disease in Europe. Currently, the microarray technology is developing rapidly and there are many efforts to apply it to infectious diseases diagnostics. In order to arrive at an assay system useful for high throughput analysis of samples from aseptic meningo-encephalitis cases the authors developed a combined multiplex ligation-dependent probe amplification and flow-through microarray assay for the detection of European Bunyaviruses. These results show that this combined assay indeed is highly sensitive, and specific for the accurate detection of multiple viruses

Large scale analysis of positional effects of single-base mismatches on microarray gene expression data

<p>Abstract</p> <p>Background</p> <p>Affymetrix GeneChips utilize 25-mer oligonucleotides probes linked to a silica surface to detect targets in solution. Mismatches due to single nucleotide polymorphisms (SNPs) can affect the hybridization between probes and targets. Previous research has indicated that binding between probes and targets strongly depends on the positions of these mismatches. However, there has been substantial variability in the effect of mismatch type across studies.</p> <p>Methods</p> <p>By taking advantage of naturally occurring mismatches between rhesus macaque transcripts and human probes from the Affymetrix U133 Plus 2 GeneChip, we collected the largest 25-mer probes dataset with single-base mismatches at each of the 25 positions on the probe ever used in this type of analysis.</p> <p>Results</p> <p>A mismatch at the center of a probe led to a greater loss in signal intensity than a mismatch at the ends of the probe, regardless of the mismatch type. There was a slight asymmetry between the ends of a probe: effects of mismatches at the 3' end of a probe were greater than those at the 5' end. A cross study comparison of the effect of mismatch types revealed that results were not in good agreement among different reports. However, if the mismatch types were consolidated to purine or pyrimidine mismatches, cross study conclusions could be generated.</p> <p>Conclusion</p> <p>The comprehensive assessment of the effects of single-base mismatches on microarrays provided in this report can be useful for improving future versions of microarray platform design and the corresponding data analysis algorithms.</p

Experimental free energy measurements of kinetic molecular states using fluctuation theorems

Recent advances in non-equilibrium statistical mechanics and single molecule technologies make it possible to extract free energy differences from irreversible work measurements in pulling experiments. To date, free energy recovery has been focused on native or equilibrium molecular states, whereas free energy measurements of kinetic states (i.e. finite lifetime states that are generated dynamically and are metastable) have remained unexplored. Kinetic states can play an important role in various domains of physics, such as nanotechnology or condensed matter physics. In biophysics, there are many examples where they determine the fate of molecular reactions: protein and peptide-nucleic acid binding, specific cation binding, antigen-antibody interactions, transient states in enzymatic reactions or the formation of transient intermediates and non-native structures in molecular folders. Here we demonstrate that it is possible to obtain free energies of kinetic states by applying extended fluctuation relations. This is shown by using optical tweezers to mechanically unfold and refold DNA structures exhibiting intermediate and misfolded kinetic states.Comment: main paper (16 pages, 5 figures) and supplementary information (22 pages, 14 figures

Controlling the temperature sensitivity of DNA-mediated colloidal interactions through competing linkages

We propose a new strategy to improve the self-assembly properties of DNA-functionalised colloids. The problem that we address is that DNA-functionalised colloids typically crystallize in a narrow temperature window, if at all. The underlying reason is the extreme sensitivity of DNA-mediated interactions to temperature or other physical control parameters. We propose to widen the window for colloidal crystallization by exploiting the competition between DNA linkages with different nucleotide sequences, which results in a temperature-dependent switching of the dominant bond type. Following such a strategy, we can decrease the temperature dependence of DNA-mediated self assembly to make systems that can crystallize in a wider temperature window than is possible with existing systems of DNA functionalised colloids. We report Monte Carlo simulations that show that the proposed strategy can indeed work in practice for real systems and specific, designable DNA sequences. Depending on the length ratio of the different DNA constructs, we find that the bond switching is either energetically driven (equal length or `symmetric' DNA) or controlled by a combinatorial entropy gain (`asymmetric' DNA), which results from the large number of possible binding partners for each DNA strand. We provide specific suggestions for the DNA sequences with which these effects can be achieved experimentally