Analysis of conformational parameters in nucleic acid fragments. I. Single crystals of nucleosides and nucleotides

A study has been undertaken of conformational parameters in single crystal structures of nucleosides and nucleotides using the techniques of helical conformational analysis. A "quasi-helix" was generated from the geometry of base-paired structures, using published data extracted from the Cambridge structural database. A total of 54 base-pairs were found in these structures, for each of which were calculated hydrogen bond parameters, propeller twist, buckle and C1'-C1' separations. These were analysed according to various classifications. Propeller twists are found to show a wide range of values and the magnitude of twist appears to be unrelated to hydrogen bond parameters or C1'-C1' separation. The values of the buckle parameter vary over a smaller range of values and are unrelated to propeller twist magnitude. There is found to be a greater tendency to form homo-base-pairs among compounds containing adenine bases.     

Analysis of conformational parameters in nucleic acid fragments. II. Co-crystal complexes of nucleic acid bases.

Studies have been made of conformational parameters in co-crystal complexes and compounds of nucleic acid bases in which there is the possibility of formation of hetero-base-pairs. Using published data extracted from the Cambridge structural database, a total of 37 base-pairs were found, of which 25 were hetero-pairs and 12 homo-pairs. These base-pairs were subject to analysis to reveal hydrogen bond parameters, propeller twist, buckle and C1'-C1' separation (or a similar parameter if C1' atoms were not present). Hetero-pairs were found to show larger twists than homo-pairs, the magnitude of twist being unrelated to hydrogen bond parameters or buckle value. The propeller twisting is less pronounced in these nucleic acid bases than in nucleosides, but still has a significant magnitude. Propeller twisting in hetero-pairs is found to be larger than in homo-pairs. Hetero-pairs appear to be formed preferentially in competitive situations.     

The propeller DNA conformation of poly(dA).poly(dT).

Physical properties of the DNA duplex, poly(dA).poly(dT) differ considerably from the alternating copolymer poly(dAT). A number of molecular models have been used to describe these structures obtained from fiber X-ray diffraction data. The recent solutions of single crystal DNA dodecamer structures with segments of oligo-A.oligo-T have revealed the presence of a high propeller twist in the AT regions which is stabilized by the formation of bifurcated (three-center) hydrogen bonds on the floor of the major groove, involving the N6 amino group of adenine hydrogen bonding to two O4 atoms of adjacent thymine residues on the opposite strand. Here we show that it is possible to incorporate the features of the single crystal analysis, specifically high propeller twist, bifurcated hydrogen bonds, and a narrow minor groove, as well as the close interstrand NMR signal between adenine HC2 and ribose HC1' of the opposite strand, into a model that is fully compatible with the diffraction data obtained from poly(dA).poly(dT).     

Studies on the cross stand hydrogen bonds in DNA double helices.

A systematic analysis has been carried out to examine all the stereochemically possible bifurcated hydrogen bonds including those of cross strand type between propeller twisted base pairs in DNA double helices by stereochemical considerations involving base pairs alone and by molecular mechanics studies on dimer and trimer duplexes. The results show that there are limited number of combinations of adjacent base pairs that would facilitate bifurcated cross-strand hydrogen bond (CSH). B-type helices concomitant with negative propeller twist seem to be more favored for the occurrence of CSH than canonical A-type helices because of slide in the latter. The results also demonstrate that helices with appropriate sequences may possess continuous run of these propeller twist driven cross strand hydrogen bonds indicating that they may in fact be considered as yet another general structural feature of DNA helices.     

Molecular dynamics simulations on parallel and antiparallel C.G*G triplexes

Molecular dynamics (MD) studies have been carried out on the Hoogsteen hydrogen bonded parallel and the reverse Hoogsteen hydrogen bonded antiparallel C.G*G triplexes. Earlier, the molecular mechanics studies had shown that the parallel structure was energetically more favourable than the antiparallel structure. To characterize the structural stability of the two triplexes and to investigate whether the antiparallel structure can transit to an energetically more favourable structure, due to the local fluctuations in the structure during the MD simulation, the two structures were subjected to 200ps of constant temperature vacuum MD simulations at 300K. Initially no constraints were applied to the structures and it was observed that for the antiparallel triplex, the structure showed a large root mean square deviation from the starting structure within the first 12ps and the N4-H41--O6 hydrogen bond in the WC duplex got distorted due to a high propeller twist and a moderate increase in the opening angle in the basepairs. Starting from an initial value of 30 degrees , helical twist of the average structure from this simulation had a value of 36 degrees , while the parallel structure stabilized at a twist of 33 degrees. In spite of the hydrogen bond distortions in the antiparallel triplex, it was energetically comparable to the parallel triplex. To examine the structural characteristics of an undistorted structure, another MD simulation was performed on the antiparallel triplex by constraining all the hydrogen bonds. This structure stabilized at an average twist of 33 degrees. In the course of the dynamics though the energy of the molecule - compared to the initial structure - improved, it did not become comparable to the parallel structure. Energy minimization studies performed in the presence of explicit water and counterions also showed the two structures to be equally favourable energetically. Together these results indicate that the parallel C.G*G triplex with Hoogsteen hydrogen bonds also represents a stereochemically and energetically favourable structure for this class of triplexes.     

Crystal structure of a B-DNA dodecamer containing inosine, d(CGCIAATTCGCG), at 2.4 A resolution and its comparison with other B-DNA dodecamers.

The crystal structure of the dodecamer, d(CGCIAATTCGCG), has been determined at 2.4 A resolution by molecular replacement, and refined to an R-factor of 0.174. The structure is isomorphous with that of the B-DNA dodecamer, d(CGCGAATTCGCG), in space group P2(1)2(1)2(1) with cell dimensions of a = 24.9, b = 40.4, and c = 66.4 A. The initial difference Fourier maps clearly indicated the presence of inosine instead of guanine. The structure was refined with 44 water molecules, and compared to the parent dodecamer. Overall the two structures are very similar, and the I:C forms Watson-Crick base pairs with similar hydrogen bond geometry to the G:C base pairs. The propeller twist angle is low for I4:C21 and relatively high for the I16:C9 base pair (-3.2 degrees compared to -23.0 degrees), and the buckle angles alter, probably due to differences in the contacts with symmetry related molecules in the crystal lattice. The central base pairs of d(CGCIAATTCGCG) show the large propeller twist angles, and the narrow minor groove that characterize A-tract DNA, although I:C base pairs cannot form the major groove bifurcated hydrogen bonds that are possible for A:T base pairs.     

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

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

Analysis of conformational parameters in nucleic acid fragments. IV. Intercalating drug complexes of very short chain oligonucleotides.

Studies have been made of base-pairing conformational parameters in single crystal structures of very short chain oligonucleotide structures complexed with drug molecules, using data extracted from the Cambridge structural database. The planar portion of the drug has a tendency to intercalate between two bases, utilising strong stacking interactions to stabilise the configuration. The effect of the existence of a formative backbone is seen in the high occurrence of standard base-pairing schemes and a consistent C1'-C1' separation, although the choice of compounds studied does tend to emphasise complementary pairing. In addition to the modulation of the general magnitude which is reduced from that in uncomplexed oligonucleotides, there appears to be some correlation of propeller twist value with the presence of planar groups sandwiching a base-pair. The average twist in such sandwiched pairs is lower than in any other group studied to date.     

Sequence effects on the propeller twist of base pairs in DNA helices.

Molecular mechanics calculations are performed on all the ten base pair steps (duplex dimers) and also a number of trimer and tetramer duplexes comprising them, in an attempt to systematically examine the possible base sequence effects on the magnitudes of propeller twists of base pairs at a given step. The analysis reveals that though propeller twist is a base pair property, it behaves very much like other base step parameters such as slide, roll, twist etc., Hence, it may be necessary to monitor the nature and variation of magnitudes of pt at a step. Calculations performed on 45 out of the 136 unique tetramer combinations involving all the ten unique base steps show that the difference in magnitudes of propeller twists of the base pairs of a given step has been found to be either steep or moderate depending on base pairs that flank the base step. These observations compare very well with the available experimental data. Tetramer sequences, wherein a base pair of a base step repeats in the same direction, exhibit a relatively steep difference in propeller twist at the step. Tetramers other than these exhibit moderate difference in propeller twist. Such sequences are broadly classified as type-I and type-II respectively. Practically all the tetrads considered in the study, excepting those with GT step and a few involving CG and GC steps, conform to the above classification.     

Studies On The Cross Strand Hydrogen Bonds In DNA Double Helices

Abstract

A systematic analysis has been carried out to examine all the stereochemically possible bifurcated hydrogen bonds including those of cross strand type between propeller twisted base pairs in DNA double helices by stereochemical considerations involving base pairs alone and by molecular mechanics studies on dimer and trimer duplexes. The results show that there are limited number of combinations of adjacent base pairs that would facilitate bifurcated cross- strand hydrogen bond (CSH). B-type helices concomitant with negative propeller twist seem to be more favored for the occurrence of CSH than canonical A-type helices because of slide in the latter. The results also demonstrate that helices with appropriate sequences may possess continuous run of these propeller twist driven cross strand hydrogen bonds indicating that they may infact be considered as yet another general structural feature of DNA helices.     

New insights into the structure of An tracts and B'-B' bends in DNA

Energy calculations suggest that the currently available NOE distance constraints for An tracts in DNA are incapable of distinguishing between structures with a narrowed minor groove arising from a large propeller twist with a small inclination or from a small propeller twist with a large negative inclination. Furthermore, analysis of published data, together with energy estimations, strongly argue against bifurcated hydrogen bonding between A and T residues being the cause of the anomalous structural properties of An tracts. A conformational analysis of the B'-B' junction has been performed in which a single variable base pair has been inserted between two regions of B' structure. We have calculated low-energy structures for AnGAn,AnCAn,AnTAn,AnCTn, and TnCAn duplexes, where the An and Tn tracts were fixed in the anomalous B' conformation. Upon optimization, all these structures were found to contain a pronounced roll-like bending into the major groove at the site of the insertion. The important factors in the formation of these B'-B' bends are the destruction of the B' conformation and the concomitant widening of the minor groove at the junction region in order to reduce minor groove interstrand base clashes and improve interstrand stacking energy. If the B' conformation has strong negative inclination, the improved intrastrand stacking energy also contributes to the bending. In calculations of duplexes with An and Tn tracts in the B conformation instead of B', the bending disappears.     

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)₂ 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.     

High propeller twist and unusual hydrogen bonding patterns from the MD simulation of (dG)6.(dC)6

A molecular dynamics (MD) study of (dG)6.(dC)6 including counter ions and 292 water molecules was made. The hydrogen bonding pattern and propeller twist angles for the mini-helix are reported as averages for times spanning 21-30, 31-40, 41-50, and 51-60 ps. The propeller twist angles range from 18 degrees to 38 degrees. Bifurcated and interstrand neighboring base (twisted) hydrogen bonding patterns were found.     

A stable antiparallel cytosine-thymine base pair occurring only at the end of a duplex.

Ultraviolet hyperchromicity experiments indicate that in DNA duplex formation, a C-T mismatch is destabilizing in the center of a duplex, but behaves as a stable base pair at the terminus of a duplex. The C-T base pair is thought to contain two hydrogen bonds, but has thermodynamic parameters (delta Ho and delta Go of dissociation) that are similar to a G-C base pair. AMBER molecular mechanics calculations were performed to study the possible structural properties of DNA duplexes with central and terminal C-T combinations. These calculations also indicate that a central C-T pair destabilizes a duplex, while terminal C-T forms a stable base pair. Hydrogen bonding between cytosine and thymine occurs only in the energy-minimized structures when the helix diameter decreases and the propeller twist angle between the bases increases. These changes are found to occur only at the end of a duplex in the calculations, which may explain the experimental results.     

Groove Width and Depth of B-DNA Structures Depend on Local Variation in Slide

Abstract

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

Sequence Effects On The Propeller Twist Of Base Pairs In DNA Helices

Abstract

Molecular mechanics calculations are performed on all the ten base pair steps (duplex dimers) and also a number of trimer and tetramer duplexes comprising them, in an attempt to systematically examine the possible base sequence effects on the magnitudes of propeller twists of base pairs at a given step. The analysis reveals that though propeller twist is a base pair property, it behaves very much like other base step parameters such as slide, roll, twist etc., Hence, it may be necessary to monitor the nature and variation of magnitudes of pt at a step. Calculations performed on 45 out of the 136 unique tetramer combinations involving all the ten unique base steps show that the difference in magnitudes of propeller twists of the base pairs of a given step has been found to be either steep or moderate depending on base pairs that flank the base step. These observations compare very well with the available experimental data. Tetramer sequences, wherein abase pair of a base step repeats in the same direction, exhibit a relatively steep difference in propeller twist at the step. Tetramers other than these exhibit moderate difference in propeller twist Such sequences are broadly classified as type-I and type-II respectively. Practically all the tetrads considered in the study, excepting those with GT step and a few involving CG and GC steps, conform to the above classification.     

A molecular dynamics study of the effect of G.T mispairs on the conformation of DNA in solution

The effect of G.T mispair incorporation into a double-helical environment was examined by molecular dynamics simulation. The 60-ps simulations performed on the two hexanucleotide duplexes d (G3C3)2 and d(G3TC2)2 included 10 Na+ counterions and first hydration shell waters. The resulting backbone torsional angle trajectories were analyzed to select time spans representative of conformational domains. The average backbone angles and helical parameters of the last time span for both duplexes are reported. During the simulation the hexamers retained B-type DNA structures that differed from typical A- or B-DNA forms. The overall helical structures for the two duplexes are vary similar. The presence of G.T mispairs did not alter the overall helical structure of the oligonucleotide duplex. Large propeller twist and buckle angles were obtained for both duplexes. The purine/pyrimidine crossover step showed a large decrease in propeller twist in the normal duplex but not in the mismatch duplex. Upon the formation of wobble mispairs in the mismatched duplex, the guanines moved into the minor groove and the thymines moved into the major groove. This helped prevent purine/purine clash and created a deformation in the relative orientation of the glycosidic bonds. It also exposed the free O4 of the thymines in the major groove and N2 of the guanines in the minor groove to interactions with solvent and counterions. These factors seemed to contribute to the apparently higher rigidity of the mismatched duplex during the simulation.     

Analysis of local helix geometry in three B-DNA decamers and eight dodecamers

Local variations in B-DNA helix structure are compared among three decamers and eight dodecamers, which contain examples of all ten base-pair step types. All pairwise combinations of helix parameters are compared by linear regression analysis, in a search for internal relationships as well as correlations with base sequence. The primary conclusions are: (1) Three-center hydrogen bonds between base-pairs occur frequently in the major groove at C-C, C-A, A-A and A-C steps, but are less convincing at C-C and C-T steps in the minor groove. The requirements for large base-pair propeller are (1) that the base-pair should be A.T rather than G.C, and (2) that it be involved in a major groove three-center hydrogen bond with the following base-pair. Either condition alone is insufficient. Hence, a large propeller is expected at the leading base-pair of A-A and A-C steps, but not at A-T, T-A, C-A or C-C steps. (2) A systematic and quantitative linkage exists between helix variables twist, rise, cup and roll, of such strength that the rise between base-pairs can hardly be described as an independent variable at all. Two typical patterns of behavior are observed at steps from one base-pair to the next: high twist profile (HTP), characterized by high twist, low rise, positive cup and negative roll, and low twist profile (LTP), marked by low twist, high rise; negative cup and positive roll. Examples of HTP are steps G-C, G-A and Y-C-A-R, where Y is pyrimidine and R is purine. Examples of LTP steps are C-G, G-G, A-G and C-A steps other than Y-C-A-R. (3) The minor groove is especially narrow across the two base-pairs of the following steps: A-T, T-A, A-A and G-A. (4) In general, base step geometry cannot be correlated solely with the bases that define the step in question; the two flanking steps also must be taken into account. Hence, local helix structure must be studied in the context, not of two base-pairs: A-B, but of four: x-A-B-y. Calladine's rules, although too simple in detail, were correct in defining the length of sequence over which a given perturbation is expressed. Whereas ten different two-base steps are possible, allowing for the identity of complementary sequences, there are 136 different four-base steps. Only 33 of these 136 four-base steps are represented in the decamer and dodecamer structures solved to date, and hence it is premature to try to set up detailed structural algorithms. (5) The sugar-phosphate backbone chains of B-DNA place strong limits on sequence-induced structural variation, damping down most variables within four or five base-pairs, and preventing purine-purine anti-anti mismatches from causing bulges in the double helix. Hence, although short-range sequence-induced deformations (or deformability) are observed, long-range deformations propagated down the helix are not to be expected.     

Determining local conformational variations in DNA. Nuclear magnetic resonance structures of the DNA duplexes d(CGCCTAATCG) and d(CGTCACGCGC) generated using back-calculation of the nuclear Overhauser effect spectra, a distance geometry algorithm and constrained molecular dynamics

Two-dimensional nuclear magnetic resonance (n.m.r.) spectroscopy and a variety of computational techniques have been used to generate three-dimensional structures of the two DNA duplexes d(CGCCTAATCG) and d(CGTCACGCGC). The central six base-pairs in these two decamers contain all ten dinucleotide pairs in DNA and thus, represent a model system for investigating how the local structure of DNA varies with base sequence. Resonance assignments were made for the non-exchangeable base protons and most of the C-1'-C-4' sugar protons in both decamers. Three-dimensional structures were generated using a distance geometry algorithm and these initial structures were refined by optimizing the fit of back-calculated spectra against the experimental two-dimensional nuclear Overhauser effect (NOE) spectra. This back-calculation procedure consists of calculating NOE cross relaxation rates for a given structure by solution of the Bloch equations, and directly accounts for spin diffusion effects. Use of this refinement procedure eliminates some assumptions that have been invoked when generating structures of DNA oligomers from n.m.r. data. Constrained energy minimization and constrained quenched molecular dynamics calculation were also performed on both decamers to help generate energetically favorable structures consistent with the experimental data. Analysis of the local conformational parameters of helical twist, helical rise, propeller twist, displacement and the alpha, beta, gamma, epison and zeta backbone torsion angles in these structures shows that these parameters span a large range of values relative to the X-ray data of nucleic acids. However, the glycosidic and pseudorotation angles are quite well defined in these structures. The implications that these results have for determination of local structural variations of DNA in solution, such as those predicted by Callidine's rules, are discussed. Our results differ significantly from some previous studies on determining local conformations of nucleic acids and comparisons with these studies are made.