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

Helix morphology changes in B-DNA induced by spontaneous B(I)<==>B(II) substrate interconversion

Investigations of spontaneous, i.e. not forced, B-DNA's B(I)<==>B(II) substate transitions are carried out on the d(CGCGAATTCGCG)2 EcoRI dodecamer sequence using Molecular Dynamics Simulations. Analysis of the resulting transition processes with respect to the backbone angles reveals concerted changes not only for backbone angles epsilon, zeta, and beta, but also for the 5'-delta and 5'-chi angles. For alpha and delta inside the interconverting base step, a change is seen in short lived B(II) conformers. With respect to base morphology distinct changes are observed for buckle, propeller twist, shift, roll and twist, as well as x-displacement and tip. The base mainly involved in the changes is identified as the base preceding the interconverting phosphate. Altogether single B(I)<==>B(II) interconversions result only in local distortions represented by the larger spread of most parameters. Comparison of the atomic positional fluctuations derived from the simulation with those obtained from the static X-ray structure results in striking similarities.

The 2D [31P] spin-echo-difference constant-time [13C, 1H]-HMQC experiment for simultaneous determination of 3J(H3'P) and 3J(C4'P) in 13C-labeled nucleic acids and their protein complexes

A two-dimensional [31P] spin-echo-difference constant-time [13C, 1H]-HMQC experiment (2D [31P]-sedct-[13C, 1H]-HMQC) is introduced for measurements of 3J(C4'P) and 3J(H3'P) scalar couplings in large 13C-labeled nucleic acids and in DNA-protein complexes. This experiment makes use of the fact that 1H-13C multiple-quantum coherences in macromolecules relax more slowly than the corresponding 13C single-quantum coherences. 3J(C4'P) and 3J(H3'P) are related via Karplus-type functions with the phosphodiester torsion angles beta and epsilon, respectively, and their experimental assessment therefore contributes to further improved quality of NMR solution structures. Data are presented for a uniformly 13C, 15N-labeled 14-base-pair DNA duplex, both free in solution and in a 17-kDa protein-DNA complex.

Structure of the trp RNA-binding attenuation protein, TRAP, bound to RNA

The trp RNA-binding attenuation protein (TRAP) regulates expression of the tryptophan biosynthetic genes of several bacilli by binding single-stranded RNA. The binding sequence is composed of eleven triplet repeats, predominantly GAG, separated by two or three non-conserved nucleotides. Here we present the crystal structure of a complex of TRAP and a 53-base single-stranded RNA containing eleven GAG triplets, revealing that each triplet is accommodated in a binding pocket formed by beta-strands. In the complex, the RNA has an extended structure without any base-pairing and binds to the protein mostly by specific protein-base interactions. Eleven binding pockets on the circular TRAP 11-mer form a belt with a diameter of about 80 A. This simple but elegant mechanism of arresting the RNA segment by encircling it around a protein disk is applicable to both transcription, when TRAP binds the nascent RNA, and to translation, when TRAP binds the same sequence within a non-coding leader region of the messenger RNA.

Modeling helix-turn-helix protein-induced DNA bending with knowledge-based distance restraints

A crucial element of many gene functions is protein-induced DNA bending. Computer-generated models of such bending have generally been derived by using a presumed bending angle for DNA. Here we describe a knowledge-based docking strategy for modeling the structure of bent DNA recognized by a major groove-inserting alpha-helix of proteins with a helix-turn-helix (HTH) motif. The method encompasses a series of molecular mechanics and dynamics simulations and incorporates two experimentally derived distance restraints: one between the recognition helix and DNA, the other between respective sites of protein and DNA involved in chemical modification-enabled nuclease scissions. During simulation, a DNA initially placed at a distance was "steered" by these restraints to dock with the binding protein and bends. Three prototype systems of dimerized HTH DNA binding were examined: the catabolite gene activator protein (CAP), the phage 434 repressor (Rep), and the factor for inversion stimulation (Fis). For CAP-DNA and Rep-DNA, the root mean square differences between model and x-ray structures in nonhydrogen atoms of the DNA core domain were 2.5 A and 1.6 A, respectively. An experimental structure of Fis-DNA is not yet available, but the predicted asymmetrical bending and the bending angle agree with results from a recent biochemical analysis.

Crystal structure of the topoisomerase II poison 9-amino-[N-(2-dimethylamino)ethyl]acridine-4-carboxamide bound to the DNA hexanucleotide d(CGTACG)2

The structure of the complex formed between d(CGTACG)(2) and the antitumor agent 9-amino-[N-(2-dimethylamino)ethyl]acridine-4-carboxamide has been solved to a resolution of 1.6 A using X-ray crystallography. The complex crystallized in space group P6(4) with unit cell dimensions a = b = 30.2 A and c = 39.7 A, alpha = beta = 90 degrees, gamma = 120 degrees. The asymmetric unit contains a single strand of DNA, 1. 5 drug molecules, and 29 water molecules. The final structure has an overall R factor of 19.3%. A drug molecule intercalates between each of the CpG dinucleotide steps with its side chain lying in the major groove, and the protonated dimethylamino group partially occupies positions close to ( approximately 3.0 A) the N7 and O6 atoms of guanine G2. A water molecule forms bridging hydrogen bonds between the 4-carboxamide NH and the phosphate group of the same guanine. Sugar rings adopt the C2'-endo conformation except for cytosine C1 which moves to C3'-endo, thereby preventing steric collision between its C2' methylene group and the intercalated acridine ring. The intercalation cavity is opened by rotations of the main chain torsion angles alpha and gamma at guanines G2 and G6. Intercalation perturbs helix winding throughout the hexanucleotide compared to B-DNA, steps 1 and 2 being unwound by 8 degrees and 12 degrees, respectively, whereas the central TpA step is overwound by 17 degrees. An additional drug molecule, lying with the 2-fold axis in the plane of the acridine ring, is located at the end of each DNA helix, linking it to the next duplex to form a continuously stacked structure. The protonated N,N-dimethylamino group of this "end-stacked" drug hydrogen bonds to the N7 atom of guanine G6. In both drug molecules, the 4-carboxamide group is internally hydrogen bonded to the protonated N-10 atom of the acridine ring. The structure of the intercalated complex enables a rationalization of the known structure-activity relationships for inhibition of topoisomerase II activity, cytotoxicity, and DNA-binding kinetics for 9-aminoacridine-4-carboxamides.

DNA stretching and compression: large-scale simulations of double helical structures

Computer-simulated elongation and compression of A - and B -DNA structures beyond the range of thermal fluctuations provide new insights into high energy "activated" forms of DNA implicated in biochemical processes, such as recombination and transcription. All-atom potential energy studies of regular poly(dG).poly(dC) and poly(dA).poly(dT) double helices, stretched from compressed states of 2.0 A per base-pair step to highly extended forms of 7.0 A per residue, uncover four different hyperfamilies of right-handed structures that differ in mutual base-pair orientation and sugar-phosphate backbone conformation. The optimized structures embrace all currently known right-handed forms of double-helical DNA identified in single crystals as well as non-canonical forms, such as the original "Watson-Crick" duplex with trans conformations about the P-O5' and C5'-C4' backbone bonds. The lowest energy minima correspond to canonical A and B -form duplexes. The calculations further reveal a number of unusual helical conformations that are energetically disfavored under equilibrium conditions but become favored when DNA is highly stretched or compressed. The variation of potential energy versus stretching provides a detailed picture of dramatic conformational changes that accompany the transitions between various families of double-helical forms. In particular, the interchanges between extended canonical and non-canonical states are reminiscent of the cooperative transitions identified by direct stretching experiments. The large-scale, concerted changes in base-pair inclination, brought about by changes in backbone and glycosyl torsion angles, could easily give rise to the observed sharp increase in force required to stretch single DNA molecules more than 1.6-1.65 times their canonical extension. Our extended duplexes also help to tie together a number of previously known structural features of the RecA-DNA complex and offer a self-consistent stereochemical model for the single-stranded/duplex DNA recognition brought in register by recombination proteins. The compression of model duplexes, by contrast, yields non-canonical structures resembling the deformed steps in crystal complexes of DNA with the TATA-box binding protein (TBP). The crystalline TBP-bound DNA steps follow the calculated compression-elongation pattern of an unusual "vertical" duplex with base planes highly inclined with respect to the helical axis, exposed into the minor groove, and accordingly accessible for recognition.Significantly, the double helix can be stretched by a factor of two and compressed roughly in half before its computed internal energy rises sharply. The energy profiles show that DNA extension-compression is related not only to the variation of base-pair Rise but also to concerted changes of Twist, Roll, and Slide. We suggest that the high energy "activated" forms calculated here are critical for DNA processing, e.g. nucleo-protein recognition, DNA/RNA synthesis, and strand exchange.

The mechanics of minor groove width variation in DNA, and its implications for the accommodation of ligands

In duplex DNA, groove width and depth are salient structural features that may influence the binding of drugs and proteins. These features are affected by movement of the bases, which for example may enforce groove compression or expansion through a rolling action of the adjacent base-pairs. Moreover, the sugar-phosphate backbone can also undergo limited movement, independently of the bases, which will affect the groove shape. We have examined how the movement of the sugar-phosphate backbone may affect the minor groove width for a fixed base geometry. In agreement with earlier studies, the sugar-phosphate backbone is found to have a certain degree of conformational flexibility in A and B-like helices, and we note a comparable freedom even in the highly curved TATA element of the TATA-binding protein/DNA complex. Phosphate mobility is highly anisotropic in all cases with favoured directions that can significantly change the groove width, independent of any changes in base geometry. We describe how the movement of the sugar-phosphate backbone may affect the accommodation of drugs and proteins in the minor groove, and we present a co-ordinate scheme which emphasises the groove adjustments associated with ligand binding. The observations have implications for the related problem of how cognate molecules are accommodated in the major groove.

Protein-DNA interactions: A structural analysis

A detailed analysis of the DNA-binding sites of 26 proteins is presented using data from the Nucleic Acid Database (NDB) and the Protein Data Bank (PDB). Chemical and physical properties of the protein-DNA interface, such as polarity, size, shape, and packing, were analysed. The DNA-binding sites shared common features, comprising many discontinuous sequence segments forming hydrophilic surfaces capable of direct and water-mediated hydrogen bonds. These interface sites were compared to those of protein-protein binding sites, revealing them to be more polar, with many more intermolecular hydrogen bonds and buried water molecules than the protein-protein interface sites. By looking at the number and positioning of protein residue-DNA base interactions in a series of interaction footprints, three modes of DNA binding were identified (single-headed, double-headed and enveloping). Six of the eight enzymes in the data set bound in the enveloping mode, with the protein presenting a large interface area effectively wrapped around the DNA.A comparison of structural parameters of the DNA revealed that some values for the bound DNA (including twist, slide and roll) were intermediate of those observed for the unbound B-DNA and A-DNA. The distortion of bound DNA was evaluated by calculating a root-mean-square deviation on fitting to a canonical B-DNA structure. Major distortions were commonly caused by specific kinks in the DNA sequence, some resulting in the overall bending of the helix. The helix bending affected the dimensions of the grooves in the DNA, allowing the binding of protein elements that would otherwise be unable to make contact. From this structural analysis a preliminary set of rules that govern the bending of the DNA in protein-DNA complexes, are proposed.

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.

A hairpin conformation for the 3' overhang of Oxytricha nova telomeric DNA

The solution secondary structure of the Oxytricha nova telomeric 3' overhang, d(T4G4)2, has been investigated by Raman spectroscopy, hydrogen-deuterium exchange kinetics and gel electrophoresis. The electrophoretic mobility of d(T4G4)2 in non-denaturing gels indicates a highly compact conformation, consistent with a hairpin secondary structure. Raman markers show that the d(T4G4)2 hairpin contains equal numbers of C2'-endo/syn and C2'-endo/anti deoxyguanosine conformers, as well as G.G base-pairs of the Hoogsteen type. The hydrogen-deuterium exchange kinetics of d(T4G4)2, monitored by time-resolved Raman spectroscopy, reveal two kinetically distinct classes of guanine imino (N1H) protons. The more slowly exchanging fraction (kN1H(1)=4.6x10(-3) min-1), which represents 50% of N1H groups, is attributed to Hoogsteen-paired residues. The more rapidly exchanging fraction (kN1H(2)>/=0.3 min-1) is attributable to solvent-exposed residues. Raman dynamic probe of the kinetics of guanine C8H-->C8(2)H exchange in d(T4G4)2 reveals modest retardation vis-à-vis dGMP, which rules out quadruplex formation by the telomeric repeat and confirms an ordered secondary structure consistent with a Hoogsteen-paired hairpin. Similar Raman, hydrogen-isotope exchange and electrophoretic mobility experiments on the related telomeric model, dT6(T4G4)2, also reveal a hairpin stabilized by Hoogsteen G.G pairs. Presence of the 5' thymidine tail preceding the Oxytricha telomeric repeat has no apparent effect on the hairpin secondary structure. We propose a molecular model for the hairpin conformation of the Oxytricha nova telomeric repeat and consider its possible roles in mechanisms of telomeric DNA interaction in vitro and telomere function in vivo.

Sequence-dependent DNA structure: the role of the sugar-phosphate backbone

A detailed analysis of the coupling between the conformational properties of the sugar-phosphate backbone and the base stacking interactions in dinucleotide steps of double helical DNA is described. In X-ray crystal structures of oligonucleotides, the backbone shows one major degree of freedom, consisting of the torsion angles chi, delta, zeta and the pseudorotation phase angle, P. The remaining torsion angles (beta, epsilon, alpha and gamma) comprise two less important degrees of freedom. The base stacking interactions show three degrees of freedom: slide-roll-twist, shift-tilt, and rise (which is more or less constant). Coupling is observed between the base and backbone degrees of freedom. The major base stacking mode, slide-roll-twist, is coupled to the major backbone mode, chi-P-delta-zeta. The secondary base stacking mode, shift-tilt, is coupled to epsilon and zeta and to a lesser extent to the chi-P-delta-zeta mode. We show that the length of the backbone, C, given by the same strand C1'-C1' separation, is an excellent single parameter descriptor for the conformation of the backbone and the way in which it is coupled to the base stacking geometry. The slide-roll-twist motion relates to changes in the mean backbone length, C, and the shift-tilt motion to the difference between the lengths of the two backbone strands, DeltaC. We use this observation to develop a simple virtual bond model which describes the coupling of the backbone conformations and the base stacking geometry. A semi-flexible bond is used to connect the same strand C1'-C1' atoms. Analysis of the X-ray crystal structure database, simple geometric considerations and model building experiments all show that this bond is flexible with respect to slide, shift and propeller but rigid with respect to the other 14 local base stacking parameters. Using this simple model for the backbone in conjunction with potential energy calculations of the base stacking interactions, we show that it is possible to predict accurately the values of these 14 base step parameters, given values of slide, shift and propeller. We also show that the base step parameters fall into three distinct groups: roll, tilt and rise are determined solely by the base stacking interactions and are independent of the backbone; twist is insensitive to the base stacking interactions and is determined solely by the constraints of a relatively rigid fixed length backbone; slide and shift are the primary degrees of freedom and cannot be predicted accurately at the dinucleotide level because they are influenced by the conformations of neighbouring steps in a sequence. We have found that the context effect on slide is mediated by the chi torsion angles while the context effect on shift results from a BI to BII transition in the backbone. We have therefore reduced the dimensionality of the dinucleotide step problem to two parameters, slide and shift.

Structural equilibrium of DNA represented with different force fields

We have recently indicated preliminary evidence of different equilibrium average structures with the CHARMM and AMBER force fields in explicit solvent molecular dynamics simulations on the DNA duplex d(C5T5) . d(A5G5) (Feig, M. and B.M. Pettitt, 1997, Experiment vs. Force Fields: DNA conformation from molecular dynamics simulations. J. Phys. Chem. B. (101:7361-7363). This paper presents a detailed comparison of DNA structure and dynamics for both force fields from extended simulation times of 10 ns each. Average structures display an A-DNA base geometry with the CHARMM force field and a base geometry that is intermediate between A- and B-DNA with the AMBER force field. The backbone assumes B form on both strands with the AMBER force field, while the CHARMM force field produces heterogeneous structures with the purine strand in A form and the pyrimidine strand in dynamical equilibrium between A and B conformations. The results compare well with experimental data for the cytosine/guanine part but fail to fully reproduce an overall B conformation in the thymine/adenine tract expected from crystallographic data, particularly with the CHARMM force field. Fluctuations between A and B conformations are observed on the nanosecond time scale in both simulations, particularly with the AMBER force field. Different dynamical behavior during the first 4 ns indicates that convergence times of several nanoseconds are necessary to fully establish a dynamical equilibrium in all structural quantities on the time scale of the simulations presented here.