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. © 1993 John Wiley & Sons, Inc.     

Conformational Flexibility of the Furanose Ring in DNA and its Dipole Moment

Abstract

The influence of conformational rearrangement of the furanose ring in DNA on its dipole moment is studied. The dipole moment of the deoxyribose molecule as a function of its puckered state is calculated by the quantum-mechanical method using the MINDO/3 approximation. The values of the dipole moment and its components are obtained at various magnitudes of the pseudorotation phase angle. The C3′-endo = C2′-endo conformational transition of deoxyribose is shown to be accompanied by the change in the dipole moment up to 3D. The results obtained are used to explain the structural properties of the DNA hydration shell.     

Conformational flexibility of the furanose ring in DNA and its dipole moment

The influence of conformational rearrangement of the furanose ring in DNA on its dipole moment is studied. The dipole moment of the deoxyribose molecule as a function of its puckered state is calculated by the quantum-mechanical method using the MINDO/3 approximation. The values of the dipole moment and its components are obtained at various magnitudes of the pseudorotation phase angle. The C3'-endo in equilibrium C2'-endo conformational transition of deoxyribose is shown to be accompanied by the change in the dipole moment up to 3D. The results obtained are used to explain the structural properties of the DNA hydration shell.     

Study of furanose ring flexibility in polynucleotide chains using Raman spectra analysis

The bands within the range of 800–850 cm^?1 of Raman spectra of polynucleotides sensitive to the change in conformation of sugar–phosphate backbone are analyzed theoretically. The bands are interpreted as the appearance of a quasi-local deoxyribose vibrational mode whose frequency is dependent on the ring puckering. The localization region of the vibrational mode is pointed out. The theory establishes a relationship between the observed spectral intensity and the population of deoxyribose conformational states described in the framework of the pseudorotation concept. The approach developed allowed one to describe the band shapes and their temperature behavior, and to determine the pseudorotation potential of deoxyribose in the helix B -form of A · T containing polynucleotides. Using the analysis of Raman spectra of DNA fibers in water–ethanol mixture the deoxyribose flexibility during the B -A transition is investigated in terms of the population of conformers and effective potential. It is shown that N- and S-type deoxyribose conformers are populated in the DNA B -form (those of the S-type are preferable), whereas N-type conformers are primarily populated in the DNA A -form.     

Structure of DNA hexamer sequence d-CGATCG by two-dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics

Solution conformation of self-complementary DNA duplex d-CGATCG, containing 5' d-CpG 3' site for intercalation of anticancer drug, daunomycin and adriamycin, has been investigated by nuclear magnetic resonance (NMR) spectroscopy. Complete resonance assignments of all the protons (except some H5'/H5" protons) have been obtained following standard procedures based on double quantum filtered correlation spectroscopy (dQF COSY) and two-dimensional nuclear Overhauser effect (NOE) spectra. Analysis of sums of coupling constants in one-dimensional NMR spectra, cross peak patterns in dQF COSY spectra and inter proton distances shows that the DNA sequence assumes a conformation close to the B-DNA family. The deoxyribose sugar conformation is in dynamic equilibrium with predominantly S-type conformer and a minor N-type conformer with N<-->S equilibrium varying with temperature. At 325 K, the mole fraction of the N-conformer increases for some of the residues by approximately 9%. Using a total of 10 spin-spin coupling constants and 112 NOE intensities, structural refinement has been carried out using Restrained Molecular Dynamics (rMD) with different starting structures, potential functions and rMD protocols. It is observed that pseudorotation phase angle of deoxyribose sugar for A3 and T4 residues is approximately 180 degrees and approximately 120 degrees, respectively while all other residues are close to C2'endo-conformation. A large propeller twist (approximately -18 degrees) and smallest twist angle (approximately 31 degrees) at A3pT4 step, in the middle of the sequence, a wider (12 A) and shallower (3.0 A) major groove with glycosidic bond rotation as high anti at both the ends of hexanucleotide are observed. The structure shows base-sequence dependent variations and hence strong local structural heterogeneity, which may have implications in ligand binding.     

DNA minor groove electrostatic potential: influence of sequence-specific transitions of the torsion angle gamma and deoxyribose conformations

The structural adjustments of the sugar-phosphate DNA backbone (switching of the γ angle (O5′–C5′–C4′–C3′) from canonical to alternative conformations and/or C2′-endo → C3′-endo transition of deoxyribose) lead to the sequence-specific changes in accessible surface area of both polar and non-polar atoms of the grooves and the polar/hydrophobic profile of the latter ones. The distribution of the minor groove electrostatic potential is likely to be changing as a result of such conformational rearrangements in sugar-phosphate DNA backbone. Our analysis of the crystal structures of the short free DNA fragments and calculation of their electrostatic potentials allowed us to determine: (1) the number of classical and alternative γ angle conformations in the free B-DNA; (2) changes in the minor groove electrostatic potential, depending on the conformation of the sugar-phosphate DNA backbone; (3) the effect of the DNA sequence on the minor groove electrostatic potential. We have demonstrated that the structural adjustments of the DNA double helix (the conformations of the sugar-phosphate backbone and the minor groove dimensions) induce changes in the distribution of the minor groove electrostatic potential and are sequence-specific. Therefore, these features of the minor groove sizes and distribution of minor groove electrostatic potential can be used as a signal for recognition of the target DNA sequence by protein in the implementation of the indirect readout mechanism.     

S-N transition of the saccharide ring in B-form DNA

A conformational transition of a single deoxyribose was analyzed in B-form trimers dA3:dT3 and dG3:dC3, both in the purine and pyrimidine chains. The main results were obtained for the duplexes with frozen ends, which could be extended by regular double helixes. The geometry of the central sugar ring in the duplexes may strongly deviate from the regular conformation. When deoxyribose changed its conformation in the central pyrimidine, the energy increase was proved to be less significant in comparison with that for purine. In the case of Thy, a decrease in pseudorotation angle P from 140 to 80 degrees causes the energy increase of 0.5 kcal/mol only, the barrier being 1.2 kcal/mol. The energy profile for Cyt has several local minima. The results of calculations were compared with numerous experimental data, they help to explain some NMR data. A perturbation of the duplex AAA:TTT structure caused by the thymine sugar ring transition, produces 5 degrees bend of the DNA axis directed toward adenines. We also investigated the influence of such conformational disturbance on the neighbouring base pairs, in particular the transition in the trimers with unfrozen ends.     

8 Using pseudorotation as a reaction coordinate in free energy simulations of nucleic acids

Backbone sugar groups are central components of nucleic acids. The conformations of the ribose/deoxyribose can be elegantly described using the concept of pseudorotation (Altona and Sundaralingam, 1972), and are dominated by the C2′- and C3′-endo conformers. The free energy barrier of the transition between these two major puckering modes can be probed by NMR relaxation experiments (Johnson and Hoogstraten, 2008), but an atomic picture of the transition path per se is only available for several truncated nucleoside analogues (Brameld & Goddard III, 1999). Here, we implemented a new free energy simulation method for Molecular Dynamics simulations using pseudorotation as the reaction coordinate (Cremer and Pople, 1975). This allowed us to compute the free energy landscape of a complete pseudorotation cycle. The free energy landscape revealed not only the relative stability of C2′- and C3′-endo conformers, but also the main transition path and its free energy barrier. As a validation of our new approach, we calculated free energy surface of the pseudorotation of guanosine monophosphate. The free energy surface revealed that the C2′-endo conformation is ?1?kcal/mol that is more stable and the free energy barrier for the transition is 4.5–5?kcal/mol. These are in excellent agreement with previous NMR measurements (Zhang et al., 2012; Röder et al., 1975). We have further applied this method to other systems that are important in pre-biotic chemistry, including an RNA duplex with unique 2′, 5′-phosphodiester linkages.     

Accurate representation of B-DNA double helical structure with implicit solvent and counterions

High-resolution nuclear magnetic resonance (NMR) and crystallographic data have been taken to refine the force field used in the torsion angle space nucleic acids molecular mechanics program DUPLEX. The population balance deduced from NMR studies of two carcinogen-modified DNA conformers in equilibrium was used to fine tune a sigmoidal, distance-dependent dielectric function so that reasonable relative energies could be obtained. In addition, the base-pair and backbone geometry from high-resolution crystal structures of the Dickerson-Drew dodecamer was used to re-evaluate the deoxyribose pseudorotation profile and the Lennard-Jones nonbonded energy terms. With a modified dielectric function that assumes a very steep distance-dependent form, a deoxyribose pseudorotation profile with reduced energy barriers between C2'- and C3'-endo minima, and a shift of the Lennard-Jones potential energy minimum to a distance approximately 0.4 A greater than the sum of the van der Waals' radii, the sequence-dependent conformational features of the Dickerson-Drew dodecamer in both the solid state and the aqueous liquid crystalline phase are well reproduced. The robust performance of the revised force field, in conjunction with its efficiency through implicit treatment of solvent and counterions, provides a valuable tool for elucidating conformations and structure-function relationships of DNA, including those of molecules modified by carcinogens and other ligands.     

Conformational studies of d-(AAAAATTTTT)2 using constraints from nuclear overhauser effects and from quantitative analysis of the cross-peak fine structures in two-dimensional 1H nuclear magnetic resonance spectra

The conformation at the dA-dT junction in d-(AAAAATTTTT)2 was investigated by using a variety of phase-sensitive two-dimensional nuclear magnetic resonance experiments at 500 MHz for detailed studies of the deoxyribose ring puckers. Conformational constraints were collected from two-dimensional nuclear Overhauser enhancement spectra recorded with short mixing times and from quantitative simulations of the cross-peaks in two-dimensional correlated spectra. Overall, the decamer duplex adopts a conformation of the B-DNA type, and for dA4 and dA5 the pseudorotation phase angle P is in the standard range 150-180 degrees. The deoxyribose puckers for the other nucleotides deviate significantly from the standard B-DNA structure. Spectrum simulations assuming either static deviations from standard B-DNA or a simple two-state dynamic equilibrium between the C2'-endo and C3'-endo forms of the deoxyribose were used to analyze the experimental data. It was thus found that the ring pucker for dT6 deviates from the regular C2'-endo form of B-DNA by a static distortion, with the pseudorotation phase angle P in the range 100-130 degrees, and a similar value of P is indicated for dT7. For the peripheral base pairs dynamic distortions of the C2'-endo form of the deoxyribose were found. In agreement with recent papers on related duplexes containing (dA)n tracts, we observed prominent nuclear Overhauser effects between adenine-2H and deoxyribose-1'H, which could be largely due to pronounced propeller twisting as observed in the crystal structures of (dA)n-containing compounds.     

Intrinsic conformational properties of deoxyribonucleosides: implicated role for cytosine in the equilibrium among the A, B, and Z forms of DNA.

Structural properties of biomolecules are dictated by their intrinsic conformational energetics in combination with environmental contributions. Calculations using high-level ab initio methods on the deoxyribonucleosides have been performed to investigate the influence of base on the intrinsic conformational energetics of nucleosides. Energy minima in the north and south ranges of the deoxyribose pseudorotation surfaces have been located, allowing characterization of the influence of base on the structures and energy differences between those minima. With all bases, chi values associated with the south energy minimum are lower than in canonical B-DNA, while chi values associated with the north energy minimum are close to those in canonical A-DNA. In deoxycytidine, chi adopts an A-DNA conformation in both the north and south energy minima. Energy differences between the A and B conformations of the nucleosides are <0.5 kcal/mol in the present calculations, except with deoxycytidine, where the A form is favored by 2.3 kcal/mol, leading the intrinsic conformational energetics of GC basepairs to favor the A form of DNA by 1.5 kcal/mol as compared with AT pairs. This indicates that the intrinsic conformational properties of cytosine at the nucleoside level contribute to the A form of DNA containing predominately GC-rich sequences. In the context of a B versus Z DNA equilibrium, deoxycytidine favors the Z form over the B form by 1.6 kcal/mol as compared with deoxythymidine, suggesting that the intrinsic conformational properties of cytosine also contribute to GC-rich sequences occurring in Z DNA with a higher frequency than AT-rich sequences. Results show that the east pseudorotation energy barrier involves a decrease in the furanose amplitude and is systematically lower than the inversion barrier, with the energy differences influenced by the base. Energy barriers going from the south (B form) sugar pucker to the east pseudorotation barrier are lower in pyrimidines as compared with purines, indicating that the intrinsic conformational properties associated with base may also influence the sugar pseudorotational population distribution seen in DNA crystal structures and the kinetics of B to A transitions. The present work provides evidence that base composition, in addition to base sequence, can influence DNA conformation.     

cis-Platinum induced distortions in DNA. Conformational analysis of d(GpCpG) and cis-pt(NH3)2[d(GpCpG)], studied by 500-MHz NMR

Proton NMR studies at 500 MHz in aqueous solution were carried out on the G-G chelated deoxytrinucleosidediphosphate platinum complex cis-Pt(NH3)2[d(GpCpG], on the uncoordinated trinucleotide d(GpCpG) and on the constituent monomers cis-Pt(NH3)2[d(Gp)]2, cis-Pt(NH3)2[d(pG)]2, d(Gp), d(pCp) and d(pG). Complete NMR spectral assignments are given and chemical shifts and coupling constants are analysed to obtain an impression of the detailed structure of d(GpCpG) and the distortion of the structure due to chelation with [cis-Pt(NH3)2]2+. Platination of the guanosine monophosphates affects the sugar conformational equilibrium to favour the N conformation of the deoxyribose ring. This feature is also apparent in ribose mononucleotides and is possibly caused by an increased anomeric effect. In cis-Pt(NH3)2[d(pG)]2 the phase angle of pseudorotation of the S-type sugar ring is 20 degrees higher than in 'free' d(pG) which might be an indication for an ionic interaction between the positive platinum and the negatively charged phosphate. It appears that d(GpCpG) reverts from a predominantly random coil to a normal right-handed B-DNA-like single-helical structure at lower temperatures, whereas the conformational features of cis-Pt(NH3)2[d(GpCpG)] are largely temperature-independent. In the latter compound much conformational freedom along the backbone angles is seen. The cytosine protons and deoxyribose protons exhibit almost no shielding effect as should normally be exerted by the guanine bases in stacking positions. This is interpreted in terms of a 'turning away' of the cytosine residue from both chelating guanines. Conformational features of cis-Pt(NH3)2[d(GpCpG)[ are compared with the 'bulge-out' of the ribose-trinucleotide m6(2)ApUpm6(2)A.     

Effect of trehalose on the contributions to the dipole potential of lipid monolayers

The dipole potential and the area changes induced by trehalose on dimyristoyl phosphatidylcholine (DMPC), 1,2-di-O-tetradecyl-sn-glycero-3-phosphocholine (dietherPC), dimyristoyl phosphatidylethanolamine (DMPE), 1,2-di-O-tetradecyl-sn-glycero-3-phosphoethanolamine (dietherPE) monolayers have been studied at different temperatures. The insertion of trehalose into DMPC monolayers in the fluid and gel states requires of the presence of carbonyl groups. The area increase observed at 0.15M trehalose is congruent with the decrease in the dipole potential. However, in dietherPC, in which trehalose does not affect the area, a decrease in the dipole potential is also observed. This is interpreted as a result of the displacement of water from the phosphate groups exposed to the aqueous phase. In DMPE, trehalose also decreases the dipole potential without affecting the area of saturated monolayers and in dietherPE no effect on dipole potential and area was observed. It is concluded that the spacer effect of trehalose depends on the specific interaction with CO, which is modulated by the strength of the interaction of the PO groups with lateral NH groups. However, it is not the only contribution to the dipole potential decrease.     

The rate of Qx→Qy relaxation in bacteriochlorophylls of reaction centers from Rhodobacter sphaeroides determined by kinetics of the ultrafast carotenoid bandshift

Transient absorption changes induced by excitation of isolated reaction centers (RCs) from Rhodobacter sphaeroides with 600 nm laser pulses of 20 fs (full width at half maximum) were monitored in the wavelength region of 420–560 nm. The spectral features of the spectrum obtained are characteristic for an electrochromic band shift of the single carotenoid (Car) molecule spheroidene, which is an integral constituent of these RCs. This effect is assigned to an electrochromic bandshift of Car due to the local electric field of the dipole moment formed by electronic excitation of bacteriochlorophyll (BChl) molecule(s) in the neighborhood of Car. Based on the known distances between the pigments, the monomeric BChl (B_B) in the inactive B-branch is inferred to dominate this effect. The excitation of B_B at 600 nm leads to a transition into the S₂ state (Q_x band), which is followed by rapid internal conversion to the S₁ state (Q_y band), thus leading to a change of strength and orientation of the dipole moment, i.e., of the electric field acting on the Car molecule. Therefore, the time course of the electrochromic bandshift reflects the rate of the internal conversion from S₂ to S₁ of B_B. The evaluation of the kinetics leads to a value of 30 fs for this relaxation process. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.     

The Interaction of Dipole Modifiers with Polyene-Sterol Complexes

Recently, we showed that the effect of dipole modifiers (flavonoids and styrylpyridinium dyes) on the conductance of single amphotericin B (AmB) channels in sterol-containing lipid bilayers primarily resulted from changes in the membrane dipole potential. The present study examines the effect of dipole modifiers on the AmB multi-channel activity. The addition of phloretin to cholesterol-containing membranes leads to a significant increase in the steady-state AmB-induced transmembrane current. Quercetin significantly decreases and RH 421 increases the current through ergosterol-containing bilayers. Other tested flavonoids and styrylpyridinium dyes do not affect the channel-forming activity of AmB independently on the sterol composition of the bilayers. The effects obtained in these trials may instead be attributed to the direct interaction of dipole modifiers with AmB/sterol complexes and not to the effect of dipole potential changes. The presence of double bonds in the Δ7 and Δ22 positions of sterol molecules, the number of conjugated double bonds and amino sugar residues in polyene molecules, and the conformation and adsorption plane of dipole modifiers are important factors impacting this interaction.     

Stereochemical studies on nucleic acid analogues. I. Conformations of alpha-nucleosides and alpha-nucleotides: interconversion of sugar puckers via O4'-exo.

The preferred conformations of ribo and deoxyribo alpha-nucleosides and alpha-nucleotides, the stereoisomers of the naturally occurring beta-isomers, are worked out by minimizing the conformational energy as a function of all the major parameters including the sugar ring conformations along the pseudorotation path. The results of the studies bring out certain distinct conformational features that are at variance with their beta counterparts. The range of glycosyl conformations are found to be not only quite restricted here but favor predominantly the anti conformation. The syn glycosyl conformation for the entire region of P values are found to be energetically less favorable, with the barrier to anti in equilibrium with syn interconversion being higher especially in alpha-ribonucleosides. The energetically preferred range of pseudorotation phase angles (P) is also considerably restricted and P values corresponding to the C1'-exo range of sugars are highly unfavorable for alpha-nucleosides, in sharp contrast to the broad range of sugar ring conformations favored by beta-isomers. While both trans congruent to 180 degrees and skew congruent to 270 degrees conformations around the C3'-O3' (phi') bond are favored for alpha-3'-nucleotides with deoxyribose sugars, ribose sugars seem to favor only the skew values of phi'. Most interestingly and in sharp contrast to beta-stereoisomers, an energy barrier is encountered at P values corresponding to O4'-endo sugars. This suggests that the possible sugar pucker interconversion between C2'-endo/C3'-exo and C3'-endo/C2'-exo in alpha-anomers could take place only through the O4'-exo region. Likewise the possible path of anti in equilibrium with syn interconversion in alpha-nucleosides is not via high anti, in sharp contrast to beta-nucleosides. These observations should be borne in mind while assigning the sugar ring conformers in alpha-nucleosides and those containing them from nmr investigations. Comparison of the results with beta-anomers seem to suggest on the whole a lack of conformational variability or the restricted nature of alpha-stereoisomers. This could be one of the reasons for its nonselection in the naturally occurring nucleic acids.     

Three-state models of furanose pseudorotation

A general procedure is described to treat the pseudorotation of the furanose ring in terms of a three-state conformational equilibrium. In addition to the principal n (C3'-endo) and s (C2'-endo) puckering domains, the unusual e (01'-endo) intermediate is included in the analysis. Each of these three conformational categories is represented by a blend of five closely related puckered forms rather than by a single rotational isomeric state. Using this model together with experimentally measured nmr coupling constants, the puckering populations of various nucleic acid analogs are estimated. The conventional two-state n/s equilibria is confirmed in ordinary ribose and deoxyribose systems. The e domain, however, is found to be of major importance in several chemically modified furanoses including certain pyrimidine deoxynucleosides damaged by radiation and various nucleosides and nucleotides forced by bulky substituents on the base into unusual syn glycosyl arrangements. The "free" pseudorotation of these modified systems is not detected by conventional two-state puckering analyses.     

Oligonucleotide conformations. III. Comparative optical and thermodynamic studies of uridylyl-3'-5'-nucleosides containing ribose, deoxyribose, or 2'-deoxy-2'-fluororibose in the uridine moiety.

The synthesis of uridylyl-3′-5′-nucleosides containing ribose, deoxyribose, or 2′-fluoro-2′-deoxyribose in the uridine-3′-bound moiety and adenosine, guanosine, cytidine or uridine in the 5′-nucleoside is reported. The temperature dependence of the circular dichroism of these dinucleoside phosphates in 0.06 M phosphate buffer at pH 7 was analyzed by the two-state model and the oscillating dimer model. From the former, apparent thermodynamic parameters were determined by means of an iterative computer method. The comparison between the three different dinucleoside phosphates in each series indicated that the fluororiboside and the riboside resembled each other and were more stacked than the analogue containing deoxyribose. It further appeared that the similarity between the fluororiboside and the riboside is influenced by the nature of the neighboring 5′-bound base. The interaction between the 3′-bound sugar moiety and the 5′-bound base is evoked as a possible stabilization mechanism.     

NMR relaxation and the structure of a synthetic DNA poly(dA-dT).poly(dA-dT)

1H-1H and 31P-1H nuclear Overhauser effects and 31P NMR spin-lattice relaxation times were measured for a synthetic DNA poly(dA-dT).poly(dA-dT) in a low-salt aqueous solution. The results have shown that all bases in the double helix are anti-orientated with respect to deoxyribose residues and that the sugar-phosphate backbone has an alternating architecture.