BII nucleotides in the B and C forms of natural-sequence polymeric DNA: A new model for the C form of DNA

Experimental detection of conformational transitions between forms of DNA: problems and prospects

Under different conditions, the DNA double helix can take different geometric forms. Of the large number of its conformations, in addition to the "canonical" B form, the A, C, and Z forms are widely known, and the D, Hoogsteen, and X forms are less known. DNA locally takes the A, C, and Z forms in the cell, in complexes with proteins. We compare different methods for detecting non-canonical DNA conformations: X-ray, IR, and Raman spectroscopy, linear and circular dichroism in both the infrared and ultraviolet regions, as well as NMR (measurement of chemical shifts and their anisotropy, scalar and residual dipolar couplings and inter-proton distances from NOESY (nuclear Overhauser effect spectroscopy) data). We discuss the difficulties in applying these methods, the problems of theoretical interpretation of the experimental results, and the prospects for reliable identification of non-canonical DNA conformations.

C-B-A Test of DNA Force Fields

The DNA duplex may be locally strongly bent in complexes with proteins, for example, with polymerases or in a nucleosome. At such bends, the DNA helix is locally in the noncanonical forms A (with a narrow major groove and a large amount of north sugars) or C (with a narrow minor groove and a large share of BII phosphates). To model the formation of such complexes by molecular dynamics methods, the force field is required to reproduce these conformational transitions for a naked DNA. We analyzed the available experimental data on the B-C and B-A transitions under the conditions easily implemented in modeling: in an aqueous NaCl solution. We selected six DNA duplexes which conformations at different salt concentrations are known reliably enough. At low salt concentrations, poly(GC) and poly(A) are in the B-form, classical and slightly shifted to the A-form, respectively. The duplexes ATAT and GGTATACC have a strong and salt concentration dependent bias toward the A-form. The polymers poly(AC) and poly(G) take the C- and A-forms, respectively, at high salt concentrations. The reproduction of the behavior of these oligomers can serve as a test for the balance of interactions between the base stacking and the conformational flexibility of the sugar-phosphate backbone in a DNA force field. We tested the AMBER bsc1 and CHARMM36 force fields and their hybrids, and we failed to reproduce the experiment. In all the force fields, the salt concentration dependence is very weak. The known B-philicity of the AMBER force field proved to result from the B-philicity of its excessively strong base stacking. In the CHARMM force field, the B-form is a result of a fragile balance between the A-philic base stacking (especially for G:C pairs) and the C-philic backbone. Finally, we analyzed some recent simulations of the LacI-, SOX-4-, and Sac7d-DNA complex formation in the framework of the AMBER force field.

C-DNA may facilitate homologous DNA pairing

Recombination-independent homologous pairing represents a prominent yet largely enigmatic feature of chromosome biology. As suggested by studies in the fungus Neurospora crassa, this process may be based on the direct pairing of homologous DNA molecules. Theoretical search for the DNA structures consistent with those genetic results has led to an all-atom model in which the B-DNA conformation of the paired double helices is strongly shifted toward C-DNA. Coincidentally, C-DNA also features a very shallow major groove that could permit initial homologous contacts without atom-atom clashes. The hereby conjectured role of C-DNA in homologous pairing should encourage the efforts to discover its biological functions and may also clarify the mechanism of recombination-independent recognition of DNA homology.

Solid-State NMR Dipolar and Chemical Shift Anisotropy Recoupling Techniques for Structural and Dynamical Studies in Biological Systems

With the development of NMR methodology and technology during the past decades, solid-state NMR (ssNMR) has become a particularly important tool for investigating structure and dynamics at atomic scale in biological systems, where the recoupling techniques play pivotal roles in modern high-resolution MAS NMR. In this review, following a brief introduction on the basic theory of recoupling in ssNMR, we highlight the recent advances in dipolar and chemical shift anisotropy recoupling methods, as well as their applications in structural determination and dynamical characterization at multiple time scales (i.e., fast-, intermediate-, and slow-motion). The performances of these prevalent recoupling techniques are compared and discussed in multiple aspects, together with the representative applications in biomolecules. Given the recent emerging advances in NMR technology, new challenges for recoupling methodology development and potential opportunities for biological systems are also discussed.

Variability of the DNA Backbone Geometry in DNA-Protein Complexes: Experimental Data Analysis

We have analyzed and compared the available experimental data (PDB) on the backbone geometry of the DNA in solution (NMR), in crystals (X-rays), and in complexes with proteins (X-rays and cryo-electron microscopy). The deoxyribose (pseudorotational angle τ₀) and ε/ζ (BI-BII transition in phosphates) flexibilities are practically the same in the four samples. The α/γ mobility is minimal in crystalline DNA: on the histograms, there is one canonical and one noncanonical t/t peak. The α/γ mobility increases in DNA solutions (three more noncanonical peaks) and is maximal in DNA-protein complexes (another additional peak). On a large amount of data, we have confirmed that the three main degrees of freedom of the sugar-phosphate backbone are "orthogonal": changes in any of the angles τ₀, (ζ-ε), and (γ-α) occur, as a rule, at a constant (usually canonical) value of any other. In the DNA-protein complexes, none of the geometrical parameters commonly used to distinguish the A and B forms of DNA, except for Zp and its simpler analog Zp', show an unambiguous correlation with τ₀. Proteins, binding to DNA, in 59% of cases change the local shape of the helix up to the characteristic of the A-form without switching the deoxyribose conformation from south to north. However, we have found simple local characteristics of one nucleotide that correlate with the angles τ₀ and (ζ-ε). These are the angles C3'C1'N* and C4'C3'P(2), respectively. They are orthogonal in DNA-protein complexes exactly as the pair τ₀ and (ζ-ε). Most characteristics of DNA in complexes with proteins are the same in X-ray and in cryo-EM data, except for the histogram for the angle τ₀. We offer a possible explanation for this difference. We also discuss the artifacts on the ε/ζ histogram for DNA in solutions caused by the currently used NMR refinement protocols.

Recombination-independent recognition of DNA homology for meiotic silencing in Neurospora crassa

The pairing of homologous chromosomes represents a critical step of meiosis in nearly all sexually reproducing species. In many organisms, pairing involves chromosomes that remain apparently intact. The mechanistic nature of homology recognition at the basis of such pairing is unknown. Using "meiotic silencing by unpaired DNA" (MSUD) as a model process, we demonstrate the existence of a cardinally different approach to DNA homology recognition in meiosis. The main advantage of MSUD over other experimental systems lies in its ability to identify any relatively short DNA fragment lacking a homologous allelic partner. Here, we show that MSUD does not rely on the canonical mechanism of meiotic recombination, yet it is promoted by REC8, a conserved component of the meiotic cohesion complex. We also show that certain patterns of interspersed homology are recognized as pairable during MSUD. Such patterns need to be colinear and must contain short tracts of sequence identity spaced apart at 21 or 22 base pairs. By using these periodicity values as a guiding parameter in all-atom molecular modeling, we discover that homologous DNA molecules can pair by forming quadruplex-based contacts with an interval of 2.5 helical turns. This process requires right-handed plectonemic coiling and additional conformational changes in the intervening double-helical segments. Our results 1) reconcile genetic and biophysical evidence for the existence of direct homologous double-stranded DNA (dsDNA)-dsDNA pairing, 2) identify a role for this process in initiating RNA interference, and 3) suggest that chromosomes can be cross-matched by a precise mechanism that operates on intact dsDNA molecules.

Ammonium Gemini Surfactants Form Complexes with Model Oligomers of siRNA and dsDNA

Dimeric cationic surfactants (gemini-type) are a group of amphiphilic compounds with potential use in gene therapy as effective carriers for nucleic acid transfection (i.e., siRNA, DNA, and plasmid DNA). Our studies have shown the formation of lipoplexes composed of alkanediyl-α,ω-bis[(oxymethyl)dimethyldodecylammonium] chlorides and selected 21-base-pair nucleic acid (dsDNA and siRNA) oligomers. To examine the structure and physicochemical properties of these systems, optical microscopy, circular dichroism spectroscopy (CD), small-angle X-ray scattering of synchrotron radiation (SR-SAXS), and agarose gel electrophoresis (AGE) were used. The lengths of spacer groups of the studied surfactants had a significant influence on the surfactants’ complexing properties. The lowest charge ratio (p/n) at which stable lipoplexes were observed was 1.5 and the most frequently occurring microstructure of these lipoplexes were cubic and micellar phases for dsDNA and siRNA, respectively. The cytotoxicity tests on HeLa cells indicated the non-toxic concentration of surfactants to be at approximately 10 µM. The dicationic gemini surfactants studied form complexes with siRNA and dsDNA oligomers; however, the complexation process is more effective towards siRNA. Therefore these systems could be applied as transfection systems for therapeutic nucleic acids.

Single-Stranded DNA Oligonucleotides Retain Rise Coordinates Characteristic of Double Helices

The structures of single-stranded DNA oligonucleotides from dimeric to hexameric sequences have been thoroughly investigated. Computations performed at the density functional level of theory including dispersion forces and solvation show that single-stranded helices adopt conformations very close to crystallographic B-DNA, with rise coordinates amounting up to 3.3 Å. Previous results, suggesting that single strands should be shorter than double helices, largely originated from the incompleteness of the adopted basis set. Although sensible deviations with respect to standard B-DNA are predicted, computations indicate that sequences rich in stacked adenines are the most ordered ones, favoring the B-DNA pattern and inducing regular arrangements also on flanking nucleobases. Several structural properties of double helices rich in adenine are indeed already reflected by the corresponding single strands.

Comparative analysis of inosine-substituted duplex DNA by circular dichroism and X-ray crystallography

Leveraging structural biology tools, we report the results of experiments seeking to determine if the different mechanical properties of DNA polymers with base analog substitutions can be attributed, at least in part, to induced changes from classical B-form DNA. The underlying hypothesis is that different inherent bending and twisting flexibilities may characterize non-canonical B-DNA, so that it is inappropriate to interpret mechanical changes caused by base analog substitution as resulting simply from 'electrostatic' or 'base stacking' influences without considering the larger context of altered helical geometry. Circular dichroism spectra of inosine-substituted oligonucleotides and longer base-substituted DNAs in solution indicated non-canonical helical conformations, with the degree of deviation from a standard B-form geometry depending on the number of I⋅C pairs. X-ray diffraction of a highly inosine-substituted DNA decamer crystal (eight I⋅C and two A⋅T pairs) revealed an A-tract-like conformation with a uniformly narrow minor groove, reduced helical rise, and the majority of sugars adopting a C1'-exo (southeastern) conformation. This contrasts with the standard B-DNA geometry with C2'-endo sugar puckers (south conformation). In contrast, the crystal structure of a decamer with only four I⋅C pairs has a geometry similar to that of the reference duplex with eight G⋅C and two A⋅T pairs. The unique crystal geometry of the inosine-rich duplex is noteworthy given its unusual CD signature in solution and the altered mechanical properties of some inosine-containing DNAs.

Microscopic structure of the polymer-induced liquid precursor for calcium carbonate

Many biomineral crystals form complex non-equilibrium shapes, often via transient amorphous precursors. Also in vitro crystals can be grown with non-equilibrium morphologies, such as thin films or nanorods. In many cases this involves charged polymeric additives that form a polymer-induced liquid precursor (PILP). Here, we investigate the CaCO3 based PILP process with a variety of techniques including cryoTEM and NMR. The initial products are 30-50 nm amorphous calcium carbonate (ACC) nanoparticles with ~2 nm nanoparticulate texture. We show the polymers strongly interact with ACC in the early stages, and become excluded during crystallization, with no liquid-liquid phase separation detected during the process. Our results suggest that "PILP" is actually a polymer-driven assembly of ACC clusters, and that its liquid-like behavior at the macroscopic level is due to the small size and surface properties of the assemblies. We propose that a similar biopolymer-stabilized nanogranular phase may be active in biomineralization.

The study of complexation between dicationic surfactants and the DNA duplex using structural and spectroscopic methods

Dicationic (also known as gemini or dimeric) bis-alkylimidazolium surfactants belong to a group of non-viral transfection systems proposed for the successful introduction of different types of nucleic acids (i.e., siRNA, DNA oligomers, and plasmid DNA) into living cells.

The intrinsic mechanics of B-DNA in solution characterized by NMR

Experimental characterization of the structural couplings in free B-DNA in solution has been elusive, because of subtle effects that are challenging to tackle. Here, the exploitation of the NMR measurements collected on four dodecamers containing a substantial set of dinucleotide sequences provides new, consistent correlations revealing the DNA intrinsic mechanics. The difference between two successive residual dipolar couplings (ΔRDCs) involving C6/8-H6/8, C3′-H3′ and C4′-H4′ vectors are correlated to the ³¹P chemical shifts (δP), which reflect the populations of the BI and BII backbone states. The δPs are also correlated to the internucleotide distances (D_inter) involving H6/8, H2′ and H2″ protons. Calculations of NMR quantities on high resolution X-ray structures and controlled models of DNA enable to interpret these couplings: the studied ΔRDCs depend mostly on roll, while D_inter are mainly sensitive to twist or slide. Overall, these relations demonstrate how δP measurements inform on key inter base parameters, in addition to probe the BI↔BII backbone equilibrium, and shed new light into coordinated motions of phosphate groups and bases in free B-DNA in solution. Inspection of the 5′ and 3′ ends of the dodecamers also supplies new information on the fraying events, otherwise neglected.

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.

Lipoplexes of dicationic gemini surfactants with DNA: Structural features of DNA compaction and transfection efficiency

The internal structure of DNA lipoplexes with hydroxyethylated alkylammonium gemini surfactants (GS) with high transfection activity was studied by circular dichroism. It was shown that the efficiency of transfection of HEK293T cells with the pEGFP-N1 circular plasmid was different from zero only in the region of existence of chiral supramolecular DNA-GS complexes and reaches a maximum at concentrations at which the spontaneous aggregation of components is observed.

Simulations Meet Experiment to Reveal New Insights into DNA Intrinsic Mechanics

The accurate prediction of the structure and dynamics of DNA remains a major challenge in computational biology due to the dearth of precise experimental information on DNA free in solution and limitations in the DNA force-fields underpinning the simulations. A new generation of force-fields has been developed to better represent the sequence-dependent B-DNA intrinsic mechanics, in particular with respect to the BI ↔ BII backbone equilibrium, which is essential to understand the B-DNA properties. Here, the performance of MD simulations with the newly updated force-fields Parmbsc0εζOLI and CHARMM36 was tested against a large ensemble of recent NMR data collected on four DNA dodecamers involved in nucleosome positioning. We find impressive progress towards a coherent, realistic representation of B-DNA in solution, despite residual shortcomings. This improved representation allows new and deeper interpretation of the experimental observables, including regarding the behavior of facing phosphate groups in complementary dinucleotides, and their modulation by the sequence. It also provides the opportunity to extensively revisit and refine the coupling between backbone states and inter base pair parameters, which emerges as a common theme across all the complementary dinucleotides. In sum, the global agreement between simulations and experiment reveals new aspects of intrinsic DNA mechanics, a key component of DNA-protein recognition.

DNA Backbone BI/BII Distribution and Dynamics in E2 Protein-Bound Environment Determined by Molecular Dynamics Simulations

BI and BII conformational substates in the DNA backbone typify canonical B-form DNA. The BI and BII substates are important for structural variation of DNA and have been implicated in protein-nucleic acid recognition mechanisms. Recent refinements have been made to nucleic acid force fields employed in molecular dynamics simulations that demonstrate a better ability to model the BI and BII states, leading to overall improved modeling of DNA structure and dynamics. These force field improvements have yet to be significantly demonstrated in the context of a protein-DNA system extended to long time scales. Our plan was to run molecular dynamics simulations of a well-studied protein-DNA system (E2-DNA) into the microsecond time scale and determine the ability of the force field to populate BII states in the DNA backbone consistent with dinucleotide steps crystallized in the BII conformation. The results showed that the dinucleotide steps in the E2-DNA complex with the highest BII populations from simulation trajectories corresponded to the dinucleotide steps crystallized in the BII state and that decoy BI and BII states converge to the same results within approximately one microsecond.

Solid-state NMR studies of nucleic acid components

Recent applications of solid-state NMR spectroscopy to studies of nucleic acids and their components.

NMR studies of DNA support the role of pre-existing minor groove variations in nucleosome indirect readout

We investigated how the intrinsic sequence-dependent properties probed via the phosphate linkages (BI ↔ BII equilibrium) influence the preferred shape of free DNA, and how this affects the nucleosome formation. First, this exploits NMR solution studies of four B-DNA dodecamers that together cover 39 base pairs of the 5' half of the sequence 601, of special interest for nucleosome formation. The results validate our previous prediction of a systematic, general sequence effect on the intrinsic backbone BII propensities. NMR provides new evidence that the backbone behavior is intimately coupled to the minor groove width. Second, application of the backbone behavior predictions to the full sequence 601 and other relevant sequences demonstrates that alternation of intrinsic low and high BII propensities, coupled to intrinsic narrow and wide minor grooves, largely coincides with the sinusoidal variations of the DNA minor groove width observed in crystallographic structures of the nucleosome. This correspondence is much poorer with low affinity sequences. Overall, the results indicate that nucleosome formation involves an indirect readout process implicating pre-existing DNA minor groove conformations. It also illustrates how the prediction of the intrinsic structural DNA behavior offers a powerful framework to gain explanatory insight on how proteins read DNA.

DNA conformational transitions induced by supercoiling control transcription in chromatin

Regulation of transcription in eukaryotes is considered in the light of recent findings demonstrating the presence of negative and positive superhelical tension in chromatin. This tension induces conformational transitions in DNA duplex. Particularly, the transition into A-form renders DNA accessible and waylaying for initiation of transcription producing RNA molecules long known to belong to the A-conformation. Competition between conformational transitions in various DNA sequences for the energy of elastic spring opens a possibility for understanding of fine tuning of transcription at a distance.

Complete chemical shift assignment of the ssDNA in the filamentous bacteriophage fd reports on its conformation and on its interface with the capsid shell

The fd bacteriophage is a filamentous virus consisting of a circular single-stranded DNA (ssDNA) wrapped by thousands of copies of a major coat protein subunit (the capsid). The coat protein subunits are mostly α-helical and curved, and are arranged in the capsid in consecutive pentamers related by a translation along the main viral axis and a rotation of ~36° (C5S2 symmetry). The DNA is right-handed and helical, but information on its structure and on its interface with the capsid is incomplete. We present here an approach for assigning the DNA nucleotides and studying its interactions with the capsid by magic-angle spinning solid-state NMR. Capsid contacts with the ssDNA are obtained using a two-dimensional (13)C-(13)C correlation experiment and a proton-mediated (31)P-(13)C polarization transfer experiment, both acquired on an aromatic-unlabeled phage sample. Our results allow us to map the residues that face the interior of the capsid and to show that the ssDNA-capsid interactions are sustained mainly by electrostatic interactions between the positively charged lysine side chains and the phosphate backbone. The use of natural abundance aromatic amino acids in the growth media facilitated the complete assignment of the four nucleotides and the observation of internucleotide contacts. Using chemical shift analysis, our study shows that structural features of the deoxyribose carbons reporting on the sugar pucker are strikingly similar to those observed recently for the Pf1 phage. However, the ssDNA-protein interface is different, and chemical shift markers of base pairing are different. This experimental approach can be utilized in other filamentous and icosahedral bacteriophages, and also in other biomolecular complexes involving structurally and functionally important DNA-protein interactions.

5-Methylation of cytosine in CG:CG base-pair steps: a physicochemical mechanism for the epigenetic control of DNA nanomechanics

van der Waals density functional theory is integrated with analysis of a non-redundant set of protein-DNA crystal structures from the Nucleic Acid Database to study the stacking energetics of CG:CG base-pair steps, specifically the role of cytosine 5-methylation. Principal component analysis of the steps reveals the dominant collective motions to correspond to a tensile "opening" mode and two shear "sliding" and "tearing" modes in the orthogonal plane. The stacking interactions of the methyl groups globally inhibit CG:CG step overtwisting while simultaneously softening the modes locally via potential energy modulations that create metastable states. Additionally, the indirect effects of the methyl groups on possible base-pair steps neighboring CG:CG are observed to be of comparable importance to their direct effects on CG:CG. The results have implications for the epigenetic control of DNA mechanics.

Simulating DNA by molecular dynamics: aims, methods, and validation

The structure and dynamics of the B-DNA double helix involves subtle sequence-dependent effects which are decisive for its function, but difficult to characterize. These structural and dynamic effects can be addressed by simulations of DNA sequences in explicit solvent. Here, we present and discuss the state-of-art of B-DNA molecular dynamics simulations with the major force fields in use today. We explain why a critical analysis of the MD trajectories is required to assess their reliability, and estimate the value and limitations of these models. Overall, simulations of DNA bear great promise towards deciphering the structural and physical subtleties of this biopolymer, where much remains to be understood.

Mg2+ in the major groove modulates B-DNA structure and dynamics

This study investigates the effect of Mg(2+) bound to the DNA major groove on DNA structure and dynamics. The analysis of a comprehensive dataset of B-DNA crystallographic structures shows that divalent cations are preferentially located in the DNA major groove where they interact with successive bases of (A/G)pG and the phosphate group of 5'-CpA or TpG. Based on this knowledge, molecular dynamics simulations were carried out on a DNA oligomer without or with Mg(2+) close to an ApG step. These simulations showed that the hydrated Mg(2+) forms a stable intra-strand cross-link between the two purines in solution. ApG generates an electrostatic potential in the major groove that is particularly attractive for cations; its intrinsic conformation is well-adapted to the formation of water-mediated hydrogen bonds with Mg(2+). The binding of Mg(2+) modulates the behavior of the 5'-neighboring step by increasing the BII (ε-ζ>0°) population of its phosphate group. Additional electrostatic interactions between the 5'-phosphate group and Mg(2+) strengthen both the DNA-cation binding and the BII character of the 5'-step. Cation binding in the major groove may therefore locally influence the DNA conformational landscape, suggesting a possible avenue for better understanding how strong DNA distortions can be stabilized in protein-DNA complexes.

DNA information: from digital code to analogue structure

The digital linear coding carried by the base pairs in the DNA double helix is now known to have an important component that acts by altering, along its length, the natural shape and stiffness of the molecule. In this way, one region of DNA is structurally distinguished from another, constituting an additional form of encoded information manifest in three-dimensional space. These shape and stiffness variations help in guiding and facilitating the DNA during its three-dimensional spatial interactions. Such interactions with itself allow communication between genes and enhanced wrapping and histone-octamer binding within the nucleosome core particle. Meanwhile, interactions with proteins can have a reduced entropic binding penalty owing to advantageous sequence-dependent bending anisotropy. Sequence periodicity within the DNA, giving a corresponding structural periodicity of shape and stiffness, also influences the supercoiling of the molecule, which, in turn, plays an important facilitating role. In effect, the super-helical density acts as an analogue regulatory mode in contrast to the more commonly acknowledged purely digital mode. Many of these ideas are still poorly understood, and represent a fundamental and outstanding biological question. This review gives an overview of very recent developments, and hopefully identifies promising future lines of enquiry.

Characterization of the geometry and topology of DNA pictured as a discrete collection of atoms

The structural and physical properties of DNA are closely related to its geometry and topology. The classical mathematical treatment of DNA geometry and topology in terms of ideal smooth space curves was not designed to characterize the spatial arrangements of atoms found in high-resolution and simulated double-helical structures. We present here new and rigorous numerical methods for the rapid and accurate assessment of the geometry and topology of double-helical DNA structures in terms of the constituent atoms. These methods are well designed for large DNA datasets obtained in detailed numerical simulations or determined experimentally at high-resolution. We illustrate the usefulness of our methodology by applying it to the analysis of three canonical double-helical DNA chains, a 65-bp minicircle obtained in recent molecular dynamics simulations, and a crystallographic array of protein-bound DNA duplexes. Although we focus on fully base-paired DNA structures, our methods can be extended to treat the geometry and topology of melted DNA structures as well as to characterize the folding of arbitrary molecules such as RNA and cyclic peptides.

Study on the adsorption mechanism of DNA with mesoporous silica nanoparticles in aqueous solution

Among the numerous adsorption strategies for DNA adsorption into mesopores, the salt-solution-induced adsorption method has a great application potential in nucleic acids science; thus, it is important to understand the adsorption mechanism. This work demonstrates the mechanistic aspects underlying the adsorption behaviors of DNA with mesoporous silica nanoparticles (MSNs) in aqueous solution. The driving forces for the adsorption process can be categorized into three parts: the shielded electrostatic force, the dehydration effect, and the intermolecular hydrogen bonds. Compared to the adsorption behaviors of DNA with a solid silica nanosphere, we find some unique features for DNA adsorption into the mesopores, such as increasing the salt concentration or decreasing the pH value can promote DNA adsorption into the mesoporous silica. Further analysis indicates that the entrance of DNA into mesopores is probably controlled by the Debye length in solution and DNA can generate direct and indirect hydrogen bonds in the pores with different diameters. The following desorption study depicts that such types of hydrogen bonds result in different energy barriers for the desorption process. In summary, our study depicts the mechanism of DNA adsorption within mesopores in aqueous solution and sets the stage for formulating MSNs as carriers of nucleic acids.

Working the kinks out of nucleosomal DNA

Condensation of DNA in the nucleosome takes advantage of its double-helical architecture. The DNA deforms at sites where the base pairs face the histone octamer. The largest so-called kink-and-slide deformations occur in the vicinity of arginines that penetrate the minor groove. Nucleosome structures formed from the 601 positioning sequence differ subtly from those incorporating an AT-rich human α-satellite DNA. Restraints imposed by the histone arginines on the displacement of base pairs can modulate the sequence-dependent deformability of DNA and potentially contribute to the unique features of the different nucleosomes. Steric barriers mimicking constraints found in the nucleosome induce the simulated large-scale rearrangement of canonical B DNA to kink-and-slide states. The pathway to these states shows nonharmonic behavior consistent with bending profiles inferred from AFM measurements.

Understanding the sequence-dependence of DNA groove dimensions: implications for DNA interactions

Background

The B-DNA major and minor groove dimensions are crucial for DNA-protein interactions. It has long been thought that the groove dimensions depend on the DNA sequence, however this relationship has remained elusive. Here, our aim is to elucidate how the DNA sequence intrinsically shapes the grooves.

Methodology/Principal Findings

The present study is based on the analysis of datasets of free and protein-bound DNA crystal structures, and from a compilation of NMR ³¹P chemical shifts measured on free DNA in solution on a broad range of representative sequences. The ³¹P chemical shifts can be interpreted in terms of the BI↔BII backbone conformations and dynamics. The grooves width and depth of free and protein-bound DNA are found to be clearly related to the BI/BII backbone conformational states. The DNA propensity to undergo BI↔BII backbone transitions is highly sequence-dependent and can be quantified at the dinucleotide level. This dual relationship, between DNA sequence and backbone behavior on one hand, and backbone behavior and groove dimensions on the other hand, allows to decipher the link between DNA sequence and groove dimensions. It also firmly establishes that proteins take advantage of the intrinsic DNA groove properties.

Conclusions/Significance

The study provides a general framework explaining how the DNA sequence shapes the groove dimensions in free and protein-bound DNA, with far-reaching implications for DNA-protein indirect readout in both specific and non specific interactions.

Bacterial melanin interacts with double-stranded DNA with high affinity and may inhibit cell metabolism in vivo

Melanin has been found to interact with a number of molecules including metal ions, antibiotics and proteins. In this study, we showed how melanin from bacteria can interact with double-stranded DNA. Investigation using capillary electrophoresis, various spectroscopic techniques and circular dichroism found that melanin interacts with DNA by intercalating between the base pairs of DNA. And this was further supported by simulating different forms of melanin docking to oligonucleotides. Transmission electron microscopy of recombinant Escherichia coli producing melanin suggested the interaction in vivo. Furthermore, we showed how the cytoplasmic localization of melanin may provide a novel function in inhibiting cellular metabolism using microcalorimetry. The implications of the interaction in prokaryotes and eukaryotes were discussed.

DNA structures from phosphate chemical shifts

For B-DNA, the strong linear correlation observed by nuclear magnetic resonance (NMR) between the (31)P chemical shifts (deltaP) and three recurrent internucleotide distances demonstrates the tight coupling between phosphate motions and helicoidal parameters. It allows to translate deltaP into distance restraints directly exploitable in structural refinement. It even provides a new method for refining DNA oligomers with restraints exclusively inferred from deltaP. Combined with molecular dynamics in explicit solvent, these restraints lead to a structural and dynamical view of the DNA as detailed as that obtained with conventional and more extensive restraints. Tests with the Jun-Fos oligomer show that this deltaP-based strategy can provide a simple and straightforward method to capture DNA properties in solution, from routine NMR experiments on unlabeled samples.

Intrinsic flexibility of B-DNA: the experimental TRX scale

B-DNA flexibility, crucial for DNA–protein recognition, is sequence dependent. Free DNA in solution would in principle be the best reference state to uncover the relation between base sequences and their intrinsic flexibility; however, this has long been hampered by a lack of suitable experimental data. We investigated this relationship by compiling and analyzing a large dataset of NMR ³¹P chemical shifts in solution. These measurements reflect the BI ↔ BII equilibrium in DNA, intimately correlated to helicoidal descriptors of the curvature, winding and groove dimensions. Comparing the ten complementary DNA dinucleotide steps indicates that some steps are much more flexible than others. This malleability is primarily controlled at the dinucleotide level, modulated by the tetranucleotide environment. Our analyses provide an experimental scale called TRX that quantifies the intrinsic flexibility of the ten dinucleotide steps in terms of Twist, Roll, and X-disp (base pair displacement). Applying the TRX scale to DNA sequences optimized for nucleosome formation reveals a 10 base-pair periodic alternation of stiff and flexible regions. Thus, DNA flexibility captured by the TRX scale is relevant to nucleosome formation, suggesting that this scale may be of general interest to better understand protein-DNA recognition.

Web 3DNA--a web server for the analysis, reconstruction, and visualization of three-dimensional nucleic-acid structures

The w3DNA (web 3DNA) server is a user-friendly web-based interface to the 3DNA suite of programs for the analysis, reconstruction, and visualization of three-dimensional (3D) nucleic-acid-containing structures, including their complexes with proteins and other ligands. The server allows the user to determine a wide variety of conformational parameters in a given structure—such as the identities and rigid-body parameters of interacting nucleic-acid bases and base-pair steps, the nucleotides comprising helical fragments, etc. It is also possible to build 3D models of arbitrary nucleotide sequences and helical types, customized single-stranded and double-helical structures with user-defined base-pair parameters and sequences, and models of DNA ‘decorated’ at user-defined sites with proteins and other molecules. The visualization component offers unique, publication-quality representations of nucleic-acid structures, such as ‘block’ images of bases and base pairs and stacking diagrams of interacting nucleotides. The w3DNA web server, located at http://w3dna.rutgers.edu, is free and open to all users with no login requirement.

DNA sequence-directed organization of chromatin: structure-based computational analysis of nucleosome-binding sequences

The folding of DNA on the nucleosome core particle governs many fundamental issues in eukaryotic molecular biology. In this study, an updated set of sequence-dependent empirical "energy" functions, derived from the structures of other protein-bound DNA molecules, is used to investigate the extent to which the architecture of nucleosomal DNA is dictated by its underlying sequence. The potentials are used to estimate the cost of deforming a collection of sequences known to bind or resist uptake in nucleosomes along various left-handed superhelical pathways and to deduce the features of sequence contributing to a particular structural form. The deformation scores reflect the choice of template, the deviations of structural parameters at each step of the nucleosome-bound DNA from their intrinsic values, and the sequence-dependent "deformability" of a given dimer. The correspondence between the computed scores and binding propensities points to a subtle interplay between DNA sequence and nucleosomal folding, e.g., sequences with periodically spaced pyrimidine-purine steps deform at low cost along a kinked template whereas sequences that resist deformation prefer a smoother spatial pathway. Successful prediction of the known settings of some of the best-resolved nucleosome-positioning sequences, however, requires a template with "kink-and-slide" steps like those found in high-resolution nucleosome structures.

The electrostatic origin of low-hydration polymorphism in DNA

In recent years, significant progress has been made towards uncovering the physical mechanisms of low-hydration polymorphism in double-helical DNA. The effect appears to be mechanistically similar in different biological systems, and it is due to the ability of water to form spanning H-bonded networks around biomacromolecules via a quasi-two-dimensional percolation transition. In the case of DNA, disintegration of the spanning H-bonded network leads to electrostatic condensation of DNA strands because, below the percolation threshold, water loses its high dielectric permittivity, whereas the concentration of neutralizing counterions becomes high. In this Concept article arguments propose that this simple electrostatic mechanism represents the universal origin of low-hydration polymorphism in DNA.

Ion dynamics and water percolation effects in DNA polymorphism

The dynamics of ions and water at the surface of DNA are studied by computer simulations in a wide range of hydrations involving the zone of low-hydration polymorphism in DNA. The long-range mobility of ions exhibits a stepwise increase at three distinct hydration levels. The first of them is close to the midpoint of the water percolation transition as well as the midpoint of the transition between A- and B-DNA forms. It coincides with the onset of the dissociation of ion pairs on the DNA surface probably caused by the increase in the water dielectric permittivity due to the appearance of the spanning hydrogen-bonding network. The other two steps are attributed to the formation of percolating water layers on the surface of DNA accompanied by the progressive escape of ions from the DNA surface. The results agree with earlier experimental data and further corroborate the suggested universal mechanism of the low hydration polymorphism in DNA including intraduplex electrostatic condensation close to the water percolation threshold.

Quantification of DNA BI/BII backbone states in solution. Implications for DNA overall structure and recognition

The backbone states of B-DNA influence its helical parameters, groove dimensions, and overall curvature. Therefore, detection and fine characterization of these conformational states are desirable. Using routine NMR experiments on a nonlabeled B-DNA oligomer and analyzing high-resolution X-ray structures, we investigated the relationship between interproton distances and backbone conformational states. The three H2'i-H6/8i+1, H2' 'i-H6/8i+1, and H6/8i-H6/8i+1 sequential distances were found cross-correlated and linearly coupled to epsilon-zeta values in X-ray structures and 31P chemical shifts (deltaP) in NMR that reflect the interconversion between the backbone BI (epsilon-zeta < 0 degrees ) and BII (epsilon-zeta > 0 degrees) states. These relationships provide a detailed check of the NMR data consistency and the possibility to extend the set of restraints for structural refinement through various extrapolations. Furthermore, they allow translation of deltaP in terms of BI/BII ratios. Also, comparison of many published deltaP in solution to crystal data shows that the impact of sequence on the BI/BII propensities is similar in both environments and is therefore an intrinsic and general property of B-DNA. This quantification of the populations of BI and BII is of general interest because these sequence-dependent backbone states act on DNA overall structure, a key feature for DNA-protein-specific recognition.

Influence of the Spacer of Cationic Gemini Amphiphiles on the Hydration of Lipoplexes

The impact of the length of gemini surfactant spacer on complexation and condensation of calf thymus DNA by cationic mixed phospholipid/gemini liposomes was investigated by monitoring the conformational changes of DNA by circular dichroism and the lipid hydration level by the emission characteristics of the fluorescent probe laurdan included in the lipid bilayer. The length of the spacer was shown to influence, on one hand, the hydration level and the organization of the corresponding liposomes and, on the other, the variation of lipid hydration level and the DNA conformation upon complexation. In fact, in correspondence with the longest spacer we observed more hydrated liposomes, probably organized in domains, a higher extent of dehydration promoted by the addition of DNA, and a minor extent of DNA conformational change. The physicochemical features of lipoplexes were shown to depend on the [cationic headgroup]/[DNA single base] ratio.

The DNA structure responds differently to physiological concentrations of K(+) or Na(+)

The influence of monovalent cations on DNA conformation and readout is an open question. This NMR study of DNA with either Na(+) or K(+) at physiological concentrations shows that the nature of the cation affects the (31)P chemical shifts (deltaP) and the sequential distances H2'(i)-H6/8(i+1), H2"(i)-H6/8(i+1), and H6/8(i)-H6/8(i+1). The deltaP and distance variations ascertain that the nature of the cation affects the DNA overall structure, i.e. both the conformational equilibria between the backbone BI (epsilon-zeta <0 degrees ) and BII (epsilon-zeta >0 degrees ) states and the helical parameters, via their strong mechanical coupling. These results reveal that Na(+) and K(+) interactions with DNA are different and sequence-dependent. These ions modulate the overall intrinsic properties of DNA, and possibly its packaging and readout.

Water Clustering and Percolation in Low Hydration DNA Shells

The hydrogen-bonded networks of water at the surface of a model DNA molecule are analyzed. At low hydrations, only small water clusters are attached to the DNA surface, whereas, at high hydrations, it is homogeneously covered by a spanning water network. The spanning water network is formed via a percolation transition at an intermediate hydration number of about 15 water molecules per nucleotide, which is very close to the midpoint of polymorphic transitions between A- and B-forms of the double helix. The percolation transition can occur in both A- and B-DNA hydration shells with nearly identical percolation thresholds. However, the mechanism of the percolation transition in A- and B-DNA is qualitatively different in regard to the roles played by the two opposite grooves of the double helix. Free ions can shift the percolation threshold by preventing some water molecules from hydrogen bond networking. The results corroborate the suggested relationship between water percolation and the low hydration polymorphism in DNA.

Structure of the alpha-homo-DNA:RNA duplex and the function of twist and slide to catalogue nucleic acid duplexes

High-resolution NMR studies of an alpha-homo-DNA:RNA duplex reveal the formation of a right-handed parallel-oriented helix. It differs significantly from a standard A- or B-type helix by a small twist value (26.2 degrees ), which leads to a helical pitch of 13.7 base pairs per helical turn, a negative inclination (-1.78 Angstrom) and a large x displacement (5.90 Angstrom). The rise (3.4 Angstrom) is similar to that found in B-DNA. The solution of this new helix structure has stimulated us to develop a mathematical and geometrical model based on slide and twist parameters to describe nucleic acid duplexes. All existing duplexes can be positioned within this landscape, which can be used to understand the helicalization process.