| Nucleic Acids Research | Pages |
Theoretical calculations, synthesis and base pairing properties of oligonucleotides containing 8-amino-2[prime]-deoxyadenosine
Introduction
Materials And Methods
Molecular modeling
Theoretical calculations
Chemicals
8-Amino-2[prime]-deoxy-N2,N8-bis(dimethylaminomethyliden) adenosine (2)
8-Amino-2[prime]-deoxy-5[prime]-O-dimethoxytrityl-N2,N8-bis(dimethyl-aminomethyliden) adenosine-3[prime]-O-(2-cyanoethyl)-N,N-diisopropyl phosphoramidite (4)
Deprotection studies
Oligonucleotide synthesis
HPLC conditions
Melting experiments
Results And Discussion
QM calculations
MD and MD-TI calculations
Preparation of the phosphoramidite derivative and deprotection studies
Oligonucleotide synthesis
Base pairing properties
Conclusions
Acknowledgements
References
Theoretical calculations, synthesis and base pairing properties of oligonucleotides containing 8-amino-2[prime]-deoxyadenosine
ABSTRACT
INTRODUCTION
Several years ago, oligonucleotides that could bind on homopurine-homopyrimidine sequences of double-stranded DNA by forming triple helices were described. The formation of nucleic acid triple helices opens the way to designing sequence-specific DNA-binding molecules which may have therapeutic and diagnostic uses such as antigene therapy (1,2). Due to these potential uses of triple helices, a lot of effort has been put into the design and preparation of modified oligonucleotides in order to enhance triple helix stability (1,2). One of the most successful modifications was to replace natural bases with some modified bases such as 5-methyl-cytidine (3,4), 5-bromouracil (5,6), 5-aminouracil (7), N4-spermine-5-methylcytidine (8) or 5-methyl-2,6(1H,3H)-pyrimidinedione (9). Due to the lack of structural information, all chemical modifications were made after screening a large number of derivatives.
A detailed picture of the DNA triplex in solution, including the specific pattern of hydration which may be used for the design of new bases with triplex stabilization properties, was recently obtained by molecular dynamics simulations (10). Average triplex structures in Shields et al. (10) displayed the existence of a small groove [the minor-Major (mM) groove] in the surroundings of the C8 of adenines. Such a narrow groove has a defined pattern of hydration and a region of large electronegative molecular interaction potential (see fig. 3 in ref. 10). These structural features suggest that a notable increase in the stability of the triplex can be achieved by introducing a polar group near C8, which will replace highly ordered water molecules in the mM groove. The stabilization should increase if the group displays a positive electrostatic potential so that it can interact with the negative potential of the DNA in the mM groove. Analysis of a d(A·T·T) triad suggests that the introduction of an amino group at C8 might lead to the formation of an additional A·T (Hoogsteen) hydrogen bond (H-bond), which should add extra stabilization to the triple helix (Fig.
Figure 1. Schematic representation of A·T·T and M·T·T triads. The grooves generated in the structure of a triplex DNA are indicated following Shields et al. (10). The Insight-II computer package (13) was used for all the modeling. The Gaussian 94 (14) program was used for quantum mechanical (QM) calculations. Molecular dynamics (MD) and MD-thermodynamic integration (TI) calculations were made using AMBER5.0 (15). All calculations used the ORIGIN-2000 parallel supercomputer of the Centre de Supercomputació de Catalunya, as well as work stations in our laboratory.
MATERIALS AND METHODS
Molecular modeling
Theoretical calculations
QM calculations. The first stage in the theoretical study was to determine the suitability of the H-bond pattern shown in Figure
A standard B-DNA d(A·T)10 duplex and a DNA triplex of sequence d(A·T·T)10 were generated as follows. The duplex structure was constructed using Arnotts standard fibber parameters (19), as implemented in Insight-II (13). For the triplex structure a MD-averaged structure of the d(A·T·T)10 triplex was used as template (10). In both cases the fifth adenine was substituted by 8-aminoadenine using Insight-II, while the rest of the bases remained unchanged in the model.
Duplex and triplex structures were immersed in a water bath and, in order to obtain electroneutral systems, Na+ ions were added (10). The duplex system contains 3211 water molecules and 18 Na+ ions while the triplex system has 4207 water molecules and 27 Na+ ions, for a total of 10 177 and 13 607 atoms, respectively. 5000 cycles of minimization were applied to both systems prior to the heating and relaxation process during 130 ps. The heating and relaxation process consisted of seven stages: (i) 10 ps of MD (T = 100 K, P = 1 atm) with restraints of 100 kcal/mol Å2 in DNA and sodium; (ii) 10 ps of MD (T = 100 K, P = 1 atm) with restraints of 100 kcal/mol Å2 in DNA; (iii) 10 ps of MD (T = 300 K, P = 1 atm) with restraints of 100 kcal/mol Å2 in DNA; (iv) 25 ps of MD (T = 300 K, P = 1 atm) with restraints of 50 kcal/mol Å2 in DNA; (v) 25 ps of MD (T = 300 K, P = 1 atm) with restraints of 25 kcal/mol Å2 in DNA; (vi) 25 ps of MD(T = 300 K, P = 1 atm) with restraints of 10 kcal/mol Å2 in DNA; (vii) 25 ps of MD (T = 300 K, P = 1 atm) with restraints of 5 kcal/mol Å2 in DNA. After these stages the system was equilibrated [the starting model for triplex DNA had already been equilibrated for 1 ns (see ref. 10), so 100 ps of free MD were considered enough to equilibrate the system fully] through 100 (triplex) and 300 (duplex) ps of non-restrained MD at constant pressure (1 atm) and temperature (298°C).
Equilibrated duplex and triplex structures were the starting point for two MD-TI studies (20). In MD-TI calculations the free energy change relating to the M->A mutation in a given environment is computed as the reversible work necessary to mutate an 8-aminoadenine molecule into an A molecule. The mutation coordinate is then defined as the mixing parameter ([lambda]) which determines the force field parameters defining the mutated residue: Vl = [lambda]VM + (1 + [lambda])VA. These calculations give us a first approximation to the relative stability of the triplex arising from the substitution of a central adenine by 8-aminoadenine under physiological conditions. Simulations were always performed in the M->A direction using 21 windows. Each window consisted of 10 ps of equilibration and 10 ps of averaging for a total of 420 ps in each mutation study. Not only inter-group contributions to the free energy difference, but also the non-bonded intra-group contributions and the bond PMF corrections (21) were considered.
All simulations (equilibration MD and MD-TI) used periodic boundary conditions and the particle-mesh Ewal approach (22) to introduce long-range electrostatic effects. All bond lengths were frozen using SHAKE (23), which allows the use of a 2 fs time step for integration. In all cases the DNA structures were represented using AMBER-95 force field parameters (24,25), while water was represented by means of the well-known TIP3P model (24,25). Atomic changes for 8-aminoadenine were determined (data available on request) using the standard RESP strategy (26).
Chemicals
Protected phosphoramidites and DNA synthesis reagents were purchased from Perkin Elmer-Applied Biosystems (USA). Nucleosides were purchased from Pharma-Waldorf (Germany). Anhydrous solvents were obtained from SDS (France). Phosphodiesterase and alkaline phosphatase were obtained from Boehringer Mannheim (Germany). The rest of the chemicals were from Fluka (Switzerland) and Aldrich (USA) and were used without further purification. 8-Amino-2[prime]-deoxyadenosine (1) was prepared following previously described protocols (27,28).
8-Amino-2[prime]-deoxy-N2,N8-bis(dimethylaminomethyliden) adenosine (2)
8-Amino-2[prime]-deoxyadenosine (27,28) (1) (3.25 g, 12.2 mmol) was suspended in 200 ml of methanol and 8.2 ml of N,N-dimethyl-formamide dimethylacetal (61 mmol) were added. The mixture was stirred overnight at room temperature. The reaction mixture was concentrated to dryness and the residue was dissolved in 10% ethanol (EtOH) in dichloromethane (DCM) and purified on a silica gel column eluted with 15% EtOH in DCM. Yield 3.52 g (9.35 mmol, 76%). TLC (10% EtOH in DCM) Rf 0.29. UV (H2O) [lambda]max 267, 334 nm. 1H-NMR (DMSO-d6) [delta]H: 8.81 (1H, s), 8.66 (1H, s), 8.17 (1H, s), 6.57 (1H, t), 5.55 (1H, m), 5.25 (2H, m), 4.46 (1H, m), 3.86 (1H, m), 3.5-3.6 (2H, m), 3.2-3.05 (13H, 3s and 1m), 2.05 (1H, m). 13C-NMR (DMSO-d6) [delta]C: 157.6, 157.1, 156.0, 155.6, 151.6, 148.4, 124.7, 87.8, 82.9, 71.7, 62.6, 40.4 (under solvent peak), 37.1, 34.5, 34.4. Analysis of C16H24N8O3.½H2O: predicted, C 49.86, H 6.54, N 29.07; found, C 49.52, H 6.57, N 28.88.
8-Amino-2[prime]-deoxy-5[prime]-O-dimethoxytrityl-N2,N8-bis(dimethyl-aminomethyliden) adenosine-3[prime]-O-(2-cyanoethyl)-N,N-diisopropyl phosphoramidite (4)
Compound 2 (1.0 g, 2.65 mmol) was dissolved in 50 ml of dry pyridine and 1.08 g of dimethoxytrityl chloride (3.18 mmol) were added. The mixture was stirred for 2 h at room temperature. To the reaction mixture 1 ml of methanol was added and the solution was concentrated to dryness. The residue was dissolved in DCM and washed with 5% NaCO3H and saturated NaCl aqueous solution. The organic phase was dried with anhydrous Na2SO4 and concentrated to dryness. The residue was purified by chromato-graphy on silica gel. The column was packed with silica gel using a 1% triethylamine solution in DCM. The product was eluted with a 0-4% methanol gradient in DCM. Compound 3 was obtained as a white foam (1.5 g, 2.21 mmol, 83%). TLC (5% EtOH/DCM)Rf 0.6.1H-NMR (Cl3CD) [delta]H: 8.71 (1H, s), 8.68 (1H, s), 8.14 (1H, s), 7.5-7.1 (9H, m), 6.8-6.6 (5H, m), 4.85 (1H, m), 4.18 (1H, m), 3.69 (6H, s), 3.4 (2H, m), 3.16, 3.07, 3.02, 2.93 (13H, 4 s and 1m), 2.2 (1H, m). 13C-NMR (Cl3CD) [delta]C: 158.1, 157.5, 156.8, 156.7, 155.3, 152.1, 149.2, 144.8, 136.0, 129.9, 128.0, 127.4, 126.4, 124.8, 112.8, 85.9, 85.3, 82.4, 73.2, 64.4, 55.0, 40.8, 40.6, 36.4, 34.8, 34.6.
Compound 3 (1.16 g, 1.7 mmol) and diisopropylethylamine (0.9 ml, 5.1 mmol) were dissolved in 20 ml of dry acetonitrile. The solution was cooled with an ice bath and 0.57 ml of 2-cyanoethoxy-N,N-diisopropylamino-chlorophosphine (2.55 mmol) were added. After 1 h of magnetic stirring at room temperature, the mixture was concentrated to dryness. The residue was dissolved in DCM and washed with 5% NaCO3H aqueous solution. The organic phase was dried with anhydrous Na2SO4 and concentrated to dryness. The residue was purified on silica gel. The column was packed using 1% triethylamine solution in DCM/hexane 9:1 and the product was eluted with DCM/hexane 9:1. Compound 4 was obtained as a white foam. Yield 1.4 g (1.59 mmol, 93%). TLC (DCM/hexane 9:1) Rf 0.8. 1H-NMR as described in Kawai et al. (12). 31P-NMR (Cl3CD) [delta]P: 144.96. Analysis of C46H59N10O6P: predicted, C 62.86, H 6.77, N 15.93; found, C 62.58, H 6.80, N 15.64.
Deprotection studies
Aliquots of dinucleotide supports containing 8-amino-deoxyadenosine (5[prime]-MT-3[prime]) were treated with concentrated ammonia either at room temperature or at 55°C. At different intervals, the ammonia solutions were filtered, concentrated to dryness and analyzed by HPLC. HPLC conditions are described below. A 30 min linear gradient from 0 to 60% B was used. Elution at 260 nm was followed.
Dimer MT. HPLC: retention time 17.2 min; UVmax 269 nm. Mass spectra found 570.4 (expected for C20H27N8O10P: 570.4).
Dimer MT containing a N6,N8-bis-dimethylaminomethyliden (dmf) group. HPLC: retention time 18.8 min. UVmax 270, 305 nm. Mass spectra found 625.3 (expected for C23H32N9O10P: 625.5).
Dimer MT containing two dmf groups. HPLC: retention time: 20 min. UVmax 266, 335 nm. Mass spectra found 679.3 (expected for C26H37N10O10P: 680.5).
Oligonucleotide synthesis
Oligonucleotides were prepared on an automatic DNA synthesizer using standard 2-cyanoethyl phosphoramidites and the modified phosphoramidite 4. Complementary pentadecamers and undecamers containing natural bases were also prepared using commercially available chemicals and following standard protocols. After the assembly of the sequences, the oligonucleotide supports were treated with 32% aqueous ammonia at 55°C for 16 h. Ammonia solutions were concentrated to dryness and the products were purified by reverse phase HPLC. Dimer A and pentamer B were prepared without the last DMT group (DMT-off protocol) on 1 µmol scale. The remaining oligonucleotides were synthesized on 0.2 µmol scale and with the last DMT group at the 5[prime]-end (DMT-on protocol) to help reverse phase purification. All purified products presented a major peak which was collected and analyzed by snake venom phosphodiesterase and alkaline phosphatase digestion followed by HPLC analysis of the nucleosides (HPLC conditions B). The retention times of nucleosides obtained after enzyme digestion were: dC, 4.6 min; dG, 8.4 min; T, 9.2 min; dA, 17.6 min; 8-amino-2[prime]-deoxyadenosine, 21.2 min. Mass spectra: sequence B, found 1501.8 (expected for C49H63N21O27P4 1501.9); sequence C, found 4579 (expected for C146H185N59O86P14 4575.5). Yield (OD units at 260 nm after HPLC purification, 0.2 µmol): sequence B, 6.3 OD; sequence C, 10.4 OD; sequence D, 16.8 OD; sequence E, 16.9 OD; sequence F, 13.5 OD; sequence G, 9.6 OD.
HPLC conditions
HPLC solutions were as follows: solvent A, 5% ACN in 100 mM triethylammonium acetate (pH 6.5); solvent B, 70% ACN in 100 mM triethylammonium acetate (pH 6.5). For analytical runs the following conditions were used: column, Nucleosil 120C18, 250 × 4 mm; flow rate, 1 ml/min; condition A, a 40 min linear gradient from 0 to 75% B; condition B, a 20 min linear gradient from 0 to 20% B. For preparative runs the following conditions were used: column, PRP-1 (Hamilton), 250 × 10 mm; flow rate, 3 ml/min; 30 min linear gradient from 10 to 80% B (DMT-on) or a 30 min linear gradient from 0 to 50% B (DMT-off)
Melting experiments
Melting experiments of pentadecamer duplexes were carried out by mixing equimolar amounts of two pentadecamer strands dissolved in a solution containing 0.15 M NaCl, 0.05 N Tris-HCl buffer, pH 7.5. Duplexes were annealed by slow cooling from 80 to 4°C. UV absorption spectra and melting curves (absorbance versus temperature) were recorded in 1 cm pathlength cells using a Varian Cary 13 spectrophotometer having a temperature controller with a programmed temperature increase of 0.5°C/min. Melts were run on duplex concentrations of 4 µM at 260 nm.
Melting experiments with triple helix were performed as follows. Solutions of equimolar amounts of the hairpin oligo-nucleotide (h26) and the 11mer (s11) were mixed in the appropriate buffer. The solutions were heated to 80°C, allowed to cool slowly to room temperature and then samples were kept in the refrigerator overnight. UV absorption spectra and melting experiments (absorbance versus temperature) were recorded in 1 cm pathlength cells using a spectrophotometer having a temperature controller with a programmed temperature increase of 0.5°C/min. Melts were run on duplex concentrations of 4 µM at 270 nm.
RESULTS AND DISCUSSION
QM calculations
The first stage in the theoretical study was to determine the suitability of the H-bond pattern shown in Figure
A
![]() |
B
![]() |
Figure 2. B3LYP/6-31G(d) optimized geometries for Hoogsteen A·T and M·T pairs in the vacuum.
The stabilization energy determined for the A·T Hoogsteen pairing is -12.7 kcal/mol (after distortion and BSSE corrections) and a value of -17.9 kcal/mol is reached for the M·T Hoogsteen pairing (after distortion and BSSE corrections). The extra H-bond present in the M·T (Hoogsteen) dimer is very strong and it seems to contribute to the stability of the purine-pyrimidine (Hoogsteen) dimer. Remarkably, the difference in stability between the A·T (Hoogsteen) and M·T (Hoogsteen) dimers (-5.2 kcal/mol from Fig.
M·T (Hoogsteen) binding energy determined from B3LYP/6-31G(d) calculations cannot be compared directly with experimental data to verify the correctness of the theoretical estimate. However, an early study of heterodimerization of adenine and thymine in the gas phase suggests (19) an enthalpy of dimerization of -13 kcal/mol, which is consistent with our theoretical estimation of -12.7 kcal/mol and gives additional support to the B3LYP/6-31G(d) calculations presented in this paper.
B3LYP/6-31G(d) calculations were also performed to investigate the role of the extra 8-amino group in the Watson-Crick binding of adenine and thymine. Geometry optimization (data not shown) is not able to find any significant change in dimer geometry. Thus, for the Watson-Crick A:T dimer, the H-bonds are 2.95 Å for the N6-O4 and 2.88 Å for the N1-N3 bond. Almost the same values (2.96 and 2.87 Å) are found for the M·T Watson-Crick dimer.
Dimerization energy for the Watson-Crick pairing of A·T is -12.1 kcal/mol, which compares with -12.7 kcal/mol found for the most stable Hoogsteen pairing. The Watson-Crick pairing of M·T has a stabilization energy of -11.7 kcal/mol, clearly less than the Hoogsteen binding of M·T (-17.9 kcal/mol), reflecting the fact that the Hoogsteen binding is stabilized for M·T by one more H-bond than the Watson-Crick binding. Comparison of M·T and A·T Watson-Crick pairing reveals that the presence of the amino group at position 8 leads to a slight decrease (0.4 kcal/mol) in the stability of the pairing, due probably to a decrease in the acidity of the purine N6 group.
MD and MD-TI calculations
Quantum mechanical calculations suggest that the replacement of A by M can lead to the formation of an extra H-bond which stabilizes in vacuo the purine·pyrimidine Hoogsteen H-bond. However, 8-aminoadenine behavior in solution cannot be completely described with this method because water greatly hinders the formation of H-bonds. For a more direct estimate, MD-TI calculations were performed. MD-TI calculations can provide a direct estimate of the influence of the replacement of an A by M on triplex stability, by computing the differences in the reversible work necessary to replace an M by A in a duplex and in a triplex DNA (Fig.
Figure 3. Schematic representation of the thermodynamic cycle used to compute the stabilization of the triplex due to the replacement of an adenine by 8-aminoadenine. Note that the magnitude in a physical sense is [Delta]G(A) - [Delta]G(B), while the magnitude computed from simulations is [Delta]G(1) - [Delta]G(2). Change of H8 by an amino group does not lead to any significant difference in duplex or triplex geometry. The Watson-Crick and Hoogsteen H-bonds were preserved during all the simulations. The average M·T (WC) and M·T (Hoogsteen) H-bond distances found in the MD are close to the optimum values determined from gas phase calculations. The amino group is placed in the mM groove of the triplex and seems well hydrated. Mutations were done in the 8-aminoadenine ([lambda] = 1) to adenine ([lambda] = 0) direction. The free energy profiles for the duplex and triplex were smooth, without apparent discontinuity. When the free energy values were computed with the first 10 ps of each window (equilibration part) or with the last 10 ps (averaging part), almost the same free energy profiles were obtained (to <0.1 kcal/mol), confirming the suitability of the sampling protocol used to compute the free energy profile (Fig. Figure 4. Free energy profile for the mutation from 8-aminoadenine ([lambda] = 1) to adenine ([lambda] = 0) in the triplex and duplex systems. Comparison of free profiles in the duplex and triplex simulations clearly demonstrates that the change from M is 2.7 kcal/mol more favorable in a duplex DNA than in a triplex DNA (Figs Figure 5. Differential (triplex-duplex) free energy profile for the mutation of 8-aminoadenine ([lambda] = 1) into adenine ([lambda] = 0). Moreover, our molecular dynamic simulation results suggest that the change of a single A residue to an 8-aminoadenine in d(A)10 does not introduce dramatic changes in the structural properties of the helix. The 8-amino group does not induce any change of conformation in the glycosidic bond and the amino group interacts nicely with water and the ribose moiety (O4[prime] and O5[prime]). These results suggest that the A->M mutation is conservative in terms of DNA duplex structure. These results confirm QM calculations noted above, which demonstrate the similar stability of Watson-Crick A·T and M·T pairings. In a similar manner, the introduction of the 8-amino group into a step of the d(A·T·T)10 triplex does not lead to significant changes in the structure of the triplex. MD simulations during 500 ps allowed us to define an average structure (for details see ref. 10) which is very similar to that obtained for the pure d(A·T·T)10 triplex. As in the duplex, the extra amino group seems reasonably solvated and establishes several H-bonding interactions with the ribose moiety and, to a lesser extent, with the phosphate group. In summary, MD suggests that a single A->M mutation in a poly(A) tract does not lead to dramatic alterations in the structure of the triplex.
Preparation of the phosphoramidite derivative and deprotection studies
The synthesis of the phosphoramidite derivative of M is illustrated in Figure
Figure 6. Preparation of 8-amino-2[prime]-deoxyadenosine phosphoramidite. The stability of the dmf groups to ammonia was analyzed on the dimer 5[prime]-MT-3[prime]. The dimer was assembled on controlled pore glass using standard phosphoramidite protocols. Aliquots of dinucleotide supports were treated with concentrated ammonia at room temperature and at 55°C. HPLC analysis of the resulting products followed by mass spectrometry analysis of the products showed complete deprotection at 55°C in <3 h. At room temperature only one dmf group was removed rapidly and the dimer containing a single dmf group was slowly converted to the unprotected dimer. Other minor products due to exchange of the dmf groups to acetyl groups during the capping reaction were also observed. All these intermediate compounds were not detected when deprotection was performed at 55°C. Therefore, for the rest of the sequences containing M, all deprotection reactions were carried out overnight at 55°C.
Oligonucleotide synthesis
Oligonucleotide sequences containing M ranging from 5 to 26 bases (Table 1) were prepared using phosphoramidite chemistry on an automatic DNA synthesizer. Coupling efficiency of 8-aminoadenine phosphoramidite 4 was similar to the commercially available phosphoramidites. After standard deprotection with concentrated ammonia, the products were purified by reverse phase HPLC using the DMT-on and DMT-off protocols. In all cases a major peak was obtained and collected. The purified oligonucleotides were analyzed by mass spectrometry (electrospray and MALDI-TOF) and by enzyme digestion followed by HPLC analysis of the resulting nucleosides. In all cases the purified oligonucleotides were obtained in good yield and had the correct nucleoside composition and the expected mass.
Table 1.
| Sequence (5[prime]->3[prime]) | ||
| A | MT | |
| B | AGMCT | |
| C | GCAATGGAMCCTCTA | |
| D | h261 | GAAGGMGGAGATTTTTCTCCTCCTTC |
| E | h262 | GAAGGMGGMGATTTTTCTCCTCCTTC |
| F | h263 | GAMGGMGGMGATTTTTCTCCTCCTTC |
| G | h264 | GMMGGMGGMGMTTTTTCTCCTCCTTC |
Base pairing properties
The hybridization properties of 8-aminoadenine were measured spectrophotometrically on a duplex formed by two complementary pentadecamers in which M is paired with the four natural bases (Table 2). 8-Aminoadenine forms the strongest base pairs with T as expected. This base pair is slightly less stable (1.6°C) than the natural A·T. Compared with adenine, 8-aminoadenine base pairs are between 1.6 and 6.4°C less stable than the A·T base pair, while the same base pairs with A were from 0 to 9.6°C less stable than the A·T base pair. There is a tendency to level the melting temperatures of mismatches that comes from a small decrease in the stability of the M·T (1.6°C compared to A·T) base pair and an increase in the M·C (3°C compared to A·C) and M·A (2.9°C compared to A·A) base pairs. This tendency to similar melting temperatures for the different mismatches may assist in the design of degenerate primers or probes. This possibility was investigated with primers containing multiple 8-aminoadenine mismatches for PCR reactions. Unfortunately, no special properties were found in the oligonucleotides carrying M as against the same primers carrying A (data not shown), so precluding its use as a universal base.
Table 2.
| Xa | Ya | ||
| M | A | G | |
| A | 52.6 | 49.7 | 54.2 |
| C | 54.5 | 51.5 | 60.9 |
| T | 57.7 | 59.3 | 51.3 |
| G | 54.9 | 54.8 | 56.8 |
QM calculations suggest that the Watson-Crick pairing of M·T is slightly less favorable than A·T. This and the large desolvation penalty upon binding of 8-aminoadenine justify that duplexes containing M are less stable than those with A. Preliminary NMR and circular dichroism analysis of duplexes containing M·T base pairs show no major changes in the duplex structure on the introduction of 8-aminoadenine (data not shown). In summary, neither experimental nor theoretical results suggest that changes of adenine to M cause any dramatic effect on duplex structure.
Triple helix stabilization properties of M were investigated using a triple helix model formed by a self-complementary hairpin of 26 bases (h26) and an all-pyrimidine single-strand oligonucleotide (s11) described previously (3). Sequential replacement of A by M in triple helix results in a stabilization of 5-8°C per substitution in the range from pH 5.5 to 7.0, except when the number of M substitutions goes from three to five residues (Table 3). In any case, the observation of triple helical structures at neutral pH becomes possible only when the M-modified base is used. The pH dependence of triplex stability is due to the presence of G·C·C+ triads. These results corroborate those of Kawai et al. (12) obtained on a three-strand triple helix. The degree of stabilization and its additive character are consistent with the existence of an extra hydrogen bond between the 8-amino group of 8-aminoadenine and the 2-keto group of T (Fig.
CONCLUSIONS
In summary, we predicted by theoretical calculations that the introduction of an amino group at position 8 of adenine would broadly stabilize triple helix formation, due to the combined effect of the gain in one Hoogsteen purine-pyrimidine H-bond and the ability of the amino group to be integrated into the spine of hydration located in the mM groove of the triplex structure. Oligonucleotides carrying M residues were synthesized using protected phosphoramidite 4. As we predicted, these oligonucleotides form very stable triple helices. As a matter of fact, the degree of stabilization is one of the highest reported for a non-natural base analog and it could be further increased by adding other reported stabilizing molecules such as intercalating compounds (1,2). The high degree of stabilization obtained by this simple substitution is especially relevant to the development of new applications based on triple helix formation such as structural studies, DNA-based diagnostic tools and antigene therapy (1,2). However, it should not be overlooked that applications using DNA from natural sources such as control of gene expression (1) may not benefit from the use of 8-aminoadenine. Since the modified purine is on the target strand, the direct applicability of this modification in certain applications may be limited. One possible way to overcome this problem may be the preparation of oligonucleotides containing the purine Watson-Crick strand (with 8-aminoadenine) linked with the pyrimidine Hoogsteen strand. If the polarity of the strands is correct, these oligonucleotides may bind the pyrimidine Watson-Crick strand and leave the natural purine Watson-Crick strand unpaired due to the exceptionally high stability of 8-aminoadenine triplexes. Experiments along these lines are currently being undertaken.
Table 3.
![]() |
|||||
| Sequence | Mb | pH 5.5 | pH 6.0 | pH 6.5 | pH 7.0 |
| h26:s11 | 0 | 40 | ~20 | - | - |
| h261:s11 | 1 | 48.3 | 30.9 | 24 | - |
| h262:s11 | 2 | 56 | 41.8 | 32.4 | 22.2 |
| h263:s11 | 3 | 65.2 | 48.7 | 39.6 | 26.9 |
| h264:s11 | 5 | 68.3 | 51.9 | 43.5 | 32.2 |
Our results demonstrate the use of state-of-the-art quantum mechanics, molecular dynamics and thermodynamic integration techniques in the analysis and design of polynucleotidic sequences. Overall, our results suggest that introduction of amino groups at position 8 of purines may lead to new products with interesting triple helix stabilization properties. Just recently, oligonucleotides carrying 8-aminoguanine have also been shown to form a very stable triple helix (29), which indicates that, indeed, replacement at position 8 of purines by amino groups is key to triple helix stabilization.
See supplementary material available in NAR Online (40 kb PDF file).
ACKNOWLEDGEMENTS
We thank M. Wiersma for valuable suggestions. We also thank Dr Matthias Mann and his group for mass spectroscopic measurements. This research was partially supported by the Centre de Supercomputació de Catalunya (CESCA, Molecular Recognition Project) and the Spanish Dirección General de Investigación Científica y Técnica (DGICYT, PB96-1005).
REFERENCES
This article has been cited by other articles:
This page is run by Oxford University Press, Great Clarendon Street, Oxford OX2 6DP, as part of the OUP Journals
Comments and feedback: jnl.info{at}oup.co.uk
Last modification: 10 Apr 1999
Copyright©Oxford University Press, 1999.
![]()
CiteULike
Connotea
Del.icio.us What's this?
![]()
![]()

![]()
![]()
![]()
A. Avino, M. Frieden, J. C. Morales, B. Garcia de la Torre, R. Guimil Garcia, F. Azorin, J. L. Gelpi, M. Orozco, C. Gonzalez, and R. Eritja
Properties of triple helices formed by parallel-stranded hairpins containing 8-aminopurines
Nucleic Acids Res.,
June 15, 2002;
30(12):
2609 - 2619.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
E. Cubero, R. Guimil-Garcia, F. J. Luque, R. Eritja, and M. Orozco
The effect of amino groups on the stability of DNA duplexes and triplexes based on purines derived from inosine
Nucleic Acids Res.,
June 15, 2001;
29(12):
2522 - 2534.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
R. Soliva, R. Guimil Garcia, J. R. Blas, R. Eritja, J. L. Asensio, C. Gonzalez, F. J. Luque, and M. Orozco
DNA-triplex stabilizing properties of 8-aminoguanine
Nucleic Acids Res.,
November 15, 2000;
28(22):
4531 - 4539.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
V. S. Rana and K. N. Ganesh
Recognition of 5-aminouracil (U#) in the central strand of a DNA triplex: orientation selective binding of different third strand bases
Nucleic Acids Res.,
March 1, 2000;
28(5):
1162 - 1169.
[Abstract]
[Full Text]
[PDF]
![]()
This Article ![]()
![]()
Abstract
![]()
Print PDF (110K)
![]()
Supplementary Data
![]()
Alert me when this article is cited
![]()
Alert me if a correction is posted
![]()
Services ![]()
![]()
Email this article to a friend
![]()
Similar articles in this journal
![]()
Similar articles in ISI Web of Science
![]()
Similar articles in PubMed
![]()
Alert me to new issues of the journal
![]()
Add to My Personal Archive
![]()
Download to citation manager
![]()
Search for citing articles in:
ISI Web of Science (8)
![]()
Request Permissions ![]()
Commercial Re-use Guidelines
for Open Access NAR Content
![]()
Google Scholar ![]()
![]()
Articles by Garcia, R. G.
![]()
Articles by Orozco, M.
![]()
Search for Related Content
![]()
PubMed ![]()
![]()
PubMed Citation
![]()
Articles by Garcia, R. G.
![]()
Articles by Orozco, M.
![]()
Social Bookmarking ![]()
![]()
What's this?


