Rational design of DNA sequences for nanotechnology, microarrays and molecular computers using Eulerian graphs

Nucleic acids are molecules of choice for both established and emerging nanoscale technologies. These technologies benefit from large functional densities of ‘DNA processing elements’ that can be readily manufactured. To achieve the desired functionality, polynucleotide sequences are currently designed by a process that involves tedious and laborious filtering of potential candidates against a series of requirements and parameters. Here, we present a complete novel methodology for the rapid rational design of large sets of DNA sequences. This method allows for the direct implementation of very complex and detailed requirements for the generated sequences, thus avoiding ‘brute force’ filtering. At the same time, these sequences have narrow distributions of melting temperatures. The molecular part of the design process can be done without computer assistance, using an efficient ‘human engineering’ approach by drawing a single blueprint graph that represents all generated sequences. Moreover, the method eliminates the necessity for extensive thermodynamic calculations. Melting temperature can be calculated only once (or not at all). In addition, the isostability of the sequences is independent of the selection of a particular set of thermodynamic parameters. Applications are presented for DNA sequence designs for microarrays, universal microarray zip sequences and electron transfer experiments.     

Deoxyribozymes that recode sequence information

Allosteric nucleic acid ligases have been used previously to transform analyte-binding into the formation of oligonucleotide templates that can be amplified and detected. We have engineered binary deoxyribozyme ligases whose two components are brought together by bridging oligonucleotide effectors. The engineered ligases can ‘read’ one sequence and then ‘write’ (by ligation) a separate, distinct sequence, which can in turn be uniquely amplified. The binary deoxyribozymes show great specificity, can discriminate against a small number of mutations in the effector, and can read and recode DNA information with high fidelity even in the presence of excess obscuring genomic DNA. In addition, the binary deoxyribozymes can read non-natural nucleotides and write natural sequence information. The binary deoxyribozyme ligases could potentially be used in a variety of applications, including the detection of single nucleotide polymorphisms in genomic DNA or the identification of short nucleic acids such as microRNAs.     

Molecular barcodes detect redundancy and contamination in hairpin-bisulfite PCR

PCR amplification of limited amounts of DNA template carries an increased risk of product redundancy and contamination. We use molecular barcoding to label each genomic DNA template with an individual sequence tag prior to PCR amplification. In addition, we include molecular ‘batch-stamps’ that effectively label each genomic template with a sample ID and analysis date. This highly sensitive method identifies redundant and contaminant sequences and serves as a reliable method for positive identification of desired sequences; we can therefore capture accurately the genomic template diversity in the sample analyzed. Although our application described here involves the use of hairpin-bisulfite PCR for amplification of double-stranded DNA, the method can readily be adapted to single-strand PCR. Useful applications will include analyses of limited template DNA for biomedical, ancient DNA and forensic purposes.     

Can Natural Proteins Designed with ‘Inverted’ Peptide Sequences Adopt Native-Like Protein Folds?

We have carried out a systematic computational analysis on a representative dataset of proteins of known three-dimensional structure, in order to evaluate whether it would possible to ‘swap’ certain short peptide sequences in naturally occurring proteins with their corresponding ‘inverted’ peptides and generate ‘artificial’ proteins that are predicted to retain native-like protein fold. The analysis of 3,967 representative proteins from the Protein Data Bank revealed 102,677 unique identical inverted peptide sequence pairs that vary in sequence length between 5–12 and 18 amino acid residues. Our analysis illustrates with examples that such ‘artificial’ proteins may be generated by identifying peptides with ‘similar structural environment’ and by using comparative protein modeling and validation studies. Our analysis suggests that natural proteins may be tolerant to accommodating such peptides.     

A versatile method for the preparation of conjugates of peptides with DNA/PNA/analog by employing chemo-selective click reaction in water

The specific 1,3 dipolar Hüisgen cycloaddition reaction known as ‘click-reaction’ between azide and alkyne groups is employed for the synthesis of peptide–oligonucleotide conjugates. The peptide nucleic acids (PNA)/DNA and peptides may be appended either by azide or alkyne groups. The cycloaddition reaction between the azide and alkyne appended substrates allows the synthesis of the desired conjugates in high purity and yields irrespective of the sequence and functional groups on either of the two substrates. The versatile approach could also be employed to generate the conjugates of peptides with thioacetamido nucleic acid (TANA) analog. The click reaction is catalyzed by Cu (I) in either water or in organic medium. In water, ~3-fold excess of the peptide-alkyne/azide drives the reaction to completion in 2 h with no side products.     

Amino acid and nucleotide recurrence in aligned sequences: synonymous substitution patterns in association with global and local base compositions

The tendency for repetitiveness of nucleotides in DNA sequences has been reported for a variety of organisms. We show that the tendency for repetitive use of amino acids is widespread and is observed even for segments conserved between human and Drosophila melanogaster at the level of >50% amino acid identity. This indicates that repetitiveness influences not only the weakly constrained segments but also those sequence segments conserved among phyla. Not only glutamine (Q) but also many of the 20 amino acids show a comparable level of repetitiveness. Repetitiveness in bases at codon position 3 is stronger for human than for D.melanogaster, whereas local repetitiveness in intron sequences is similar between the two organisms. While genes for immune system-specific proteins, but not ancient human genes (i.e. human homologs of Escherichia coli genes), have repetitiveness at codon bases 1 and 2, repetitiveness at codon base 3 for these groups is similar, suggesting that the human genome has at least two mechanisms generating local repetitiveness. Neither amino acid nor nucleotide repetitiveness is observed beyond the exon boundary, denying the possibility that such repetitiveness could mainly stem from natural selection on mRNA or protein sequences. Analyses of mammalian sequence alignments show that while the ‘between gene’ GC content heterogeneity, which is linked to ‘isochores’, is a principal factor associated with the bias in substitution patterns in human, ‘within gene’ heterogeneity in nucleotide composition is also associated with such bias on a more local scale. The relationship amongst the various types of repetitiveness is discussed.     

Use of locked nucleic acid oligonucleotides to add functionality to plasmid DNA

The available reagents for the attachment of functional moieties to plasmid DNA are limiting. Most reagents bind plasmid DNA in a non-sequence- specific manner, with undefined stoichiometry, and affect DNA charge and delivery properties or involve chemical modifications that abolish gene expression. The design and ability of oligonucleotides (ODNs) containing locked nucleic acids (LNAs) to bind supercoiled, double-stranded plasmid DNA in a sequence-specific manner are described for the first time. The main mechanism for LNA ODNs binding plasmid DNA is demonstrated to be by strand displacement. LNA ODNs are more stably bound to plasmid DNA than similar peptide nucleic acid (PNA) ‘clamps’ for procedures such as particle-mediated DNA delivery (gene gun). It is shown that LNA ODNs remain associated with plasmid DNA after cationic lipid-mediated transfection into mammalian cells. LNA ODNs can bind to DNA in a sequence-specific manner so that binding does not interfere with plasmid conformation or gene expression. Attachment of CpG-based immune adjuvants to plasmid by ‘hybrid’ phosphorothioate–LNA ODNs induces tumour necrosis factor-α production in the macrophage cell line RAW264.7. This observation exemplifies an important new, controllable methodology for adding functionality to plasmids for gene delivery and DNA vaccination.     

Localization and dynamics of small circular DNA in live mammalian nuclei

While genomic DNA, packaged into chromatin, is known to be locally constrained but highly dynamic in the nuclei of living cells, little is known about the localization and dynamics of small circular DNA molecules that invade cells by virus infection, application of gene therapy vectors or experimental transfection. To address this point, we have created traceable model substrates by direct labeling of plasmid DNA with fluorescent peptide nucleic acids, and have investigated their fate after microinjection into living cells. Here, we report that foreign DNA rapidly undergoes interactions with intranuclear structural sites that strongly reduce its mobility and restrict the DNA to regions excluding nucleoli and nuclear bodies such as PML bodies. The labeled plasmids partially co-localize with SAF-A, a well characterized marker protein for the nuclear ‘scaffold’ or ‘matrix’, and are resistant towards extraction by detergent and, in part, elevated salt concentrations. We show that the localization and the low mobility of plasmids is independent of the plasmid sequence, and does not require the presence of either a scaffold attachment region (SAR) DNA element or a functional promoter.     

Genome sequencing and analysis of two early-flowering cherry (Cerasus × kanzakura) varieties, ‘Kawazu-zakura’ and ‘Atami-zakura’

To gain genetic insights into the early-flowering phenotype of ornamental cherry, also known as sakura, we determined the genome sequences of two early-flowering cherry (Cerasus × kanzakura) varieties, ‘Kawazu-zakura’ and ‘Atami-zakura’. Because the two varieties are interspecific hybrids, likely derived from crosses between Cerasus campanulata (early-flowering species) and Cerasus speciosa, we employed the haplotype-resolved sequence assembly strategy. Genome sequence reads obtained from each variety by single-molecule real-time sequencing (SMRT) were split into two subsets, based on the genome sequence information of the two probable ancestors, and assembled to obtain haplotype-phased genome sequences. The resultant genome assembly of ‘Kawazu-zakura’ spanned 519.8 Mb with 1,544 contigs and an N50 value of 1,220.5 kb, while that of ‘Atami-zakura’ totalled 509.6 Mb with 2,180 contigs and an N50 value of 709.1 kb. A total of 72,702 and 69,528 potential protein-coding genes were predicted in the genome assemblies of ‘Kawazu-zakura’ and ‘Atami-zakura’, respectively. Gene clustering analysis identified 2,634 clusters uniquely presented in the C. campanulata haplotype sequences, which might contribute to its early-flowering phenotype. Genome sequences determined in this study provide fundamental information for elucidating the molecular and genetic mechanisms underlying the early-flowering phenotype of ornamental cherry tree varieties and their relatives.     

Sensitive isothermal detection of nucleic-acid sequence by primer generation–rolling circle amplification

A simple isothermal nucleic-acid amplification reaction, primer generation–rolling circle amplification (PG–RCA), was developed to detect specific nucleic-acid sequences of sample DNA. This amplification method is achievable at a constant temperature (e.g. 60°C) simply by mixing circular single-stranded DNA probe, DNA polymerase and nicking enzyme. Unlike conventional nucleic-acid amplification reactions such as polymerase chain reaction (PCR), this reaction does not require exogenous primers, which often cause primer dimerization or non-specific amplification. Instead, ‘primers’ are generated and accumulated during the reaction. The circular probe carries only two sequences: (i) a hybridization sequence to the sample DNA and (ii) a recognition sequence of the nicking enzyme. In PG–RCA, the circular probe first hybridizes with the sample DNA, and then a cascade reaction of linear rolling circle amplification and nicking reactions takes place. In contrast with conventional linear rolling circle amplification, the signal amplification is in an exponential mode since many copies of ‘primers’ are successively produced by multiple nicking reactions. Under the optimized condition, we obtained a remarkable sensitivity of 84.5 ymol (50.7 molecules) of synthetic sample DNA and 0.163 pg (~60 molecules) of genomic DNA from Listeria monocytogenes, indicating strong applicability of PG–RCA to various molecular diagnostic assays.     

Use of extremely short F?rster resonance energy transfer probes in real-time polymerase chain reaction

Described in the article is a new approach for the sequence-specific detection of nucleic acids in real-time polymerase chain reaction (PCR) using fluorescently labeled oligonucleotide probes. The method is based on the production of PCR amplicons, which fold into dumbbell-like secondary structures carrying a specially designed ‘probe-luring’ sequence at their 5′ ends. Hybridization of this sequence to a complementary ‘anchoring’ tail introduced at the 3′ end of a fluorescent probe enables the probe to bind to its target during PCR, and the subsequent probe cleavage results in the florescence signal. As it has been shown in the study, this amplicon-endorsed and guided formation of the probe-target duplex allows the use of extremely short oligonucleotide probes, up to tetranucleotides in length. In particular, the short length of the fluorescent probes makes possible the development of a ‘universal’ probe inventory that is relatively small in size but represents all possible sequence variations. The unparalleled cost-effectiveness of the inventory approach is discussed. Despite the short length of the probes, this new method, named Angler real-time PCR, remains highly sequence specific, and the results of the study indicate that it can be effectively used for quantitative PCR and the detection of polymorphic variations.     

Haloarcula hispanica CRISPR authenticates PAM of a target sequence to prime discriminative adaptation

The prokaryotic immune system CRISPR/Cas (Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR-associated genes) adapts to foreign invaders by acquiring their short deoxyribonucleic acid (DNA) fragments as spacers, which guide subsequent interference to foreign nucleic acids based on sequence matching. The adaptation mechanism avoiding acquiring ‘self’ DNA fragments is poorly understood. In Haloarcula hispanica, we previously showed that CRISPR adaptation requires being primed by a pre-existing spacer partially matching the invader DNA. Here, we further demonstrate that flanking a fully-matched target sequence, a functional PAM (protospacer adjacent motif) is still required to prime adaptation. Interestingly, interference utilizes only four PAM sequences, whereas adaptation-priming tolerates as many as 23 PAM sequences. This relaxed PAM selectivity explains how adaptation-priming maximizes its tolerance of PAM mutations (that escape interference) while avoiding mis-targeting the spacer DNA within CRISPR locus. We propose that the primed adaptation, which hitches and cooperates with the interference pathway, distinguishes target from non-target by CRISPR ribonucleic acid guidance and PAM recognition.     

Inhibition of MDR1 gene expression by chimeric HNA antisense oligonucleotides

Hexitol nucleic acids (HNAs) are nuclease resistant and provide strong hybridization to RNA. However, there is relatively little information on the biological properties of HNA antisense oligonucleotides. In this study, we compared the antisense effects of a chimeric HNA ‘gapmer’ oligonucleotide comprising a phosphorothioate central sequence flanked by 5′ and 3′ HNA sequences to conventional phosphorothioate oligonucleotides and to a 2′-O-methoxyethyl (2′-O-ME) phosphorothioate ‘gapmer’. The antisense oligomers each targeted a sequence bracketing the start codon of the message of MDR1, a gene involved in multi-drug resistance in cancer cells. Antisense and control oligonucleotides were delivered to MDR1-expressing cells using transfection with the cationic lipid Lipofectamine 2000. The anti-MDR1 HNA gapmer was substantially more potent than a phosphorothioate oligonucleotide of the same sequence in reducing expression of P-glycoprotein, the MDR1 gene product. HNA and 2′-O-ME gapmers displayed similar potency, but a pure HNA antisense oligonucleotide (lacking the phosphorothioate ‘gap’) was ineffective, indicating that RNase H activity was likely required. Treatment with anti-MDR1 HNA gapmer resulted in increased cellular accumulation of the drug surrogate Rhodamine 123 that correlated well with the reduced cell surface expression of P-glycoprotein. Thus, HNA gapmers may provide a valuable additional tool for antisense-based investigations and therapeutic approaches.     

Random multi-recombinant PCR for the construction of combinatorial protein libraries

Development of a new methodology to create protein libraries, which enable the exploration of global protein space, is an exciting challenge. In this study we have developed random multi-recombinant PCR (RM-PCR), which permits the shuffling of several DNA fragments without homologous sequences. In order to evaluate this methodology, we applied it to create two different combinatorial DNA libraries. For the construction of a ‘random shuffling library’, RM-PCR was used to shuffle six DNA fragments each encoding 25 amino acids; this affords many different fragment sequences whose every position has an equal probability to encode any of the six blocks. For the construction of the ‘alternative splicing library’, RM-PCR was used to perform different alternative splicings at the DNA level, which also yields different block sequences. DNA sequencing of the RM-PCR products in both libraries revealed that most of the sequences were quite different, and had a long open reading frame without a frame shift or stop codon. Furthermore, no distinct bias among blocks was observed. Here we describe how to use RM-PCR for the construction of combinatorial DNA libraries, which encode protein libraries that would be suitable for selection experiments in the global protein space.     

The crystal structure of an ‘All Locked’ nucleic acid duplex

‘Locked nucleic acids’ (LNAs) are known to introduce enhanced bio- and thermostability into natural nucleic acids rendering them powerful tools for diagnostic and therapeutic applications. We present the 1.9 Å X-ray structure of an ‘all LNA’ duplex containing exclusively modified β-d-2′-O-4′C-methylene ribofuranose nucleotides. The helix illustrates a new type of nucleic acid geometry that contributes to the understanding of the enhanced thermostability of LNA duplexes. A notable decrease of several local and overall helical parameters like twist, roll and propeller twist influence the structure of the LNA helix and result in a widening of the major groove, a decrease in helical winding and an enlarged helical pitch. A detailed structural comparison to the previously solved RNA crystal structure with the corresponding base pair sequence underlines the differences in conformation. The surrounding water network of the RNA and the LNA helix shows a similar hydration pattern.     

An in vitro selection scheme for oligonucleotide probes to discriminate between closely related DNA sequences

Using an in vitro selection, we have obtained oligonucleotide probes with high discriminatory power against multiple, similar nucleic acid sequences, which is often required in diagnostic applications for simultaneous testing of such sequences. We have tested this approach, referred to as iterative hybridizations, by selecting probes against six 22-nt-long sequence variants representing human papillomavirus, (HPV). We have obtained probes that efficiently discriminate between HPV types that differ by 3–7nt. The probes were found effective to recognize HPV sequences of the type 6, 11, 16, 18 and a pair of type 31 and 33, either when immobilized on a solid support or in a reverse configuration, as well to discriminate HPV types from the clinical samples. This methodology can be extended to generate diagnostic kits that rely on nucleic acid hybridization between closely related sequences. In this approach, instead of adjusting hybridization conditions to the intended set of probe–target pairs, we ‘adjust’, through in vitro selection, the probes to the conditions we have chosen. Importantly, these conditions have to be ‘relaxed’, allowing the formation of a variety of not fully complementary complexes from which those that efficiently recognize and discriminate intended from non-intended targets can be readily selected.     

Real-time investigation of nucleic acids phosphorylation process using molecular beacons

Phosphorylation of nucleic acids is an indispensable process to repair strand interruption of nucleic acids. We have studied the process of phosphorylation using molecular beacon (MB) DNA probes in real-time and with high selectivity. The MB employed in this method is devised to sense the product of a ‘phosphorylation–ligation’ coupled enzyme reaction. Compared with the current assays, this novel method is convenient, fast, selective, highly sensitive and capable of real-time monitoring in a homogenous solution. The preference of T4 polynucleotide kinase (T4 PNK) has been investigated using this approach. The results revealed that a single-stranded oligonucleotide containing guanine at the 5′ termini is most preferred, while those utilizing cytosine in this location are least preferred. The preference of (T)₉ was reduced greatly when phosphoryl was modified at the 5′ end, implying that T4 PNK could discern the phosphorylated/unphosphorylated oligonucleotides. The increase of oligonucleotide DNA length leads to an enhancement in preference. A fast and accurate method for assaying the kinase activity of T4 PNK has been developed with a wide linear detection range from 0.002 to 4.0 U/ml in 3 min. The effects of certain factors, such as NTP, ADP, (NH₄)₂SO₄ and Na₂HPO₄, on phosphorylation have been investigated. This novel approach enables us to investigate the interactions between proteins and nucleic acids in a homogenous solution, such as those found in DNA repair or in drug development.