Different families of double-stranded conformations of DNA as revealed by computer calculations

Pseudoknot structures with conserved base triples in telomerase RNAs of ciliates

Telomerase maintains the integrity of telomeres, the ends of linear chromosomes, by adding G-rich repeats to their 3′-ends. Telomerase RNA is an integral component of telomerase. It contains a template for the synthesis of the telomeric repeats by the telomerase reverse transcriptase. Although telomerase RNAs of different organisms are very diverse in their sequences, a functional non-template element, a pseudoknot, was predicted in all of them. Pseudoknot elements in human and the budding yeast Kluyveromyces lactis telomerase RNAs contain unusual triple-helical segments with AUU base triples, which are critical for telomerase function. Such base triples in ciliates have not been previously reported. We analyzed the pseudoknot sequences in 28 ciliate species and classified them in six different groups based on the lengths of the stems and loops composing the pseudoknot. Using miniCarlo, a helical parameter-based modeling program, we calculated 3D models for a representative of each morphological group. In all cases, the predicted structure contains at least one AUU base triple in stem 2, except for that of Colpidium colpoda, which contains unconventional GCG and AUA triples. These results suggest that base triples in a pseudoknot element are a conserved feature of all telomerases.

A triple helix within a pseudoknot is a conserved and essential element of telomerase RNA

Telomerase copies a short template within its integral telomerase RNA onto eukaryotic chromosome ends, compensating for incomplete replication and degradation. Telomerase action extends the proliferative potential of cells, and thus it is implicated in cancer and aging. Nontemplate regions of telomerase RNA are also crucial for telomerase function. However, they are highly divergent in sequence among species, and their roles are largely unclear. Using in silico three-dimensional modeling, constrained by mutational analysis, we propose a three-dimensional model for a pseudoknot in telomerase RNA of the budding yeast Kluyveromyces lactis. Interestingly, this structure includes a U-A.U major-groove triple helix. We confirmed the triple-helix formation in vitro using oligoribonucleotides and showed that it is essential for telomerase function in vivo. While triplex-disrupting mutations abolished telomerase function, triple compensatory mutations that formed pH-dependent G-C.C(+) triples restored the pseudoknot structure in a pH-dependent manner and partly restored telomerase function in vivo. In addition, we identified a novel type of triple helix that is formed by G-C.U triples, which also partly restored the pseudoknot structure and function. We propose that this unusual structure, so far found only in telomerase RNA, provides an essential and conserved telomerase-specific function.

NMR structure of the full-length linear dimer of stem-loop-1 RNA in the HIV-1 dimer initiation site

The packaging signal of HIV-1 RNA contains a stem-loop structure, SL1, which serves as the dimerization initiation site for two identical copies of the genome and is important for packaging of the RNA genome into the budding virion and for overall infectivity. SL1 spontaneously dimerizes via a palindromic hexanucleotide sequence in its apical loop, forming a metastable kissing dimer form. Incubation with nucleocapsid protein causes this form to refold to a thermodynamically stable mature linear dimer. Here, we present an NMR structure of the latter form of the full-length SL1 sequence of the Lai HIV-1 isolate. The structure was refined using nuclear Overhauser effect and residual dipolar coupling data. The structure presents a symmetric homodimer of two RNA strands of 35 nucleotides each; it includes five stems separated by four internal loops. The central palindromic stem is surrounded by two symmetric adenine-rich 1-2 internal loops, A-bulges. All three adenines in each A-bulge are stacked inside the helix, consistent with the solution structures of shorter SL1 constructs determined previously. The outer 4-base pair stems and, proximal to them, purine-rich 1-3 internal loops, or G-bulges, are the least stable parts of the molecule. The G-bulges display high conformational variability in the refined ensemble of structures, despite the availability of many structural restraints for this region. Nevertheless, most conformations share a similar structural motif: a guanine and an adenine from opposite strands form a GA mismatch stacked on the top of the neighboring stem. The two remaining guanines are exposed, one in the minor groove and another in the major groove side of the helix, consistent with secondary structure probing data for SL1. These guanines may be recognized by the nucleocapsid protein, which binds tightly to the G-bulge in vitro.

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

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

Pulling chromatin fibers: computer simulations of direct physical micromanipulations

A low-resolution molecular model, which combines the known mechanical properties of protein-free DNA with the accumulating picture of chromatosome structure, has been developed to account for the stretching of single chromatin fibers by an imposed external force. Force-extension characteristics of sets of chains accumulated by Monte Carlo sampling are consistent with recently observed findings in the non-destructive regime (<20 pN imposed force), where the structure of the chromatosome remains intact. The correspondence between simulation and the relaxation phase of the experiment limits the equilibrium entry-exit angle of linker DNA on the chromatosome to W=50(+/-10) degrees and the effective DNA linker length to L(eff)=40(+/-5) bp. The computed force-extension characteristics are relatively insensitive to other parameters of the model, precluding their accurate estimation. The introduction of an attractive potential between closely spaced nucleosomes reproduces the added initial resistance of single fibers to extension at high salt conditions. The consideration of elastic linkers also improves the fitting of assorted classical measurements of unstressed chromatin structure in solution. The overall picture of chromatin that emerges is an irregular, fluctuating, three-dimensional, zig-zag structure with intact, mechanically stable chromatosome units and deformable linkers. The modeled fiber undergoes large-scale configurational rearrangements without significant perturbation of the constituent chromatosome beads, collapsing into a highly condensed form in response to small (<2kT) inter-nucleosomal attractions.

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.

Electrostatic interaction between helical macromolecules in dense aggregates: an impetus for DNA poly- and meso-morphism

DNA exhibits a surprising multiplicity of structures when it is packed into dense aggregates. It undergoes various polymorphous transitions (e.g., from the B to A form) and mesomorphous transformations (from hexagonal to orthorhombic or monoclinic packing, changes in the mutual alignment of nearest neighbors, etc). In this report we show that such phenomena may have their origin in the specific helical symmetry of the charge distribution on DNA surface. Electrostatic interaction between neighboring DNA molecules exhibits strong dependence on the patterns of molecular surface groups and adsorbed counter-ions. As a result, it is affected by such structural parameters as the helical pitch, groove width, the number of base pairs per helical turn, etc. We derive expressions which relate the energy of electrostatic interaction with these parameters and with the packing variables characterizing the axial and azimuthal alignment between neighboring macromolecules. We show, in particular, that the structural changes upon the B-to-A transition reduce the electrostatic energy by approximately kcal/mol per base pair, at a random adsorption of counter ions. Ion binding into the narrow groove weakens or inverts this effect, stabilizing B-DNA, as it is presumably the case in Li+-DNA assemblies. The packing symmetry and molecular alignment in DNA aggregates are shown to be affected by the patterns of ion binding.

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.

DNA bending by small, mobile multivalent cations

We propose a purely electrostatic mechanism by which small, mobile, multivalent cations can induce DNA bending. A multivalent cation binds at the entrance to the B-DNA major groove, between the two phosphate strands, electrostatically repelling sodium counterions from the neighboring phosphates. The unscreened phosphates on both strands are strongly attracted to the groove-bound cation. This leads to groove closure, accompanied by DNA bending toward the cationic ligand. We explicitly treat the dynamic character of the cation-DNA interaction using an adiabatic approximation, noting that DNA bending is much slower than the diffusion of nonspecifically bound, mobile cations. We make semiquantitative estimates of the free energy components of bending-electrostatic (with a sigmoidal distance-dependent dielectric function), elastic, and entropic cation localization-and find that the equilibrium state is bent B-DNA stabilized with a self-localized cation. This is a bending polaron, formation of which should be critically dependent on the strength of electrostatic interaction and the concentration of highly mobile cations available for self-localization. We predict that the resultant bend will be large (approximately 20-40 degrees), smooth (because it is spread over 6 bp), and infrequent. The stability of such a bend can be variable, from transient to highly stable (static) bending, observable with standard curvature-measuring techniques. We further predict that this bending mechanism will have an unusual sequence dependence: sequences with less binding specificity will be more bent, unless the specific binding site is in the major groove.

Structure and conformation of helical nucleic acids: analysis program (SCHNAaP)

We present a new versatile program, SCHNAaP, for the analysis of double-helical nucleic acid structures. The program uses mathematically rigorous and fully reversible procedures for calculating the structural parameters: the Cambridge University Engineering Department Helix computation Scheme (CEHS) is used to determine the local helical parameters and an analogous procedure is used to determine the global helical parameters. These parameters form a complete set that conforms to the "Cambridge Accord" on definitions and nomenclature of nucleic acid structure parameters. In addition to the two standard Watson-Crick base-pairs, the program handles mismatched base-pairs and chemically modified bases. An analysis of the sugar-phosphate backbone conformation is included. Standardized base-stacking diagrams of each dinucleotide step with reference to the mid-step triad are generated. Structures are classified as one of the four polymorphic families, A/B, Z, W or R, although W- and R-DNA (two types of hypothetical structure) have yet to be observed experimentally.

The study of possible A and B conformations of alternating DNA using a new program for conformational analysis of duplexes (CONAN)

A new program, CONAN has been designed for CONformational ANalysis of oligonucleotide duplexes with natural and modified bases. It allows to model both regular DNA fragments with different types of symmetry and irregular ones including bends, junctions, mismatched pairs and base lesions. Computations and minimization of the energy are performed in a space of internal structural variables chosen to build start structure easier and conveniently analyze the results obtained. These internal structural variables determine mutual base-base and base-sugar arrangement and sugar puckering. The analytical closure procedure is applied both to sugar rings and to backbone fragments between adjacent sugars. For more effective energy minimization, analytical gradient is calculated. The CONAN was applied to the search for low-energy conformations of poly(dA-dT).poly(dA-dT) and poly(dG-dC).poly(dG-dC) duplexes. Extended regions of low-energy A and B conformations are revealed and characterized. These regions contain structures with different relative values of helical twist, tau, for pur-pyr and pyr-pur steps, namely, conformations with tau (pur-pyr) > tau (pyr-pur) and with tau (pur-pyr) < tau (pyr-pur). Two types of sugar puckering were found for B-form low-energy conformations, the first type with all C2'-endo sugar residues and the second one - with C2'-endo purines and O1'-endo pyrimidines. The calculated conformations are compared with X-ray diffraction data for crystals and fibers and NMR data for solution.

Structure of d(T)n.d(A)n.d(T)n: the DNA triple helix has B-form geometry with C2'-endo sugar pucker

The polynucleotide helix d(T)n.d(A)n.d(T)n is the only deoxypolynucleotide triple helix for which a structure has been published, and it is generally assumed as the structural basis for studies of DNA triplexes. The helix has been assigned to an A-form conformation with C3'-endo sugar pucker by Arnott and Selsing [1974; cf. Arnott et al. (1976)]. We show here by infrared spectroscopy in D2O solution that the helix is instead B-form and that the sugar pucker is in the C2'-endo region. Distamycin A, which binds only to B-form and not to A-form helices, binds to the triple helix without displacement of the third strand, as demonstrated by CD spectroscopy and gel electrophoresis. Molecular modeling shows that a stereochemically satisfactory structure can be build using C2'-endo sugars and a displacement of the Watson-Crick base-pair center from the helix axis of 2.5 A. Helical constraints of rise per residue (h = 3.26 A) and residues per turn (n = 12) were taken from fiber diffraction experiments of Arnott and Selsing (1974). The conformational torsion angles are in the standard B-form range, and there are no short contacts. In contrast, we were unable to construct a stereochemically allowed model with A-form geometry and C3'-endo sugars. Arnott et al. (1976) observed that their model had short contacts (e.g., 2.3 A between the phosphate-dependent oxygen on the A strand and O2 in the Hoogsteen-paired thymine strand) which are generally known to be outside the allowed range.(ABSTRACT TRUNCATED AT 250 WORDS)

Solution structure of [d(ATGAGCGAATA)]2. Adjacent G:A mismatches stabilized by cross-strand base-stacking and BII phosphate groups

The solution structure of a rather unusual B-form duplex [d(ATGAGCGAATA)]2 has been determined using two-dimensional nuclear magnetic resonance (2D-NMR) and distance geometry methods. This sequence forms a stable ten base-pair B-form duplex with 3' overhangs and two pairs of adjacent G:A mismatches paired via a sheared hydrogen-bonding scheme. All non-exchangeable protons, including the stereo-specific H-5'S/H-5'R of the 3G and 7G residues, were assigned by 2D-NMR. The phosphorus spectrum was assigned using heteronuclear correlation with H-3' and H-4' reasonances. The complete assignments reveal several unusual nuclear Overhauser enhancements (NOEs) and unusual chemical shifts for the neighboring G:A mismatch pairs and their adjacent nucleotides. Inter-proton distances were derived from time-dependent NOEs and used to generate initial structures, which were further refined by iterative back-calculation of the two-dimensional nuclear Overhauser enhancement spectra; 22 final structures were calculated from the refined distance bounds. All these final structures exhibit fully wound helical structures with small penalty values against the refined distance bounds and small pair-wise root-mean-square deviation values (typically 0.5 A to 0.9 A). The two helical strands exchange base stacking at both of the two G:A mismatch sites, resulting in base stacking down each side rather than down each strand of the twisted duplex. Very large twist angles (77 degrees) were found at the G:A mismatch steps. All the final structures were found to have BII phosphate conformations at the adjacent G:A mismatch sites, consistent with observed downfield 31P chemical shifts and Monte-Carlo conformational search results. Our results support the hypothesis that 31P chemical shifts are related to backbone torsion angles. These BII phosphate conformations in the adjacent G:A mismatch step suggest that hydrogen bonding of the G:A pair G-NH2 to a nearby phosphate oxygen atom is unlikely. The unusual structure of the duplex may be stabilized by strong interstrand base stacking as well as intrastrand stacking, as indicated by excellent base overlap within the mismatch stacks.

Systematic study of nuclear Overhauser effects vis-à-vis local helical parameters, sugar puckers, and glycosidic torsions in B DNA: insensitivity of NOE to local transitions in B DNA oligonucleotides due to internal structural compensations

A method has been developed to solve structures of DNA oligomers in solution from the experimental NOE data. The method is a combination of two approaches: (1) full matrix NOESY simulations and (2) conformational calculations of DNA double helix based on generalized helical parameters. The process of the refinement of a solution structure does not involve NMR-derived interproton distance constraints; rather it consists of a direct fitting of a structure to the experimental NOE data, a weighted sum of energy, and R factor being under minimization. A helical parameters-based generation of DNA forms makes it possible to organize the search for the optimal structure more effectively, systematically varying starting conformations. The method has been used to calculate a structure for the self-complementary DNA hexamer GGATCC, which is consistent with the available experimental data. The structure belongs to the B family of forms, although the local structural heterogeneity is very strong. Sugar puckers vary from O4'-exo to C3'-exo; helical steps are open with different magnitudes toward the minor groove. Next, we have addressed the question of how uniquely the structure is defined by the existing NMR data. Different structural parameters have been systematically varied, and their effect on individual NOE's and the R factor has been studied. Two energetically conjugated parameters, sugar puckers and glycosidic angles, can be determined very reliably, because of the strong dependences of the intraresidue H6/H8 to H2'/H2''/H3' NOE's. In contrast, the local helical conformation of DNA and the geometry of base pairs proved to be underdetermined by the existing NOE information, because the effect of any helical parameter on interproton distances can be compensated by the concerted changes in other parameters.

Are time-averaged restraints necessary for nuclear magnetic resonance refinement? A model study for DNA

A recently suggested method for refinement of structural data obtained from two-dimensional nuclear magnetic resonance experiments using molecular dynamics (MD) is explored. In this method, the time-averaged values of the appropriate internal co-ordinates of the molecule, calculated from the MD trajectory, are driven by restraints towards the experimental target values. This contrasts with most refinement procedures currently in use, where restraints are applied based on the instantaneous values of the appropriate co-ordinates. Both refinement methods are applied to the EcoRI restriction site DNA hexamer d(GAATTC)2, using target nuclear Overhauser enhancement distances derived from a one nanosecond unrestrained MD simulation of this structure. The resulting refined structures are compared to the results of the unrestrained MD trajectory, which serves as our "experimental" data. We show that although both methods can yield an average structure with the correct gross morphology, the new method allows both a much more realistic picture of inherent flexibility, and reproduces fine conformational detail better, such as sequence dependency. We also analyze the very long MD trajectory generated here (longer than any previously reported for a DNA oligomer), and find that significantly shorter simulations, typical of those frequently performed, may not yield acceptably reliable values for certain structural parameters.

Four-stranded DNA helices: conformational analysis of regular poly(dT).poly(dA).poly(dA).poly(dT) helices with various types of base binding

The paper presents results obtained in conformational analysis of homopolymeric four-stranded poly(dT).poly(dA).poly(dA).poly(dT) DNA helices in which the pairs of strands with identical bases are parallel and have a two-fold symmetry axis. All possible models of base binding to yield a symmetric complex have been considered. The dihedral angles of sugar-phosphate backbones and helix parameters, which are consistent with the minima of conformational energy for four-stranded DNAs, have been determined using the results of optimization of conformational energy calculated at atom-atom approximation. Potential energy is shown to depend on the structure of base complexes and on the mutual orientation of unlike strands. Possible biological functions of four-stranded helices are discussed.

Block-units method for conformational calculations of large nucleic acid chains. II. The two-hierarchical approach and its application to conformational arrangement of the unusual T psi C loop of rabbit tRNA(Val)

The two-level hierarchical methodology is suggested for conformational calculations of large fragments of nucleic acids. The method of the first level is intended for performing a fast screening of the conformational phase space. The high-level method may be used to refine structurally important conformations. The method of the first level is the block-units method, which has been developed specially for these purposes (see part I). It has been shown that the block-units method allows the satisfactory calculation of the structure parameters of the optimal conformations of polynucleotides. The results of the conformational rearrangement calculations of the T psi C loop of the tRNA(Phe) after modification of its sequence are represented.

Structure calculations for single-stranded DNA complexed with the single-stranded DNA binding protein GP32 of bacteriophage T4: a remarkable DNA structure

In this study it is established by calculation which regular conformations single-stranded DNA and RNA can adopt in the complex with the single-stranded DNA binding protein GP32 of bacteriophage T4. In order to do so, information from previous experiments about base orientations and the length and diameter of the complexes is used together with knowledge about bond lengths and valence angles between chemical bonds. It turns out that there is only a limited set of similar conformations which are in agreement with experimental data. The arrangement of neighboring bases is such that there is ample space for aromatic residues of the protein to partly intercalate between the bases, which is in agreement with a previously proposed model for the binding domain of the protein [Prigodich, R. V., Shamoo, Y., Williams, K. R., Chase, J. W., Konigsberg, W. H., & Coleman, J. E. (1986) Biochemistry 25, 3666-3671]. Both C2'endo and C3'endo sugar conformations lead to calculated DNA conformations that are consistent with experimental data. The orientation of the O2' atoms of the sugars in RNA can explain why the binding affinity of GP32 for polyribonucleotides is lower than for polydeoxyribonucleotides.

Parallel double helices of DNA. Conformational analysis of regular helices with the second order symmetry axis

We have performed a conformational analysis of DNA double helices with parallel directed backbone strands connected with the second order symmetry axis being at the same time the helix axis. The calculations were made for homopolymers poly(dA).poly(dA), poly(dC).poly(dC), poly(dG) poly(dG), and poly(dT).poly(dT). All possible variants of hydrogen bonding of base pairs of the same name were studied for each polymer. The maps of backbone chain geometrical existence were constructed. Conformational and helical parameters corresponding to local minima of conformational energy of "parallel" DNA helices, calculated at atom-atom approximation, were determined. The dependence of conformational energy on the base pair and on the hydrogen bond type was analysed. Two major conformational advantageous for "parallel" DNA's do not depend much on the hydrogen-bonded base pair type were indicated. One of them coincided with the conformational region typical for "antiparallel" DNA, in particular for the B-form DNA. Conformational energy of "parallel" DNA depends on the base pair type and for the most part is similar to the conformational energy of "antiparallel" B-DNA.

The structure of the aflatoxin B1-DNA adduct at N7 of guanine. Theoretical intercalation and covalent adduct models

Two theoretical models are proposed for the conformational structure of both intercalated and covalent adduct complexes of aflatoxin B1, designated AFB1, with N7 of guanine of DNA. The covalent adduct model requires the DNA to kink a minimum of 39 degrees about the covalent site of the C8 and N7 atoms comprising the bond of the covalent complex. The preference of AFB1 for specific G bases within a sequence of GC content followed that of experimental studies with the added feature that for binding to the third G base of a tetramer sequence from the 3'-end, the AFB1 displayed enhanced binding at the 3' site of the targeted guanosine. Binding of AFB1 to the second G base of a tetramer sequence from the 3'-end leads to preference for a 5' site of the targeted guanosine. Inhibition of AFB1's interaction with the targeted DNA in the presence of intercalated ethidium bromide is explained by these proposed models.

Nucleic acid model building: the multiple backbone solutions associated with a given base morphology

A constrained model building procedure is used to generate nucleic acid structures of the familiar A-, B-, and Z-DNA duplexes. Attention is focused upon the multiple structural solutions associated with the arrangements of nucleic acid base pairs rather than the optimum sugar-phosphate structure. The glycosyl (chi) and sugar torsions (both the ring puckering and the exocyclic C5'-C4' (psi) torsion) are treated as independent variables and the resulting O3'...O5' distances are used as closure determinants. When such distances conform to the known geometry of phosphate chemical bonding, an intervening phosphorus atom with correct C-O-P valence angles can be located. Four sequential torsion angles--phi', omega', omega and phi--about the C3'-O3'-P-O5'-C5' bonds are then obtained as dependent variables. The resulting structures are categorized in terms of conformation, ranked in potential energy, and analyzed for torsional correlations. The numerical results are quite interesting with implications regarding nucleic acid models constructed to fit less than ideal experimental data. The multiple solutions to the problem are useful for comprehending the conformational complexities of the local sugar-phosphate backbone and for understanding the transitions between different helical forms. According to these studies, unique characterization of a nucleic acid duplex involves more than the determination of its base pair morphology, its sugar puckering preferences, or its groove binding features.

Conformational studies of nucleic acids: IV. The conformational energetics of oligonucleotides: d(ApApApA) and ApApApA

Utilizing a new method for modeling furanose pseudorotation (D. A. Pearlman and S.-H. Kim, J. Biomol. Struct. Dyn. 3, 85 (1985)) and the empirical multiple correlations between nucleic acid torsion angles we derived in the previous report (D. A. Pearlman and S.-H. Kim, previous paper in this issue), we have made an energetic examination of the entire conformational spaces available to two nucleic acid oligonucleotides: d(ApApApA) and ApApApA. The energies are calculated using a semi-empirical potential function. From the resulting body of data, energy contour map pairs (one for the DNA molecule, one for the RNA structure) have been created for each of the 21 possible torsion angle pairs in a nucleotide repeating unit. Of the 21 pairs, 15 have not been reported previously. The contour plots are different from those made earlier in that for each point in a particular angle-angle plot, the remaining five variable torsion angles are rotated to the values which give a minimum energy at this point. The contour maps are overall quite consistent with the experimental distribution of oligonucleotide data. A number of these maps are of particular interest: delta (C5'-C4'-C3'-O3')-chi (O4'-C1'-N9-C4), where the energetic basis for an approximately linear delta-chi correlation can be seen: zeta (C3'-O3'-P-O5')-delta, in which the experimentally observed linear correlation between zeta and delta in DNA(220 degrees less than zeta less than 280 degrees) is clearly predicted; zeta-epsilon (C4'-C3'-O3'-P), which shows that epsilon increases with decreasing zeta less than 260 degrees; alpha (O3'-P-O5'-C5')-gamma (O5'-C5'-C4'-C3') where a clear linear correlation between these angles is also apparent, consistent with experiment; and several others. For the DNA molecule studied here, the sugar torsion delta is predicted to be the most flexible, while for the RNA molecule, the greatest amount of flexibility is expected to reside in alpha and gamma. Both the DNA and RNA molecules are predicted to be highly polymorphic. Complete energy minimization has been performed on each of the minima found in the energy searches and the results further support this prediction. Possible pathways for B-form to A-form DNA interconversion suggested by the results of this study are discussed. The results of these calculations support use of the new sugar modeling technique and torsion angle correlations in future conformational studies of nucleic acids.

A unique four-stranded model of a homologous recombination intermediate

This paper proposes a model of four-stranded DNA synapsis during recombination between homologous segments of two DNA duplexes. The proposed intermediate is one of only two known models having relative chain orientations about the synaptic junction that are consistent with recent topological results on the integrative recombination of bacteriophage lambda. This model has the advantage of providing a mechanism for recognition of sequence homology between duplexes through specific hydrogen-bond formation; other models are discussed in comparison. The new model is based on an alternative family of DNA structures having chain directions opposite to those of the Watson-Crick family of structures. Idealized coordinates for generating both right- and left-handed forms of these alternative structures are presented for further study.

The flexibility of the nucleic acids: (I). "SIR", a novel approach to the variation of polymer geometry in constrained systems

A novel and powerful methodology is developed which allows the alteration of molecular structures subjected to constraints and its application to polynucleotides with mononucleotide repeat symmetry, including the treatment of the flexible sugar rings is described. In contrast to procedures proposed by other authors, the constraints are formulated as differential equations which are linear with respect to the differentials of the geometrical variables. These equations can be solved easily by stepwise numerical integration involving sucessive infinitesimal rotations (SIR). Moreover, these equations define a set of independent curvilinear coordinates which can be used directly as the parameters of the energy functional in an energy minimisation procedure. This methodology allows the scanning of the full configurational space of a complex macromolecule, with direct access to the helicoidal variables in the case of periodic systems. Through this approach many problems involving biomacromolecular conformation, which would otherwise be intractable, may be studied with considerable ease.

Three-state diagram for DNA

Experimental phase diagrams (A form, B form, Coil) were built in the coordinates (a, alcohol fraction: T, temperature) for the natural DNAs and poly d(A-T). The main parameter of the B-A transition,-cooperativity length, nu o, was estimated by the slopes of the branches A-B, A-Coil, B-Coil in the vicinity of the triple point: nu o = 10-20 base pairs, which corresponds to the energy for the B/A junction of 1.2-1.8 kcal/mol. We discovered two new effects which are due to the coexistence of the three different conformations in one polymeric molecule: an increase in the melting temperature above that for the 'ideal' triple point (i.e. for the case of the ideal phase transitions); a widening of the melting curve within the B-A transition range.

The conformation of the complex of the helix destabilizing protein GP32 of bacteriophage T4 and single stranded DNA

The conformation of single stranded polynucleotides is changed specifically upon binding of the helix destabilizing protein of bacteriophage T4 (GP32). On the basis of circular dichroism (CD) and absorption experiments it is shown that denaturing conditions and the binding of oligopeptides can not induce the altered conformation. On the contrary, according to the current CD and absorption theory, the optical properties of the complex can be explained by a specific, regular conformation, characterized by an appreciable tilt of the bases (less than or equal to -10 degrees) and either a small rotation per base or a small helix diameter. This conformation agrees nicely with the increase of the base-base distance in the complex as determined in solution by electric field induced birefringence measurements. Our calculations show that also the model proposed by Alma (Ph.D. Thesis Catholic University Nijmegen, The Netherlands (1982)) for the complex of the helix destabilizing protein of bacteriophage fd, in which the helix diameter is large and the bases are almost parallel to the helix axis, would agree with the CD- and absorption spectra of the GP32-complex. For the latter protein this model would have to be modified with regard to the axial increment of the bases which is much larger in the GP32-complexes.

Flexibility of 3',5' deoxyribonucleoside diphosphates

In 3',5' deoxyribonucleoside diphosphates, in addition to the nature of the base and the sugar puckering, there are six single bond rotations. However, from the analysis of crystal structure data on the constituents of nucleic acids, only three rotational angles, that are about glycosyl bond, about C4'-C5' and about C3'-O3' bonds, are flexible. For a given sugar puckering and a base, potential energy calculations using non-bonded, electrostatic and torsional functions were carried out by varying the three torsion angles. The energies are represented as isopotential energy surfaces. Since the availability of the real-time color graphics, it is possible to analyse these isopotential energy surfaces. The calculations were carried out for C3' exo and C3' endo puckerings for deoxyribose and also for four bases. These calculations throw more light not only on the allowed regions for the three rotational angles but also on the relationships among them. The dependence of base and the puckering of the sugar on these rotational angles and thereby the flexibility of the 3',5' deoxyribonucleoside diphosphates is discussed. From our calculations, it is now possible to follow minimum energy path for interconversion among various conformers.

Sequence-dependent anisotropic flexibility of B-DNA. A conformational study

Bending flexibility of the six tetrameric duplexes was investigated d(AAAA):d(TTTT), d(AATT)2, d(TTAA)2, d(GGGG):d(CCCC), d(GGCC)2 and d(CCGG)2,. The tetramers were extended in the both directions by regular double helices. The stiffness of the B-DNA double helix when bent into the both grooves proved to be less than that in the perpendicular direction by an order of magnitude. Such an anisotropy is a property of the sugar-phosphate backbone structure. The calculated fluctuations of the DNA bending along the dyad axis, 5-7 degree, are in agreement with experimental value of the DNA persistence length. Anisotropy of the double helix is sequence-dependent: most easily bent into the minor groove are the tetramers with purine-pyrimidine dimer (RY) in the middle. In contrast, YR dinucleotides prefer bending into the major groove. Moreover, they have an equilibrium bend of 6-12 degree into this groove. The above inequality is caused by stacking interaction of the bases. The bend in the central dimer is distributed to some extent between the adjacent links, though the main fraction of the bend remains within the central link. Variation of the sugar-phosphate geometry in the bent helix is inessential, so that DNA remains within the B-family of forms: namely, when the helical axis is bent by 20 degree. the backbone dihedral angles vary by no more than 15 degree. The obtained results are in accord with x-ray structure of the B-DNA dodecamer; they further substantiate our early model of DNA wrapping in the nucleosome by means of "mini-kinks" separated by a half-pitch of the double helix, i.e. by 5-6 b.p. Sequence-dependent anisotropy of DNA presumably dictates the three-dimensional structure of DNA in solution as well. We have found that nonrandom allocation of YR dimers leads to the systematic bends in equilibrium structure of certain DNA fragments.

A study of DNA melting in concentrated water-alcohol solutions

The DNA melting profiles with high resolution have been studied for conditions corresponding to the B and A conformations of DNA in water-alcohol solutions. The melting profiles of the A-form and B-form DNA, their mean melting temperatures and melting range width were found to differ. DNA was shown to be heterogeneous in respect of the B-A transition, the GC-rich regions more readily converting into the A form than AT-rich ones. The presence of boundaries between the A and B sections within the transition zone did not smooth off the fine structure of melting profiles.

Aspects of specific DNA-protein interaction; local bending of DNA molecules by in-register binding of the oligopeptide antibiotic distamycin

Distamycin A binds strongly/moderately to DNA below/above r = 0.08 molecules bound per DNA phosphate. Titration viscometric measurements for high and low molecular weight DNA yielded the relative changes of DNA persistence length, delta alpha/alpha 0, and contour length, delta L/L0 (Nucl. Acids Res. 7 (1979) 1375). delta L/L0 is negligible/positive in the range below/above r = 0.09 at 0.2 M Na+! A two-line covering of the small groove by ligand molecules explains the increase of contour length. The characteristic delta alpha/alpha 0 drop is quantitatively interpreted by local DNA bending (kinking). The underlying theoretical basis is presented in two appendices and applied, in a third one, to literature data for the DNA-actinomycin system. The angle gamma of local DNA bending as induced by complex formation with different distamycin derivatives is presented for DNA species of different base composition. By means of appropriate model, a length mismatch per dinucleotide of (0.03(2) +/- 0.01)nm [or (0.04(3) +/- 0.01)nm] was derived from the experimentally obtained bending angle per dinucleotide of (1.6 +/- 0.4)0 [or (2.1 +/- 0.05)0], independent of DNA base composition and distamycin chain length.

Conformational transitions in closed circular DNA molecules. I. Topological and energetical considerations

A theory of conformational transitions in closed circular DNA as a function of topological linking number of the molecule (alpha) is elaborated taking into account topological and energetical considerations. The theory predicts a step-like dependence of a number of superhelical turns in DNA molecules (tau) on delta alpha. Thus, the number of superhelical turns tau = delta alpha for small values of delta alpha. For a large delta alpha (when conformational begin to occur) tau = delta alpha - sigma phi ij, where sigma phi ij is the total angle of conformational transitions for a given delta alpha. This prediction is in good agreement with published data on the dependence of the sedimentation coefficient of circular DNA molecules on their topological linking number. The results also allow to explain the disagreement between a number of titratable superhelical turns in circular DNA molecules and a number of supercoiles seen on electron micrographs for molecules with sufficiently large delta alpha.

Circular dichroism anisotrophy of DNA with different modifications at N7 of guanine

The complexex DNA-Ag1+, DNA-Cu1+, protonated DNA and DNA methylated at N7 of guanine were oriented by pumping the solutions through a multicapillary cell in the direction of a light beam. The CD components along the DNA axis, delta epsilon parallel, and normal to it, 2 delta epsilon perpendicular, were calculated from the CD spectra of the oriented samples by the method of Chung and Holzwarth, (1975) J. Mol. Biol. 92, 449--466. It was shown that in most cases, except that of the protonated DNA, the degree of orientation was only slightly less than that for pure DNA. This demonstrated the absence of aggregation and of appreciable denaturation. In all cases the modifications of DNA give rise to a negative component 2 delta epsilon perpendicular, whose magnitude increased as the extent of modification increased. From both the CD spectra of non-oriented samples and the absorption spectra, an inference is drawn that Ag1+ and Cu1+ are attached to the same site as CH3 groups i.e., to the N7 atom of guanine. Proton transfer along the H-bond from the N1 atom of G to the N3 atom of the complementary cytosine is suggested to be a result of the modifications, although the case of H+-DNA may differ from the others. Based on the CD spectra for the anisotropic components, delta epsilon parallel and 2 delta epsilon perpendicular, it is proposed that ligand binding is accompanied by winding of the DNA helix.

Statistical mechanics of supercoils and the torsional stiffness of the DNA double helix

The distribution of closed unknotted polymer chains over the writhing number is calculated by the Monte-Carlo method. For circular duplex DNA the variance of the distribution equals approximately half the observed variance of equilibrium distribution over the linking number. The balance which arises from fluctuations in DNA twisting makes it possible to estimate the torsional stiffness of the double helix.