Critical effect of the N2 amino group on structure, dynamics, and elasticity of DNA polypurine tracts

Structural basis of DNA recognition by BEN domain proteins reveals a role for oligomerization in unmethylated DNA selection by BANP

The BEN domain is a newly discovered type of DNA-binding domain that exists in a variety of species. There are nine BEN domain-containing proteins in humans, and most have been shown to have chromatin-related functions. NACC1 preferentially binds to CATG motif-containing sequences and functions primarily as a transcriptional coregulator. BANP and BEND3 preferentially bind DNA bearing unmethylated CpG motifs, and they function as CpG island-binding proteins. To date, the DNA recognition mechanism of quite a few of these proteins remains to be determined. In this study, we solved the crystal structures of the BEN domains of NACC1 and BANP in complex with their cognate DNA substrates. We revealed the details of DNA binding by these BEN domain proteins and unexpectedly revealed that oligomerization is required for BANP to select unmethylated CGCG motif-containing DNA substrates. Our study clarifies the controversies surrounding DNA recognition by BANP and demonstrates a new mechanism by which BANP selects unmethylated CpG motifs and functions as a CpG island-binding protein. This understanding will facilitate further exploration of the physiological functions of the BEN domain proteins in the future.

Temperature-dependent elasticity of DNA, RNA, and hybrid double helices

Nucleic acid double helices in their DNA, RNA, and DNA-RNA hybrid form play a fundamental role in biology and are main building blocks of artificial nanostructures, but how their properties depend on temperature remains poorly understood. Here, we report thermal dependence of dynamic bending persistence length, twist rigidity, stretch modulus, and twist-stretch coupling for DNA, RNA, and hybrid duplexes between 7°C and 47°C. The results are based on all-atom molecular dynamics simulations using different force field parameterizations. We first demonstrate that unrestrained molecular dynamics can reproduce experimentally known mechanical properties of the duplexes at room temperature. Beyond experimentally known features, we also infer the twist rigidity and twist-stretch coupling of the hybrid duplex. As for the temperature dependence, we found that increasing temperature softens all the duplexes with respect to bending, twisting, and stretching. The relative decrease of the stretch moduli is 0.003-0.004/°C, similar for all the duplex variants despite their very different stretching stiffness, whereas RNA twist stiffness decreases by 0.003/°C, and smaller values are found for the other elastic moduli. The twist-stretch couplings are nearly unaffected by temperature. The stretching, bending, and twisting stiffness all include an important entropic component. Relation of our results to the two-state model of DNA flexibility is discussed. Our work provides temperature-dependent elasticity of nucleic acid duplexes at the microsecond scale relevant for initial stages of protein binding.

Harnessing non-Watson-Crick's base pairing to enhance CRISPR effectors cleavage activities and enable gene editing in mammalian cells

Genomic DNA of the cyanophage S-2L virus is composed of 2-aminoadenine (Z), thymine (T), guanine (G), and cytosine (C), forming the genetic alphabet ZTGC, which violates Watson-Crick base pairing rules. The Z-base has an extra amino group on the two position that allows the formation of a third hydrogen bond with thymine in DNA strands. Here, we explored and expanded applications of this non-Watson-Crick base pairing in protein expression and gene editing. Both ZTGC-DNA (Z-DNA) and ZUGC-RNA (Z-RNA) produced in vitro show detectable compatibility and can be decoded in mammalian cells, including Homo sapiens cells. Z-crRNA can guide CRISPR-effectors SpCas9 and LbCas12a to cleave specific DNA through non-Watson-Crick base pairing and boost cleavage activities compared to A-crRNA. Z-crRNA can also allow for efficient gene and base editing in human cells. Together, our results help pave the way for potential strategies for optimizing DNA or RNA payloads for gene editing therapeutics and give insights to understanding the natural Z-DNA genome.

Self-Consistent Parameterization of DNA Residues for the Non-Polarizable AMBER Force Fields

Fixed-charge (non-polarizable) forcefields are accurate and computationally efficient tools for modeling the molecular dynamics of nucleic acid polymers, particularly DNA, well into the µs timescale. The continued utility of these forcefields depends in part on expanding the residue set in step with advancing nucleic acid chemistry and biology. A key step in parameterizing new residues is charge derivation which is self-consistent with the existing residues. As atomic charges are derived by fitting against molecular electrostatic potentials, appropriate structural models are critical. Benchmarking against the existing charge set used in current AMBER nucleic acid forcefields, we report that quantum mechanical models of deoxynucleosides, even at a high level of theory, are not optimal structures for charge derivation. Instead, structures from molecular mechanics minimization yield charges with up to 6-fold lower RMS deviation from the published values, due to the choice of such an approach in the derivation of the original charge set. We present a contemporary protocol for rendering self-consistent charges as well as optimized charges for a panel of nine non-canonical residues that will permit comparison with literature as well as studying the dynamics of novel DNA polymers.

DNA methylation cues in nucleosome geometry, stability and unwrapping

Cytosine methylation at the 5-carbon position is an essential DNA epigenetic mark in many eukaryotic organisms. Although countless structural and functional studies of cytosine methylation have been reported, our understanding of how it influences the nucleosome assembly, structure, and dynamics remains obscure. Here, we investigate the effects of cytosine methylation at CpG sites on nucleosome dynamics and stability. By applying long molecular dynamics simulations on several microsecond time scale, we generate extensive atomistic conformational ensembles of full nucleosomes. Our results reveal that methylation induces pronounced changes in geometry for both linker and nucleosomal DNA, leading to a more curved, under-twisted DNA, narrowing the adjacent minor grooves, and shifting the population equilibrium of sugar-phosphate backbone geometry. These DNA conformational changes are associated with a considerable enhancement of interactions between methylated DNA and the histone octamer, doubling the number of contacts at some key arginines. H2A and H3 tails play important roles in these interactions, especially for DNA methylated nucleosomes. This, in turn, prevents a spontaneous DNA unwrapping of 3-4 helical turns for the methylated nucleosome with truncated histone tails, otherwise observed in the unmethylated system on several microseconds time scale.

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

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

Interaction of Thymine DNA Glycosylase with Oxidised 5-Methyl-cytosines in Their Amino- and Imino-Forms

Thymine DNA Glycosylase (TDG) is an enzyme of the base excision repair mechanism and removes damaged or mispaired bases from DNA via hydrolysis of the glycosidic bond. Specificity is of high importance for such a glycosylase, so as to avoid the damage of intact DNA. Among the substrates reported for TDG are mispaired uracil and thymine but also formyl-cytosine and carboxyl-cytosine. Methyl-cytosine and hydroxylmethyl-cytosine are, in contrast, not processed by the TDG enzyme. We have in this work employed molecular dynamics simulations to explore the conformational dynamics of DNA carrying a formyl-cytosine or carboxyl-cytosine and compared those to DNA with the non-cognate bases methyl-cytosine and hydroxylmethyl-cytosine, as amino and imino tautomers. Whereas for the mispairs a wobble conformation is likely decisive for recognition, all amino tautomers of formyl-cytosine and carboxyl-cytosine exhibit the same Watson–Crick conformation, but all imino tautomers indeed form wobble pairs. The conformational dynamics of the amino tautomers in free DNA do not exhibit differences that could be exploited for recognition, and also complexation to the TDG enzyme does not induce any alteration that would indicate preferable binding to one or the other oxidised methyl-cytosine. The imino tautomers, in contrast, undergo a shift in the equilibrium between a closed and a more open, partially flipped state, towards the more open form upon complexation to the TDG enzyme. This stabilisation of the more open conformation is most pronounced for the non-cognate bases methyl-cytosine and hydroxyl-cytosine and is thus not a likely mode for recognition. Moreover, calculated binding affinities for the different forms indicate the imino forms to be less likely in the complexed DNA. These findings, together with the low probability of imino tautomers in free DNA and the indifference of the complexed amino tautomers, suggest that discrimination of the oxidised methyl-cytosines does not take place in the initial complex formation.

Molecular mechanism of methyl-dependent and spatial-specific DNA recognition of c-Jun homodimer

DNA methylation is important in regulation of gene expression and normal development because it alters the interplay between protein and DNA. Experiments have shown that a single 5-methylcytosine at different CpG sites (mCpG) might have different effects on specific recognition, but the atomistic origin and dynamic details are largely unclear. In this work, we investigated the mechanism of monomethylation at different CpG sites in the cognate motif and the cooperativity of full methylation. By constructing four models of c-Jun/Jun protein binding to the 5[Formula: see text]-XGAGTCA-3[Formula: see text] (X represents C or methylated C) motif, we characterized the dynamics of the contact interface using the all-atom molecular dynamics method. Free energy analysis of MM/GBSA suggests that regardless of whether the C12pG13 site of the bottom strand is methylated, the effects from mC25 of the top strand are dominant and can moderately enhance the binding by [Formula: see text] 31 kcal/mol, whereas mC12 showed a relatively small contribution, in agreement with the experimental data. Remarkably, we found that this spatial-specific influence was induced by different regulatory rules. The influence of the mC25 site is mainly mediated by steric hindrance. The additional methyl group leads to the conformational changes in nearby residues and triggers an obvious structural bending in the protein, which results in the formation of a new T-Asn-C triad that enhances the specific recognition of TCA half-sites. The substitution of the methyl group at the mC12 site of the bottom strand breaks the original H-bonds directly. Such changes in electrostatic interactions also lead to the remote allosteric effects of protein by multifaceted interactions but have negligible contributions to binding. Although these two influence modes are different, they can both fine-tune the local environment, which might produce remote allosteric effects through protein-protein interactions. Further analysis reveals that the discrepancies in these two modes are primarily due to their location. Moreover, when both sites are methylated, the major determinant of binding specificity depends on the context and the location of the methylation site, which is the result of crosstalk and cooperativity.

The 2-Amino Group of 8-Aza-7-deaza-7-bromopurine-2,6-diamine and Purine-2,6-diamine as Stabilizer for the Adenine-Thymine Base Pair in Heterochiral DNA with Strands in Anomeric Configuration

Stabilization of DNA is beneficial for many applications in the fields of DNA therapeutics, diagnostics, and materials science. Now, this phenomenon is studied on heterochiral DNA, an autonomous DNA recognition system with complementary strands in α-D and β-D configuration showing parallel strand orientation. The 12-mer heterochiral duplexes were constructed from anomeric (α/β-D) oligonucleotide single-strands. Purine-2,6-diamine and 8-aza-7-deaza-7-bromopurine-2,6-diamine 2'-deoxyribonucleosides having the capability to form tridentate base pairs with dT were used to strengthen the stability of the dA-dT base pair. Tm data and thermodynamic values obtained from UV melting profiles indicated that the 8-aza-7-deaza 2'-deoxyribonucleoside decorated with a bromo substituent is so far the most efficient stabilizer for heterochiral DNA. Compared with that, the stabilizing effect of the purine-2,6-diamine 2'-deoxyribonucleoside is low. Global changes of helix structures were identified by circular dichroism (CD) spectra during melting.

Structural dynamics of double-stranded DNA with epigenome modification

Modification of cytosine plays an important role in epigenetic regulation of gene expression and genome stability. Cytosine is converted to 5-methylcytosine (5mC) by DNA methyltransferase; in turn, 5mC may be oxidized to 5-hydroxymethylcytosine (5hmC) by ten-eleven translocation enzyme. The structural flexibility of DNA is known to affect the binding of proteins to methylated DNA. Here, we have carried out a semi-quantitative analysis of the dynamics of double-stranded DNA (dsDNA) containing various epigenetic modifications by combining data from imino 1H exchange and imino 1H R1ρ relaxation dispersion NMR experiments in a complementary way. Using this approach, we characterized the base-opening (kopen) and base-closing (kclose) rates, facilitating a comparison of the base-opening and -closing process of dsDNA containing cytosine in different states of epigenetic modification. A particularly striking result is the increase in the kopen rate of hemi-methylated dsDNA 5mC/C relative to unmodified or fully methylated dsDNA, indicating that the Watson-Crick base pairs undergo selective destabilization in 5mC/C. Collectively, our findings imply that the epigenetic modulation of cytosine dynamics in dsDNA mediates destabilization of the GC Watson-Crick base pair to allow base-flipping in living cells.

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

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

Systematic analysis of the binding behaviour of UHRF1 towards different methyl- and carboxylcytosine modification patterns at CpG dyads

The multi-domain protein UHRF1 is essential for DNA methylation maintenance and binds DNA via a base-flipping mechanism with a preference for hemi-methylated CpG sites. We investigated its binding to hemi- and symmetrically modified DNA containing either 5-methylcytosine (mC), 5-hydroxymethylcytosine (hmC), 5-formylcytosine (fC), or 5-carboxylcytosine (caC). Our experimental results indicate that UHRF1 binds symmetrically carboxylated and hybrid methylated/carboxylated CpG dyads in addition to its previously reported substrates. Complementary molecular dynamics simulations provide a possible mechanistic explanation of how the protein could differentiate between modification patterns. First, we observe different local binding modes in the nucleotide binding pocket as well as the protein's NKR finger. Second, both DNA modification sites are coupled through key residues within the NKR finger, suggesting a communication pathway affecting protein-DNA binding for carboxylcytosine modifications. Our results suggest a possible additional function of the hemi-methylation reader UHRF1 through binding of carboxylated CpG sites. This opens the possibility of new biological roles of UHRF1 beyond DNA methylation maintenance and of oxidised methylcytosine derivates in epigenetic regulation.

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.

Nanomechanics of Diaminopurine-Substituted DNA

2,6-diaminopurine (DAP) is a nucleobase analog of adenine. When incorporated into double-stranded DNA (dsDNA), it forms three hydrogen bonds with thymine. Rare in nature, DAP substitution alters the physical characteristics of a DNA molecule without sacrificing sequence specificity. Here, we show that in addition to stabilizing double-strand hybridization, DAP substitution also changes the mechanical and conformational properties of dsDNA. Thermal melting experiments reveal that DAP substitution raises melting temperatures without diminishing sequence-dependent effects. Using a combination of atomic force microscopy (AFM), magnetic tweezer (MT) nanomechanical assays, and circular dichroism spectroscopy, we demonstrate that DAP substitution increases the flexural rigidity of dsDNA yet also facilitates conformational shifts, which manifest as changes in molecule length. DAP substitution increases both the static and dynamic persistence length of DNA (measured by AFM and MT, respectively). In the static case (AFM), in which tension is not applied to the molecule, the contour length of DAP-DNA appears shorter than wild-type (WT)-DNA; under tension (MT), they have similar dynamic contour lengths. At tensions above 60 pN, WT-DNA undergoes characteristic overstretching because of strand separation (tension-induced melting) and spontaneous adoption of a conformation termed S-DNA. Cyclic overstretching and relaxation of WT-DNA at near-zero loading rates typically yields hysteresis, indicative of tension-induced melting; conversely, cyclic stretching of DAP-DNA showed little or no hysteresis, consistent with the adoption of the S-form, similar to what has been reported for GC-rich sequences. However, DAP-DNA overstretching is distinct from GC-rich overstretching in that it happens at a significantly lower tension. In physiological salt conditions, evenly mixed AT/GC DNA typically overstretches around 60 pN. GC-rich sequences overstretch at similar if not slightly higher tensions. Here, we show that DAP-DNA overstretches at 52 pN. In summary, DAP substitution decreases the overall stability of the B-form double helix, biasing toward non-B-form DNA helix conformations at zero tension and facilitating the B-to-S transition at high tension.

Dispersed dynamics of solvation in G-quadruplex DNA: comparison of dynamic Stokes shifts of probes in parallel and antiparallel quadruplex structures

G-quadruplex DNA (GqDNA) structures play an important role in many specific cellular functions and are promising anti-tumor targets for small molecules (ligands). Here, we measured the dynamic Stokes shift of a ligand (Hoechst) bound to parallel c-Myc (mPu22) GqDNA over five decades of time from 100 fs to 10 ns, and compared it with the previously reported dynamics of DAPI bound to antiparallel human telomeric (hTelo22) GqDNA (Pal et al 2015 J. Phys. Chem. Lett. 6 1754). Stokes shift data from fluorescence up-conversion and time-correlated single photon counting experiments was combined to cover the broad dynamic range. The results show that the solvation dynamics of Hoechst in parallel mPu22 GqDNA follow a power law relaxation, added to fast 2 ps exponential relaxation, from 100 fs to 10 ns, with only a subtle difference of power law exponents in the two ligand-GqDNA systems (0.06 in Hoechst-mPu22 compared to 0.16 in DAPI-hTelo22). We measured steady-state fluorescence spectra and time-resolved anisotropy decays which confirm the tight binding of Hoechst to parallel mPu22 with a binding constant of ~1  ×  10⁵ M^-1. The molecular docking of Hoechst in parallel GqDNA followed by a 50 ns molecular dynamics (MD) simulation on a Hoechst-GqDNA complex reveals that Hoechst binds to one of the outer G-tetrads by end-stacking near G13 and G4, which is different from the binding site of DAPI inside a groove of antiparallel hTelo22 GqDNA. Reconciling previous experimental and simulation results, we assign the 2 ps component to the hydration dynamics of only weakly perturbed water near mPu22 and the power law relaxation to the coupled motion of water and DNA (i.e. DNA backbone, unpaired bases and loops connecting G-tetrads) which come near the Hoechst inside parallel GqDNA.

Effect of C5-Methylation of Cytosine on the UV-Induced Reactivity of Duplex DNA: Conformational and Electronic Factors

C5-methylation of cytosines is strongly correlated with UV-induced mutations detected in skin cancers. Mutational hot-spots appearing at TCG sites are due to the formation of pyrimidine cyclobutane dimers (CPDs). The present study, performed for the model DNA duplex (TCGTA)3·(TACGA)3 and the constitutive single strands, examines the factors underlying the effect of C5-methylation on pyrimidine dimerization at TCG sites. This effect is quantified for the first time by quantum yields ϕ. They were determined following irradiation at 255, 267, and 282 nm and subsequent photoproduct analysis using HPLC coupled to mass spectrometry. C5-methylation leads to an increase of the CPD quantum yield up to 80% with concomitant decrease of that of pyrimidine(6-4) pyrimidone adducts (64PPs) by at least a factor of 3. The obtained ϕ values cannot be explained only by the change of the cytosine absorption spectrum upon C5-methylation. The conformational and electronic factors that may affect the dimerization reaction are discussed in light of results obtained by fluorescence spectroscopy, molecular dynamics simulations, and quantum mechanical calculations. Thus, it appears that the presence of an extra methyl on cytosine affects the sugar puckering, thereby enhancing conformations of the TC step that are prone to CPD formation but less favorable to 64PPs. In addition, C5-methylation diminishes the amplitude of conformational motions in duplexes; in the resulting stiffer structure, ππ* excitations may be transferred from initially populated exciton states to reactive pyrimidines giving rise to CPDs.

Mechanical properties of base-modified DNA are not strictly determined by base stacking or electrostatic interactions

This work probes the mystery of what balance of forces creates the extraordinary mechanical stiffness of DNA to bending and twisting. Here we explore the relationship between base stacking, functional group occupancy of the DNA minor and major grooves, and DNA mechanical properties. We study double-helical DNA molecules substituting either inosine for guanosine or 2,6-diaminopurine for adenine. These DNA variants, respectively, remove or add an amino group from the DNA minor groove, with corresponding changes in hydrogen-bonding and base stacking energy. Using the techniques of ligase-catalyzed cyclization kinetics, atomic force microscopy, and force spectroscopy with optical tweezers, we show that these DNA variants have bending persistence lengths within the range of values reported for sequence-dependent variation of the natural DNA bases. Comparison with seven additional DNA variants that modify the DNA major groove reveals that DNA bending stiffness is not correlated with base stacking energy or groove occupancy. Data from circular dichroism spectroscopy indicate that base analog substitution can alter DNA helical geometry, suggesting a complex relationship among base stacking, groove occupancy, helical structure, and DNA bend stiffness.

Understanding the structural and dynamic consequences of DNA epigenetic modifications: computational insights into cytosine methylation and hydroxymethylation

We report a series of molecular dynamics (MD) simulations of up to a microsecond combined simulation time designed to probe epigenetically modified DNA sequences. More specifically, by monitoring the effects of methylation and hydroxymethylation of cytosine in different DNA sequences, we show, for the first time, that DNA epigenetic modifications change the molecule's dynamical landscape, increasing the propensity of DNA toward different values of twist and/or roll/tilt angles (in relation to the unmodified DNA) at the modification sites. Moreover, both the extent and position of different modifications have significant effects on the amount of structural variation observed. We propose that these conformational differences, which are dependent on the sequence environment, can provide specificity for protein binding.

E. coli Gyrase Fails to Negatively Supercoil Diaminopurine-Substituted DNA

Type II topoisomerases modify DNA supercoiling, and crystal structures suggest that they sharply bend DNA in the process. Bacterial gyrases are a class of type II topoisomerases that can introduce negative supercoiling by creating a wrap of DNA before strand passage. Isoforms of these essential enzymes were compared to reveal whether they can bend or wrap artificially stiffened DNA. Escherichia coli gyrase and human topoisomerase IIα were challenged with normal DNA or stiffer DNA produced by polymerase chain reaction reactions in which diaminopurine (DAP) replaced adenine deoxyribonucleotide triphosphates. On single DNA molecules twisted with magnetic tweezers to create plectonemes, the rates or pauses during relaxation of positive supercoils in DAP-substituted versus normal DNA were distinct for both enzymes. Gyrase struggled to bend or perhaps open a gap in DAP-substituted DNA, and segments of wider DAP DNA may have fit poorly into the N-gate of the human topoisomerase IIα. Pauses during processive activity on both types of DNA exhibited ATP dependence consistent with two pathways leading to the strand-passage-competent state with a bent gate segment and a transfer segment trapped by an ATP-loaded and latched N-gate. However, E. coli DNA gyrase essentially failed to negatively supercoil 35% stiffer DAP DNA.

SiRNA-induced mutation in HIV-1 polypurine tract region and its influence on viral fitness

Converting single-stranded viral RNA into double stranded DNA for integration is an essential step in HIV-1 replication. Initial polymerization of minus-strand DNA is primed from a host derived tRNA, whereas subsequent plus-strand synthesis requires viral primers derived from the 3' and central polypurine tracts (3' and cPPTs). The 5' and 3' termini of these conserved RNA sequence elements are precisely cleaved by RT-associated RNase H to generate specific primers that are used to initiate plus-strand DNA synthesis. In this study, siRNA wad used to produce a replicative HIV-1 variant contained G(-1)A and T(-16)A substitutions within/adjacent to the 3'PPT sequence. Introducing either or both mutations into the 3'PPT region or only the G(-1)A substitution in the cPPT region of NL4-3 produced infectious virus with decreased fitness relative to the wild-type virus. In contrast, introducing the T(-16)A or both mutations into the cPPT rendered the virus(es) incapable of replication, most likely due to the F185L integrase mutation produced by this nucleotide substitution. Finally, the effects of G(-1)A and T(-16)A mutations on cleavage of the 3'PPT were examined using an in vitro RNase H cleavage assay. Substrate containing both mutations was mis-cleaved to a greater extent than either wild-type substrate or substrate containing the T(-16)A mutation alone, which is consistent with the observed effects of the equivalent nucleotide substitutions on the replication fitness of NL4-3 virus. In conclusion, siRNA targeting of the HIV-1 3'PPT region can substantially suppress virus replication, and this selective pressure can be used to generate infectious virus containing mutations within or near the HIV-1 PPT. Moreover, in-depth analysis of the resistance mutations demonstrates that although virus containing a G(-1)A mutation within the 3'PPT is capable of replication, this nucleotide substitution shifts the 3'-terminal cleavage site in the 3'PPT by one nucleotide (nt) and significantly reduces viral fitness.

Mechanical Model of DNA Allostery

The importance of allosteric effects in DNA is becoming increasingly appreciated, but the underlying mechanisms remain poorly understood. In this work, we propose a general modeling framework to study DNA allostery. We describe DNA in a coarse-grained manner by intra-base pair and base pair step coordinates, complemented by groove widths. Quadratic deformation energy is assumed, yielding linear relations between the constraints and their effect. Model parameters are inferred from standard unrestrained, explicit-solvent molecular dynamics simulations of naked DNA. We applied the approach to study minor groove binding of diamidines and pyrrole-imidazole polyamides. The predicted DNA bending is in quantitative agreement with experiment and suggests that diamidine binding to the alternating TA sequence brings the DNA closer to the A-tract conformation, with potentially important functional consequences. The approach can be readily applied to other allosteric effects in DNA and generalized to model allostery in various molecular systems.

Roles of the amino group of purine bases in the thermodynamic stability of DNA base pairing

The energetic aspects of hydrogen-bonded base-pair interactions are important for the design of functional nucleotide analogs and for practical applications of oligonucleotides. The present study investigated the contribution of the 2-amino group of DNA purine bases to the thermodynamic stability of oligonucleotide duplexes under different salt and solvent conditions, using 2'-deoxyriboinosine (I) and 2'-deoxyribo-2,6-diaminopurine (D) as non-canonical nucleotides. The stability of DNA duplexes was changed by substitution of a single base pair in the following order: G•C > D•T ≈ I•C > A•T > G•T > I•T. The apparent stabilization energy due to the presence of the 2-amino group of G and D varied depending on the salt concentration, and decreased in the water-ethanol mixed solvent. The effects of salt concentration on the thermodynamics of DNA duplexes were found to be partially sequence-dependent, and the 2-amino group of the purine bases might have an influence on the binding of ions to DNA through the formation of a stable base-paired structure. Our results also showed that physiological salt conditions were energetically favorable for complementary base recognition, and conversely, low salt concentration media and ethanol-containing solvents were effective for low stringency oligonucleotide hybridization, in the context of conditions employed in this study.

FeON-FeOFF: the Helicobacter pylori Fur regulator commutates iron-responsive transcription by discriminative readout of opposed DNA grooves

Most transcriptional regulators bind nucleotide motifs in the major groove, although some are able to recognize molecular determinants conferred by the minor groove of DNA. Here we report a transcriptional commutator switch that exploits the alternative readout of grooves to mediate opposite output regulation for the same input signal. This mechanism accounts for the ability of the Helicobacter pylori Fur regulator to repress the expression of both iron-inducible and iron-repressible genes. When iron is scarce, Fur binds to DNA as a dimer, through the readout of thymine pairs in the major groove, repressing iron-inducible transcription (FeON). Conversely, on iron-repressible elements the metal ion acts as corepressor, inducing Fur multimerization with consequent minor groove readout of AT-rich inverted repeats (FeOFF). Our results provide first evidence for a novel regulatory paradigm, in which the discriminative readout of DNA grooves enables to toggle between the repression of genes in a mutually exclusive manner.

Exploring polymorphisms in B-DNA helical conformations

The traditional mesoscopic paradigm represents DNA as a series of base-pair steps whose energy response to equilibrium perturbations is elastic, with harmonic oscillations (defining local stiffness) around a single equilibrium conformation. In addition, base sequence effects are often analysed as a succession of independent XpY base-pair steps (i.e. a nearest-neighbour (NN) model with only 10 unique cases). Unfortunately, recent massive simulations carried out by the ABC consortium suggest that the real picture of DNA flexibility may be much more complex. The paradigm of DNA flexibility therefore needs to be revisited. In this article, we explore in detail one of the most obvious violations of the elastic NN model of flexibility: the bimodal distributions of some helical parameters. We perform here an in-depth statistical analysis of a very large set of MD trajectories and also of experimental structures, which lead to very solid evidence of bimodality. We then suggest ways to improve mesoscopic models to account for this deviation from the elastic regime.

Can We Accurately Describe the Structure of Adenine Tracts in B-DNA? Reference Quantum-Chemical Computations Reveal Overstabilization of Stacking by Molecular Mechanics

Sequence-dependent local variations of helical parameters, structure, and flexibility are crucial for molecular recognition processes involving B-DNA. A-tracts, i.e., stretches of several consecutive adenines in one strand that are in phase with the DNA helical repeat, mediate significant DNA bending. During the past few decades, there have been intense efforts to understand the sequence dependence of helical parameters in DNA. Molecular dynamics (MD) simulations can provide valuable insights into the molecular mechanism behind the relationship between sequence and structure. However, although recent improvements in empirical force fields have helped to capture many sequence-dependent B-DNA properties, several problems remain, such as underestimation of the helical twist and suspected underestimation of the propeller twist in A-tracts. Here, we employ reference quantum mechanical (QM) calculations, explicit solvent MD, and bioinformatics to analyze the underestimation of propeller twisting of A-tracts in simulations. Although we did not identify a straightforward explanation, we discovered two imbalances in the empirical force fields. The first was overestimation of stacking interactions accompanied by underestimation of base-pairing energy, which we attribute to anisotropic polarizabilities that are not reflected by the isotropic force fields. This may lead to overstacking with potentially important consequences for MD simulations of nucleic acids. The second observed imbalance was steric clash between A(N1) and T(N3) nitrogens of AT base pairs in force-field descriptions, resulting in overestimation of the AT pair stretch in MD simulations. We also substantially extend the available set of benchmark estimated CCSD(T)/CBS data for B-DNA base stacking and provide a code that allows the generation of diverse base-stacking geometries suitable for QM computations with predefined intra- and interbase pair parameters.

Molecular Dynamics Simulations of a Nucleosome and Free DNA

All atom molecular dynamics simulations (10ns) of a nucleosome and of its 146 basepairs of DNA free in solution have been conducted. DNA helical parameters (Roll, Tilt, Twist, Shift, Slide, Rise) were extracted from each trajectory to compare the conformation, effective force constants, persistence length measures, and fluctuations of nucleosomal DNA to free DNA. The conformation of DNA in the nucleosome, as determined by helical parameters, is found to be largely within the range of thermally accessible values obtained for free DNA. DNA is found to be less flexible on the nucleosome than when free in solution, however such measures are length scale dependent. A method for disassembling and reconstructing the conformation and dynamics of the nucleosome using Fourier analysis is presented. Long length variations in the conformation of nucleosomal DNA are identified other than those associated with helix repeat. These variations are required to create a proposed tetrasome conformation or to qualitatively reconstruct the 1.75 turns of the nucleosome's superhelix. Reconstruction of free DNA using selected long wavelength variations in conformation can produce either a left-handed or a right-handed superhelix. The long wavelength variations suggest 146 basepairs is a natural length of DNA to wrap around the histone core.

A duplex DNA model with regular inter-base-pair hydrogen bonds

It is well known that base-pair stacking is the main factor in stabilizing DNA duplex and plays an important role in determining DNA sequence-dependence. What is the dominant force in base-pair stacking? This fundamental biological question remains a challenging problem. Here, based on recent studies about the non-planarity of amino groups on DNA bases, we propose a new duplex DNA model, in which all base amino groups are non-planar and participate in forming regular inter-base-pair hydrogen bonds (IBP H-bonds). This model implies that IBP H-bonds are the dominant force stabilizing base-pair stacking and play a crucial role in determining the geometry and physical properties of sequence-dependent twisted stacking between adjacent base pairs. The model presents a new insight into the link, through regular IBP H-bonds, between base-sequence, fine structure and physical properties at dinucleotide step level, and provides an attractively concise, uniform and quantitative interpretation for various experimentally observed DNA sequence-dependent properties in terms of regular IBP H-bonds. It would provide a new approach to understanding the dynamics and underlying mechanisms of DNA sequence-dependent biological processes, sequence-structure-property relationships, DNA strand separation during replication and transcriptions, etc.

A measure of bending in nucleic acids structures applied to A-tract DNA

A method is proposed to measure global bending in DNA and RNA structures. It relies on a properly defined averaging of base-fixed coordinate frames, computes mean frames of suitably chosen groups of bases and uses these mean frames to evaluate bending. The method is applied to DNA A-tracts, known to induce considerable bend to the double helix. We performed atomistic molecular dynamics simulations of sequences containing the A(4)T(4) and T(4)A(4) tracts, in a single copy and in two copies phased with the helical repeat. Various temperature and salt conditions were investigated. Our simulations indicate bending by roughly 10 degrees per A(4)T(4) tract into the minor groove, and an essentially straight structure containing T(4)A(4), in agreement with electrophoretic mobility data. In contrast, we show that the published NMR structures of analogous sequences containing A(4)T(4) and T(4)A(4) tracts are significantly bent into the minor groove for both sequences, although bending is less pronounced for the T(4)A(4) containing sequence. The bending magnitudes obtained by frame averaging are confirmed by the analysis of superhelices composed of repeated tract monomers.

Non-additive interactions of nucleobases in model dinucleotide steps occurring in B-DNA crystals

Non-additivity of base-base interactions in all ten possible model dinucleotide steps were analyzed on MP2/aug-cc-pvDZ quantum chemistry level. Conformations of four nucleobases exactly matched to ones occurring in B-DNA crystals. In most of thw 162 analyzed tetramers both three- and four-body contributions are negligible except for d(GpG) steps. However, in these dinucleotides both contributions are always of opposite signs and in all cases the sum of all non-additive part of intermolecular interactions do not exceed 2.6 kcal mol(-1). This stands for less than 5% of the overall binding energy of dinucleotide steps. Also replacements of guanine with 8-oxoguanine in d(GpG) systems introduces non-additivity of the same magnitude as for canonical dinucleotides. It is observed linear relationships between values of total binding energy obtained in the tetramer basis set and estimated energy exclusively in dimers basis sets with assumption of pairwise additivities. For all analyzed dinucleotides steps there are also linear correlations between amount of non-additive contributions and magnitude of pairs interactions. Based on differences in electrostatic contribution to the total binding energy of four nucleobases and polarity of dinucleotide steps three distinct classes of dinucleotide steps were identified.

Comparison of intrinsic stacking energies of ten unique dinucleotide steps in A-RNA and B-DNA duplexes. Can we determine correct order of stability by quantum-chemical calculations?

High level ab initio methods have been used to study stacking interactions in ten unique base pair steps both in A-RNA and in B-DNA duplexes. The protocol for selection of geometries based on molecular dynamics (MD) simulations is proposed, and its suitability is demonstrated by comparison with stacking in steps at fiber diffraction geometries. It is shown that fiber diffraction geometries are not sufficiently accurate for interaction energy calculations. In addition, the protocol for selection of geometries based on MD simulations allows for the evaluation of the variability of the intrinsic stacking energies along the MD trajectories. The uncertainty in stacking energies (difference between the most and least stable geometry) due to the dynamical nature of systems can be, in some cases, as large as 3.0 kcal x mol(-1), which is almost 50% of the actual sequence dependence of base stacking energies (the energy difference between the most and least stable sequences). Thus, assessing the relative magnitude of the gas phase stacking energy using a single geometry for each sequence is insufficient to obtain an unambiguous order of gas phase stacking energies in canonical double helices. Though the ordering of ten unique dinucleotide steps cannot be definitive, some general conclusions were drawn. The stacking energies of base pair steps in A-RNA are more evenly separated compared to B-DNA, and their ordering is less sensitive to the dynamics of the system compared to be B-DNA. The most stable step both in B-DNA and A-RNA is the GC/GC [corrected] step that is well separated from the second most stable step CG/CG. [corrected] Also the least stable step (the CC/GG step) is well separated from the rest of the structures. The calculations further show that B-DNA stacking is favorable only marginally (on average by 1.14 kcal x mol(-1) per base pair step) over A-RNA stacking, and this difference vanishes after subtracting the stabilizing van der Waals effect of the thymine 5-methyl group that is absent in RNA. Basically, no correlation between the sequence dependence of gas phase stacking energies and the sequence dependence of DeltaG degrees(37) free energies used in nearest-neighbor models was found either for B-DNA or for A-RNA. This reflects the complexity of the balance of forces that are responsible for the sequence dependence of thermodynamics stability of nucleic acids, which masks the effect of the intrinsic interactions between the stacked base pairs.

Nature and magnitude of aromatic stacking of nucleic acid bases

This review summarises recent advances in quantum chemical calculations of base-stacking forces in nucleic acids. We explain in detail the very complex relationship between the gas-phase base-stacking energies, as revealed by quantum chemical (QM) calculations, and the highly variable roles of these interactions in nucleic acids. This issue is rarely discussed in quantum chemical and physical chemistry literature. We further extensively discuss methods that are available for base-stacking studies, complexity of comparison of stacking calculations with gas phase experiments, balance of forces in stacked complexes of nucleic acid bases, and the relation between QM and force field descriptions. We also review all recent calculations on base-stacking systems, including details analysis of the B-DNA stacking. Specific attention is paid to the highest accuracy QM calculations, to the decomposition of the interactions, and development of dispersion-balanced DFT methods. Future prospects of computational studies of base stacking are discussed.

Towards a molecular dynamics consensus view of B-DNA flexibility

We present a systematic study of B-DNA flexibility in aqueous solution using long-scale molecular dynamics simulations with the two more recent versions of nucleic acids force fields (CHARMM27 and parmbsc0) using four long duplexes designed to contain several copies of each individual base pair step. Our study highlights some differences between pambsc0 and CHARMM27 families of simulations, but also extensive agreement in the representation of DNA flexibility. We also performed additional simulations with the older AMBER force fields parm94 and parm99, corrected for non-canonical backbone flips. Taken together, the results allow us to draw for the first time a consensus molecular dynamics picture of B-DNA flexibility.

Computational modeling of the chromatin fiber

The packing of the genomic DNA in the living cell is essential for its biological function. While individual aspects of the genome architecture, such as DNA and nucleosome structure or the arrangement of chromosome territories are well studied, much information is missing for a unified description of cellular DNA at all its structural levels. Computer modeling can contribute to such a description. We present here some typical approaches to models of the chromatin fiber, including different amounts of detail in the description of the local nucleosome structure. The main results from our simulations are that the physical properties of the chromatin fiber can be well described by a simplified model consisting of cylinder-like nucleosomes connected by flexible DNA segments, with a geometry determined by the bending and twisting angles between nucleosomes. Randomness in the local geometry - such as random absence of linker histone H1 - leads to a dramatic increase in the chromatin fiber flexibility. Furthermore, we show that chromatin is much more flexible to bending than to stretching, and that the structure of the chromatin fiber favors the formation of sharp bends.

In vitro analysis of the effects of mutations in the G-tract of the human immunodeficiency virus type 1 polypurine tract on RNase H cleavage specificity

The recognition and precise cleavage of the polypurine tract (PPT) of the human immunodeficiency virus type 1 (HIV-1) is an essential step in HIV-1 reverse transcription. The accurate cleavage, and the subsequent removal, of the PPT by the RNase H activity of HIV-1 RT defines the left end of the double-stranded viral DNA genome, the substrate for integration into the host genome. Previous analyses have shown that mutations in the 3'-end (G-tract) of the PPT cause alterations in RNase H cleavage specificity. In particular, mutations at positions 2 and 5 in the G-tract increased the frequency of retention of PPT sequences in the 2-LTR circle junction. To better understand why these mutations affected PPT cleavage in vivo, we analyzed the cleavage of PPT substrates in vitro that contained altered sequences and unusual base substitutions. Our results, herein, confirm that mutations at positions 2 and 5 of the G-tract do significantly alter the cleavage specificity at the PPT/U3 junction, and further suggest that the miscleavages observed in vivo were due to an improper generation of the PPT primer, as opposed to its improper removal. Finally, our results point to the structure of the PPT, rather than the base-specific contacts between the PPT and HIV-1 RT, as the primary determinants of RNase H cleavage specificity at the PPT/U3 junction.

Nature of base stacking: reference quantum-chemical stacking energies in ten unique B-DNA base-pair steps

Base-stacking energies in ten unique B-DNA base-pair steps and some other arrangements were evaluated by the second-order Møller-Plesset (MP2) method, complete basis set (CBS) extrapolation, and correction for triple (T) electron-correlation contributions. The CBS(T) calculations were compared with decade-old MP2/6-31G*(0.25) reference data and AMBER force field. The new calculations show modest increases in stacking stabilization compared to the MP2/6-31G*(0.25) data and surprisingly large sequence-dependent variation of stacking energies. The absolute force-field values are in better agreement with the new reference data, while relative discrepancies between quantum-chemical (QM) and force-field values increase modestly. Nevertheless, the force field provides good qualitative description of stacking, and there is no need to introduce additional pair-additive electrostatic terms, such as distributed multipoles or out-of-plane charges. There is a rather surprising difference of about 0.1 A between the vertical separation of base pairs predicted by quantum chemistry and derived from crystal structures. Evaluations of different local arrangements of the 5'-CG-3' step indicate a sensitivity of the relative stacking energies to the level of calculation. Thus, describing quantitative relations between local DNA geometrical variations and stacking may be more complicated than usually assumed. The reference calculations are complemented by continuum-solvent assessment of solvent-screening effects.

Minor groove deformability of DNA: a molecular dynamics free energy simulation study

The conformational deformability of nucleic acids can influence their function and recognition by proteins. A class of DNA binding proteins including the TATA box binding protein binds to the DNA minor groove, resulting in an opening of the minor groove and DNA bending toward the major groove. Explicit solvent molecular dynamics simulations in combination with the umbrella sampling approach have been performed to investigate the molecular mechanism of DNA minor groove deformations and the indirect energetic contribution to protein binding. As a reaction coordinate, the distance between backbone segments on opposite strands was used. The resulting deformed structures showed close agreement with experimental DNA structures in complex with minor groove-binding proteins. The calculated free energy of minor groove deformation was approximately 4-6 kcal mol(-1) in the case of a central TATATA sequence. A smaller equilibrium minor groove width and more restricted minor groove mobility was found for the central AAATTT and also a significantly ( approximately 2 times) larger free energy change for opening the minor groove. The helical parameter analysis of trajectories indicates that an easier partial unstacking of a central TA versus AT basepair step is a likely reason for the larger groove flexibility of the central TATATA case.

Polymer chain models of DNA and chromatin

Many properties of the genome in the cell nucleus can be understood by modeling DNA and chromatin as a flexible polymer chain. This article introduces into current models for such a coarse-grained description and reviews some recent results from our own group. Examples given are the unrolling of DNA from the histone core and the response of the 30 nm chromatin fiber to mechanical stretching.

X-ray diffraction "fingerprinting" of DNA structure in solution for quantitative evaluation of molecular dynamics simulation

Solution state x-ray diffraction fingerprinting is demonstrated as a method for experimentally assessing the accuracy of molecular dynamics (MD) simulations. Fourier transforms of coordinate data from MD simulations are used to produce reciprocal space "fingerprints" of atomic pair distance correlations that are characteristic of the ensemble and are the direct numerical analogues of experimental solution x-ray diffraction (SXD). SXD experiments and MD simulations were carried out to test the ability of experiment and simulation to resolve sequence-dependent modifications in helix conformation for B-form DNA. SXD experiments demonstrated that solution-state poly(AT) and poly(A)-poly(T) duplex DNA sequences exist in ensembles close to canonical B-form and B'-form structures, respectively. In contrast, MD simulations analyzed in terms of SXD fingerprints are shown to deviate from experiment, most significantly for poly(A)-poly(T) duplex DNA. Compared with experiment, MD simulation shortcomings were found to include both mismatches in simulated conformer structures and number population within the ensembles. This work demonstrates an experimental approach for quantitatively evaluating MD simulations and other coordinate models to simulate biopolymer structure in solution and suggests opportunities to use solution diffraction data as experimental benchmarks for developing supramolecular force fields optimized for a range of in situ applications.

Towards an understanding of DNA recognition by the methyl-CpG binding domain 1

CpG methylation determines a variety of biological functions of DNA. The methylation signal is interpreted by proteins containing a methyl-CpG binding domain (MBDs). Based on the NMR structure of MBD1 complexed with methylated DNA we analysed the recognition mode by means of molecular dynamics simulations. As the protein is monomeric and recognizes a symmetrically methylated CpG step, the recognition mode is an asymmetric one. We find that the two methyl groups do not contribute equally to the binding energy. One methyl group is associated with the major part of the binding energy and the other one nearly does not contribute at all. The contribution of the two cytosine methyl groups to binding energy is calculated to be -3.6 kcal/mol. This implies a contribution of greater than two orders of magnitude to the binding constant. The conserved amino acid Asp32 is known to be essential for DNA binding by MBD1, but so far no direct contact with DNA has been observed. We detected a direct DNA base contact to Asp32. This could be the main reason for the importance of this amino acid. MBD contacts DNA exclusively in the major groove, the minor groove is reserved for histone contacts. We found a deformation of the minor groove shape due to complexation by MBD1, which indicates an information transfer between the major and the minor groove.