Chiral phosphorothioate analogues of B-DNA. The crystal structure of Rp-d[Gp(S)CpGp(S)CpGp(S)C]

Structural dynamics of therapeutic nucleic acids with phosphorothioate backbone modifications

Antisense oligonucleotides (ASOs) offer ground-breaking possibilities for selective pharmacological intervention for any gene product-related disease. Therapeutic ASOs contain extensive chemical modifications that improve stability to enzymatic cleavage and modulate binding affinity relative to natural RNA/DNA. Molecular dynamics (MD) simulation can provide valuable insights into how such modifications affect ASO conformational sampling and target binding. However, force field parameters for chemically modified nucleic acids (NAs) are still underdeveloped. To bridge this gap, we developed parameters to allow simulations of ASOs with the widely applied phosphorothioate (PS) backbone modification, and validated these in extensive all-atom MD simulations of relevant PS-modified NA systems representing B-DNA, RNA, and DNA/RNA hybrid duplex structures. Compared to the corresponding natural NAs, single PS substitutions had marginal effects on the ordered DNA/RNA duplex, whereas substantial effects of phosphorothioation were observed in single-stranded RNA and B-DNA, corroborated by the experimentally derived structure data. We find that PS-modified NAs shift between high and low twist states, which could affect target recognition and protein interactions for phosphorothioated oligonucleotides. Furthermore, conformational sampling was markedly altered in the PS-modified ssRNA system compared to that of the natural oligonucleotide, indicating sequence-dependent effects on conformational preference that may in turn influence duplex formation.

Nucleic acid-protein interfaces studied by MAS solid-state NMR spectroscopy

Solid-state NMR (ssNMR) has become a well-established technique to study large and insoluble protein assemblies. However, its application to nucleic acid-protein complexes has remained scarce, mainly due to the challenges presented by overlapping nucleic acid signals. In the past decade, several efforts have led to the first structure determination of an RNA molecule by ssNMR. With the establishment of these tools, it has become possible to address the problem of structure determination of nucleic acid-protein complexes by ssNMR. Here we review first and more recent ssNMR methodologies that study nucleic acid-protein interfaces by means of chemical shift and peak intensity perturbations, direct distance measurements and paramagnetic effects. At the end, we review the first structure of an RNA-protein complex that has been determined from ssNMR-derived intermolecular restraints.

From Angstroms to Nanometers: Measuring Interatomic Distances by Solid-State NMR

Internuclear distances represent one of the main structural constraints in molecular structure determination using solid-state NMR spectroscopy, complementing chemical shifts and orientational restraints. Although a large number of magic-angle-spinning (MAS) NMR techniques have been available for distance measurements, traditional ¹³C and ¹⁵N NMR experiments are inherently limited to distances of a few angstroms due to the low gyromagnetic ratios of these nuclei. Recent development of fast MAS triple-resonance ¹⁹F and ¹H NMR probes has stimulated the design of MAS NMR experiments that measure distances in the 1-2 nm range with high sensitivity. This review describes the principles and applications of these multiplexed multidimensional correlation distance NMR experiments, with an emphasis on ¹⁹F- and ¹H-based distance experiments. Representative applications of these long-distance NMR methods to biological macromolecules as well as small molecules are reviewed.

Sulfur-centered hemi-bond radicals as active intermediates in S-DNA phosphorothioate oxidation

Phosphorothioate (PS) modifications naturally appear in bacteria and archaea genome and are widely used as antisense strategy in gene therapy. But the chemical effects of PS introduction as a redox active site into DNA (S-DNA) is still poorly understood. Herein, we perform time-resolved spectroscopy to examine the underlying mechanisms and dynamics of the PS oxidation by potent radicals in free model, in dinucleotide, and in S-oligomer. The crucial sulphur-centered hemi-bonded intermediates -P–S∴S–P- were observed and found to play critical roles leading to the stable adducts of -P–S–S–P-, which are backbone DNA lesion products. Moreover, the oxidation of the PS moiety in dinucleotides d[G_PSG], d[A_PSA], d[G_PSA], d[A_PSG] and in S-oligomers was monitored in real-time, showing that PS oxidation can compete with adenine but not with guanine. Significantly, hole transfer process from A^+• to PS and concomitant -P–S∴S–P- formation was observed, demonstrating the base-to-backbone hole transfer unique to S-DNA, which is different from the normally adopted backbone-to-base hole transfer in native DNA. These findings reveal the distinct backbone lesion pathway brought by the PS modification and also imply an alternative -P–S∴S–P-/-P–S–S–P- pathway accounting for the interesting protective role of PS as an oxidation sacrifice in bacterial genome.

Symmetry in Nucleic-Acid Double Helices

In nature and in the test tube, nucleic acids occur in many different forms. Apart from single-stranded, coiled molecules, DNA and RNA prefer to form helical arrangements, in which the bases are stacked to shield their hydrophobic surfaces and expose their polar edges. Focusing on double helices, we describe the crucial role played by symmetry in shaping DNA and RNA structure. The base pairs in nucleic-acid double helices display rotational pseudo-symmetry. In the Watson–Crick base pairs found in naturally occurring DNA and RNA duplexes, the symmetry axis lies in the base-pair plane, giving rise to two different helical grooves. In contrast, anti-Watson–Crick base pairs have a dyad axis perpendicular to the base-pair plane and identical grooves. In combination with the base-pair symmetry, the syn/anti conformation of paired nucleotides determines the parallel or antiparallel strand orientation of double helices. DNA and RNA duplexes in nature are exclusively antiparallel. Watson–Crick base-paired DNA or RNA helices display either right-handed or left-handed helical (pseudo-) symmetry. Genomic DNA is usually in the right-handed B-form, and RNA double helices adopt the right-handed A-conformation. Finally, there is a higher level of helical symmetry in superhelical DNA in which B-form double strands are intertwined in a right- or left-handed sense.

Submolecular Structure and Orientation of Oligonucleotide Duplexes Tethered to Gold Electrodes Probed by Infrared Reflection Absorption Spectroscopy: Effect of the Electrode Potentials

Unique electronic and ligand recognition properties of the DNA double helix provide basis for DNA applications in biomolecular electronic and biosensor devices. However, the relation between the structure of DNA at electrified interfaces and its electronic properties is still not well understood. Here, potential-driven changes in the submolecular structure of DNA double helices composed of either adenine-thymine (dAdT)₂₅ or cytosine-guanine (dGdC)₂₀ base pairs tethered to the gold electrodes are for the first time analyzed by in situ polarization modulation infrared reflection absorption spectroscopy (PM IRRAS) performed under the electrochemical control. It is shown that the conformation of the DNA duplexes tethered to gold electrodes via the C₆ alkanethiol linker strongly depends on the nucleic acid sequence composition. The tilt of purine and pyrimidine rings of the complementary base pairs (dAdT and dGdC) depends on the potential applied to the electrode. By contrast, neither the conformation nor orientation of the ionic in character phosphate-sugar backbone is affected by the electrode potentials. At potentials more positive than the potential of zero charge (pzc), a gradual tilting of the double helix is observed. In this tilted orientation, the planes of the complementary purine and pyrimidine rings lie ideally parallel to each other. These potentials do not affect the integral stability of the DNA double helix at the charged interface. At potentials more negative than the pzc, DNA helices adopt a vertical to the gold surface orientation. Tilt of the purine and pyrimidine rings depends on the composition of the double helix. In monolayers composed of (dAdT)₂₅ molecules the rings of the complementary base pairs lie parallel to each other. By contrast, the tilt of purine and pyrimidine rings in (dGdC)₂₀ helices depends on the potential applied to the electrode. Such potential-induced mobility of the complementary base pairs can destabilize the helix structure at a submolecular level. These pioneer results on the potential-driven changes in the submolecular structure of double stranded DNA adsorbed on conductive supports contribute to further understanding of the potential-driven sequence-specific electronic properties of surface-tethered oligonucleotides.

Phosphorothioate DNA Stabilized Fluorescent Gold and Silver Nanoclusters

Unmodified single-stranded DNA has recently gained popularity for the templated synthesis of fluorescent noble metal nanoclusters (NCs). Bright, stable, and biocompatible clusters have been developed primarily through optimization of DNA sequence. However, DNA backbone modifications have not yet been investigated. In this work, phosphorothioate (PS) DNAs are evaluated in the synthesis of Au and Ag nanoclusters, and are employed to successfully template a novel emitter using T₁₅ DNA at neutral pH. Mechanistic studies indicate a distinct UV-dependent formation mechanism that does not occur through the previously reported thymine N3. The positions of PS substitution have been optimized. This is the first reported use of a T₁₅ template at physiological pH for AgNCs.

Theoretical study on the relationship between Rp-phosphorothioation and base-step in S-DNA: based on energetic and structural analysis

Phosphorothioation (PT), previously used in synthetic antisense drugs to arrest the transcription or translation process, is also a novel physiological modification in bacteria DNAs. In the previous study, we reported that Rp-phosphorothioation (Rp-PT) destabilizes B-type helix significantly, using a quantum-mechanics-based energy scoring function developed with a dinucleotide model ( Zhang et al. J. Phys. Chem. B , 2012 , 116 , 10639 - 10648 ). A consequent question surfaces in the field of the phosphorothioated DNA (S-DNA) research: does the endogenous chemical modification interact with the base sequence in the bacterial genomes, e.g., in terms of the most common structure of the B-type helix? In this work, we carried out further energetic analysis on the backbone relative energies calculated with the scoring function according to 16 groups of base-step classifications. Moreover, we conducted molecular dynamics simulations of the B-helical structure with the different base-pair steps, to investigate the detailed structural changes upon the O-/S-substitution. As a result, the Rp-PT modification definitively enhances the stiffness of the backbone and differentiates backbone stability as an interaction with base-steps. Furthermore, certain exceptional sequences such as GT and CC were highlighted in the structural analysis of the sulfur local contacts and relative orientation of double strands, indicating that Rp-PT can cross-talk with particular base-steps. The special effects between the phosphorothioation and base-step may be related to the conservative consensus observed highly frequently in bacterial genomes.

Formation of S-Cl phosphorothioate adduct radicals in dsDNA S-oligomers: hole transfer to guanine vs disulfide anion radical formation

In phosphorothioate-containing dsDNA oligomers (S-oligomers), one of the two nonbridging oxygen atoms in the phosphate moiety of the sugar-phosphate backbone is replaced by sulfur. In this work, electron spin resonance (ESR) studies of one-electron oxidation of several S-oligomers by Cl2(•-) at low temperatures are performed. Electrophilic addition of Cl2(•-) to phosphorothioate with elimination of Cl(-) leads to the formation of a two-center three-electron σ(2)σ*(1)-bonded adduct radical (-P-S-̇Cl). In AT S-oligomers with multiple phosphorothioates, i.e., d[ATATAsTsAsT]2, -P-S-̇Cl reacts with a neighboring phosphorothioate to form the σ(2)σ*(1)-bonded disulfide anion radical ([-P-S-̇S-P-](-)). With AT S-oligomers with a single phosphorothioate, i.e., d[ATTTAsAAT]2, reduced levels of conversion of -P-S-̇Cl to [-P-S-̇S-P-](-) are found. For guanine-containing S-oligomers containing one phosphorothioate, -P-S-̇Cl results in one-electron oxidation of guanine base but not of A, C, or T, thereby leading to selective hole transfer to G. The redox potential of -P-S-̇Cl is thus higher than that of G but is lower than those of A, C, and T. Spectral assignments to -P-S-̇Cl and [-P-S-̇S-P-](-) are based on reaction of Cl2(•-) with the model compound diisopropyl phosphorothioate. The results found for d[TGCGsCsGCGCA]2 suggest that [-P-S-̇S-P-](-) undergoes electron transfer to the one-electron-oxidized G, healing the base but producing a cyclic disulfide-bonded backbone with a substantial bond strength (50 kcal/mol). Formation of -P-S-̇Cl and its conversion to [-P-S-̇S-P-](-) are found to be unaffected by O2, and this is supported by the theoretically calculated electron affinities and reduction potentials of [-P-S-S-P-] and O2.

Intermediate rate atomic trajectories of RNA by solid-state NMR spectroscopy

Many RNAs undergo large conformational changes in response to the binding of proteins and small molecules. However, when RNA functional dynamics occur in the nanosecond-microsecond time scale, they become invisible to traditional solution NMR relaxation methods. Residual dipolar coupling methods have revealed the presence of extensive nanosecond-microsecond domain motions in HIV-1 TAR RNA, but this technique lacks information on the rates of motions. We have used solid-state deuterium NMR to quantitatively describe trajectories of key residues in TAR by exploiting the sensitivity of this technique to motions that occur in the nanosecond-microsecond regime. Deuterium line shape and relaxation data were used to model motions of residues within the TAR binding interface. The resulting motional models indicate two functionally essential bases within the single-stranded bulge sample both the free and Tat-bound conformations on the microsecond time scale in the complete absence of the protein. Thus, our results strongly support a conformational capture mechanism for recognition: the protein does not induce a new RNA structure, but instead captures an already-populated conformation.

Crystallographic studies of chemically modified nucleic acids: a backward glance

Chemically modified nucleic acids (CNAs) are widely explored as antisense oligonucleotide or small interfering RNA (siRNA) candidates for therapeutic applications. CNAs are also of interest in diagnostics, high-throughput genomics and target validation, nanotechnology and as model systems in investigations directed at a better understanding of the etiology of nucleic acid structure, as well as the physicochemical and pairing properties of DNA and RNA, and for probing protein-nucleic acid interactions. In this article, we review research conducted in our laboratory over the past two decades with a focus on crystal-structure analyses of CNAs and artificial pairing systems. We highlight key insights into issues ranging from conformational distortions as a consequence of modification to the modulation of pairing strength, and RNA affinity by stereoelectronic effects and hydration. Although crystal structures have only been determined for a subset of the large number of modifications that were synthesized and analyzed in the oligonucleotide context to date, they have yielded guiding principles for the design of new analogs with tailor-made properties, including pairing specificity, nuclease resistance, and cellular uptake. And, perhaps less obviously, crystallographic studies of CNAs and synthetic pairing systems have shed light on fundamental aspects of DNA and RNA structure and function that would not have been disclosed by investigations solely focused on the natural nucleic acids.

Conformations of the sugar-phosphate backbone in helical DNA crystal structures

The DNA backbone geometry was analyzed for 96 crystal structures of oligodeoxynucleotides. The ranges and mean values of the torsion angles for the observed subclasses of the A-, B-, and Z-DNA conformational types were determined by analyzing distributions of the torsion angles and scattergrams relating pairs of angles.

Hydration dependent dynamics in RNA

The essential role played by local and collective motions in RNA function has led to a growing interest in the characterization of RNA dynamics. Recent investigations have revealed that even relatively simple RNAs experience complex motions over multiple time scales covering the entire ms-ps motional range. In this work, we use deuterium solid-state NMR to systematically investigate motions in HIV-1 TAR RNA as a function of hydration. We probe dynamics at three uridine residues in different structural environments ranging from helical to completely unrestrained. We observe distinct and substantial changes in (2)H solid-state relaxation times and lineshapes at each site as hydration levels increase. By comparing solid-state and solution state (13)C relaxation measurements, we establish that ns-micros motions that may be indicative of collective dynamics suddenly arise in the RNA as hydration reaches a critical point coincident with the onset of bulk hydration. Beyond that point, we observe smaller changes in relaxation rates and lineshapes in these highly hydrated solid samples, compared to the dramatic activation of motion occurring at moderate hydration.

Sequence preference for BI/BII conformations in DNA: MD and crystal structure data analysis

Deciphering sequence information from sugar-phosphate backbone is finely tuned through the conformational substates of DNA. BII conformation, one of the conformational substates of B-DNA, is known to play a key role in DNA-protein recognition. BI and BII are identified by the epsilon-zeta difference, which is negative in BI and positive in BII. Our analysis of MD and crystal structures shows that BII conformation is sequence specific and dinucleotides GC, CG, CA, TG, TA show high preference to take up BII conformation, while TT, TC, CT, CC dinucleotides rarely take up this conformation. Significant changes were observed in the dinucleotide parameters viz. twist, roll, and slide for the steps having BII conformation. Interestingly, the magnitude of variation in the dinucleotide parameters is seen to depend mainly on two factors, the magnitude of epsilon-zeta difference and the presence or absence of BII conformation in the second strand, across the WC base-paired dinucleotide step. Based on these two factors, the conformational substate of a dinucleotide step can be further classified as BI.BI (BI conformation in both strands), BI.BII (BI conformation in one strand and BII conformation in the other), and BII.BII (BII conformation in both strands). The occurrence of BII in both strands was found to be quite rare and thus, it can be concluded that BI.BI and BI.BII hybrid steps are more favorable than a BII.BII step. In conformity with the sequence preference seen for dinucleotides in each strand, BII.BII combination of backbone conformation was observed only for GC, CG, CA, and TG containing dinucleotide steps. We further classified BII.BII step as strong BII and weak BII depending on the magnitude of the average epsilon-zeta difference. The dinucleotide steps which belong to the category of strong BII, have large twist, high positive slide and negative roll values, while those in the weak BII group have roll, twist, and slide values similar to that of hybrid BI.BII steps. This conformational property could be contributing to the groove opening/closing and thus can modulate protein-DNA interaction.

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.

Structure of the U6 RNA intramolecular stem-loop harboring an S(P)-phosphorothioate modification

Phosphorothioate-substitution experiments are often used to elucidate functionally important metal ion-binding sites on RNA. All previous experiments with S(P)-phosphorothioate-substituted RNAs have been done in the absence of structural information for this particular diastereomer. Yeast U6 RNA contains a metal ion-binding site that is essential for spliceosome function and includes the pro-S(P) oxygen 5' of U(80). S(P)-phosphorothioate substitution at this location creates spliceosomes dependent on thiophilic ions for the first step of splicing. We have determined the solution structure of the U(80) S(P)-phosphorothioate-substituted U6 intramolecular stem-loop (ISL), and also report the refined NMR structure of the unmodified U6 ISL. Both structures were determined with inclusion of (1)H-(13)C residual dipolar couplings. The precision of the structures with and without phosphorothioate (RMSD = 1.05 and 0.79 A, respectively) allows comparison of the local and long-range structural effect of the modification. We find that the U6-ISL structure is unperturbed by the phosphorothioate. Additionally, the thermodynamic stability of the U6 ISL is dependent on the protonation state of the A(79)-C(67) wobble pair and is not affected by the adjacent phosphorothioate. These results indicate that a single S(P)-phosphorothioate substitution can be structurally benign, and further validate the metal ion rescue experiments used to identify the essential metal-binding site(s) in the spliceosome.

The conformer substates of nonoriented B-type DNA in double helical poly(dG-dC)

A nonoriented hydrated film of poly(dG-dC) with ?20 water molecules per nucleotide (called B* by Loprete and Hartman (Biochem. 32, 4077-4082 (1993)) was studied by Fourier transform infrared (FT-IR) spectroscopy either as equilibrated sample between 290 and 270 K or, after quenching into the glassy state, as nonequilibrated film isothermally at 200 and 220 K. IR spectral changes on isothermal relaxation at 200 and 220 K, caused by interconversion of two conformer substates, are revealed by difference spectra. Comparison with difference curves obtained in the same manner from two classical B-DNA forms, namely the d(CGCGAATTCGCG)(2) dodecamer and polymeric NaDNA from salmon testes, revealed that the spectral changes on B(I)-to-B(II) interconversion in the classical B-DNA forms are very similar to those in the B*-form, and that the spectroscopic differences between the B(I) and B(II) features from classical B-DNA and those from the modified B*-form are minor. Nonexponential kinetics of the B(I)-->B(II) transition in the B*-form of poly(dG-dC) at 200 K showed that the structural relaxation time is about three times of that in the classical B-DNA forms (approximately equal to 30 versus approximately equal to 10 min at 200 K). The unexpected reversal of conformer substates interconversion (that is B(II)-->B(I) transition on cooling from 290 K and B(I)-->B(II) transition on isothermal relaxation at 200 K) observed for classical B-DNA occurs also in the modified B*-form. We therefore conclude that restructuring of hydration shells rules the low-temperature dynamics of the B*-form via its two conformer substates in the same manner reported for classical B-DNA by Pichler et al. (J. Phys. Chem. B 106, 3263-3274 (2002)).

Crystal studies of B-DNA: the answers and the questions

The history of DNA crystallography is reviewed and is followed by discussion of the methods used for structure determination. The features of B-DNA molecular and crystal structures are described.

Synthesis and NMR characterization of diastereomeric CPSMeG derivatives

Synthesis and stereochemical characterization of enantiomerically pure CPSMeG derivatives by NMR methods are reported. The effect of 5'-dimethoxytrityl on the conformational properties is described. It was found that in P-diastereomers the conformational differences about the C3'-O3' bond, as discernible from the delta J values, are enhanced by the presence of this protecting group.

Phosphorothioate oligodeoxynucleotides: what is their origin and what is unique about them?

The development of nucleoside phosphorothioates is described in its historical context. Examples of the interaction of phosphorothioate groups, present either in oligodeoxynucleotides or in DNA, with nucleases are presented. The structural features responsible for the resistance of the phosphorothioates toward degradation by nucleases are discussed, as are the possible reasons for the high-affinity interaction of phosphorothioate oligodeoxynucleotides with certain proteins.

Sequence effects on energetic and structural properties of phosphorothioate DNA: a molecular modelling study

Phosphorothioate (PS) oligonucleotides constitute a new class of potent drugs, resulting from the replacement of one anionic oxygen of the phosphodiester backbone by one sulphur atom. This replacement confers chirality to the phosphorus atom (PSS or PSR) and alters the energetic, structural and biological properties of B-DNA. These properties were assessed by molecular mechanics calculations on a set of regular sequences, d(YR)8.d(YR)8 and d(RR)8.d(YY)8 (R, purine; Y, pyrimidine). Results indicated: (i) destabilisation of both the PS(R)and the PSS oligomers, the loss of total energy being mainly due to a variation in the electrostatic term; (ii) an additional chirality effect, due to van der Waals and backbone angle energies, larger for PSS oligomers than for PSR oligomers; (iii) a clear sequence effect on stability, particularly from the base immediately preceding the PS group. Even though the PS group alters the stability of oligomers, it does not significantly modify the conformation. Altogether, our molecular modelling data parallel the available experimental data. Our results reveal that sequence effects on the energetic properties of PS oligomers are local and additive. Therefore, studies of the set of the 10 unique double-stranded modified dinucleotide steps included in regular oligomers could be used to predict the behaviour of any double-stranded PS-DNA.

Unexpected BII conformer substate population in unoriented hydrated films of the d(CGCGAATTCGCG)2 dodecamer and of native B-DNA from salmon testes

Conformational substates of B-DNA had been observed so far in synthetic oligonucleotides but not in naturally occurring highly polymeric B-DNA. Our low-temperature experiments show that native B-DNA from salmon testes and the d(CGCGAATTCGCG)2 dodecamer have the same BI and BII substates. Nonequilibrium distribution of conformer population was generated by quenching hydrated unoriented films to 200 K, and isothermal structural relaxation toward equilibrium by interconversion of substates was followed by Fourier transform infrared spectroscopy. BI interconverts into BII on isothermal relaxation at 200 K, whereas on slow cooling from ambient temperature, BII interconverts into BI. Our estimation of the dodecamer's BI-to-BII conformer substate population by curve resolution of the symmetrical stretching vibration of the ionic phosphate is 2.4 +/- 0.5 to 1 at 200 K, and it is 1.3 +/- 0.5 to 1 between 270 and 290 K. Pronounced spectral changes upon BI-to-BII interconversion are consistent with base destacking coupled with migration of water from ionic phosphate toward the phosphodiester and sugar moieties. Nonspecific interaction of proteins with the DNA backbone could become specific by induced-fit-type interactions with either BI or BII backbone conformations. This suggests that the BI-to-BII substate interconversion could be a major contributor to the protein recognition process.

Molecular dynamic simulations of environment and sequence dependent DNA conformations: the development of the BMS nucleic acid force field and comparison with experimental results

Molecular dynamic (MD) simulations using the BMS nucleic acid force field produce environment and sequence dependent DNA conformations that closely mimic experimentally derived structures. The parameters were initially developed to reproduce the potential energy surface, as defined by quantum mechanics, for a set of small molecules that can be used as the building blocks for nucleic acid macromolecules (dimethyl phosphate, cyclopentane, tetrahydrofuran, etc.). Then the dihedral parameters were fine tuned using a series of condensed phase MD simulations of DNA and RNA (in zero added salt, 4M NaCl, and 75% ethanol solutions). In the tuning process the free energy surface for each dihedral was derived from the MD ensemble and fitted to the conformational distributions and populations observed in 87 A- and B-DNA x-ray and 17 B-DNA NMR structures. Over 41 nanoseconds of MD simulations are presented which demonstrate that the force field is capable of producing stable trajectories, in the correct environments, of A-DNA, double stranded A-form RNA, B-DNA, Z-DNA, and a netropsin-DNA complex that closely reproduce the experimentally determined and/or canonical DNA conformations. Frequently the MD averaged structure is closer to the experimentally determined structure than to the canonical DNA conformation. MD simulations of A- to B- and B- to A-DNA transitions are also shown. A-DNA simulations in a low salt environment cleanly convert into the B-DNA conformation and converge into the RMS space sampled by a low salt simulation of the same sequence starting from B-DNA. In MD simulations using the BMS force field the B-form of d(GGGCCC)2 in a 75% ethanol solution converts into the A-form. Using the same methodology, parameters, and conditions the A-form of d(AAATTT)2 correctly converts into the B-DNA conformation. These studies demonstrate that the force field is capable of reproducing both environment and sequence dependent DNA structures. The 41 nanoseconds (nsec) of MD simulations presented in this paper paint a global picture which suggests that the DNA structures observed in low salt solutions are largely due to the favorable internal energy brought about by the nearly uniform screening of the DNA electrostatics. While the conformations sampled in high salt or mixed solvent environments occur from selective and asymmetric screening of the phosphate groups and DNA grooves, respectively, brought about by sequence induced ion and solvent packing.

Hydration of the phosphate group in double-helical DNA

Water distributions around phosphate groups in 59 B-, A-, and Z-DNA crystal structures were analyzed. It is shown that the waters are concentrated in six hydration sites per phosphate and that the positions and occupancies of these sites are dependent on the conformation and type of nucleotide. The patterns of hydration that are characteristic of the backbone of the three DNA helical types can be attributed in part to the interactions of these hydration sites.

Structure of a stereoregular phosphorothioate DNA/RNA duplex

In this work, we present the first NMR solution structure of a DNA/RNA hybrid containing stereoregular Rp-phosphorothioate modifications of all DNA backbone linkages. The complex of the enzymatically synthesized phosphorothioate DNA octamer (all-Rp)-d(GCGTCAGG) and its complementary RNA r(CCUGACGC) was found to adopt an overall conformation within the A-form family. Most helical parameters and the sugar puckers of the DNA strand assume values intermediate between A- and B-form. The close structural similarity with the unmodified DNA/RNA hybrid of the same sequence may explain why both the natural and the sulfur-substituted complex can be recognized and digested by ribonuclease H.

Crystal structure of d(GCGCGCG) with 5'-overhang G residues

The crystal structure of the DNA heptamer d(GCGCGCG) has been solved at 1.65 A resolution by the molecular replacement method and refined to an R-value of 0.184 for 3598 reflections. The heptamer forms a Z-DNA d(CGCGCG)2 with 5'-overhang G residues instead of an A-DNA d(GCGCGC)2 with 3'-overhang G residues. The overhang G residues from parallel strands of two adjacent duplexes form a trans reverse Hoogsteen G x G basepair that stacks on the six Z-DNA basepairs to produce a pseudocontinuous helix. The reverse Hoogsteen G x G basepair is unusual in that the displacement of one G base relative to the other allows them to participate in a bifurcated (G1)N2 . . . N7(G8) and an enhanced (G8)C8-H . . . O6(G1) hydrogen bond, in addition to the two usual hydrogen bonds. The 5'-overhang G residues are anti and C2'-endo while the 3'-terminal G residues are syn and C2'-endo. The conformations of both G residues are different from the syn/C3'-endo for the guanosine in a standard Z-DNA. The two cobalt hexammine ions bind to the phosphate groups in both GpC and CpG steps in Z(I) and Z(II) conformations. The water structure motif is similar to the other Z-DNA structures.

Crystal structure of the self-complementary 5'-purine start decamer d(GCGCGCGCGC) in the Z-DNA conformation. I

Alternating self-complementary oligonucleotides starting with a 5'-pyrimidine usually form left-handed Z-DNA; however, with a 5'-purine start sequence they form the right-handed A-DNA. Here we report the crystal structure of the decamer d(GCGCGCGCGC) with a 5'-purine start in the Z-DNA form. The decamer crystallizes in the hexagonal space group P6(5)22, unit cell dimensions a = b = 18.08 and c = 43.10 A, with one of the following four dinucleotide diphosphates in the asymmetric unit: d(pGpC)/d(GpCp)/d(pCpG)/d(CpGp). The molecular replacement method, starting with d(pGpC) of the isomorphous Z-DNA hexamer d(araC-dG)3 without the 2'-OH group of arabinose, was used in the structure analysis. The method gave the solution only after the sugar-phosphate conformation of the GpC step was manipulated. The refinement converged to a final R value of 18.6% for 340 unique reflections in the resolution range 8.0-1.9 A. A result of the sequence alternation is the alternation in the nucleotide conformation; guanosine is C3'-endo, syn, and cytidine is C2'-endo, anti. The CpG step phosphodiester conformation is the same as ZI or ZII, whereas that of the GpC step phosphodiester is "intermediate" in the sense that zeta (O3'-P bond) is the same as ZII but alpha (P-O5' bond) is the same as ZI. The duplexes generated from the dinucleotide asymmetric unit are stacked one on top of the other in the crystal to form an infinite pseudocontinuous helix. This renders it a quasi-polymerlike structure that has assumed the Z-DNA conformation further strengthened by the long inner Z-forming stretch d(CG)4. An interesting feature of the structure is the presence of water strings in both the major and the minor grooves. In the minor groove the cytosine carbonyl oxygen atoms of the GpC and CpG steps are cross-bridged by water molecules that are not themselves hydrogen bonded but are enclosed by the water rings in the mouth of the minor groove. In the major groove three independent water molecules form a zigzagging continuous water string that runs throughout the duplex.

Crystal structure of the B-DNA hexamer d(CTCGAG): model for an A-to-B transition

The crystal structure of the B-DNA hexamer d(CTCGAG) has been solved at 1.9 A resolution by iterative single isomorphous replacement, using the brominated derivative d(CG5BrCGAG), and refined to an R-factor of 18.6% for 120 nonhydrogen nucleic acid atoms and 32 water molecules. Although the central four base pairs form a typical B-form helix, several parameters suggest a transition to an A-like conformation at the termini. Based on this observation, a B-to-A transition was modeled, maintaining efficient base stacking across the junction. The wide minor groove (approximately 6.9 A) is reminiscent of that in the side-by-side double drug-DNA complexes and hosts a double spine of hydration. The global helix axes of the pseudo-continuous helices are at an acute angle of 60 degrees. The pseudocontinuous stacking is reinforced by the minor groove water structure extending between the two duplexes. The crossover point of two pairs of stacked duplexes is at the stacking junction, unlike that observed in the B-DNA decamers and dodecamers. This arrangement may have implications for the structure of a four-way DNA junction. The duplexes are arranged around a large (approximately 20 A diameter) channel centered on a 6(2) screw axis.

Influence of the phosphate backbone on the recognition and hydrolysis of DNA by the EcoRV restriction endonuclease. A study using oligodeoxynucleotide phosphorothioates

A set of phosphorothioate-containing oligonucleotides based on pGACGATATCGTC, a self-complementary dodecamer that contains the EcoRV recognition sequence (GATATC), has been prepared. The phosphorothioate group has been individually introduced at the central nine phosphate positions and the two diastereomers produced at each site separated and purified. The Km and Vmax values found for each of these modified DNA molecules with the EcoRV restriction endonuclease have been determined and compared with those seen for the unmodified all-phosphate-containing dodecamer. This has enabled an evaluation of the roles that both of the non-esterified oxygen atoms in the individual phosphates play in DNA binding and hydrolysis by the endonuclease. The results have also been compared with crystal structures of the EcoRV endonuclease, complexed with an oligodeoxynucleotide, to allow further definition of phosphate group function during substrate binding and turnover. For further study, see the related article "Probing the Indirect Readout of the Restriction Enzyme EcoRV: Mutational Analysis of Contacts to the DNA Backbone" (Wenz, A., Jeltsch, A., and Pingoud, A. (1996) J. Biol. Chem. 271, 5565-5573).

NMR investigations of duplex stability of phosphorothioate and phosphorodithioate DNA analogues modified in both strands

Duplex formation from the self-complementary 12mer d(CGCGAATTCGCG) (Dickerson dodecamer) in which all phosphodiester linkages were replaced by phosphorothioate or phosphorodithioate linkages was studied using variable-temperature 1H and 31P NMR spectroscopy. Melting temperatures of the dodecamer, measured spectrophotometrically, showed significant decrease upon sulfur substitution (Tm 49 degrees C for the phosphorothioate and 21 degrees C for the phosphorodithioate, compared with 68 degrees C for the unmodified oligomer, in 1 M salt). Hyperchromicity observed upon melting of the dithioate was surprisingly low. NOESY spectra of the monothioate showed a cross-peak pattern characteristic for a right-handed duplex. Imino proton resonances of the duplex, shown by the mono- and the dithioate, were similar to those of the parent compound. In spite of monophasic melting curves, temperature dependence of the imino proton resonances and phosphorus resonances of the phosphorodithioate indicated heterogeneity with respect to base-pairing, compatible with the presence of a hairpin loop. Relaxation times (T1) of the imino protons in the phosphorothioate, determined by the saturation recovery method, were considerably shorter than in the unmodified oligomer. Base-pair lifetimes in the unmodified Dickerson dodecamer, determined by catalyst-dependent changes in relaxation rates of imino protons, were in the range of 2-30 ms at 20 degrees C. Strongly reduced base-pair lifetimes were found in the phosphorothioate analogue.

Nucleic acid crystallography: a view from the nucleic acid database

What are the future directions of the field of nucleic acid crystallography? Although there have been many duplex structures determined, the sample is still relatively small. This is especially true if one wants to derive enough information about the relationships between sequence and structure. Indeed, there are data for all the possible 10 dimer steps, but for some steps it is very limited. If the structural code resides in trimers or tetrad steps then there is simply not enough data to do meaningful statistical analyses. So the first direction that needs to be explored is the determination of more structures with more varied sequences. The other noticeable thing about the data is the shortness of the strands. While it is probably true that attempts to crystallize very long sequences will not meet with success, the idea of crystallizing sequences engineered to fit together via sticky ends such as has been done for the CAP-DNA complex (Schultz et al., 1990) should give data about the behavior of much longer stretches of DNA. The question of the effects of environment on the structure of DNA continues to be a very important one to address since DNA is rarely alone. The preliminary data we have analysed from the current sample shows that the conformation of some steps are very sensitive to packing type. Numerous studies of the hydration around DNA shows that there is a real synergy between the hydration structure and the base conformation. More data will allow further quantitation of these observations. RNA structure is the next very exciting frontier. The emerging structures of duplexes with internal loops, the two hammerhead ribozyme structures and the group I intron ribozyme have given us a glimpse of the complexity and elegance of this class of molecules. With the technology now in place to allow the determination of the structures of these molecules, the expectation is that now we will see a large increase in the number of these structures in the NDB.

An A-DNA triplet code: thermodynamic rules for predicting A- and B-DNA

The ability to predict macromolecular conformations from sequence and thermodynamic principles has long been coveted but generally has not been achieved. We show that differences in the hydration of DNA surfaces can be used to distinguish between sequences that form A- and B-DNA. From this, a "triplet code" of A-DNA propensities was derived as energetic rules for predicting A-DNA formation. This code correctly predicted > 90% of A- and B-DNA sequences in crystals and correlates with A-DNA formation in solution. Thus, with our previous studies on Z-DNA, we now have a single method to predict the relative stability of sequences in the three standard DNA duplex conformations.

Analysis of local helix bending in crystal structures of DNA oligonucleotides and DNA-protein complexes

Sequence-dependent bending of the helical axes in 112 oligonucleotide duplex crystal structures resident in the Nucleic Acid Database have been analyzed and compared with the use of bending dials, a computer graphics tool. Our analysis includes structures of both A and B forms of DNA and considers both uncomplexed forms of the double helix as well as those bound to drugs and proteins. The patterns in bending preferences in the crystal structures are analyzed by base pair steps, and emerging trends are noted. Analysis of the 66 B-form structures in the Nucleic Acid Database indicates that uniform trends within all pyrimidine-purine and purine-pyrimidine steps are not necessarily observed but are found particularly at CG and GC steps of dodecamers. The results support the idea that AA steps are relatively straight and that larger roll bends occur at or near the junctions of these A-tracts with their flanking sequences. The data on 16 available crystal structures of protein-DNA complexes indicate that the majority of the DNA bends induced via protein binding are sharp localized kinks. The analysis of the 30 available A-form DNA structures indicates that these structures are also bent and show a definitive preference for bending into the deep major groove over the shallow minor groove.

A knowledge-based model of DNA hydration

The aqueous hydration of DNA is an important aspect of its structure, which is of direct relevance to mechanisms of radiation damage. We have made a quantitative analysis of solvent interactions within hydrogen bonding distance of polar atoms of oligonucleotides using 12 B-DNA oligonucleotide crystal structures. The distribution of water molecules around the four bases, the sugar residues and the phosphate groups were generated and analysed both qualitatively and quantitatively. These data have then been used in a knowledge-based method to generate the likely hydration sites around a canonical B-DNA conformation in order to generate models of use in track studies of radiation damage.

Comparative studies of high resolution Z-DNA crystal structures. Part 1: Common hydration patterns of alternating dC-dG

The water structure in three crystal forms of the left-handed Z-DNA hexamer [d(CGCGCG)]2 has been analyzed. Several common motifs have been found in the first hydration shells. On the convex surface, the major groove of the left-handed conformation, water molecules bridge the guanine O-6 keto groups at GpC steps. Cytosine N-4 nitrogens of opposite strands are hydrated by tandem water molecules. At the bottom of the minor groove, a string of water molecules connects the cytosine O-2 keto groups. Across the minor groove guanine N-2 nitrogens are bridged to phosphate oxygens of cytosine and guanine residues by one or two water molecules. In contrast to the very regular geometry of the water structure around the bases, the arrangement of water molecules between phosphate groups appears to be less ordered. However, there is a strong correlation between the interphosphate distances and the number of water molecules or ions which link the phosphate groups. In all three structures various ions, such as sodium and magnesium ions, as well as the protonated amino and imino groups of the polycation spermine displace and replace water molecules in the first hydration shell. Nevertheless, the analysis reveals that numerous first hydration shell water molecules in Z-DNA crystals can be regarded as part of the DNA structure. Their positions and thermal parameters are generally independent of changes in the local crystallographic environment.

Binding of two distamycin A molecules in the minor groove of an alternating B-DNA duplex

Here we report the 1.8A X-ray structure of a 2:1 drug-DNA complex between distamycin A and an alternating B-DNA octamer duplex d(ICICICIC)2. The two distamycin A molecules are bound side by side with dyad symmetry in an antiparallel orientation in the expanded minor groove. The amides of each drug molecule are hydrogen bonded to the minor groove base atoms of only one DNA strand. The complex not only shows binding of two drug molecules, but the DNA duplex also exhibits striking low-high alternations in the helical twist angles, the sugar puckering and the phosphate conformations, providing the basis for a new model for an alternating B-DNA with a dinucleotide repeat.

Molecular and crystal structure of Sp-thymidin-3'-yl 4-thiothymidin-5'-yl methylphosphonate

The molecular structure of one diastereomer of the dinucleoside methylphosphonate Tp(Me)sT (1) has been determined by X-ray diffraction methods. The crystal asymmetric unit contains one molecule of 1 and one methanol in an orthorhombic unit cell of dimensions a = 13.241(4), b = 13.844(3), c = 14.944(7) A, space group P2(1)2(1)2(1). Both pyrimidine bases in 1 are oriented anti relative to the 2'-deoxyribose rings, and the sugar conformations are 2E and 2(3)T in the 4-thiothymidine and thymidine moieties, respectively. The deoxyribose-phosphonate backbone has an extended conformation with the bases completely unstacked and almost parallel. The absolute configuration at the phosphorus center in 1 is Sp.

Conformational sub-states in B-DNA

Theoretical studies of the sequence-dependent conformation of B-DNA have been carried out using Jumna, a helicoidal co-ordinate minimization algorithm. The results obtained for a series of six oligomers with repetitive sequences show that, with the exception of the homopolymers (dA)n.(dT)n and (dG)n.(dC)n, all sequences can adopt a variety of conformations characterized by considerable changes in helicoidal parameters and also in sugar puckers which adopt C(2')-endo (falling into 2 classes) or, in the case of pyrimidine nucleotides, O(1')-endo forms. These studies lead to an improved understanding of the role of base sequence on DNA conformation and point to a number of interesting correlations between the various structural parameters describing the double helix.