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.     

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.     

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.     

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.     

Thermodynamics of transfer ribonucleic acids: the effect of sodium on the thermal unfolding of yeast tRNAPhe.

The thermal unfolding of yeast phenylalanine-specific tRNA (tRNA^Phe) has been calorimetrically investigated at several salt concentrations in the absence of magnesium. Application of the deconvolution theory of macromolecular conformational transitions allows calculation of the thermodynamic parameters of unfolding. It is demonstrated that the unfolding of tRNA^Phe occurs in a sequential fashion and that four separate transitions or five macromolecular thermodynamic states exist in the temperature range 8–72°C under the experimental conditions of these studies (0.067–0.52M Na⁺). The enthalpy and entropy changes between states and the relative population of each state as a function of temperature and salt concentration have been obtained. Sodium stabilizes the low-temperature conformations of tRNA^Phe. The increase in the melting temperatures of each transition is shown to be linearly dependent on the logarithm of sodium concentration. These results allow calculation of the “phase” diagram for the transitions as a function of salt concentration.     

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.     

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.     

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.     

Force-induced melting of the DNA double helix. 2. Effect of solution conditions

In this paper, we consider the implications of the general theory developed in the accompanying paper, to interpret experiments on DNA overstretching that involve variables such as solution temperature, pH, and ionic strength. We find the DNA helix-coil phase boundary in the force-temperature space. At temperatures significantly below the regular (zero force) DNA melting temperature, the overstretching force, f(ov)(T), is predicted to decrease nearly linearly with temperature. We calculate the slope of this dependence as a function of entropy and heat-capacity changes upon DNA melting. Fitting of the experimental f(ov)(T) dependence allows determination of both of these quantities in very good agreement with their calorimetric values. At temperatures slightly above the regular DNA melting temperature, we predict stabilization of dsDNA by moderate forces, and destabilization by higher forces. Thus the DNA stretching curves, f(b), should exhibit two rather than one overstretching transitions: from single stranded (ss) to double stranded (ds) and then back at the higher force. We also predict that any change in DNA solution conditions that affects its melting temperature should have a similar effect on DNA overstretching force. This result is used to calculate the dependence of DNA overstretching force on solution pH, f(ov)(pH), from the known dependence of DNA melting temperature on pH. The calculated f(ov)(pH) is in excellent agreement with its experimental determination (M. C. Williams, J. R. Wenner, I. Rouzina, and V. A. Bloomfield, Biophys. J., accepted for publication). Finally, we quantitatively explain the measured dependence of DNA overstretching force on solution ionic strength for crosslinked and noncrosslinked DNA. The much stronger salt dependence of f(ov) in noncrosslinked DNA results from its lower linear charge density in the melted state, compared to crosslinked or double-stranded overstretched S-DNA.     

Electrophoretic Karyotype and Linkage Groups of the Amoeboflagellate Naegleria gruberi

ABSTRACT. We have constructed a molecular karyotype for two strains of Naegleria gruberi using pulsed field gel electrophoresis. Each strain has about 23 chromosomes, considerably more than any previous estimate. These chromosomes range in size from 400 kilobasepairs to over 2,000 kilobasepairs. In Naegleria , construction of the DNA karyotype depends on assessment of the anomalous electrophoretic mobility of the circular ribosomal RNA genes. We have determined the chromosomal locations of an identified unique gene (flagellar calmodulin) and four identified multigene families (α- and β-tubulin, actin, ubiquitin), as well as three differentially expressed genes of unknown functions. The ca. 12 actin genes are dispersed over at least seven chromosomes, whereas the majority of the more than eight α-tubulin genes are confined to a single chromosome. The ubiquitin genes are found on five chromosomes in one strain and seven in the other and the β-tubulin genes are on three or four. Our observations provide a foundation for molecular genetic studies in this organism.     

Structure and Dynamics of Double Helical DNA in Torsion Angle Hyperspace: A Molecular Mechanics Approach

Abstract

Analysis of the conformational space populated by the torsion angles and the correlation between the conformational energy and the sequence of DNA are important for fully understanding DNA structure and function. Presence of seven variable torsion angles about single covalent bonds in DNA main chain puts a big challenge for such analysis. We have carried out restrained energy minimization studies for four representative dinucleosides, namely d(ApA):d(TpT), d(CpG):d(CpG), d(GpC):d(GpC) and d(CpA):d(TpG) to determine the energy hyperspace of DNA in context to the values of the torsion angles and the structural properties of the DNA conformations populating the favorable regions of this energy hyperspace. The torsion angles were manipulated by constraining their values at the reference points and then performing energy minimization. The energy minima obtained on the potential energy contour plots mostly correspond to the conformations populated in crystal structures of DNA. Some novel favorable conformations that are not present in crystal structure data are also found. The plots also suggest few low energy routes for conformational transitions or the associated energy barrier heights. Analyses of base pairing and stacking possibility reveal structural changes accompanying these transitions as well as the flexibility of different base steps towards variations in different torsion angles.     

Electrophoretic karyotype and linkage groups of the amoeboflagellate Naegleria gruberi

We have constructed a molecular karyotype for two strains of Naegleria gruberi using pulsed field gel electrophoresis. Each strain has about 23 chromosomes, considerably more than any previous estimate. These chromosomes range in size from 400 kilobasepairs to over 2,000 kilobasepairs. In Naegleria, construction of the DNA karyotype depends on assessment of the anomalous electrophoretic mobility of the circular ribosomal RNA genes. We have determined the chromosomal locations of an identified unique gene (flagellar calmodulin) and four identified multigene families (alpha- and beta-tubulin, actin, ubiquitin), as well as three differentially expressed genes of unknown functions. The ca. 12 actin genes are dispersed over at least seven chromosomes, whereas the majority of the more than eight alpha-tubulin genes are confined to a single chromosome. The ubiquitin genes are found on five chromosomes in one strain and seven in the other and the beta-tubulin genes are on three or four. Our observations provide a foundation for molecular genetic studies in this organism.     

Cyanelle DNA from Cyanophora paradoxa exists in two forms due to intramolecular recombination

The DNA of cyanelles which are described as endosymbiotic cyanobacteria from Cyanophora paradoxa (strain LB555UTEX) is equivalent to 127 kilobasepairs. It is characterized by two inverted repeat segments, 10 kilobasepairs in size, which are separated from each other by long single-copy DNA segments of unequal size. This morphology of the chromosome is also found in the chloroplast DNA of most higher plants and some green algae. The cyanelle DNA exists in two forms of circular molecules which differ only in the orientation of the two single-copy DNA segments relative to each other. This is likely due to intramolecular recombination within the inverted repeat segments.     

THE MORPHOLOGY OF THE SWELLING PROCESS IN RAT LIVER MITOCHONDRIA

Through the use of combined spectrophotometric and electron microscope techniques, large amplitude swelling of rat liver mitochondria has been described as an ordered sequence of ultrastructural transitions. Prior to the actual swelling, mitochondria undergo two major conformational changes: condensed to twisted form and twisted to orthodox form. This sequence is independent of (a) the nature of swelling agents and (b) the time of onset of swelling. Agents that delay the onset of swelling act to increase the duration of the twisted conformation. Agents that prevent extensive swelling hold mitochondria in intermediate conformations. Gross swelling, immediately preceded by a decrease in electron opacity of the matrix, involves the rupture of the outer membrane and expansion of the inner compartment of the mitochondrion.     

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.     

Nanomechanics of Diaminopurine-Substituted DNA

2,6-diaminopurine (DAP) is a nucleobase analog of adenine. When incorporated into double-stranded DNA (dsDNA), it forms three hydrogen bonds with thymine. Rare in nature, DAP substitution alters the physical characteristics of a DNA molecule without sacrificing sequence specificity. Here, we show that in addition to stabilizing double-strand hybridization, DAP substitution also changes the mechanical and conformational properties of dsDNA. Thermal melting experiments reveal that DAP substitution raises melting temperatures without diminishing sequence-dependent effects. Using a combination of atomic force microscopy (AFM), magnetic tweezer (MT) nanomechanical assays, and circular dichroism spectroscopy, we demonstrate that DAP substitution increases the flexural rigidity of dsDNA yet also facilitates conformational shifts, which manifest as changes in molecule length. DAP substitution increases both the static and dynamic persistence length of DNA (measured by AFM and MT, respectively). In the static case (AFM), in which tension is not applied to the molecule, the contour length of DAP-DNA appears shorter than wild-type (WT)-DNA; under tension (MT), they have similar dynamic contour lengths. At tensions above 60 pN, WT-DNA undergoes characteristic overstretching because of strand separation (tension-induced melting) and spontaneous adoption of a conformation termed S-DNA. Cyclic overstretching and relaxation of WT-DNA at near-zero loading rates typically yields hysteresis, indicative of tension-induced melting; conversely, cyclic stretching of DAP-DNA showed little or no hysteresis, consistent with the adoption of the S-form, similar to what has been reported for GC-rich sequences. However, DAP-DNA overstretching is distinct from GC-rich overstretching in that it happens at a significantly lower tension. In physiological salt conditions, evenly mixed AT/GC DNA typically overstretches around 60 pN. GC-rich sequences overstretch at similar if not slightly higher tensions. Here, we show that DAP-DNA overstretches at 52 pN. In summary, DAP substitution decreases the overall stability of the B-form double helix, biasing toward non-B-form DNA helix conformations at zero tension and facilitating the B-to-S transition at high tension.     

Extension of torsionally stressed DNA by external force.

Metropolis Monte Carlo simulation was used to study the elasticity of torsionally stressed double-helical DNA. Equilibrium distributions of DNA conformations for different values of linking deficit, external force, and ionic conditions were simulated using the discrete wormlike chain model. Ionic conditions were specified in terms of DNA effective diameter, i.e., hard-core radius of the model chain. The simulations show that entropic elasticity of the double helix depends on how much it is twisted. For low amounts of twisting (less than about one turn per twist persistence length) the force versus extension is nearly the same as in the completely torsionally relaxed case. For more twisting than this, the molecule starts to supercoil, and there is an increase in the force needed to realize a given extension. For sufficiently large amounts of twist, the entire chain is plectonemically supercoiled at low extensions; a finite force must be applied to obtain any extension at all in this regime. The simulation results agree well with the results of recent micromanipulation experiments.