Conformational analysis of the single-helical DNA fragment d(T-A-A-T) in aqueous solution. The combined use of NMR proton chemical shifts and coupling constants obtained at 300 MHz and 500 MHz

Development of Force Field Parameters for the Simulation of Single- and Double-Stranded DNA Molecules and DNA-Protein Complexes

Although molecular dynamics (MD) simulations have been used extensively to study the structural dynamics of proteins, the role of MD simulation in studies of nucleic acid based systems has been more limited. One contributing factor to this disparity is the historically lower level of accuracy of the physical models used in such simulations to describe interactions involving nucleic acids. By modifying nonbonded and torsion parameters of a force field from the Amber family of models, we recently developed force field parameters for RNA that achieve a level of accuracy comparable to that of state-of-the-art protein force fields. Here we report force field parameters for DNA, which we developed by transferring nonbonded parameters from our recently reported RNA force field and making subsequent adjustments to torsion parameters. We have also modified the backbone charges in both the RNA and DNA parameter sets to make the treatment of electrostatics compatible with our recently developed variant of the Amber protein and ion force field. We name the force field resulting from the union of these three parameter sets (the new DNA parameters, the revised RNA parameters, and the existing protein and ion parameters) DES-Amber. Extensive testing of DES-Amber indicates that it can describe the thermal stability and conformational flexibility of single- and double-stranded DNA systems with a level of accuracy comparable to or, especially for disordered systems, exceeding that of state-of-the-art nucleic acid force fields. Finally, we show that, in certain favorable cases, DES-Amber can be used for long-timescale simulations of protein–nucleic acid complexes.

Oligonucleotides conjugated to natural lipids: synthesis of phosphatidyl-anchored antisense oligonucleotides

Antisense oligonucleotides are promising therapeutic agents against a variety of diseases. Effective delivery of these molecules is critical in view of their clinical application. Despite the richness of synthetic strategies addressed to the lipophilic modification of oligodeoxynucleotides (ODNs) for enhancing their pharmacokinetic behavior and trans-membrane delivery, the phosphatidyl group (1,2-di-O-acyl-sn-glycero-3-phosphoryl) has been never used as the lipophilic moiety of lipid-ODN conjugates. The present paper reports a general procedure for synthesizing 5'-phosphatidyl-ODNs. By this procedure, phosphatidyl conjugates of a VEGF antisense-ODN have been prepared, which differ in the fatty acid composition of their phosphatidyl moiety. These new lipid-ODN conjugates, which have been characterized on the basis of their physicochemical properties, showed an improved resistance to exonucleases and were able to lower the VEGF-mRNA expression in human SH-SY5Y neuroblastoma cells more effectively than the relevant free antisense-ODN did.

Impact of C5-cytosine methylation on the solution structure of d(GAAAACGTTTTC)2. An NMR and molecular modelling investigation

The solution structures of d(GAAAACGTTTTC)2 and of its methylated derivative d(GAAAAMe5CGTTTTC)2 have been determined by NMR and molecular modelling in order to examine the impact of cytosine methylation on the central CpG conformation. Detailed 1H NMR and 31P NMR investigation of the two oligomers includes quantitative NOESY, 2D homonuclear Hartmann-Hahn spectroscopy, double-quantum-filtered COSY and heteronuclear 1H-31P correlation. Back-calculations of NOESY spectra and simulations of double-quantum-filtered COSY patterns were performed to gain accurate information on interproton distances and sugar phase angles. Molecular models under experimental constraints were generated by energy minimization by means of the molecular mechanics program JUMNA. The MORASS software was used to iteratively refine the structures obtained. After methylation, the oligomer still has a B-DNA conformation. However, there are differences in the structural parameters and the thermal stability as compared to the unmethylated molecule. Careful structural analysis shows that after methylation CpG departs from the usual conformation observed in other ACGT tetramers with different surroundings. Subtle displacements of bases, sugars and backbone imposed by the steric interaction of the two methyl groups inside the major groove are accompanied by severe pinching of the minor groove at the C-G residues.

Thermodynamics of melting of the circular dumbbell d<pCGC-TT-GCG-TT>

The conformational behavior of DNA minihairpin loops is sensitive to the directionality of the base pair that closes the loop. Especially tailored circular dumbbells, consisting of a stem of three Watson-Crick base pairs capped on each side with a minihairpin loop, serve as excellent model compounds by means of which deeper insight is gained into the relative stability and melting properties of hairpin loops that differ only in directionality of the closing pair: C-G vs G-C. For this reason the thermodynamic properties of the circular DNA decamers 5'-d<pCGC-TT-GCG-TT>-3' (I) and reference compounds 5'-d<pGGC-TT-GCC-TT>-3' (II) and 5'-d(GCG-TC-CGC)-3' (III) are studied by means of nmr spectroscopy. Molecules I and II adopt dumbbell structures closed on both sides by a two-membered hairpin loop. At low temperature I consists of a mixture of two slowly exchanging forms, denoted L2L2 and L2L4. The low-temperature L2L2 form is the fully intact minihairpin structure with three Watson-Crick C-G base pairs. The high-temperature form, L2L4, contains a partially disrupted closing G-C base pair in the 5'-GTTC-3' loop, with the cytosine base placed in a syn orientation. The opposite 5'-CTTG-3' loop remains stable. A study of the noncircular hairpin structure III shows similar conformational behavior for the 5'-GTTC-3' loop as found in I; a syn orientation for C(6) and two slowly exchanging imino proton signals for G(3). The melting point Tm of II was estimated to lie above 365 K. The Tm value of the duplex stem and the 5'-CTTG-3' loop of the L2L4 form of I is 352 +/- 2 K. The delta H0 is calculated as -89 +/- 10 kJ/mol. The Tm value determined for the individual residues of the 5'-GTTC-3' loop lies 4 degrees-11 degrees lower. The enthalpy delta H0 of melting the thymine residues in the 5'-GTTC-3' loop is calculated to be -61 +/- 7 kJ/mol. Thermodynamic data of the equilibrium between the slowly exchanging two- and four-membered loop conformers of I reveal an upper limit for delta H0 of +30 kJ/mol in going from a two-membered to a four-membered loop, in agreement with the enthalpy difference of +28 kJ/mol between the two loops at the Tm midpoint. For hairpin III the upper limit for delta H0 in going from a two-membered to a four-membered loop amounts to +/- 21 kJ/mol.(ABSTRACT TRUNCATED AT 400 WORDS)

A Monte Carlo method for generating structures of short single-stranded DNA sequences

A Monte Carlo method has been developed for generating the conformations of short single-stranded DNAs from arbitrary starting states. The chain conformers are constructed from energetically favorable arrangements of the constituent mononucleotides. Minimum energy states of individual dinucleotide monophosphate molecules are identified using a torsion angle minimizer. The glycosyl and acyclic backbone torsions of the dimers are allowed to vary, while the sugar rings are held fixed in one of the two preferred puckered forms. A total of 108 conformationally distinct states per dimer are considered in this first stage of minimization. The torsion angles within 5 kcal/mole of the global minimum in the resulting optimized states are then allowed to vary by +/- 10 degrees in an effort to estimate the breadth of the different local minima. The energies of a total of 2187 (3(7)) angle combinations are examined per local conformational minimum. Finally, the energies of all dinucleotide conformers are scaled so that the populations of differently puckered sugar rings in the theoretical sample match those found in nmr solution studies. This last step is necessitated by limitations in the theoretical methods to predict DNA sugar puckering accurately. The conformer populations of the individual acyclic torsion angles in the composite dimer ensembles are found to be in good agreement with the distributions of backbone conformations deduced from nmr coupling constants and the frequencies of glycosyl conformations in x-ray crystal structures, suggesting that the low energy states are reasonable. The low energy dimer forms (consisting of 150-325 conformational states per dimer step) are next used as variables in a Monte Carlo algorithm, which generates the conformations of single-stranded d(CXnG) chains, where X = A, T and n = 3, 4, 5. The oligonucleotides are built sequentially from the 5' end of the chain using random numbers to select the conformations of overlapping dimer units. The simulations are very fast, involving a total of 10(6) conformations per chain sequence. The potential errors in the buildup procedure are minimized by taking advantage of known rotational interdependences in the sugar-phosphate backbone. The distributions of oligonucleotide conformations are examined in terms of the magnitudes, positions, and orientations of the end-to-end vectors of the chains. The differences in overall flexibility and extension of the oligomers are discussed in terms of the conformations of the constituent dinucleotide steps, while the general methodology is discussed and compared with other nucleic acid model building techniques.

Conformational consequences of the incorporation of arabinofuranosylcytidine in DNA. An NMR study of the DNA fragments d(CGCTAGCG) and d(CGaCTAGCG) in solution

The self-complementary octamers d(CGCTAGCG) and d(CGaCTAGCG) (aC, arabinofuranosylcytidine) were studied by means of NMR spectroscopy. It is shown that d(CGaCTAGCG), under suitable conditions of oligonucleotide concentration, ionic strength and temperature, exclusively adopts a hairpin structure. However, under the same experimental conditions (5 mM DNA, no added salt, 295 K) d(CGCTAGCG) mainly adopts a B-DNA-type duplex. At lower temperatures (less than or equal to 290 K) the hairpin form of d(CGaCTAGCG) occurs in slow exchange with an intact B-DNA-type duplex. When the DNA concentration of d(CGCTAGCG) is dramatically reduced (less than or equal to 0.5 mM) the hairpin form becomes highly favoured at the expense of the dimer. Moreover, proton-chemical-shift considerations indicate that the structural features of the hairpin structure of d(CGCTAGCG) mimic, in part, those of the modified octamer d(CGaCTAGCG), i.e. a loop comprising only the two central residues with the thymine located into the minor groove (Pieters, J. M. L., de Vroom, E., van der Marel, G. A., van Boom, J. H., Koning, T. M. G., Kaptein, R. and Altona, C. unpublished results). Thermodynamic analysis of d(CGCTAGCG) yields an average Tmd value of 342 K (1 M DNA) and a delta Hod value of -266 kJ/mol for the dimer/coil transition and an average Tmh value of 321 K and delta Hoh - 102 kJ/mol for the hairpin/coil equilibrium. For the duplex/coil equilibrium of d(CGaCTAGCG) an average Tmd value of 336 K (1 M DNA) and delta Hod value of -253 kJ/mol are deduced. The hairpin/coil transition of d(CGaCTAGCG) is characterized by a delta Hoh value of -104 kJ/mol and an average Tmh value of 331 K. It is concluded that incorporation of an arabinofuranosylcytidine in the octamer d(CGaCTAGCG) results in stabilization of the hairpin form, whereas the dimer is destablized by two aC.dG base pairs.

Conformational and thermodynamic consequences of the introduction of a nick in duplexed DNA fragments: an NMR study augmented by biochemical experiments

NMR studies were carried out on various equimolar mixtures consisting of a combination of oligomers: d(ACGGCT) (I). d(pACGGCT) (Ia), d(TGCAGT) (II), d(AGCCGTACTGCA) (III), d(TGCAGTACGGCT) (IV). It is shown that I + II + III (MI) and Ia + II + III (M2) form stable duplexes with nicks in the centre of the respective double helices. A close analysis of the NOESY experiments of M1 and M2 revealed that these fragments form B-DNA type duplex structures. A comparison of the chemical-shift data of the nicked duplexes with those of the intact duplex of III + IV (M3) demonstrated that only small local distortions occur when a nick is introduced. The chemical-shift profiles of M1 and M3 were used to obtain the thermodynamic data for the duplex/coil transitions. The profiles of M1 were analysed by means of a new thermodynamic model (TRIDUP). From the calculated thermodynamic data of M1 and M3 it is concluded that the melting behaviour of M1 occurs cooperatively. A ligation experiment demonstrated that the relatively small substrate (M2) was almost completely joined after an overnight incubation at 14 degrees C.

Thermodynamics of the various forms of the dodecamer d(ATTACCGGTAAT) and of its constituent hexamers from proton nmr chemical shifts and UV melting curves: three-state and four-state thermodynamic models

Chemical shifts of base and H1' protons of the single-stranded hexamers d(ATTACC) and d(GGTAAT), of the 1:1 mixtures of these complementary hexamers, and of the self-complementary dodecamer d(ATTACCGGTAAT) were measured at various temperatures in aqueous solution. Four different sample concentrations were used in the case of the dodecamer and of the mixture of the complementary hexamers; the individual hexamers were measured at two different DNA concentrations. Absorbance temperature profiles at five different NaCl concentrations were measured for the dodecamer in order to quantify the effect of the ionic strength on the duplex formation. Under suitable conditions of nucleotide concentration, temperature, and ionic strength, the dodecamer adopts either a B-DNA duplex or a hairpin-loop structure. Chemical shift vs temperature profiles, constructed for all samples, were used to obtain thermodynamic parameters either for the various stacking interactions in the single strands or for the duplex or the hairpin-loop formation. In the analysis of the duplex formation of the hexamers, a two-state approach appeared too simple, because systematic deviations were revealed. Therefore, a new three-state model (DUPSTAK) was developed. In order to investigate the magnitude of error arising from the use of the two-state approach in cases where the DUPSTAK model appears more appropriate, a series of test calculations was made. The magnitude of error in the enthalpy and in the entropy of duplex melting is found to depend linearly upon the actual melting temperature and not upon the individual delta Hd degrees and delta Sd degrees values. Thermodynamic analysis of the chemical shift vs temperature profiles in D2O solution (no added salt) yields an average Tmd value of 341 K (1M DNA) and delta Hd degrees of - 121 kJ.mol-1 for the dimer/random-coil transition of the hexamer duplex d(ATTACC).d(GGTAAT). For the duplex in equilibrium random-coil transition of the 12-mer d(ATTACCGGTAAT) an average Tmd value of 336 K (1M DNA) and delta Hd degrees of -372 kJ.mol-1 are found. The hairpin/random-coil transition of d(ATTACCGGTAAT) is characterized by a rather large delta Hh degrees value, -130 kJ.mol-1, and an average Tmh value of 304 K.

Influence of the base sequence on the conformational behaviour of DNA polynucleotides in solution

NMR studies were carried out on samples of the non-self-complementary tetramers d(C-A-C-A), d(T-G-T-G), d(G-A-G-A) and d(T-C-T-C) and of 1:1 mixtures of the complementary tetramers d(C-A-C-A) X d(T-G-T-G) and d(G-A-G-A) X d(T-C-T-C) at two DNA concentrations and of the self-complementary octamers d(C-A-C-A-T-G-T-G) and d(G-A-G-A-T-C-T-C). Assignments, based upon one-dimensional NOE and homonuclear-decoupling and two-dimensional correlated and NOE spectroscopies are given of the resonances of most of the base and sugar protons. Chemical shift vs temperature profiles, constructed for all samples, yielded insight into the temperature- and concentration-dependent conformational behaviour of the compounds and were used to obtain thermodynamic parameters pertaining to the stacked-single-strand----random-coil and duplex----random-coil equilibria. Vicinal proton-proton couplings were analyzed in terms of the conformation of the deoxyribose rings in the single-stranded tetramers and duplexed octamers. The NOE patterns, chemical shift profiles, imino-proton resonances and coupling data revealed that the compounds adopt B-DNA-like structures. The ratio duplexed/stacked-single-strand/random coil depends upon external conditions as well as upon base sequence. The thermodynamic data indicate that: in terms of single-helical stacking, the R-R steps (Tm 321-328 K) appear more stable than the Y-R or R-Y steps (Tm 308-316 K) and the Y-Y steps score least (Tm 290-300 K), and the duplexes consisting of alternating, d(Y-R)n, strands are more stable, in terms of delta H degrees, compared to the d(R-R)n X d(Y-Y)n duplexes. The analyses of the couplings demonstrated that the sugars of the single-stranded tetramers and duplexed octamers occur as a blend of N- and S-type conformers, with a preference for the S-type (C2'-endo) sugar conformation: upon duplex formation, no significant shift in the N-type/S-type ratio was observed. The fraction S-type sugar conformation of a given residue, %S, in the stacked-single strands was found to depend upon the nature of its own base and that of the adjacent residues: sugars in an R-R stretch display high values of %S (90-100), whereas those in Y-Y stretches show relatively low values (approximately equal to 65).(ABSTRACT TRUNCATED AT 400 WORDS)

Conformational analysis of the deoxyribofuranose ring in DNA by means of sums of proton-proton coupling constants: a graphical method

A graphical method is presented for the conformational analysis of the sugar ring in DNA fragments by means of proton-proton couplings. The coupling data required for this analysis consist of sums of couplings, which are referred to as sigma 1' (= J1'2' + J1'2''), sigma 2' (= J1'2' + J2'3' + J2'2''), sigma 2'' (= J1'2'' + J2''3' + J2'2'') and sigma 3' (= J2'3' + J2''3' + J3'4'). These sums of couplings correspond to the distance between the outer peaks of the H1', H2', H2'' and H3' [31P] resonances, respectively, (except for sigma 2' and sigma 2'' in the case of a small chemical shift difference between the H2' and H2'' resonances) and can often be obtained from 1H-NMR spectra via first-order measurement, obviating the necessity of a computer-assisted simulation of the fine structure of these resonances. Two different types of graphs for the interpretation of the coupling data are discussed: the first type of graph serves to probe as to whether or not the sugar ring occurs as a single conformer, and if so to analyze the coupling data in terms of the geometry of this sugar ring. In cases where the sugar ring does not occur as a single conformer, but as a blend of N- and S-type sugar puckers, the second type of graph is used to analyze the coupling data in terms of the geometry and population of the most abundant form. It is shown that the latter type of analysis can be carried out on the basis of experimental values for merely sigma 1',sigma 2' and sigma 2'', without any assumptions or restrictions concerning a relation between the geometry of the N- and S-type conformer. In addition, the question is discussed as to how insight can be gained into the conformational purity of the sugar ring from the observed fine structure of the H1' resonance. Finally, a comparison is made between experimental coupling data reported for single-stranded and duplex DNA fragments and covalent RNA-DNA hybrids on the one hand and the predicted couplings and sums of couplings presented in this paper on the other hand.

Influence of 2-aminoadenosine, A', on the conformational behaviour of d(T-A'-T-A'). A one-dimensional proton NMR study at 300 MHz and 500 MHz

Proton NMR studies at 300 MHz and 500 MHz in aqueous solution were carried out on the modified self-complementary tetranucleoside trisphosphate d(T-A'-T-A'), in which d(A') represents 2-aminodeoxyadenosine. The NMR spectra were observed at two sample concentrations over the temperature range 2-70 degrees C. Assignments, based on homonuclear decoupling and nuclear Overhauser enhancement (NOE) experiments, are given. The concentration dependence of the chemical shift vs temperature profiles was used to extract information concerning duplex formation. The proton-proton and proton-phosphorus coupling constants were obtained at four temperatures and yielded accurate conformational data on the sugar ring and on the back-bone angles beta, gamma and delta. From the observed line broadening, shift profiles, NOEs, and from the observation of imino proton resonances, it is concluded that the compound exists as a miniduplex at low temperature. Comparison of these observations with similar observations on the parent compound d(T-A-T-A) indicate that substitution of dA by dA' increases the tendency towards duplex formation. At low temperature the compound adopts a stacked B-DNA type structure: destacking occurs on raising the temperature. The sequence-dependent sugar ring geometry [Mellema, J.-R., Pieters, J. M. L., van der Marel, G. A., van Boom, J. H., Haasnoot, C. A. G. & Altona, C. (1984) Eur. J. Biochem. 143, 285-301] is present in this molecule. The conformational parameters of d(T-A'-T-A') and d(T-A-T-A) are quite similar, thus substitution of dA by dA' has no measurable influence on the geometry of the sugar ring or on the backbone angles beta, gamma and delta.

cis-diamminedichloroplatinum(II) induced distortion of a single and double stranded deoxydecanucleosidenonaphosphate studied by nuclear magnetic resonance

The structural distortion of a single- and a double-stranded decadeoxynucleotide upon binding of cis-PtCl2(NH3)2 was studied by 1H-NMR. After selective platination of d(T-C-T-C-G-G-T-C-T-C) (I) at the central d(-GpG-) site (resulting in I-Pt), several non-exchangeable base protons as well as H1', H2', H2" and H3' protons could be assigned by means of conventional NMR double-resonance techniques. Addition of the complementary decamer strand to I and I-Pt yielded the double-stranded III and III-Pt, respectively. All non-exchangeable base, H1', and most of the H2' and H2" protons in the two double stranded compounds could be assigned using 2D-chemical shift correlation (COSY) and nuclear Overhauser enhancement (NOESY) techniques. The double stranded compound III appears to adopt a B-DNA like structure. Comparison of NOEs and proton-proton coupling constants in the d(-GpG-).cisPt part in I-Pt and III-Pt reveals that their structure displays large similarity. Significant chemical shift changes (i.e. larger than 0.1 ppm) between III and III-Pt are restricted to the central four base pairs. It follows that the outer three base pairs, located on either side of the central four base pairs in III-Pt are likely to adopt a regular B-DNA type helix. The observed large upfield and downfield chemical shifts in the d(-CpGpG-) part of III with respect to III-Pt can be rationalized by describing the distortion of the double helix as a kink. A discussion of the observed physical effects upon platination of a double-stranded oligonucleotide is presented.

Sequence-dependent structural variation in single-helical DNA. Proton NMR studies of d(T-A-T-A) and d(A-T-A-T) in aqueous solution

The two deoxyribotetranucleoside triphosphates d(T-A-T-A) and d(A-T-A-T) were investigated in aqueous solution by one- and two-dimensional proton NMR at 300 and 500 MHz. It is demonstrated that both compounds occur predominantly in the single-helical form. Accurate coupling constants are obtained by computer simulation of several 500-MHz spectra. The data are interpreted in terms of N and S pseudorotational ranges. The geometry of the major S-type conformers displays a clear sequence dependence, as expressed by variation of the endocyclic backbone angle delta (C5'-C4'-C3'-O3'). A simple sum rule is proposed to predict delta variation in single-helical DNA fragments. Comparisons are made with other sequence-dependent geometries as observed in a double-helical B-DNA fragment in the crystalline state. Furthermore, one- and two-dimensional nuclear Overhauser effect (NOE) spectroscopy was carried out on d(T-A-T-A). An inventory is made of the observed intra- and inter-residue NOEs. The NOE data confirm the presence of a highly stacked single-helical conformation of d(T-A-T-A) in solution. No indications are found for the formation of a bulge-out structure as observed for analogous alternating purine-pyrimidine oligoribonucleotides.