Bending modes of DNA directly addressed by cryo-electron microscopy of DNA minicircles

We use cryo-electron microscopy (cryo-EM) to study the 3D shapes of 94-bp-long DNA minicircles and address the question of whether cyclization of such short DNA molecules necessitates the formation of sharp, localized kinks in DNA or whether the necessary bending can be redistributed and accomplished within the limits of the elastic, standard model of DNA flexibility. By comparing the shapes of covalently closed, nicked and gapped DNA minicircles, we conclude that 94-bp-long covalently closed and nicked DNA minicircles do not show sharp kinks while gapped DNA molecules, containing very flexible single-stranded regions, do show sharp kinks. We corroborate the results of cryo-EM studies by using Bal31 nuclease to probe for the existence of kinks in 94-bp-long minicircles.     

Kinking occurs during molecular dynamics simulations of small DNA minicircles

Recent experiments on minicircle formation suggest that a conformational mechanism other than smooth deformation may be playing a role in enhancing DNA flexibility. Both local base unpairing and kink formation have been suggested as possible explanations. Although kinks within isolated DNA were proposed 30 years ago, they have, until now, only been observed within DNA complexed with proteins. In order to test how DNA behaves in the strong bending regime, we have carried out molecular dynamics simulations of a 94 base pair minicircle in explicit solvent with two different linking numbers, corresponding to a torsionally relaxed state and a positively supercoiled state. The simulations suggest that sharp kinks can indeed arise in small minicircles. The relaxed minicircle is generally associated with a single kink, while two kinks occur with the supercoiled state. No evidence is seen of base unpaired regions.     

Atomistic simulations reveal bubbles, kinks and wrinkles in supercoiled DNA

Although DNA is frequently bent and supercoiled in the cell, much of the available information on DNA structure at the atomistic level is restricted to short linear sequences. We report atomistic molecular dynamics (MD) simulations of a series of DNA minicircles containing between 65 and 110 bp which we compare with a recent biochemical study of structural distortions in these tight DNA loops. We have observed a wealth of non-canonical DNA structures such as kinks, denaturation bubbles and wrinkled conformations that form in response to bending and torsional stress. The simulations show that bending alone is sufficient to induce the formation of kinks in circles containing only 65 bp, but we did not observe any defects in simulations of larger torsionally relaxed circles containing 110 bp over the same MD timescales. We also observed that under-winding in minicircles ranging in size from 65 to 110 bp leads to the formation of single stranded bubbles and wrinkles. These calculations are used to assess the ability of atomistic MD simulations to determine the structure of bent and supercoiled DNA.     

Kinking the double helix by bending deformation

DNA bending and torsional deformations that often occur during its functioning inside the cell can cause local disruptions of the regular helical structure. The disruptions created by negative torsional stress have been studied in detail, but those caused by bending stress have only been analyzed theoretically. By probing the structure of very small DNA circles, we determined that bending stress disrupts the regular helical structure when the radius of DNA curvature is smaller than 3.5 nm. First, we developed an efficient method to obtain covalently closed DNA minicircles. To detect structural disruptions in the minicircles we treated them by single-strand-specific endonucleases. The data showed that the regular DNA structure is disrupted by bending deformation in the 64–65-bp minicircles, but not in the 85–86-bp minicircles. Our results suggest that strong DNA bending initiates kink formation while preserving base pairing.     

Theoretical Analysis of Disruptions in DNA Minicircles

Under sufficient bending stress, which appears in DNA minicircles and small DNA loops, the double helix experiences local disruptions of its regular structure. We developed a statistical-mechanical treatment of the disruptions in DNA minicircles, studied experimentally by Du et al. The model of disruptions used in our Monte Carlo simulation of minicircle conformations specifies these conformations by three parameters: DNA bend angle at the disruption, θ_d; local DNA unwinding caused by the disruption; and the free energy associated with the disruption in the unstressed double helix, G_d. The model is applicable to any structural type of disruption, kinks or opening of single basepairs. The simulation shows that accounting for both torsional and bending deformation associated with the disruptions is very important for proper analysis. We obtained a relationship between values of G_d and θ_d under which the simulation results are compatible with the experimental data. The relationship suggests that the free energy of basepair opening, which includes flipping out both bases, is significantly higher than the generally accepted value. The model is also applied to the analysis of j-factors of very short DNA fragments.     

Replication of kinetoplast DNA in isolated kinetoplasts from Crithidia fasciculata. Identification of minicircle DNA replication intermediates

The kinetoplast DNA (kDNA) of trypanosomes is comprised of thousands of DNA minicircles and 20-50 maxicircles catenated into a single network. We show that kinetoplasts isolated from the trypanosomatid species Crithidia fasciculata incorporate labeled nucleotides and support minicircle DNA replication in a manner which mimics two characteristics of minicircle replication in vivo: 1) the minicircles are replicated as free molecules and subsequently reattached to the kDNA network, and 2) a replication intermediate having a structure consistent with a highly gapped minicircle species is generated. In addition, a class of minicircle DNA replication intermediates is observed containing discontinuities at specific sites within each of the newly synthesized DNA strands. By using a strain of C. fasciculata possessing nearly homogenous minicircles, we were able to map the discontinuities to two small regions situated 180 degrees apart on the minicircle. Each region has two sites at which a discontinuity can occur, one on each strand and separated by approximately 100 base pairs. These sites may represent origins of minicircle DNA replication.     

Orientation of DNA Minicircles Balances Density and Topological Complexity in Kinetoplast DNA

Kinetoplast DNA (kDNA), a unique mitochondrial structure common to trypanosomatid parasites, contains thousands of DNA minicircles that are densely packed and can be topologically linked into a chain mail-like network. Experimental data indicate that every minicircle in the network is, on average, singly linked to three other minicircles (i.e., has mean valence 3) before replication and to six minicircles in the late stages of replication. The biophysical factors that determine the topology of the network and its changes during the cell cycle remain unknown. Using a mathematical modeling approach, we previously showed that volume confinement alone can drive the formation of the network and that it induces a linear relationship between mean valence and minicircle density. Our modeling also predicted a minicircle valence two orders of magnitude greater than that observed in kDNA. To determine the factors that contribute to this discrepancy we systematically analyzed the relationship between the topological properties of the network (i.e., minicircle density and mean valence) and its biophysical properties such as DNA bending, electrostatic repulsion, and minicircle relative position and orientation. Significantly, our results showed that most of the discrepancy between the theoretical and experimental observations can be accounted for by the orientation of the minicircles with volume exclusion due to electrostatic interactions and DNA bending playing smaller roles. Our results are in agreement with the three dimensional kDNA organization model, initially proposed by Delain and Riou, in which minicircles are oriented almost perpendicular to the horizontal plane of the kDNA disk. We suggest that while minicircle confinement drives the formation of kDNA networks, it is minicircle orientation that regulates the topological complexity of the network.     

Structural Mechanics of DNA Wrapping in the Nucleosome

Experimental X-ray crystal structures and a database of calculated structural parameters of DNA octamers were used in combination to analyse the mechanics of DNA bending in the nucleosome core complex. The 1kx5 X-ray crystal structure of the nucleosome core complex was used to determine the relationship between local structure at the base-step level and the global superhelical conformation observed for nucleosome-bound DNA. The superhelix is characterised by a large curvature (597°) in one plane and very little curvature (10°) in the orthogonal plane. Analysis of the curvature at the level of 10-step segments shows that there is a uniform curvature of 30° per helical turn throughout most of the structure but that there are two sharper kinks of 50° at ± 2 helical turns from the central dyad base pair. The curvature is due almost entirely to the base-step parameter roll. There are large periodic variations in roll, which are in phase with the helical twist and account for 500° of the total curvature. Although variations in the other base-step parameters perturb the local path of the DNA, they make minimal contributions to the total curvature. This implies that DNA bending in the nucleosome is achieved using the roll-slide-twist degree of freedom previously identified as the major degree of freedom in naked DNA oligomers. The energetics of bending into a nucleosome-bound conformation were therefore analysed using a database of structural parameters that we have previously developed for naked DNA oligomers. The minimum energy roll, the roll flexibility force constant and the maximum and minimum accessible roll values were obtained for each base step in the relevant octanucleotide context to account for the effects of conformational coupling that vary with sequence context. The distribution of base-step roll values and corresponding strain energy required to bend DNA into the nucleosome-bound conformation defined by the 1kx5 structure were obtained by applying a constant bending moment. When a single bending moment was applied to the entire sequence, the local details of the calculated structure did not match the experiment. However, when local 10-step bending moments were applied separately, the calculated structure showed excellent agreement with experiment. This implies that the protein applies variable bending forces along the DNA to maintain the superhelical path required for nucleosome wrapping. In particular, the 50° kinks are constraints imposed by the protein rather than a feature of the 1kx5 DNA sequence. The kinks coincide with a relatively flexible region of the sequence, and this is probably a prerequisite for high-affinity nucleosome binding, but the bending strain energy is significantly higher at these points than for the rest of the sequence. In the most rigid regions of the sequence, a higher strain energy is also required to achieve the standard 30° curvature per helical turn. We conclude that matching of the DNA sequence to the local roll periodicity required to achieve bending, together with the increased flexibility required at the kinks, determines the sequence selectivity of DNA wrapping in the nucleosome.     

Strong bending of the DNA double helix

During the past decade, the issue of strong bending of the double helix has attracted a lot of attention. Here, we overview the major experimental and theoretical developments in the field sorting out reliably established facts from speculations and unsubstantiated claims. Theoretical analysis shows that sharp bends or kinks have to facilitate strong bending of the double helix. It remains to be determined what is the critical curvature of DNA that prompts the appearance of the kinks. Different experimental and computational approaches to the problem are analyzed. We conclude that there is no reliable evidence that any anomalous behavior of the double helix happens when DNA fragments in the range of 100 bp are circularized without torsional stress. The anomaly starts at the fragment length of about 70 bp when sharp bends or kinks emerge in essentially every molecule. Experimental data and theoretical analysis suggest that kinks may represent openings of isolated base pairs, which had been experimentally detected in linear DNA molecules. The calculation suggests that although the probability of these openings in unstressed DNA is close to 10⁻⁵, it increases sharply in small DNA circles reaching 1 open bp per circle of 70 bp.     

Role of p38 in replication of Trypanosoma brucei kinetoplast DNA

Trypanosomes have an unusual mitochondrial genome, called kinetoplast DNA, that is a giant network containing thousands of interlocked minicircles. During kinetoplast DNA synthesis, minicircles are released from the network for replication as theta-structures, and then the free minicircle progeny reattach to the network. We report that a mitochondrial protein, which we term p38, functions in kinetoplast DNA replication. RNA interference (RNAi) of p38 resulted in loss of kinetoplast DNA and accumulation of a novel free minicircle species named fraction S. Fraction S minicircles are so underwound that on isolation they become highly negatively supertwisted and develop a region of Z-DNA. p38 binds to minicircle sequences within the replication origin. We conclude that cells with RNAi-induced loss of p38 cannot initiate minicircle replication, although they can extensively unwind free minicircles.     

The structure of replicating kinetoplast DNA networks

Kinetoplast DNA (kDNA), the mitochondrial DNA of Crithidia fasciculata and related trypanosomatids, is a network containing approximately 5,000 covalently closed minicircles which are topologically interlocked. kDNA synthesis involves release of covalently closed minicircles from the network, and, after replication of the free minicircles, reattachment of the nicked or gapped progeny minicircles to the network periphery. We have investigated this process by electron microscopy of networks at different stages of replication. The distribution of nicked and closed minicircles is easily detectable either by autoradiography of networks radiolabeled at endogenous nicks by nick translation or by twisting the covalently closed minicircles with intercalating dye. The location of newly synthesized minicircles within the network is determined by autoradiography of network is determined by autoradiography of networks labeled in vivo with a pulse of [3H]thymidine. These studies have clarified structural changes in the network during replication, the timing of repair of nicked minicircles after replication, and the mechanism of division of the network.     

Structural organization of kinetoplast DNA and its compactization in a model system in vitro

Studies on compactization and decompactization of the genome are of great importance for elucidation of structural mechanisms taking part in the regulation of gene activity. Kinetoplast DNA (kpDNA) is a convenient model for studies of compactization processes. KpDNA represents unique structure ("network"), consisting of catenated circular molecules of two types: minicircles (900 b.p.) and maxicircles (40 000 b.p.). The compactization process of kpDNA in vitro caused by interaction with synthetic peptide-dansylhydraside trivaline was studied. It was shown that at the initial stages the hairpins are observed on minicircles as if triple rings are being organized. The formation of hairpin is probably favoured by the presence in the minicircles of bent DNA, a specific nucleotide sequence causing rigid bending of the DNA helix. The hairpin does not make contact with the neighbouring DNA segment to form a triple ring, because the sizes of minicircles are too small. The minicircles compactization is finished with a complete collapse of the minicircles with the formation of rod-like structures. The catenation causes branching of rod-like structures. As a result of their intermolecular interaction, the branched rod-like structures become thicker. The process is completed with formation of the compact network, its diameter being 3-6 times smaller compared to the initial one.     

Short-range template switching in great ape genomes explored using pair hidden Markov models

Many complex genomic rearrangements arise through template switch errors, which occur in DNA replication when there is a transient polymerase switch to an alternate template nearby in three-dimensional space. While typically investigated at kilobase-to-megabase scales, the genomic and evolutionary consequences of this mutational process are not well characterised at smaller scales, where they are often interpreted as clusters of independent substitutions, insertions and deletions. Here we present an improved statistical approach using pair hidden Markov models, and use it to detect and describe short-range template switches underlying clusters of mutations in the multi-way alignment of hominid genomes. Using robust statistics derived from evolutionary genomic simulations, we show that template switch events have been widespread in the evolution of the great apes’ genomes and provide a parsimonious explanation for the presence of many complex mutation clusters in their phylogenetic context. Larger-scale mechanisms of genome rearrangement are typically associated with structural features around breakpoints, and accordingly we show that atypical patterns of secondary structure formation and DNA bending are present at the initial template switch loci. Our methods improve on previous non-probabilistic approaches for computational detection of template switch mutations, allowing the statistical significance of events to be assessed. By specifying realistic evolutionary parameters based on the genomes and taxa involved, our methods can be readily adapted to other intra- or inter-species comparisons.     

Intramitochondrial location and dynamics of Crithidia fasciculata kinetoplast minicircle replication intermediates

Kinetoplast DNA, the mitochondrial DNA of Crithidia fasciculata, is organized into a network containing 5,000 topologically interlocked minicircles. This network, situated within the mitochondrial matrix, is condensed into a disk-shaped structure located near the basal body of the flagellum. Fluorescence in situ hybridization revealed that before their replication, minicircles are released vectorially from the network face nearest the flagellum. Replication initiates in the zone between the flagellar face of the disk and the mitochondrial membrane (we term this region the kinetoflagellar zone [KFZ]). The replicating minicircles then move to two antipodal sites that flank the disk-shaped network. In later stages of replication, the number of free minicircles increases, accumulating transiently in the KFZ. The final replication events, including primer removal, repair of many of the gaps, and reattachment of the progeny minicircles to the network periphery, are thought to take place within the antipodal sites.     

The majority of minicircle DNA in Crithidia fasciculata strain CF-C1 is of a single class with nearly homogeneous DNA sequence

DNA minicircles found within the kinetoplast of the trypanosomatid Crithidia fasciculata, like those of most other kinetoplastid species, are heterogeneous in sequence. The pattern of minicircle DNA fragments generated by cleavage of kinetoplast DNA with various restriction enzymes has been used to demonstrate this heterogeneity. Here we describe a strain of Crithidia fasciculata in which more than 90% of the DNA minicircles exhibit a common pattern of restriction enzyme cleavage sites. A map of cleavage sites within this major minicircle DNA class is presented for seven restriction enzymes with hexanucleotide recognition sequences. Sequence homogeneity at an even finer level is reflected in minicircle DNA digestion patterns generated by restriction enzymes with tetranucleotide recognition sites. Partial DNA sequence analysis of multiple clones from the major minicircle class shows nearly complete homogeneity at the nucleotide level. The existence of a near homogeneous complement of DNA minicircles in Crithidia should facilitate the study of their replication in this organism.     

Characterization of the molecular components in kinetoplast- mitochondrial DNA of Trypanosoma equiperdum. Comparative study of the dyskinetoplastic and wild strains

The structure of the kinetoplast DNA of Trypanosoma equiperdum has been studied and compared to the structure of the circular mitochondrial DNA extracted from a dyskinetoplastic strain of T. equiperdum. In T. equiperdum wild type, the kinetoplast DNA constitutes approximately 6% of the total cellular DNA and is composed of approximately 3,000 supercoiled minicircles of 6.4 x 10(5) daltons and approximately 50 circular supercoiled molecules of 15.4 x 10(6) daltons topologically interlocked; The buoyant density in CsCl of the minicircles is 1.691 g/cm 3. The large circles have a buoyant density of 1.684 g/cm 3, are homogeneous in size and are selectively cleaved by several restriction endonucleases which do not cleave the minicircles. The cleavage sites of six different restriction endonucleases have been mapped on the large circle. The minicircles are cleaved by two other restriction endonucleases, and their cleavage sites have been mapped. The mitochondrial DNA extracted from the dyskinetoplastic strain of T. equiperdum represents 7% of the total DNA of the cell and is composed of supercoiled circles, heterogeneous in size, and topologically associated in catenated oligomers. Its buoyant density in CsCl is 1.688 g/cm 3. These molecules are not cleaved by any of the eight restriction endonucleases tested. The reassociation kinetics of in vitro labeled kDNA minicircles and large circles has been studied. The results indicate that the minicircles as well as the large circles are homogeneous in sequence and that the circular DNA of the dyskinetoplastic strain has no sequence in common with the kDNA of the wild strain.     

Selection for arsenite resistance causes reversible changes in minicircle composition and kinetoplast organization in Leishmania mexicana.

Certain minor minicircle sequence classes in the kinetoplast DNA (kDNA) networks of arsenite- or tunicamycin-resistant Leishmania mexicana amazonensis variants whose nuclear DNA is amplified appear to be preferentially selected to replicate (S. T. Lee, C. Tarn, and K. P. Chang, Mol. Biochem. Parasitol. 58:187-204, 1993). These sequences replace the predominant wild-type minicircle sequences to become dominant species in the kDNA network. The switch from wild-type-specific to variant-specific minicircles takes place rapidly within the same network, the period of minicircle dominance changes being defined as the transition period. To investigate the structural organization of the kDNA networks during this transition period, we analyzed kDNA from whole arsenite-resistant Leishmania parasites by dot hybridization with sequence-specific DNA probes and by electron-microscopic examination of isolated kDNA networks in vitro. Both analyses concluded that during the switch of dominance the predominant wild-type minicircle class was rapidly lost and that selective replication of variant-specific minicircles subsequently filled the network step by step. There was a time during the transition when few wild-type- or variant-specific minicircles were present, leaving the network almost empty and exposing a species of thick, long, fibrous DNA which seemed to form a skeleton for the network. Both minicircles and maxicircles were found to attach to these long DNA fibrils. The nature of the long DNA fibrils is not clear, but they may be important in providing a framework for the network structure and a support for the replication of minicircles and maxicircles.