Effect of a neutralized phosphate backbone on the minor groove of B-DNA: molecular dynamics simulation studies

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.

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

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

Mono- and bis-pyrazolophthalazines: Design, synthesis, cytotoxic activity, DNA/HSA binding and molecular docking studies

In an attempt to find new potent cytotoxic compounds, several mono- and bis-pyrazolophthalazines 4a-m and 6a-h were synthesized through an efficient, one-pot, three- and pseudo five-component synthetic approach. All derivatives were evaluated for their in vitro cytotoxic activities against four human cancer cell lines of A549, HepG2, MCF-7, and HT29. Compound 4e showed low toxicity against normal cell lines (MRC-5 and MCF 10A, IC₅₀ > 200 µM) and excellent cytotoxic activity against A549 cell line with IC₅₀ value of 1.25 ± 0.19 µM, which was 1.8 times more potent than doxorubicin (IC₅₀ = 2.31 ± 0.13 µM). In addition, compound 6c exhibited remarkable cytotoxic activity against A549 and MCF-7 cell lines (IC₅₀ = 1.35 ± 0.12 and 0.49 ± 0.01 µM, respectively), more than two-fold higher than that of doxorubicin. The binding properties of the best active mono- and bis-pyrazolophthalazine (4e and 6c) with HSA and DNA were fully evaluated by various techniques including UV-Vis absorption, circular dichroism (CD), Zeta potential and dynamic light scattering analyses indicating interaction of the compounds with the secondary structure of HSA and significant change of DNA conformation, presumably via a groove binding mechanism. Additionally, molecular docking and site-selective binding studies confirmed the fundamental interaction of compounds 4e and 6c with base pairs of DNA. Compounds 4e and 6c showed promising features to be considered as potential lead structures for further studies in cancer therapy.

Structural effects of Diamond nanoparticles and Paclitaxel combination on calf thymus DNA

The combination effects of nanodiamonds (NDs) and Paclitaxel (PTX) on the DNA structure were examined. The UV-Visible, steady-state and time-resolved fluorescence spectroscopy, CD, viscosity and zeta potential results showed that PTX + NDs could form a complex via groove binding mechanism. The values of binding constants, ΔG° and ΔH° and ΔS° values showed that PTX + NDs interact strongly with DNA and the hydrophobic force plays main role in this interaction. The ΔG25ο and T_m study indicated the instability of DNA in presence of PTX + NDs. This study demonstrated that NDs could enhance the effect of PTX on DNA structure as well as its affinity and binding to DNA.

Graph-Theoretic Analysis of Monomethyl Phosphate Clustering in Ionic Solutions

All-atom molecular dynamics simulations combined with graph-theoretic analysis reveal that clustering of monomethyl phosphate dianion (MMP^2-) is strongly influenced by the types and combinations of cations in the aqueous solution. Although Ca²⁺ promotes the formation of stable and large MMP^2- clusters, K⁺ alone does not. Nonetheless, clusters are larger and their link lifetimes are longer in mixtures of K⁺ and Ca²⁺. This "synergistic" effect depends sensitively on the Lennard-Jones interaction parameters between Ca²⁺ and the phosphorus oxygen and correlates with the hydration of the clusters. The pronounced MMP^2- clustering effect of Ca²⁺ in the presence of K⁺ is confirmed by Fourier transform infrared spectroscopy. The characterization of the cation-dependent clustering of MMP^2- provides a starting point for understanding cation-dependent clustering of phosphoinositides in cell membranes.

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

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

DNA translocation by human uracil DNA glycosylase: role of DNA phosphate charge

Human DNA repair glycosylases must encounter and inspect each DNA base in the genome to discover damaged bases that may be present at a density of <1 in 10 million normal base pairs. This remarkable example of specific molecular recognition requires a reduced dimensionality search process (facilitated diffusion) that involves both hopping and sliding along the DNA chain. Despite the widely accepted importance of facilitated diffusion in protein-DNA interactions, the molecular features of DNA that influence hopping and sliding are poorly understood. Here we explore the role of the charged DNA phosphate backbone in sliding and hopping by human uracil DNA glycosylase (hUNG), which is an exemplar that efficiently locates rare uracil bases in both double-stranded DNA and single-stranded DNA. Substitution of neutral methylphosphonate groups for anionic DNA phosphate groups weakened nonspecific DNA binding affinity by 0.4-0.5 kcal/mol per substitution. In contrast, sliding of hUNG between uracil sites embedded in duplex and single-stranded DNA substrates persisted unabated when multiple methylphosphonate linkages were inserted between the sites. Thus, a continuous phosphodiester backbone negative charge is not essential for sliding over nonspecific DNA binding sites. We consider several alternative mechanisms for these results. A model consistent with previous structural and nuclear magnetic resonance dynamic results invokes the presence of open and closed conformational states of hUNG. The open state is short-lived and has weak or nonexistent interactions with the DNA backbone that are conducive for sliding, and the populated closed state has stronger interactions with the phosphate backbone. These data suggest that the fleeting sliding form of hUNG is a distinct weakly interacting state that facilitates rapid movement along the DNA chain and resembles the transition state for DNA dissociation.

Effects of an unnatural base pair replacement on the structure and dynamics of DNA and neighboring water and ions

Incorporating small molecule probes into biomolecular systems to report on local structure and dynamics is a powerful strategy that underlies a wide variety of experimental techniques, including fluorescence, electron paramagnetic resonance (EPR), and Forster resonance energy transfer (FRET) measurements. When an unnatural probe is inserted into a protein or DNA, the degree to which the presence of the probe has perturbed the local structure and dynamics it was intended to study is always an important concern. Here, molecular dynamics (MD) simulations are used to systematically study the effect of replacing a DNA base pair with a fluorescent probe, coumarin 102 deoxyriboside, at six unique sites along an A-tract DNA dodecamer. While the overall structure of the DNA oligonucleotide remains intact, replacement of A*T base pairs leads to widespread structural and dynamic perturbations up to four base pairs away from the probe site, including widening of the minor groove and increased DNA flexibility. New DNA conformations, not observed in the native sequence, are sometimes found in the vicinity of the probe and its partner abasic site analog. Strong correlations are demonstrated between DNA surface topology and water mobility.

The energetic contribution of induced electrostatic asymmetry to DNA bending by a site-specific protein

DNA bending can be promoted by reducing the net negative electrostatic potential around phosphates on one face of the DNA, such that electrostatic repulsion among phosphates on the opposite face drives bending toward the less negative surface. To provide the first assessment of energetic contribution to DNA bending when electrostatic asymmetry is induced by a site-specific DNA binding protein, we manipulated the electrostatics in the EcoRV endonuclease-DNA complex by mutation of cationic side chains that contact DNA phosphates and/or by replacement of a selected phosphate in each strand with uncharged methylphosphonate. Reducing the net negative charge at two symmetrically located phosphates on the concave DNA face contributes -2.3 kcal mol(-1) to -0.9 kcal mol(-1) (depending on position) to complex formation. In contrast, reducing negative charge on the opposing convex face produces a penalty of +1.3 kcal mol(-1). Förster resonance energy transfer experiments show that the extent of axial DNA bending (about 50°) is little affected in modified complexes, implying that modification affects the energetic cost but not the extent of DNA bending. Kinetic studies show that the favorable effects of induced electrostatic asymmetry on equilibrium binding derive primarily from a reduced rate of complex dissociation, suggesting stabilization of the specific complex between protein and markedly bent DNA. A smaller increase in the association rate may suggest that the DNA in the initial encounter complex is mildly bent. The data imply that protein-induced electrostatic asymmetry makes a significant contribution to DNA bending but is not itself sufficient to drive full bending in the specific EcoRV-DNA complex.

On the molecular basis of uracil recognition in DNA: comparative study of T-A versus U-A structure, dynamics and open base pair kinetics

Uracil (U) can be found in DNA as a mismatch paired either to adenine (A) or to guanine (G). Removal of U from DNA is performed by a class of enzymes known as uracil–DNA–glycosylases (UDG). Recent studies suggest that recognition of U–A and U–G mismatches by UDG takes place via an extra-helical mechanism. In this work, we use molecular dynamics simulations to analyze the structure, dynamics and open base pair kinetics of U–A base pairs relative to their natural T–A counterpart in 12 dodecamers. Our results show that the presence of U does not alter the local conformation of B-DNA. Breathing dynamics and base pair closing kinetics are only weakly dependent on the presence of U versus T, with open T–A and U–A pairs lifetimes in the nanosecond timescale. Additionally, we observed spontaneous base flipping in U–A pairs. We analyze the structure and dynamics for this event and compare the results to available crystallographic data of open base pair conformations. Our results are in agreement with both structural and kinetic data derived from NMR imino proton exchange measurements, providing the first detailed description at the molecular level of elusive events such as spontaneous base pair opening and flipping in mismatched U–A sequences in DNA. Based on these results, we propose that base pair flipping can occur spontaneously at room temperature via a 3-step mechanism with an open base pair intermediate. Implications for the molecular basis of U recognition by UDG are discussed.

Intensities of DNA ion-phosphate modes in the low-frequency Raman spectra

The Raman intensities of counterion vibrations with respect to the phosphate groups of the double-helix backbone (ion-phosphate modes) in the low-frequency spectra (< 200 cm(-1)) of B -DNA with different alkali metal counterions have been calculated using the model for DNA conformational vibrations and the valence-optic approach. The results have showed that the spectra of DNA with heavy counterions (Rb(+) and Cs(+)) differ from the spectra of DNA with light counterions (Na(+) and K(+)). The calculated spectra of DNA with heavy counterions are characterized by intensive modes of ion-phosphate vibrations that form one united band near 115 cm(-1). Ion-phosphate modes in the spectra of DNA with light counterions are characterized by higher frequencies (near 180 cm(-1)) and much lower intensity. Our calculations explain why the modes of ion-phosphate vibrations are observed in Cs-DNA spectra rather than in Na-DNA. The determined sensitivity of the intensities of DNA low-frequency spectra to the counterion type proves the existence of the ion-phosphate modes.

Site-specific DNA structural and dynamic features revealed by nucleotide-independent nitroxide probes

In site-directed spin labeling, a covalently attached nitroxide probe containing a chemically inert unpaired electron is utilized to obtain information on the local environment of the parent macromolecule. Studies presented here examine the feasibility of probing local DNA structural and dynamic features using a class of nitroxide probes that are linked to chemically substituted phosphorothioate positions at the DNA backbone. Two members of this family, designated as R5 and R5a, were attached to eight different sites of a dodecameric DNA duplex without severely perturbing the native B-form conformation. Measured X-band electron paramagnetic resonance (EPR) spectra, which report on nitroxide rotational motions, were found to vary depending on the location of the label (e.g., duplex center vs termini) and the surrounding DNA sequence. This indicates that R5 and R5a can provide information on the DNA local environment at the level of an individual nucleotide. As these probes can be attached to arbitrary nucleotides within a nucleic acid sequence, they may provide a means to "scan" a given DNA molecule in order to interrogate its local structural and dynamic features.

Conformational tuning of magnetic interactions in metal-DNA complexes

The alignment of Cu(2+) ions along a modified DNA helix is studied with either hydroxypyridone (H) or bis(salicylaldehyde)ethylenediamine (S-en) metalated base pairs (MBPs). The conformational motion of H-MBP leads to the interlinking of the H-MBPs by an extended Cu-O network that is ferromagnetic, whereas the conformational freezing of the S-en-MBP leads to an ordered pairwise-stacked arrangement that is weakly antoferrimagnetic.

Ability of polymerase eta and T7 DNA polymerase to bypass bulge structures

DNA misalignment occurs in homopolymer tracts during replication and can lead to frameshift mutations. Polymerase (pol) recognition of primer-templates containing bulge structures and the transmission of a bulge through a polymerase binding site or replication complex are important components of frameshift mutagenesis. In this report, we describe the interaction of the catalytic core of pol eta with primer-templates containing bulge structures by single round primer extension. We found that pol eta could stabilize a frayed primer terminus, which enhances its ability to extend primer-templates containing bulges. Based on methylphosphonate-DNA mapping, pol eta interacts with the single strand template but not appreciably with the template strand of the DNA stem greater than two nucleotides from the primer terminus. These latter characteristics, combined with the ability to stabilize a frayed primer terminus, may explain why primer-templates containing template bulges are extended so efficiently by pol eta. Although pol eta could accommodate large bulges and continue synthesis without obstruction, bulge structures in the template, but not in the primer, caused termination of the T7 DNA replication complex. Terminations occurred when the template bulge neared the helix-loop-helix domain of the polymerase thumb. Terminations were not observed, however, when bulge structures approached the site of interaction of the DNA with the extended thumb and thioredoxin. At low temperature, however, terminations did occur at this site.

Accelerated molecular dynamics: a promising and efficient simulation method for biomolecules

Many interesting dynamic properties of biological molecules cannot be simulated directly using molecular dynamics because of nanosecond time scale limitations. These systems are trapped in potential energy minima with high free energy barriers for large numbers of computational steps. The dynamic evolution of many molecular systems occurs through a series of rare events as the system moves from one potential energy basin to another. Therefore, we have proposed a robust bias potential function that can be used in an efficient accelerated molecular dynamics approach to simulate the transition of high energy barriers without any advance knowledge of the location of either the potential energy wells or saddle points. In this method, the potential energy landscape is altered by adding a bias potential to the true potential such that the escape rates from potential wells are enhanced, which accelerates and extends the time scale in molecular dynamics simulations. Our definition of the bias potential echoes the underlying shape of the potential energy landscape on the modified surface, thus allowing for the potential energy minima to be well defined, and hence properly sampled during the simulation. We have shown that our approach, which can be extended to biomolecules, samples the conformational space more efficiently than normal molecular dynamics simulations, and converges to the correct canonical distribution.

Nuclear aggregates of polyamines are supramolecular structures that play a crucial role in genomic DNA protection and conformation

In a previous study we showed that natural polyamines interact in the nuclear environment with phosphate groups to form molecular aggregates [nuclear aggregates of polyamines (NAPs)] with estimated molecular mass values of 8000, 4800 and 1000 Da. NAPs were found to interact with genomic DNA, influence its conformation and interfere with the action of nucleases. In the present work, we demonstrated that NAPs protect naked genomic DNA from DNase I, whereas natural polyamines (spermine, spermidine and putrescine) fail to do so. In the context of DNA protection, NAPs induced noticeable changes in DNA conformation, which were revealed by temperature-dependent modifications of DNA electrophoretic properties. In addition, we presented, for NAPs, a structural model of polyamine aggregation into macropolycyclic compounds. We believe that NAPs are the sole biological forms by which polyamines efficiently protect genomic DNA against DNase I, while maintaining its dynamic structure.

Mechanisms of coupling between DNA recognition specificity and catalysis in EcoRI endonuclease

Proteins that bind to specific sites on DNA often do so in order to carry out catalysis or specific protein-protein interaction while bound to the recognition site. Functional specificity is enhanced if this second function is coupled to correct DNA site recognition. To analyze the structural and energetic basis of coupling between recognition and catalysis in EcoRI endonuclease, we have studied stereospecific phosphorothioate (PS) or methylphosphonate (PMe) substitutions at the scissile phosphate GpAATTC or at the adjacent phosphate GApATTC in combination with molecular-dynamics simulations of the catalytic center with bound Mg2+. The results show the roles in catalysis of individual phosphoryl oxygens and of DNA distortion and suggest that a "crosstalk ring" in the complex couples recognition to catalysis and couples the two catalytic sites to each other.

The contribution of phosphate-phosphate repulsions to the free energy of DNA bending

DNA bending is important for the packaging of genetic material, regulation of gene expression and interaction of nucleic acids with proteins. Consequently, it is of considerable interest to quantify the energetic factors that must be overcome to induce bending of DNA, such as base stacking and phosphate–phosphate repulsions. In the present work, the electrostatic contribution of phosphate–phosphate repulsions to the free energy of bending DNA is examined for 71 bp linear and bent-form model structures. The bent DNA model was based on the crystallographic structure of a full turn of DNA in a nucleosome core particle. A Green's function approach based on a linear-scaling smooth conductor-like screening model was applied to ascertain the contribution of individual phosphate–phosphate repulsions and overall electrostatic stabilization in aqueous solution. The effect of charge neutralization by site-bound ions was considered using Monte Carlo simulation to characterize the distribution of ion occupations and contribution of phosphate repulsions to the free energy of bending as a function of counterion load. The calculations predict that the phosphate–phosphate repulsions account for ∼30% of the total free energy required to bend DNA from canonical linear B-form into the conformation found in the nucleosome core particle.

Structural selectivity of aromatic diamidines

Competition dialysis was used to study the interactions of 13 substituted aromatic diamidine compounds with 13 nucleic acid structures and sequences. The results show a striking selectivity of these compounds for the triplex structure poly dA:(poly dT)2, a novel aspect of their interaction with nucleic acids not previously described. The triplex selectivity of selected compounds was confirmed by thermal denaturation studies. Triplex selectivity was found to be modulated by the location of amidine substituents on the core phenyl-furan-phenyl ring scaffold. Molecular models were constructed to rationalize the triplex selectivity of DB359, the most selective compound in the series. Its triplex selectivity was found to arise from optimal ring stacking on base triplets, along with proper positioning of its amidine substituents to occupy the minor and the major-minor grooves of the triplex. New insights into the molecular recognition of nucleic acid structures emerged from these studies, adding to the list of available design principles for selectively targeting DNA and RNA.

A molecular dynamics simulation study of polyamine? and sodium?DNA. Interplay between polyamine binding and DNA structure

Four different molecular dynamics (MD) simulations have been performed for infinitely long ordered DNA molecules with different counterions, namely the two natural polyamines spermidine(3+) (Spd3+) and putrescine(2+) (Put2+), the synthetic polyamine diaminopropane(2+) (DAP2+), and the simple monovalent cation Na+. All systems comprised a periodical hexagonal cell with three identical DNA decamers, 15 water molecules per nucleotide, and counterions balancing the DNA charge. The simulation setup mimics the DNA state in oriented DNA fibers, previously studied using NMR and other experimental methods. In this paper the interplay between polyamine binding and local DNA structure is analyzed by investigating how and if the minor groove width of DNA depends on the presence and dynamics of the counterions. The results of the MD simulations reveal principal differences in the polyamine-DNA interactions between the natural [spermine(4+), Spd3+, Put2+] and the synthetic (DAP2+) polyamines.

Effect of phosphorothioate chirality on the grooves of DNA double helices: a molecular dynamics study

Phosphorothioate oligonucleotides (PS-ODNs) have gained considerable attention in drug therapy, primarily as potent antisense or antigene oligomers, which bind to specific DNA or mRNA sequences and lead to transcriptional or translational arrest. These are obtained by substituting one of the anionic oxygen of the phosphate group by a sulfur atom, which introduces chirality to the phosphorus atom of the DNA backbone. In this molecular dynamics simulation study, structural parameters like groove widths, environmental parameters like hydration or cation binding, and electrostatic energy surfaces of both the chiral forms of DNA/PS-DNA duplexes were assessed and compared with that of a normal DNA. Results indicate that, PS-S form with its sulfur atoms facing the minor groove has a widened minor groove, while the scenario is reverse for the PS-R form. Further analysis reveals the existence of several factors like large van der Waals radius of sulfur and the effect it has on its neighboring hydration pattern along with the net electrostatic environment, influencing such structural alterations. This also indicates, for the first time, the effect of absolute phosphorothioate chirality on the global structure of a DNA/PS-DNA hybrid that otherwise resembles a regular B-DNA structure.

Daunomycin Intercalation Stabilizes Distinct Backbone Conformations of DNA

Daunomycin is a widely used antibiotic of the anthracycline family. In the present study we reveal the structural properties and important intercalator-DNA interactions by means of molecular dynamics. As most of the X-ray structures of DNA-daunomycin intercalated complexes are short hexamers or octamers of DNA with two drug molecules per doublehelix we calculated a self complementary 14-mer oligodeoxyribonucleotide duplex d(CGCGCGATCGCGCG)2 in the B-form with two putative intercalation sites at the 5'-CGA-3' step on both strands. Consequently we are able to look at the structure of a 1:1 complex and exclude crystal packing effects normally encountered in most of the X-ray crystallographic studies conducted so far. We performed different 10 to 20 ns long molecular dynamics simulations of the uncomplexed DNA structure, the DNA-daunomycin complex and a 1:2 complex of DNA-daunomycin where the two intercalator molecules are stacked into the two opposing 5'-CGA-3' steps. Thereby--in contrast to X-ray structures--a comparison of a complex of only one with a complex of two intercalators per doublehelix is possible. The chromophore of daunomycin is intercalated between the 5'-CG-3' bases while the daunosamine sugar moiety is placed in the minor groove. We observe a flexibility of the dihedral angle at the glycosidic bond, leading to three different positions of the ammonium group responsible for important contacts in the minor groove. Furthermore a distinct pattern of BI and BII around the intercalation site is induced and stabilized. This indicates a transfer of changes in the DNA geometry caused by intercalation to the DNA backbone.

DNA-thumb interactions and processivity of T7 DNA polymerase in comparison to yeast polymerase eta

The replicative polymerase of bacteriophage T7 is structurally and mechanistically well characterized. The crystal structure of T7 DNA polymerase or gene 5 protein complexed to its processivity factor, Escherichia coli thioredoxin, a primer-template, and a dideoxynucleotide reveals how this enzyme interacts with the 3'-end of the primer-template, but does not show how thioredoxin confers processivity to the polymerase. In the crystal structure highly conserved amino acids Asn(335) and Ser(338) of the thumb subdomain of T7 DNA polymerase are seen to interact with phosphates 7 and 8 of the DNA template strand. Results with a mutant T7 DNA polymerase in which aliphatic residues are substituted for these amino acids and experiments with different length and methylphosphonate-modified primer-templates demonstrate that these interactions are essential for processive synthesis and d(A.T)(n) tract bypass. Our data with methylphosphonate-modified DNA suggests that thioredoxin confers processivity to T7 DNA polymerase in part by causing an interaction with the phosphate backbone or minor groove of DNA. Residues Asn(335) and Ser(338) may also function with a nearby helix-loop-helix motif located at residues 339-372 to enclose the DNA during processive synthesis. Our results suggest that this structure must be held close to the DNA by ionic interactions to function. These interactions also allow for DNA sliding but physically block the passage of a 3T bulge in the template. In contrast, yeast polymerase eta, a polymerase that non-mutagenically repairs cis-syn thymidine dimers, allows the same bulge to slide past its thumb subdomain during synthesis. A relaxed thumb interaction with the DNA could account for the notably low processivity of polymerase eta.

Altered structural fluctuations in duplex RNA versus DNA: a conformational switch involving base pair opening

DNA and RNA are known to have different structural properties. In the present study, molecular dynamics (MD) simulations on a series of RNA and DNA duplexes indicate differential structural flexibility for the two classes of oligonucleotides. In duplex RNA, multiple base pairs experienced local opening events into the major groove on the nanosecond time scale, while such events were not observed in the DNA simulations. Three factors are indicated to be responsible for the base opening events in RNA: solvent-base interactions, 2'OH(n)-O4'(n+1) intra-strand hydrogen bonding, and enhanced rigid body motion of RNA at the nucleoside level. Water molecules in the major groove of RNA contribute to initiation of base pair opening. Stabilization of the base pair open state is due to a 'conformational switch' comprised of 2'OH(n)-O4'(n+1) hydrogen bonding and a rigid body motion of the nucleoside moiety in RNA. This rigid body motion is associated with decreased flexibility of the glycosyl linkage and sugar moieties in A-form structures. The observed opening rates in RNA are consistent with the imino proton exchange experiments for AU base pairs, although not for GC base pairs, while structural and flexibility changes associated with the proposed conformational switch are consistent with survey data of RNA and DNA crystal structures. The possible relevance of base pair opening events in RNA to its many biological functions is discussed.

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

The fine structure of the DNA double helix and a number of its physical properties depend upon nucleotide sequence. This includes minor groove width, the propensity to undergo the B-form to A-form transition, sequence-directed curvature, and cation localization. Despite the multitude of studies conducted on DNA, it is still difficult to appreciate how these fundamental properties are linked to each other at the level of nucleotide sequence. We demonstrate that several sequence-dependent properties of DNA can be attributed, at least in part, to the sequence-specific localization of cations in the major and minor grooves. We also show that effects of cation localization on DNA structure are easier to understand if we divide all DNA sequences into three principal groups: A-tracts, G-tracts, and generic DNA. The A-tract group of sequences has a peculiar helical structure (i.e., B*-form) with an unusually narrow minor groove and high base-pair propeller twist. Both experimental and theoretical studies have provided evidence that the B*-form helical structure of A-tracts requires cations to be localized in the minor groove. G-tracts, on the other hand, have a propensity to undergo the B-form to A-form transition with increasing ionic strength. This property of G-tracts is directly connected to the observation that cations are preferentially localized in the major groove of G-tract sequences. Generic DNA, which represents the vast majority of DNA sequences, has a more balanced occupation of the major and minor grooves by cations than A-tracts or G-tracts and is thereby stabilized in the canonical B-form helix. Thus, DNA secondary structure can be viewed as a tug of war between the major and minor grooves for cations, with A-tracts and G-tracts each having one groove that dominates the other for cation localization. Finally, the sequence-directed curvature caused by A-tracts and G-tracts can, in both cases, be explained by the cation-dependent mismatch of A-tract and G-tract helical structures with the canonical B-form helix of generic DNA (i.e., a cation-dependent junction model).

The role of minor groove functional groups in DNA hydration

Here we describe the crystal structure of modified [d(CGCGAATTCGCG)]2 refined to 2.04 A. The modification, which affects only the two thymines at the central ApT step, involves isosteric removal of the 2-keto oxygen atoms and substitution of the N1 nitrogen with carbon. The crystal structure reveals the ability of this modified thymine to effectively base pair with adenine in [d(CGCGAAtTCGCG)]2. The structure also suggests that the minor groove 'spine of hydration' is destabilized but essentially intact.