Dynamics of the B-A transition of DNA double helices

Thermodynamic and kinetic properties of a single base pair in A-DNA and B-DNA

Double stranded DNA can adopt different forms, the so-called A-, B-, and Z-DNA, which play different biological roles. In this work, the thermodynamic and the kinetic parameters for the base-pair closing and opening in A-DNA and B-DNA were calculated by all-atom molecular dynamics simulations at different temperatures. The thermodynamic parameters of the base pair in B-DNA were in good agreement with the experimental results. The free energy barrier of breaking a single base stack results from the enthalpy increase ΔH caused by the disruption of hydrogen bonding and base-stacking interactions, as well as water and base interactions. The free energy barrier of base pair closing comes from the unfavorable entropy loss ΔS caused by the restriction of torsional angles and hydration. It was found that the enthalpy change ΔH and the entropy change ΔS for the base pair in A-DNA are much larger than those in B-DNA, and the transition rates between the opening and the closing state for the base pair in A-DNA are much slower than those in B-DNA. The large difference of the enthalpy and entropy change for forming the base pair in A-DNA and B-DNA results from different hydration in A-DNA and B-DNA. The hydration pattern observed around DNA is an accompanying process for forming the base pair, rather than a follow-up of the conformation.

Kinetics of the Conformational Transformation between B- and A-Forms in the Drew-Dickerson Dodecamer

Some DNA sequences in crystals and in complexes with proteins can exist in the forms intermediate between the B- and A-DNA. Based on this, it was implied that the B-to-A transition for any DNA molecule should go through these intermediate forms also in kinetics. More precisely, the helix parameter Slide has to change first, and the molecule should take the E-form. After that, the Roll parameter changes. In the present work, we simulated the kinetics of the B-A transition in the Drew-Dickerson dodecamer, a known B-philic DNA oligomer. We used the "sugar" coarse-grained model that reproduces ribose flexibility, preserves sequence specificity, employs implicit water and explicit ions, and offers the possibility to vary friction. As the control parameter of the transition, we chose the volume available for a counterion and considered the change from a large to a small volume. In the described system, the B-to-A conformational transformation proved to correspond to a first-order phase transition. The molecule behaves like a small cluster in the region of such a transition, jumping between the A- and B-forms in a wide range of available volumes. The viscosity of the solvent does not affect the midpoint of the transition but only the overall mobility of the system. All helix parameters change synchronously on average, we have not observed the sequence "Slide first, Roll later" in kinetics, and the E-DNA is not a necessary step for the transition between the B- and A-forms in the studied system. So, the existence of the intermediate DNA forms requires specific conditions, shifting the common balance of interactions: certain nucleotide sequence in specific solution or/and the interaction with some protein.

Role of protein interactions in stabilizing canonical DNA features in simulations of DNA in crowded environments

Background

Cellular environments are highly crowded with biological macromolecules resulting in frequent non-specific interactions. While the effect of such crowding on protein structure and dynamics has been studied extensively, very little is known how cellular crowding affects the conformational sampling of nucleic acids.

Results

The effect of protein crowding on the conformational preferences of DNA (deoxyribonucleic acid) is described from fully atomistic molecular dynamics simulations of systems containing a DNA dodecamer surrounded by protein crowders. From the simulations, it was found that DNA structures prefer to stay in B-like conformations in the presence of the crowders. The preference for B-like conformations results from non-specific interactions of crowder proteins with the DNA sugar-phosphate backbone. Moreover, the simulations suggest that the crowder interactions narrow the conformational sampling to canonical regions of the conformational space.

Conclusions

The overall conclusion is that crowding effects may stabilize the canonical features of DNA that are most important for biological function. The results are complementary to a previous study of DNA in reduced dielectric environments where reduced dielectric environments alone led to a conformational shift towards A-DNA. Such a shift was not observed here suggested that the reduced dielectric response of cellular environments is counteracted by non-specific interactions with protein crowders under in vivo conditions.

Electronic supplementary material

The online version of this article (10.1186/s13628-018-0048-y) contains supplementary material, which is available to authorized users.

AMOEBA Polarizable Atomic Multipole Force Field for Nucleic Acids

The AMOEBA polarizable atomic multipole force field for nucleic acids is presented. Valence and electrostatic parameters were determined from high-level quantum mechanical data, including structures, conformational energy, and electrostatic potentials, of nucleotide model compounds. Previously derived parameters for the phosphate group and nucleobases were incorporated. A total of over 35 μs of condensed-phase molecular dynamics simulations of DNA and RNA molecules in aqueous solution and crystal lattice were performed to validate and refine the force field. The solution and/or crystal structures of DNA B-form duplexes, RNA duplexes, and hairpins were captured with an average root-mean-squared deviation from NMR structures below or around 2.0 Å. Structural details, such as base pairing and stacking, sugar puckering, backbone and χ-torsion angles, groove geometries, and crystal packing interfaces, agreed well with NMR and/or X-ray. The interconversion between A- and B-form DNAs was observed in ethanol-water mixtures at 328 K. Crystal lattices of B- and Z-form DNA and A-form RNA were examined with simulations. For the RNA tetraloop, single strand tetramers, and HIV TAR with 29 residues, the simulated conformational states, 3 J-coupling, nuclear Overhauser effect, and residual dipolar coupling data were compared with NMR results. Starting from a totally unstacked/unfolding state, the rCAAU tetranucleotide was folded into A-form-like structures during ∼1 μs molecular dynamics simulations.

Kinetics of the B-A transition of DNA: analysis of potential contributions to a reaction barrier

Because of open problems in the relation between results obtained by relaxation experiments and molecular dynamics simulations on the B-A transition of DNA, relaxation measurements of the B-A dynamics have been extended to a wider range of conditions. Field-induced reaction effects are measured selectively by the magic angle technique using a novel cell construction preventing perturbations from cell window anisotropy. The kinetics was recorded for the case of poly[d(AT)] up to the salt concentration limit of 4.4 mM, where aggregation does not yet interfere. Now experimental data on the B-A dynamics are available for poly[d(AT)] at salt concentrations of 0.18, 0.73, 2.44 and 4.4 mM. In all cases, a spectrum of time constants is found, ranging from ~ 10 μs up to components approaching ~ 1 ms. The relatively small dependence of these data on the salt concentration indicates that electrostatic effects on the kinetics are not as strong as may be expected. The ethanol content at the transition center is a linear function of the logarithm of the salt concentration, and the slope is close to that expected from polyelectrolyte theory. The B-A transition dynamics was also measured in D₂O at a salt concentration of 2.4 mM: the center of the transition is found at 20.0 mol/l H₂O and at 20.1 mol/l D₂O with an estimated accuracy of ± 0.1 mol/l; the spectrum of time constants at the respective transition centers is very similar. The experimental results are discussed regarding the data obtained by molecular dynamics simulations.

Comparison of the Binding Geometry of Free-Base and Hexacoordinated Cationic Porphyrins to A- and B-Form DNA

Although the transition from B-DNA to the A-form is essential for many biological concerns, the properties of this transition have not been resolved. The B to A equilibrium can be analyzed conveniently because of the significant changes in circular dichroism (CD) and absorption spectrum. CD and linear dichroism (LD) methods were used to examine the binding of water-soluble meso-tetrakis(N-methylpyridinium-4-yl)porphyrin (TMPyP) and its derivatives, Co-TMPyP, with B- and A-calf thymus DNA. B- to A-transitions occurred when the physiological buffer was replaced with a water-ethanol mixture (∼80 v/v %), and the fluorescence emission spectra of TMPyP bound to DNA showed a different pattern under ethanol-water conditions and water alone. The featureless broad emission bands of TMPyP were split into two peaks near at 658 and 715 nm in the presence of DNA under an aqueous solution. In the case of an ethanol-water system, however, the emission bands are split in two peaks near at 648 and 708 nm and 656 and 715 nm with and without DNA, respectively. This may be due to a change in the solution polarity. On the basis of the CD and LD data, TMPyP interacts with B-DNA via intercalation at a low ratio under a low ionic strength, 1 mM sodium phosphate. On the other hand, the interaction with A-DNA (80 v/v % ethanol-water system) occurs in a nonintercalating manner. This difference might be because the structural conformations, such as the groove of A-DNA, are not as deep as in B-DNA and the bases are much more tilted. In the case of Co-TMPyP, porphyrin binds preferably via an outside self-stacking mode with B- and A-DNA.

Contraction and Expansion of Stimuli-Responsive DNA Bonds in Flexible Colloidal Crystals

DNA surface ligands can be used as programmable "bonds" to control the arrangement of nanoparticles into crystalline superlattices. Here, we study the intrinsic responsiveness of these DNA bonds to changes in local dielectric constant (εr) as a new approach to dynamically modulate superlattice structure. Remarkably, ethanol (EtOH) addition can be used to controllably tune DNA bond length from 16 to 3 nm and to increase bond stability by >40 °C, while retaining long-range order and crystal habit. Interestingly, we find that these structural changes, which involve the expansion and contraction of crystals by up to 75% in volume, occur in a cooperative fashion once a critical percentage of EtOH is reached. These results provide a facile and robust approach to create stimuli-responsive lattices, to access high volume fractions, and to improve thermal stability.

Boundary conditions for free A-DNA in solution and the relation of local to global DNA structures at reduced water activity

Because of repeated claims that A-DNA cannot exist without aggregation or condensation, the state of DNA restriction fragments with 84–859 bp has been analyzed in aqueous solutions upon reduction of the water activity. Rotational diffusion times τ^d measured by electric dichroism at different water activities with a wide variation of viscosities are normalized to values τ^c at the viscosity of water, which indicate DNA structures at a high sensitivity. For short helices (chain lengths \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\ell} $$\end{document}ℓ ≤ persistence length p), cooperative formation of A-DNA is reflected by the expected reduction of the hydrodynamic length; the transition to the A-form is without aggregation or condensation upon addition of ethanol at monovalent salt ≤1 mM. The aggregation boundary, indicated by a strong increase of τ^c, is shifted to higher monovalent salt (≥4 mM) when ethanol is replaced by trifluoroethanol. The BA transition is not indicated anymore by a cooperative change of τ^c for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\ell} $$\end{document}ℓ » p; τ^c values for these long chains decrease upon reduction of the water activity continuously over the full range, including the BA transition interval. This suggests a non-cooperative BC transition, which induces DNA curvature. The resulting wide distribution of global structures hides changes of local length during the BA transition. Free A-DNA without aggregation/condensation is found at low-salt concentrations where aggregation is inhibited and/or very slow. In an intermediate range of solvent conditions, where the A-form starts to aggregate, a time window remains that can be used for analysis of free A-DNA in a quasi-equilibrium state.

A reversible B–A transition of DNA duplexes induced by synthetic cationic copolymers

Reversible B–A transitions of DNA duplexes were induced by synthetic cationic and anionic polymers.

Single-Walled Carbon Nanotubes Modulate the B- to A-DNA Transition

We study the conformational equilibrium between B-to-A forms of ds-DNA adsorbed onto a single-walled carbon nanotube (SWNT) using free energy profile calculations based on all-atom molecular dynamics simulations. The potential of mean force (PMF) of the B-to-A transition of ds-DNA in the presence of an uncharged (10,0) carbon nanotube for two dodecamers with poly-AT or poly-GC sequences is calculated as a function of a root-mean-square-distance (ΔRMSD) difference metric for the B-to-A transition. The calculations reveal that in the presence of a SWNT DNA favors B-form DNA significantly in both poly-GC and poly-AT sequences. Furthermore, the poly-AT DNA:SWNT complex shows a higher energy penalty for adopting an A-like conformation than poly-GC DNA:SWNT by several kcal/mol. The presence of a SWNT on either poly-AT or poly-GC DNA affects the PMF of the transition such that the B form is favored by as much as 10 kcal/mol. In agreement with published data, we find a potential energy minimum between A and B-form DNA at ΔRMSD ≈ -1.5 Å and that the presence of the SWNT moves this minimum by as much as ΔRMSD = 3 Å.

Conformational preferences of DNA in reduced dielectric environments

The effect of reduced dielectric environments on the conformational sampling of DNA was examined through molecular dynamics simulations. Different dielectric environments were used to model one aspect of cellular environments. Implicit solvent based on the Generalized Born methodology was used to reflect different dielectric environments in the simulations. The simulation results show a tendency of DNA structures to favor noncanonical A-like conformations rather than canonical A- and B-forms as a result of the reduced dielectric environments. The results suggest that the reduced dielectric response in cellular environments may be sufficient to enhance the sampling of A-like DNA structures compared to dilute solvent conditions.

Sequence dependent free energy profiles of localized B- to A-form transition of DNA in water

DNA carries an inherent polymorphism, which surfaces under various external conditions. While B-form remains predominant under normal physiological conditions for most of the DNA sequences, low humidity and increased ion concentration cause B- to A-form transition. Certain proteins and molecules also sometimes cause local deformation of the DNA to the specific A-form. Previous experimental and computational studies focused on the overall B- to A-form transition. Here for the first time we investigated thermodynamics and mechanism of B- to A-form transition in water for various DNA sequences at a local dinucleotide base pair level. We introduced a new reaction coordinate Zp', based on the unique order parameter Zp, to drive B- to A-form transition locally and thereby calculate free energy profiles for the same for all the ten different dinucleotide steps embedded in a twelve base pair DNA. Results show that the trend of "A" and "B" philicity observed in experiment is preserved even at this local dinucleotide level, indicating its localized origin. Higher free energy cost obtained here is attributed to the cost of creating B∕A junctions along with formation of B->A transition at dimer level. We find that while water energetically stabilizes A-form for all the ten different dinucleotide steps to various extents, entropy acts against it. Therefore, we find that the stability of B-form DNA in water is entropic in origin. Mechanism of the conversion appears to be triggered by Slide; however, backbone parameters change concertedly.

DNA manipulation by means of insulator-based dielectrophoresis employing direct current electric fields

Electrokinetic techniques offer a great potential for biological particle manipulation. Among these, dielectrophoresis (DEP) has been successfully utilized for the concentration of bioparticles. Traditionally, DEP is performed employing microelectrodes, an approach with attractive characteristics but expensive due to microelectrode fabrication costs. An alternative is insulator-based DEP, a method where non-uniform electric fields are created with arrays of insulating structures. This study presents the concentration of linear DNA particles (pET28b) employing a microchannel, with an array of cylindrical insulating structures and direct current electric fields. Results showed manipulation of DNA particles with a combination of electroosmotic, electrophoretic, and dielectrophoretic forces. Employing suspending media with conductivity of 104 muS/cm and pH of 11.15, under applied fields between 500 and 1500 V/cm, DNA particles were observed to be immobilized due to negative dielectrophoretic trapping. The observation of DNA aggregates that occurred at higher applied fields, and dispersed once the field was removed is also included. Finally, concentration factors varying from 8 to 24 times the feed concentration were measured at 2000 V/cm after concentration time-periods of 20-40 s. The results presented here demonstrate the potential of insulator-based DEP for DNA concentration, and open the possibility for fast DNA manipulation for laboratory and large-scale applications.

Phase transition between A and B forms of DNA: a free-energy perspective

We study the structural transition from B form of DNA to A form of DNA using group theoretic methods. The transition is not of the order-disorder type and hence to construct a Landau kind of theory for the transition we define a higher symmetry and relevant order parameters. We also discuss the issue of all the conformations, observed experimentally during the course of transition, being fundamentally different or not.

Lesion-induced DNA weak structural changes detected by pulsed EPR spectroscopy combined with site-directed spin labelling

Double electron-electron resonance (DEER) was applied to determine nanometre spin-spin distances on DNA duplexes that contain selected structural alterations. The present approach to evaluate the structural features of DNA damages is thus related to the interspin distance changes, as well as to the flexibility of the overall structure deduced from the distance distribution. A set of site-directed nitroxide-labelled double-stranded DNA fragments containing defined lesions, namely an 8-oxoguanine, an abasic site or abasic site analogues, a nick, a gap and a bulge structure were prepared and then analysed by the DEER spectroscopic technique. New insights into the application of 4-pulse DEER sequence are also provided, in particular with respect to the spin probes' positions and the rigidity of selected systems. The lesion-induced conformational changes observed, which were supported by molecular dynamics studies, confirm the results obtained by other, more conventional, spectroscopic techniques. Thus, the experimental approaches described herein provide an efficient method for probing lesion-induced structural changes of nucleic acids.

Spectroscopic methods for the physical characterization and formulation of nonviral gene delivery systems

Currently, with the exception of naked DNA formulations, most pharmaceutical preparations of plasmid DNA employ some type of polycationic delivery vector such as synthetic cationic polymers and lipids to enhance delivery. A number of biophysical techniques are readily available for the structural characterization of plasmid DNA within synthetic gene delivery complexes. Here we present applications of ultraviolet (UV) absorption, circular dichroism (CD), infrared (IR), and fluorescence spectroscopies as well as dynamic light scattering to the structural analysis of the oligonucleotide component of nonviral gene delivery vectors. We also illustrate this approach for the investigation of the formulation of lipoplex and polyplex-based gene delivery systems. To summarize such data, we show how the macromolecular complexes can be represented as vectors in a highly dimensional space in which the components of the vector consist of normalized values of experimental parameters measured as a function of different solution conditions such as pH and ionic strength.

Spectroscopic and molecular dynamics evidence for a sequential mechanism for the A-to-B transition in DNA

The A-to-B form transition has been examined in three DNA duplexes, d(CGCGAATTCGCG)(2), d(CGCGAATTGCGC), and d(CGCAAATTTCGC), using circular dichroism spectroscopy, ultraviolet resonance Raman (UVRR) spectroscopy, and molecular dynamics (MD) simulation. Circular dichroism spectra confirm that these molecules adopt the A form under conditions of reduced water activity. UVRR results, obtained under similar conditions, suggest that the transition involves a series of intermediate forms between A and B. Cooperative and distinct transitions were observed for the bases and the sugars. Independent MD simulations on d(CGCGAATTCGCG)(2) show a spontaneous change from the A to B form in aqueous solution and describe a kinetic model that agrees well with UVRR results. Based on these observations, we predict that the mechanism of the transition involves a series of A/B hybrid forms and is sequential in nature, similar to previous crystallographic studies of derivatized duplexes. A simulation in which waters were restrained in the major groove of B DNA shows a rapid, spontaneous change from B to A at reduced water activity. These results indicate that a quasiergodic sampling of the solvent distribution may be a problem in going from B to A at reduced water activity in the course of an MD simulation.

A Novel Major Groove Binding Site in B-form DNA for Ethidium Cation

A stably-bound external binding site for ethidium cation in the major groove of B-form DNA is proposed. This complex is stabilized by hydrogen bonding between this ligand and the nucleophilic centers O6 and N7 of guanine, both of which are accessible via the major groove. This binding site is not the same as the well-characterized electrostatically-stabilized external binding site, but rather is seen to be a covalently bound complex which is stabilized by two hydrogen bonds between the ethidium ligand and guanine in the double stranded (ds) B-form DNA. This site [(1), R. Monaco, F. Hasheer. J Biomol Struct Dyn 10, 675 (1993)] can only exist at very low occupancy ratios. The existence of this binding site leads directly to the expectation that there will exist particular mechanistic steps along the pathway of interaction between ethidium and ds B-DNA at low and high ligand concentrations that involve this binding mode. This would not only explain observations published recently [for example, see (2-6), W. Wilson, I. Lopp. Biopolymers 18, 3025 (1979); L. Wakelin, M. Waring. J Mol Biol 144, 183-214 (1980); A. Karpetyan, N. Mehrabian, G. Terzikian, A. Antonian, P. Vardevanian, M. Frank-Kamenetshii. Proceedings of the 10th Conversation, SUNY Albany, 275 (1998); P. Vardevanyan, A. Antonyan, G. Manukyan, A. Karapetyan. Experimental and Molecular Medicine 33, 205 (2001); P. Vardevanyan, A. Antonyan, L. Minasbekan, A. Karapetyan. Proceedings of the 2002 Miami Nature Biotechnology Winter Symposium, 2(S1), 144 (2002)] but also give insight into discrepancies reported in the literature over the years by different workers studying the mechanism of interaction between ethidium and DNA. In this paper this novel binding interaction is discussed, and it is shown how the elucidation of this interaction led to the proposal of two distinct mechanisms of intercalation between ds B-DNA and ethidium cation for high and low concentrations of ligand. Modeling studies show the stability, configuration, and relative energies of this outside binding site. It is expected that this externally bound complex between ethidium cation and ds B-form DNA will be experimentally detectable using fluorescent polarization and/or linear and circular dichroism spectroscopic studies [(7, 8) E. Tuite, U. Sehlstedt, P. Hagmar, B. Norden, M. Takahashi. Euro J Biochem 243, 482-492 (1997); T. Hard. Biopolymers 26, 613-618 (1987)].

Simulation of adsorption of DNA on carbon nanotubes

We report molecular dynamics simulations of DNA adsorption on a single-walled carbon nanotube (SWNT) in an aqueous environment. We have modeled a DNA segment with 12 base pairs (Dickerson dodecamer) and a (8,8) SWNT in water, with counterions to maintain total charge neutrality. Simulations show that DNA binds to the external surface of an uncharged or positively charged SWNT on a time scale of a few hundred picoseconds. The hydrophobic end groups of DNA are attracted to the hydrophobic SWNT surface of uncharged SWNTs, while the hydrophilic backbone of DNA does not bind to the uncharged SWNT. The binding mode of DNA to charged SWNTs is qualitatively different from uncharged SWNTs. The phosphodiester groups of the DNA backbone are attracted to a positively charged SWNT surface while DNA does not adsorb on negatively charged SWNTs. There is no evidence for canonical double-stranded DNA wrapping around either charged or uncharged SWNTs on the very short time scales of the simulations. The adsorption process appears to have negligible effect on the internal stacking structure of the DNA molecule but significantly affects the A to B form conversion of A-DNA. The adsorption of A-DNA onto an uncharged SWNT inhibits the complete relaxation of A-DNA to B-DNA within the time scale of the simulations. In contrast, binding of the A-DNA onto a positively charged SWNT may promote slightly the A to B conversion.

Theoretical study of large conformational transitions in DNA: the B<-->A conformational change in water and ethanol/water

We explore here the possibility of determining theoretically the free energy change associated with large conformational transitions in DNA, like the solvent-induced B<-->A conformational change. We find that a combination of targeted molecular dynamics (tMD) and the weighted histogram analysis method (WHAM) can be used to trace this transition in both water and ethanol/water mixture. The pathway of the transition in the A-->B direction mirrors the B-->A pathway, and is dominated by two processes that occur somewhat independently: local changes in sugar puckering and global rearrangements (particularly twist and roll) in the structure. The B-->A transition is found to be a quasi-harmonic process, which follows closely the first spontaneous deformation mode of B-DNA, showing that a physiologically-relevant deformation is in coded in the flexibility pattern of DNA.

The nature of "unusual" electro-optical transients observed for DNA

Unusual electro-optical transients have been observed for many different polymers and colloidal systems. These effects provoked serious confusion, because a simple-minded interpretation can be completely misleading. The case of double helical DNA is of particular interest, because DNA has been studied in more detail than other systems and because of its biological function. DNA is subject to bending, which implies a loss of symmetry. Due to its high charge density, non-symmetric conformations must have a non-symmetric distribution of charges leading to a torque of considerable magnitude in the presence of external electric fields. The dipole moment describing this torque must be calculated in a coordinate system with its origin at the center of diffusion. The resulting dipole values are in the range of thousands of Debye units. Because the new dipole type is analogous to but not identical with permanent dipoles, the notation "quasi-permanent" dipole is suggested. Application of this concept, using commonly accepted parameters for DNA and established procedures for calculation of electro-optical transients, leads to "unusual" transients. Thus, these transients must be expected from well-known parameters of DNA double helices. The influence of the quasi-permanent dipole moment may be amplified considerably by hydrodynamic coupling. This effect has been demonstrated for the case of smoothly bent rods. Both model calculations and experiments illustrate the danger of getting data that may be completely misleading. For example, depending on pulse amplitudes and/or pulse lengths, electro-optical decays may be accelerated artificially due to superposition of decay components with opposite amplitudes. Experiments show that unusual transients and apparent acceleration effects disappear, when high frequency sine pulses are used for the electro-optical analysis of DNA. Electro-optical effects depend upon the internal dynamics of the object under investigation. In general, the dynamics of DNA bending was assumed to be fast compared to rotational diffusion. Because stacking rearrangements in single stranded nucleic acids are relatively slow and recently the dynamics of the B-A transition was observed in the time range >1 micros, it is likely that there are also relatively slow rearrangements between bending conformers. Bending transitions are expected to be relatively fast, when there are no activation barriers in the bending pathway, and may be slow, when activation barriers must be passed between bending conformers.

Carbon nanotubes selective destabilization of duplex and triplex DNA and inducing B-A transition in solution

Single-walled carbon nanotubes (SWNTs) have been considered as the leading candidate for nanodevice applications ranging from gene therapy and novel drug delivery to membrane separations. The miniaturization of DNA-nanotube devices for biological applications requires fully understanding DNA-nanotube interaction mechanism. We report here, for the first time, that DNA destabilization and conformational transition induced by SWNTs are sequence-dependent. Contrasting changes for SWNTs binding to poly[dGdC]:poly[dGdC] and poly[dAdT]:poly[dAdT] were observed. For GC homopolymer, DNA melting temperature was decreased 40°C by SWNTs but no change for AT-DNA. SWNTs can induce B–A transition for GC-DNA but AT-DNA resisted the transition. Our circular dichroism, competitive binding assay and triplex destabilization studies provide direct evidence that SWNTs induce DNA B–A transition in solution and they bind to the DNA major groove with GC preference.

The dynamics of the B-A transition of natural DNA double helices

The dynamics of the B-A transition of DNA double helices with different GC contents and various chain lengths has been characterized by an electric field pulse technique. The field-induced B-A reaction is separated from orientation effects using the magic angle technique. Amplitudes reflecting the B-A reaction are observed selectively in the limited range of ethanol contents, where CD spectra demonstrate the B-A transition. The maximum amplitude appears at 1-2% higher ethanol content than the center of the B-A transition observed by CD because electric field pulses induce a relatively large perturbation from the A- toward the B-form. The relaxation curves measured after pulse termination reflect a spectrum of up to three relaxation processes. For DNA's with approximately 50% GC, the main part of the amplitude ( approximately 75%) is associated with time constants of approximately 2 micros, and another major component appears with time constants of 50-100 micros. These relaxation effects have been observed for DNA samples with 859, 2629, 7160, and 48501 bp. The time constant associated with the main amplitude increases with decreasing GC content from approximately 2 micros at 50% GC to approximately 3 mus at 41% GC and approximately 10 micros at 0% GC at the center of the B-A transition. Model calculations on the kinetics of cooperative linear Ising lattices predict the appearance of a distinct maximum of the mean relaxation time at the center of the transition. The absence of such maximum in our experimental data indicates a low cooperativity of the B-A transition with a nucleation parameter of approximately 0.1. The rate of the B-A transition is lower by approximately 3 orders of magnitude than that predicted by molecular dynamics simulations.

The physical determinants of the DNA conformational landscape: an analysis of the potential energy surface of single-strand dinucleotides in the conformational space of duplex DNA

A multivariate analysis of the backbone and sugar torsion angles of dinucleotide fragments was used to construct a 3D principal conformational subspace (PCS) of DNA duplex crystal structures. The potential energy surface (PES) within the PCS was mapped for a single-strand dinucleotide model using an empirical energy function. The low energy regions of the surface encompass known DNA forms and also identify previously unclassified conformers. The physical determinants of the conformational landscape are found to be predominantly steric interactions within the dinucleotide backbone, with medium-dependent backbone-base electrostatic interactions serving to tune the relative stability of the different local energy minima. The fidelity of the PES to duplex DNA properties is validated through a correspondence to the conformational distribution of duplex DNA crystal structures and the reproduction of observed sequence specific propensities for the formation of A-form DNA. The utility of the PES is demonstrated through its succinct and accurate description of complex conformational processes in simulations of duplex DNA. The study suggests that stereochemical considerations of the nucleic acid backbone play a role in determining conformational preferences of DNA which is analogous to the role of local steric interactions in determining polypeptide secondary structure.

Thermal difference spectra: a specific signature for nucleic acid structures

We show that nucleic acid structures may be conveniently and inexpensively characterized by their UV thermal difference spectra. A thermal difference spectrum (TDS) is obtained for a nucleic acid by simply recording the ultraviolet absorbance spectra of the unfolded and folded states at temperatures above and below its melting temperature (T_m). The difference between these two spectra is the TDS. The TDS has a specific shape that is unique for each type of nucleic acid structure, a conclusion that is based on a comparison of >900 spectra from 200 different sequences. The shape of the TDS reflects the subtleties of base stacking interactions that occur uniquely within each type of nucleic acid structure. TDS provides a simple, inexpensive and rapid method to obtain structural insight into nucleic acid structures, which is applicable to both DNA and RNA from short oligomers to polynucleotides. TDS complements circular dichroism as a tool for the structural characterization of nucleic acids in solution.

Free energy landscape of A-DNA to B-DNA conversion in aqueous solution

The interconversion between the well-characterized A- and B-forms of DNA is a structural transition for which the intermediate states and the free energy difference between the two endpoints are not known precisely. In the present study, the difference between the Root Mean Square Distance (RMSD) from canonical A-form and B-form DNA is used as an order parameter to characterize this free energy difference using umbrella sampling molecular dynamics (MD) simulations with explicit solvent. The constraint imposed along this order parameter allows relatively unrestricted evolution of the intermediate structures away from both canonical A- and B-forms. The free energy difference between the A- and B-forms for the hexamer DNA sequence CTCGAG in aqueous solution is conservatively estimated to be at least 2.8 kcal/mol. A continuum of intermediate structures with no well-defined local minima links the two forms. The absence of any major barriers in the free energy surface is consistent with spontaneous conversion of the A-form DNA to B-form DNA in unconstrained simulations. The extensive sampling in the MD simulations (>0.1 mus) also allowed quantitative energetic characterization of local backbone conformational variables such as sugar pseudorotation angles and BI/BII state equilibria and their dependence on base identity. The absolute minimum in the calculated free energy profile corresponds closely to the crystal structure of the hexamer sequence, indicating that the present method has the potential to identify the most stable state for an arbitrary DNA sequence in water.

The B- to A-DNA transition and the reorganization of solvent at the DNA surface

DNA geometry depends on relative humidity. Using the CHARMM22 force field to push B-DNA to A-DNA, a molecular dynamics simulation of a mixed-sequence 24-basepair DNA double-stranded oligomer, starting from B-DNA, was carried out to explore both the mechanism of the transition and the evolution of hydration patterns on the surface of DNA. Over the 11-ns trajectory, the transition recapitulates the slide-first, roll-later mechanism, is opposed by DNA electrostatics, and is favored by an increasing amount of condensed sodium ions. Hydration was characterized by counting the hydrogen bonds between water and DNA, and by the number of water bridges linking two DNA atoms. The number of hydrogen bonds between water and DNA remains constant during the transition, but there is a 40% increase in the number of water bridges, in agreement with the principle of economy of hydration. Water bridges emerge as delicate sensors of both structure and dynamics of DNA. Both local flexibility and the frustration of the water network on the surface of DNA probably account for the low populations and short residence times of the bridges, and for the lubricant role of water in ligand-DNA interactions.