Mediation of the A/B-DNA helix transition by G-tracts in the crystal structure of duplex CATGGGCCCATG

The crystal structure of the DNA dodecamer duplex CATGGGCCCATG lies on a structural continuum along the transition between A- and B-DNA. The dodecamer possesses the normal vector plot and inclination values typical of B-DNA, but has the crystal packing, helical twist, groove width, sugar pucker, slide and x-displacement values typical of A-DNA. The structure shows highly ordered water structures, such as a double spine of water molecules against each side of the major groove, stabilizing the GC base pairs in an A-like conformation. The different hydration of GC and AT base pairs provides a physical basis for solvent-dependent facilitation of the A↔B helix transition by GC base pairs. Crystal structures of CATGGGCCCATG and other A/B-DNA intermediates support a ‘slide first, roll later’ mechanism for the B→A helix transition. In the distribution of helical parameters in protein–DNA crystal structures, GpG base steps show A-like properties, reflecting their innate predisposition for the A conformation.     

Invariant and variable base stacking geometries in B-DNA and A-DNA

We calculated the interatomic distances between all couples of non-hydrogen atoms belonging to the neighboring Watson-Crick base pairs in the available crystal structures of DNA. Their standard deviations revealed remarkably large differences in the variability of the base stacking geometries of the particular steps. In line with experimental studies in solution, (CpA)-(TpG) and (TpA).(TpA) were identified as the most variable or flexible steps in the crystal structures of B-DNA. On the other hand, base stacking geometries of the (ApT).(ApT) steps were the most invariant, which was very surprising because all three steps composed only of C and G were much more flexible. This finding suggests that conformational stability of DNA and the rigidity have different origins. Furthermore, the nucleotide sequence dependence of the flexibility was almost reversed in A-DNA because the most flexible steps in B-DNA were the least flexible in A-DNA. The most invariant steps of B-DNA were variable in A-DNA. The (ApT).(ApT) step was a notable exception to this rule because it belonged to the most rigid steps in both B-DNA and A-DNA. The present results are fully consistent with the properties that poly(dA-dT).poly(dA-dT), poly(dA).poly(dT), poly(dAdC).poly(dG-dT) and poly(dA-dG).poly(dC-dT) exhibit in solution.     

The post-SCF quantum chemistry characteristics of inter- and intra-strand stacking interactions in d(CpG) and d(GpC) steps found in B-DNA, A-DNA and Z-DNA crystals

The energies of intra- and inter-strand stacking interactions in model d(GpC) and d(CpG) two-base-pair steps were estimated by MP2/aug-cc-pVDZ single point calculations corrected for basis superposition errors. The stacked two-nucleobase pairs were constructed using experimental values of base pair and base step parameters taken from Nucleic Acid Database (). Three distinct polymorphic forms were analysed, namely A-, B- and Z-DNA. The applied methodology enables statistical analysis of structural and energetic diversities. The structural relationships between polymorphic forms are quite complex and depend on the sequence of pairs. The variability of parameters such as shift and tilt is almost the same irrespective of the polymorphic form and sequence of steps analysed. In contrast, shift and twist distributions easily discriminate all three polymorphic forms of DNA. Interestingly, despite significant structural diversities, the energies of the most frequent energy ranges are comparable irrespective of the polymorphic form and base sequence. There was observed compensation of inter- and intra-strand interactions, especially for d(GpC) and d(CpG) steps found in A- and B-DNA. Thus, among many other roles, these pairs act as a kind of energetic buffer, balancing the double helix. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Solution structure of the dodecamer d-(CATGGGCC-CATG)2 is B-DNA. Experimental and molecular dynamics study

The DNA duplex d-(CATGGGCCCATG)2 has been studied in solution by FTIR, NMR and CD. The experimental approaches have been complemented by series of large-scale unrestrained molecular dynamics simulation with explicit inclusion of solvent and counterions. Typical proton-proton distances extracted from the NMR spectra and the CD spectra are completely in agreement with slightly modified B-DNA. By molecular dynamics simulation, starting from A-type sugar pucker, a spontaneous repuckering to B-type sugar pucker was observed. Both experimental and theoretical approaches suggest for the dodecamer d-(CATGGGCCCATG)2 under solution conditions puckering of all 2'-deoxyribose residues in the south conformation (mostly C2'-endo) and can exclude significant population of sugars in the north conformation (C3'-endo). NMR, FTIR and CD data are in agreement with a B-form of the dodecamer in solution. Furthermore, the duplex shows a cooperative B-A transition in solution induced by addition of trifluorethanol. This contrasts a recently published crystal structure of the same oligonucleotide found as an intermediate between B- and A-DNA where 23 out of 24 sugar residues were reported to adopt the north (N-type) conformation (C3'-endo) like in A-DNA (Ng, H. L., Kopka, M. L. and Dickerson, R. E., Proc. Natl. Acad. Sci. U S A 97, 2035-2039 (2000)). The simulated structures resemble standard B-DNA. They nevertheless show a moderate shift towards A-type stacking similar to that seen in the crystal, despite the striking difference in sugar puckers between the MD and X-ray structures. This is in agreement with preceding MD reports noticing special stacking features of G-tracts exhibiting a tendency towards the A-type stacking supported by the CD spectra also reflecting the G-tract stacking. MD simulations reveal several noticeable local conformational variations, such as redistribution of helical twist and base pair roll between the central GpC steps and the adjacent G-tract segments, as well as a substantial helical twist variability in the CpA(TpG) steps combined with a large positive base pair roll. These local variations are rather different from those seen in the crystal.     

Energy hyperspace for stacking interaction in AU/AU dinucleotide step: Dispersion‐corrected density functional theory study

Double helical structures of DNA and RNA are mostly determined by base pair stacking interactions, which give them the base sequence‐directed features, such as small roll values for the purine–pyrimidine steps. Earlier attempts to characterize stacking interactions were mostly restricted to calculations on fiber diffraction geometries or optimized structure using ab initio calculations lacking variation in geometry to comment on rather unusual large roll values observed in AU/AU base pair step in crystal structures of RNA double helices. We have generated stacking energy hyperspace by modeling geometries with variations along the important degrees of freedom, roll, and slide, which were chosen via statistical analysis as maximally sequence dependent. Corresponding energy contours were constructed by several quantum chemical methods including dispersion corrections. This analysis established the most suitable methods for stacked base pair systems despite the limitation imparted by number of atom in a base pair step to employ very high level of theory. All the methods predict negative roll value and near‐zero slide to be most favorable for the purine–pyrimidine steps, in agreement with Calladine's steric clash based rule. Successive base pairs in RNA are always linked by sugar–phosphate backbone with C3′‐endo sugars and this demands C1′–C1′ distance of about 5.4 Å along the chains. Consideration of an energy penalty term for deviation of C1′–C1′ distance from the mean value, to the recent DFT‐D functionals, specifically ωB97X‐D appears to predict reliable energy contour for AU/AU step. Such distance‐based penalty improves energy contours for the other purine–pyrimidine sequences also. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 107–120, 2014.     

Characterization of the base stacking interactions in DNA by means of Lennard-Jones empirical potentials

Three empirical potentials of the Lennard-Jones type taken from literature were used to calculate van der Waals contributions to the base-pair couples stacking energies in B-DNA and A-DNA type double helical conformations. The information obtained can be summarized as follows: (1) Purine-pyrimidine and purine-purine (pyrimidine-pyrimidine in the complementary strand) sequences preferred right-handed helical arrangement, whereas pyrimidine-purine sequences favoured left-handed (C-G) or unwound (T-A) stacking geometry; in the latter case this only held for B- but not A-DNA (the C-G sequence was not studied in A-DNA owing to difficulties (see below) with the G amino group in B-DNA); (2) Positive propeller twist of base-pairs was stable in both B- and A-DNA; the thymine methyl group promoted the propeller and this effect was strongest in the A-T step; (3) Tilt of base pairs occurred around zero in B-DNA and between 15-20 degrees C in A-DNA, in agreement with the experimental observations; (4) Vertical separation of base pairs was optimal within 0.33-0.34 nm for B-DNA and around 0.29 nm for A-DNA using the 9-6 potential. The 12-6 potential gave similar results with B-DNA as the 9-6 potential if, however, base pairs were separated by 0.35-0.36 nm; (5) The calculated effect of the guanine amino group was substantially stronger than expected on the basis of data derived from X-ray diffraction studies of oligonucleotide single crystals; (6) In comparison with the 9-6 potential, the 12-6 potential provided more strict energy minima. In summary, the empirical potentials reproduce, at least semiquantitatively, many but not all DNA properties; this should be taken into account whenever the potentials are used for prediction purposes.     

Base pair buckling can eliminate the interstrand purine clash at the CpG steps in B-DNA caused by the base pair propeller twisting

Results of calculations using various empirical potentials suggest that base pair buckling, which commonly occurs in DNA crystal structures, is sufficient to eliminate the steric clash at CpG steps in B-DNA, originating from the base pair propeller twisting. The buckling is formed by an inclination of cytosines while deviations of guanines from a plane perpendicular to the double helix axis are unfavorable. The buckling is accompanied by an increased vertical separation of the base pair centers but the buckled arrangement of base pairs is at least as stable as when the vertical separation is normal and buckle zero. In addition, room is created by the increased vertical separation for the bases to propeller twist as is observed in DNA crystal structures. Further stabilization of base stacking is introduced into the buckled base pair arrangement by roll opening the base pairs into the double helix minor groove. The roll may lead to the double helix bending and liberation of guanines from the strictly perpendicular orientation to the double helix axis. The liberated guanines further contribute to the base pair buckling and stacking improvement. This work also suggests a characteristic very stable DNA structure promoted by nucleotide sequences in which runs of purines follow runs of pyrimidine bases.     

The stacking interactions of bipyridine complexes: the influence of the metal ion type on the strength of interactions

The strength of the stacking interactions in the bipy complexes of nickel, palladium, and platinum, [M(CN)₂ bipy]₂ (M?=?Ni, Pd, Pt), was calculated using the ωB97xD/def2-TZVP method. The results show that for all considered geometries, interactions are the strongest for platinum, and weakest for nickel complexes, as a result of higher dispersion contributions of platinum over the palladium and nickel complexes. It was also shown that strength of interactions considerably rises with an increase of the stacking overlap area. As a consequence of the favorable electrostatic term, the strength of interactions also rises when metal atom and cyano ligands are involved in the overlap with bipy ligand. The strongest interaction was calculated in the platinum complex, for the geometry that has overlap of metal and cyano ligands with bipy ligand with an energy of -39.80 kcal mol^-1. The energies for similar geometries of palladium and nickel complexes are -34.60 and -32.45 kcal mol^-1. These energies, remarkably, exceed the strength of the stacking interactions between organic aromatic molecules. These results can be of importance in all systems with stacking interactions, from materials to biomolecules.     

Finding and visualizing nucleic acid base stacking

Base stacking is one of the primary factors stabilizing nucleic acid structure. Yet, methods for locating stacking interactions in DNA and RNA are rare and methods for displaying stacking are rarer still. We present here simple, automated procedures to search nucleic acid molecules for base-base and base-oxygen stacking and to display these interactions graphically in a manner that readily conveys both the location and the quality of the interaction. The method makes no a priori assumptions about relative base positions when searching for stacking, nor does it rely on empirical energy functions. This is a distinct advantage for two reasons. First, the relative contributions of the forces stabilizing stacked bases are unknown. Second, the electrostatic and hydrophobic components of base stacking are both poorly defined by existing potential energy functions.     

Index to Authors,Volume 3

Abstract

Results of calculations using various empirical potentials suggest that base pair buckling, which commonly occurs in DNA crystal structures, is sufficient to eliminate the steric clash at CpG steps in B-DNA, originating from the base pair propeller twisting. The buckling is formed by an inclination of cytosines while deviations of guanines from a plane perpendicular to the double helix axis are unfavorable. The buckling is accompanied by an increased vertical separation of the base pair centers but the buckled arrangement of base pairs is at least as stable as when the vertical separation is normal and buckle zero. In addition, room is created by the increased vertical separation for the bases to propeller twist as is observed in DNA crystal structures. Further stabilization of base stacking is introduced into the buckled base pair arrangement by roll opening the base pairs into the double helix minor groove. The roll may lead to the double helix bending and liberation of guanines from the strictly perpendicular orientation to the double helix axis. The liberated guanines further contribute to the base pair buckling and stacking improvement. This work also suggests a characteristic very stable DNA structure promoted by nucleotide sequences in which runs of purines follow runs of pyrimidine bases.     

Double helix conformation, groove dimensions and ligand binding potential of a G/C stretch in B-DNA.

The self-complementary DNA fragment CCGGCGCCGG crystallizes in the rhombohedral space group R3 with unit cell parameters a = 54.07 A and c = 44.59 A. The structure has been determined by X-ray diffraction methods at 2.2 A resolution and refined to an R value of 16.7%. In the crystal, the decamer forms B-DNA double helices with characteristic groove dimensions: compared with B-DNA of random sequence, the minor groove is wide and deep and the major groove is rather shallow. Local base pair geometries and stacking patterns are within the range commonly observed in B-DNA crystal structures. The duplex bears no resemblance to A-form DNA as might have been expected for a sequence with only GC base pairs. The shallow major groove permits an unusual crystal packing pattern with several direct intermolecular hydrogen bonds between phosphate oxygens and cytosine amino groups. In addition, decameric duplexes form quasi-infinite double helices in the crystal by end-to-end stacking. The groove geometries and accessibilities of this molecule as observed in the crystal may be important for the mode of binding of both proteins and drug molecules to G/C stretches in DNA.     

Hybrid simulation approach incorporating microscopic interaction along with rigid body degrees of freedom for stacking between base pairs

Stacking interaction between the aromatic heterocyclic bases plays an important role in the double helical structures of nucleic acids. Considering the base as rigid body, there are total of 18 degrees of freedom of a dinucleotide step. Some of these parameters show sequence preferences, indicating that the detailed atomic interactions are important in the stacking. Large variants of non‐canonical base pairs have been seen in the crystallographic structures of RNA. However, their stacking preferences are not thoroughly deciphered yet from experimental results. The current theoretical approaches use either the rigid body degrees of freedom where the atomic information are lost or computationally expensive all atom simulations. We have used a hybrid simulation approach incorporating Monte‐Carlo Metropolis sampling in the hyperspace of 18 stacking parameters where the interaction energies using AMBER‐parm99bsc0 and CHARMM‐36 force‐fields were calculated from atomic positions. We have also performed stacking energy calculations for structures from Monte‐Carlo ensemble by Dispersion corrected density functional theory. The available experimental data with Watson–Crick base pairs are compared to establish the validity of the method. Stacking interaction involving A:U and G:C base pairs with non‐canonical G:U base pairs also were calculated and showed that these structures were also sequence dependent. This approach could be useful to generate multiscale modeling of nucleic acids in terms of coarse‐grained parameters where the atomic interactions are preserved. This method would also be useful to predict structure and dynamics of different base pair steps containing non Watson–Crick base pairs, as found often in the non‐coding RNA structures. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 212–226, 2016.     

Theoretical studies on protein-nucleic acid interactions. III. Stacking of aromatic amino acids with bases and base pairs of nucleic acids

Stacking of aromatic amino acids tryptophan (Trp), tyrosine (Tyr), phenylalanine (Phe), and histidine (His) with bases and base pairs of nucleic acids has been studied. Stacking energies of the amino acid–base (or base pair) complexes have been calculated by second-order perturbation theory. Our results show that, in general, the predominant contribution to the total stacking energy comes from the dispersion terms. In these cases, repulsion energy is greater than the sum of electrostatic and polarization energies. In contrast to this, interaction of histidine with the bases and base pairs is largely Coulombic in nature. The complexes of guanine with aromatic amino acids are more stable than the corresponding complexes of adenine. Among pyrimidines, cytosine forms the most stable complexes with the aromatic amino acids. The G · C base pair has the highest affinity with aromatic amino acids among various sets of base pairs. Optimized geometries of the stacked complexes show that the aromatic moieties overlap only partially. The heteroatom of one residue generally overlaps with the other aromatic moiety. There is a considerable degree of configurational freedom in the stacked geometries. The role of stacking in specific recognition of base sequences by proteins is discussed.     

Molecular Mechanics Studies of Dinucleoside Methylphosphonates: Influence of Methylphosphonate and its Chirality on the Phosphodiester Conformation

Abstract

Molecular mechanics studies are performed on single stranded as well as base paired forms of dinucleoside methylphosphonates comprising different base sequences for both the Sand R-isomers of methylphosphonate (MP). S-MP produces noticeable distortions in the geometry, locally at the phosphate center, and this enables the stereochemical feasibility of compact g^? g^? phosphodiester. Besides, it tends to perturb the conformations around the P- 03′ and glycosyl bonds, causing minor variations in stacking interactions. In single stranded dinucleosides, the gain in adjacent base stacking interaction energies seems to be sufficient to overcome the barrier to P-03′ bond rotation arising due to S-MP…sugar interaction, and this results in transition to a compact phosphodiester (B_I-type) from an initial extended phosphodiester (B_II-type) conformation. Such a thing seems rather difficult in base pair constrained duplexes. Dinucleosides with R-MP behave analogous to normal phosphate duplexes as the methyl group is away from the sugar. It is found that dinucleoside methylphosphonates are energetically less favoured than the corresponding dinucleoside phosphates mainly due to the depletion of contributions from electrostatic attractive interactions involving the base and sugar with the methylphosphonate consequent to the nonionic nature of the latter. Neither S-MP nor R-MP seem to significantly alter the stereochemistry of duplex structure.     

A free energy analysis of nucleic acid base stacking in aqueous solution.

This paper reports a theoretical study of the free energy contributions to nucleic acid base stacking in aqueous solution. Electrostatic interactions are treated by using the finite difference Poisson-Boltzmann method and nonpolar effects are treated with explicit calculation of van der Waals interactions and/or free energy-surface area relationships. Although for some pairs of bases there is a favorable Coulombic interaction in the stacked conformation, generally the net effect of electrostatic interactions is to oppose stacking. This result is caused by the loss of favorable base-solvent electrostatic interactions, that accompany the partial removal of polar atoms from water in the stacked conformation. Nonpolar interactions, involving the hydrophobic effect and enhancement of van der Waals interactions caused by close-packing, drive stacking. The calculations qualitatively reproduce the experimental dependence of stacking free energy on purine-pyrimidine composition.     

Interaction Between the Guanine Amino Group and the Adenine Six Membered Ring Stabilizes the Unusual Conformation of the CpA Step in B-DNA

Abstract

The CpA step is dramatically overwound in several B-DNA oligonucleotide crystal structures and its AT pair is substantially shifted towards the cytosine of the preceding base pair and towards the minor groove. We show using a geometrical analysis of the crystal data and empirical potential calculations that a strong interaction between the guanine amino group and the adenine six membered ring is responsible for the unique conformational properties of the CpA step.     

Structural studies of DNA fragments: the G.T wobble base pair in A, B and Z DNA; the G.A base pair in B-DNA

The crystal structures of five double helical DNA fragments containing non-Watson-Crick complementary base pairs are reviewed. They comprise four fragments containing G.T base pairs: two deoxyoctamers d(GGGGCTCC) and d(GGGGTCCC) which crystallise as A type helices; a deoxydodecamer d(CGCGAATTTGCG) which crystallises in the B-DNA conformation; and the deoxyhexamer d(TGCGCG), which crystallises as a Z-DNA helix. In all four duplexes the G and T bases form wobble base pairs, with bases in the major tautomer forms and hydrogen bonds linking N1 of G with O2 of T and O6 of G with N3 of T. The X-ray analyses establish that the G.T wobble base pair can be accommodated in the A, B or Z double helix with minimal distortion of the global conformation. There are, however, changes in base stacking in the neighbourhood of the mismatched bases. The fifth structure, d(CGCGAATTAGCG), contains the purine purine mismatch G.A where G is in the anti and A in the syn conformation. The results represent the first direct structure determinations of base pair mismatches in DNA fragments and are discussed in relation to the fidelity of replication and mismatch recognition.