Probing the Elasticity of DNA on Short Length Scales by Modeling Supercoiling under Tension

The wormlike-chain (WLC) model is widely used to describe the energetics of DNA bending. Motivated by recent experiments, alternative, so-called subelastic chain models were proposed that predict a lower elastic energy of highly bent DNA conformations. Until now, no unambiguous verification of these models has been obtained because probing the elasticity of DNA on short length scales remains challenging. Here we investigate the limits of the WLC model using coarse-grained Monte Carlo simulations to model the supercoiling of linear DNA molecules under tension. At a critical supercoiling density, the DNA extension decreases abruptly due to the sudden formation of a plectonemic structure. This buckling transition is caused by the large energy required to form the tightly bent end-loop of the plectoneme and should therefore provide a sensitive benchmark for model evaluation. Although simulations based on the WLC energetics could quantitatively reproduce the buckling measured in magnetic tweezers experiments, the buckling almost disappears for the tested linear subelastic chain model. Thus, our data support the validity of a harmonic bending potential even for small bending radii down to 3.5 nm.     

Conformational Effects of DNA Stretching

DNA is an extensible molecule, and an extended conformation of DNA is involved in some biological processes. We have examined the effect of elongation stress on the conformational properties of DNA base pairs by conformational analysis. The calculations show that stretching does significantly affect the conformational properties and flexibilities of base pairs. In particular, we have found that the propeller twist in base pairs reverses its sign upon stretching. The energy profile analysis indicates that electrostatic interactions make a major contribution to the stabilization of the positive-propeller-twist configuration in stretched DNA. This stretching also results in a monotonic decrease in the helical twist angle, tending to unwind the double helix. Fluctuations in most variables initially increase upon stretching, because of unstacking of base pairs, but then the fluctuations decrease as DNA is stretched further, owing to the formation of specific interactions between base pairs induced by the positive propeller twist. Thus, the stretching of DNA has particularly significant effects upon DNA flexibility. These changes in both the conformation and flexibility of base pairs probably have a role in functional interactions with proteins.     

Conformational effects of DNA stretching

DNA is an extensible molecule, and an extended conformation of DNA is involved in some biological processes. We have examined the effect of elongation stress on the conformational properties of DNA base pairs by conformational analysis. The calculations show that stretching does significantly affect the conformational properties and flexibilities of base pairs. In particular, we have found that the propeller twist in base pairs reverses its sign upon stretching. The energy profile analysis indicates that electrostatic interactions make a major contribution to the stabilization of the positive-propeller-twist configuration in stretched DNA. This stretching also results in a monotonic decrease in the helical twist angle, tending to unwind the double helix. Fluctuations in most variables initially increase upon stretching, because of unstacking of base pairs, but then the fluctuations decrease as DNA is stretched further, owing to the formation of specific interactions between base pairs induced by the positive propeller twist. Thus, the stretching of DNA has particularly significant effects upon DNA flexibility. These changes in both the conformation and flexibility of base pairs probably have a role in functional interactions with proteins.     

Probing the elastic limit of DNA bending

Sharp bending of double-stranded DNA (dsDNA) plays an essential role in genome structure and function. However, the elastic limit of dsDNA bending remains controversial. Here, we measured the opening rates of small dsDNA loops with contour lengths ranging between 40 and 200 bp using single-molecule Fluorescence Resonance Energy Transfer. The relationship of loop lifetime to loop size revealed a critical transition in bending stress. Above the critical loop size, the loop lifetime changed with loop size in a manner consistent with elastic bending stress, but below it, became less sensitive to loop size, indicative of softened dsDNA. The critical loop size increased from ∼60 bp to ∼100 bp with the addition of 5 mM magnesium. We show that our result is in quantitative agreement with the kinkable worm-like chain model, and furthermore, can reproduce previously reported looping probabilities of dsDNA over the range between 50 and 200 bp. Our findings shed new light on the energetics of sharply bent dsDNA.     

Probing the Flexibility of Large Conformational Changes in Protein Structures through Local Perturbations

Protein conformational changes and dynamic behavior are fundamental for such processes as catalysis, regulation, and substrate recognition. Although protein dynamics have been successfully explored in computer simulation, there is an intermediate-scale of motions that has proven difficult to simulate—the motion of individual segments or domains that move independently of the body the protein. Here, we introduce a molecular-dynamics perturbation method, the Rotamerically Induced Perturbation (RIP), which can generate large, coherent motions of structural elements in picoseconds by applying large torsional perturbations to individual sidechains. Despite the large-scale motions, secondary structure elements remain intact without the need for applying backbone positional restraints. Owing to its computational efficiency, RIP can be applied to every residue in a protein, producing a global map of deformability. This map is remarkably sparse, with the dominant sites of deformation generally found on the protein surface. The global map can be used to identify loops and helices that are less tightly bound to the protein and thus are likely sites of dynamic modulation that may have important functional consequences. Additionally, they identify individual residues that have the potential to drive large-scale coherent conformational change. Applying RIP to two well-studied proteins, Dihdydrofolate Reductase and Triosephosphate Isomerase, which possess functionally-relevant mobile loops that fluctuate on the microsecond/millisecond timescale, the RIP deformation map identifies and recapitulates the flexibility of these elements. In contrast, the RIP deformation map of α-lytic protease, a kinetically stable protein, results in a map with no significant deformations. In the N-terminal domain of HSP90, the RIP deformation map clearly identifies the ligand-binding lid as a highly flexible region capable of large conformational changes. In the Estrogen Receptor ligand-binding domain, the RIP deformation map is quite sparse except for one large conformational change involving Helix-12, which is the structural element that allosterically links ligand binding to receptor activation. RIP analysis has the potential to discover sites of functional conformational changes and the linchpin residues critical in determining these conformational states.     

Conformational Properties of the TATA-Box Binding Sequence of DNA

Abstract

Nanosecond scale molecular dynamics simulations in water demonstrate that the DNA oligomer, GCGTATATAAAACGC, which contains a target site for the TATA-box binding protein (TBP), has an intrinsic preference to adopt an A-like conformation in the region of the TATA-box and undergoes bending related to that seen within in the TBP complex. This result is obtained from two independent simulations using different starting structures. In line with earlier suggestions of Guzikevich-Guerstein and Shakked, these simulations imply that an A-DNA conformation may be an important intermediate step in forming the strongly distorted DNA observed within the crystallographically determined complex with TBP. These results also support modeling studies by Lebrun et al. which suggest that the TBP binding mechanism can be broken down into a backbone transition to an A-like form coupled with a mechanical distortion which locally stretches and unwinds the DNA.     

Conformational constraints in nuclear DNA.

We have investigated DNA superstructure in a wide range of nuclei of higher cells by gently lysing cells to release structures that resemble nuclei but are depleted of nuclear proteins. The sedimentation properties of these structures, which we call nucleoids, have been examined in sucrose gradients containing the intercalating agent, ethidium. The sedimentation rate of nucleoids derived from the growing cells of mammals, birds, amphibians and insects varies in the manner characteristic of circular and superhelical molecules of DNA. These characteristic changes in sedimentation rate are abolished by irradiating the nucleoids with low doses of gamma-rays, a procedure known to introduce single-strand scissions into DNA. We have also investigated by similar means DNA superstructure in nucleoids derived from a variety of different chick cells. Nucleoids derived from adult hen erythrocytes differ from the other nucleoids studied in that their sedimentation rate does not vary in the manner characteristic of supercoiled DNA.     

Length Mutations in Human Mitochondrial DNA

By high-resolution, restriction mapping of mitochondrial DNAs purified from 112 human individuals, we have identified 14 length variants caused by small additions and deletions (from about 6 to 14 base pairs in length). Three of the 14 length differences are due to mutations at two locations within the D loop, whereas the remaining 11 occur at seven sites that are probably within other noncoding sequences and at junctions between coding sequences. In five of the nine regions of length polymorphism, there is a sequence of five cytosines in a row, this sequence being comparatively rare in coding DNA. Phylogenetic analysis indicates that, in most of the polymorphic regions, a given length mutation has arisen several times independently in different human lineages. The average rate at which length mutations have been arising and surviving in the human species is estimated to be many times higher for noncoding mtDNA than for noncoding nuclear DNA. The mystery of why vertebrate mtDNA is more prone than nuclear DNA to evolve by point mutation is now compounded by the discovery of a similar bias toward rapid evolution by length mutation.     

Distributions of circular DNA according its topological states

A distinctive feature of closed circular DNA molecules is their particular topological state, which cannot be altered by any conformational rearrangement short of breaking at least one strand. This topological constraint opens unique possibilities for experimental studies of the distributions of topological states created in different ways. Primarily, the equilibrium distributions of topological properties are considered in the review. It is described how such distributions can be obtained and measured experimentally, and how they can be computed. Comparison of the calculated and measured equilibrium distributions over the linking number of complementary strands, equilibrium fractions of knots and links formed by circular molecules has provided much valuable information about the properties of the double helix. Study of the steady-state fraction of knots and links created by type II DNA topoisomerases has revealed a surprising property of the enzymes: their ability to reduce these fractions considerably below the equilibrium level.     

Conformational changes in DNA and soft modes

In this paper we explore the applicability of the soft mode approach to study the conformational transitions of DNA. It is believed that the A-B conformation change is a first order transition. Soft mode theories only apply to the initial stages of a first order transition. However the mode softening in such a case can be the initiating factor which ultimately leads to the transition. The first order transition is, then, a breakdown of what otherwise would have been a true second order transition. The mode softening is causally connected to onset of the transition. We use the eigenvectors obtained from lattice dynamics calculations to identify the softmode. We use the eigenvector projections to form a force constant matrix that is required to drive a mode soft. We explore the methods by which this force constant matrix can be formed. We suggest that the breaking of specific "water bridges" between phosphate groups in the two single strands can drive the conformation change.     

Probing the electrostatic shielding of DNA with capillary electrophoresis

The free solution mobility of a 20-bp double-stranded DNA oligomer has been measured in diethylmalonate (DM) and Tris-acetate buffers, with and without added NaCl or TrisCl. DM buffers have the advantage that the buffering ion is anionic, so the cation composition in the solution can be varied at will. The results indicate that the free solution mobility of DNA decreases linearly with the logarithm of ionic strength when the ionic strength is increased by increasing the buffer concentration. The mobility also decreases linearly with the logarithm of ionic strength when NaCl is added to NaDM buffer or TrisCl is added to TrisDM buffer. Nonlinear effects are observed if the counterion in the added salt differs from the counterion in the buffer. The dependence of the mobility on ionic strength cannot be predicted using the Henry, Debye-Hückel-Onsager, or Pitts equations for electrophoresis. However, the mobilities observed in all buffer and buffer/salt solutions can be predicted within approximately 20% by the Manning equation for electrophoresis, using no adjustable parameters. The results suggest that the electrostatic shielding of DNA is determined not only by the relative concentrations of the various ions in the solution, but also by their equivalent conductivities.     

Conformational Changes in DNA and Soft Modes

Abstract

In this paper we explore the applicability of the soft mode approach to study the conformational transitions of DNA. It is believed that the A-B conformation change is a first order transition. Soft mode theories only apply to the initial stages of a first order transition. However the mode softening in such a case can be the initiating factor which ultimately leads to the transition. The first order transition is, then, a breakdown of what otherwise would have been a true second order transition. The mode softening is causally connected to onset of the transition. We use the eigenvectors obtained from lattice dynamics calculations to identify the softmode. We use the eigenvector projections to form a force constimi matrix that is required to drive a mode soft. We explore the methods by which this force constant matrix can be formed. We suggest that the breaking of specific “water bridges” between phosphate groups in the two single strands can drive the conformation change.     

Sequence Recognition of DNA by Protein-Induced Conformational Transitions

The binding of proteins to specific sequences of DNA is an important feature of virtually all DNA transactions. Proteins recognize specific DNA sequences using both direct readout (sensing types and positions of DNA functional groups) and indirect readout (sensing DNA conformation and deformability). Previously we showed that the P22 c2 repressor N-terminal domain (P22R NTD) forces the central non-contacted 5'-ATAT-3' sequence of the DNA operator into the B′ state, a state known to affect DNA hydration, rigidity and bending. Usually the B′ state, with a narrow minor groove and a spine of hydration, is reserved for A-tract DNA (TpA steps disrupt A-tracts). Here, we have co-crystallized P22R NTD with an operator containing a central 5′-ACGT-3′ sequence in the non-contacted region. C·G base pairs have not previously been observed in the B′ state and are thought to prevent it. However, P22R NTD induces a narrow minor groove and a spine of hydration to 5'-ACGT-3'. We observe that C·G base pairs have distinctive destabilizing and disordering effects on the spine of hydration. It appears that the reduced stability of the spine results in a higher energy cost for the B to B′ transition. The differential effect of DNA sequence on the barrier to this transition allows the protein to sense the non-contacted DNA sequence.     

Conformational studies of antisense DNA by PFG NMR

Pulsed field gradient diffusion constant measurements were used to resolve the ambiguity in determining the conformational states of single-stranded DNA dodecanucleotides (d1s, d4s and d5s). For d1s and d5s, because of the spectral symmetry conventional NMR analyses cannot differentiate whether they are hairpins or homo-duplexes. However, the diffusion constants of these sequences at 300 K are 1.4 times greater than those of the comparison complementary duplexes. This result agrees well with what is expected for D_hairpin>/D_duplex based on classic liquid-phase translational diffusion models and the Einstein–Stokes equation, confirming that d1s and d5s form hairpins. d4s did not show a structured spectral pattern, but its diffusion constant measurement suggests that this sequence may not be a random coil. The DNA sequences studied contain chemically modified backbone linkages and are potential antisense agents for gene regulation. The knowledge of their diffusion constants, in combination with conventional NMR analysis and other biophysical spectroscopic measurements, provides new insights into the relationships of chemical structure and conformational preference of antisense oligonucleotides and their analogs.     

DNA Stable-Isotope Probing (DNA-SIP)

DNA stable-isotope probing (DNA-SIP) is a powerful technique for identifying active microorganisms that assimilate particular carbon substrates and nutrients into cellular biomass. As such, this cultivation-independent technique has been an important methodology for assigning metabolic function to the diverse communities inhabiting a wide range of terrestrial and aquatic environments. Following the incubation of an environmental sample with stable-isotope labelled compounds, extracted nucleic acid is subjected to density gradient ultracentrifugation and subsequent gradient fractionation to separate nucleic acids of differing densities. Purification of DNA from cesium chloride retrieves labelled and unlabelled DNA for subsequent molecular characterization (e.g. fingerprinting, microarrays, clone libraries, metagenomics). This JoVE video protocol provides visual step-by-step explanations of the protocol for density gradient ultracentrifugation, gradient fractionation and recovery of labelled DNA. The protocol also includes sample SIP data and highlights important tips and cautions that must be considered to ensure a successful DNA-SIP analysis.Download video file.^{(159M, mp4)}