Sugar-phosphate Backbone Conformation Research Articles

Why Are Left-Handed G-Quadruplexes Scarce?

G-quadruplexes (G4s) are nucleic acid structures crucial for the regulation of gene expression and genome maintenance. While they hold promise as nanodevice components, achieving desired G4 folds requires understanding the interplay between stability and structural properties, like helicity. Although right-handed G4 structures dominate the experimental data, the molecular basis for this preference over left-handed helicity is unclear. To address this, we employ all-atom molecular dynamics simulations and quantum chemical methods. Our results reveal that right-handed G4s exhibit greater thermodynamic and kinetic stability as a result of favorable sugar-phosphate backbone conformations in guanine tracts. Moreover, while hydrogen-bonding patterns influence helicity-specific G4 loop conformations, they minimally affect stability differences. We also elucidate the strong correlation between helicity and the strand progression direction, essential for G4 structures. These findings deepen our understanding of G4s, providing molecular-level insights into their structural and energetic preferences, which could inform the design of novel nanodevices.

8-Oxoguanine Affects DNA Backbone Conformation in the EcoRI Recognition Site and Inhibits Its Cleavage by the Enzyme.

8-oxoguanine is one of the most abundant and impactful oxidative DNA lesions. However, the reasons underlying its effects, especially those not directly explained by the altered base pairing ability, are poorly understood. We report the effect of the lesion on the action of EcoRI, a widely used restriction endonuclease. Introduction of 8-oxoguanine inside, or adjacent to, the GAATTC recognition site embedded within the Drew—Dickerson dodecamer sequence notably reduced the EcoRI activity. Solution NMR revealed that 8-oxoguanine in the DNA duplex causes substantial alterations in the sugar—phosphate backbone conformation, inducing a BI→BII transition. Moreover, molecular dynamics of the complex suggested that 8-oxoguanine, although does not disrupt the sequence-specific contacts formed by the enzyme with DNA, shifts the distribution of BI/BII backbone conformers. Based on our data, we propose that the disruption of enzymatic cleavage can be linked with the altered backbone conformation and dynamics in the free oxidized DNA substrate and, possibly, at the protein—DNA interface.

Design, synthesis, and duplex-stabilizing properties of conformationally constrained tricyclic analogues of LNA.

The design, synthesis and biophysical evaluation of two highly-constrained tricyclic analogues of locked nucleic acid (LNA), which restrict rotation around the C4'-C5'-exocyclic bond (torsion angle γ) and enhance hydrophobicity in the minor groove and along the major groove, are reported. A structural model that provides insights into the sugar-phosphate backbone conformations required for efficient hybridization to complementary nucleic acids is also presented.

Understanding the Sequence Preference of Recurrent RNA Building Blocks using Quantum Chemistry: The Intrastrand RNA Dinucleotide Platform.

Folded RNA molecules are shaped by an astonishing variety of highly conserved noncanonical molecular interactions and backbone topologies. The dinucleotide platform is a widespread recurrent RNA modular building submotif formed by the side-by-side pairing of bases from two consecutive nucleotides within a single strand, with highly specific sequence preferences. This unique arrangement of bases is cemented by an intricate network of noncanonical hydrogen bonds and facilitated by a distinctive backbone topology. The present study investigates the gas-phase intrinsic stabilities of the three most common RNA dinucleotide platforms - 5'-GpU-3', ApA, and UpC - via state-of-the-art quantum-chemical (QM) techniques. The mean stability of base-base interactions decreases with sequence in the order GpU > ApA > UpC. Bader's atoms-in-molecules analysis reveals that the N2(G)…O4(U) hydrogen bond of the GpU platform is stronger than the corresponding hydrogen bonds in the other two platforms. The mixed-pucker sugar-phosphate backbone conformation found in most GpU platforms, in which the 5'-ribose sugar (G) is in the C2'-endo form and the 3'-sugar (U) in the C3'-endo form, is intrinsically more stable than the standard A-RNA backbone arrangement, partially as a result of a favorable O2'…O2P intra-platform interaction. Our results thus validate the hypothesis of Lu et al. (Lu Xiang-Jun, et al. Nucleic Acids Res. 2010, 38, 4868-4876), that the superior stability of GpU platforms is partially mediated by the strong O2'…O2P hydrogen bond. In contrast, ApA and especially UpC platform-compatible backbone conformations are rather diverse and do not display any characteristic structural features. The average stabilities of ApA and UpC derived backbone conformers are also lower than those of GpU platforms. Thus, the observed structural and evolutionary patterns of the dinucleotide platforms can be accounted for, to a large extent, by their intrinsic properties as described by modern QM calculations. In contrast, we show that the dinucleotide platform is not properly described in the course of atomistic explicit-solvent simulations. Our work also gives methodological insights into QM calculations of experimental RNA backbone geometries. Such calculations are inherently complicated by rather large data and refinement uncertainties in the available RNA experimental structures, which often preclude reliable energy computations.

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

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

Conformational Characteristics and Correlations in Crystal Structures of Nucleic Acid Oligonucleotides: Evidence for Sub-states

Sugar phosphate backbone conformations are a structural element inextricably involved in a complete understanding of specific recognition nucleic acid ligand interactions, from early stage discrimination of the correct target to complexation per se, including any structural adaptation on binding. The collective results of high resolution DNA, RNA and protein/DNA crystal structures provide an opportunity for an improved and enhanced statistical analysis of standard and unusual sugar-phosphate backbone conformations together with corresponding dinucleotide sequence effects as a basis for further exploration of conformational effects on binding. In this study, we have analyzed the conformations of all relevant crystal structures in the nucleic acids data base, determined the frequency distribution of all possible ϵ, ζ, α, β and γ backbone angle arrangements within four nucleic acid categories (A-RNA and A-DNA, free and bound B-DNA) and explored the relationships between backbone angles, sugar puckers and selected helical parameters. The trends in the correlations are found to be similar regardless of the nucleic acid category. It is interesting that specific structural effects exhibited by the different unusual backbone sub-states are in some cases contravariant. Certain α/γ changes are accompanied by C3′ endo (north) sugars, small twist angles and positive values of base pair roll, and favor a displacement of nucleotide bases towards the minor groove compared to that of canonical B form structures. Unusual ϵ/ζ combinations occur with C2′ (south) sugars, high twist angles, negative values of base pair roll, and base displacements towards the major groove. Furthermore, any unusual backbone correlates with a reduced dispersion of equilibrium structural parameters of the whole double helix, as evidenced by the reduced standard deviations of almost all con- formational parameters. Finally, a strong sequence effect is displayed in the free oligomers, but reduced somewhat in the ligand bound forms. The most variable steps are GpA and CpA, and, to a lesser extent, their partners TpC and TpG. The results provide a basis for considering if the variable and non-variable steps within a biological active sequence precisely determine morphological structural features as the curvature direction, the groove depth, and the accessibility of base pair for non covalent associations.

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations.

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)(2) and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar-phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving 'lock and key' recognition.

Direct and indirect readout in mutant Met repressor–operator complexes

Background: The methionine repressor, MetJ, represses the transcription of genes involved in methionine biosynthesis by binding to arrays of two to five adjacent copies of an eight base-pair ‘metbox’ sequence. Naturally occurring operators differ from the consensus sequence to a greater extent as the number of metboxes increases. MetJ, while accommodating this sequence variation in natural operators, is very sensitive to particular base changes, even where bases are not directly contacted in the crystal structure of a complex formed between the repressor and consensus operator. Results: Here we report the high-resolution structure of a MetJ mutant, Q44K, bound to the consensus operator sequence (Q44Kwt19) and two related sequences containing mutations at sites believed to be important for indirect readout at non-contacted bases. The overall structure of the Q44Kwt19 complex is very similar to the wild-type complex, but there are small variations in sugar–phosphate backbone conformation and direct contacts to the DNA bases. The mutant complexes show a mixture of direct and indirect readout of sequence variations, with differences in direct contacts and DNA conformation. Conclusions: Comparison of the wild-type and mutant repressor–operator complexes shows that the repressor makes sufficiently strong interactions with the sugar–phosphate backbone to accommodate some variation in operator sequence with minor changes in direct bases contacts. The reduction in repressor affinity for the two mutant repressor complexes can be partially attributed to a loss in direct contacts to the DNA. In one case, however, the replacement of a flexible TA base-step leads to an unfavourable DNA conformation that reduces the stability of the repressor–operator complex.

Crystal structures of two plasmid copy control related RNA duplexes: An 18 base pair duplex at 1.20 A resolution and a 19 base pair duplex at 1.55 A resolution.

The structures of two RNA duplexes, whose sequences correspond to portions of the ColE1 plasmid copy control RNA I and RNA II, have been determined. Crystals containing the 18mers 5'-CA CCGUUGGUAGCGGUGC-3' and 5'-CACCGCUACCAACGGUGC-3' diffract to 1.20 A resolution while those containing the 19mers 5'-GCACCGUUGGUAGCGGUGC-3' and 5'-GCACCGCUACCAACGGUGC-3' diffract to 1.55 A resolution. Both duplexes are standard A form, with Watson-Crick base pairing throughout. Use of anisotropic atomic displacement factors in refinement of the 1.20 A structure dramatically improved refinement statistics, resulting in a final R(free) of 15.0% and a crystallographic R-factor of 11.6%. Perhaps surprisingly, these crystals of the 18 base pair RNA exhibit a 36-fold static disorder, resulting in a structure with a single sugar-phosphate backbone conformation and an averaged base composition at each residue. Since the sugar-phosphate backbone structure is identical in the 36 different nucleotides that are superimposed, there can be no sequence-dependent variation in the structure. The average ribose pucker amplitude is 45.8 degrees for the 18 base pair structure and 46.4 degrees for the 19 base pair structure; these values are respectively 19% and 20% larger than the average pucker amplitude reported from nucleoside crystal structures. A standard RNA water structure, based on analysis of the hydration of these crystal structures and that of the TAR RNA stem [Ippolito, J. A., and Steitz, T. A. (1998) Proc. Natl. Acad. Sci. U.S.A. 95, 9819-9824], has been derived, which has allowed us to predict water positions in lower resolution RNA crystal structures. We report a new RNA packing motif, in which three pro-S(p) phosphate oxygens interact with an ammonium ion.

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 Å per base-pair step to highly extended forms of 7.0 Å 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 Information Transmission from the Nucleobase Drives the Sugar-Phosphate Backbone Conformation in the Nucleotide Wire

We herein present our results showing that a simple mononucleotide block actually acts as a wire. This is evidenced by the fact that any change of the electronic character of the nucleobase is transmitted to drive the sugar conformation, which in turn dictates cooperatively the preferred orientation of the phosphate backbone.

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.

Influence of electrostatic interactions on the sugar-phosphate backbone conformation in DNA.

The influence of sugar ring flexibility in DNA on the mechanism of the B <===> A conformational transition is studied. The dipole moment of the deoxyribose as a function of its puckered states is calculated by the quantum-mechanical method using the MINDO/3 approximation. The interaction of the sugar dipole with the neighbor molecular groups in polynucleotide chains is estimated. The sugar dipole interaction with phosphate groups and counterions is shown to be strong and capable to deform the pseudorotation potential of deoxyribose. The effective pseudorotation potential of sugar ring in the B and A helices is obtained. The results are used to explain the behavior of Raman bands in the region of sugar-phosphate vibrations. The mechanism of the effect of electrostatic forces on the sugar-phosphate backbone conformation, which is essential for the B <===> A and other structure transitions, is offered.

Distance‐dependent dielectric constants and their application to double‐helical DNA

Detailed studies of structures of biological macromolecules, even in simplified models, involve many costly and time-consuming calculations. Any thorough methods require sampling of an extremely large conformation and momentum space. Calculations of electrostatic interactions, which depend on many physical factors, such as the details of solvent, solvent accessibility in macromolecules, and molecular polarizability, are always developed in a compromise between more rigorous, detailed models and the need for immediate application to complicated biological systems. In this paper, a middle ground is taken between the more exact theoretical models and the simplest constant values for the dielectric constant. The effects of solvent, counterions, and molecular polarizability are incorporated through a set of adjustable parameters that should be determined from experimental conditions. Several previous forms for the dielectric function are compared with the new ones. The present methods use Langevin functions to span the region of dielectric constant between bulk solvent and cavity values. Application of such dielectric models to double-helical DNA is important because base-stacking preferences were previously demonstrated [A. Sarai, J. Mazur, R. Nussinov, and R. L. Jernigan (1988) Biochemistry, vol. 27, pp. 8498-8502] to be sensitive to the electrostatic formulation. Here we find that poly(dG).poly(dC) can be A form for high screening and B form for low screening. By contrast, poly(dA).poly(dT) can only take the B form. Base stacking is more sensitive to the form of the dielectric function than are the sugar-phosphate backbone conformations. Also in B form, the backbone conformations are not so affected by the base types as in A form.

Structural variation in d(CTCTAGAG). Implications for protein-DNA interactions

Single-crystal X-ray diffraction techniques have been used to characterize the structure of the self-complementary DNA oligomer d(CTCTAGAG). The structure was refined to an R factor of 14.7% using data to 2.15-A resolution. The tetragonal unit cell, space group P4(3)2(1)2, has dimensions a = 42.53 and c = 24.33 A. The asymmetric unit consists of a single strand or four base pairs. Two strands, related by a crystallographic dyad axis, coil about each other to form a right-handed duplex. This octamer duplex has a mean helix rotation of 32 degrees, 11.3 base pairs per turn, an average rise of 3.1 A, C3'-endo furanose conformations, a shallow minor groove, and a deep major groove. Such averaged parameters suggest classification of the octamer as a member of the A-DNA family. However, the global parameters tend to mask variations in conformational parameters observed at the level of the base pairs. In particular, the central TpA (= TpA) step displays extensive interstrand purine-purine overlap and an unusual sugar-phosphate backbone conformation. These structural features may be directly related to certain sequence-specific protein-DNA interactions involving nucleases and repressors.

NMR study of stacking interactions and conformational adjustments in the dinucleotide-carcinogen adduct 2'-deoxycytidylyl-(3----5)-2'-deoxy-8-(N-fluoren-2-ylac etamido)guanosi ne.

The conformation and dynamics of the dinucleotide d-CpG modified at the C(8) position of the guanine ring by the carcinogen 2-(acetylamino)fluorene has been investigated by high-field 1H NMR spectroscopy. A two-state analysis of chemical shift data has enabled estimation of the extent of intramolecular stacking in aqueous solution as a function of temperature. The stacking, which is mostly fluorene-cytosine, is virtually complete in the low-temperature range. The 500-MHz 1H NMR spectrum consists of two subspectra near ambient temperatures due to a 14.3 +/- 0.3 kcal/mol barrier to internal rotation about the amide bond in the stacked form. A large barrier to internal rotation about the guanyl-nitrogen bond at C(8) has also been ascertained, but separate NMR subspectra were not detected due to the predominance of one of the torsional diastereomers (alpha' = 90 degrees) in the fully stacked state. Problems of self-association and chemical exchange were identified and overcome to enable analysis of the sugar-phosphate backbone conformation utilizing coupling constants. For the exocyclic C(4')-C(5') bond of the deoxyguanosine moiety, there is a high gauche+ (gamma = 60 degrees) conformer population, which is uncommon for a purine nucleotide with a syn orientation about the glycosyl bond. The gauche- conformation (gamma = 300 degrees), which is normally present in syn purine nucleotides in solution, was not detected. The exocyclic C(5')-O(5') torsion of the deoxy-guanosine moiety remains near the classical energy minimum (beta = 180 degrees) in the major stacked conformations. The sugar ring of the deoxycytidine moiety is predominantly in the C2'-endo conformation, while the deoxyguanosine ring is a mixture of conformations, one of which appears to be unusually puckered. The results support intercalation models of modified DNA and suggest a looped-out structure, with the modified guanine being the first base in the loop. Such structures could explain the relatively rapid rate of repair and the frame-shift mutations of this type of adduct.

A Proposal For A Specific Double-Helical Structure In Which The Polynucleotide Strands Intercalate Instead of Forming Base-pairs

Double helices, since the discovery of the DNA structure by Watson and Crick, represent the single most important secondary structural form of nucleic acids. The secondary structures of a variety of polynucleotide helices have now been well characterised with hydrogen-bonded base-pairs as building blocks. We wish to propose here the possibility, in a specific case, of a double stranded helical structure without any base-pair, but having a repeat unit of two nucleotides with their bases stacked through intercalation. The proposal comes from the initial models we have built for poly(dC) using the stacking patterns found in the crystal structures of 5'-dCMPNa2 which crystallises in two forms depending on the degree of hydration. These structures have pairs of nucleotides with the cytosine rings partially overlapping and separated by 3.3A. Using these as repeat units one could generate a model for poly(dC) with parallel strands, having a turn angle of 30 degrees and a base separation of 6.6A along each strand. Both right and left handed models with these parameters can be built in a smooth fashion without any obviously unreasonable stereochemical contacts. The helix diameter is about 13.5A, much smaller than that of normal helices with base-pair repeats. The changes in the sugar-phosphate backbone conformation in the present models compared to normal duplexes only reflect the torsional flexibility available for extension of polynucleotide chains as manifested by the crystal structures of drug-inserted oligonucleotide complexes. Intercalation proposed here could have some structural relevance elsewhere, for instance to the base-mismatched regions on the double helix and the packing of noncomplementary single strands as found in the filamentous bacteriophage Pf1.

Conformational studies on guanosine nucleotides and polynucleotides. The effect of the base on the glycosyl and backbone conformations.

AbstractPotential energy calculations have been carried out on guanosine 5′‐monophosphate and guanosine 3′,5′‐diphosphate by simultaneous variation of the rotation angles χ and ψ to assess the influence of the base on the glycosyl and the sugar–phosphate backbone conformations and their possible interdependence. The results obtained predict that guanine containing nucleotides tend to favor the syn glycosyl conformation in sharp contrast to the anti conformation found for adenine and the common pyrimidine nucleotides. This characteristic behavior of guanine nucleotides is largely a result of increased non‐bonded van der Waals' and electrostatic attractive interactions between the 5′‐phosphate and the amino group of the base. As seen previously [N. Yathindra and M. Sundaralingam, Biopolymers, 12, 297 (1973)] changes in the glycosyl conformation mainly affects the C(4′)–C(5′) bond rotations of the backbone. The most favored (χ,ψ) conformer corresponding to the global minimum both in 5′‐GMP and 3′,5′‐GDP is syn‐gg, (190°,80°) for C(2′)‐endo rings and (180°,50°) for C(3′)‐endo rings. It is noteworthy that the syn angle is different from the usually observed range of 210–260°. The next significant minimum in the (χ,ψ) energy surface corresponds to the syn‐gt conformational combination (the syn angle now in the normal range), indicating that the occurrence of the normal syn conformer induces a distortion in the backbone c(4′)‐c(5′) bond conformation from the preferred gg. These studies reveal the significant influence of the base not only on the glycosyl conformation but also on the backbone conformations of polynucleotides.