A standard reference frame for the description of nucleic acid base-pair geometry

The Contribution of Backbone Electrostatic Repulsion to DNA Mechanical Properties is Length-Scale-Dependent

The mechanics of DNA bending is crucially related to many vital biological processes. Recent experiments reported anomalous flexibility for DNA on short length scales, calling into doubt the validity of the harmonic worm-like chain (WLC) model in this region. In the present work, we systematically probed the bending dynamics of DNA at different length scales. In contrast to the remarkable deviation from the WLC description for DNA duplexes of less than three helical turns, our atomistic studies indicate that the neutral "null isomer" behaves in accord with the ideal elastic WLC and exhibits a uniform decay for the directional correlation of local bending. The backbone neutralization weakens the anisotropy in the effective bending preference and the helical periodicity of bend correlation that have previously been observed for normal DNA. The contribution of electrostatic repulsion to stretching cooperativity and the mechanical properties of DNA strands is length-scale-dependent: the phosphate neutralization increases the stiffness of DNA below two helical turns, but it is decreased for longer strands. We find that DNA rigidity is largely determined by base pair stacking, with electrostatic interactions contributing only around 10% of the total persistence length.

cgDNAweb: a web interface to the cgDNA sequence-dependent coarse-grain model of double-stranded DNA

The sequence-dependent statistical mechanical properties of fragments of double-stranded DNA is believed to be pertinent to its biological function at length scales from a few base pairs (or bp) to a few hundreds of bp, e.g. indirect read-out protein binding sites, nucleosome positioning sequences, phased A-tracts, etc. In turn, the equilibrium statistical mechanics behaviour of DNA depends upon its ground state configuration, or minimum free energy shape, as well as on its fluctuations as governed by its stiffness (in an appropriate sense). We here present cgDNAweb, which provides browser-based interactive visualization of the sequence-dependent ground states of double-stranded DNA molecules, as predicted by the underlying cgDNA coarse-grain rigid-base model of fragments with arbitrary sequence. The cgDNAweb interface is specifically designed to facilitate comparison between ground state shapes of different sequences. The server is freely available at cgDNAweb.epfl.ch with no login requirement.

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.

Comparison of DNA and RNA substrate effects on TET2 structure

Ten-eleven translocation (TET) enzymes can perform the stepwise oxidation of 5-methylcytosine (5mC) to 5-carboxylcytosine on both single-stranded (ss) and double-stranded (ds) DNA and RNA. It has been established that TET2 has a preference for ds DNA substrates, but it can catalyze the oxidation reaction on both ssDNA and RNA. The reasons for this substrate preference have been investigated for only a substrate 5mC ribonucleotide in a DNA strand, but not other nucleic acid configurations (Biochemistry58 (2019) 411). We performed molecular dynamics simulations on TET2 with various ss and ds substrates in order to better understand the structural and dynamical reasons for TET2's preference to act on ds DNA. Our simulations show that substrates that have a ribonucleotide experience several disruptions in their overall backbone shape, hydrogen bonding character, and non-bonded interactions. These differences appear to lead to the instability of ribonucleotide in the active site, and provide further rational for TET2's experimental behavior.

Effects of Noncanonical Base Pairing on RNA Folding: Structural Context and Spatial Arrangements of G·A Pairs

Noncanonical base pairs play important roles in assembling the three-dimensional structures critical to the diverse functions of RNA. These associations contribute to the looped segments that intersperse the canonical double-helical elements within folded, globular RNA molecules. They stitch together various structural elements, serve as recognition elements for other molecules, and act as sites of intrinsic stiffness or deformability. This work takes advantage of new software (DSSR) designed to streamline the analysis and annotation of RNA three-dimensional structures. The multiscale structural information gathered for individual molecules, combined with the growing number of unique, well-resolved RNA structures, makes it possible to examine the collective features deeply and to uncover previously unrecognized patterns of chain organization. Here we focus on a subset of noncanonical base pairs involving guanine and adenine and the links between their modes of association, secondary structural context, and contributions to tertiary folding. The rigorous descriptions of base-pair geometry that we employ facilitate characterization of recurrent geometric motifs and the structural settings in which these arrangements occur. Moreover, the numerical parameters hint at the natural motions of the interacting bases and the pathways likely to connect different spatial forms. We draw attention to higher-order multiplexes involving two or more G·A pairs and the roles these associations appear to play in bridging different secondary structural units. The collective data reveal pairing propensities in base organization, secondary structural context, and deformability and serve as a starting point for further multiscale investigations and/or simulations of RNA folding.

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.

Overtwisting induces polygonal shapes in bent DNA

By combining analytical results and simulations of various coarse-grained models, we investigate the minimal energy shape of DNA minicircles which are torsionally constrained by an imposed over or undertwist. We show that twist-bend coupling, a cross interaction term discussed in the recent DNA literature, induces minimal energy shapes with a periodic alternation of parts with high and low curvature resembling rounded polygons. We briefly discuss the possible experimental relevance of these findings. We finally show that the twist and bending energies of minicircles are governed by renormalized stiffness constants, rather than the bare ones. This has important consequences for the analysis of experiments involving circular DNA meant to determine DNA elastic constants.

Impact of N-(2-aminoethyl) Glycine Unit on Watson-Crick Base Pairs

We aim to understand the structure and stability of the backbone tailored Watson-Crick base pairs, Guanine-Cytosine (GC), Adenine-Thymine (AT) and Adenine-Uracil (AU) by incorporating N-(2-aminoethyl) glycine units (linked by amide bonds) at the purine and pyrimidine sites of the nucleobases. Density functional theory (DFT) is employed in which B3LYP/6-311++G∗ ∗ level of theory has been used to optimize all the structures. The peptide attached base pairs are compared with the natural deoxyribose nucleic acid (DNA)/ribonucleic acid (RNA) base pairs and the calculations are carried out in both the gas and solution phases. The structural propensities of the optimized base pairs are analyzed using base pair geometries, hydrogen bond distances and stabilization energies and, compared with the standard reference data. The structural parameters were found to correlate well with the available data. The addition of peptide chain at the back bone of the DNA/RNA base pairs results only with a minimal distortion and hence does not alter the structural configuration of the base pairs. Also enhanced stability of the base pairs is spotted while adding peptidic chain at the purine site rather than the pyrimidine site of the nucleobases. The stability of the complexes is further interpreted by considering the hydrogen bonded N–H stretching frequencies of the respective base pairs. The discrimination in the interaction energies observed in both gas and solution phases are resulted due to the existence of distinct lowest unoccupied molecular orbitals (LUMO) in the solution phase. The reactivity of the base pairs is also analyzed through the in-depth examinations on the highest occupied molecular orbital (HOMO)-LUMO orbitals.

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.

Evidence for Hidden Involvement of N3-Protonated Guanine in RNA Structure and Function

Charged nucleobases have been found to occur in several known RNA molecules and are considered essential for their structure and function. The mechanism of their involvement is however not yet fully understood. Revelation of the role of N7-protonated guanine, in modulating the geometry and stability of noncanonical base pairs formed through its unprotonated edges [Watson-Crick (WC) and sugar], has triggered the need to evaluate the feasibility of similar roles of other protonated nucleobases [Halder et al., Phys Chem Chem Phys, 2015, 17, 26249]. In this context, N3 protonation of guanine makes an interesting case as its influence on the charge distribution of the WC edge is similar to that of N7 protonation, though its thermodynamic cost of protonation is significantly higher. In this work, we have carried out structural bioinformatics analyses and quantum mechanics-based calculations to show that N3 protonation of guanine may take place in a cellular environment, at least in the G:C W:W Trans and G:G W:H Cis base pairs. Our results provide a reasonable starting point for future investigations in order to address the larger mechanistic question.

Effect of the sulfur atom on S2 and S4 positions of the uracil ring in different DNA:RNA hybrid microhelixes with three nucleotide base pairs

The effect of the sulphur atom on the uracil ring was analyzed in different DNA:RNA microhelixes with three nucleotide base-pairs, including uridine, 2-thiouridine, 4-thiouridine, 2,4-dithiouridine, cytidine, adenosine and guanosine. Distinct backbone and helical parameters were optimized at different density functional (DFT) levels. The Watson-Crick pair with 2-thiouridine appears weaker than with uridine, but its interaction with water molecules appears easier. Two types of microhelixes were found, depending on the H-bond of H2' hydroxyl atom: A-type appears with the ribose ring in ³ E-envelope C_3' -endo, and B-type in ² E-envelope C_2' -endo. B-type is less common but it is more stable and with higher dipole-moment. The sulphur atoms significantly increase the dipole-moment of the microhelix, as well as the rise and propeller twist parameters. Simulations with four Na atoms H-bonded to the phosphate groups, and further hydration with explicit water molecules were carried out. A re-definition of the numerical value calculation of several base-pair and base-stacking parameters is suggested.

Trinucleotide Base Pair Stacking Free Energy for Understanding TF-DNA Recognition and the Functions of SNPs

Single nucleotide polymorphisms (SNPs) affect base pair stacking, which is the primary factor for maintaining the stability of DNA. However, the mechanism of how SNPs lead to phenotype variations is still unclear. In this work, we connected SNPs and base pair stacking by a 3-mer base pair stacking free energy matrix. The SNPs with large base pair stacking free energy differences led to phenotype variations. A molecular dynamics (MD) simulation was then applied. Our results showed that base pair stacking played an important role in the transcription factor (TF)-DNA interaction. Changes in DNA structure mainly originate from TF-DNA interactions, and with the increased base pair stacking free energy, the structure of DNA approaches its free type, although its binding affinity was increased by the SNP. In addition, quantitative models using base pair stacking features revealed that base pair stacking can be used to predict TF binding specificity. As such, our work combined knowledge from bioinformatics and structural biology and provided a new understanding of the relationship between SNPs and phenotype variations. The 3-mer base pair stacking free energy matrix is useful in high-throughput screening of SNPs and predicting TF-DNA binding affinity.

6'-Fluoro[4.3.0]bicyclo nucleic acid: synthesis, biophysical properties and molecular dynamics simulations

Here we report on the synthesis, biophysical properties and molecular modeling of oligonucleotides containing unsaturated 6'-fluoro[4.3.0]bicyclo nucleotides (6'F-bc4,3-DNA). Two 6'F-bc4,3 phosphoramidite building blocks (T and C) were synthesized starting from a previously described [3.3.0]bicyclic sugar. The conversion of this sugar to a gem-difluorinated tricyclic intermediate via difluorocarbene addition followed either by a NIS-mediated or Vorbrüggen nucleosidation yielded in both cases the β-tricyclic nucleoside as major anomer. Subsequent desilylation and cyclopropane ring opening of these tricyclic intermediates afforded the unsaturated 6'F-bc4,3 nucleosides. The successful incorporation of the corresponding phosphoramidite building blocks into oligonucleotides was achieved with tert-butyl hydroperoxide as oxidation agent. Thermal melting experiments of the modified duplexes disclosed a destabilizing effect versus DNA and RNA complements, but with a lesser degree of destabilization versus complementary DNA (ΔT m/mod = -1.5 to -3.7 °C). Molecular dynamics simulation on the nucleoside and oligonucleotide level revealed the preference of the C1'-exo/C2'-endo alignment of the furanose ring. Moreover, the simulation of duplexes with complementary RNA disclosed a DNA/RNA-type duplex structure suggesting that this modification might be a substrate for RNase H.

Selectivity and Promiscuity in TET-Mediated Oxidation of 5-Methylcytosine in DNA and RNA

Enzymes of the ten-eleven translocation (TET) family add diversity to the repertoire of nucleobase modifications by catalyzing the oxidation of 5-methylcytosine (5mC). TET enzymes were initially found to oxidize 5-methyl-2'-deoxycytidine in genomic DNA, yielding products that contribute to epigenetic regulation in mammalian cells, but have since been found to also oxidize 5-methylcytidine in RNA. Considering the different configurations of single-stranded (ss) and double-stranded (ds) DNA and RNA that coexist in a cell, defining the scope of TET's preferred activity and the mechanisms of substrate selectivity is critical to better understand the enzymes' biological functions. To this end, we have systematically examined the activity of human TET2 on DNA, RNA, and hybrid substrates in vitro. We found that, while ssDNA and ssRNA are well tolerated, TET2 is most proficient at dsDNA oxidation and discriminates strongly against dsRNA. Chimeric and hybrid substrates containing mixed DNA and RNA character helped reveal two main features by which the enzyme discriminates between substrates. First, the identity of the target nucleotide alone is the strongest reactivity determinant, with a preference for 5-methyldeoxycytidine, while both DNA or RNA are relatively tolerated on the rest of the target strand. Second, while a complementary strand is not required for activity, DNA is the preferred partner, and complementary RNA diminishes reactivity. Our biochemical analysis, complemented by molecular dynamics simulations, provides support for an active site optimally configured for dsDNA reactivity but permissive for various nucleic acid configurations, suggesting a broad range of plausible roles for TET-mediated 5mC oxidation in cells.

The AT Interstrand Cross-Link: Structure, Electronic Properties, and Influence on Charge Transfer in dsDNA

The interaction of chemical and physical agents with genetic material can lead to almost 80 different DNA damage formations. The targeted intentional DNA damage by radiotherapy or chemotherapy is a front-line anticancer therapy. An interstrand cross-link can result from ionization radiation or specific chemical agents, such as trans-/cisplatin activity. Here, the influence of the adenine and thymidine (AT) interstrand linkage, the covalent bond between the adenine N6 and thymidine C5 methylene group, on the isolated base pair as well as double-stranded DNA (dsDNA) was taken into quantum mechanical/molecular mechanical (QM/MM) consideration at the m062x/6-31+G*:UFF level of theory in the aqueous phase. All the results presented in this article, for the first time, show that an AT-interstrand cross-link (ICL) changes the positive and negative charge migration process due to a higher activation energy forced by the cross-link’s presence. However, the final radical cation destination in cross-linked DNA is left in the same place as in a native double-stranded-deoxyoligonucleotide. Additionally, the direction of the radical anion transfer was found to be opposite to that of native dsDNA. Therefore, it can be postulated that the appearance of the AT-ICL does not disturb the hole migration in the double helix, with subsequent effective changes in the electron migration process.

Interactive Visualization of RNA and DNA Structures

The analysis and visualization of nucleic acids (RNA and DNA) is playing an increasingly important role due to their fundamental importance for all forms of life and the growing number of known 3D structures of such molecules. The great complexity of these structures, in particular, those of RNA, demands interactive visualization to get deeper insights into the relationship between the 2D secondary structure motifs and their 3D tertiary structures. Over the last decades, a lot of research in molecular visualization has focused on the visual exploration of protein structures while nucleic acids have only been marginally addressed. In contrast to proteins, which are composed of amino acids, the ingredients of nucleic acids are nucleotides. They form structuring patterns that differ from those of proteins and, hence, also require different visualization and exploration techniques. In order to support interactive exploration of nucleic acids, the computation of secondary structure motifs as well as their visualization in 2D and 3D must be fast. Therefore, in this paper, we focus on the performance of both the computation and visualization of nucleic acid structure. We present a ray casting-based visualization of RNA and DNA secondary and tertiary structures, which enables for the first time real-time visualization of even large molecular dynamics trajectories. Furthermore, we provide a detailed description of all important aspects to visualize nucleic acid secondary and tertiary structures. With this, we close an important gap in molecular visualization.

Representation of DNA environment: Spiral staircase distribution function

In the present study, we investigated the local structure of DNA and its environment using a new visualization technique. The spiral staircase distribution function (SSDF) is determined as two-dimensional density distribution of atoms of water and ligands in local reference frames linked with each base pair of poly-DNA molecule, either GC or AT. This property of SSDF provides opportunity to study sequence-specific binding of ions, peptides, and other agents derived from a molecular dynamics computer simulation. The spatial structure of double-stranded DNA environment in water solution containing either Mg²⁺ or Na⁺ ions was investigated using of SSDF. The distributions of ions around GC and AT base pairs are shown separately. It is observed that Mg²⁺ ions interact with DNA atoms by means of the layer of water molecules and penetrate into the major groove only. Na⁺ ions have a direct contact with DNA atoms and penetrate both into the major and minor grooves of DNA. © 2018 Wiley Periodicals, Inc.

The temperature dependence of the helical twist of DNA

DNA is the carrier of all cellular genetic information and increasingly used in nanotechnology. Quantitative understanding and optimization of its functions requires precise experimental characterization and accurate modeling of DNA properties. A defining feature of DNA is its helicity. DNA unwinds with increasing temperature, even for temperatures well below the melting temperature. However, accurate quantitation of DNA unwinding under external forces and a microscopic understanding of the corresponding structural changes are currently lacking. Here we combine single-molecule magnetic tweezers measurements with atomistic molecular dynamics and coarse-grained simulations to obtain a comprehensive view of the temperature dependence of DNA twist. Experimentally, we find that DNA twist changes by ΔTw(T) = (-11.0 ± 1.2)°/(°C·kbp), independent of applied force, in the range of forces where torque-induced melting is negligible. Our atomistic simulations predict ΔTw(T) = (-11.1 ± 0.3)°/(°C·kbp), in quantitative agreement with experiments, and suggest that the untwisting of DNA with temperature is predominantly due to changes in DNA structure for defined backbone substates, while the effects of changes in substate populations are minor. Coarse-grained simulations using the oxDNA framework yield a value of ΔTw(T) = (-6.4 ± 0.2)°/(°C·kbp) in semi-quantitative agreement with experiments.

Analogues of P and Z as Efficient Artificially Expanded Genetic Information System

To artificially expand the genetic information system and to realize artificial life, it is necessary to discover new functional DNA bases that can form stable duplex DNA and participate in error-free replication. It is recently proposed that the 2-amino-imidazo[1,2- a]-1,3,5-triazin-4(8 H)one (P) and 6-amino-5-nitro-2(1 H)-pyridone (Z) would form a base pair complex, which is more stable than that of the normal G-C base pair and would produce an unperturbed duplex DNA. Here, by using quantum chemical calculations in aqueous medium, it is shown that the P and Z molecules can be modified with the help of electron-withdrawing and -donating substituents mainly found in B-DNA to generate new bases that can produce even more stable base pairs. Among the various bases studied, P3, P4, Z3, and Z5 are found to produce base pairs, which are about 2-15 kcal/mol more stable than the P-Z base pair. It is further shown that these base pairs can be stacked onto the G-C and A-T base pairs to produce stable dimers. The consecutive stacking of these base pairs is found to yield even more stable dimers. The influence of charge penetration effects and backbone atoms in stabilizing these dimers are also discussed. It is thus proposed that the P3, P4, Z3, and Z5 would form promiscuous artificial genetic information system and can be used for different biological applications. However, the evaluations of the dynamical effects of these bases in DNA-containing several nucleotides and the efficacy of DNA polymerases to replicate these bases would provide more insights.

Monte Carlo simulation algorithm for B-DNA

Understanding the structure-function relationship of biomolecules containing DNA has motivated experiments aimed at determining molecular structure using methods such as small-angle X-ray and neutron scattering (SAXS and SANS). SAXS and SANS are useful for determining macromolecular shape in solution, a process which benefits by using atomistic models that reproduce the scattering data. The variety of algorithms available for creating and modifying model DNA structures lack the ability to rapidly modify all-atom models to generate structure ensembles. This article describes a Monte Carlo algorithm for simulating DNA, not with the goal of predicting an equilibrium structure, but rather to generate an ensemble of plausible structures which can be filtered using experimental results to identify a sub-ensemble of conformations that reproduce the solution scattering of DNA macromolecules. The algorithm generates an ensemble of atomic structures through an iterative cycle in which B-DNA is represented using a wormlike bead-rod model, new configurations are generated by sampling bend and twist moves, then atomic detail is recovered by back mapping from the final coarse-grained configuration. Using this algorithm on commodity computing hardware, one can rapidly generate an ensemble of atomic level models, each model representing a physically realistic configuration that could be further studied using molecular dynamics. © 2016 Wiley Periodicals, Inc.

New Insights into the Role of DNA Shape on Its Recognition by p53 Proteins

The tumor suppressor p53 acts as a transcription factor recognizing diverse DNA response elements (REs). Previous structural studies of p53-DNA complexes revealed non-canonical Hoogsteen geometry of A/T base pairs at conserved CATG motifs leading to changes in DNA shape and its interface with p53. To study the effects of DNA shape on binding characteristics, we designed REs with modified base pairs "locked" into either Hoogsteen or Watson-Crick form. Here we present crystal structures of these complexes and their thermodynamic and kinetic parameters, demonstrating that complexes with Hoogsteen base pairs are stabilized relative to those with all-Watson-Crick base pairs. CATG motifs are abundant in p53REs such as GADD45 and p53R2 related to cell-cycle arrest and DNA repair. The high-resolution structures of these complexes validate their propensity to adopt the unique Hoogsteen-induced structure, thus providing insights into the functional role of DNA shape and broadening the mechanisms that contribute to DNA recognition by proteins.

Sub-picosecond proton tunnelling in deformed DNA hydrogen bonds under an asymmetric double-oscillator model

We present a model of proton tunnelling across DNA hydrogen bonds, compute the characteristic tunnelling time (CTT) from donor to acceptor and discuss its biological implications. The model is a double oscillator characterised by three geometry parameters describing planar deformations of the H bond, and a symmetry parameter representing the energy ratio between ground states in the individual oscillators. We discover that some values of the symmetry parameter lead to CTTs which are up to 40 orders of magnitude smaller than a previous model predicted. Indeed, if the symmetry parameter is sufficiently far from its extremal values of 1 or 0, then the proton's CTT under any physically realistic planar deformation is guaranteed to be below one picosecond, which is a biologically relevant time-scale. This supports theories of links between proton tunnelling and biological processes such as spontaneous mutation.

How a single 5-methylation of cytosine regulates the recognition of C/EBPβ transcription factor: a molecular dynamic simulation study

CpG methylation can regulate gene expression by altering the specific binding of protein and DNA. In order to understand how a single 5mC regulates protein-DNA interactions, we have compared the structures and dynamics of CEBP/βprotein-DNA complexes before and after methylation, and the results indicate that even a single 5mC can regulate protein-DNA recognition by steric-hindrance effect of methyl group and changing the hydrogen bond interactions. The interactions between the methyl group, mCpG motif, and the conserved residue arginine make the protein read out the variation of local environment, which further enhances the specific recognition and affects the base pair stacking. The stacking interactions can propagate along the backbone of DNA and lead to long-range allosteric effects, including obvious conformational variations for DNA base pairs.

Modeling DNA oxidation in water

A novel set of hole-site energies and electronic coupling parameters to be used, in the framework of the simplest tight-binding approximation, for predicting DNA hole trapping efficiencies and rates of hole transport in oxidized DNA is proposed. The novel parameters, significantly different from those previously reported in the literature, have been inferred from reliable density functional calculations, including both the sugar-phosphate ionic backbone and the effects of the aqueous environment. It is shown that most of the experimental oxidation free energies of DNA tracts and of oligonucleotides available from photoelectron spectroscopy and voltammetric measurements are reproduced with great accuracy, without the need for introducing sequence dependent parameters.

DEPENDENCE OF HOMO–LUMO GAP OF DNA BASE PAIR STEPS ON TWIST ANGLE: A DENSITY FUNCTIONAL APPROACH

Electronic structure calculations of all 10 unique base pair (bp) steps have been calculated to study the interaction energies of the bp steps, density of states (DOS), projected density of states (pDOS) using the density functional theory (DFT). Plane wave basis with ultrasoft pseudo-potential method has been used within the local density approximation (LDA) for the exchange correlation functional. Electron densities of the bp steps corresponding to HOMO and LUMO level have been calculated to understand the difference in stacking energies of the bp steps. The variation of HOMO–LUMO gap (g) of all possible bp steps on twist angle has been studied. We have observed that out of the 10 bp steps, 4 purine–purine bp steps (d(AA), d(GG), d(AG) and d(GA)), show significant variation of [Formula: see text] on twist angle. The observed variation on twist angle of d(AA) bp step has been explained by the calculated DOS and electron densities.

A unified approach for quantifying and interpreting DNA shape readout by transcription factors

Transcription factors (TFs) interpret DNA sequence by probing the chemical and structural properties of the nucleotide polymer. DNA shape is thought to enable a parsimonious representation of dependencies between nucleotide positions. Here, we propose a unified mathematical representation of the DNA sequence dependence of shape and TF binding, respectively, which simplifies and enhances analysis of shape readout. First, we demonstrate that linear models based on mononucleotide features alone account for 60-70% of the variance in minor groove width, roll, helix twist, and propeller twist. This explains why simple scoring matrices that ignore all dependencies between nucleotide positions can partially account for DNA shape readout by a TF Adding dinucleotide features as sequence-to-shape predictors to our model, we can almost perfectly explain the shape parameters. Building on this observation, we developed a post hoc analysis method that can be used to analyze any mechanism-agnostic protein-DNA binding model in terms of shape readout. Our insights provide an alternative strategy for using DNA shape information to enhance our understanding of how cis-regulatory codes are interpreted by the cellular machinery.

A novel form of RNA double helix based on G·U and C·A+ wobble base pairing

Wobble base pairs are critical in various physiological functions and have been linked to local structural perturbations in double-helical structures of nucleic acids. We report a 1.38-Å resolution crystal structure of an antiparallel octadecamer RNA double helix in overall A conformation, which includes a unique, central stretch of six consecutive wobble base pairs (W helix) with two G·U and four rare C·A+ wobble pairs. Four adenines within the W helix are N1-protonated and wobble-base-paired with the opposing cytosine through two regular hydrogen bonds. Combined with the two G·U pairs, the C·A+ base pairs facilitate formation of a half turn of W-helical RNA flanked by six regular Watson-Crick base pairs in standard A conformation on either side. RNA melting experiments monitored by differential scanning calorimetry, UV and circular dichroism spectroscopy demonstrate that the RNA octadecamer undergoes a pH-induced structural transition which is consistent with the presence of a duplex with C·A+ base pairs at acidic pH. Our crystal structure provides a first glimpse of an RNA double helix based entirely on wobble base pairs with possible applications in RNA or DNA nanotechnology and pH biosensors.

Expanding the repertoire of DNA shape features for genome-scale studies of transcription factor binding

Uncovering the mechanisms that affect the binding specificity of transcription factors (TFs) is critical for understanding the principles of gene regulation. Although sequence-based models have been used successfully to predict TF binding specificities, we found that including DNA shape information in these models improved their accuracy and interpretability. Previously, we developed a method for modeling DNA binding specificities based on DNA shape features extracted from Monte Carlo (MC) simulations. Prediction accuracies of our models, however, have not yet been compared to accuracies of models incorporating DNA shape information extracted from X-ray crystallography (XRC) data or Molecular Dynamics (MD) simulations. Here, we integrated DNA shape information extracted from MC or MD simulations and XRC data into predictive models of TF binding and compared their performance. Models that incorporated structural information consistently showed improved performance over sequence-based models regardless of data source. Furthermore, we derived and validated nine additional DNA shape features beyond our original set of four features. The expanded repertoire of 13 distinct DNA shape features, including six intra-base pair and six inter-base pair parameters and minor groove width, is available in our R/Bioconductor package DNAshapeR and enables a comprehensive structural description of the double helix on a genome-wide scale.

Consequences of Mg2+ binding on the geometry and stability of RNA base pairs

Quantum chemical calculations reveal the role of magnesium in stabilizing the geometries of intrinsically unstable RNA base pairs.

A DFT study of 2-aminopurine-containing dinucleotides: prediction of stacked conformations with B-DNA structure

The fluorescence properties of dinucleotides incorporating 2-aminopurine (2AP) suggest that the simplest oligonucleotides adopt conformations similar to those found in duplex DNA. However, there is a lack of structural data for these systems. We report a density functional theory (DFT) study of the structures of 2AP-containing dinucleotides (deoxydinucleoside monophosphates), including full geometry optimisation of the sugar-phosphate backbone. Our DFT calculations employ the M06-2X functional for reliable treatment of dispersion interactions and include implicit aqueous solvation. Dinucleotides with 2AP in the 5'-position and each of the natural bases in the 3'-position are examined, together with the analogous 5'-adenine-containing systems. Computed structures are compared in detail with typical B-DNA base-step parameters, backbone torsional angles and sugar pucker, derived from crystallographic data. We find that 2AP-containing dinucleotides adopt structures that closely conform to B-DNA in all characteristic parameters. The structures of 2AP-containing dinucleotides closely resemble those of their adenine-containing counterparts, demonstrating the fidelity of 2AP as a mimic of the natural base. As a first step towards exploring the conformational heterogeneity of dinucleotides, we also characterise an imperfectly stacked conformation and one in which the bases are completely unstacked.

Parallels between DNA and collagen - comparing elastic models of the double and triple helix

Multi-stranded helices are widespread in nature. The interplay of polymeric properties with biological function is seldom discussed. This study probes analogies between structural and mechanical properties of collagen and DNA. We modeled collagen with Eulerian rotational and translational parameters of adjacent rungs in the triple-helix ladder and developed statistical potentials by extracting the dispersion of the parameters from a database of atomic-resolution structures. The resulting elastic model provides a common quantitative way to describe collagen deformations upon interacting with integrins or matrix metalloproteinase and DNA deformations upon protein binding. On a larger scale, deformations in Type I collagen vary with a periodicity consistent with the D-periodic banding of higher-order fibers assemblies. This indicates that morphologies of natural higher-order collagen packing might be rooted in the characteristic deformation patterns.

Electronic Structures of LNA Phosphorothioate Oligonucleotides

Important oligonucleotides in anti-sense research have been investigated in silico and experimentally. This involves quantum mechanical (QM) calculations and chromatography experiments on locked nucleic acid (LNA) phosphorothioate (PS) oligonucleotides. iso-potential electrostatic surfaces are essential in this study and have been calculated from the wave functions derived from the QM calculations that provide binding information and other properties of these molecules. The QM calculations give details of the electronic structures in terms of e.g., energy and bonding, which make them distinguish or differentiate between the individual PS diastereoisomers determined by the position of sulfur atoms. Rules are derived from the electronic calculations of these molecules and include the effects of the phosphorothioate chirality and formation of electrostatic potential surfaces. Physical and electrochemical descriptors of the PS oligonucleotides are compared to the experiments in which chiral states on these molecules can be distinguished. The calculations demonstrate that electronic structure, electrostatic potential, and topology are highly sensitive to single PS configuration changes and can give a lead to understanding the activity of the molecules.

LNA effects on DNA binding and conformation: from single strand to duplex and triplex structures

The anti-gene strategy is based on sequence-specific recognition of double-strand DNA by triplex forming (TFOs) or DNA strand invading oligonucleotides to modulate gene expression. To be efficient, the oligonucleotides (ONs) should target DNA selectively, with high affinity. Here we combined hybridization analysis and electrophoretic mobility shift assay with molecular dynamics (MD) simulations to better understand the underlying structural features of modified ONs in stabilizing duplex- and triplex structures. Particularly, we investigated the role played by the position and number of locked nucleic acid (LNA) substitutions in the ON when targeting a c-MYC or FXN (Frataxin) sequence. We found that LNA-containing single strand TFOs are conformationally pre-organized for major groove binding. Reduced content of LNA at consecutive positions at the 3'-end of a TFO destabilizes the triplex structure, whereas the presence of Twisted Intercalating Nucleic Acid (TINA) at the 3'-end of the TFO increases the rate and extent of triplex formation. A triplex-specific intercalating benzoquinoquinoxaline (BQQ) compound highly stabilizes LNA-containing triplex structures. Moreover, LNA-substitution in the duplex pyrimidine strand alters the double helix structure, affecting x-displacement, slide and twist favoring triplex formation through enhanced TFO major groove accommodation. Collectively, these findings should facilitate the design of potent anti-gene ONs.

Analyzing DNA curvature and its impact on the ionic environment: application to molecular dynamics simulations of minicircles

We propose a method for analyzing the magnitude and direction of curvature within nucleic acids, based on the curvilinear helical axis calculated by Curves+. The method is applied to analyzing curvature within minicircles constructed with varying degrees of over- or under-twisting. Using the molecular dynamics trajectories of three different minicircles, we are able to quantify how curvature varies locally both in space and in time. We also analyze how curvature influences the local environment of the minicircles, notably via increased heterogeneity in the ionic distributions surrounding the double helix. The approach we propose has been integrated into Curves+ and the utilities Canal (time trajectory analysis) and Canion (environmental analysis) and can be used to study a wide variety of static and dynamic structural data on nucleic acids.

DSSR-enhanced visualization of nucleic acid structures in Jmol

Sophisticated and interactive visualizations are essential for making sense of the intricate 3D structures of macromolecules. For proteins, secondary structural components are routinely featured in molecular graphics visualizations. However, the field of RNA structural bioinformatics is still lagging behind; for example, current molecular graphics tools lack built-in support even for base pairs, double helices, or hairpin loops. DSSR (Dissecting the Spatial Structure of RNA) is an integrated and automated command-line tool for the analysis and annotation of RNA tertiary structures. It calculates a comprehensive and unique set of features for characterizing RNA, as well as DNA structures. Jmol is a widely used, open-source Java viewer for 3D structures, with a powerful scripting language. JSmol, its reincarnation based on native JavaScript, has a predominant position in the post Java-applet era for web-based visualization of molecular structures. The DSSR-Jmol integration presented here makes salient features of DSSR readily accessible, either via the Java-based Jmol application itself, or its HTML5-based equivalent, JSmol. The DSSR web service accepts 3D coordinate files (in mmCIF or PDB format) initiated from a Jmol or JSmol session and returns DSSR-derived structural features in JSON format. This seamless combination of DSSR and Jmol/JSmol brings the molecular graphics of 3D RNA structures to a similar level as that for proteins, and enables a much deeper analysis of structural characteristics. It fills a gap in RNA structural bioinformatics, and is freely accessible (via the Jmol application or the JSmol-based website http://jmol.x3dna.org).

Martini Coarse-Grained Force Field: Extension to RNA

RNA has an important role not only as the messenger of genetic information but also as a regulator of gene expression. Given its central role in cell biology, there is significant interest in studying the structural and dynamic behavior of RNA in relation to other biomolecules. Coarse-grain molecular dynamics simulations are a key tool to that end. Here, we have extended the coarse-grain Martini force field to include RNA after our recent extension to DNA. In the same way DNA was modeled, the tertiary structure of RNA is constrained using an elastic network. This model, therefore, is not designed for applications involving RNA folding but rather offers a stable RNA structure for studying RNA interactions with other (bio)molecules. The RNA model is compatible with all other Martini models and opens the way to large-scale explicit-solvent molecular dynamics simulations of complex systems involving RNA.

[KCl] Dependence of B-DNA Groove Bending Anisotropy

The energetics of B-DNA bending toward the major and minor grooves were quantified by free energy simulations at four different KCl concentrations. Increased [KCl] led to more flexible DNA, with persistence lengths that agreed well with experimental values. At all salt concentrations, major groove bending was preferred, although preferences for major and minor groove bending were similar for the A-tract containing sequence. Since the phosphate repulsions and DNA internal energy favored minor groove bending, the preference for major groove bending was thought to originate from differences in solvation. Water in the minor groove was tighter bound than water in the major groove, and harder to displace than major groove water, which favored the compression of the major groove upon bending. Higher [KCl] decreased the persistence length for both major and minor groove bending but did not greatly affect the free energy spacing between the minor and major groove bending curves. For sequences without A-tracts, salt affected major and minor bending to nearly the same degree, and did not change the preference for major groove bending. For the A-tract containing sequence, an increase in salt concentration decreased the already small energetic difference between major and minor groove bending. Since salts did not significantly affect the relative differences in bending energetics and hydration, it is likely that the increased bending flexibilities upon salt increase are simply due to screening.

Role of Pseudoisocytidine Tautomerization in Triplex-Forming Oligonucleotides: In Silico and in Vitro Studies

Pseudoisocytidine (^ΨC) is a synthetic cytidine analogue that can target DNA duplex to form parallel triplex at neutral pH. Pseudoisocytidine has mainly two tautomers, of which only one is favorable for triplex formation. In this study, we investigated the effect of sequence on ^ΨC tautomerization using λ-dynamics simulation, which takes into account transitions between states. We also performed in vitro binding experiments with sequences containing ^ΨC and furthermore characterized the structure of the formed triplex using molecular dynamics simulation. We found that the neighboring methylated or protonated cytidine promotes the formation of the favorable tautomer, whereas the neighboring thymine or locked nucleic acid has a poor effect, and consecutive ^ΨC has a negative influence. The deleterious effect of consecutive ^ΨC in a triplex formation was confirmed using in vitro binding experiments. Our findings contribute to improving the design of ^ΨC-containing triplex-forming oligonucleotides directed to target G-rich DNA sequences.

How Does Mg2+ Modulate the RNA Folding Mechanism: A Case Study of the G:C W:W Trans Basepair

Reverse Watson-Crick G:C basepairs (G:C W:W Trans) occur frequently in different functional RNAs. This is one of the few basepairs whose gas-phase-optimized isolated geometry is inconsistent with the corresponding experimental geometry. Several earlier studies indicate that through post-transcriptional modification, direct protonation, or coordination with Mg²⁺, accumulation of positive charge near N7 of guanine can stabilize the experimental geometry. Interestingly, recent studies reveal significant variation in the position of putatively bound Mg²⁺. This, in conjunction with recently raised doubts regarding some of the Mg²⁺ assignments near the imino nitrogen of guanine, is suggestive of the existence of multiple Mg²⁺ binding modes for this basepair. Our detailed investigation of Mg²⁺-bound G:C W:W Trans pairs occurring in high-resolution RNA crystal structures shows that they are found in 14 different contexts, eight of which display Mg²⁺ binding at the Hoogsteen edge of guanine. Further examination of occurrences in these eight contexts led to the characterization of three different Mg²⁺ binding modes: 1) direct binding via N7 coordination, 2) direct binding via O6 coordination, and 3) binding via hydrogen-bonding interaction with the first-shell water molecules. In the crystal structures, the latter two modes are associated with a buckled and propeller-twisted geometry of the basepair. Interestingly, respective optimized geometries of these different Mg²⁺ binding modes (optimized using six different DFT functionals) are consistent with their corresponding experimental geometries. Subsequent interaction energy calculations at the MP2 level, and decomposition of its components, suggest that for G:C W:W Trans , Mg²⁺ binding can fine tune the basepair geometries without compromising with their stability. Our results, therefore, underline the importance of the mode of binding of Mg²⁺ ions in shaping RNA structure, folding and function.

Crystal structure of d(CCGGGGTACCCCGG)2 at 1.4 Å resolution

The X-ray crystal structure of the DNA tetradecamer sequence d(CCGGGGTACCCCGG)2 is reported at 1.4 Å resolution in the tetragonal space group P41212. The sequence was designed to fold as a four-way junction. However, it forms an A-type double helix in the presence of barium chloride. The metal ion could not be identified in the electron-density map. The crystallographic asymmetric unit consists of one A-type double helix with 12 base pairs per turn, in contrast to 11 base pairs per turn for canonical A-DNA. A large number of solvent molecules have been identified in both the grooves of the duplex and around the backbone phosphate groups.

Consecutive non-natural PZ nucleobase pairs in DNA impact helical structure as seen in 50 μs molecular dynamics simulations

Z

Little is known about the influence of multiple consecutive 'non-standard' ( , 6-amino-5-nitro-2(1H)-pyridone, and , 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one) nucleobase pairs on the structural parameters of duplex DNA. nucleobase pairs follow standard rules for Watson-Crick base pairing but have rearranged hydrogen bonding donor and acceptor groups. Using the X-ray crystal structure as a starting point, we have modeled the motions of a DNA duplex built from a self-complementary oligonucleotide (5΄-CTTATPPPZZZATAAG-3΄) in water over a period of 50 μs and calculated DNA local parameters, step parameters, helix parameters, and major/minor groove widths to examine how the presence of multiple, consecutive nucleobase pairs might impact helical structure. In these simulations, the -containing DNA duplex exhibits a significantly wider major groove and greater average values of stagger, slide, rise, twist and h-rise than observed for a 'control' oligonucleotide in which nucleobase pairs are replaced by . The molecular origins of these structural changes are likely associated with at least two differences between and . First, the electrostatic properties of differ from in terms of density distribution and dipole moment. Second, differences are seen in the base stacking of pairs in dinucleotide steps, arising from energetically favorable stacking of the nitro group in with π-electrons of the adjacent base.

Comparison of X-ray crystal structures of a tetradecamer sequence d(CCCGGGTACCCGGG)2 at 1.7 Å resolution

We present here a comparison of three different X-ray crystal structures of DNA tetradecamer sequence d(CCCGGGTACCCGGG)₂ all at about 1.7 Å resolution. The sequence was designed as an attempt to form a DNA four-way junction with A-type helical arms. However, in the presence of zinc, magnesium, and in the absence of any metal ion, it does not take up the junction structure, but forms an A-type double helix. This allowed us to study possible conformational changes in the double helix due to the presence of metal ions. Upon addition of the zinc ion, there is a change in the space group from P4₁2₁2 to P4₁. The overall conformation of the duplex remains the same. There are small changes in the interaction of the metal ions with the DNA. In the zinc-bound structure, there are two zinc ions that show direct interaction with the N7 atoms of terminal G₁₃ bases at either end of the molecule. There are small changes in the interhelical contacts. The consequence of these differences is to break some of the symmetry and change the space group.

The tilt-dependent potential of mean force of a pair of DNA oligomers from all-atom molecular dynamics simulations

Electrostatic interactions between DNA molecules have been extensively studied experimentally and theoretically, but several aspects (e.g. its role in determining the pitch of the cholesteric DNA phase) still remain unclear. Here, we performed large-scale all-atom molecular dynamics simulations in explicit water and 150 mM sodium chloride, to reconstruct the potential of mean force (PMF) of two DNA oligomers 24 base pairs long as a function of their interaxial angle and intermolecular distance. We find that the potential of mean force is dominated by total DNA charge, and not by the helical geometry of its charged groups. The theory of homogeneously charged cylinders fits well all our simulation data, and the fit yields the optimal value of the total compensated charge on DNA to  ≈65% of its total fixed charge (arising from the phosphorous atoms), close to the value expected from Manning's theory of ion condensation. The PMF calculated from our simulations does not show a significant dependence on the handedness of the angle between the two DNA molecules, or its size is on the order of [Formula: see text]. Thermal noise for molecules of the studied length seems to mask the effect of detailed helical charge patterns of DNA. The fact that in monovalent salt the effective interaction between two DNA molecules is independent on the handedness of the tilt may suggest that alternative mechanisms are required to understand the cholesteric phase of DNA.