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.     

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.     

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.     

Modeling DNA Thermodynamics under Torsional Stress

Negatively twisted DNA is essential to many biological functions. Due to torsional stress, duplex DNA can have local, sequence-dependent structural defects. In this work, a thermodynamic model of DNA was built to qualitatively predict the local sequence-dependent mechanical instabilities under torsional stress. The results were compared to both simulation of a coarse-grained model and experiment results. By using the Kirkwood superposition approximation, we built an analytical model to represent the free energy difference ΔW of a hydrogen-bonded basepair between the B-form helical state and the basepair opened (or locally melted) state, within a given sequence under torsional stress. We showed that ΔW can be well approximated by two-body interactions with its nearest-sequence-neighbor basepairs plus a free energy correction due to long-range correlations. This model is capable of rapidly predicting the position and thermodynamics of local defects in a given sequence. The result qualitatively matches with an in vitro experiment for a long DNA sequence (>4000 basepairs). The 12 parameters used in this model can be further quantitatively refined when more experimental data are available.     

NMR structure of the DNA decamer duplex containing double T*G mismatches of cis-syn cyclobutane pyrimidine dimer: implications for DNA damage recognition by the XPC-hHR23B complex

The cis–syn cyclobutane pyrimidine dimer (CPD) is a cytotoxic, mutagenic and carcinogenic DNA photoproduct and is repaired by the nucleotide excision repair (NER) pathway in mammalian cells. The XPC–hHR23B complex as the initiator of global genomic NER binds to sites of certain kinds of DNA damage. Although CPDs are rarely recognized by the XPC–hHR23B complex, the presence of mismatched bases opposite a CPD significantly increased the binding affinity of the XPC–hHR23B complex to the CPD. In order to decipher the properties of the DNA structures that determine the binding affinity for XPC–hHR23B to DNA, we carried out structural analyses of the various types of CPDs by NMR spectroscopy. The DNA duplex which contains a single 3′ T·G wobble pair in a CPD (CPD/GA duplex) induces little conformational distortion. However, severe distortion of the helical conformation occurs when a CPD contains double T·G wobble pairs (CPD/GG duplex) even though the T residues of the CPD form stable hydrogen bonds with the opposite G residues. The helical bending angle of the CPD/GG duplex was larger than those of the CPD/GA duplex and properly matched CPD/AA duplex. The fluctuation of the backbone conformation and significant changes in the widths of the major and minor grooves at the double T·G wobble paired site were also observed in the CPD/GG duplex. These structural features were also found in a duplex that contains the (6–4) adduct, which is efficiently recognized by the XPC–hHR23B complex. Thus, we suggest that the unique structural features of the DNA double helix (that is, helical bending, flexible backbone conformation, and significant changes of the major and/or minor grooves) might be important factors in determining the binding affinity of the XPC–hHR23B complex to DNA.     

DNA sequence and methylation prescribe the inside-out conformational dynamics and bending energetics of DNA minicircles

Eukaryotic genome and methylome encode DNA fragments’ propensity to form nucleosome particles. Although the mechanical properties of DNA possibly orchestrate such encoding, the definite link between ‘omics’ and DNA energetics has remained elusive. Here, we bridge the divide by examining the sequence-dependent energetics of highly bent DNA. Molecular dynamics simulations of 42 intact DNA minicircles reveal that each DNA minicircle undergoes inside-out conformational transitions with the most likely configuration uniquely prescribed by the nucleotide sequence and methylation of DNA. The minicircles’ local geometry consists of straight segments connected by sharp bends compressing the DNA’s inward-facing major groove. Such an uneven distribution of the bending stress favors minimum free energy configurations that avoid stiff base pair sequences at inward-facing major grooves. Analysis of the minicircles’ inside-out free energy landscapes yields a discrete worm-like chain model of bent DNA energetics that accurately account for its nucleotide sequence and methylation. Experimentally measuring the dependence of the DNA looping time on the DNA sequence validates the model. When applied to a nucleosome-like DNA configuration, the model quantitatively reproduces yeast and human genomes’ nucleosome occupancy. Further analyses of the genome-wide chromatin structure data suggest that DNA bending energetics is a fundamental determinant of genome architecture.     

Direct measurements of the nucleosome-forming preferences of periodic DNA motifs challenge established models

Several periodic motifs have been implicated in facilitating the bending of DNA around the histone core of the nucleosome. For example, di-nucleotides AA/TT/TA and GC at ∼10-bp periods, but offset by 5 bp, are found with higher-than-expected occurrences in aligned nucleosomal DNAs in vitro and in vivo. Additionally, regularly oscillating period-10 trinucleotide motifs non-T, A/T, G and their complements have been implicated in the formation of regular nucleosome arrays. The effects of these periodic motifs on nucleosome formation have not been systematically tested directly by competitive reconstitution assays. We show that, in general, none of these period-10 motifs, except TA, in certain sequence contexts, facilitates nucleosome formation. The influence of periodic TAs on nucleosome formation is appreciable; with some of the 200-bp DNAs out-competing bulk nucleosomal DNA by more than 400-fold. Only the nucleotides immediately flanking TA influence its nucleosome-forming ability. Period-10 TA, when flanked by a pair of permissive nucleotides, facilitates DNA bending through compression of the minor groove. The free energy change for nucleosome formation decreases linearly with the number of consecutive TAs, up to eight. We suggest how these data can be reconciled with previous findings.     

An Evaluation of the Crystal Structure of C-terminal Truncated Apolipoprotein A-I in Solution Reveals Structural Dynamics Related to Lipid Binding

Apolipoprotein (apo) A-I mediates many of the anti-atherogenic functions attributed to high density lipoprotein. Unfortunately, efforts toward a high resolution structure of full-length apoA-I have not been fruitful, although there have been successes with deletion mutants. Recently, a C-terminal truncation (apoA-I^Δ185–243) was crystallized as a dimer. The structure showed two helical bundles connected by a long, curved pair of swapped helical domains. To compare this structure to that existing under solution conditions, we applied small angle x-ray scattering and isotope-assisted chemical cross-linking to apoA-I^Δ185–243 in its dimeric and monomeric forms. For the dimer, we found evidence for the shared domains and aspects of the N-terminal bundles, but not the molecular curvature seen in the crystal. We also found that the N-terminal bundles equilibrate between open and closed states. Interestingly, this movement is one of the transitions proposed during lipid binding. The monomer was consistent with a model in which the long shared helix doubles back onto the helical bundle. Combined with the crystal structure, these data offer an important starting point to understand the molecular details of high density lipoprotein biogenesis.     

Structure of discontinuities in kinetoplast DNA-associated minicircles during S phase in Crithidia fasciculata

Kinetoplast DNA (kDNA) is a novel form of mitochondrial DNA consisting of thousands of interlocked minicircles and 20–30 maxicircles. The minicircles replicate free of the kDNA network but nicks and gaps in the newly synthesized strands remain at the time of reattachment to the kDNA network. We show here that the steady-state population of replicated, network-associated minicircles only becomes repaired to the point of having nicks with a 3′OH and 5′deoxyribonucleoside monophosphate during S phase. These nicks represent the origin/terminus of the strand and occur within the replication origins (oriA and oriB) located 180° apart on the minicircle. Minicircles containing a new L strand have a single nick within either oriA or oriB but not in both origins in the same molecule. The discontinuously synthesized H strand contains single nicks within both oriA and oriB in the same molecule implying that discontinuities between the H-strand Okazaki fragments become repaired except for the fragments initiated within the two origins. Nicks in L and H strands at the origins persist throughout S phase and only become ligated as a prelude to network division. The failure to ligate these nicks until just prior to network division is not due to inappropriate termini for ligation.     

Mammalian chromatin substructure studies with the calcium–magnesium endonuclease and two-dimensional polyacrylamide-gel electrophoresis

The basic regularity of chromatin substructure that has been reported in rat liver chromatin (Hewish & Burgoyne, 1973b) was also detected in mouse chromatin. The regular series of DNA fragments produced by the action of Ca–Mg endonuclease on rat chromatin were studied further. The smallest single-stranded class has a molecular weight of approx. 45000–63000 and the smallest double-stranded class has a molecular weight of approx. 120000–150000. Studies of the substructure of the DNA fragments produced by the Ca–Mg endonuclease have shown that the regular series of double-stranded fragments have regular series of single-stranded fragments within them. It was concluded that the regular series of double-stranded fragments was probably a consequence of the regular series of single-stranded fragments. Digestion time-courses are presented for mouse and rat nuclear DNA.     

DNA phosphate crowding correlates with protein cationic side chain density and helical curvature in protein/DNA crystal structures

Sequence-specific binding of proteins to their DNA targets involves a complex spectrum of processes that often induce DNA conformational variation in the bound complex. The forces imposed by protein binding that cause the helical deformations are intimately interrelated and difficult to parse or rank in importance. To investigate the role of electrostatics in helical deformation, we quantified the relationship between protein cationic residue density (Cpc) and DNA phosphate crowding (Cpp). The correlation between Cpc and Cpp was then calculated for a subset of 58 high resolution protein–DNA crystal structures. Those structures containing strong Cpc/Cpp correlation (>±0.25) were likely to contain DNA helical curvature. Further, the correlation factor sign predicted the direction of helical curvature with positive (16 structures) and negative (seven structures) correlation containing concave (DNA curved toward protein) and convex (DNA curved away from protein) curvature, respectively. Protein–DNA complexes without significant Cpc/Cpp (36 structures) correlation (-0.25<0<0.25) tended to contain DNA without significant curvature. Interestingly, concave and convex complexes also include more arginine and lysine phosphate contacts, respectively, whereas linear complexes included essentially equivalent numbers of Lys/Arg phosphate contacts. Together, these findings suggest an important role for electrostatic interactions in protein–DNA complexes involving helical curvature.     

Helical coherence of DNA in crystals and solution

The twist, rise, slide, shift, tilt and roll between adjoining base pairs in DNA depend on the identity of the bases. The resulting dependence of the double helix conformation on the nucleotide sequence is important for DNA recognition by proteins, packaging and maintenance of genetic material, and other interactions involving DNA. This dependence, however, is obscured by poorly understood variations in the stacking geometry of the same adjoining base pairs within different sequence contexts. In this article, we approach the problem of sequence-dependent DNA conformation by statistical analysis of X-ray and NMR structures of DNA oligomers. We evaluate the corresponding helical coherence length—a cumulative parameter quantifying sequence-dependent deviations from the ideal double helix geometry. We find, e.g. that the solution structure of synthetic oligomers is characterized by 100–200 Å coherence length, which is similar to ~150 Å coherence length of natural, salmon-sperm DNA. Packing of oligomers in crystals dramatically alters their helical coherence. The coherence length increases to 800–1200 Å, consistent with its theoretically predicted role in interactions between DNA at close separations.     

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.     

Insights into structure, dynamics and hydration of locked nucleic acid (LNA) strand-based duplexes from molecular dynamics simulations

Locked nucleic acid (LNA) is a chemically modified nucleic acid with its sugar ring locked in an RNA-like (C3′-endo) conformation. LNAs show extraordinary thermal stabilities when hybridized with DNA, RNA or LNA itself. We performed molecular dynamics simulations on five isosequential duplexes (LNA–DNA, LNA–LNA, LNA–RNA, RNA–DNA and RNA–RNA) in order to characterize their structure, dynamics and hydration. Structurally, the LNA–DNA and LNA–RNA duplexes are found to be similar to regular RNA–DNA and RNA–RNA duplexes, whereas the LNA–LNA duplex is found to have its helix partly unwound and does not resemble RNA–RNA duplex in a number of properties. Duplexes with an LNA strand have on average longer interstrand phosphate distances compared to RNA–DNA and RNA–RNA duplexes. Furthermore, intrastrand phosphate distances in LNA strands are found to be shorter than in DNA and slightly shorter than in RNA. In case of induced sugar puckering, LNA is found to tune the sugar puckers in partner DNA strand toward C3′-endo conformations more efficiently than RNA. The LNA–LNA duplex has lesser backbone flexibility compared to the RNA–RNA duplex. Finally, LNA is less hydrated compared to DNA or RNA but is found to have a well-organized water structure.     

The structure of metallo-DNA with consecutive thymine–HgII–thymine base pairs explains positive entropy for the metallo base pair formation

We have determined the three-dimensional (3D) structure of DNA duplex that includes tandem Hg^II-mediated T–T base pairs (thymine–Hg^II–thymine, T–Hg^II–T) with NMR spectroscopy in solution. This is the first 3D structure of metallo-DNA (covalently metallated DNA) composed exclusively of ‘NATURAL’ bases. The T–Hg^II–T base pairs whose chemical structure was determined with the ¹⁵N NMR spectroscopy were well accommodated in a B-form double helix, mimicking normal Watson–Crick base pairs. The Hg atoms aligned along DNA helical axis were shielded from the bulk water. The complete dehydration of Hg atoms inside DNA explained the positive reaction entropy (ΔS) for the T–Hg^II–T base pair formation. The positive ΔS value arises owing to the Hg^II dehydration, which was approved with the 3D structure. The 3D structure explained extraordinary affinity of thymine towards Hg^II and revealed arrangement of T–Hg^II–T base pairs in metallo-DNA.     

kDNA minicircles of the major sequence class of C. fasciculata contain a single region of bent helix widely separated from the two origins of replication

The major sequence class of Crithidia fasciculata minicircles is shown to have a single region of bent helical DNA widely separated from the two replication origins located 180 degrees apart on the minicircle map. The position of the bend in the DNA has been mapped both by gel electrophoretic methods and by direct electron microscopic observation of the DNA. This sequence directed bending is apparently the result of homopolymeric dA X dT tracts 4-6 base pairs long repeated in phase with the helix screw. The region of the bend contains nineteen such homopolymeric tracts in a region of about 200 base pairs with sixteen of the tracts oriented in the same direction.     

Recognition and incision of site-specifically modified C8 guanine adducts formed by 2-aminofluorene, N-acetyl-2-aminofluorene and 1-nitropyrene by UvrABC nuclease

Nucleotide excision repair plays a crucial role in removing many types of DNA adducts formed by UV light and chemical carcinogens. We have examined the interactions of Escherichia coli UvrABC nuclease proteins with three site-specific C8 guanine adducts formed by the carcinogens 2-aminofluorene (AF), N-acetyl-2-acetylaminofluorene (AAF) and 1-nitropyrene (1-NP) in a 50mer oligonucleotide. Similar to the AF and AAF adducts, the 1-NP-induced DNA adduct contains an aminopyrene (AP) moiety covalently linked to the C8 position of guanine. The dissociation constants for UvrA binding to AF–, AAF– and AP–DNA adducts, determined by gel mobility shift assay, are 33 ± 9, 8 ± 2 and 23 ± 9 nM, respectively, indicating that the AAF adduct is recognized much more efficiently than the other two. Incision by UvrABC nuclease showed that AAF–DNA was cleaved ~2-fold more efficiently than AF– or AP–DNA (AAF > AF ≈ AP), even though AP has the largest molecular size in this group. However, an opened DNA structure of six bases around the adduct increased the incision efficiency for AF–DNA (but not for AP–DNA), making it equivalent to that for AAF–DNA. These results are consistent with a model in which DNA damage recognition by the E.coli nucleotide excision repair system consists of two sequential steps. It includes recognition of helical distortion in duplex DNA followed by recognition of the type of nucleotide chemical modification in a single-stranded region. The difference in incision efficiency between AF– and AAF–DNA adducts in normal DNA sequence, therefore, is a consequence of their difference in inducing structural distortions in DNA. The results of this study are discussed in the light of NMR solution structures of these DNA adducts.