B-DNA Structure Research Articles

What history tells us IX. Z-DNA: when nature is not opportunistic

In December 1979, Alexander Rich and his collaborators described in Nature the structure of a left-handed double helical DNA fragment (Wang et al 1979). This structure differed signifi cantly from the B-DNA structure described 26 years before by Watson and Crick (fi gure 1). It was named ZDNA because the deoxyribose-phosphate backbone follows a zig-zag course instead of the regular one found in the Bstructure. But the last letter of the alphabet was also chosen to indicate that the new form was highly different from the previously described A and B structures. This discovery was not totally unexpected: since Watson and Crick, many other models had been proposed for DNA, including left-handed models, just like Z-DNA. Nevertheless, it was a big shock, and it attracted a lot of attention. From the beginning, Z-DNA was a structure in search of a function. From 1981 to 1983, preliminary results suggested that Z-DNA might exist in physiological conditions, and play a role in the regulation of transcription. No obvious major experimental support was further obtained. In fact, negative results were reported in the following years and interest in Z-DNA progressively declined. The strange fate of Z-DNA raises interesting historical issues. The fi rst concerns the relation between structure and function, a hallmark of molecular biology, in which the structure of macromolecules explains their function. What happens when the determination of the structure precedes the characterization of the potential functions it could fulfi l? This new structure, and its potential roles, fi tted quite well with the expectations of the numerous researchers working on the regulation of gene transcription in those years. But the expectations may have been an illusion, supported by an overly opportunistic vision of nature in which what is possible and what exists have a (too) strong tendency to coincide. Series

Synthetic metallomolecules as agents for the control of DNA structure

This tutorial review summarises B-DNA structure and metallomolecule binding modes and illustrates some DNA structures induced by molecules containing metallic cations. The effects of aquated metal ions, cobalt amines, ruthenium octahedral metal complexes, metallohelicates and platinum complexes such as cis-platin are discussed alongside the techniques of NMR, X-ray crystallography, gel electrophoresis, circular dichroism, linear dichroism and molecular dynamics. The review will be of interest to people interested in both DNA structure and roles of metallomolecules in biological systems.

Charge Localization in Stacked Radical Cation DNA Base Pairs and the Benzene Dimer Studied by Self-Interaction Corrected Density-Functional Theory

The incomplete cancellation of the electron self-interaction can be a serious shortcoming of density-functional theory especially when treating odd-electron systems. In this work, several popular and potentially viable correction schemes are applied in order to characterize the electronic structure of stacked molecular pairs, consisting of a neutral molecule and adjacent radical cation, as a function of separation distance. The unphysical sharing of the positive charge between adjacent molecules separated by 6-7 A is corrected for by applying a new empirical scheme proposed by VandeVondele and Sprik [Phys. Chem. Chem. Phys. 2005, 7, 1363] with a unique choice of parameters. This method is subsequently applied to characterize the electronic structure of two neighboring guanines excised from a canonical Arnott B-DNA structure and will be used in future investigations of certain model DNA fibers.

Recognition of B-DNA by Neomycin−Hoechst 33258 Conjugates

Recent developments have indicated that aminoglycoside binding is limited not to RNA but to nucleic acids that, like RNA, adopt conformations similar to the A-form. We have further sought to expand the utility of aminoglycoside binding to B-DNA structures by conjugating neomycin, an aminoglycoside antibiotic, with the B-DNA minor groove binding ligand Hoechst 33258. Described herein are novel neomycin-Hoechst 33258 conjugates developed for exploring B-DNA groove recognition. We have varied the two reported conjugates in linker length and composition in an effort to improve our understanding of the spatial differences that define B-DNA binding. Spectroscopic studies such as ultraviolet (UV) melting, isothermal fluorescence titrations, differential scanning calorimetry (DSC), and circular dichroism (CD) together illustrate the mode of binding by such conjugates. Both conjugates exhibit enhanced thermal stabilization of A.T rich duplexes when compared to Hoechst 33258.

Watson–Crick Base‐Pairing Properties of Nucleic Acid Analogues with Stereocontrolled α and β Torsion Angles (α,β‐D‐CNAs)

Constrained nucleic acids (α,β-D-CNAs) have an ethylene bridge that locks the α and β backbone torsion angles of DNA in the canonical (gauche(−),trans) conformation, as in B-form DNA (see picture), or in the noncanonical (gauche(+),trans) conformation. Canonical α,β-D-CNAs enhance the duplex-forming ability of oligonucleotides with complementary DNA, while noncanonical α,β-D-CNAs stabilize distorted B-DNA structures. Supporting information for this article is available on the WWW under http://www.wiley-vch.de/contents/jc_2002/2006/z504475_s.pdf or from the author. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.

DNA Interaction with Human Serum Albumin Studied by Affinity Capillary Electrophoresis and FTIR Spectroscopy

The question addressed in this study is how does the protein-DNA complexation affect the structure and dynamics of DNA and protein in aqueous solution. We examined the interaction of calf-thymus DNA with human serum albumin (HSA) in aqueous solution at physiological conditions, using constant DNA concentration of 12.5 mM (phosphate) and various HSA contents 0.25 to 2% or 0.04 to 0.3 mM. Affinity capillary electrophoresis and FTIR spectroscopic methods were used to determine the protein binding mode, the association constant, sequence preference, and the biopolymer secondary structural changes in the HSA-DNA complexes. Spectroscopic evidence showed two types of HSA-DNA complexes with strong binding of K(1) = 4.5 x 10(5) M(-1) and weak binding of K(2) = 6.10 x 10(4) M(-1). The two major binding sites were located on the G-C bases and the backbone PO(2) group. The protein-DNA interaction stabilizes the HSA secondary structure. A minor alteration of B-DNA structure was observed, while no major protein conformational changes occurred.

Sequence-Dependent Conformational Energy of DNA Derived from Molecular Dynamics Simulations: Toward Understanding the Indirect Readout Mechanism in Protein−DNA Recognition

Sequence dependence of DNA conformation plays a crucial role in its recognition by proteins and ligands. To clarify the relationship between sequence and conformation, it is necessary to quantify the conformational energy and specificity of DNA. Here, we make a systematic analysis of dodecamer DNA structures including all the 136 unique tetranucleotide sequences at the center by molecular dynamics simulations. Using a simplified conformational model with six parameters to describe the geometry of adjacent base pairs and harmonic potentials along these coordinates, we estimated the equilibrium conformational parameters and the harmonic potentials of mean force for the central base-pair steps from many trajectories of the simulations. This enabled us to estimate the conformational energy and the specificity for any given DNA sequence and structure. We tested our method by using sequence-structure threading to estimate the conformational energy and the Z-score as a measure of specificity for many B-DNA and A-DNA crystal structures. The average Z-scores were negative for both kinds of structures, indicating that the potential of mean force from the simulation is capable of predicting sequence specificity for the crystal structures and that it may be used to study the sequence specificity of both types of DNA. We also estimated the positional distribution of conformational energy and Z-score within DNA and showed that they are strongly position dependent. This analysis enabled us to identify particular conformations responsible for the specificity. The presented results will provide an insight into the mechanisms of DNA sequence recognition by proteins and ligands.

Interaction of Arsenic Trioxide As2O3with DNA and RNA

Arsenic salts have been used for centuries to treat a variety of medical conditions ranging from infectious disease to cancer. More recently, trivalent arsenic trioxide was found to exhibit high antitumor activity towards hematological malignancies. Even though much is known about antitumor activity and DNA damage by As2O3, there has been no report on the interaction of arsenic trioxide with isolated DNA or RNA. Therefore, it was of interest to examine the interaction of As2O3 with DNA and RNA in aqueous solution at physiological pH. FTIR and UV-visible difference spectroscopic methods were used to characterize the nature of drug-DNA and drug-RNA interactions and to determine the As binding site, the binding constant, the sequence selectivity, the helix stability, and the biopolymer secondary structure in the As2O3-polynucleotide complexes in vitro. The FTIR spectroscopic studies were conducted with As2O3-polynucleotide (phosphate) ratios of 1/40, 1/20, 1/10, and 1/5, with a final DNA (P) or RNA (P) concentration of 6.25 mmol/l. Spectroscopic results showed As2O3 binds to DNA and RNA at G-C, A-T, and A-U bases, and no interaction with the backbone PO2 group. As2O3-DNA and -RNA adducts showed one type of binding with overall binding constant of K(As2O3-DNA) = 1.24 x 10(5) M(-1) and K(As2O3-RNA) = 2.60 x 10(5) M(-1). The As2O3-polynucleotide complexation is associated with a partial biopolymer aggregation and no major alterations of B-DNA or A-RNA structure.

MOLECULAR DYNAMICS SIMULATIONS OF DNA DIMERS BASED ON REPLICA-EXCHANGE UMBRELLA SAMPLING I: TEST OF SAMPLING EFFICIENCY

In order to elucidate the stacking-unstacking process of DNA dimers, we have performed molecular dynamics simulations based on replica-exchange umbrella sampling (REUS), which is one of powerful conformational sampling techniques. We studied four DNA dimers composed of the adenine and thymine bases in both the 5′ and the 3′ positions (dApdA, dApdT, dTpdA, and dTpdT). We examined the time series of the distance between the glycosidic nitrogen atoms, root-mean-square deviations from A-DNA and B-DNA, various backbone and glycosidic torsion angles, and the pseudorotation phase angles as functions of the simulation time step. All these time series imply that the present simulation has indeed sampled a very wide conformational space. The results for the backbone and glycosidic torsion angles and pseudorotation phase angles imply that B-DNA structures are the dominant motif of the stacked dimers, while a small population of A-DNA also exists in the stacked states.

Effects of sequence and structure in the separation of nucleic acids using ion pair reverse phase liquid chromatography

Ion pair reverse phase high performance liquid chromatography on non-porous alkylated poly(styrene-divinylbenzene) particles enables the high resolution separation of double stranded DNA fragments. To further understand the separation mechanisms involved in ion pair reverse phase liquid chromatography we have analysed the effects of curved or “bent” DNA fragments with respect to their separation using both gel electrophoresis and ion pair reverse phase liquid chromatography. Size dependent separations of curved DNA fragments that migrate anomalously during gel electrophoresis were observed using ion pair reverse phase liquid chromatography. To further study the sequence effect and resulting changes in hydrophobicity of the duplex DNA, PCR fragments were generated that contain uracil in place of thymine. The resulting fragments were shown to elute with shorter retention times, demonstrating that sequence-specific effects can alter the retention of duplex DNA. The study was extended to the investigation of non-canonical B-DNA structures (Holliday junctions) under various chromatographic conditions, demonstrating that the coaxial stacking of the helices in such structures, in the presence of magnesium causes a change in retention.

Structure of B-DNA with Cations Tethered in the Major Groove,

Here, we describe the 1.6-A X-ray structure of the DDD (Dickerson-Drew dodecamer), which has been covalently modified by the tethering of four cationic charges. This modified version of the DDD, called here the DDD(4+), is composed of [d(CGCGAAXXCGCG)](2), where X is effectively a thymine residue linked at the 5 position to an n-propyl-amine. The structure was determined from crystals soaked with thallium(I), which has been broadly used as a mimic of K(+) in X-ray diffraction experiments aimed at determining positions of cations adjacent to nucleic acids. Three of the tethered cations are directed radially out from the DNA. The radially directed tethered cations do not appear to induce structural changes or to displace counterions. One of the tethered cations is directed in the 3' direction, toward a phosphate group near one end of the duplex. This tethered cation appears to interact electrostatically with the DNA. This interaction is accompanied by changes in helical parameters rise, roll, and twist and by a displacement of the backbone relative to a control oligonucleotide. In addition, these interactions appear to be associated with displacement of counterions from the major groove of the DNA.

Taxol Anticancer Activity and DNA Binding

The interaction of taxol with DNA has major biological importance since it is shown the presence of higher concentration of taxol in the nucleus, than in the human lung tumor cell. Therefore, in this report we examine the interaction of taxol with calf-thymus DNA in aqueous solution at physiological pH, using constant DNA concentration (25 or 1.25 mM phosphate) and various taxol/DNA (phosphate) ratios 1/200 to 1/2. Capillary electrophoresis and Fourier transform infrared (FTIR) difference spectroscopic methods are used to characterize the nature of drug-DNA interaction and to determine the taxol binding site, the binding constant, sequence selectivity, helix stability and biopolymer secondary structure in the taxol-DNA complexes in vitro. Structural analysis showed that taxol is an external DNA binder with no affinity towards DNA intercalation. The major target of taxol is A-T, G-C bases and the backbone PO(2) groups. Two bindings were observed for taxol-DNA complexes with K(1)= 1.4 x 10(4) M(-1) and K(2)=3.5 X 10(3) M(-1). The taxol-DNA interaction is associated with a partial helix stabilization and no major alterations of B-DNA structures.

Oligonucleotides containing 7-propynyl-7-deazaguanine: synthesis and base pair stability

Oligonucleotides incorporating the propynyl derivative of 7-deaza-2′-deoxyguanosine ( 1 ) were synthesized by solid-phase oligonucleotide synthesis. As building blocks the phosphoramidites 7a, b were prepared. The incorporation of 1 into oligonucleotides exerts a positive effect on the DNA duplex stability. The duplex stabilization by 1 was higher than that of 7-iodo-7-deaza-2′-deoxyguanosine ( 2b ). The stabilizing effect of the 7-propynyl group introduced in the 7-deazapurines is similar to that reported for 8-aza-7-deazapurines. From CD spectra it was deduced that the B-DNA structure is not significantly altered by compound 1 .

Knowledge-based prediction of DNA atomic structure from nucleic sequence.

A simple knowledge-based method for DNA atomic structure prediction from nucleic sequence is presented. We used free B-DNA crystal structures to estimate the distribution of trinucleotide base pairs and tetranucleotide base-pair steps conformational coordinates. We used these distributions as a basis to predict the 3D position of the non-hydrogen atoms of the nucleic bases of any arbitrary DNA sequence of any length. The only constraint imposed was that the structure is a B-DNA one with Watson-Crick complementary base pairs. The method was tested on not seen DNA structures with sequence lengths varying from 6bp to 12bp. The obtained predictions have RMSE around 0.5 A for the translational conformational coordinates, and around 5 degrees for the rotational. For the estimation of the nucleic base non-hydrogen atom coordinates the RMSE is around 1.1 A. The knowledge-based method outperformed a technique based on genetic algorithms in the prediction of B-DNA structures.

Taxol interaction with DNA and RNA — Stability and structural features

Taxol (paclitaxel) is an anticancer drug that interacts with microtubule proteins in a manner that catalyzes their formation from tubulin and stabilizes the resulting structures. However, in the human lung tumor cell, the concentration of paclitaxel is highest in the nucleus. Therefore, it was of interest to examine the interaction of taxol with DNA and RNA in aqueous solution at physiological pH. Capillary electrophoresis and Fourier transform infrared (FTIR) difference spectroscopic methods were used to characterize the nature of drug–DNA and drug–RNA interactions and to determine the taxol binding site, the binding constant, the sequence selectivity, the helix stability, and the biopolymer secondary structure in the taxol–polynucleotide complexes in vitro. The FTIR spectroscopic studies were conducted with taxol/polynucleotide (phosphate) ratios of 1/80, 1/40, 1/20, 1/10, 1/4, and 1/2 with a final DNA(P) or RNA(P) concentration of 12.5 mmol/L, and capillary electrophoresis was performed after incubation of taxol with polynucleotides at ratios of 1/200 to 1/12 with a final polynucleotide concentration of 1.25 mmol/L. Taxol was shown to bind to DNA and RNA at G–C, A–T, or A–U bases and the backbone PO2group. Two types of binding were observed for taxol–DNA with K1 = 1.3 × 104L mol–1and K2 = 3.5 × 103L mol–1, whereas taxol–RNA complexes showed one type of binding with K = 1.3 × 104L mol–1. The taxol–polynucleotide complexation is associated with a partial helix stabilization and no major alterations of B-DNA or A-RNA structure. Key words: DNA, RNA, taxol, binding site, binding constant, conformation, helix stability, electrophoresis, FTIR spectroscopy.

Local sequential minimization of double stranded B-DNA using Monte Carlo annealing.

A software algorithm has been developed to investigate the folding process in B-DNA structures in vacuum under a simple and accurate force field. This algorithm models linear double stranded B-DNA sequences based on a local, sequential minimization procedure. The original B-DNA structures were modeled using initial nucleotide structures taken from the Brookhaven database. The models contain information at the atomic level allowing one to investigate as accurately as possible the structure and characteristics of the resulting DNA structures. A variety of DNA sequences and sizes were investigated containing coding and non-coding, random and real, homogeneous or heterogeneous sequences in the range of 2 to 40 base pairs. The force field contains terms such as angle bend, Lennard-Jones, electrostatic interactions and hydrogen bonding which are set up using the Dreiding II force field and defined to account for the helical parameters such as twist, tilt and rise. A close comparison was made between this local minimization algorithm and a global one (previously published) in order to find out advantages and disadvantages of the different methods. From the comparison, this algorithm gives better and faster results than the previous method, allowing one to minimize larger DNA segments. DNA segments with a length of 40 bases need approximately 4 h, while 2.5 weeks are needed with the previous method. After each minimization the angles between phosphate-oxygen-carbon A1, the oxygen-phosphate-oxygen A2 and the average helical twists were calculated. From the generated fragments it was found that the bond angles are A1=150 degrees +/-2 degrees and A2=130 degrees +/-10 degrees, while the helical twist is 36.6 degrees +/-2 degrees in the A strand and A1=150 degrees +/-6 degrees and A2=130+/-6 degrees with helical twist 39.6 degrees +/-2 degrees in the B strand for the DNA segment with the same sequence as the Dickerson dodecamer.

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.

DNA fine structure and dynamics in crystals and in solution: the impact of BI/BII backbone conformations.

Sugar-phosphate backbone conformations are an important structural element for a complete understanding of specific recognition in nucleic acid-protein interactions. They can be involved both in early stages of target discrimination and in structural adaptation upon binding. In the first part of this study, we have analyzed high-resolution structures of double-stranded B-DNA either isolated or bound to proteins, and explored the impact of both the standard BI and the unusual BII phosphate backbone conformations on neighboring sugar puckers and on selected helical parameters. Correlations are found to be similar for free and bound DNA, and in both categories, the possible facing backbone conformations (BI.BI, BI.BII, and BII.BII) define well-characterized substates in the B-DNA conformational space. Notably, BII.BII steps are characterized by specific, and sequence-independent, structural effects involving reduced standard deviations for almost all conformational parameters. In the second part of this work, we analyze four 10 ns molecular dynamics simulations in explicit solvent on the DNA targets of NF-kappaB and bovine papillomavirus E2 proteins, highlighting the multiplicity of backbone dynamical behavior. These results show sequence effects on the percentages of BI and BII conformers, the preferential state of facing backbones, the occurrence of coupled transitions. The backbone states can consequently be seen as a mechanism for transmitting information from the bases to the phosphate groups and thus for modulating the overall structural properties of the target DNA.

Full-strand Minimization of Single and Double Stranded B-DNA Using Monte-Carlo Annealing

A new algorithm is developed which implements the Monte Carlo annealing method and is used to model and minimize linear, single and double stranded B-DNA sequences. The preliminary B-DNA structures are modeled using initial structures obtained from the Brookhaven database. The model contains structural details at the atomic level and is therefore more elaborate and accurate than the pseudo-atomic and elastomechanical models. The minimization concerns the entire chain length and not only local nucleotide complexes. A variety of DNA sequences (coding or non-coding, random or real, homogeneous or heterogeneous) are investigated in the range of 20-40 base pairs. The potential energy function is written in terms of a set of internal coordinates defined to account for the helical parameters such as twist, tilt and rise, which are important parameters for the description of the global shape of any type of DNA or RNA molecule. The force field used is composed of a limited number of bonded and non-bonded interactions such as bond stretch, angle bend and Lennard-Jones interactions with the Dreiding II force field parameter set used for these interactions. From the minimized structures the angles between Phosphate-Oxygen-Carbon " A 1 " and Oxygen-Phosphate-Oxygen " A 2 " and the average helical twist were calculated. For single strands it is shown that the bond angles are A 1 =107 - 1° and A 2 =122 - 1°, while the helical twist is 37.8 - 1°. For double stranded DNA our model predicts the helical twist of h =35.5 - 2° well, in the A strand, while the prediction is less accurate, h =47 - 2° for the complementary strand B. The average values for the angles A 1 and A 2 are 130 - 1° and 150 - 3° for strand A and 102 - 4° and 123 - 5° for strand B. The reason for this discrepancy is attributed to the different conformations initially adopted by the sugars in the A and B strands.

Electrostatic energy analysis of 8-oxoguanine DNA lesion—molecular dynamics study

One nanosecond molecular dynamics (MD) simulation was performed for two DNA segments each composed of 30 base pairs. In one DNA segment the native guanines at nucleotides positions 17 and 19 were replaced with two 8-oxoguanines (8-oxoG) (8-oxoG is mutagenic DNA oxo-lesion). The analysis of results was focused on the electrostatic energy that is supposed to be significant factor causing the disruption of DNA base stacking in DNA duplex and may also serve as a signal toward the repair enzyme informing the presence of the lesion. The repulsive interaction between 8-oxoG and the entire DNA molecule was observed, which caused the extrahelical position of 8-oxoG (position 19). The repulsive electrostatic interaction between both 8-oxoG lesions contributed to the flipping out of one 8-oxoG and to the local instability of the lesioned DNA region. The electrostatic potential at the surface of DNA close to the lesions has more negative value than the same region on the native DNA. This electrostatic potential may signal presence of the lesion to the repair enzyme. In the simulation of native DNA segment, no significant structural changes were observed and B-DNA structure was well preserved throughout the MD simulation.