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.     

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.     

Cooperative kinking at distant sites in mechanically stressed DNA

In cells, DNA is routinely subjected to significant levels of bending and twisting. In some cases, such as under physiological levels of supercoiling, DNA can be so highly strained, that it transitions into non-canonical structural conformations that are capable of relieving mechanical stress within the template. DNA minicircles offer a robust model system to study stress-induced DNA structures. Using DNA minicircles on the order of 100 bp in size, we have been able to control the bending and torsional stresses within a looped DNA construct. Through a combination of cryo-EM image reconstructions, Bal31 sensitivity assays and Brownian dynamics simulations, we have been able to analyze the effects of biologically relevant underwinding-induced kinks in DNA on the overall shape of DNA minicircles. Our results indicate that strongly underwound DNA minicircles, which mimic the physical behavior of small regulatory DNA loops, minimize their free energy by undergoing sequential, cooperative kinking at two sites that are located about 180° apart along the periphery of the minicircle. This novel form of structural cooperativity in DNA demonstrates that bending strain can localize hyperflexible kinks within the DNA template, which in turn reduces the energetic cost to tightly loop 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.     

A unique endonuclease from Crithidia fasciculata which recognizes a bend in the DNA helix. Specificity of the cleavage reaction

The introduction of a single nick in DNA circles by Crithidia fasciculata nicking enzyme (Shlomai, J., and Linial, M. (1986) J. Biol. Chem. 261, 16219-16225) requires the presence of a bent structure in the DNA helix. However, the sequence directing the local bending of the DNA helix is not per se a preferred site for nicking by the enzyme. No extensive sequence specificity is involved in defining the cleavage site for C. fasciculata nicking enzyme in the duplex circular DNA substrate. However, the abundance of A and T residues is significantly high at both the 3' and the 5' termini generated at the nicked site. Nicking of the sequence-directed bent fragment from C. fasciculata kinetoplast DNA minicircles correlates with the periodicity determined by the unique nucleotide distribution in the bent sequence, reflected in its thermodynamic parameters. Occurrence of nicking is best correlated with the predicted minima of the melting temperature and delta G profiles, as well as with A and T dinucleotide sequences at the nicked site, in both the supercoiled and the relaxed sequence-directed bent DNA substrates. The potential role of the bend-dependent nicking reaction in the replication of kinetoplast DNA minicircles is discussed.     

Visualization of the bent helix in kinetoplast DNA by electron microscopy

Kinetoplast DNA minicircles from the trypanosomatid Crithidia fasciculata contain a segment of approximately 200 bp which is probably more highly bent than any other DNA previously studied. Electron microscopy (EM) of relaxed minicircles (2.5 kb) revealed 200-300 bp loops within the larger circles, and the loops could also be detected on full-length linear molecules. Examination by EM of a 219 bp cloned fragment which contains the bent helix revealed that up to 70% of the molecules appeared circular whether or not the ends were cohesive. In contrast, a 207 bp fragment from pBR322 showed no circles and the fragments in general appeared much straighter than the kinetoplast fragments. Treatment of the 219 bp bent kinetoplast fragment with the drug distamycin caused a striking reduction in curvature.     

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.     

A computational study of nucleosomal DNA flexibility

Molecular dynamics simulations of the nucleosome core particle and its isolated DNA free in solution are reported. The simulations are based on the implicit solvent methodology and provide insights into the nature of large-scale structural fluctuations and flexibility of the nucleosomal DNA. In addition to the kinked regions previously identified in the x-ray structure of the nucleosome, the simulations support the existence of a biochemically identified distorted region of the DNA. Comparison of computed relative free energies shows that formation of the kinks is associated with little, if any, energy cost relative to a smooth, ideal conformation of the DNA superhelix. Isolated nucleosomal DNA is found to be considerably more flexible than expected for a 147 bp stretch of DNA based on its canonical persistence length of 500 A. Notably, the significant bending of the DNA observed in our simulations occurs without breaking of Watson-Crick bonds. The computed relative stability of bent conformations is sensitive to the ionic strength of the solution in the physiological range; the sensitivity suggests possible experiments that might provide further insights into the structural origins of the unusual flexibility of the DNA.     

The chromosomal protein MC1 from the archaebacterium Methanosarcina sp. CHTI 55 induces DNA bending and supercoiling.

We have investigated the effect on the DNA structure of protein MC1, a basic and small polypeptide (Mr 10700) representing the major chromosomal protein in Methanosarcinaceae. The ability of protein MC1 to strongly favour cyclization upon polymerization of short DNA fragments by T4 DNA ligase indicates that protein MC1 mediates DNA bending. Several negatively supercoiled topoisomers of minicircles were obtained with DNA fragments of 203 and 146 bp, their distribution depends upon the amount of protein MC1 complexed with DNA. In addition, protein MC1 can induce a compaction of a nicked plasmid.     

3D reconstruction and comparison of shapes of DNA minicircles observed by cryo-electron microscopy

We use cryo-electron microscopy to compare 3D shapes of 158 bp long DNA minicircles that differ only in the sequence within an 18 bp block containing either a TATA box or a catabolite activator protein binding site. We present a sorting algorithm that correlates the reconstructed shapes and groups them into distinct categories. We conclude that the presence of the TATA box sequence, which is believed to be easily bent, does not significantly affect the observed shapes.     

DNA supercoiling changes the spacing requirement of two lac operators for DNA loop formation with lac repressor.

We have used a gel retardation assay to investigate the influence of DNA supercoiling on loop formation between lac repressor and two lac operators. A series of 15 DNA minicircles of identical size (452 bp) was constructed carrying two lac operators at distances ranging from 153 to 168 bp. Low positive or negative supercoiling (sigma = +/- 0.023) changed the spacing between the two lac operators required for the formation of the most stable loops. This reveals the presence of altered double helical repeats (ranging from 10.3 to 10.7 bp) in supercoiled DNA minicircles. At elevated negative supercoiling (sigma = -0.046) extremely stable loops were formed at all operator distances tested, with a slight spacing periodicity remaining. After relaxation of minicircle-repressor complexes with topoisomerase I one superhelical turn was found to be constrained in those minicircles which carry operators at distances corresponding to a non-integral number of helical turns. This indicates that DNA loop formation can define local DNA domains with altered topological properties of the DNA helix.     

The sequence-directed bent structure in kinetoplast DNA is recognized by an enzyme from Crithidia fasciculata

Crithidia fasciculata nicking enzyme (Shlomai, J., and Linial, M. (1986) J. Biol. Chem. 261, 16219-16225) interrupts a single phosphodiester bond in duplex DNA circles from various sources, only in their supercoiled form, but not following their relaxation by DNA topoisomerases. However, this requirement for DNA substrate supercoiling was not observed using the natural kinetoplast DNA as a substrate. Relaxed kinetoplast DNA minicircles, either free or topologically linked, were efficiently nicked by the enzyme. Furthermore, bacterial plasmids, containing a unit length kinetoplast DNA minicircle insert, were used as substrates for nicking in their relaxed form. This capacity to activate a relaxed DNA topoisomer as a substrate for nicking is an intrinsic property of the sequence-directed bend, naturally present in kinetoplast DNA. The 211-base pair fragment of the bent region from C. fasciculata kinetoplast DNA could support the nicking of a relaxed DNA substrate in a reaction dependent upon the DNA helix curvature.     

Cationic metals promote sequence-directed DNA bending

A DNA segment of approximately 200 base pairs (bp) from Crithidia fasciculata kinetoplast minicircles was previously shown by electron microscopy (EM) to bend into a small circle due to its unique nucleotide sequence containing repeated blocks of 4-6 A's. When this segment was flanked by 207 bp of plasmid DNA on one side and 460 bp on the other, the resulting 890-bp DNA was found to appear either relatively straight or extremely bent as visualized by EM. The bend was located one-third the distance from one end. The fraction of molecules with the most extreme bend increased from approximately 2% to 50-60% following incubation of the DNA with increasing concentrations of Zn2+, Co2+, Ba2+, and Mn2+. These observations suggest that sequence-directed bending in DNA is an inducible and not a static phenomenon. Possible roles of transitions between the bent and straight conformations in the control of gene expression are discussed.     

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.     

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.     

Pankinetoplast DNA structure in a primitive bodonid flagellate, Cryptobia helicis.

The mitochondrial DNA (mtDNA) of a primitive kinetoplastid flagellate Cryptobia helicis is composed of 4.2 kb minicircles and 43 kb maxicircles. 85% and 6% of the minicircles are in the form of supercoiled (SC) and relaxed (OC) monomers, respectively. The remaining minicircles (9%) constitute catenated oligomers composed of both the SC and OC molecules. Minicircles contain bent helix and sequences homologous to the minicircle conserved sequence blocks. Maxicircles encode typical mitochondrial genes and are not catenated. The mtDNA, which we describe with the term 'pankinetoplast DNA', is spread throughout the mitochondrial lumen, where it is associated with multiple electron-lucent loci. There are approximately 8400 minicircles per pankinetoplast-mitochondrion, with the pan-kDNA representing approximately 36% of the total cellular DNA. Based on the similarity of the C.helicis minicircles to plasmids, we present a theory on the formation of the kDNA network.     

'Empty' minicircles and petB/atpA and psbD/psbE (cytb559 alpha) genes in tandem in Amphidinium carterae plastid DNA

Amphidinium carterae minicircle chloroplast DNA was separated from total DNA by centrifugation through a sucrose/NaCl gradient. Sequences of minicircles with psbA and 23S rRNA contained a common region of 67 bp. Primers designed from this generated numerous polymerase chain reaction products of 1.5-2.6 kb. These contained psaA and psaB as one gene/circle, and petB/atpA and psbD/psbE as two genes/circle. 'Empty' minicircles of 1.7-2.5 kb containing no identifiable genes or parts of genes were more abundant than gene-containing circles. From 15 minicircles a minimum common region of 48 bp was identified, with little identity to that from other dinoflagellate minicircles.     

DNA gyrase can supercoil DNA circles as small as 174 base pairs.

DNA gyrase introduces negative supercoils into closed-circular DNA using the free energy of ATP hydrolysis. Consideration of steric and thermodynamic aspects of the supercoiling reaction indicates that there should be a lower limit to the size of DNA circle which can be supercoiled by gyrase. We have investigated the supercoiling reaction of circles from 116-427 base pairs (bp) in size and have determined that gyrase can supercoil certain relaxed isomers of circles as small as 174 bp, dependent on the final superhelix density of the supercoiled product. Furthermore, this limiting superhelical density (-0.11) is the same as that determined for the supercoiling of plasmid pBR322. We also find that although circles in the range 116-152 bp cannot be supercoiled, they can nevertheless be relaxed by gyrase when positively supercoiled. These data suggest that the conformational changes associated with the supercoiling reaction can be carried out by gyrase in a circle as small as 116 bp. We discuss these results with respect to the thermodynamics of DNA supercoiling and steric aspects of the gyrase mechanism.     

Characterizing a Histidine Switch Controlling pH-Dependent Conformational Changes of the Influenza Virus Hemagglutinin

During the fusion of the influenza virus to the host cell, bending of the HA2 chain of hemagglutinin into a hairpin-shaped structure in a pH-dependent manner facilitates the fusion of the viral envelope and the endosomal membrane. To characterize the structural and dynamical responses of the hinge region of HA2 to pH changes and examine the role of a conserved histidine in this region (the hinge histidine), we have performed an extensive set of molecular dynamics (MD) simulations of 26-residue peptides encompassing the hinge regions of several hemagglutinin subtypes under both neutral and low pH conditions, modeled by the change of the protonation state of the hinge histidine. More than 70 sets of MD simulations (collectively amounting to 25.1 μs) were performed in both implicit and explicit solvents to study the effect of histidine protonation on structural dynamics of the hinge region. In both explicit and implicit solvent simulations, hinge bending was consistently observed upon the protonation of the histidine in all the simulations starting with an initial straight helical conformation, whereas the systems with a neutral histidine retained their primarily straight conformation throughout the simulations. Conversely, the MD simulations starting from an initially bent conformation resulted in the formation of a straight helical structure upon the neutralization of the hinge histidine, whereas the bent structure was maintained when the hinge histidine remained protonated. Finally, mutation of the hinge histidine to alanine abolishes the bending response of the peptide altogether. A molecular mechanism based on the interaction of the hinge histidine with neighboring acidic residues is proposed to be responsible for its role in controlling the conformation of the hinge. We propose that this might present a common mechanism for pH-controlled structural changes in helical structures when histidines act as the pH sensor.