DNA bending and sequence-dependent backbone conformation NMR and computer experiments

DNA polymorphism in crystals: three stable conformations for the decadeoxynucleotide d(GCATGCATGC)

High-resolution structures of DNA fragments determined using X-ray crystallography or NMR have provided descriptions of a veritable alphabet of conformations. They have also shown that DNA is a flexible molecule, with some sequences capable of adopting two different structures. Here, the first example is presented of a DNA fragment that can assume three different and distinct conformations in crystals. The decanucleotide d(GCATGCATGC) was previously reported to assume a single-stranded double-fold structure. In one of the two crystal structures described here the decamer assumes both the double-fold conformation and, simultaneously, the more conventional B-type double-helical structure. In the other crystal the sequence assumes the A-type double-helical conformation. These results, taken together with CD spectra, which were recorded as the decamer was titrated against four metal ions and spermine, indicate that the molecule may exist as a mixed population of structures in solution. Small differences in the environmental conditions, such as the concentration of metal ion, may decide which of these crystallizes out. The results also support the idea that it may be possible for DNA to change its structure to suit the binding requirements of proteins or drugs.

A non-Karplus effect: evidence from phosphorus heterocycles and DFT calculations of the dependence of vicinal phosphorus-hydrogen NMR coupling constants on lone-pair conformation

In contrast to literature reports of a Karplus-type curve that correlates (3)J(PH) with phosphorus-hydrogen dihedral angle, a recently reported glycine-derived 1,3,2-oxazaphospholidine (7c) has two hydrogen atoms on the ring with identical PNCH dihedral angles but measured coupling constants of ∼6 and 1.5 Hz. DFT calculations were in accord with these values and suggested that the smaller coupling constant is negative. Experimental evidence of the opposite signs of these coupling constants was obtained by analysis of the ABX NMR spectrum of the new glycine-derived N-p-toluenesulfonyl phosphorus heterocycle 6c. DFT calculations on 6c and on Me(2)NPCl(2) and t-BuPCl(2) were also in accord with NMR data and allowed confirmation of unusual features including a lone pair effect on (3)J(PH), the negative coupling constant, temperature-dependent chemical shifts due to rotation about the sulfonamide S-N bond, and vicinal phosphorus-hydrogen coupling constants over 40 Hz. Calculation of phosphorus-hydrogen coupling constants both as a function of PYCH dihedral angle θ (Y = O, N, C) and lone pair-PYC dihedral angle ω shows similar θ,ω surfaces for (3)J(PH) with a range of (3)J(PH) from -4.4 to +51 Hz and demonstrates the large non-Karplus effect of lone-pair conformation on vicinal phosphorus-hydrogen coupling constants.

Probing phosphorylation by non-mammalian isoprenoid biosynthetic enzymes using (1)H-(31)P-(31)P correlation NMR spectroscopy

The biogenesis of isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP) is accomplished by the methylerythritol phosphate (MEP) pathway in plants, bacteria and parasites, making it a potential target for the development of anti-infective agents and herbicides. The biosynthetic enzymes comprising this pathway catalyze intriguing chemical transformations on diphosphate scaffolds, offering an opportunity to generate novel analogs in this synthetically challenging compound class. Such a biosynthetic approach to generating new diphosphate analogs may involve transformation through discrete diphosphate species, presenting unique challenges in structure determination and characterization of unnatural enzyme-generated diphosphate products produced in tandem. We have developed (1)H-(31)P-(31)P correlation NMR spectroscopy techniques for the direct characterization of crude MEP pathway enzyme products at low concentrations (200 microM to 5 mM) on a room temperature (non-cryogenic) NMR probe. Coupling the 100% natural abundance of the (31)P nucleus with the high intrinsic sensitivity of proton NMR, (1)H-(31)P-(31)P correlation spectroscopy is particularly useful for characterization of unnatural diphosphate enzyme products in the MEP pathway. As proof of principle, we demonstrate the rapid characterization of natural enzyme products of the enzymes IspD, E and F in tandem enzyme incubations. In addition, we have characterized several unnatural enzyme products using this technique, including new products of cytidyltransferase IspD bearing erythritol, glycerol and ribose components. The results of this study indicate that IspD may be a useful biocatalyst and highlight (1)H-(31)P-(31)P correlation spectroscopy as a valuable tool for the characterization of other unnatural products in non-mammalian isoprenoid biosynthesis.

Local and global effects of strong DNA bending induced during molecular dynamics simulations

DNA bending plays an important role in many biological processes, but its molecular and energetic details as a function of base sequence remain to be fully understood. Using a recently developed restraint, we have studied the controlled bending of four different B-DNA oligomers using molecular dynamics simulations. Umbrella sampling with the AMBER program and the recent parmbsc0 force field yield free energy curves for bending. Bending 15-base pair oligomers by 90 degrees requires roughly 5 kcal mol(-1), while reaching 150 degrees requires of the order of 12 kcal mol(-1). Moderate bending occurs mainly through coupled base pair step rolls. Strong bending generally leads to local kinks. The kinks we observe all involve two consecutive base pair steps, with disruption of the central base pair (termed Type II kinks in earlier work). A detailed analysis of each oligomer shows that the free energy of bending only varies quadratically with the bending angle for moderate bending. Beyond this point, in agreement with recent experiments, the variation becomes linear. An harmonic analysis of each base step yields force constants that not only vary with sequence, but also with the degree of bending. Both these observations suggest that DNA is mechanically more complex than simple elastic rod models would imply.

DNA and its counterions: a molecular dynamics study

The behaviour of mobile counterions, Na+ and K+, was analysed around a B-DNA double helix with the sequence CCATGCGCTGAC in aqueous solution during two 50 ns long molecular dynamics trajectories. The movement of both monovalent ions remains diffusive in the presence of DNA. Ions sample the complete space available during the simulation time, although individual ions sample only about one-third of the simulation box. Ions preferentially sample electronegative sites around DNA, but direct binding to DNA bases remains a rather rare event, with highest site occupancy values of <13%. The location of direct binding sites depends greatly on the nature of the counterion. While Na+ binding in both grooves is strongly sequence-dependent with the preferred binding site in the minor groove, K+ mainly visits the major groove and binds close to the centre of the oligomer. The electrostatic potential of an average DNA structure therefore cannot account for the ability of a site to bind a given cation; other factors must also play a role. An extensive analysis of the influence of counterions on DNA conformation showed no evidence of minor groove narrowing upon ion binding. A significant difference between the conformations of the double helix in the different simulations can be attributed to extensive alpha/gamma transitions in the phosphate backbone during the simulation with Na+. These transitions, with lifetimes over tens of nanoseconds, however, appear to be correlated with ion binding to phosphates. The ion-specific conformational properties of DNA, hitherto largely overlooked, may play an important role in DNA recognition and binding.

Effect of phosphorothioate chirality on the grooves of DNA double helices: a molecular dynamics study

Phosphorothioate oligonucleotides (PS-ODNs) have gained considerable attention in drug therapy, primarily as potent antisense or antigene oligomers, which bind to specific DNA or mRNA sequences and lead to transcriptional or translational arrest. These are obtained by substituting one of the anionic oxygen of the phosphate group by a sulfur atom, which introduces chirality to the phosphorus atom of the DNA backbone. In this molecular dynamics simulation study, structural parameters like groove widths, environmental parameters like hydration or cation binding, and electrostatic energy surfaces of both the chiral forms of DNA/PS-DNA duplexes were assessed and compared with that of a normal DNA. Results indicate that, PS-S form with its sulfur atoms facing the minor groove has a widened minor groove, while the scenario is reverse for the PS-R form. Further analysis reveals the existence of several factors like large van der Waals radius of sulfur and the effect it has on its neighboring hydration pattern along with the net electrostatic environment, influencing such structural alterations. This also indicates, for the first time, the effect of absolute phosphorothioate chirality on the global structure of a DNA/PS-DNA hybrid that otherwise resembles a regular B-DNA structure.

Protein-DNA hydrophobic recognition in the minor groove is facilitated by sugar switching

Information readout in the DNA minor groove is accompanied by substantial DNA deformations, such as sugar switching between the two conformational domains, B-like C2'-endo and A-like C3'-endo. The effect of sugar puckering on the sequence-dependent protein-DNA interactions has not been studied systematically, however. Here, we analyzed the structural role of A-like nucleotides in 156 protein-DNA complexes solved by X-ray crystallography and NMR. To this end, a new algorithm was developed to distinguish interactions in the minor groove from those in the major groove, and to calculate the solvent-accessible surface areas in each groove separately. Based on this approach, we found a striking difference between the sets of amino acids interacting with B-like and A-like nucleotides in the minor groove. Polar amino acids mostly interact with B-nucleotides, while hydrophobic amino acids interact extensively with A-nucleotides (a hydrophobicity-structure correlation). This tendency is consistent with the larger exposure of hydrophobic surfaces in the case of A-like sugars. Overall, the A-like nucleotides aid in achieving protein-induced fit in two major ways. First, hydrophobic clusters formed by several consecutive A-like sugars interact cooperatively with the non-polar surfaces in proteins. Second, the sugar switching occurs in large kinks promoted by direct protein contact, predominantly at the pyrimidine-purine dimeric steps. The sequence preference for the B-to-A sugar repuckering, observed for pyrimidines, suggests that the described DNA deformations contribute to specificity of the protein-DNA recognition in the minor groove.

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 epsilon, zeta, alpha, beta and gamma 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 alpha/gamma 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 epsilon/zeta 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 conformational 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.

Alpha/gamma transitions in the B-DNA backbone

In the crystal structures of protein complexes with B-DNA, alpha and gamma DNA backbone torsion angles often exhibit non-canonical values. It is not known if these alternative backbone conformations are easily accessible in solution and can contribute to the specific recognition of DNA by proteins. We have analysed the coupled transition of the alpha and gamma torsion angles within the central GpC step of a B-DNA dodecamer by computer simulations. Five stable or metastable non-canonical alpha/gamma sub-states are found. The most favourable pathway from the canonical alpha/gamma structure to any unusual form involves a counter-rotation of alpha and gamma, via the trans conformation. However, the corresponding free energy indicates that spontaneous flipping of the torsions is improbable in free B-DNA. This is supported by an analysis of the available high resolution crystallographic structures showing that unusual alpha/gamma states are only encountered in B-DNA complexed to proteins. An analysis of the structural consequences of alpha/gamma transitions shows that the non-canonical backbone geometry influences essentially the roll and twist values and reduces the equilibrium dispersion of structural parameters. Our results support the hypothesis that unusual alpha/gamma backbones arise during protein-DNA complexation, assisting the fine structural adjustments between the two partners and playing a role in the overall complexation free energy.

DNA, RNA and hybrid RNA-DNA oligomers of identical sequence: structural and dynamic differences

A 27-mer sequence was synthesised as DNA duplex (DD), RNA duplex (RR), and RNA-DNA (RD) hybrid in order to characterise their structural and dynamic features. The hydrodynamic radius (Rh) and the rise (b) values of the three samples were consistent with the conformations predicted by CD analysis. The value of the torsional constant (alpha) of the samples containing RNA was approximately twice that of the DD sample and followed the order: DD < RD < RR. The same order was observed in the thermodynamic stability and in the reduction of the electrophoretic mobility. gamma-Ray footprinting analysis was carried out to resolve the individual strand conformation in the hybrid. The RNA strand preserved its conformation, while the DNA strand showed local deformations mainly at TA and TG steps.