A unified model for the origin of DNA sequence-directed curvature

DNA polymerase ε and δ variants drive mutagenesis in polypurine tracts in human tumors

Alterations in the exonuclease domain of DNA polymerase ε cause ultramutated cancers. These cancers accumulate AGA>ATA transversions; however, their genomic features beyond the trinucleotide motifs are obscure. We analyze the extended DNA context of ultramutation using whole-exome sequencing data from 524 endometrial and 395 colorectal tumors. We find that G>T transversions in POLE-mutant tumors predominantly affect sequences containing at least six consecutive purines, with a striking preference for certain positions within polypurine tracts. Using this signature, we develop a machine-learning classifier to identify tumors with hitherto unknown POLE drivers and validate two drivers, POLE-E978G and POLE-S461L, by functional assays in yeast. Unlike other pathogenic variants, the E978G substitution affects the polymerase domain of Pol ε. We further show that tumors with POLD1 drivers share the extended signature of POLE ultramutation. These findings expand the understanding of ultramutation mechanisms and highlight peculiar mutagenic properties of polypurine tracts in the human genome.

DNA major versus minor groove occupancy of monomeric and dimeric crystal violet derivatives. Toward structural correlations

Six monomeric (1a-1f) and five dimeric (2a-2e) derivatives of the triphenylmethane dye crystal violet (CV) have been prepared. Evaluation of the binding of these compounds to CT DNA by competitive fluorescent intercalator displacement (FID) assays, viscosity experiments, and UV and CD spectroscopy suggest that monomeric derivative 1a and dimeric derivative 2d likely associate with the major groove of DNA, while dimeric derivatives 2a and 2e likely associate with the minor groove of DNA. Additional evidence for the groove occupancy assignments of these derivatives was obtained from ITC experiments and from differential inhibition of DNA cleavage by the major groove binding restriction enzyme BamHI, as revealed by agarose gel electrophoresis. The data indicate that major groove ligands may be optimally constructed from dye units containing a sterically bulky 3,5-dimethyl-N,N-dimethylaniline group; furthermore, the groove-selectivity of olefin-tethered dimer 2d suggests that stereoelectronic interactions (n → π*) between the ligand and DNA are also an important design consideration in the crafting of major-groove binding ligands.

Unusual similarity of DNA solvation dynamics in high-salinity crowding with divalent cations of varying concentrations

Double-stranded DNA bears the highest linear negative charge density (2e- per base-pair) among all biopolymers, leading to strong interactions with cations and dipolar water, resulting in the formation of a dense 'condensation layer' around DNA. Interactions involving proteins and ligands binding to DNA are primarily governed by strong electrostatic forces. Increased salt concentrations impede such electrostatic interactions - a situation that prevails in oceanic species due to their cytoplasm being enriched with salts. Nevertheless, how these interactions' dynamics are affected in crowded hypersaline environments remains largely unexplored. Here, we employ steady-state and time-resolved fluorescence Stokes shifts (TRFSS) of a DNA-bound ligand (DAPI) to investigate the static and dynamic solvation properties of DNA in the presence of two divalent cations, magnesium (Mg2+), and calcium (Ca2+) at varying high to very-high concentrations of 0.15 M, 1 M and 2 M. We compare the results to those obtained in physiological concentrations (0.15 M) of monovalent Na+ ions. Combining data from fluorescence femtosecond optical gating (FOG) and time-correlated single photon counting (TCSPC) techniques, dynamic fluorescence Stokes shifts in DNA are analysed over a broad range of time-scales, from 100 fs to 10 ns. We find that while divalent cation crowding strongly influences the DNA stability and ligand binding affinity to DNA, the dynamics of DNA solvation remain remarkably similar across a broad range of five decades in time, even in a high-salinity crowded environment with divalent cations, as compared to the physiological concentration of the Na+ ion. Steady-state and time-resolved data of the DNA-groove-bound ligand are seemingly unaffected by ion-crowding in hypersaline solution, possibly due to ions being mostly displaced by the DNA-bound ligand. Furthermore, the dynamic coupling of cations with nearby water may possibly contribute to a net-neutral effect on the overall collective solvation dynamics in DNA, owing to the strong anti-correlation of their electrostatic interaction energy fluctuations. Such dynamic scenarios may persist within the cellular environment of marine life and other biological cells that experience hypersaline conditions.

Conformational Morphing by a DNA Analogue Featuring 7-Deazapurines and 5-Halogenpyrimidines and the Origins of Adenine-Tract Geometry

Several efforts are currently directed at the creation and cellular implementation of alternative genetic systems composed of pairing components that are orthogonal to the natural dA/dT and dG/dC base pairs. In an alternative approach, Watson-Crick-type pairing is conserved, but one or all of the four letters of the A, C, G, and T alphabet are substituted by modified components. Thus, all four nucleobases were altered to create halogenated deazanucleic acid (DZA): dA was replaced by 7-deaza-2'-deoxyadenosine (dzA), dG by 7-deaza-2'-deoxyguanosine (dzG), dC by 5-fluoro-2'-deoxycytidine (FdC), and dT by 5-chloro-2'-deoxyuridine (CldU). This base-pairing system was previously shown to retain function in Escherichia coli. Here, we analyze the stability, hydration, structure, and dynamics of a DZA Dickerson-Drew Dodecamer (DDD) of sequence 5'-FdC-dzG-FdC-dzG-dzA-dzA-CldU-CldU-FdC-dzG-FdC-dzG-3'. Contrary to similar stabilities of DDD and DZA-DDD, osmotic stressing revealed a dramatic loss of hydration for the DZA-DDD relative to that for the DDD. The parent DDD 5'-d(CGCGAATTCGCG)-3' features an A-tract, a run of adenosines uninterrupted by a TpA step, and exhibits a hallmark narrow minor groove. Crystal structures─in the presence of RNase H─and MD simulations show increased conformational plasticity ("morphing") of DZA-DDD relative to that of the DDD. The narrow dzA-tract minor groove in one structure widens to resemble that in canonical B-DNA in a second structure. These changes reflect an indirect consequence of altered DZA major groove electrostatics (less negatively polarized compared to that in DNA) and hydration (reduced compared to that in DNA). Therefore, chemical modifications outside the minor groove that lead to collapse of major groove electrostatics and hydration can modulate A-tract geometry.

Monovalent cation localization in DNA A-tracts with different sequences

The free solution mobilities of 26-base pair (bp) DNA oligomers containing A-tracts with and without internal ApT steps have been measured by capillary electrophoresis, using the mobility of a 26-bp random-sequence oligomer as a reference. The background electrolytes (BGEs) contained mixtures of Li+ and tetrapropylammonium (TPA+ ) ions, keeping the total cation concentration constant at 0.3 M. The mobility ratios equaled 1.00 in 0.3 M TPA+ , indicating that the A-tract and reference oligomers had the same B-form conformation in this BGE. With increasing [Li+ ], the mobility ratio decreased as Li+ ions became localized in the A-tract minor groove, suggesting that the A-tract was now in the B* conformation. If the A-tract contained an internal ApT step and the oligomer contained less than ∼50% A + T, the mobility ratio reached a reduced plateau value that remained constant as the [Li+ ] increased to 0.3 M. However, for A-tracts without an internal ApT step and for A-tracts embedded in oligomers containing more than 50% A + T, the mobility ratios increased again at high [Li+ ], eventually reaching a plateau value of 1.00. Hence, DNA A-tracts in solution appear to exist as mixtures of the B and B* conformations, with the fractional concentration of each conformer depending on the [Li+ ], the A-tract sequence, and the total A + T content of the oligomer.

Cracking the Floral Quartet Code: How Do Multimers of MIKCC-Type MADS-Domain Transcription Factors Recognize Their Target Genes?

MADS-domain transcription factors (MTFs) are involved in the control of many important processes in eukaryotes. They are defined by the presence of a unique and highly conserved DNA-binding domain, the MADS domain. MTFs bind to double-stranded DNA as dimers and recognize specific sequences termed CArG boxes (such as 5'-CC(A/T)₆GG-3') and similar sequences that occur hundreds of thousands of times in a typical flowering plant genome. The number of MTF-encoding genes increased by around two orders of magnitude during land plant evolution, resulting in roughly 100 genes in flowering plant genomes. This raises the question as to how dozens of different but highly similar MTFs accurately recognize the cis-regulatory elements of diverse target genes when the core binding sequence (CArG box) occurs at such a high frequency. Besides the usual processes, such as the base and shape readout of individual DNA sequences by dimers of MTFs, an important sublineage of MTFs in plants, termed MIKC^C-type MTFs (M^C-MTFs), has evolved an additional mechanism to increase the accurate recognition of target genes: the formation of heterotetramers of closely related proteins that bind to two CArG boxes on the same DNA strand involving DNA looping. M^C-MTFs control important developmental processes in flowering plants, ranging from root and shoot to flower, fruit and seed development. The way in which M^C-MTFs bind to DNA and select their target genes is hence not only of high biological interest, but also of great agronomic and economic importance. In this article, we review the interplay of the different mechanisms of target gene recognition, from the ordinary (base readout) via the extravagant (shape readout) to the idiosyncratic (recognition of the distance and orientation of two CArG boxes by heterotetramers of M^C-MTFs). A special focus of our review is on the structural prerequisites of M^C-MTFs that enable the specific recognition of target genes.

Insights into DNA solvation found in protein-DNA structures

The proteins that bind double-helical DNA present various microenvironments that sense and/or induce signals in the genetic material. The high-resolution structures of protein-DNA complexes reveal the nature of both the microenvironments and the conformational responses in DNA and protein. Complex networks of interactions within the structures somehow tie the protein and DNA together and induce the observed spatial forms. Here we show how the cumulative buildup of amino acid atoms around the sugars, phosphates, and bases in different protein-DNA complexes produces a binding cloud around the double helix and how different types of atoms fill that cloud. Rather than focusing on the principles of molecular binding and recognition suggested by the arrangements of amino acids and nucleotides in the macromolecular complexes, we consider the proteins in contact with DNA as organized solvents. We describe differences in the mix of atoms that come in closest contact with DNA, subtle sequence-dependent features in the microenvironment of the sugar-phosphate backbone, a direct link between the localized buildup of ionic species and the electrostatic potential surfaces of the DNA bases, and sites of atomic buildup above and below the basepair planes that transmit the unique features of the base environments along the chain backbone. The inferences about solvation that can be drawn from the survey provide new stimuli for improvement of nucleic acid force fields and fresh ideas for exploration of the properties of DNA in solution.

Analysis of nucleoid-associated protein-binding regions reveals DNA structural features influencing genome organization in Mycobacterium tuberculosis

Nucleoid-associated proteins (NAPs) maintain bacterial nucleoid configuration through their architectural properties of DNA bending, wrapping, and bridging. However, the contribution of DNA structural alterations to DNA-NAP recognition at the genomic scale remains unresolved. Present work dissects the DNA sequence, shape and altered structural preferences at a genomic scale for six NAPs in Mycobacterium tuberculosis. Results suggest narrower minor groove width (MGW) and higher DNA rigidity are marked for the binding sites of EspR and Lsr2, while mIHF, MtHU and NapM have heterogeneous DNA structural predilections. In contrast, WhiB4-DNA-binding sites were characterized by wider MGW, highly deformable and less curved DNA. This work provides systematic insight into NAP-mediated genome organization as a function of DNA structural features.

Examining the Effects of Netropsin on the Curvature of DNA A-Tracts Using Electrophoresis

A-tracts are sequences of repeated adenine bases that, under the proper conditions, are capable of mediating DNA curvature. A-tracts occur naturally in the regulatory regions of many organisms, yet their biological functions are not fully understood. Orienting multiple A-tracts together constructively or destructively in a phase has the potential to create different shapes in the DNA helix axis. One means of detecting these molecular shape differences is from altered DNA mobilities measured using electrophoresis. The small molecule netropsin binds the minor groove of DNA, particularly at AT-rich sequences including A-tracts. Here, we systematically test the hypothesis that netropsin binding eliminates the curvature of A-tracts by measuring the electrophoretic mobilities of seven 98-base pair DNA samples containing different numbers and arrangements of centrally located A-tracts under varying conditions with netropsin. We find that netropsin binding eliminates the mobility difference between the DNA fragments with different A-tract arrangements in a concentration-dependent manner. This work provides evidence for the straightening of A-tracts upon netropsin binding and illustrates an artificial approach to re-sculpt DNA shape.

Ionic Atmosphere Effect on the Absorption Spectrum of a Flavoprotein: A Reminder to Consider Solution Ions

This study utilizes the FMN-dependent NADH:quinone oxidoreductase from Pseudomonas aeruginosa PAO1 to investigate the effect of introducing an active site negative charge on the flavin absorption spectrum both in the absence and presence of a long-range electrostatic potential coming from solution ions. There were no observed changes in the flavin UV-visible spectrum when an active site tyrosine (Y277) becomes deprotonated in vitro. These results could only be reproduced computationally using average solvent electrostatic configuration (ASEC) QM/MM simulations that include both positive and negative solution ions. The same calculations performed with minimal ions to neutralize the total protein charge predicted that deprotonating Y277 would significantly alter the flavin absorption spectrum. Analyzing the distribution of solution ions indicated that the ions reorganize around the protein surface upon Y277 deprotonation to cancel the effect of the tyrosinate on the flavin absorption spectrum. Additional biochemical experiments were performed to test this hypothesis.

Allostery through DNA drives phenotype switching

Allostery is a pervasive principle to regulate protein function. Growing evidence suggests that also DNA is capable of transmitting allosteric signals. Yet, whether and how DNA-mediated allostery plays a regulatory role in gene expression remained unclear. Here, we show that DNA indeed transmits allosteric signals over long distances to boost the binding cooperativity of transcription factors. Phenotype switching in Bacillus subtilis requires an all-or-none promoter binding of multiple ComK proteins. We use single-molecule FRET to demonstrate that ComK-binding at one promoter site increases affinity at a distant site. Cryo-EM structures of the complex between ComK and its promoter demonstrate that this coupling is due to mechanical forces that alter DNA curvature. Modifications of the spacer between sites tune cooperativity and show how to control allostery, which allows a fine-tuning of the dynamic properties of genetic circuits.

Most insights on DNA-mediated allostery upon transcription factor (TF) binding were either based on artificial promoters or found to be short-ranged. Here authors use single-molecule FRET and cryo-EM to show that Bacillus subtilis bacteria utilize long-range allostery in a stochastic and reversible phenotype switch.

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA's diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

DNA-binding properties of the MADS-domain transcription factor SEPALLATA3 and mutant variants characterized by SELEX-seq

We studied the DNA-binding profile of the MADS-domain transcription factor SEPALLATA3 and mutant variants by SELEX-seq. DNA-binding characteristics of SEPALLATA3 mutant proteins lead us to propose a novel DNA-binding mode. MIKC-type MADS-domain proteins, which function as essential transcription factors in plant development, bind as dimers to a 10-base-pair AT-rich motif termed CArG-box. However, this consensus motif cannot fully explain how the abundant family members in flowering plants can bind different target genes in specific ways. The aim of this study was to better understand the DNA-binding specificity of MADS-domain transcription factors. Also, we wanted to understand the role of a highly conserved arginine residue for binding specificity of the MADS-domain transcription factor family. Here, we studied the DNA-binding profile of the floral homeotic MADS-domain protein SEPALLATA3 by performing SELEX followed by high-throughput sequencing (SELEX-seq). We found a diverse set of bound sequences and could estimate the in vitro binding affinities of SEPALLATA3 to a huge number of different sequences. We found evidence for the preference of AT-rich motifs as flanking sequences. Whereas different CArG-boxes can act as SEPALLATA3 binding sites, our findings suggest that the preferred flanking motifs are almost always the same and thus mostly independent of the identity of the central CArG-box motif. Analysis of SEPALLATA3 proteins with a single amino acid substitution at position 3 of the DNA-binding MADS-domain further revealed that the conserved arginine residue, which has been shown to be involved in a shape readout mechanism, is especially important for the recognition of nucleotides at positions 3 and 8 of the CArG-box motif. This leads us to propose a novel DNA-binding mode for SEPALLATA3, which is different from that of other MADS-domain proteins known.

An artificial cationic oligosaccharide combined with phosphorothioate linkages strongly improves siRNA stability

Small interfering RNAs (siRNAs) are potential tools for gene-silencing therapy, but their instability is one of the obstacles in the development of siRNA-based drugs. To improve siRNA stability, we synthesised a double-stranded RNA-binding cationic oligodiaminogalactose 4mer (ODAGal4) and investigated here its characteristics for siRNA stabilisation in vitro. ODAGal4 improved the resistance of various siRNAs against serum degradation. The effect of ODAGal4 on siRNA stabilisation was further amplified by introduction of modified nucleotides into the siRNA. In particular, a combination of ODAGal4 and incorporation of phosphorothioate linkages into the siRNA prominently prevented degradation by serum. The half-lives of fully phosphorothioate-modified RNA duplexes with ODAGal4 were more than 15 times longer than those of unmodified siRNAs without ODAGal4; this improvement in serum stability was superior to that observed for other chemical modifications. Serum degradation assays of RNAs with multiple chemical modifications showed that ODAGal4 preferentially improves the stability of RNAs with phosphorothioate modification among chemical modifications. Furthermore, melting temperature analysis showed that ODAGal4 greatly increases the thermal stability of phosphorothioate RNAs. Importantly, ODAGal4 did not interrupt gene-silencing activity of all the RNAs tested. Collectively, these findings demonstrate that ODAGal4 is a potent stabiliser of siRNAs, particularly nucleotides with phosphorothioate linkages, representing a promising tool in the development of gene-silencing therapies.

Understanding the paradoxical mechanical response of in-phase A-tracts at different force regimes

A-tracts are A:T rich DNA sequences that exhibit unique structural and mechanical properties associated with several functions in vivo. The crystallographic structure of A-tracts has been well characterized. However, the mechanical properties of these sequences is controversial and their response to force remains unexplored. Here, we rationalize the mechanical properties of in-phase A-tracts present in the Caenorhabditis elegans genome over a wide range of external forces, using single-molecule experiments and theoretical polymer models. Atomic Force Microscopy imaging shows that A-tracts induce long-range (∼200 nm) bending, which originates from an intrinsically bent structure rather than from larger bending flexibility. These data are well described with a theoretical model based on the worm-like chain model that includes intrinsic bending. Magnetic tweezers experiments show that the mechanical response of A-tracts and arbitrary DNA sequences have a similar dependence with monovalent salt supporting that the observed A-tract bend is intrinsic to the sequence. Optical tweezers experiments reveal a high stretch modulus of the A-tract sequences in the enthalpic regime. Our work rationalizes the complex multiscale flexibility of A-tracts, providing a physical basis for the versatile character of these sequences inside the cell.

Employing LiCl salt gradient in the wild-type α-hemolysin nanopore to slow down DNA translocation and detect methylated cytosine

In this research, we demonstrate a label-free detection, biological nanopore-based method to distinguish methylated cytosine (mC) from naked cytosine (C) in sample mixtures containing both C and mC at a prolonged translocation duration. Using a 15-fold increase in LiCl salt concentration going from a cis to trans chamber, we increased the translocation dwell time of ssDNA by over 5-fold and the event capture rate by 6-fold in comparison with the symmetric concentration of 1.0 M KCl (control). This is a consequence of counter-ion binding and effective lowering of the overall charge of DNA, which in turn lessens the electrophoretic drive of the system and slows the translocation velocity. Moreover, salt gradients can create a large electric field that will funnel ions and polymers towards the pore, increasing the capture rate and translocation dwell time of DNA. As a result, in 0.2 M-3.0 M LiCl solution, ssDNA achieved a prolonged dwell time of 52 μs per nucleotide and a capture rate of 60 ssDNA per second. Importantly, lowering the translocation speed of ssDNA enhances the resulting resolution, allowing 5'-mC to be distinguished from C without using methyl-specific labels. We successfully distinguished 5'-mC from C when mixed together at ratios of 1 : 1, 3 : 7 and 7 : 3. The distribution of current blockade amplitudes of all mixtures adopted bimodal shapes, with peak-to-peak ratios coarsely corresponding to the mixture composition (e.g. the density and distribution of events shifted in correspondence with changes in 18b-0mC and 18-2mC concentration ratios in the mixture).

Synthesis and DNA binding profile of monomeric, dimeric, and trimeric derivatives of crystal violet

Monomeric, dimeric, and trimeric derivatives of the triphenylmethane dye crystal violet (1a-1f) have been synthesized for the purpose of evaluating their affinity and sequence selectivity for duplex DNA. Competitive ethidum displacement assays indicate that 1a-1f have apparent association constants for CT DNA in the range of 1.80-16.2 × 107 M-1 and binding site sizes of 10-14 bp. Viscosity experiments performed on ligand 1f confirmed that these dyes associate with duplex DNA by a non-intercalative mode of binding. Circular dichroism and competition binding studies of the tightest binding ligand 1e with known major and minor groove binding molecules suggest that these dye derivatives likely occupy the major groove of DNA. Data from the binding of 1e to polynucleotides indicate close to an order of magnitude preference for associating with AT rich homopolymers over GC rich homopolymers, suggesting a shape-selective match of the sterically bulky ligand with DNA containing a wider major groove.

Characterization of the binding interactions between EvaGreen dye and dsDNA

Understanding the dsDNA·EG binding interaction is important because the EvaGreen (EG) dye is increasingly used in real-time quantitative polymerase chain reaction, high resolution melting analysis, and routine quantification of DNA. In this work, a binding isotherm for the interactions of EG with duplex DNA (poly-dA₁₇·poly-dT₁₇) has been determined from the absorption and fluorescence spectra of the EG and dsDNA·EG complex. The isotherm has a sigmoidal shape and can be modeled with the Hill equation, indicating positive cooperativity for the binding interaction. A Scatchard plot of the binding data yields a concave-down curve in agreement with the Hill analysis of the binding isotherm for a positive cooperative binding interaction. Analysis of the Scatchard plot with the modified McGhee and von Hippel model for a finite one-dimensional homogeneous lattice and nonspecific binding of ligands to duplex DNA yields the intrinsic binding constant, the number of lattice sites occluded by a bound ligand, and the cooperativity parameter of 3.6 × 10⁵ M^-1, 4.0, and 8.1, respectively. The occluded site size of 4 indicates that moieties of the EG intercalate into the adjacent base pairs of the duplex DNA with a gap of 1 intercalation site between EG binding sites, as expected for a bifunctional molecule. Interestingly, at high [EG]/[base pair], the intercalation is disrupted. A model is proposed based on the fluorescence spectrum where the formation of anti-parallel stacked chains of EGs bound externally to the duplex DNA occur at these high ratios.

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.

The floral homeotic protein SEPALLATA3 recognizes target DNA sequences by shape readout involving a conserved arginine residue in the MADS-domain

SEPALLATA3 of Arabidopsis thaliana is a MADS-domain transcription factor (TF) and a key regulator of flower development. MADS-domain proteins bind to sequences termed 'CArG-boxes' [consensus 5'-CC(A/T)₆ GG-3']. Because only a fraction of the CArG-boxes in the Arabidopsis genome are bound by SEPALLATA3, more elaborate principles have to be discovered to better understand which features turn CArG-boxes into genuine recognition sites. Here, we investigate to what extent the shape of the DNA is involved in a 'shape readout' that contributes to the binding of SEPALLATA3. We determined in vitro binding affinities of SEPALLATA3 to DNA probes that all contain the CArG-box motif, but differ in their predicted DNA shape. We found that binding affinity correlates well with a narrow minor groove of the DNA. Substitution of canonical bases with non-standard bases supports the hypothesis of minor groove shape readout by SEPALLATA3. Analysis of mutant SEPALLATA3 proteins further revealed that a highly conserved arginine residue, which is expected to contact the DNA minor groove, contributes significantly to the shape readout. Our studies show that the specific recognition of cis-regulatory elements by a plant MADS-domain TF, and by inference probably also of other TFs of this type, heavily depends on shape readout mechanisms.

An overview of recent advances in duplex DNA recognition by small molecules

As the carrier of genetic information, the DNA double helix interacts with many natural ligands during the cell cycle, and is amenable to such intervention in diseases such as cancer biogenesis. Proteins bind DNA in a site-specific manner, not only distinguishing between the geometry of the major and minor grooves, but also by making close contacts with individual bases within the local helix architecture. Over the last four decades, much research has been reported on the development of small non-natural ligands as therapeutics to either block, or in some cases, mimic a DNA–protein interaction of interest. This review presents the latest findings in the pursuit of novel synthetic DNA binders. This article provides recent coverage of major strategies (such as groove recognition, intercalation and cross-linking) adopted in the duplex DNA recognition by small molecules, with an emphasis on major works of the past few years.

Co(II) and Ni(II) binding of the Escherichia coli transcriptional repressor RcnR orders its N terminus, alters helix dynamics, and reduces DNA affinity

RcnR, a transcriptional regulator in Escherichia coli, derepresses the expression of the export proteins RcnAB upon binding Ni(II) or Co(II). Lack of structural information has precluded elucidation of the allosteric basis for the decreased DNA affinity in RcnR's metal-bound states. Here, using hydrogen-deuterium exchange coupled with MS (HDX-MS), we probed the RcnR structure in the presence of DNA, the cognate metal ions Ni(II) and Co(II), or the noncognate metal ion Zn(II). We found that cognate metal binding altered flexibility from the N terminus through helix 1 and modulated the RcnR-DNA interaction. Apo-RcnR and RcnR-DNA complexes and the Zn(II)-RcnR complex exhibited similar ²H uptake kinetics, with fast-exchanging segments located in the N terminus, in helix 1 (residues 14-24), and at the C terminus. The largest difference in ²H incorporation between apo- and Ni(II)- and Co(II)-bound RcnR was observed in helix 1, which contains the N terminus and His-3, and has been associated with cognate metal binding. ²H uptake in helix 1 was suppressed in the Ni(II)- and Co(II)-bound RcnR complexes, in particular in the peptide corresponding to residues 14-24, containing Arg-14 and Lys-17. Substitution of these two residues drastically affected DNA-binding affinity, resulting in rcnA expression in the absence of metal. Our results suggest that cognate metal binding to RcnR orders its N terminus, decreases helix 1 flexibility, and induces conformational changes that restrict DNA interactions with the positively charged residues Arg-14 and Lys-17. These metal-induced alterations decrease RcnR-DNA binding affinity, leading to rcnAB expression.

Functional Mechanisms of Microsatellite DNA in Eukaryotic Genomes

Microsatellite repeat DNA is best known for its length mutability, which is implicated in several neurological diseases and cancers, and often exploited as a genetic marker. Less well-known is the body of work exploring the widespread and surprisingly diverse functional roles of microsatellites. Recently, emerging evidence includes the finding that normal microsatellite polymorphism contributes substantially to the heritability of human gene expression on a genome-wide scale, calling attention to the task of elucidating the mechanisms involved. At present, these are underexplored, but several themes have emerged. I review evidence demonstrating roles for microsatellites in modulation of transcription factor binding, spacing between promoter elements, enhancers, cytosine methylation, alternative splicing, mRNA stability, selection of transcription start and termination sites, unusual structural conformations, nucleosome positioning and modification, higher order chromatin structure, noncoding RNA, and meiotic recombination hot spots.

Large-scale chromosome folding versus genomic DNA sequences: A discrete double Fourier transform technique

Using state-of-the-art techniques combining imaging methods and high-throughput genomic mapping tools leaded to the significant progress in detailing chromosome architecture of various organisms. However, a gap still remains between the rapidly growing structural data on the chromosome folding and the large-scale genome organization. Could a part of information on the chromosome folding be obtained directly from underlying genomic DNA sequences abundantly stored in the databanks? To answer this question, we developed an original discrete double Fourier transform (DDFT). DDFT serves for the detection of large-scale genome regularities associated with domains/units at the different levels of hierarchical chromosome folding. The method is versatile and can be applied to both genomic DNA sequences and corresponding physico-chemical parameters such as base-pairing free energy. The latter characteristic is closely related to the replication and transcription and can also be used for the assessment of temperature or supercoiling effects on the chromosome folding. We tested the method on the genome of E. coli K-12 and found good correspondence with the annotated domains/units established experimentally. As a brief illustration of further abilities of DDFT, the study of large-scale genome organization for bacteriophage PHIX174 and bacterium Caulobacter crescentus was also added. The combined experimental, modeling, and bioinformatic DDFT analysis should yield more complete knowledge on the chromosome architecture and genome organization.

Contributions of Sequence to the Higher-Order Structures of DNA

One of the critical unanswered questions in genome biophysics is how the primary sequence of DNA bases influences the global properties of very-long-chain molecules. The local sequence-dependent features of DNA found in high-resolution structures introduce irregularities in the disposition of adjacent residues that facilitate the specific binding of proteins and modulate the global folding and interactions of double helices with hundreds of basepairs. These features also determine the positions of nucleosomes on DNA and the lengths of the interspersed DNA linkers. Like the patterns of basepair association within DNA, the arrangements of nucleosomes in chromatin modulate the properties of longer polymers. The intrachromosomal loops detected in genomic studies contain hundreds of nucleosomes, and given that the simulated configurations of chromatin depend on the lengths of linker DNA, the formation of these loops may reflect sequence-dependent information encoded within the positioning of the nucleosomes. With knowledge of the positions of nucleosomes on a given genome, methods are now at hand to estimate the looping propensities of chromatin in terms of the spacing of nucleosomes and to make a direct connection between the DNA base sequence and larger-scale chromatin folding.

Photophysical Characterization of Enhanced 6-Methylisoxanthopterin Fluorescence in Duplex DNA

The structure and dynamic motions of bases in DNA duplexes and other constructs are important for understanding mechanisms of selectivity and recognition of DNA-binding proteins. The fluorescent guanine analogue, 6-methylisoxanthopterin 6-MI, is well suited to this purpose as it exhibits an unexpected 3- to 4-fold increase in relative quantum yield upon duplex formation when incorporated into the following sequences: ATFAA, AAFTA, or ATFTA (where F represents 6-MI). To better understand some of the factors leading to the 6-MI fluorescence increase upon duplex formation, we characterized the effect of local sequence and structural perturbations on 6-MI photophysics through temperature melts, quantum yield measurements, fluorescence quenching assays, and fluorescence lifetime measurements. By examining 21 sequences we have determined that the duplex-enhanced fluorescence (DEF) depends on the composition of bases adjacent to 6-MI and the presence of adenines at locations n ± 2 from the probe. Investigation of duplex stability and local solvent accessibility measurements support a model in which the DEF arises from a constrained geometry of 6-MI in the duplex, which remains H-bonded to cytosine, stacked with adjacent bases and inaccessible to quenchers. Perturbation of DNA structure through the introduction of an unpaired base 3' to 6-MI or a mismatched basepair increases 6-MI dynamic motion leading to fluorescence quenching and a reduction in quantum yield. Molecular dynamics simulations suggest the enhanced fluorescence results from a greater degree of twist at the X-F step relative to the quenched duplexes examined. These results point to a model where adenine residues located at n ± 2 from 6-MI induce a structural geometry with greater twist in the duplex that hinders local motion reducing dynamic quenching and producing an increase in 6-MI fluorescence.

Shape readout of AT-rich DNA by carbohydrates

Gene expression can be altered by small molecules that target DNA; sequence as well as shape selectivities are both extremely important for DNA recognition by intercalating and groove-binding ligands. We have characterized a carbohydrate scaffold (1) exhibiting DNA "shape readout" properties. Thermodynamic studies with 1 and model duplex DNAs demonstrate the molecule's high affinity and selectivity towards B* form (continuous AT-rich) DNA. Isothermal titration calorimetry (ITC), circular dichroism (CD) titration, ultraviolet (UV) thermal denaturation, and Differential Scanning Calorimetry were used to characterize the binding of 1 with a B* form AT-rich DNA duplex d[5'-G2 A6 T6 C2 -3']. The binding constant was determined using ITC at various temperatures, salt concentrations, and pH. ITC titrations were fit using a two-binding site model. The first binding event was shown to have a 1:1 binding stoichiometry and was predominantly entropy-driven with a binding constant of approximately 10(8) M(-1) . ITC-derived binding enthalpies were used to obtain the binding-induced change in heat capacity (ΔCp ) of -225 ± 19 cal/mol·K. The ionic strength dependence of the binding constant indicated a significant electrolytic contribution in ligand:DNA binding, with approximately four to five ion pairs involved in binding. Ligand 1 displayed a significantly higher affinity towards AT-tract DNA over sequences containing GC inserts, and binding experiments revealed the order of binding affinity for 1 with DNA duplexes: contiguous B* form AT-rich DNA (d[5'-G2 A6 T6 C2 -3']) >B form alternate AT-rich DNA (d[5'-G2 (AT)6 C2- 3']) > A form GC-rich DNA (d[5'-A2 G6 C6 T2 -3']), demonstrating the preference of ligand 1 for B* form DNA.

Flanking A·T basepairs destabilize the B(∗) conformation of DNA A-tracts

Capillary electrophoresis has been used to characterize the interaction of monovalent cations with 26-basepair DNA oligomers containing A-tracts embedded in flanking sequences with different basepair compositions. A 26-basepair random-sequence oligomer was used as the reference; lithium and tetrabutylammonium (TBA(+)) ions were used as the probe ions. The free solution mobilities of the A-tract and random-sequence oligomers were identical in solutions containing <∼ 100 mM cation. At higher cation concentrations, the A-tract oligomers migrated faster than the reference oligomer in TBA(+) and slower than the reference in Li(+). Hence, cations of different sizes can interact very differently with DNA A-tracts. The increased mobilities observed in TBA(+) suggest that the large hydrophobic TBA(+) ions are preferentially excluded from the vicinity of the A-tract minor groove, increasing the effective net charge of the A-tract oligomers and increasing the mobility. By contrast, Li(+) ions decrease the mobility of A-tract oligomers because of the preferential localization of Li(+) ions in the narrow A-tract minor groove. Embedding the A-tracts in AT-rich flanking sequences markedly alters preferential interactions of monovalent cations with the B(∗) conformation. Hence, A-tracts embedded in genomic DNA may or may not interact preferentially with monovalent cations, depending on the relative number of A · T basepairs in the flanking sequences.

Organization of DNA in a bacterial nucleoid

BACKGROUND

It is unclear how DNA is packaged in a bacterial cell in the absence of nucleosomes. To investigate the initial level of DNA condensation in bacterial nucleoid we used in vivo DNA digestion coupled with high-throughput sequencing of the digestion-resistant fragments. To this end, we transformed E. coli cells with a plasmid expressing micrococcal nuclease. The nuclease expression was under the control of AraC repressor, which enabled us to perform an inducible digestion of bacterial nucleoid inside a living cell.

RESULTS

Analysis of the genomic localization of the digestion-resistant fragments revealed their non-random distribution. The patterns observed in the distribution of the sequenced fragments indicate the presence of short DNA segments protected from the enzyme digestion, possibly because of interaction with DNA-binding proteins. The average length of such digestion-resistant segments is about 50 bp and the characteristic repeat in their distribution is about 90 bp. The gene starts are depleted of the digestion-resistant fragments, suggesting that these genomic regions are more exposed than genomic sequences on average. Sequence analysis of the digestion-resistant segments showed that while the GC-content of such sequences is close to the genome-wide value, they are depleted of A-tracts as compared to the bulk genomic DNA or to the randomized sequence of the same nucleotide composition.

CONCLUSIONS

Our results suggest that DNA is packaged in the bacterial nucleoid in a non-random way that facilitates interaction of the DNA binding factors with regulatory regions of the genome.

A new general model for predicting melting thermodynamics of complementary and mismatched B-form duplexes containing locked nucleic acids: application to probe design for digital PCR detection of somatic mutations

Advances in real-time polymerase chain reaction (PCR), as well as the emergence of digital PCR (dPCR) and useful modified nucleotide chemistries, including locked nucleic acids (LNAs), have created the potential to improve and expand clinical applications of PCR through their ability to better quantify and differentiate amplification products, but fully realizing this potential will require robust methods for designing dual-labeled hydrolysis probes and predicting their hybridization thermodynamics as a function of their sequence, chemistry, and template complementarity. We present here a nearest-neighbor thermodynamic model that accurately predicts the melting thermodynamics of a short oligonucleotide duplexed either to its perfect complement or to a template containing mismatched base pairs. The model may be applied to pure-DNA duplexes or to duplexes for which one strand contains any number and pattern of LNA substitutions. Perturbations to duplex stability arising from mismatched DNA:DNA or LNA:DNA base pairs are treated at the Gibbs energy level to maintain statistical significance in the regressed model parameters. This approach, when combined with the model's accounting of the temperature dependencies of the melting enthalpy and entropy, permits accurate prediction of T(m) values for pure-DNA homoduplexes or LNA-substituted heteroduplexes containing one or two independent mismatched base pairs. Terms accounting for changes in solution conditions and terminal addition of fluorescent dyes and quenchers are then introduced so that the model may be used to accurately predict and thereby tailor the T(m) of a pure-DNA or LNA-substituted hydrolysis probe when duplexed either to its perfect-match template or to a template harboring a noncomplementary base. The model, which builds on classic nearest-neighbor thermodynamics, should therefore be of use to clinicians and biologists who require probes that distinguish and quantify two closely related alleles in either a quantitative PCR or dPCR assay. This potential is demonstrated by using the model to design allele-specific probes that completely discriminate and quantify clinically relevant mutant alleles (BRAF V600E and KIT D816V) in a dPCR assay.

Recent findings on the mechanisms involved in tenofovir resistance

Since its approval for clinical use in 2001, tenofovir (TFV) has become one of the most frequently prescribed nucleotide analogues used in combination with other antiretroviral agents against HIV-1 infection. Although reverse transcriptase inhibitors (RTIs) including TFV have been shown to be highly potent with reasonable safety profiles in the clinic, drug resistance hinders the effectiveness of current therapies and even causes treatment failure. Therefore, understanding the resistance mechanisms of RT and exploring the potential antiviral synergy between the different RTIs in combination therapies against the resistance mechanisms would greatly improve the long-term efficacy of existing and future regimens. We have studied the pyrophosphorolytic removal of TFV, a major resistance mechanism that RT utilizes, from two different viral sequences and observed interesting outcomes associated with the sequence context. Furthermore, addition of efavirenz, a non-nucleoside RTI, inhibits this removal process confirming the synergistic antiviral effects. This article highlights our recently published work on the viral sequence context contributing to the study of anti-HIV drug resistance in conjunction with the benefits of combining various RTIs that may have been neglected previously.

Specific minor groove solvation is a crucial determinant of DNA binding site recognition

The DNA sequence preferences of nearly all sequence specific DNA binding proteins are influenced by the identities of bases that are not directly contacted by protein. Discrimination between non-contacted base sequences is commonly based on the differential abilities of DNA sequences to allow narrowing of the DNA minor groove. However, the factors that govern the propensity of minor groove narrowing are not completely understood. Here we show that the differential abilities of various DNA sequences to support formation of a highly ordered and stable minor groove solvation network are a key determinant of non-contacted base recognition by a sequence-specific binding protein. In addition, disrupting the solvent network in the non-contacted region of the binding site alters the protein's ability to recognize contacted base sequences at positions 5–6 bases away. This observation suggests that DNA solvent interactions link contacted and non-contacted base recognition by the protein.

Close encounters with DNA

Over the past ten years, the all-atom molecular dynamics method has grown in the scale of both systems and processes amenable to it and in its ability to make quantitative predictions about the behavior of experimental systems. The field of computational DNA research is no exception, witnessing a dramatic increase in the size of systems simulated with atomic resolution, the duration of individual simulations and the realism of the simulation outcomes. In this topical review, we describe the hallmark physical properties of DNA from the perspective of all-atom simulations. We demonstrate the amazing ability of such simulations to reveal the microscopic physical origins of experimentally observed phenomena. We also discuss the frustrating limitations associated with imperfections of present atomic force fields and inadequate sampling. The review is focused on the following four physical properties of DNA: effective electric charge, response to an external mechanical force, interaction with other DNA molecules and behavior in an external electric field.

Interactions with polynucleotides and antitumor activity of amidino and imidazolinyl substituted 2-phenylbenzothiazole mesylates

Based on previously reported antiproliferative activity screening, four most promising disubstituted 2-phenylbenzothiazole hydrochlorides were chosen for detailed study. Water solubility, as well as liphophilicity/hydrophilicity balance of organic core were modified by conversion to mesylate salts. For purpose of structure/activity studies their structures were determined by X-ray structure analysis. Detailed analysis of interactions of new compounds with double stranded (ds-) DNA/RNA by UV/Vis and CD titrations, thermal melting and viscometry experiments revealed that most of studied compounds intercalate into ds-RNA but bind into minor groove of AT-DNA, and agglomerate along GC-DNA. Furthermore, compounds also interact with ss-RNA, but only amino-imidazolinyl 2-phenylbenzothiazole, 4b displayed well defined orientation and dominant binding mode (by induced CD signals) with poly A and poly G. Besides, in vitro investigations revealed moderate to high antiproliferative activity of benzothiazoles against seven human cancer cell lines, while in some cases (HTC 116, SW620, MIA PaCa-2) high correlation between the type of the amidino group and cytotoxic activity was observed.

Sequence-specific DNA interactions with calixarene-based langmuir monolayers

The interactions of an amphiphilic calixarene, namely p-guanidino-dodecyloxy-calix[4]arene, 1, self-assembled as Langmuir monolayers, with short double stranded DNA, were investigated by surface pressure-area (π-A) isotherms, surface ellipsometry and Brewster angle microscopy (BAM). Three DNA 30mers were used as models, poly(AT), poly(GC) and a random DNA sequence with 50% of G:C base pairs. The interactions of these model DNA duplexes with 1-based Langmuir monolayers were studied by measuring compression isotherms using increasing DNA concentrations (10(-6), 10(-5), 10(-4), and 5 × 10(-4) g L(-1)) in the aqueous subphase. The isotherms of 1 showed an expansion of the monolayer with, interestingly, significant differences depending on the duplex DNA sequence studied. Indeed, the interactions of 1-based monolayers with poly(AT) led to an expansion of the monolayer that was significantly more pronounced that for monolayers on subphases of poly(GC) and the random DNA sequence. The structure and thickness of 1-based Langmuir monolayers were investigated by BAM and surface ellipsometry that showed differences in thickness and structure between a monolayer formed on pure water or on a DNA subphase, with here again relevant dissimilarities depending on the DNA composition.

DNA damage by histone radicals in nucleosome core particles

Although DNA binding proteins shield the genetic material from diffusible reactive oxygen species by reacting with them, the resulting protein (peroxyl) radicals can oxidize the bound DNA. To explore this possible DNA damage by protein radicals, histone H4 proteins containing an azoalkane radical precursor at defined sites were prepared. Photolysis of a nucleosome core particle containing the modified protein produces DNA damage that is consistent with selective C4′-oxidation. The nucleotide(s) damaged is highly dependent on proximity to the protein radical. These experiments provide insight into the effects of oxidative stress on protein-bound DNA, revealing an additional layer of complexity concerning nucleic acid damage.

MD and NMR analyses of choline and TMA binding to duplex DNA: on the origins of aberrant sequence-dependent stability by alkyl cations in aqueous and water-free solvents

It has been known for decades that alkylammonium ions, such as tetramethyl ammonium (TMA), alter the usual correlation between DNA GC-content and duplex stability. In some cases it is even possible for an AT-rich duplex to be more stable than a GC-rich duplex of the same length. There has been much speculation regarding the origin of this aberration in sequence-dependent DNA duplex stability, but no clear resolution. Using a combination of molecular dynamics simulations and NMR spectroscopy we demonstrate that choline (2-hydroxy-N,N,N-trimethylethanaminium) and TMA are preferentially localized in the minor groove of DNA duplexes at A·T base pairs and these same ions show less pronounced localization in the major groove compared to what has been demonstrated for alkali and alkali earth metal ions. Furthermore, free energy calculations show that single-stranded GC-rich sequences exhibit more favorable solvation by choline than single-stranded AT-rich sequences. The sequence-specific nature of choline and TMA binding provides a rationale for the enhanced stability of AT-rich sequences when alkyl-ammonium ions are used as the counterions of DNA. Our combined theoretical and experimental study provides one of the most detailed pictures to date of cations localized along DNA in the solution state, and provides insights that go beyond understanding alkyl-ammonium ion binding to DNA. In particular, because choline and TMA bind to DNA in a manner that is found to be distinct from that previously reported for Na(+), K(+), Mg(2+), and Ca(2+), our results reveal the important but underappreciated role that most other cations play in sequence-specific duplex stability.

DNA A-tracts are not curved in solutions containing high concentrations of monovalent cations

The intrinsic curvature of seven 98 bp DNA molecules containing up to four centrally located A6-tracts has been measured by gel and capillary electrophoresis as a function of the number and arrangement of the A-tracts. At low cation concentrations, the electrophoretic mobility observed in polyacrylamide gels and in free solution decreases progressively with the increasing number of phased A-tracts, as expected for DNA molecules with increasingly curved backbone structures. Anomalously slow electrophoretic mobilities are also observed for DNA molecules containing two pairs of phased A-tracts that are out of phase with each other, suggesting that out-of-phase distortions of the helix backbone do not cancel each other out. The mobility decreases observed for the A-tract samples are due to curvature, not cation binding in the A-tract minor groove, because identical free solution mobilities are observed for a molecule with four out-of-phase A-tracts and one with no A-tracts. Surprisingly, the curvature of DNA A-tracts is gradually lost when the monovalent cation concentration is increased to ∼200 mM, regardless of whether the cation is a hydrophilic ion like Na+, NH4+, or Tris+ or a hydrophobic ion like tetrabutylammonium. The decrease in A-tract curvature with increasing ionic strength, along with the known decrease in A-tract curvature with increasing temperature, suggests that DNA A-tracts are not significantly curved under physiological conditions.

Control of DNA minor groove width and Fis protein binding by the purine 2-amino group

The width of the DNA minor groove varies with sequence and can be a major determinant of DNA shape recognition by proteins. For example, the minor groove within the center of the Fis–DNA complex narrows to about half the mean minor groove width of canonical B-form DNA to fit onto the protein surface. G/C base pairs within this segment, which is not contacted by the Fis protein, reduce binding affinities up to 2000-fold over A/T-rich sequences. We show here through multiple X-ray structures and binding properties of Fis–DNA complexes containing base analogs that the 2-amino group on guanine is the primary molecular determinant controlling minor groove widths. Molecular dynamics simulations of free-DNA targets with canonical and modified bases further demonstrate that sequence-dependent narrowing of minor groove widths is modulated almost entirely by the presence of purine 2-amino groups. We also provide evidence that protein-mediated phosphate neutralization facilitates minor groove compression and is particularly important for binding to non-optimally shaped DNA duplexes.

Unusually high-affinity Mg(2+) binding at the AU-rich sequence within the antiterminator hairpin of a Mg(2+) riboswitch

Mg(2+)-Responsive riboswitches represent a fascinating example of bifunctional RNAs that sense Mg(2+) ions with high selectivity and autonomously regulate the expression of Mg(2+)-transporter proteins. The mechanism of the mgtA riboswitch is scarcely understood, and a detailed structural analysis is called for to study how this RNA can selectively recognize Mg(2+) and respond by switching between two alternative stem loop structures. In this work, we investigated the structure and Mg(2+)-binding properties of the lower part of the antiterminator loop C from the mgtA riboswitch of Yersinia enterocolitica by solution NMR and report a discrete Mg(2+)-binding site embedded in the AU-rich sequence. At the position of Mg(2+) binding, the helical axis exhibits a distinct kink accompanied by a widening of the major groove, which accommodates the Mg(2+)-binding pocket. An unusually large overlap between two adenine residues on the opposite strands suggests that the bending may be sequence-induced by strong stacking interactions, enabling Mg(2+) to bind at this so-far not described metal-ion binding site.

Indirect readout of DNA sequence by p22 repressor: roles of DNA and protein functional groups in modulating DNA conformation

The repressor of bacteriophage P22 (P22R) discriminates between its various DNA binding sites by sensing the identity of non-contacted base pairs at the center of its binding site. The "indirect readout" of these non-contacted bases is apparently based on DNA's sequence-dependent conformational preferences. The structures of P22R-DNA complexes indicate that the non-contacted base pairs at the center of the binding site are in the B' state. This finding suggests that indirect readout and therefore binding site discrimination depend on P22R's ability to either sense and/or impose the B' state on the non-contacted bases of its binding sites. We show here that the affinity of binding sites for P22R depends on the tendency of the central bases to assume the B'-DNA state. Furthermore, we identify functional groups in the minor groove of the non-contacted bases as the essential modulators of indirect readout by P22R. In P22R-DNA complexes, the negatively charged E44 and E48 residues are provocatively positioned near the negatively charged DNA phosphates of the non-contacted nucleotides. The close proximity of the negatively charged groups on protein and DNA suggests that electrostatics may play a key role in the indirect readout process. Changing either of two negatively charged residues to uncharged residues eliminates the ability of P22R to impose structural changes on DNA and to recognize non-contacted base sequence. These findings suggest that these negatively charged amino acids function to force the P22R-bound DNA into the B' state and therefore play a key role in indirect readout by P22R.

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.

Understanding the molecular mechanism of sequence dependent tenofovir removal by HIV-1 reverse transcriptase: differences in primer binding site versus polypurine tract

Tenofovir (TFV) is a nucleotide reverse transcriptase inhibitor (NtRTI) that is often administered as first-line therapy against human immunodeficiency virus type-1 (HIV-1) infection and acts as a chain terminator when incorporated into viral DNA. However, HIV-1 reverse transcriptase (RT) excises TFV in the presence of either ATP or pyrophosphate, which is an important drug resistance mechanism that would interfere with the effective treatment. Previous studies have shown conflicting results on excision efficiencies for TFV-terminated primer-templates derived from either primer binding site (PBS) or polypurine tract (PPT) sequences. To provide mechanistic insight into the variation in TFV removal from both sequences that are vital for the HIV-1 life cycle, we compared the efficiencies of removal reaction in response to sequence dependence via utilizing blocked PBS and PPT primer-templates. We found an enhanced TFV excision with PPT sequence over PBS sequence through ATP-mediated removal and a subsequent incorporation of ATP into the unblocked primers. Furthermore, the rate of pyrophosphorolytic excision of TFV from PPT sequence was 21-fold higher than that for the PBS sequence. However, the addition of efavirenz, nonnucleoside reverse transcriptase inhibitor (NNRTI), to the removal reaction effectively inhibits the TFV excision from both primers by forming a stable complex that would leave TFV inaccessible for excision. These results illuminate the degree of primer-template sequence contribution on TFV removal as well as increase our understanding of the molecular mechanism for the beneficial effects of widely used combinations of antiretroviral regimens in the context of synergistic antiviral activity and drug resistance.

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.

B-DNA structure is intrinsically polymorphic: even at the level of base pair positions

Increasingly exact measurement of single crystal X-ray diffraction data offers detailed characterization of DNA conformation, hydration and electrostatics. However, instead of providing a more clear and unambiguous image of DNA, highly accurate diffraction data reveal polymorphism of the DNA atomic positions and conformation and hydration. Here we describe an accurate X-ray structure of B-DNA, painstakingly fit to a multistate model that contains multiple competing positions of most of the backbone and of entire base pairs. Two of ten base-pairs of CCAGGCCTGG are in multiple states distinguished primarily by differences in slide. Similarly, all the surrounding ions are seen to fractionally occupy discrete competing and overlapping sites. And finally, the vast majority of water molecules show strong evidence of multiple competing sites. Conventional resolution appears to give a false sense of homogeneity in conformation and interactions of DNA. In addition, conventional resolution yields an average structure that is not accurate, in that it is different from any of the multiple discrete structures observed at high resolution. Because base pair positional heterogeneity has not always been incorporated into model-building, even some high and ultrahigh-resolution structures of DNA do not indicate the full extent of conformational polymorphism.

Energies and pressures in viruses: contribution of nonspecific electrostatic interactions

We summarize some aspects of electrostatic interactions in the context of viruses. A simplified but, within well defined limitations, reliable approach is used to derive expressions for electrostatic energies and the corresponding osmotic pressures in single-stranded RNA viruses and double-stranded DNA bacteriophages. The two types of viruses differ crucially in the spatial distribution of their genome charge which leads to essential differences in their free energies, depending on the capsid size and total charge in a quite different fashion. Differences in the free energies are trailed by the corresponding characteristics and variations in the osmotic pressure between the inside of the virus and the external bathing solution.

Neomycin-neomycin dimer: an all-carbohydrate scaffold with high affinity for AT-rich DNA duplexes

A dimeric neomycin-neomycin conjugate 3 with a flexible linker, 2,2'-(ethylenedioxy)bis(ethylamine), has been synthesized and characterized. Dimer 3 can selectively bind to AT-rich DNA duplexes with high affinity. Biophysical studies have been performed between 3 and different nucleic acids with varying base composition and conformation by using ITC (isothermal calorimetry), CD (circular dichroism), FID (fluorescent intercalator displacement), and UV (ultraviolet) thermal denaturation experiments. A few conclusions can be drawn from this study: (1) FID assay with 3 and polynucleotides demonstrates the preference of 3 toward AT-rich sequences over GC-rich sequences. (2) FID assay and UV thermal denaturation experiments show that 3 has a higher affinity for the poly(dA)·poly(dT) DNA duplex than for the poly(dA)·2poly(dT) DNA triplex. Contrary to neomycin, 3 destabilizes poly(dA)·2poly(dT) triplex but stabilizes poly(dA)·poly(dT) duplex, suggesting the major groove as the binding site. (3) UV thermal denaturation studies and ITC experiments show that 3 stabilizes continuous AT-tract DNA better than DNA duplexes with alternating AT bases. (4) CD and FID titration studies show a DNA binding site size of 10-12 base pairs/drug, depending upon the structure/sequence of the duplex for AT-rich DNA duplexes. (5) FID and ITC titration between 3 and an intramolecular DNA duplex [d(5'-A(12)-x-T(12)-3'), x = hexaethylene glycol linker] results in a binding stoichiometry of 1:1 with a binding constant ∼10(8) M(-1) at 100 mM KCl. (6) FID assay using 3 and 512 hairpin DNA sequences that vary in their AT base content and placement also show a higher binding selectivity of 3 toward continuous AT-rich than toward DNA duplexes with alternate AT base pairs. (7) Salt-dependent studies indicate the formation of three ion pairs during binding of the DNA duplex d[5'-A(12)-x-T(12)-3'] and 3. (8) ITC-derived binding constants between 3 and DNA duplexes have the following order: AT continuous, d[5'-G(3)A(5)T(5)C(3)-3'] > AT alternate, d[5'-G(3)(AT)(5)C(3)-3'] > GC-rich d[5'-A(3)G(5)C(5)T(3)-3']. (9) 3 binds to the AT-tract-containing DNA duplex (B* DNA, d[5'-G(3)A(5)T(5)C(3)-3']) with 1 order of magnitude higher affinity than to a DNA duplex with alternating AT base pairs (B DNA, d[5'-G(3)(AT)(5)C(3)-3']) and with almost 3 orders of magnitude higher affinity than a GC-rich DNA (A-form, d[5'-A(3)G(5)C(5)T(3)-3']).