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.

Computational Simulation of Adapter Length-Dependent LASSO Probe Capture Efficiency.

Multiplexed cloning of long DNA sequences is a valuable technique in many biotechnology applications, such as long-read genome sequencing and the creation of open reading frame (ORF) libraries. Long-adapter single-stranded oligonucleotide (LASSO) probes have shown promise as a tool to clone long DNA fragments. LASSO probes are molecular inversion probes (MIP) engineered with an adapter region of user-defined length, flanked between template-specific probe sequences. Herein, we demonstrate that the adapter length is a key feature of LASSO that influences the efficiency of gene capture and cloning. Furthermore, we applied a model based on Monte Carlo molecular simulation in order to study the relationship between the long-adapter length of LASSO and capture enrichment. Our results suggest that the adapter length is a factor that contributes to the free energy of target-probe interaction, thereby determining the efficiency of capture. The results indicate that LASSOs with extremely long adapters cannot capture the targets well. They also suggest that targets of different lengths may prefer adapters of different lengths.

Assembly and Assessment of DNA Scaffolded Vaccines.

Vaccines play an important role in preventing many life-threatening infectious diseases. To meet the demand of vaccination for treating a wide range of diseases, rational vaccine design has been recognized as a desirable and necessary strategy for development of safe and effective vaccines. DNA nanostructures are advantageous in the design and construction of synthetic vaccines, owing to their robust self-assembly, programmability, and precision control in complex organization, as well as their intrinsic adjuvant activity. Here, we describe a modular assembly of DNA scaffolded vaccine complex, composing of a model antigen, streptavidin, and adjuvant, CpG oligonucleotide. The DNA-assembled vaccines were found to elicit strong antigen-specific antibody responses, but causing little or no adverse reactions. Conceivably, this vaccine platform can be further optimized for improved immunogenicity and extended to the construction of various subunit vaccines.

TopDown real-time gene synthesis.

This chapter introduces a simple, cost-effective TopDown one-step gene synthesis method, which is suitable for the sequence assembly of fairly long DNA. This method can be distinguished from conventional gene synthesis methods by two key features: (1) the melting temperature of the outer primers is designed to be ∼8°C lower than that of the assembly oligonucleotides, and (2) different annealing temperatures are utilized to selectively control the efficiencies of oligonucleotide assembly and full-length template amplification. This method eliminates the interference between polymerase chain reactions (PCR) assembly and amplification in one-step gene synthesis. Additionally, the TopDown gene synthesis has been combined with the LCGreen I DNA fluorescence dye in a real-time gene synthesis approach for investigating the stepwise efficiency and kinetics of PCR-based gene synthesis. The obtained real-time fluorescence signals are compared with gel electrophoresis results to optimize gene synthesis conditions.

DNA synthesis security.

It is generally assumed that genetic engineering advances will, inevitably, facilitate the misapplication of biotechnology toward the production of biological weapons. Unexpectedly, however, some of these very advances in the areas of DNA synthesis and sequencing may enable the implementation of automated and nonintrusive safeguards to avert the illicit applications of biotechnology. In the case of DNA synthesis, automated DNA screening tools could be built into DNA synthesizers in order to block the synthesis of hazardous agents. In addition, a comprehensive safety and security regime for dual-use genetic engineering research could include nonintrusive monitoring of DNA sequencing. This is increasingly feasible as laboratories outsource this service to just a few centralized sequencing factories. The adoption of automated, nonintrusive monitoring and surveillance of the DNA synthesis and sequencing pipelines may avert many risks associated with dual-use biotechnology. Here, we describe the historical background and current challenges associated with dual-use biotechnologies and propose strategies to address these challenges.

Comparative assessment of fungal cellobiohydrolase I richness and composition in cDNA generated using oligo(dT) primers or random hexamers.

Understanding soil fungal distribution and activities, particularly at the level of gene expression, is important in unveiling mechanisms regulating their activities in situ. Recent identification of fungal genes involved in carbon cycling has provided the foundation for developing reverse-transcriptase PCR assays to monitor spatiotemporal gene expression patterns in soils and other complex microbial systems. The polyadenylated 3' ends of eukaryotic mRNA transcripts enables the use of oligo(dT) primers for cDNA synthesis, but this can result in the overrepresentation of the 3' end of transcripts in cDNA pools. In an effort to increase the uniformity of transcripts represented in cDNA pools, random hexamers have been used. The use of both priming methods is abundant in the literature, but we do not know how these methods perform relative to each other. We performed comparative richness and compositional analyses of the fungal glycosyl hydrolase family 7 cellobiohydrolase I gene cbhI amplified from soil cDNAs that had been generated using either oligo(dT) primers or random hexamers. Our results demonstrate that similar cbhI richness and composition were recovered using both approaches. Richness estimates and compositional profiles of cbhI sequence libraries generated from random hexamer-primed cDNA were more variable than from libraries generated from oligo(dT) primed cDNA. However, our overall results indicate that, on average, comparable richness and composition were recovered from soil cDNAs when either priming method was used.

Homogeneous insulin and C-Peptide sensors for rapid assessment of insulin and C-peptide secretion by the islets.

OBJECTIVE: Glucose-stimulated islet insulin or C-peptide secretion experiments are a fundamental tool for studying and assessing islet function. The goal of this work was to develop Ab-based fluorescent homogenous sensors that would allow rapid, inexpensive, near-instantaneous determinations of insulin and C-peptide levels in biological samples. RESEARCH DESIGN AND METHODS: Our approach was to use molecular pincer design (Heyduk et al., Anal Chem 2008;80:5152-5159) in which a pair of antibodies recognizing nonoverlapping epitopes of the target are modified with short fluorochrome-labeled complementary oligonucleotides and are used to generate a fluorescence energy transfer (FRET) signal in the presence of insulin or C-peptide. RESULTS: The sensors were capable of detecting insulin and C-peptide with high specificity and with picomolar concentration detection limits in times as short as 20 min. Insulin and C-peptide levels determined with the FRET sensors showed outstanding correlation with determinations performed under the same conditions with enzyme-linked immunosorbent assay. Most importantly, the sensors were capable of rapid and accurate determinations of insulin and C-peptide secreted from human or rodent islets, verifying their applicability for rapid assessment of islet function. CONCLUSIONS: The homogeneous, rapid, and uncomplicated nature of insulin and C-peptide FRET sensors allows rapid assessment of ß-cell function and could enable point-of-care determinations of insulin and C-peptide.

ImmunoAT method: An initial assessment for the detection of abnormal isoforms of prion protein in formalin-fixed and paraffin-embedded tissues.

The AT-tailing method is a labelling technique that utilises oligo(dA-dT)-dependent signal amplification. In this study, a new immunohistochemical application of the immunoAT method was developed. This method uses an oligo(dA-dT)-conjugated primary antibody (direct immunoAT method) or an oligo(dA-dT)-conjugated secondary antibody (indirect immunoAT method). Fifteen-base oligo(dA-dT)-conjugated antibodies (IgG-ATs) were prepared in advance by conjugating maleimide-activated oligo(dA-dT) to IgG via free sulfhydryl residues that had been introduced on the surface of IgG using Traut's reagent. Following the reaction with the target antigen and the IgG-AT, oligo(dA-dT) was elongated by DeltaTth DNA polymerase in the presence of dATP, dTTP and biotinylated dUTP, consequently labelling the antigen-antibody complex with a large amount of biotin. To initially evaluate the immunoAT method, the presence or absence of prion protein (PrP(sc)) was determined in formalin-fixed and paraffin-embedded sections of the medulla oblongata of cattle which had been under active surveillance for bovine spongiform encephalopathy. Sections were examined using direct and indirect immunoAT methods and the EnVision+ system (Dako) under conditions that were identical except for the differing IgG-AT and AT-tailing methods. PrP(sc) detection was consistent using all three methods. The clearest signals were obtained using the indirect immunoAT method, suggesting significant potential for this method.

Comment on "Remeasuring the double helix".

Mathew-Fenn et al. (Reports, 17 October 2008, p. 446) reported unexpected distance fluctuations in short end-labeled DNA constructs and interpreted them as evidence for cooperative DNA stretching modes. We show that when accounting for a subtle linker leverage effect, their data can be understood within standard noncooperative DNA elasticity.

Comparative assessment of plasmid and oligonucleotide DNA substrates in measurement of in vitro base excision repair activity.

Mammalian base excision repair (BER) is mediated through at least two subpathways designated 'single-nucleotide' (SN) and 'long-patch' (LP) BER (2-nucleotides long/more repair patch). Two forms of DNA substrate are generally used for in vitro BER assays: oligonucleotide- and plasmid-based. For plasmid-based BER assays, the availability of large quantities of substrate DNA with a specific lesion remains the limiting factor. Using sequence-specific endonucleases that cleave only one strand of DNA on a double-stranded DNA substrate, we prepared large quantities of plasmid DNA with a specific lesion. We compared the kinetic features of BER using plasmid and oligonucleotide substrates containing the same lesion and strategic restriction sites around the lesion. The K(m) for plasmid DNA substrate was slightly higher than that for the oligonucleotide substrate, while the V(max) of BER product formation for the plasmid and oligonucleotide substrates was similar. The catalytic efficiency of BER with the oligonucleotide substrate was slightly higher than that with the plasmid substrate. We conclude that there were no significant differences in the catalytic efficiency of in vitro BER measured with plasmid and oligonucleotide substrates. Analysis of the ratio of SN BER to LP BER was addressed using cellular extracts and a novel plasmid substrate.

An assessment of the role of DNA adenine methyltransferase on gene expression regulation in E coli.

N6-Adenine methylation is an important epigenetic signal, which regulates various processes, such as DNA replication and repair and transcription. In gamma-proteobacteria, Dam is a stand-alone enzyme that methylates GATC sites, which are non-randomly distributed in the genome. Some of these overlap with transcription factor binding sites. This work describes a global computational analysis of a published Dam knockout microarray alongside other publicly available data to throw insights into the extent to which Dam regulates transcription by interfering with protein binding. The results indicate that DNA methylation by DAM may not globally affect gene transcription by physically blocking access of transcription factors to binding sites. Down-regulation of Dam during stationary phase correlates with the activity of TFs whose binding sites are enriched for GATC sites.

NMR assessment of the global shape of a non-labelled DNA dodecamer containing a tandem of G-T mismatches.

We have carried out a solution study of two non-labelled self-complementary DNA dodecamers d(GACTGTACAGTC)2 and d(GACTGTGCAGTC)2 by NMR, the second sequence composed of two G-T mismatches. Structures were determined using distances extracted from NOE effects alone or using both NOE and RDC constraints, measured in three different liquid crystalline media. We ensured that our data on the influence of the mesogen on the DNA structures, and the way in which the RDCs were incorporated as constraints in the protocol refinement, were consistent. We also tested the influence of different sets of RDCs and the best means of optimizing the calculation of D(a) and R. Resolution and accuracy of the ten best energy final structures were compared. The addition of a small set of RDC constraints significantly improves the final determined structures. We took advantage of the specificity of the RDC, i.e. it contains orientational information, and explored the global shape of the DNA duplexes; it was found that the duplexes do not have a large curvature. For the G-T base pair, we observed, in this particular sequence (tandem of G-T mismatches), a new pattern of base pairing, which involved the formation of a bifurcated hydrogen bond.

High-resolution NMR structure of an AT-rich DNA sequence.

We have determined, by proton NMR and complete relaxation matrix methods, the high-resolution structure of a DNA oligonucleotide in solution with nine contiguous AT base pairs. The stretch of AT pairs, TAATTATAA x TTATAATTA, is imbedded in a 27-nucleotide stem-and-loop construct, which is stabilized by terminal GC base pairs and an extraordinarily stable DNA loop GAA (Hirao et al., 1994, Nucleic Acids Res. 22, 576-582). The AT-rich sequence has three repeated TAA x TTA motifs, one in the reverse orientation. Comparison of the local conformations of the three motifs shows that the sequence context has a minor effect here: atomic RMSD between the three TAA x TTA fragments is 0.4-0.5 A, while each fragment is defined within the RMSD of 0.3-0.4 A. The AT-rich stem also contains a consensus sequence for the Pribnow box, TATAAT. The TpA, ApT, and TpT x ApA steps have characteristic local conformations, a combination of which determines a unique sequence-dependent pattern of minor groove width variation. All three TpA steps are locally bent in the direction compressing the major groove of DNA. These bends, however, compensate each other, because of their relative position in the sequence, so that the overall helical axis is essentially straight.

Octamer-primed sequencing technology: effects of dNTP supplementation.

Octamer-primed sequencing is a directed DNA sequencing strategy that employs the use of a presynthesized octamer primer library. Together with electronic octamer sequencing technology (eOST) it provides a faster, less expensive way to obtain DNA sequence information and can be used as a valuable tool for gap closure in large-scale genomic sequencing. In this paper we discuss the effect of dGTP/TTP supplementation on octamer sequencing. We show that addition of 75 microM dGTP and 5 microM TTP can improve the sequencing success rate by increasing the length and accuracy of generated sequence information. We also discuss the effect of template base composition immediately downstream of the octamer primer on the outcome of octamer sequencing.

The 2-(N-formyl-N-methyl)aminoethyl group as a potential phosphate/thiophosphate protecting group in solid-phase oligodeoxyribonucleotide synthesis.

[reaction in text] The 2-(N-formyl-N-methyl)aminoethyl deoxyribonucleoside phosphoramidite 1 has been synthesized and used in the solid-phase synthesis of an octadecathymidylic acid as a cost-efficient monomer for potential application in the preparation of therapeutic oligonucleotides. The 2-(N-formyl-N-methyl)aminoethyl group can be cleaved from oligonucleotides according to a unique thermolytic cyclodeesterification process at pH 7.0. In addition to being cost-effective, the use of 1 simplifies oligonucleotide postsynthesis processing by eliminating the utilization of concentrated ammonium hydroxide in oligonucleotide deprotection.

Static curvature and flexibility measurements of DNA with microscopy. A simple renormalization method, its assessment by experiment and simulation.

We present the derivation of equations based on statistical polymer chain analysis and a method to quantify the average angle value of intrinsic bends and the local flexibility at a given locus on DNA fragments imaged by electron microscopy. DNA fragments of n base-pairs are considered as stiff chains of n jointed unit rigid rods. If the DNA fragments are composed of two branches A0Am and A0Bn, with, respectively, m and n base-pairs, where the standard deviations of the angle formed by two consecutive base-pairs are uniform over each branch, respectively, sigmathetaA and sigmathetaB, we show that the standard deviation of the angle AmA0Bn is: [formula: see text] where sigmatheta0 is the standard deviation of the angle at locus A0. This equation is established for small angular deviations by analysis of DNA at different scales and the validity of the methodology is controlled with the computation of the reduced chi2 statistical test. The length of the DNA fragments must be of the order of, or below, the persistence length, as determined by sets of statistics from computer simulations of DNA fragments. This is verified experimentally by a detailed analysis of the digitized contours of homogeneous linear 139 base-pair DNA fragments observed by electron microscopy. The images are compared to the reconstruction of DNA fragments from the measurements. The value found, sigma0=4.6 degrees/bp, is consistent with the well-accepted value for DNA in a plane. We discuss the relationship between the standard deviation of the measured angles and the flexibility at the base-pair level. This method is useful to quantify directly from microscopy techniques, such as electron or scanning force microscopy, the true bending angle, either intrinsic or induced by a ligand, and its associated flexibility at a given locus in any small DNA fragment.

Monte Carlo simulation of radiolytic attack to 5'-d[T4G4]4 sequence in a unimolecular quadruplex.

PURPOSE: To extend to a quadruplex the stochastic model of radiolytic attack previously applied to a quasi-random duplex DNA. MATERIALS AND METHODS: The quadruplex structure is obtained from the PDB databank (first structure from 201D entry). The probabilities of OH* radical attack at all sugar and base reaction sites are calculated using a stochastic model based on the Monte Carlo method. RESULTS: Good agreement between the calculated and experimental frank strand break (FSB) probabilities is obtained using the relative efficiencies of conversion of the C4' and C5'-centred radicals into FSB determined for the quasi-random duplex (2.8:1 respectively). Efficiencies of base radicals-to-alkali-revealed breaks (ARB) conversion are determined by fitting the calculated probabilities of base attack to the previously reported experimental probabilities for ARB. The efficiency of conversion of thymine radical into an ARB is 3.4+/-0.7 times higher than for guanine radical. CONCLUSIONS: This paper supports the calculation method and allows evaluation of the relative efficiencies of thymine and guanine radicals-to-ARB conversion.

Computational applications of DNA structural scales.

We study from a computational standpoint several different physical scales associated with structural features of DNA sequences, including dinucleotide scales such as base stacking energy and propeller twist, and trinucleotide scales such as bendability and nucleosome positioning. We show that these scales provide an alternative or complementary compact representation of DNA sequences. As an example we construct a strand invariant representation of DNA sequences. The scales can also be used to analyze and discover new DNA structural patterns, especially in combinations with hidden Markov models (HMMs). The scales are applied to HMMs of human promoter sequences revealing a number of significant differences between regions upstream and downstream of the transcriptional start point. Finally we show, with some qualifications, that such scales are by and large independent, and therefore complement each other.

Kinetic characterization of linear diffusion of the restriction endonuclease EcoRV on DNA.

We have examined the kinetic parameters of linear diffusion of EcoRV on DNA. The data were analyzed by Monte Carlo simulations in which the efficiency of recognition of EcoRV sites during linear diffusion, the efficiency of linear diffusion, and the behavior of enzymes at the ends of linear DNA is explicitly treated. The analysis of the dependence of linear diffusion on the concentrations of NaCl and MgCl2 shows that linear diffusion is maximal at 50 mM NaCl under all concentrations of MgCl2 tested and increases with increasing concentrations of Mg2+ up to 10 mM, the highest concentration used in the test. Under these conditions, EcoRV scans 2 x 10(6) bp during one binding event with a velocity of about 1.7 x 10(6) bp s-1. The enzyme tends to overlook cleavage sites at 1 mM but not at 10 mM MgCl2. This result confirms the thermodynamic finding that EcoRV does not bind very specifically to DNA in the absence of Mg2+. It demonstrates that there is a Mg2+-dependent continuous transition between a nonspecific and a specific binding mode of EcoRV to DNA. By comparing cleavage rates of linear DNA whose ends are free or blocked, we have shown that EcoRV has a very low probability to fall off at the ends of linear DNA. The enzyme rather is "reflected" and continues linear diffusion. EcoRV does not cleave oligonucleotides containing two EcoRV sites processively. Consequently, dissociation of the enzyme from the cleavage products is not preceded by a transfer to nonspecific DNA, and linear diffusion is not involved in product dissociation in EcoRV.

A computational approach to modeling nucleic acid hairpin structures.

Hairpin is a structural motif frequently observed in both RNA and DNA molecules. This motif is involved specifically in various biological functions (e.g., gene expression and regulation). To understand how these hairpin motifs perform their functions, it is important to study their structures. Compared to protein structural motifs, structures of nucleic acid hairpins are less known. Based on a set of reduced coordinates for describing nucleic acid structures and a sampling algorithm that equilibrates structures using Metropolis Monte Carlo simulation, we developed a method to model nucleic acid hairpin structures. This method was used to predict the structure of a DNA hairpin with a single-guanosine loop. The lowest energy structure from the ensemble of 200 sampled structures has a RMSD of < 1.5 A, from the structure determined using NMR. Additional constraints for the loop bases were introduced for modeling an RNA hairpin with two nucleotides in the loop. The modeled structure of this RNA hairpin has extensive base stacking and an extra hydrogen bond (between the CYT in the loop and a phosphate oxygen), as observed in the NMR structure.