Temperature Dependence of the DNA Double Helix at the Nanoscale: Structure,Elasticity, and Fluctuations

Biological organisms exist over a broad temperature range of −15°C to +120°C, where many molecular processes involving DNA depend on the nanoscale properties of the double helix. Here, we present results of extensive molecular dynamics simulations of DNA oligomers at different temperatures. We show that internal basepair conformations are strongly temperature-dependent, particularly in the stretch and opening degrees of freedom whose harmonic fluctuations can be considered the initial steps of the DNA melting pathway. The basepair step elasticity contains a weaker, but detectable, entropic contribution in the roll, tilt, and rise degrees of freedom. To extend the validity of our results to the temperature interval beyond the standard melting transition relevant to extremophiles, we estimate the effects of superhelical stress on the stability of the basepair steps, as computed from the Benham model. We predict that although the average twist decreases with temperature in vitro, the stabilizing external torque in vivo results in an increase of ∼1°/bp (or a superhelical density of Δσ ?  + 0.03) in the interval 0–100°C. In the final step, we show that the experimentally observed apparent bending persistence length of torsionally unconstrained DNA can be calculated from a hybrid model that accounts for the softening of the double helix and the presence of transient denaturation bubbles. Although the latter dominate the behavior close to the melting transition, the inclusion of helix softening is important around standard physiological temperatures.     

Elongation of repetitive DNA by DNA polymerase from a hyperthermophilic bacterium Thermus thermophilus

Short repetitive DNA sequences are believed to be one of the primordial genetic elements that served as a source of complex large DNA found in the genome of modern organisms. However, the mechanism of its expansion (increase in repeat number) during the course of evolution is unclear. We demonstrate that the DNA polymerase of the hyperthermophilic bacterium Thermus thermophilus can elongate oligoDNA with several tandem repeats to very long DNA in vitro. For instance, 48mer repetitive oligoDNA (TACATGTA)₆, which has 25% GC content and a palindromic sequence, can be elongated up to ~10 000 bases by DNA polymerase at 74°C without template DNA. OligoDNA having a different GC content or a quasi-palindromic sequence can also be elongated, but less efficiently. A spectroscopic thermal melting experiment with the oligoDNA showed that its hairpin–coil transition temperature was very close to the elongation reaction temperature (74°C), but was much higher than the temperature at which duplex oligoDNA can exist stably. Taken together, we conclude that repetitive oligoDNA with a palindromic or quasi-palindromic sequence is elongated extensively by a hyperthermophilic DNA polymerase through hairpin–coil transitions. We propose that such an elongation mechanism might have been a driving force to expand primordial short DNA.     

Polysomes from winter rye seedlings grown at low temperature : I. Size class distribution, composition, and stability

We have studied the influence of growth at low temperature on size class distribution, stability and composition of leaf cytoplasmic polysomes from rye seedlings (Secale cereale, cv Puma) grown at 5°C and at 20°C. Leaves of seedlings grown at 5°C contain 2.7 times more cytoplasmic polysomes (expressed on a DNA basis) and the polysome size class distribution is skewed toward larger polysomes. These changes were more pronounced in the free polysome fraction than in the membrane-bound fraction. The melting point of the total ribosome fraction from cold-grown leaves was decreased by 3.7°C. Electrophoresis did not reveal any difference in the rRNA or in core-ribosomal proteins (KCl nondissociable) following growth at low temperature. Some differences were noted in peripheral ribosomal proteins. This study is the first to examine the effect of growth at low and high temperatures on polysome metabolism using plants of similar developmental stage. Polysome quantity, polymerization, melting point and peripheral ribosomal proteins in rye seedlings are modified during growth at low temperature.     

Multivariate curve resolution: a powerful tool for the analysis of conformational transitions in nucleic acids

A successful application is reported of the multivariate curve resolution alternating least-squares method (MCR-ALS) for the analysis of nucleic acid melting and salt-induced transitions. Under conditions where several structures co-exist in a conformational equilibrium, MCR-ALS analysis of the UV and circular dichroism (CD) spectra at different temperatures, ionic strength and oligonucleotide concentration allows for the resolution of concentration profiles and pure spectra of the different species. The methodology is illustrated by the case of the cyclic oligonucleotide d<pTGCTCGCT>. The melting transition of this molecule at different oligonucleotide concentrations was studied at 0, 2 and 10 mM MgCl₂ by UV and CD spectroscopy. In addition, salt titration experiments were carried out at 21.0 and 54.0°C. The MCR-ALS analysis indicates that three different conformations of this molecule co-exist in solution. In agreement with previous NMR studies, these conformations were assigned to a monomeric dumbbell-like structure, a dimeric four-stranded conformation and a disordered (random coil) structure. The MCR-ALS methodology allows for a detailed analysis of how this equilibrium is affected by temperature, salt and oligonucleotide concentration.     

The conformation of double-stranded DNA inside bacteriophages depends on capsid size and shape

The packaging of double-stranded DNA into bacteriophages leads to the arrangement of the genetic material into highly-packed and ordered structures. Although modern experimental techniques reveal the most probable location of DNA inside viral capsids, the individual conformations of DNA are yet to be determined. In the current study we present the results of molecular dynamics simulations of the DNA packaging into several bacteriophages performed within the framework of a coarse-grained model. The final DNA conformations depend on the size and shape of the capsid, as well as the size of the protein portal, if any. In particular, isometric capsids with small or absent portals tend to form concentric spools, whereas the presence of a large portal favors coaxial spooling; slightly and highly elongated capsids result in folded and twisted toroidal conformations, respectively. The results of the simulations also suggest that the predominant factor in defining the global DNA arrangement inside bacteriophages is the minimization of the bending stress upon packaging.     

Recent developments in single-molecule DNA mechanics

Over the past two decades, measurements on individual stretched and twisted DNA molecules have helped define the basic elastic properties of the double helix and enabled real-time functional assays of DNA-associated molecular machines. Recently, new magnetic tweezers approaches for simultaneously measuring freely fluctuating twist and extension have begun to shed light on the structural dynamics of large nucleoprotein complexes. Related technical advances have facilitated direct measurements of DNA torque, contributing to a better understanding of abrupt structural transitions in mechanically stressed DNA. The new measurements have also been exploited in studies that hint at a developing synergistic relationship between single-molecule manipulation and structural DNA nanotechnology.     

Exploratory studies on azole carboxamides as nucleobase analogs: thermal denaturation studies on oligodeoxyribonucleotide duplexes containing pyrrole-3-carboxamide.

In order to study base pairing properties of the amide group in DNA duplexes, a nucleoside analog, 1-(2'-deoxy-beta-D-ribofuranosyl)pyrrole-3-carboxamide, was synthesized by a new route from the ester, methyl 1-(2'-deoxy-3',5'-di-O-p -toluoyl-beta-D-erythro-pentofuranosyl)pyrrole-3-carboxylate, obtained from the coupling reaction between 1-chloro-2-deoxy-3,5-di-O -toluoyl-d-erythropentofuranose and methyl pyrrole-3-carboxylate by treatment with dimethylaluminum amide. 1-(2'-Deoxy-beta-D-ribofuranosyl)pyrrole-3-carboxamide was incorporated into a series of oligodeoxyribonucleotides by solid-phase phosphoramidite technology. The corresponding oligodeoxyribonucleotides with 3-nitropyrrole in the same position in the sequence were synthesized for UV comparison of helix-coil transitions. The thermal melting studies indicate that pyrrole-3-carboxamide, which could conceptually adopt either a dA-like or a dI-like hydrogen bond conformation, pairs with significantly higher affinity to T than to dC. Pyrrole-3-carboxamide further resembles dA in the relative order of its base pairing preferences (T >dG >dA >dC). Theoretical calculations on the model compound N-methylpyrrole-3-carboxamide using density functional theory show little difference in the preference for a syntau versus anti conformation about the bond from pyrrole C3 to the amide carbonyl. The amide groups in both the minimized antitau and syntau conformations are twisted out of the plane of the pyrrole ring by 6-14 degrees. This twist may be one source of destabilization when the amide group is placed in the helix. Another contribution to the difference in stability between the base pairs of pyrrole-3-carboxamide with T and pyrrole-3-carboxamide with C may be the presence of a hydrogen bond in the former involving an acidic proton (N3-H of T).     

Histone tails cooperate to control the breathing of genomic nucleosomes

Genomic DNA is packaged in chromatin, a dynamic fiber variable in size and compaction. In chromatin, repeating nucleosome units wrap 145–147 DNA basepairs around histone proteins. Genetic and epigenetic regulation of genes relies on structural transitions in chromatin which are driven by intra- and inter-nucleosome dynamics and modulated by chemical modifications of the unstructured terminal tails of histones. Here we demonstrate how the interplay between histone H3 and H2A tails control ample nucleosome breathing motions. We monitored large openings of two genomic nucleosomes, and only moderate breathing of an engineered nucleosome in atomistic molecular simulations amounting to 24 μs. Transitions between open and closed nucleosome conformations were mediated by the displacement and changes in compaction of the two histone tails. These motions involved changes in the DNA interaction profiles of clusters of epigenetic regulatory aminoacids in the tails. Removing the histone tails resulted in a large increase of the amplitude of nucleosome breathing but did not change the sequence dependent pattern of the motions. Histone tail modulated nucleosome breathing is a key mechanism of chromatin dynamics with important implications for epigenetic regulation.     

Integration Host Factor of Mycobacterium tuberculosis,mIHF, Compacts DNA by a Bending Mechanism

The bacterial chromosomal DNA is folded into a compact structure called as ‘nucleoid’ so that the bacterial genome can be accommodated inside the cell. The shape and size of the nucleoid are determined by several factors including DNA supercoiling, macromolecular crowding and nucleoid associated proteins (NAPs). NAPs bind to different sites of the genome in sequence specific or non-sequence specific manner and play an important role in DNA compaction as well as regulation. Until recently, few NAPs have been discovered in mycobacteria owing to poor sequence similarities with other histone-like proteins of eubacteria. Several putative NAPs have now been identified in Mycobacteria on the basis of enriched basic residues or histone-like “PAKK” motifs. Here, we investigate mycobacterial Integration Host Factor (mIHF) for its architectural roles as a NAP using atomic force microscopy and DNA compaction experiments. We demonstrate that mIHF binds DNA in a non-sequence specific manner and compacts it by a DNA bending mechanism. AFM experiments also indicate a dual architectural role for mIHF in DNA compaction as well as relaxation. These results suggest a convergent evolution in the mechanism of E. coli and mycobacterial IHF in DNA compaction.     

Structural and torsional properties of the RAD51-dsDNA nucleoprotein filament

Human RAD51 is a key protein in the repair of DNA by homologous recombination. Its assembly onto DNA, which induces changes in DNA structure, results in the formation of a nucleoprotein filament that forms the basis of strand exchange. Here, we determine the structural and mechanical properties of RAD51-dsDNA filaments. Our measurements use two recently developed magnetic tweezers assays, freely orbiting magnetic tweezers and magnetic torque tweezers, designed to measure the twist and torque of individual molecules. By directly monitoring changes in DNA twist on RAD51 binding, we determine the unwinding angle per RAD51 monomer to be 45°, in quantitative agreement with that of its bacterial homolog, RecA. Measurements of the torque that is built up when RAD51-dsDNA filaments are twisted show that under conditions that suppress ATP hydrolysis the torsional persistence length of the RAD51-dsDNA filament exceeds that of its RecA counterpart by a factor of three. Examination of the filament’s torsional stiffness for different combinations of divalent ions and nucleotide cofactors reveals that the Ca²⁺ ion, apart from suppressing ATPase activity, plays a key role in increasing the torsional stiffness of the filament. These quantitative measurements of RAD51-imposed DNA distortions and accumulated mechanical stress suggest a finely tuned interplay between chemical and mechanical interactions within the RAD51 nucleoprotein filament.     

Gramicidin Channels Are Internally Gated

Gramicidin channels are archetypal molecular subjects for solid-state NMR studies and investigations of single-channel or cation conductance. Until now, the transitions between on and off conductance states have been thought, based on multichannel studies, to represent monomer ↔ dimer reactions. Here we use a single-molecule deposition method (vesicle fusion to a planar bilayer) to show that gramicidin dimer channels do not normally dissociate when conductance terminates. Furthermore, the observation of two ¹³C peaks in solid-state NMR indicates very stable dichotomous conformations for both the first and second peptide bonds in the monomers, and a two-dimensional chemical exchange spectrum with a 12-s mixing time demonstrates that the Val₁ carbonyl conformations exchange slowly, with lifetimes of several seconds. It is proposed that gramicidin channels are gated by small conformational changes in the channel near the permeation pathway. These studies demonstrate how regulation of conformations governing closed ↔ open transitions may be achieved and studied at the molecular level.     

Evidence of different G-quadruplex DNA binding with biogenic polyamines probed by electrospray ionization-quadrupole time of flight mass spectrometry,circular dichroism and atomic force microscopy

The binding properties of five G-quadruplex oligonucleotides (humtel24, k-ras32, c-myc22, c-kit1 and c-kit2) with polyamines have been investigated by electrospray ionization-quadrupole time of flight mass spectrometry, circular dichroism, melting temperature, atomic force microscopy (AFM) and molecular simulation. The MS results demonstrated that the polyamines and G-quadruplex DNA can form complexes with high affinity, and one molecule of G-quadruplex DNA can combine several molecules (1–5) of polyamines. The binding affinities of the polyamines to DNA were in the order of spermine > spermidine > putrescine. After binding with polyamines, the conformations of the G-quadruplex DNA were significantly changed, and spermine can induce the configurations of k-ras32 and c-kit1 to deviate from their G-quadruplex structures at high concentrations. In the presence of K⁺, the conformations of G-quadruplex DNA were stabilized, while polyamines can also induced alterations of their configurations. Melting temperature experiments suggested that the T_m of the DNA–polyamine complexes obviously increased both in the absence and presence of K⁺. The AFM results indicated that polyamines can induce aggregation of G-quadruplex DNA. Above results illustrated that the polyamines bound with the phosphate backbone and the base-pairs of G-quadruplex structures. Combining with the molecular simulation, the binding mode of the G-quadruplex DNA and polyamines were discussed. The results obtained would be beneficial for understanding the biological and physiological functions of polyamines and provide useful information for development of antitumor drugs.     

Spectral Analysis of Naturally Occurring Methylxanthines (Theophylline,Theobromine and Caffeine) Binding with DNA

Nucleic acids exist in a dynamic equilibrium with a number of molecules that constantly interact with them and regulate the cellular activities. The inherent nature of the structure and conformational integrity of these macromolecules can lead to altered biological activity through proper targeting of nucleic acids binding ligands or drug molecules. We studied the interaction of naturally occurring methylxanthines such as theophylline, theobromine and caffeine with DNA, using UV absorption and Fourier transform infrared (FTIR) spectroscopic methods, and especially monitored their binding affinity in the presence of Mg²⁺ and during helix-coil transitions of DNA by temperature (T_m) or pH melting profiles. The study indicates that all these molecules effectively bind to DNA in a dose dependent manner. The overall binding constants of DNA-theophylline = 3.5×10³ M⁻¹, DNA-theobromine = 1.1×10³ M⁻¹, and DNA-Caffeine = 3.8×10³ M⁻¹. On the other hand T_m/pH melting profiles showed 24–35% of enhanced binding activity of methylxanthines during helix-coil transitions of DNA rather than to its native double helical structure. The FTIR analysis divulged that theophylline, theobromine and caffeine interact with all the base pairs of DNA (A-T; G-C) and phosphate group through hydrogen bond (H-bond) interaction. In the presence of Mg²⁺, methylxanthines altered the structure of DNA from B to A-family. However, the B-family structure of DNA remained unaltered in DNA-methylxanthines complexes or in the absence of Mg²⁺. The spectral analyses indicated the order of binding affinity as “caffeine≥theophylline>theobromine” to the native double helical DNA, and “theophylline≥theobromine>caffeine to the denatured form of DNA and in the presence of divalent metal ions.     

Stretched DNA Investigated Using Molecular-Dynamics and Quantum-Mechanical Calculations

We combined atomistic molecular-dynamics simulations with quantum-mechanical calculations to investigate the sequence dependence of the stretching behavior of duplex DNA. Our combined quantum-mechanical/molecular-mechanical approach demonstrates that molecular-mechanical force fields are able to describe both the backbone and base-base interactions within the highly distorted nucleic acid structures produced by stretching the DNA from the 5′ ends, which include conformations containing disassociated basepairs, just as well as these force fields describe relaxed DNA conformations. The molecular-dynamics simulations indicate that the force-induced melting pathway is sequence-dependent and is influenced by the availability of noncanonical hydrogen-bond interactions that can assist the disassociation of the DNA basepairs. The biological implications of these results are discussed.     

A general method for manipulating DNA sequences from any organism with optical tweezers

Mechanical manipulation of single DNA molecules can provide novel information about DNA properties and protein–DNA interactions. Here we describe and characterize a useful method for manipulating desired DNA sequences from any organism with optical tweezers. Molecules are produced from either genomic or cloned DNA by PCR using labeled primers and are tethered between two optically trapped microspheres. We demonstrate that human, insect, plant, bacterial and viral sequences ranging from ~10 to 40 kilobasepairs can be manipulated. Force-extension measurements show that these constructs exhibit uniform elastic properties in accord with the expected contour lengths for the targeted sequences. Detailed protocols for preparing and manipulating these molecules are presented, and tethering efficiency is characterized as a function of DNA concentration, ionic strength and pH. Attachment strength is characterized by measuring the unbinding time as a function of applied force. An alternative stronger attachment method using an amino–carboxyl linkage, which allows for reliable DNA overstretching, is also described.     

Temperature Dependence of the DNA Double Helix at the Nanoscale: Structure,Elasticity, and Fluctuations

Biological organisms exist over a broad temperature range of −15°C to +120°C, where many molecular processes involving DNA depend on the nanoscale properties of the double helix. Here, we present results of extensive molecular dynamics simulations of DNA oligomers at different temperatures. We show that internal basepair conformations are strongly temperature-dependent, particularly in the stretch and opening degrees of freedom whose harmonic fluctuations can be considered the initial steps of the DNA melting pathway. The basepair step elasticity contains a weaker, but detectable, entropic contribution in the roll, tilt, and rise degrees of freedom. To extend the validity of our results to the temperature interval beyond the standard melting transition relevant to extremophiles, we estimate the effects of superhelical stress on the stability of the basepair steps, as computed from the Benham model. We predict that although the average twist decreases with temperature in vitro, the stabilizing external torque in vivo results in an increase of ∼1°/bp (or a superhelical density of Δσ?+0.03$\text{Δ}\sigma ?+0.03$) in the interval 0–100°C. In the final step, we show that the experimentally observed apparent bending persistence length of torsionally unconstrained DNA can be calculated from a hybrid model that accounts for the softening of the double helix and the presence of transient denaturation bubbles. Although the latter dominate the behavior close to the melting transition, the inclusion of helix softening is important around standard physiological temperatures.     

DNA melting investigated by differential scanning calorimetry and Raman spectroscopy.

Thermal denaturation of the B form of double-stranded DNA has been probed by differential scanning calorimetry (DSC) and Raman spectroscopy of 160 base pair (bp) fragments of calf thymus DNA. The DSC results indicate a median melting temperature Tm = 75.5 degrees C with calorimetric enthalpy change delta Hcal = 6.7 kcal/mol (bp), van't Hoff enthalpy change delta HVH = 50.4 kcal/mol (cooperative unit), and calorimetric entropy change delta Scal = 19.3 cal/deg.mol (bp), at the experimental conditions of 55 mg DNA/ml in 5 mM sodium cacodylate at pH 6.4. The average cooperative melting unit (nmelt) comprises 7.5 bp. The Raman signature of 160 bp DNA is highly sensitive to temperature. Analyses of several conformation-sensitive Raman bands indicate the following ranges for thermodynamic parameters of melting: 43 < delta HVH < 61 kcal/mol (cooperative unit), 75 < Tm < 80 degrees C and 6 < (nmelt) < 9 bp, consistent with the DSC results. The changes observed in specific Raman band frequencies and intensities as a function of temperature reveal that thermal denaturation is accompanied by disruption of Watson-Crick base pairs, unstacking of the bases and disordering of the B form backbone. These three types of structural change are highly correlated throughout the investigated temperature range of 20 to 93 degrees C. Raman bands diagnostic of purine and pyrimidine unstacking, conformational rearrangements in the deoxyribose-phosphate moieties, and changes in environment of phosphate groups have been identified. Among these, bands at 834 cm-1 (due to a localized vibration of the phosphodiester group), 1240 cm-1 (thymine ring) and 1668 cm-1 (carbonyl groups of dT, dG and dC), are shown by comparison with DSC results to be the most reliable quantitative indicators of DNA melting. Conversely, the intensities of Raman marker bands at 786 cm-1 (cytosine ring), 1014 cm-1 (deoxyribose ring) and 1092 cm-1 (phosphate group) are largely invariant to melting and are proposed as appropriate standards for intensity normalizations.     

Loop L5 Assumes Three Distinct Orientations during the ATPase Cycle of the Mitotic Kinesin Eg5: A TRANSIENT AND TIME-RESOLVED FLUORESCENCE STUDY*

Members of the kinesin superfamily of molecular motors differ in several key structural domains, which probably allows these molecular motors to serve the different physiologies required of them. One of the most variable of these is a stem-loop motif referred to as L5. This loop is longest in the mitotic kinesin Eg5, and previous structural studies have shown that it can assume different conformations in different nucleotide states. However, enzymatic domains often consist of a mixture of conformations whose distribution shifts in response to substrate binding or product release, and this information is not available from the “static” images that structural studies provide. We have addressed this issue in the case of Eg5 by attaching a fluorescent probe to L5 and examining its fluorescence, using both steady state and time-resolved methods. This reveals that L5 assumes an equilibrium mixture of three orientations that differ in their local environment and segmental mobility. Combining these studies with transient state kinetics demonstrates that there is a major shift in this distribution during transitions that interconvert weak and strong microtubule binding states. Finally, in conjunction with previous cryo-EM reconstructions of Eg5·microtubule complexes, these fluorescence studies suggest a model in which L5 regulates both nucleotide and microtubule binding through a set of reversible interactions with helix α3. We propose that these features facilitate the production of sustained opposing force by Eg5, which underlies its role in supporting formation of a bipolar spindle in mitosis.