The conformation of a B-DNA decamer is mainly determined by its sequence and not by crystal environment

Dissection of integrated readout reveals the structural thermodynamics of DNA selection by transcription factors

Nucleobases such as inosine have been extensively utilized to map direct contacts by proteins in the DNA groove. Their deployment as targeted probes of dynamics and hydration, which are dominant thermodynamic drivers of affinity and specificity, has been limited by a paucity of suitable experimental models. We report a joint crystallographic, thermodynamic, and computational study of the bidentate complex of the arginine side chain with a Watson-Crick guanine (Arg×GC), a highly specific configuration adopted by major transcription factors throughout the eukaryotic branches in the Tree of Life. Using the ETS-family factor PU.1 as a high-resolution structural framework, inosine substitution for guanine resulted in a sharp dissection of conformational dynamics and hydration and elucidated their role in the DNA specificity of PU.1. Our work suggests an under-exploited utility of modified nucleobases in untangling the structural thermodynamics of interactions, such as the Arg×GC motif, where direct and indirect readout are tightly integrated.

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.

The structure of the Klf4 DNA-binding domain links to self-renewal and macrophage differentiation

Krueppel-like factor 4 (Klf4) belongs to the Sp/Klf family of zinc-finger transcription factors and is indispensable for terminal maturation of epithelial tissues. Furthermore, it is part of a small set of proteins that are used to generate pluripotent embryonic stem cells from differentiated tissues. Herein, we describe that a Klf4 zinc-finger domain mutant induces self-renewal and block of maturation, while wild-type Klf4 induces terminal macrophage differentiation. Moreover, we present the crystal structure of the zinc-finger domain of Klf4 bound to its target DNA, revealing that primarily the two C-terminal zinc-finger motifs are required for site specificity. Lack of those two zinc fingers leads to deficiency of Klf4 to induce macrophage differentiation. The first zinc finger, on the other hand, inhibits the otherwise cryptic self-renewal and block of differentiation activity of Klf4. Our data show that impairing the DNA binding could potentially contribute to a monocytic leukemia.

Lipooligosaccharide-deficient Neisseria meningitidis shows altered pilus-associated characteristics

Molecular interaction between host mucosal surfaces and outer membrane components of microbes is crucial in the infection process. The outer membrane of pathogenic Neisseria contains surface molecules such as pili, PilC, and Opa and a monolayer of lipooligosaccharide (LOS), all of which are involved in the interaction with host cells. Pili mediate the initial attachment to human epithelial cells, which is followed by tight contact between bacteria and the eucaryotic cells, leading to bacterial invasion. To further examine the basis for bacterium-host cell contact, we constructed an LOS-deficient Neisseria meningitidis serogroup C mutant. LOS deficiency was without exception accompanied by altered colony opacity and morphology, which most likely represented an "on" switch for Opa540 expression, and by reduced levels of the iron-regulated proteins FetA and FbpA. We show here that LOS is essential for pilus-associated adherence but dispensable for fiber formation and twitching motility. The absence of attachment to epithelial cells could not be attributed to altered levels of piliation or defects in the pilus adhesion phenotype. Further, LOS mutants do not invade host cells and have lost the natural competence for genetic transformation.

Atomic-resolution crystal structures of B-DNA reveal specific influences of divalent metal ions on conformation and packing

Crystal structures of B-form DNA have provided insights into the global and local conformational properties of the double helix, the solvent environment, drug binding and DNA packing. For example, structures of the duplex with sequence CGCGAATTCGCG, the Dickerson-Drew dodecamer (DDD), established a unique geometry of the central A-tract and a hydration spine in the minor groove. However, our knowledge of the various interaction modes between metal ions and DNA is very limited and almost no information exists concerning the origins of the different effects on DNA conformation and packing exerted by individual metal ions. Crystallization of the DDD duplex in the presence of Mg(2+)and Ca(2+)yields different crystal forms. The structures of the new Ca(2+)-form and isomorphous structures of oligonucleotides with sequences GGCGAATTCGCG and GCGAATTCGCG were determined at a maximum resolution of 1.3 A. These and the 1.1 A structure of the DDD Mg(2+)-form have revealed the most detailed picture yet of the ionic environment of B-DNA. In the Mg(2+)and Ca(2+)-forms, duplexes in the crystal lattice are surrounded by 13 magnesium and 11 calcium ions, respectively.Mg(2+)and Ca(2+)generate different DNA crystal lattices and stabilize different end-to-end overlaps and lateral contacts between duplexes, thus using different strategies for reducing the effective repeat length of the helix to ten base-pairs. Mg(2+)crystals allow the two outermost base-pairs at either end to interact laterally via minor groove H-bonds, turning the 12-mer into an effective 10-mer. Ca(2+)crystals, in contrast, unpair the outermost base-pair at each end, converting the helix into a 10-mer that can stack along its axis. This reduction of a 12-mer into a functional 10-mer is followed no matter what the detailed nature of the 5'-end of the chain: C-G-C-G-A-ellipsis, G-G-C-G-A-ellipsis, or a truncated G-C-G-A-ellipsis Rather than merely mediating close contacts between phosphate groups, ions are at the origin of many well-known features of the DDD duplex structure. A Mg(2+)coordinates in the major groove, contributing to kinking of the duplex at one end. While Ca(2+)resides in the minor groove, coordinating to bases via its hydration shell, two magnesium ions are located at the periphery of the minor groove, bridging phosphate groups from opposite strands and contracting the groove at one border of the A-tract.

The B-DNA dodecamer at high resolution reveals a spine of water on sodium

We describe a very accurate addition (called structure X here) to the B-DNA dodecamer family of X-ray structures. Our results confirm the observation of Drew and Dickerson [(1981) J. Mol. Biol. 151, 535-556] that the spine of hydration in AT tract DNA is two layers deep. However, our results suggest that the primary spine is partially occupied by sodium ions. We suggest that many sequence-dependent features of DNA conformation are mediated by site specific binding of cations. For example, preferential localization of cations, as described here within the minor groove of structure X, is probably the structural origin of AT tract bending and groove narrowing. The secondary spine, which does not interact directly with the DNA, is as geometrically regular as the primary spine, providing a model for transmission of sequence information into solvent regions. A fully hydrated magnesium ion located in the major groove of structure X appears to pull cytosine bases partially out from the helical stack, exposing pi-systems to partial positive charges of the magnesium ion and its outer sphere. A partially ordered spermine molecule is located within the major groove of structure X. Dodecamer structures are derived from crystals of [d(CGCGAATTCGCG)]2 in space group P212121 (a = 25 A, b = 40 A, and c = 66 A). On average, those crystals diffracted to around 2.5 A resolution with 2500 unique reflections. Structure X, with the same space group, DNA sequence, and crystal form as the "Dickerson dodecamer", is refined against a complete, low-temperature, 1.4 A resolution data set, with over 11000 reflections. Structure X appears to be conformationally more ordered than previous structures, suggesting that at least a portion of the conformational heterogeneity previously attributed to DNA sequence in fact arises from experimental error.

Sequence-dependent crossed helix packing in the crystal structure of a B-DNA decamer yields a detailed model for the Holliday junction

The structure of the B-DNA decamer d(CGCAATTGCG)2 has been determined by X-ray diffraction analysis to a resolution of 2.3 A and an R-factor of 17.7%. The decamer crystallises in the monoclinic space group C2 and packs with a crossed arrangement of helices and a unique crossing contact distinct from all other decamer structures. This is believed to be a direct result of the sequence-dependent minor groove width of the duplex. Crossed helix structures of DNA are valuable starting points for modelling studies of the Holliday junction. Two unique sites are observed at the cross-over junction where strand exchange may occur. A Holliday junction model has been constructed for each case and modelled using molecular mechanics and dynamics techniques. One of these models was found to be fully consistent with the available physical data.

Molecular dynamics simulations of the bis-intercalated complexes of ditercalinium and Flexi-Di with the hexanucleotide d(GCGCGC)2: theoretical analysis of the interaction and rationale for the sequence binding specificity

The X-ray crystal structures of the complexes of ditercalinium and Flexi-Di with d(CGCG)2 have been studied by computational chemistry methods in an attempt to rationalize their distinct structural features. In addition, the complexes of these two bisintercalating drugs with d(GCGCGC)2 have been modeled and subjected to 0.5 ns of molecular dynamics simulations in explicit solvent with the aim of evaluating the relative importance of hydrogen bonding and stacking interactions in the sequence binding specificity of these compounds. According to our calculations, the electrostatic term is attractive for the stacking interactions between the pyridocarbazole chromophores of these drugs and the base pairs that make up the sandwiched GpC step. On the contrary, this energy term is repulsive for the base pairs that make up the boundaries of the bisintercalation site. This differential electrostatic binding energy component, which is shown to have a strong orientational dependence, could lie at the origin of the observed binding preferences of these drugs. In addition, both the Lennard-Jones and the electrostatic energy terms contribute to stabilizing the underwound central GpC step. The attractive electrostatic interactions between the linkers and the major groove are in concert with the stacking specificities for the sandwiched GpC step, which is thus very effectively stapled by the drugs. The hydrogen-bonding potential of the linkers, however, appears to be reduced in an aqueous medium due to competing interactions with water. Binding of either ditercalinium or Flexi-Di to d(GCGCGC)2 appears to favor the A-type conformation that this DNA molecule most likely adopts in the free state. The possible relevance of these findings to the process of bis-intercalation and to the pharmacological action of these compounds is discussed.

Sequence-dependent structural variation in B-DNA

Though fiber diffraction originally led to the belief that the structure of DNA would be a simple regular helix, X-ray crystallography of synthetic oligomers has shown that both deformability and structure depend on sequence. But the rules that determine these factors remain mysterious.

The structure of B-DNA in oriented fibers

Native, general sequence B-form DNA in uniaxially oriented fibers is a ten-fold helix with identical antiparallel strands: this is to say the molecular symmetry is 2 2 10(1). The diffraction patterns indicate that local variations, however significant, must be modest. This is true also for the lithium salt of calf thymus DNA in fibers that are polycrystalline as well as oriented. The contents of its orthorhombic unit cells are arranged with P2(1)2(1)2(1) symmetry which permits the molecular symmetry to be merely two-fold. The molecular structure of DNA in such conditions resembles, conformationally and molecularly, that of B-type DNA in oligonucleotide single crystals and in oriented polycrystalline fibers of polyoligonucleotides, and therefore provides a basis for evaluating the variations that may be due to sequence effects in polyoligonucleotides in fibers and oligonucleotides in single crystals.

Refinement of netropsin bound to DNA: bias and feedback in electron density map interpretation

The X-ray crystal structure of the complex of the B-DNA dodecamer CGCGAATTCGCG with the antitumor drug netropsin has been reexamined to locate the drug accurately for computer-based drug design. The optimum solution is with the drug centered in the AATT region of the minor groove, making three good bifurcated hydrogen bonds with adenine N3 and thymine O2 atoms along the floor of the groove. Pyrrole rings of netropsin are packed against the C2 positions of adenines, leaving no room for the amine group of guanine and, hence, providing a structural rationale for the A.T specificity of netropsin. An alternative positioning in which the drug is shifted along the minor groove by ca. one-half base pair step is rejected on the basis of free R factor calculations and the appearance of the original drug-free difference maps. Final omit maps, although of more pleasing appearance, are not a dependable means of discriminating between right and wrong structures. The shifted alternative drug position ignores potential hydrogen bonding along the floor of the groove, provides no explanation for netropsin's observed A.T specificity, and is contradicted by NMR results [Patel, D. J. (1982) Proc. Natl. Acad. Sci. U.S.A. 79, 6424].

Sequence dependent effects of CpG cytosine methylation. A joint 1H-NMR and 31P-NMR study

The impact of cytosine methylation in the central CpG step of two closely related octanucleotide duplexes d(CATCGATG)2 and d(CTTCGAAG)2 was examined by 1H-NMR and 31P-NMR experiments, and a quantitative structural analysis was performed using the NOE-derived distances, the sugar puckers and the epsilon torsion angles. The two starting oligonucleotides displayed a B-DNA conformation with, however, significant local structure differences. Although the methylated oligonucleotides retained their B-DNA conformation, different structural and thermal stability effects were observed. The magnitude of the methylation effects was to depend on the initial conformation of the CpG site, which is governed by the nature of the dinucleotide AT or TT located on the CpG flanks. As an example of sequence dependence, the methylation of CpG entailed larger conformational variation in d(CATCGATG)2 than in d(CTTCGAAG)2. In this study, the 1H and 31P chemical-shift parameters averred as extremely sensitive probes for detecting subtle conformational changes. Finally, our comparative results may aid our understanding of the structural and related biological effects produced by cytosine methylation in DNA.

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.

Complexes between osmium tetraoxide bispyridine and DNA: a molecular mechanics study

This study analyzes changes in conformation of five double-stranded DNA fragments related to binding of the single-strand selective probe osmium tetraoxide bispyridine (Os, py) to a thymine. Molecular mechanics was used to investigate four B-DNA and one A-DNA fragments including two structures containing a G.T mispair. The reactivity of a particular thymine was estimated by the difference in energies between interactions in a refined DNA fragment and the corresponding interactions in the transformed fragment with Os, py. Both calculations with and without counterions were performed and the results were in qualitative agreement with experiments. The energetically relaxed fragments with Os, py showed relatively minor global structural changes in comparison to the relaxed fragments without Os, py probe. The computed structures of fragments with Watson-Crick pairing that enable binding of Os, py had similar structural characteristics to geometries found in X-ray studies of single-base mismatches. The possible role of ions in binding Os, py to mispairs is also discussed.

Sequence changes in the pilus subunit lead to tropism variation of Neisseria gonorrhoeae to human tissue

Pili of Neisseria gonorrhoeae are correlated with increased bacterial attachment to epithelial cells and undergo both phase and antigenic variation. Phase variation of gonococcal pili can be brought about by recombination events in the pilin structural gene, pilE, or by the on/off switch in expression of PilC, a pilus biogenesis protein for which two loci exist. We have studied the binding to epithelial cell lines and to fixed tissue sections of N. gonorrhoeae MS11 derivatives and mutants carrying structurally defined PilE and PilC proteins. In situ binding studies of N. gonorrhoeae to formalin-fixed tissue sections resulted in a binding pattern similar to that obtained using viable epithelial cell lines of different origin. Piliated gonococcal clones, containing different pilE sequences, varied dramatically from one another in their efficiencies at binding to corneal and conjunctival tissue, but bound equally well to cervical and endometrial tissues. Further, the binding data suggested that PilC expression by itself, i.e. without pili, cannot confer bacterial binding and that expression of either PilC1 or PilC2 does not confer different binding properties to the bacterial cells. Possible receptors for piliated gonococci were expressed in human tissues, such as cervix, endometrium, cornea, intestine, stomach, mid-brain and meninges, but not in human kidney. Pretreatment of the target tissues with Proteinase K decreased the gonococcal binding dramatically, whereas pretreatment with neuraminidase and meta-periodate, which cleave carbon-carbon linkages between vicinal hydroxyl groups in carbohydrates, did not affect attachment of gonococci. These data argue that pilus-dependent attachment of N. gonorrhoeae to human tissue may be mediated by a eukaryotic receptor having protein characteristics, and that the pilus subunit sequence may play an important role in the interaction with human cornea.

Crystallographic studies of DNA helix structure

X-ray crystallography can reveal the three-dimensional structure of short fragments of DNA or RNA with unique precision. It provides information concerning both the global helical structure and the geometry of local features such as base-pair stacking patterns and backbone conformation. An analysis of the structures of a family of DNA decamers with related sequences, crystallizing in a number of different lattices, defines the ranges in which conformational parameters can vary in B-DNA helices and shows the correlations between them. Thus, these studies show the static structures and give insight into the mechanics of DNA helices by showing how a change of one local conformational parameter will influence others. Crystal structures are also used to assess the competing influences of nucleotide sequence and environment on the three-dimensional DNA structure. To extrapolate from DNA crystal structures to physical characteristics and function of these molecules in solution or embedded into a defined sequence context remains a major challenge.

Theoretical analysis of the base stacking in DNA: choice of the force field and a comparison with the oligonucleotide crystal structures

It follows from previous studies that changes in the base pair vertical separation (BPVS) influence the architecture of DNA much more than any other conformational parameter. This inspired us to compare BPVS in the available oligonucleotide crystal structures with the optimum values provided by nine different empirical potentials employed in the theoretical studies of DNA conformation. This comparison shows that BPVS is reproduced by three fields in all steps of the highly resolved oligonucleotide crystal structures while the remaining six empirical potentials, including AMBER, GROMOS and CHARMM, provide systematic deviations. We further find that the base pairs are poorly stacked (mostly compressed) in some other refined DNA crystal structures. Our analysis indicates that this poor stacking originates from improperly determined positions of the bases. The approach described in the present communication can be used to identify DNA structures which are not accurate enough for studies of the relationships between the base sequence and DNA conformation.

Molecular dynamics simulations of B-DNA: an analysis of the role of initial molecular configuration, randomly assigned velocity distribution, long integration times, and nonconstrained termini

Molecular dynamics simulations of three DNA sequences using the AMBER 3.0 force field were performed with implicit inclusion of water through a distance-dependent dielectric constant and solvated counterions. Simulations of the self-complementary DNA dodecamer d(CGCGAATTCGCG) were started from a regular B-DNA structure and the x-ray single crystal B-DNA structure. Although mean convergence during the 89-ps calculation was confirmed, localized differences in backbone torsionals and base-pair helicoidals were observed. A nanosecond simulation of the nonself-complementary 14 base-pair DNA d(GGCGGAATTGGCGG) indicates that most structural parameters stabilize within the first 100-200 ps, while isolated features show low-frequency oscillations throughout the calculation. The lack of harmonic constraints on the ends of the molecules was shown not to perturb the structural dynamics of the internal oligonucleotide beyond the external 2 base pairs. Comparison of three simulations of the nonself-complementary 14 base-pair DNA d(GGCGAAATTCGCGG), identical in all respects other than the assignment of initial Maxwellian atomic velocity distributions, revealed the inherent systematic variability. The three calculations result in nearly superimposable global structures, with localized variations in torsionals and helicoidals. Our results provide a basis for performing a comparative analysis of the effect of DNA sequence on localized structure.

Structure of the B-DNA decamer C-C-A-A-C-I-T-T-G-G in two different space groups: conformational flexibility of B-DNA

For the first time, the same B-DNA oligomer has been crystallized and its structure solved in two different space groups. Crystallization of C-C-A-A-C-I-T-T-G-G with Ca2+ yields monoclinic space group C2 with a = 31.87 A, b = 25.69 A, c = 34.21 A, beta = 114.1 degrees, and five base pairs per asymmetric unit. The 5026 2 sigma data to 1.3 A refine to R = 0.152 with 72 waters, one heptavalent hydrated calcium complex, and one cacodylate ion per asymmetric unit. In contrast, crystallization with Mg2+ yields trigonal space group P3(2)21 with a = b = 33.23 A, c = 94.77 A, gamma = 120 degrees, and 10 base pairs per asymmetric unit. The 1725 2 sigma data to 2.2 A refine to R = 0.164 with 36 water molecules and one octahedral magnesium complex per asymmetric unit. The monoclinic form is virtually isostructural with previously solved monoclinic decamers, including twist angles of ca. 50 degrees at C-A and T-G steps. In contrast, the trigonal structure has quite different local helix parameters, with twist angles of ca. 36 degrees at the corresponding steps. These local parameter differences can only be attributed to crystal packing, suggesting that certain sequences of B-DNA are more flexible and influenced by their surroundings than had previously been thought. Such deformability may be important for interaction of B-DNA with control proteins, where both static structure and dynamic deformability comprise components of the recognition process. The crossing of two helices at an angle of 120 degrees in the trigonal cell is a model for an antiparallel, uncrossed Holliday junction, as has been noted earlier by Timsit and Moras [Timsit, Y., & Moras, D. (1991) J. Mol. Biol. 221, 919-940] from a rhombohedral DNA dodecamer structure analysis.

An accurate method for determining the helical repeat of DNA in solution reveals differences to the crystal structures of two B-DNA decamers

Many DNA sequences have been studied by X-ray crystallography with the goal of deciphering a sequence-structure code. We have determined the helical repeats of two B-type DNA decamers in solution employing an electrophoretic method based on phasing of bent segments. The decamers contain recognition sites for the dcm methyltransferase and for the restriction nuclease NarI with a mutational hotspot. Their helical repeats are 10.59(+/- 0.05) bp and 10.52(+/- 0.03) bp, respectively, whereas crystallographic analysis yielded 10.0 bp in the solid state. This difference is greater than that for the transition between B- and A-type DNA in solution. Thus, reliable information about the polymorphism of DNA in solution must be based on both X-ray and solution data. We describe a generally applicable approach to accurately determine helical repeats of small DNA duplexes in solution.

X-ray structure of the DNase I-d(GGTATACC)2 complex at 2.3 A resolution

The crystal structure of a complex between DNase I and the self-complementary octamer duplex d(GGTATACC)2 has been solved using the molecular replacement method and refined to a crystallographic R-factor of 18.8% for all data between 6.0 and 2.3 A resolution. In contrast to the structure of the DNase I-d(GCGATCGC)2 complex solved previously, the DNA remains uncleaved in the crystal. The general architecture of the two complexes is highly similar. DNase I binds in the minor groove of a right-handed DNA duplex, and to the phosphate backbones on either side over five base-pairs, resulting in a widening of the minor groove and a concurrent bend of the DNA away from the bound enzyme. There is very little change in the structure of the DNase I on binding the substrate. Many other features of the interaction are conserved in the two complexes, in particular the stacking of a deoxyribose group of the DNA onto the side-chain of a tyrosine residue (Y76), which affects the DNA conformation and the binding of an arginine side-chain in the minor groove. Although the structures of the DNA molecules appear at first sight rather similar, detailed analysis reveals some differences that may explain the relative resistance of the d(GGTATACC)2 duplex to cleavage by DNase I: whilst some backbone parameters are characteristic of a B-conformation, the spatial orientation of the base-pairs in the d(GGTATACC)2 duplex is close to that generally observed in A-DNA. These results further support the hypothesis that the minor-groove width and depth and the intrinsic flexibility of DNA are the most important parameters affecting the interaction. The disposition of residues around the scissile phosphate group suggests that two histidine residues, H134 and H252, are involved in catalysis.

Groove width and depth of B-DNA structures depend on local variation in slide

The groove widths of DNA helix, especially minor groove width, are generally believed to be important for recognition of DNA by various types of ligands. It has been postulated earlier that large negative propeller twist, in the AT rich regions compresses the minor groove of duplex DNA. A systematic study has now been carried out by generating models with different values of local doublet and intra-basepair parameters and calculating their minor groove widths. It is found that several local doublet parameters affect the minor groove width but it depends most strongly on the local step parameters roll and slide when each parameter is considered individually. However, a detailed analysis of the various local parameters within the B-DNA family of crystal structures indicates that propeller twist and slide are most strongly correlated with the observed values of minor groove width. The groove depth is also strongly correlated with slide. Thus the local base sequence dependent variations in slide can modify both the groove width and depth and consequently determine the ligand binding properties of DNA.

Structure of poly d(A).poly d(T)

On the basis of the x-ray data from polycrystalline and well oriented fibers of the sodium salt of poly d(A).poly d(T) (Arnott et al, Nucl. Acids Res. 11, 4141-4155 (1983), a revised B'-DNA model incorporating B-like adenine and thymine strands is shown to give a much better x-ray agreement (R = 0.25) than the previously assigned model consisting of mixed sugar conformations in the two strands. The narrowing of the minor and the widening of the major grooves are promiscuous features of B'-DNA, which are common to all poly d(purine).poly d(pyrimidine) duplexes with two hydrogen bonded base-pairs and are in marked contrast with classical B-DNA. Due to modest propeller (-15 degrees), the cross strand diagonal hydrogen bonds (0.37 nm) in this duplex are not as strong as those in A,T-rich oligonucleotide crystal structures.

Structure of a B-DNA decamer with a central T-A step: C-G-A-T-T-A-A-T-C-G

The X-ray crystal structure analysis of the decamer C-G-A-T-T-A-A-T-C-G has been carried out to a resolution of 1.5 A. The crystals are space group P2(1)2(1)2(1), cell dimensions a = 38.60 A, b = 39.10 A, c = 33.07 A. The structure was solved by molecular replacement and refined with X-PLOR and NUCLSQ. The final R factor for a model with 404 DNA atoms, 108 water molecules and one magnesium hexahydrate cation is 15.7%. The double helix is essentially isostructural with C-G-A-T-C-G-A-T-C-G, with closely similar local helix parameters. The structure of the T-T-A-A center differs from that found in C-G-C-G-T-T-A-A-C-G-C-G in that the minor groove in our decamer is wide at the central T-A step rather than narrow, and the twist angle of the T-A step is small (31.1 degrees) rather than large. Whereas the tetrad model provides a convenient framework for discussing local DNA helix structure, it cannot be the entire story. The articulated helix model of DNA structure proposes that certain sequence regions of DNA show preferential twisting or bending properties, whereas other regions are less capable of deformation, in a manner that may be useful in sequence recognition by drugs and protein. Further crystal structure analyses should help to delineate the precise nature of sequence-dependent articulation in the DNA double helix.

Three-dimensional structures of bulge-containing DNA fragments

The three-dimensional structure of a DNA tridecamer d(CGCAGAATTCGCG)2 containing bulged adenine bases was determined by single crystal X-ray diffraction methods, at 120 K, to 2.6 A resolution. The structure is a B-DNA type double helix with a single duplex in the asymmetric unit. One of the bulged adenine bases loops out from the double helix, while the other stacks in to it. This is in contrast to our preliminary finding, which indicated that both adenine bases were looped out. This revised model was confirmed by the use of a covalently bound heavy-atom derivative. The conformation of the looped-out bulge hardly disrupts base stacking interactions of the bases flanking it. This is achieved by the backbone making a "loop-the-loop" curve with the extra adenine flipping over with respect to the other nucleotides in the strand. The looped-out base intercalates into the stacked-in bulge site of a symmetrically related duplex. The looped-out and stacked-in bases form an A.A reversed Hoogsteen base-pair that stacks between the surrounding base-pairs, thus stabilizing both bulges. The double helix is frayed at one end with the two "melted" bases participating in intermolecular interactions. A related structure, of the same tridecamer, after soaking the crystals with proflavin, was determined to 3.2 A resolution. The main features of this B-DNA duplex are basically similar to the native tridecamer but differ in detail especially in the conformation of the bulged-out base. Accommodation of a large perturbation such as that described here with minimal disruption of the double helix shows both the flexibility and resiliency of the DNA molecule.

Modeling DNA structures

For a molecule of biological importance, one expects a strong correlation between the three-dimensional structure and its biological function(s). Molecular simulations allow the prediction of physical properties of macromolecules. Many of these properties are closely related to the molecular structure. Model-building studies may thus supplement structurally low-resolution experimental data with detailed three-dimensional hypothetical atomic models. Such studies may give a consistent integral view of a wealth of experimental data. In most cases, such models will predict the outcome of certain experiments. Their actual results will often either confirm the model, be used for further refinement, or demand a major revision. Empirical force-fields provide a large amount of physicochemical knowledge concerning structural and other physical properties about various classes of molecules. They give reasonable bond distances and angles and prevent short van der Waals contacts. Difficulties arise for the prediction of large-scale structural elements. These are not only determined by short-range interactions, but also result from long-range electrostatic, hydration, and hydrophobic forces. In the case of the ionically driven DNA B-Z transition, the point of transition as a function of ionic strength and size can be correctly predicted (33). More must be done for a better understanding of the hydration forces (6). What type of questions will be reliably answered by a force-field? Relatively safe answers concern the local geometry of the molecules. If a conformation leads to strong distortions of bond distances or angles or to close van der Waals contacts, it can safely be rejected. Optimizing such unfavorable structures energetically may lead to structures showing how to avoid such distortions. More difficult are energetic questions: Which of two conformers is more stable, or what is the free energy of the substrate in the active site (63)? One cannot always be sure that the force-field provides the correct answers. Therefore, one should concentrate on questions that can be checked experimentally. The application of such concepts to model curved DNA (19, 60) and the DNA four-way junction (61) provides promising results. To explore the knowledge contained in the force-fields, several methods have been proposed. Taking advantage of structural symmetries may improve critically the convergence, while refining the target molecule or its building blocks. A numerically stable derivative for the torsion potential has been proposed. The optimization method of conjugated gradients (see Section II,B) is a powerful tool to find the way downhill toward a local minimum. To surmount barriers and escape local minima requires nonlocal optimization procedures.(ABSTRACT TRUNCATED AT 400 WORDS)

Modeling DNA structures: molecular mechanics and molecular dynamics

Model building studies may be used to supplement structurally low resolution experimental data with detailed three-dimensional hypothetical atomic models. Because of the strong relation between structure and function in biological molecules such models may give a consistent, integral view of a wealth of experimental data. In most cases such models will predict the outcome of certain experiments. The outcome of these experiments will often either confirm the model may be used for further refinement or even demand a major revision of the model. Coordinates obtained from X-ray fiber diffraction data or in special cases single-crystal data may provide the elements for DNA or RNA model building. Local and nonlocal optimization may be used to refine these structures and to evaluate their statistical significance as estimated by a chosen force field. Appreciable progress using nonlocal optimization procedures can only be expected if the dimensionality of the problem can be reduced sufficiently to the relevant degrees of freedom. Taking advantage of structural symmetries may critically improve the convergence while refining the target molecule or its building blocks. Monte Carlo and molecular dynamics methods allow one to calculate averaged quantities. In addition, molecular dynamics provides time evolutions of certain averages. During the simulation of certain physical properties of molecules a huge amount of data will be generated. They will provide many answers, but these answers may not always apply to the original question. So what type of questions will be reliably answered by a force field? Relatively safe answers concern the local geometry of the molecules. If a conformation leads to strong distortions of bond distances or angles or to close van der Waals contacts, this conformation can safely be rejected. Optimizing such unfavorable structures energetically may lead to structures showing how to avoid such distortions. More difficult are energetical questions: which of two conformers is more stable, or what is the free energy of the substrate in the active site? One cannot always be sure that the force field provides the correct answer. Therefore, one should pose only those questions which can be checked experimentally. Because of the many possible answers, the experiment may benefit by starting with a choice proposed by the simulation. The application of this procedure to curved DNA and the DNA four-way junction was successful.

Local supercoil-stabilized DNA structures

The DNA double helix exhibits local sequence-dependent polymorphism at the level of the single base pair and dinucleotide step. Curvature of the DNA molecule occurs in DNA regions with a specific type of nucleotide sequence periodicities. Negative supercoiling induces in vitro local nucleotide sequence-dependent DNA structures such as cruciforms, left-handed DNA, multistranded structures, etc. Techniques based on chemical probes have been proposed that make it possible to study DNA local structures in cells. Recent results suggest that the local DNA structures observed in vitro exist in the cell, but their occurrence and structural details are dependent on the DNA superhelical density in the cell and can be related to some cellular processes.