Kilobase Lengths Research Articles

The N-terminal domain of human mitochondrial helicase Twinkle has DNA-binding activity crucial for supporting processive DNA synthesis by polymerase γ

Twinkle is the ring-shaped replicative helicase within the human mitochondria with high homology to bacteriophage T7 gp4 helicase-primase. Unlike many orthologs of Twinkle, the N-terminal domain (NTD) of human Twinkle has lost its primase activity through evolutionarily acquired mutations. The NTD has no demonstrated activity thus far; its role has remained unclear. Here, we biochemically characterize the isolated NTD and C-terminal domain (CTD) with linker to decipher their contributions to full-length Twinkle activities. This novel CTD construct hydrolyzes ATP, has weak DNA unwinding activity, and assists DNA polymerase γ (Polγ)-catalyzed strand-displacement synthesis on short replication forks. However, CTD fails to promote multikilobase length product formation by Polγ in rolling-circle DNA synthesis. Thus, CTD retains all the motor functions but struggles to implement them for processive translocation. We show that NTD has DNA-binding activity, and its presence stabilizes Twinkle oligomerization. CTD oligomerizes on its own, but the loss of NTD results in heterogeneously sized oligomeric species. The CTD also exhibits weaker and salt-sensitive DNA binding compared with full-length Twinkle. Based on these results, we propose that NTD directly contributes to DNA binding and holds the DNA in place behind the central channel of the CTD like a "doorstop," preventing helicase slippages and sustaining processive unwinding. Consistent with this model, mitochondrial single-stranded DNA-binding protein (mtSSB) compensate for the NTD loss and partially restore kilobase length DNA synthesis by CTD and Polγ. The implications of our studies are foundational for understanding the mechanisms of disease-causing Twinkle mutants that lie in the NTD.

The use of resolvases T4 endonuclease VII and T7 endonuclease I in mutation detection.

Mutation and polymorphism detection is of increasing importance in the field of molecular genetics. This is reflected by the plethora of chemical, enzymatic, and physically based methods of mutation detection. The ideal method would detect mutations in large fragments of DNA and position them to single base-pair (bp) accuracy. Few methods are able to quickly screen kilobase lengths of DNA and position the mutation at the same time. The Enzyme Mismatch Cleavage (EMC) method of mutation detection is able to reliably detect nearly 100% of mutations in DNA fragments as large as 2 kb and position them to within 6 bp. This method exploits the activity of a resolvase enzyme from T4, T4 endonuclease VII, and, more recently, a second bacteriophage resolvase, T7 endonuclease I. The technique uses these enzymes to digest heteroduplex DNA formed by annealing wild-type and mutant DNA. Digestion fragments indicate the presence, and the position, of any mutations. The method is robust and reliable and much faster and cheaper than sequencing. These attributes have resulted in its increasing use in the field of mutation detection.

SIST: stress-induced structural transitions in superhelical DNA.

Supercoiling imposes stress on a DNA molecule that can drive susceptible sequences into alternative non-B form structures. This phenomenon occurs frequently in vivo and has been implicated in biological processes, such as replication, transcription, recombination and translocation. SIST is a software package that analyzes sequence-dependent structural transitions in kilobase length superhelical DNA molecules. The numerical algorithms in SIST are based on a statistical mechanical model that calculates the equilibrium probability of transition for each base pair in the domain. They are extensions of the original stress-induced duplex destabilization (SIDD) method, which analyzes stress-driven DNA strand separation. SIST also includes algorithms to analyze B-Z transitions and cruciform extrusion. The SIST pipeline has an option to use the DZCBtrans algorithm, which analyzes the competition among these three transitions within a superhelical domain. The package and additional documentation are freely available at https://bitbucket.org/benhamlab/sist_codes. dzhabinskaya@ucdavis.edu.

Theoretical Analysis of Competing Conformational Transitions in Superhelical DNA

We develop a statistical mechanical model to analyze the competitive behavior of transitions to multiple alternate conformations in a negatively supercoiled DNA molecule of kilobase length and specified base sequence. Since DNA superhelicity topologically couples together the transition behaviors of all base pairs, a unified model is required to analyze all the transitions to which the DNA sequence is susceptible. Here we present a first model of this type. Our numerical approach generalizes the strategy of previously developed algorithms, which studied superhelical transitions to a single alternate conformation. We apply our multi-state model to study the competition between strand separation and B-Z transitions in superhelical DNA. We show this competition to be highly sensitive to temperature and to the imposed level of supercoiling. Comparison of our results with experimental data shows that, when the energetics appropriate to the experimental conditions are used, the competition between these two transitions is accurately captured by our algorithm. We analyze the superhelical competition between B-Z transitions and denaturation around the c-myc oncogene, where both transitions are known to occur when this gene is transcribing. We apply our model to explore the correlation between stress-induced transitions and transcriptional activity in various organisms. In higher eukaryotes we find a strong enhancement of Z-forming regions immediately 5′ to their transcription start sites (TSS), and a depletion of strand separating sites in a broad region around the TSS. The opposite patterns occur around transcript end locations. We also show that susceptibility to each type of transition is different in eukaryotes and prokaryotes. By analyzing a set of untranscribed pseudogenes we show that the Z-susceptibility just downstream of the TSS is not preserved, suggesting it may be under selection pressure.

Multivariable Parameter Optimization Of Microfluidic Post Arrays For DNA Fractionation

Length-based fractionation of long DNA is a fundamental process in genomic analysis. Traditional methods, such as pulsed-field gel electrophoresis or capillary electrophoresis, are slow and not easily coupled to downstream analytical processes. Microfluidic devices containing arrays of micron-scale posts have previously been described for length-based fractionation of kilobase length DNA. [1] Here, we present experimental and computational studies defining and optimizing a broad array of parameters in order to maximize mass throughput and isolation of DNA > 150 kb. Parameters explored include post field geometry, electric field intensity, and field oscillation timing. Experimental studies were performed by analyzing of the mobility of single YOYO-1 intercalated DNA molecules in PDMS post arrays. Computational studies simulated kilobase length DNA as a worm-like chain model and investigated the interaction of such molecules with post arrays of various geometries. Optimization of our system results in the ability to isolate 165 kb DNA from 125 kb DNA and process 5 ng of DNA within 30 minutes of operation. This research was supported by the Department of Homeland Security, Science and Technology Directorate.[1] Huang et al.., Nature Biotechnology, 2002, 20, 1048.

The potential and challenges of nanopore sequencing

A nanopore-based device provides single-molecule detection and analytical capabilities that are achieved by electrophoretically driving molecules in solution through a nano-scale pore. The nanopore provides a highly confined space within which single nucleic acid polymers can be analyzed at high throughput by one of a variety of means, and the perfect processivity that can be enforced in a narrow pore ensures that the native order of the nucleobases in a polynucleotide is reflected in the sequence of signals that is detected. Kilobase length polymers (single-stranded genomic DNA or RNA) or small molecules (e.g., nucleosides) can be identified and characterized without amplification or labeling, a unique analytical capability that makes inexpensive, rapid DNA sequencing a possibility. Further research and development to overcome current challenges to nanopore identification of each successive nucleotide in a DNA strand offers the prospect of 'third generation' instruments that will sequence a diploid mammalian genome for approximately $1,000 in approximately 24 h.

Factors Influencing Dna Expansion in the Course of Polymerase Chain Reaction

We performed more than 3,500 polymerase chain reactions (PCRs) under various conditions with more than 400 DNA fragments of 4–150 nucleotides in length. Some of the PCRs provided expanded DNA molecules of kilobase lengths whereas others led to no expansion. Repetitiveness of the primary structure was mostly found to be necessary but not sufficient for the expansion. (A+T)-rich fragments expand better than (G+C)-rich ones and pyrimidine-rich fragments expand better than purine-rich fragments. Terminal nucleotides and the fragment length also are important for the expansion. Examples are presented when relatively small alterations of the DNA primary structure caused a dramatic change in the expansion. For example, A8T8 expanded a lot whereas T8A8 did not expand at all. The present work has implications for pathological expansions of microsatellites in the human genome as well as regarding the genome evolution in general.

Expansion during PCR of short single-stranded DNA fragments carrying nonselfcomplementary dinucleotide or trinucleotide repeats.

We performed PCR of many DNA fragments of 6-32 nucleotides in length. Some of the fragments expanded into kilobase lengths even in the absence of the complementary strand. The dramatic expansion was observed for (CA)8, (TG)8, (CA)4, (CA)6, (CA)12, (TG)4, (TG)6, (TG)12, (TC)10, (GA)10 and other single strands. Similar expansions were exhibited by related trinucleotide repeats (TTG)5, (CAA)5, (TGG)5, and (CCA)5 as well. However even small perturbations of the strict repetitive nature of the DNA primary structure substantially reduced the expansions. The expansion products had properties characteristic for normal Watson-Crick duplexes. Hence either the Taq polymerase and/or other components of the PCR buffer promote homoduplex formation of the nonselfcomplementary fragments, which is necessary to prime the synthesis of the complementary DNA strand, or the Taq polymerase is able to copy the single-stranded DNA template without any priming effect. The present observations have implications for the evolution of genomic DNA, microsatellite length polymorphism as well as the pathological expansions of trinucleotide repeats in the human genome.

The use of resolvases T4 endonuclease VII and T7 endonuclease I in mutation detection.

Mutation and polymorphism detection is of increasing importance in the field of molecular genetics. This is reflected by the plethora of chemical, enzymatic, and physically based methods of mutation detection. The ideal method would detect mutations in large fragments of DNA and position them to single base-pair (bp) accuracy. Few methods are able to quickly screen kilobase lengths of DNA and position the mutation at the same time. The Enzyme Mismatch Cleavage (EMC) method of mutation detection is able to reliably detect nearly 100% of mutations in DNA fragments as large as 2 kb and position them to within 6 bp. This method exploits the activity of a resolvase enzyme from T4, T4 endonuclease VII, and, more recently, a second bacteriophage resolvase, T7 endonuclease I. The technique uses these enzymes to digest heteroduplex DNA formed by annealing wild-type and mutant DNA. Digestion fragments indicate the presence, and the position, of any mutations. The method is robust and reliable and much faster and cheaper than sequencing. These attributes have resulted in its increasing use in the field of mutation detection.

Sieving matrix selection.

The charge density of nucleic acids is uniform and constant for nucleic acid molecules of varying length. Therefore, they typically migrate at the same velocity in response to an electrical potential in free solution. Capillary electrophoretic (CE) separation of nucleic acid species requires addition of a sieving matrix to the background electrolyte to retard the analytes in proportion to their molecular size (). This can be achieved by casting a crosslinked gel in the capillary, or by incorporating a linear gel material or polymer in the background electrolyte. Both of these approaches have been applied successfully to the separation of nucleic acid species ranging from short, single-stranded oligonucleotides to dsDNA of several kilobase lengths. The selection of the appropriate sieving matrix will depend upon the size of the molecules that are to be separated, the resolution required, and the capabilities of the analytical instrumentation. This chapter will review the various sieving matrices which have been employed for CE separation of nucleic acids, identify the appropriate matrix for specific applications, and provide practical suggestions for preparation and use of sieving media for CE.

Alpha-fetoprotein structure and function: relevance to isoforms, epitopes, and conformational variants.

Mammalian alpha-fetoprotein (AFP) is classified as a member of the albuminoid gene superfamily consisting of albumin, AFP, vitamin D (Gc) protein, and alpha-albumin. Molecular variants of AFP have long been reported in the biomedical literature. Early studies identified isoelectric pH isoforms and lectin-binding variants of AFP, which differed in their physicochemical properties, but not in amino acid composition. Genetic variants of AFP, differing in mRNA kilobase length, were later extensively described in rodent models during fetal/perinatal stages, carcinogenesis, and organ regeneration. With the advent of monoclonal antibodies in the early 1980s, multiple antigenic epitopes on native AFP were detected and categorized, culminating in the identification of six to seven major epitopes. During this period, various AFP-binding proteins and receptors were reported to inhibit certain AFP immunoreactions. Concomittantly, human and rodent AFP were cloned and the amino acid sequences of the translated proteins were divulged. Once the amino acid composition of the AFP molecule was known, enzymatic fragments could be identified and synthetic peptide segments synthesized. Following discovery of the molten globule form in 1981, the existence of transitory, intermediate forms of AFP were acknowledged and their physiological significance was realized. In the present review, the various isoforms and variants of AFP are discussed in light of their potential biological relevance.

Exact method for numerically analyzing a model of local denaturation in superhelically stressed DNA

Local denaturation, the separation at specific sites of the two strands comprising the DNA double helix, is one of the most fundamental processes in biology, required to allow the base sequence to be read both in DNA transcription and in replication. In living organisms this process can be mediated by enzymes which regulate the amount of superhelical stress imposed on the DNA. We present a numerically exact technique for analyzing a model of denaturation in superhelically stressed DNA. This approach is capable of predicting the locations and extents of transition in circular superhelical DNA molecules of kilobase lengths and specified base pair sequences. It can also be used for closed loops of DNA which are typically found in vivo to be kilobases long. The analytic method consists of an integration over the DNA twist degrees of freedom followed by the introduction of auxiliary variables to decouple the remaining degrees of freedom, which allows the use of the transfer matrix method. The algorithm implementing our technique requires ${O(N}^{2})$ operations and $O(N)$ memory to analyze a DNA domain containing N base pairs. However, to analyze kilobase length DNA molecules it must be implemented in high precision floating point arithmetic. An accelerated algorithm is constructed by imposing an upper bound M on the number of base pairs that can simultaneously denature in a state. This accelerated algorithm requires $O(\mathrm{MN})$ operations, and has an analytically bounded error. Sample calculations show that it achieves high accuracy (greater than 15 decimal digits) with relatively small values of M $(M<0.05N)$ for kilobase length molecules under physiologically relevant conditions. Calculations are performed on the superhelical $\mathrm{pBR}322$ DNA sequence to test the accuracy of the method. With no free parameters in the model, the locations and extents of local denaturation predicted by this analysis are in quantitatively precise agreement with in vitro experimental measurements. Calculations performed on the fructose-1,6-bisphosphatase gene sequence from yeast show that this approach can also accurately treat in vivo denaturation.

Released chromatin: linearized DNA for high resolution fluorescence in situ hybridization.

Free DNA was prepared from routinely harvested and fixed cells for high resolution FISH mapping using either a sodium hydroxide/ethanol mixture or 70% formamide. Hybridization signals from cosmid probes appeared as extended lines. The average length of signals on DNA prepared with sodium hydroxide was significantly greater than with formamide. A set of overlapping cosmids from the HLA class II region was used to determine how precisely the actual overlap or gap between probes can be calculated from the measured overlap or gap between their signals. Lengths of the probe signals and their known kilobase lengths were used as an internal ruler. The mean values calculated from the measured length from 30 or more signals for each probe pair showed remarkable conformity with the known kilobase lengths of their overlaps and gaps. Immediately adjacent probes could also be ordered on the released DNA. These simple procedures dramatically increase the speed with which relationships between probes can be determined during contig construction.

Theoretical analysis of heteropolymeric transitions in superhelical DNA molecules of specified sequence

This paper develops accurate and computationally tractable, theoretical methods to analyze equilibrium heteropolymeric transitions in superhelical DNAs of specified sequence and kilo-base length. Although these methods are generally applicable, their development here is focused primarily on superhelical strand separation, a heteropolymeric transition of great biological importance to which every base pair in a duplex DNA molecule is susceptible. Because the total number of states of this transition grows exponentially with molecular length, exact analytic methods that consider all possible states are not feasible for DNAs of practical interest. First, an approximate statistical mechanical analysis is developed in which all states are considered whose free energies do not exceed that of the minimum energy state by more than a specified threshold amount. An approximate partition function is constructed using these states, from which estimates of the equilibrium values of important transition parameters are calculated. Next, a density of states analysis is performed to estimate the influence of the neglected, high energy states. Its results are used to refine the computed equilibrium parameter values, correcting for the effects of the neglected states. Only the transition profile (the probability of transition of each base pair in the sequence) cannot be refined in this manner, although its overall accuracy can be assessed. Sample calculations analyzing the strand separation transition in superhelical pBR322 DNA show that the analytic methods developed here yield accurate results in a computationally feasible manner. Moreover, the predictions of this analysis agree closely with experimental results.

Sequence organization of the soybean genome

The total complexity of one constituent soybean ( Glycine max) genome is estimated to be 1.29 · 10 9 nucleotide pairs, as determined by analysis of the reassociation kinetics of sheared (0.47 kilobase) DNA. Single copy sequences are estimated to represent from 53 to 64% of the genome by analysis of hydroxyapatite binding of repetitive DNA as a function of fragment length. From 65 to 70% of these sungle copy sequences have a short period interspetion with 1.11–1.36 kilobase lengths alternating with 0.3–0.4 kilobase repetitive sequence elements. The repetitive sequences of soybean DNA are interspersed both among themselves and among single copy regions of the genome.