A-form conformational motifs in ligand-bound DNA structures

A mechanistic model of primer synthesis from catalytic structures of DNA polymerase α-primase

The mechanism by which polymerase α-primase (polα-primase) synthesizes chimeric RNA-DNA primers of defined length and composition, necessary for replication fidelity and genome stability, is unknown. Here, we report cryo-EM structures of Xenopus laevis polα-primase in complex with primed templates representing various stages of DNA synthesis. Our data show how interaction of the primase regulatory subunit with the primer 5' end facilitates handoff of the primer to polα and increases polα processivity, thereby regulating both RNA and DNA composition. The structures detail how flexibility within the heterotetramer enables synthesis across two active sites and provide evidence that termination of DNA synthesis is facilitated by reduction of polα and primase affinities for the varied conformations along the chimeric primer-template duplex. Together, these findings elucidate a critical catalytic step in replication initiation and provide a comprehensive model for primer synthesis by polα-primase.

Multistate B- to A- transition in protein-DNA Binding - How well is it described by current AMBER force fields?

When DNA interacts with a protein, its structure often undergoes a significant conformational adaptation, usually involving a transition from B-DNA towards the A-DNA form. This is not a two-state, but rather a multistate transition. The A- and B- forms differ mainly in sugar pucker (north/south) and glycosidic torsion χ (anti/high-anti). The combination of A-like pucker and B-like χ (and vice versa) represents the nature of the intermediate states between the pure A- and B- forms. Here we study how the A/B equilibrium and the A/B intermediate states at protein-DNA interfaces are modeled by current AMBER force fields. Eight diverse protein-DNA complexes and their naked (unbound) DNAs were simulated with OL15 and bsc1 force fields and an experimental combination OL15χOL3. We found that while the geometries of the A-like intermediate states agree well with the native X-ray geometries, their populations (stabilities) are significantly underestimated. Different force fields predict different propensities for A-like states growing in the order OL15 < bsc1 < OL15χOL3, yet all underestimate A-like form populations. Interestingly, the force fields seem to predict the correct sequence-dependent A-form propensity, as they predict larger populations of the A-like form in unbound DNA in those steps that acquire A-like conformations in protein-DNA complexes. The instability of A-like geometries in current force fields significantly alters the geometry of simulated protein-DNA complexes and destabilizes the binding motif, suggesting that refinement is required to improve description of protein-DNA interactions in AMBER force fields.

Synchrotron-based FTIR microspectroscopy reveals DNA methylation profile in DNA-HALO structure

Fourier transform infrared (FTIR) spectroscopy is a rapid, non-destructive and label-free technique for identifying subtle changes in all bio-macromolecules, and has been used as a method of choice for studying DNA conformation, secondary DNA structure transition and DNA damage. In addition, the specific level of chromatin complexity is introduced via epigenetic modifications forcing the technological upgrade in the analysis of such an intricacy. As the most studied epigenetic mechanism, DNA methylation is a major regulator of transcriptional activity, involved in the suppression of a broad spectrum of genes and its deregulation is involved in all non-communicable diseases. The present study was designed to explore the use of synchrotron-based FTIR analysis to monitor the subtle changes in molecule bases regarding the DNA methylation status of cytosine in the whole genome. In order to reveal the conformation-related best sample for FTIR-based DNA methylation analysis in situ, we used methodology for nuclear HALO preparations and slightly modified it to isolated DNA in HALO formations. Nuclear DNA-HALOs represent samples with preserved higher-order chromatin structure liberated of any protein residues that are closer to native DNA conformation than genomic DNA (gDNA) isolated by the standard batch procedure. Using FTIR spectroscopy we analyzed the DNA methylation profile of isolated gDNA and compared it with the DNA-HALOs. This study demonstrated the potential of FTIR microspectroscopy to detect DNA methylation marks in analyzed DNA-HALO specimens more precisely in comparison with classical DNA extraction procedures that yield unstructured whole genomic DNA. In addition, we used different cell types to assess their global DNA methylation profile, as well as defined specific infrared peaks that can be used for screening DNA methylation.

Non-canonical DNA structures in the human ribosomal DNA

Non-canonical structures (NCS) refer to the various forms of DNA that differ from the B-conformation described by Watson and Crick. It has been found that these structures are usual components of the genome, actively participating in its essential functions. The present review is focused on the nine kinds of NCS appearing or likely to appear in human ribosomal DNA (rDNA): supercoiling structures, R-loops, G-quadruplexes, i-motifs, DNA triplexes, cruciform structures, DNA bubbles, and A and Z DNA conformations. We discuss the conditions of their generation, including their sequence specificity, distribution within the locus, dynamics, and beneficial and detrimental role in the cell.

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.

Delta SARS-CoV-2 s2m Structure, Dynamics, and Entropy: Consequences of the G15U Mutation

Bioinformatic analysis of the Delta SARS-CoV-2 genome reveals a single nucleotide mutation (G15U) in the stem-loop II motif (s2m) relative to ancestral SARS-CoV-2. Despite sequence similarity, unexpected differences between SARS-CoV-2 and Delta SARS-CoV-2 s2m homodimerization experiments require the discovery of unknown structural and thermodynamic changes necessary to rationalize the data. Using our reported SARS-CoV-2 s2m model, we induced the G15U substitution and performed 3.5 microseconds of unbiased molecular dynamics simulation at 283 and 310 K. The resultant Delta s2m adopted a secondary structure consistent with our reported NMR data, resulting in significant deviations in the tertiary structure and dynamics from our SARS-CoV-2 s2m model. First, we find differences in the overall three-dimensional structure, where the characteristic 90° L-shaped kink of the SARS-CoV-2 s2m did not form in the Delta s2m resulting in a "linear" hairpin with limited bending dynamics. Delta s2m helical parameters are calculated to align closely with A-form RNA, effectively eliminating a hinge point to form the L-shape kink by correcting an upper stem defect in SARS-CoV-2 induced by a noncanonical and dynamic G:A base pair. Ultimately, the shape difference rationalizes the migration differences in reported electrophoresis experiments. Second, increased fluctuation of the Delta s2m palindromic sequence, within the terminal loop, compared to SARS-CoV-2 s2m results in an estimated increase of entropy of 6.8 kcal/mol at 310 K relative to the SARS-CoV-2 s2m. The entropic difference offers a unique perspective on why the Delta s2m homodimerizes less spontaneously, forming fewer kissing dimers and extended duplexes compared to SARS-CoV-2. In this work, both the L-shape reduction and palindromic entropic penalty provides an explanation of our reported in vitro electrophoresis homodimerization results. Ultimately, the structural, dynamical, and entropic differences between the SARS-CoV-2 s2m and Delta s2m serve to establish a foundation for future studies of the s2m function in the viral lifecycle.

Structure and Dynamics of DNA and RNA Double Helices Formed by d(CTG), d(GTC), r(CUG), and r(GUC) Trinucleotide Repeats and Associated DNA-RNA Hybrids

Myotonic dystrophy type 1 is the most frequent form of muscular dystrophy in adults caused by an abnormal expansion of the CTG trinucleotide. Both the expanded DNA and the expanded CUG RNA transcript can fold into hairpins. Co-transcriptional formation of stable RNA·DNA hybrids can also enhance the instability of repeat tracts. We performed molecular dynamics simulations of homoduplexes associated with the disease, d(CTG)n and r(CUG)n, and their corresponding r(CAG)n:d(CTG)n and r(CUG)n:d(CAG)n hybrids that can form under bidirectional transcription and of non-pathological d(GTC)n and d(GUC)n homoduplexes. We characterized their conformations, stability, and dynamics and found that the U·U and T·T mismatches are dynamic, favoring anti-anti conformations inside the helical core, followed by anti-syn and syn-syn conformations. For DNA, the secondary minima in the non-expanding d(GTC)n helices are deeper, wider, and longer-lived than those in d(CTG)n, which constitutes another biophysical factor further differentiating the expanding and non-expanding sequences. The hybrid helices are closer to A-RNA, with the A-T and A-U pairs forming two stable Watson-Crick hydrogen bonds. The neutralizing ion distribution around the non-canonical pairs is also described.

Base pair compositional variability influences DNA structural stability and tunes hydration thermodynamics and dynamics

DNA deformability and differential hydration are crucial determinants of biological processes ranging from genetic material packaging to gene expression; their associative details, however, remain inadequately understood. Herein, we report investigations of the dynamic and thermodynamic responses of the local hydration of a variety of base pair sequences. Leveraging in silico sampling and our in-house analyses, we first report the local conformational propensity of sequences that are either predisposed toward the canonical A- or B-conformations or are restrained to potential transitory pathways. It is observed that the transition from the unrestrained A-form to the B-form leads to lengthwise structural deformation. The insertion of intermittent -(CG)- base pairs in otherwise homogeneous -(AT)- sequences bears dynamical consequences for the vicinal hydration layer. Calculation of the excess (pair) entropy suggests substantially higher values of hydration water surrounding A conformations over the B- conformations. Applying the Rosenfeld approximation, we project that the diffusivity of water molecules proximal to canonical B conformation is least for the minor groove of the canonical B-conformation. We determine that structure, composition, and conformation specific groove dimension together influence the local hydration characteristics and, therefore, are expected to be important determinants of biological processes.

Stretching effects on non-adiabatic electron dynamic behavior in poly(dG)-poly(dC) DNA upon the proton irradiation

The electronic excitations caused by DNA when exposed to ion radiation is essential to DNA damage. In this paper, we investigated the energy deposition and electron excitation process of DNA with reasonable stretching range upon proton irradiation based on time-dependent density functional theory. Stretching changes the strength of hydrogen bonding between the DNA base pairs, which in turn affects the Coulomb interaction between the projectile and DNA. As a semi-flexible molecule, the way of energy deposition is weakly sensitive to the stretching rate of DNA. However, the increase of stretching rate causes the increase of charge density along the trajectory channel, sequentially resulting in an increase in proton resistance along the intruding channel. The Mulliken charge analysis indicates that the guanine base and guanine ribose are ionized, meanwhile the cytosine base and cytosine ribose are reduced at all stretching rates. In a few femtoseconds, there exists an electron flow passing through the guanine ribose, guanine, cytosine base and the cytosine ribose in turn. This electron flow increases electron transfer and DNA ionization, promoting the side chain damage of the DNA upon ion irradiation. Our results provide a theoretical insight for deciphering the physical mechanism of the early stage of the irradiation process, and are also of great significance for the study of particle beam cancer therapy in different biological tissues.

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.

Crystal structures of 'ALternative Isoinformational ENgineered' DNA in B-form

The first structural model of duplex DNA reported in 1953 by Watson & Crick presented the double helix in B-form, the form that genomic DNA exists in much of the time. Thus, artificial DNA seeking to mimic the properties of natural DNA should also be able to adopt B-form. Using a host-guest system in which Moloney murine leukemia virus reverse transcriptase serves as the host and DNA as the guests, we determined high-resolution crystal structures of three complexes including 5'-CTTBPPBBSSZZSAAG, 5'-CTTSSPBZPSZBBAAG and 5'-CTTZZPBSBSZPPAAG with 10 consecutive unnatural nucleobase pairs in B-form within self-complementary 16 bp duplex oligonucleotides. We refer to this ALternative Isoinformational ENgineered (ALIEN) genetic system containing two nucleobase pairs (P:Z, pairing 2-amino-imidazo-[1,2-a]-1,3,5-triazin-(8H)-4-one with 6-amino-5-nitro-(1H)-pyridin-2-one, and B:S, 6-amino-4-hydroxy-5-(1H)-purin-2-one with 3-methyl-6-amino-pyrimidin-2-one) as ALIEN DNA. We characterized both position- and sequence-specific helical, nucleobase pair and dinucleotide step parameters of P:Z and B:S pairs in the context of B-form DNA. We conclude that ALIEN DNA exhibits structural features that vary with sequence. Further, Z can participate in alternative stacking modes within a similar sequence context as captured in two different structures. This finding suggests that ALIEN DNA may have a larger repertoire of B-form structures than natural DNA. This article is part of the theme issue 'Reactivity and mechanism in chemical and synthetic biology'.

Desorption of plasmid DNA from anion exchangers: Salt concentration at elution is independent of plasmid size and load

Detailed studies on the sorption behavior of plasmids on anion exchangers are rare compared to proteins. In this study, we systematically compare the elution behavior of plasmid DNA on three common anion exchange resins using linear gradient and isocratic elution experiments. Two plasmids of different lengths, 8 and 20 kbp, were studied and their elution characteristics were compared to a green fluorescent protein. Using established methods for determining retention characteristics of biomolecules in ion exchange chromatography lead to remarkable results. In contrast to the green fluorescent protein, plasmid DNA consistently elutes at one characteristic salt concentration in linear gradient elution. This salt concentration was the same independent of plasmid size but differed slightly for different resins. The behavior is consistent also at preparative loadings of plasmid DNA. Thus, only a single linear gradient elution experiment is sufficient to design elution in a process scale capture step. At isocratic elution conditions, plasmid DNA elutes only above this characteristic concentration. Even at slightly lower concentrations most plasmids remain tightly bound. We hypothesize, that the desorption is accompanied by a conformational change leading to a reduced number of available negative charges for binding. This explanation is supported by structural analysis before and after elution.

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.

Understanding the self-assembly dynamics of A/T absent 'four-way DNA junctions with sticky ends' at altered physiological conditions through molecular dynamics simulations

Elucidation of structure and dynamics of alternative higher-order structures of DNA such as in branched form could be targeted for therapeutics designing. Herein, we are reporting the intrinsically dynamic and folds transitions of an unusual DNA junction with sequence d(CGGCGGCCGC)4 which self-assembles into a four-way DNA junction form with sticky ends using long interval molecular simulations under various artificial physiological conditions. The original crystal structure coordinates (PDB ID: 3Q5C) for the selected DNA junction was considered for a total of 1.1 μs molecular dynamics simulation interval, including different temperature and pH, under OPLS-2005 force field using DESMOND suite. Following, post-dynamics structure parameters for the DNA junction were calculated and analyzed by comparison to the crystal structure. We show here that the self-assembly dynamics of DNA junction is mitigated by the temperature and pH sensitivities, and discloses peculiar structural properties as function of time. From this study it can be concluded on account of temperature sensitive and pH dependent behaviours, DNA junction periodic arrangements can willingly be synthesized and redeveloped for multiple uses like genetic biomarkers, DNA biosensor, DNA nanotechnology, DNA Zipper, etc. Furthermore, the pH dis-regulation behaviour may be used to trigger the functionality of DNA made drug-releasing nanomachines.

Ethanol-Induced Condensation and Decondensation in DNA-Linked Nanoparticles: A Nucleosome-like Model for the Condensed State

Inspired by the conventional use of ethanol to induce DNA precipitation, ethanol condensation has been applied as a routine method to dynamically tune "bond" lengths (i.e., the surface-to-surface distances between adjacent nanoparticles that are linked by DNA) and thermal stabilities of colloidal crystals involving DNA-linked nanoparticles. However, the underlying mechanism of how the DNA bond that links gold nanoparticles changes in this class of colloidal crystals in response to ethanol remains unclear. Here, we conducted a series of all-atom molecular dynamic (MD) simulations to explore the free energy landscape for DNA condensation and decondensation. Our simulations confirm that DNA condensation is energetically much more favorable under 80% ethanol conditions than in pure water, as a result of ethanol's role in enhancing electrostatic interactions between oppositely charged species. Moreover, the condensed DNA adopts B-form in pure water and A-form in 80% ethanol, which indicates that the higher-order transition does not affect DNA's conformational preferences. We further propose a nucleosome-like supercoiled model for the DNA condensed state, and we show that the DNA end-to-end distance derived from this model matches the experimentally measured DNA bond length of about 3 nm in the fully condensed state for DNA where the measured length is 16 nm in water. Overall, this study provides an atomistic understanding of the mechanism underlying ethanol-induced condensation and water-induced decondensation, while our proposed nucleosome-like model allows the design of new strategies for interpreting experimental studies of DNA condensation.

RNA as a Major-Groove Ligand: RNA-RNA and RNA-DNA Triplexes Formed by GAA and UUC or TTC Sequences

Friedreich's ataxia is associated with noncanonical nucleic acid structures that emerge when GAA:TTC repeats in the first intron of the FXN gene expand beyond a critical number of repeats. Specifically, the noncanonical repeats are associated with both triplexes and R-loops. Here, we present an in silico investigation of all possible triplexes that form by attaching a third RNA strand to an RNA:RNA or DNA:DNA duplex, complementing previous DNA-based triplex studies. For both new triplexes results are similar. For a pyridimine UUC⁺ third strand, the parallel orientation is stable while its antiparallel counterpart is unstable. For a neutral GAA third strand, the parallel conformation is stable. A protonated GA⁺A third strand is stable in both parallel and antiparallel orientations. We have also investigated Na⁺ and Mg²⁺ ion distributions around the triplexes. The presence of Mg²⁺ ions helps stabilize neutral, antiparallel GAA triplexes. These results (along with previous DNA-based studies) allow for the emergence of a complete picture of the stability and structural characteristics of triplexes based on the GAA and TTC/UUC sequences, thereby contributing to the field of trinucleotide repeats and the associated unusual structures that trigger expansion.

Visualizing "Alternative Isoinformational Engineered" DNA in A- and B-Forms at High Resolution

A fundamental property of DNA built from four informational nucleotide units (GCAT) is its ability to adopt different helical forms within the context of the Watson-Crick pair. Well-characterized examples include A-, B-, and Z-DNA. For this study, we created an isoinformational biomimetic polymer, built (like standard DNA) from four informational "letters", but with the building blocks being artificial. This ALternative Isoinformational ENgineered (ALIEN) DNA was hypothesized to support two nucleobase pairs, the P:Z pair matching 2-amino-imidazo-[1,2a]-1,3,5-triazin-[8H]-4-one with 6-amino-3-5-nitro-1H-pyridin-2-one and the B:S pair matching 6-amino-4-hydroxy-5-1H-purin-2-one with 3-methyl-6-amino-pyrimidin-2-one. We report two structures of ALIEN DNA duplexes at 1.2 Å resolution and a third at 1.65 Å. All of these are built from a single self-complementary sequence (5'-CTSZZPBSBSZPPBAG) that includes 12 consecutive ALIEN nucleotides. We characterized the helical, nucleobase pair, and dinucleotide step parameters of ALIEN DNA in these structures. In addition to showing that ALIEN pairs retain basic Watson-Crick pairing geometry, two of the ALIEN DNA structures are characterized as A-form DNA and one as B-form DNA. We identified parameters that map differences effecting the transition between the two helical forms; these same parameters distinguish helical forms of isoinformational natural DNA. Collectively, our analyses suggest that ALIEN DNA retains essential structural features of natural DNA, not only its information density and Watson-Crick pairing but also its ability to adopt two canonical forms.

Saffron, Its Active Components, and Their Association with DNA and Histone Modification: A Narrative Review of Current Knowledge

Intensive screening for better and safer medications to treat diseases such as cancer and inflammatory diseases continue, and some phytochemicals have been discovered to have anti-cancer and many therapeutical activities. Among the traditionally used spices, Crocus sativus (saffron) and its principal bioactive constituents have anti-inflammatory, antioxidant, and chemopreventive properties against multiple malignancies. Early reports have shown that the epigenetic profiles of healthy and tumor cells vary significantly in the context of different epigenetic factors. Multiple components, such as carotenoids as bioactive dietary phytochemicals, can directly or indirectly regulate epigenetic factors and alter gene expression profiles. Previous reports have shown the interaction between active saffron compounds with linker histone H1. Other reports have shown that high concentrations of saffron bind to the minor groove of calf thymus DNA, resulting in specific structural changes from B- to C-form of DNA. Moreover, the interaction of crocin G-quadruplex was reported. A recent in silico study has shown that residues of SIRT1 interact with saffron bio-active compounds and might enhance SIRT1 activation. Other reports have shown that the treatment of Saffron bio-active compounds increases γH2AX, decreases HDAC1 and phosphorylated histone H3 (p-H3). However, the question that still remains to be addressed how saffron triggers various epigenetic changes? Therefore, this review discusses the literature published till 2022 regarding saffron as dietary components and its impact on epigenetic mechanisms. Novel bioactive compounds such as saffron components that lead to epigenetic alterations might be a valuable strategy as an adjuvant therapeutic drug.

B-to-A transition in target DNA during retroviral integration

Integration into host target DNA (tDNA), a hallmark of retroviral replication, is mediated by the intasome, a multimer of integrase (IN) assembled on viral DNA (vDNA) ends. To ascertain aspects of tDNA recognition during integration, we have solved the 3.5 Å resolution cryo-EM structure of the mouse mammary tumor virus (MMTV) strand transfer complex (STC) intasome. The tDNA adopts an A-like conformation in the region encompassing the sites of vDNA joining, which exposes the sugar-phosphate backbone for IN-mediated strand transfer. Examination of existing retroviral STC structures revealed conservation of A-form tDNA in the analogous regions of these complexes. Furthermore, analyses of sequence preferences in genomic integration sites selectively targeted by six different retroviruses highlighted consistent propensity for A-philic sequences at the sites of vDNA joining. Our structure additionally revealed several novel MMTV IN-DNA interactions, as well as contacts seen in prior STC structures, including conserved Pro125 and Tyr149 residues interacting with tDNA. In infected cells, Pro125 substitutions impacted the global pattern of MMTV integration without significantly altering local base sequence preferences at vDNA insertion sites. Collectively, these data advance our understanding of retroviral intasome structure and function, as well as factors that influence patterns of vDNA integration in genomic DNA.

Identification of sugar-containing natural products that interact with i-motif DNA

There are thousands of compounds shown to interact with G-quadruplex DNA, yet very few which target i-motif (iM) DNA. Previous work showed that tobramycin can interact with iM- DNA, indicating the potential for sugar-molecules to target these structures. Computational approaches indicated that the sugar-containing natural products baicalin and geniposidic acid had potential to target iM-DNA. We assessed the DNA interacting properties of these compounds using FRET-based DNA melting and a fluorescence-based displacement assay using iM-DNA structures from the human telomere and the insulin linked polymorphic region (ILPR), as well as complementary G-quadruplex and double stranded DNA. Both baicalin and geniposidic acid show promise as iM-interacting compounds with potential for use in experiments into the structure and function of i-motif forming DNA sequences and present starting points for further synthetic development of these as probes for iM-DNA.

Revisiting DNA Sequence-Dependent Deformability in High-Resolution Structures: Effects of Flanking Base Pairs on Dinucleotide Morphology and Global Chain Configuration

DNA carries more than the list of biochemical ingredients that drive the basic functions of living systems. The sequence of base pairs includes a multitude of structural and energetic signals, which determine the degree to which the long, threadlike molecule moves and how it responds to proteins and other molecules that control its processing and govern its packaging. The chemical composition of base pairs directs the spatial disposition and fluctuations of successive residues. The observed arrangements of these moieties in high-resolution protein–DNA crystal structures provide one of the best available estimates of the natural, sequence-dependent structure and deformability of the double-helical molecule. Here, we update the set of knowledge-based elastic potentials designed to describe the observed equilibrium structures and configurational fluctuations of the ten unique base-pair steps. The large number of currently available structures makes it possible to characterize the configurational preferences of the DNA base-pair steps within the context of their immediate neighbors, i.e., tetrameric context. Use of these knowledge-based potentials shows promise in accounting for known effects of sequence in long chain molecules, e.g., the degree of curvature reported in classic gel mobility studies and the recently reported sequence-dependent responses of supercoiled minicircles to nuclease cleavage.

Selective Photo-Crosslinking Detection of Methylated Cytosine in DNA Duplex Aided by a Cationic Comb-Type Copolymer

In the process of cell development and differentiation, C-5-methylation of cytosine (5-methylcytosine: 5-mC) in genome DNA is an important transcriptional regulator that switches between differentiated and undifferentiated states. Further, abnormal DNA methylations are often present in tumor suppressor genes and are associated with many diseases. Therefore, 5-mC detection technology is an important tool in the most exciting fields of molecular biology and diagnosing diseases such as cancers. In this study, we found a novel photo-crosslinking property of psoralen-conjugated oligonucleotide (Ps-Oligo) to the double-stranded DNA (ds-DNA) containing 5-mC in the presence of a cationic comb-type copolymer, poly(allylamine)-graft-dextran (PAA-g-Dex). Photo-crosslinking efficiency of Ps-Oligo to 5-mC in ds-DNA was markedly enhanced in the presence of PAA-g-Dex, permitting 5-mC-targeted crosslinking. We believe that the combination of PAA-g-Dex and Ps-Oligo will be an effective tool for detecting 5-mC in genomic DNA.

Crystal structure of an RNA/DNA strand exchange junction

Short segments of RNA displace one strand of a DNA duplex during diverse processes including transcription and CRISPR-mediated immunity and genome editing. These strand exchange events involve the intersection of two geometrically distinct helix types-an RNA:DNA hybrid (A-form) and a DNA:DNA homoduplex (B-form). Although previous evidence suggests that these two helices can stack on each other, it is unknown what local geometric adjustments could enable A-on-B stacking. Here we report the X-ray crystal structure of an RNA-5'/DNA-3' strand exchange junction at an anisotropic resolution of 1.6 to 2.2 Å. The structure reveals that the A-to-B helical transition involves a combination of helical axis misalignment, helical axis tilting and compression of the DNA strand within the RNA:DNA helix, where nucleotides exhibit a mixture of A- and B-form geometry. These structural principles explain previous observations of conformational stability in RNA/DNA exchange junctions, enabling a nucleic acid architecture that is repeatedly populated during biological strand exchange events.

Nucleosome recognition and DNA distortion by the Chd1 remodeler in a nucleotide-free state

Chromatin remodelers are ATP-dependent enzymes that reorganize nucleosomes within all eukaryotic genomes. Here we report a complex of the Chd1 remodeler bound to a nucleosome in a nucleotide-free state, determined by cryo-EM to 2.3 Å resolution. The remodeler stimulates the nucleosome to absorb an additional nucleotide on each strand at two different locations: on the tracking strand within the ATPase binding site and on the guide strand one helical turn from the ATPase motor. Remarkably, the additional nucleotide on the tracking strand is associated with a local transformation toward an A-form geometry, explaining how sequential ratcheting of each DNA strand occurs. The structure also reveals a histone-binding motif, ChEx, which can block opposing remodelers on the nucleosome and may allow Chd1 to participate in histone reorganization during transcription.

Structural basis of SNAPc-dependent snRNA transcription initiation by RNA polymerase II

RNA polymerase II (Pol II) carries out transcription of both protein-coding and non-coding genes. Whereas Pol II initiation at protein-coding genes has been studied in detail, Pol II initiation at non-coding genes, such as small nuclear RNA (snRNA) genes, is less well understood at the structural level. Here, we study Pol II initiation at snRNA gene promoters and show that the snRNA-activating protein complex (SNAPc) enables DNA opening and transcription initiation independent of TFIIE and TFIIH in vitro. We then resolve cryo-EM structures of the SNAPc-containing Pol IIpre-initiation complex (PIC) assembled on U1 and U5 snRNA promoters. The core of SNAPc binds two turns of DNA and recognizes the snRNA promoter-specific proximal sequence element (PSE), located upstream of the TATA box-binding protein TBP. Two extensions of SNAPc, called wing-1 and wing-2, bind TFIIA and TFIIB, respectively, explaining how SNAPc directs Pol II to snRNA promoters. Comparison of structures of closed and open promoter complexes elucidates TFIIH-independent DNA opening. These results provide the structural basis of Pol II initiation at non-coding RNA gene promoters.

Structural switching/polymorphism by sequential base substitution at quasi-palindromic SNP site (G → A) in LCR of human β-globin gene cluster

The human β-globin gene Locus Control Region (LCR), a dominant regulator of globin gene expression contains five tissue-specific DNase I-hypersensitive sites (HSs). A single nucleotide polymorphism (SNP) (A → G) present in HS4 region of locus control region (LCR), have shown a notable association between the G allele and the occurrence of β-thalassemia. This SNP site exhibiting a hairpin - duplex equilibrium manifested in A → B like DNA transition has previously been reported from this laboratory. Since, DNA is a dynamic and adaptable molecule, so any change of a single base within a primary DNA sequence can produce major biological consequences commonly manifested in genetic disorders such as sickle cell anemia and β-thalassemia. Herein, the differential behavior of sequential single base substitutions G → A on the quasi-palindromic sequence (d-TGGGGGCCCCA; HPG11) has been explored. A combination of native gel electrophoresis, circular dichroism (CD), and UV-thermal denaturation (Tm) techniques have been used to investigate the structural polymorphism associated with various variants of HPG11 i.e. HPG11A2 to HPG11A5. The CD spectra confirmed that all the HPG11 variants exhibit a hairpin - duplex equilibrium. Oligomer concentration dependence on CD spectra has been correlated with A → B DNA conformational transition. However, as revealed in gel electrophoresis, HPG11A2 → A5 exhibit the formation of a tetramolecular structure (four-way junction) at higher oligomer concentration. UV-melting studies also supported the melting of hairpin, duplex and four-way junction structure. This polymorphism pattern may possibly be significant for DNA-protein recognition, in the process of regulation of LCR in the β-globin gene.

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.

Accurate prediction of B-form/A-form DNA conformation propensity from primary sequence: A machine learning and free energy handshake

DNA carries the genetic code of life, with different conformations associated with different biological functions. Predicting the conformation of DNA from its primary sequence, although desirable, is a challenging problem owing to the polymorphic nature of DNA. We have deployed a host of machine learning algorithms, including the popular state-of-the-art LightGBM (a gradient boosting model), for building prediction models. We used the nested cross-validation strategy to address the issues of “overfitting” and selection bias. This simultaneously provides an unbiased estimate of the generalization performance of a machine learning algorithm and allows us to tune the hyperparameters optimally. Furthermore, we built a secondary model based on SHAP (SHapley Additive exPlanations) that offers crucial insight into model interpretability. Our detailed model-building strategy and robust statistical validation protocols tackle the formidable challenge of working on small datasets, which is often the case in biological and medical data.

•

A robust machine learning model to predict A- or B-DNA conformation

•

Outcome of machine learning model is explained with free energy values

•

Our approach works well under class imbalance and limited data constraints

The sequence in the genome of an organism encodes all the information of life. We combine a data-driven approach using machine learning (ML) and the results of free energy calculations to offer a fresh perspective on this long-standing problem of prediction of DNA conformation (A or B) from the sequence. We trained our ML model using sophisticated state-of-the art algorithms such as LightGBM along with a nested cross-validation strategy to overcome the common problems associated with data bias and overfitting when constrained by limited data size. Our study will serve the broader interest of researchers who are not only seeking accurate and reliable predictive models but also want to understand the physical and chemical origins behind the predictions.

A molecular view of DNA flexibility

DNA dynamics can only be understood by taking into account its complex mechanical behavior at different length scales. At the micrometer level, the mechanical properties of single DNA molecules have been well-characterized by polymer models and are commonly quantified by a persistence length of 50 nm (~150 bp). However, at the base pair level (~3.4 Å), the dynamics of DNA involves complex molecular mechanisms that are still being deciphered. Here, we review recent single-molecule experiments and molecular dynamics simulations that are providing novel insights into DNA mechanics from such a molecular perspective. We first discuss recent findings on sequence-dependent DNA mechanical properties, including sequences that resist mechanical stress and sequences that can accommodate strong deformations. We then comment on the intricate effects of cytosine methylation and DNA mismatches on DNA mechanics. Finally, we review recently reported differences in the mechanical properties of DNA and double-stranded RNA, the other double-helical carrier of genetic information. A thorough examination of the recent single-molecule literature permits establishing a set of general 'rules' that reasonably explain the mechanics of nucleic acids at the base pair level. These simple rules offer an improved description of certain biological systems and might serve as valuable guidelines for future design of DNA and RNA nanostructures.

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.

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

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

Function of a viral genome packaging motor from bacteriophage T4 is insensitive to DNA sequence

Many viruses employ ATP-powered motors during assembly to translocate DNA into procapsid shells. Previous reports raise the question if motor function is modulated by substrate DNA sequence: (i) the phage T4 motor exhibits large translocation rate fluctuations and pauses and slips; (ii) evidence suggests that the phage phi29 motor contacts DNA bases during translocation; and (iii) one theoretical model, the 'B-A scrunchworm', predicts that 'A-philic' sequences that transition more easily to A-form would alter motor function. Here, we use single-molecule optical tweezers measurements to compare translocation of phage, plasmid, and synthetic A-philic, GC rich sequences by the T4 motor. We observed no significant differences in motor velocities, even with A-philic sequences predicted to show higher translocation rate at high applied force. We also observed no significant changes in motor pausing and only modest changes in slipping. To more generally test for sequence dependence, we conducted correlation analyses across pairs of packaging events. No significant correlations in packaging rate, pausing or slipping versus sequence position were detected across repeated measurements with several different DNA sequences. These studies suggest that viral genome packaging is insensitive to DNA sequence and fluctuations in packaging motor velocity, pausing and slipping are primarily stochastic temporal events.

Atypical structures of GAA/TTC trinucleotide repeats underlying Friedreich's ataxia: DNA triplexes and RNA/DNA hybrids

Expansion of the GAA/TTC repeats in the first intron of the FXN gene causes Friedreich's ataxia. Non-canonical structures are linked to this expansion. DNA triplexes and R-loops are believed to arrest transcription, which results in frataxin deficiency and eventual neurodegeneration. We present a systematic in silico characterization of the possible DNA triplexes that could be assembled with GAA and TTC strands; the two hybrid duplexes [r(GAA):d(TTC) and d(GAA):r(UUC)] in an R-loop; and three hybrid triplexes that could form during bidirectional transcription when the non-template DNA strand bonds with the hybrid duplex (collapsed R-loops, where the two DNA strands remain antiparallel). For both Y·R:Y and R·R:Y DNA triplexes, the parallel third strand orientation is more stable; both parallel and antiparallel protonated d(GA+A)·d(GAA):d(TTC) triplexes are stable. Apparent contradictions in the literature about the R·R:Y triplex stability is probably due to lack of molecular resolution, since shifting the third strand by a single nucleotide alters the stability ranking. In the collapsed R-loops, antiparallel d(TTC+)·d(GAA):r(UUC) is unstable, while parallel d(GAA)·r(GAA):d(TTC) and d(GA+A)·r(GAA):d(TTC) are stable. In addition to providing new structural perspectives for specific therapeutic aims, our results contribute to a systematic structural basis for the emerging field of quantitative R-loop biology.

Stacking geometry between two sheared Watson-Crick basepairs: Computational chemistry and bioinformatics based prediction

BACKGROUND

Molecular modeling of RNA double helices is possible using most probable values of basepair parameters obtained from crystal structure database. The A:A w:wC non-canonical basepair, involving Watson-Crick edges of two Adenines in cis orientation, appears quite frequently in database. Bimodal distribution of its Shear, due to two different H-bonding schemes, introduces the confusion in assigning most the probable value. Its effect is pronounced when the A:A w:wC basepair stacks on Sheared wobble G:U W:WC basepairs.

METHODS

We employed molecular dynamics simulations of three possible double helices with GAG, UAG and GAU sequence motifs at their centers and quantum chemical calculation for non-canonical A:A w:wC basepair stacked on G:U W:WC basepair.

RESULTS

We noticed stable structures of GAG motif with specifically negative Shear of the A:A basepair but stabilities of the other motifs were not found with A:A w:wC basepairing. Hybrid DFT-D and MP2 stacking energy analyses on dinucleotide step sequences, A:A w:wC::G:U W:WC and A:A w:wC::U:G W:WC reveal that viable orientation of A:A::G:U prefers one of the H-bonding modes with negative Shear, supported by crystal structure database. The A:A::U:G dinucleotide, however, prefers structure with only positive Shear.

CONCLUSIONS

The quantum chemical calculations explain why MD simulations of GAG sequence motif only appear stable. In the cases of the GAU and UAG motifs "tug of war" situation between positive and negative Shears of A:A w:wC basepair induces conformational plasticity.

GENERAL SIGNIFICANCE

We have projected comprehensive reason behind the promiscuous nature of A:A w:wC basepair which brings occasional structural plasticity.

Orthogonal Genetic Systems

Xenobiology is an emerging area of synthetic biology that aims to safeguard genetically engineered cells by storing synthetic biology information in xeno-nucleic acid polymers (XNAs). Critical to the success of this effort is the need to establish cellular systems that can maintain an XNA chromosome in actively dividing cells. This viewpoint discusses the structural parameters of the nucleic acid backbone that should be considered when designing an orthogonal genetic system that can replicate without interference from the endogenous genome. In addition to practical value, these studies have the potential to provide new fundamental insight into the structure and function properties of unnatural nucleic acid polymers.

Web 3DNA 2.0 for the analysis, visualization, and modeling of 3D nucleic acid structures

Web 3DNA (w3DNA) 2.0 is a significantly enhanced version of the widely used w3DNA server for the analysis, visualization, and modeling of 3D nucleic-acid-containing structures. Since its initial release in 2009, the w3DNA server has continuously served the community by making commonly-used features of the 3DNA suite of command-line programs readily accessible. However, due to the lack of updates, w3DNA has clearly shown its age in terms of modern web technologies and it has long lagged behind further developments of 3DNA per se. The w3DNA 2.0 server presented here overcomes all known shortcomings of w3DNA while maintaining its battle-tested characteristics. Technically, w3DNA 2.0 implements a simple and intuitive interface (with sensible defaults) for increased usability, and it complies with HTML5 web standards for broad accessibility. Featurewise, w3DNA 2.0 employs the most recent version of 3DNA, enhanced with many new functionalities, including: the automatic handling of modified nucleotides; a set of 'simple' base-pair and step parameters for qualitative characterization of non-Watson-Crick double-helical structures; new structural parameters that integrate the rigid base plane and the backbone phosphate group, the two nucleic acid components most reliably determined with X-ray crystallography; in silico base mutations that preserve the backbone geometry; and a notably improved module for building models of single-stranded RNA, double-helical DNA, Pauling triplex, G-quadruplex, or DNA structures 'decorated' with proteins. The w3DNA 2.0 server is freely available, without registration, at http://web.x3dna.org.

DNA Conformational Changes Play a Force-Generating Role during Bacteriophage Genome Packaging

Motors that move DNA, or that move along DNA, play essential roles in DNA replication, transcription, recombination, and chromosome segregation. The mechanisms by which these DNA translocases operate remain largely unknown. Some double-stranded DNA (dsDNA) viruses use an ATP-dependent motor to drive DNA into preformed capsids. These include several human pathogens as well as dsDNA bacteriophages-viruses that infect bacteria. We previously proposed that DNA is not a passive substrate of bacteriophage packaging motors but is instead an active component of the machinery. We carried out computational studies on dsDNA in the channels of viral portal proteins, and they reveal DNA conformational changes consistent with that hypothesis. dsDNA becomes longer ("stretched") in regions of high negative electrostatic potential and shorter ("scrunched") in regions of high positive potential. These results suggest a mechanism that electrostatically couples the energy released by ATP hydrolysis to DNA translocation: The chemical cycle of ATP binding, hydrolysis, and product release drives a cycle of protein conformational changes. This produces changes in the electrostatic potential in the channel through the portal, and these drive cyclic changes in the length of dsDNA as the phosphate groups respond to the protein's electrostatic potential. The DNA motions are captured by a coordinated protein-DNA grip-and-release cycle to produce DNA translocation. In short, the ATPase, portal, and dsDNA work synergistically to promote genome packaging.

Changes in Microenvironment Modulate the B- to A-DNA Transition

B- to A-DNA transition is known to be sensitive to the macroscopic properties of the solution, such as salt and ethanol concentrations. Microenvironmental effects on DNA conformational transition have been broadly studied. Providing an intuitive picture of how DNA responds to environmental changes is, however, still needed. Analyzing the chemical equilibrium of B-to-A DNA transition at critical concentrations, employing explicit-solvent simulations, is envisioned to help understand such microenvironmental effects. In the present study, free-energy calculations characterizing the B- to A-DNA transition and the distribution of cations were carried out in solvents with different ethanol concentrations. With the addition of ethanol, the most stable structure of DNA changes from the B- to A-form, in agreement with previous experimental observation. In 60% ethanol, a chemical equilibrium is found, showing reversible transition between B- and A-DNA. Analysis of the microenvironment around DNA suggests that with the increase of ethanol concentration, the cations exhibit a significant tendency to move toward the backbone, and mobility of water molecules around the major groove and backbone decreases gradually, leading eventually to a B-to-A transition. The present results provide a free-energy view of DNA microenvironment and of the role of cation motion in the conformational transition.

DNA Crookedness Regulates DNA Mechanical Properties at Short Length Scales

Sequence-dependent DNA conformation and flexibility play a fundamental role in the specificity of DNA-protein interactions. Here we quantify the DNA crookedness: a sequence-dependent deformation of DNA that consists of periodic bends of the base pair centers chain. Using extensive 100  μs-long, all-atom molecular dynamics simulations, we found that DNA crookedness and its associated flexibility are bijective, which unveils a one-to-one relation between DNA structure and dynamics. This allowed us to build a predictive model to compute the stretch moduli of different DNA sequences from solely their structure. Sequences with very little crookedness show extremely high stretching stiffness and have been previously shown to form unstable nucleosomes and promote gene expression. Interestingly, the crookedness can be tailored by epigenetic modifications, known to affect gene expression. Our results rationalize the idea that the DNA sequence is not only a chemical code, but also a physical one that allows finely regulating its mechanical properties and, possibly, its 3D arrangement inside the cell.

Unusual stability exhibited by (AT)XN12(AT)Y motif associated with high fetal hemoglobin levels

Quasi-palindromic sequences (AT)_XN₁₂(AT)_Y present in HS2 (hypersensitive site 2) of the human β-globin locus are known to be significantly associated with increased fetal hemoglobin (HbF) levels. High HbF levels in some adults arise due to pathological conditions such as sickle cell disease and β-thalassemia. However, elevated levels of HbF are also associated with a reducing morbidity and mortality in patients with β-thalassemia and thus ameliorate the severity of the disease. Using gel-electrophoresis, ultraviolet (UV)-thermal denaturation, and circular dichroism (CD) techniques, we demonstrated that it exhibits a hairpin-duplex equilibrium. Intramolecular species (hairpin) were observed in both low and high salt concentrations in gel assay studies displaying the unusual stability of intramolecular species even at the high counter-ion concentration. The unusual stability of hairpin secondary structures was also demonstrated by the monophasic nature of the melting profiles for the oligonucleotides which persisted at low as well as high salt and oligomer concentrations. Change in CD spectra as a function of oligomer concentration indicates that the bimolecular duplex formation is selectively favored over monomolecular hairpin formation at and above 9 µM oligomer concentration. Thus, we hypothesize that imperfect inverted repeat sequence (AT)_XN₁₂(AT)_Y of HS2 of β-globin gene LCR forms the unusually stable hairpins which may result in the formation of a cruciform structure that may be recruited for binding by various nuclear proteins that could result in elevated HbF levels. Communicated by Ramaswamy H. Sarma.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Parallel Guanine Duplex and Cytosine Duplex DNA with Uninterrupted Spines of AgI-Mediated Base Pairs

Hydrogen bonding between nucleobases produces diverse DNA structural motifs, including canonical duplexes, guanine (G) quadruplexes, and cytosine (C) i-motifs. Incorporating metal-mediated base pairs into nucleic acid structures can introduce new functionalities and enhanced stabilities. Here we demonstrate, using mass spectrometry (MS), ion mobility spectrometry (IMS), and fluorescence resonance energy transfer (FRET), that parallel-stranded structures consisting of up to 20 G-Ag^I-G contiguous base pairs are formed when natural DNA sequences are mixed with silver cations in aqueous solution. FRET indicates that duplexes formed by poly(cytosine) strands with 20 contiguous C-Ag^I-C base pairs are also parallel. Silver-mediated G-duplexes form preferentially over G-quadruplexes, and the ability of Ag⁺ to convert G-quadruplexes into silver-paired duplexes may provide a new route to manipulating these biologically relevant structures. IMS indicates that G-duplexes are linear and more rigid than B-DNA. DFT calculations were used to propose structures compatible with the IMS experiments. Such inexpensive, defect-free, and soluble DNA-based nanowires open new directions in the design of novel metal-mediated DNA nanotechnology.

Strong deformations of DNA: Effect on the persistence length

Extreme deformations of the DNA double helix attracted a lot of attention during the past decades. Particularly, the determination of the persistence length of DNA with extreme local disruptions, or kinks, has become a crucial problem in the studies of many important biological processes. In this paper we review an approach to calculate the persistence length of the double helix by taking into account the formation of kinks of arbitrary configuration. The reviewed approach improves the Kratky-Porod model to determine the type and nature of kinks that occur in the double helix, by measuring a reduction of the persistence length of the kinkable DNA.

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.

Role of Inosine⁻Uracil Base Pairs in the Canonical RNA Duplexes

Adenosine to inosine (A–I) editing is the most common modification of double-stranded RNA (dsRNA). This change is mediated by adenosine deaminases acting on RNA (ADARs) enzymes with a preference of U>A>C>G for 5′ neighbor and G>C=A>U or G>C>U=A for 3′ neighbor. A–I editing occurs most frequently in the non-coding regions containing repetitive elements such as ALUs. It leads to disruption of RNA duplex structure, which prevents induction of innate immune response. We employed standard and biased molecular dynamics (MD) simulations to analyze the behavior of RNA duplexes with single and tandem inosine–uracil (I–U) base pairs in different sequence context. Our analysis showed that the I–U pairs induce changes in base pair and base pair step parameters and have different dynamics when compared with standard canonical base pairs. In particular, the first I–U pair from tandem I–U/I–U systems exhibited increased dynamics depending on its neighboring 5′ base. We discovered that UII sequence, which is frequently edited, has lower flexibility compared with other sequences (AII, GII, CII), hence it only modestly disrupts dsRNA. This might indicate that the UAA motifs in ALUs do not have to be sufficiently effective in preventing immune signaling.