Base flipping of the thymine dimer in duplex DNA

Exposure of two adjacent thymines in DNA to UV light of 260-320 nm can result in the formation of the cis,syn-cyclobutane pyrimidine dimer (CPD). The structure of DNA containing an intrahelical CPD lesion has been previously studied experimentally and computationally. However, the structure of the extrahelical, flipped-out, CPD lesion, which has been shown to be the structure that binds to the CPD repair enzyme, DNA photolyase, has yet to be reported. In this work the structure of both the flipped-in and the flipped-out CPD lesions in duplex DNA is reported. These structures were calculated using 8 ns molecular dynamics (MD) simulations. These structures are then used to define the starting and ending points for the base-flipping process for the CPD lesion. Using a complex, two-dimensional pseudodihedral coordinate, the potential of mean force (PMF) for the base-flipping process was calculcated using novel methodology. The free energy of the flipped-out CPD is roughly 6.5 kcal/mol higher than that of the flipped-in state, indicating that the barrier to flipping out is much lower for CPD than for undamaged DNA. This may indicate that the flipped-out CPD lesion may be recognized by its repair enzyme, DNA photolyase, whereas previous studies of other damaged, as well as nondamaged, bases indicate that they are recognized by enzymes in the intrahelical, flipped-in state.     

The extent of DNA deformation in DNA photolyase-substrate complexes: a solution state fluorescence study.

Cyclobutylpyrimidine dimers (CPDs) are the major UV photoproduct formed in DNA containing adjacent pyrimidines. These lesions can be repaired by DNA photolyase, a flavoprotein that utilizes blue light in a direct reversal of the cyclobutane ring. Previous studies have shown that the CPD is base flipped into the protein, with concomitant disruption of the substrate around the CPD. In this study, we use a fluorescent cytidine analog, pyrrolo-dC (PC), to probe how far base flipping propagates along the duplex. From these measurements, the degree of base destacking in the two bases flanking the adenines opposing the CPD appears to be minimal, which was consistent with the protein:substrate crystal structure. Fluorescence-detected melting temperatures for duplexes with and without a CPD were obtained, suggesting that a 5'-pyrimidine-PC-purine-3' motif is more stable than the 5'-purine-PC-pyrimidine-3' motif. This stability trend was reflected in the fluorescence intensities of ss-PC oligos but not for duplexes. The melting point depression due to the PC probe was found to be comparable to other popular fluorescent base analogs.     

Stability and structure of RNA duplexes containing isoguanosine and isocytidine.

Isoguanosine (iG) and isocytidine (iC) differ from guanosine (G) and cytidine (C), respectively, in that the amino and carbonyl groups are transposed. The thermodynamic properties of a set of iG, iC containing RNA duplexes have been measured by UV optical melting. It is found that iG-iC replacements usually stabilize duplexes, and the stabilization per iG-iC pair is sequence-dependent. The sequence dependence can be fit to a nearest-neighbor model in which the stabilities of iG--iC pairs depend on the adjacent iG--iC or G--C pairs. For 5'-CG-3'/3'-GC-5' and 5'-GG-3'/3'-CC-5' nearest neighbors, the free energy differences upon iG-iC replacement are smaller than 0.2 kcal/mol at 37 degrees C, regardless of the number of replacements. For 5'-GC-3'/3'-CG-5', however, each iG--iC replacement adds 0.6 kcal/mol stabilizing free energy at 37 degrees C. Stacking propensities of iG and iC as unpaired nucleotides at the end of a duplex are similar to those of G and C. An NMR structure is reported for r(CiGCGiCG)(2) and found to belong to the A-form family. The structure has substantial deviations from standard A-form but is similar to published NMR and/or crystal structures for r(CGCGCG)(2) and 2'-O-methyl (CGCGCG)(2). These results provide benchmarks for theoretical calculations aimed at understanding the fundamental physical basis for the thermodynamic stabilities of nucleic acid duplexes.     

Long-range cooperativity in molecular recognition of RNA by oligodeoxynucleotides with multiple C5-(1-propynyl) pyrimidines.

A heptamer composed of C5-(1-propynyl) pyrimidines (Y(p)'s) is a potent and specific antisense agent against the mRNA of SV40 large T antigen (Wagner, R. W.; Matteucci, M. D.; Grant, D.; Huang, T.; Froehler, B. C. Nat. Biotechnol. 1996, 14, 840-844). To characterize the role of the propynyl groups in molecular recognition, thermodynamic increments associated with substitutions in DNA:RNA duplexes, such as 5'-dCCUCCUU-3':3'-rGAGGAGGAAAU-5', have been measured by UV melting experiments. For nucleotides tested, an unpaired dangling end stabilizes unmodified and propynylated duplexes similarly, except that addition of a 5' unpaired rA is 1.4 kcal/mol more stabilizing on the propynylated, PODN:RNA, duplex than on the DNA:RNA duplex. Free energy increments for addition of single propynyl groups range from 0 to -4.0 kcal/mol, depending on the final number and locations of substitutions. A preliminary model for predicting the stabilities of Y(p)-containing hybrid duplexes is presented. Eliminating one amino group, and therefore a hydrogen bond, by substituting inosine (I) for guanosine (G), to give 5'-dC(p)C(p)U(p)C(p)C(p)U(p)U(p)-3':3'-rGAGIAGGAAAU-5', destabilizes the duplex by 3.9 kcal/mol, compared to 1.7 kcal/mol for the same change within the unpropynylated duplex. This 2.2 kcal/mol difference is eliminated by removing a single propynyl group three base pairs away. CD spectra suggest that single propynyl deletions within the PODN:RNA duplex have position-dependent effects on helix geometry. The results suggest long-range cooperativity between propynyl groups and provide insights for rationally programming oligonucleotides with enhanced binding and specificity. This can be exploited in developing technologies that are dependent upon nucleic acid-based molecular recognition.     

Parallel DNA double helices incorporating isoG or m5isoC bases studied by FTIR, CD and molecular modeling

FTIR spectroscopy has been used to follow the formation of parallel stranded DNA duplexes incorporating isoG or m5isoC bases and determine their base pairing scheme. The results are discussed in comparison with data concerning anti-parallel duplexes with comparable base composition and sequence. In duplexes containing A-T and isoG-C or m5isoC-G base pairs shifts of the thymine C2=O2 and C4=O4 carbonyl stretching vibrations (to lower and higher wavenumbers, respectively, when compared to their positions in classical cis Watson-Crick (WC) base pairs) reflect the formation of trans Watson-Crick A-T base pairs. All carbonyl groups of cytosines, m5isocytosines, guanines and isoguanines are found to be involved in hydrogen bonds, indicative of the formation of isoG-C and m5isoC-G base pairs with three hydrogen bonds. Molecular modeling shows that both structures form regular right handed helices with C2'endo sugar puckers. The role of the water content on the helical conformation of the parallel duplexes has been studied by FTIR and CD. It is found that a conformational transition similar to the B --> A transition observed for anti-parallel duplexes induced by a decrease of the water content of the samples can occur for these parallel duplexes. Their helical flexibility has been evidenced by FTIR studies on hydrated films by the emergence of absorption bands characteristic of A type geometry, in particular by an S-type --> N-type repuckering of the deoxyribose. All sugars in the parallel duplex with alternating d(isoG-A)/d(C-T) sequence can adopt an N-type geometry in low water content conditions. The conformational transition of the parallel hairpin duplex with alternating d(isoG-A)/d(C-T) sequence was followed by circular dichroism in water/trifluoroethanol solutions and its free energy at 0 degrees C was estimated to be 6.6 +/- 0.3 kcal mol(-1).     

Hydrophobic, Non-Hydrogen-Bonding Bases and Base Pairs in DNA

We report the properties of hydrophobic isosteres of pyrimidines and purines in synthetic DNA duplexes. Phenyl nucleosides 1 and 2 are nonpolar isosteres of the natural thymidine nucleoside, and indole nucleoside 3 is an analog of the complementary purine 2-aminodeoxyadenosine. The nucleosides were incorporated into synthetic oligodeoxynucleotides and were paired against each other and against the natural bases. Thermal denaturation experiments were used to measure the stabilities of the duplexes at neutral pH. It is found that the hydrophobic base analogs are nonselective in pairing with the four natural bases but selective for pairing with each other rather than with the natural bases. For example, compound 2 selectively pairs with itself rather than with A, T, G, or C; the magnitude of this selectivity is found to be 6.5-9.3 °C in Tm or 1.5-1.8 kcal/mol in free energy (25 °C). All possible hydrophobic pairing combinations of 1, 2, and 3 were examined. Results show that the pairing affinity depends on the nature of the pairs and on position in the duplex. The highest affinity pairs are found to be the 1-1 and 2-2 self-pairs and the 1-2 heteropair. The best stabilization occurs when the pairs are placed at the ends of duplexes rather than internally; the internal pairs may be destabilized by imperfect steric mimicry which leads to non-ideal duplex structure. In some cases the hydrophobic pairs are significantly stabilizing to the DNA duplex; for example, when situated at the end of a duplex, the 1-1 pair is more stabilizing than a T-A pair. When situated internally, the affinity of the 1-1 pair is the same as, or slightly better than, the analogous T-T mismatch pair, which is known to have two hydrogen bonds. The studies raise the possibility that hydrogen bonds may not always be required for the formation of stable duplex DNA-like structure. In addition, the results point out the importance of solvation and desolvation in natural base pairing, and lend new support to the idea that hydrogen bonds in DNA may be more important for specificity of pairing than for affinity. Finally, the study raises the possibility of using these or related base pairs to expand the genetic code beyond the natural A-T and G-C pairs.     

Hydration regulates thermodynamics of G-quadruplex formation under molecular crowding conditions

The effect of molecular crowding on the structure and stability of biomolecules has become a subject of increasing interest because it can clarify how biomolecules behave under cell-mimicking conditions. Here, we quantitatively analyzed the effects of molecular crowding on the thermodynamics of antiparallel G-quadruplex formation via Hoogsteen base pairs and of antiparallel hairpin-looped duplex (HP duplex) formation via Watson-Crick base pairs. The free energy change at 25 degrees C for G-quadruplex formation decreased from -3.5 to -5.5 kcal mol(-1) when the concentration of poly(ethylene glycol) 200 was increased from 0 to 40 wt %, whereas that of duplex formation increased from -9.8 to -6.9 kcal mol(-1). These results showed that the antiparallel G-quadruplex is stabilized under molecular crowding conditions, but that the HP duplex is destabilized. Moreover, plots of stability (ln K(obs)) of the DNA structures versus water activity (ln a(w)) demonstrated that the ln K(obs) for G-quadruplex formation decreased linearly as the ln a(w) increased, whereas that for duplex formation increased linearly with the increase in ln a(w), suggesting that the slope approximately equals the number of water molecules released or taken up during the formation of these structures. Thus, molecular crowding affects the thermodynamics of DNA structure formation by altering the hydration of the DNA. The stabilization of the DNA structures with Hoogsteen base pairs and destabilization of DNA structures with Watson-Crick base pairs under molecular crowding conditions lead to structural polymorphism of DNA sequences regulated by the state of hydration.     

The possibility of the decomposition of 2'-deoxyribose moiety of thymidine induced by the low energy electron attachment

To evaluate the possibility of the decomposition of 2-deoxyribose moiety of thymidine induced by low energy electrons (LEE) attachment, the transition states and the energy barriers of the bond breaking processes of the ribose of the nucleoside have been studied theoretically by applying the density functional theory with the double zeta basis sets (DZP++). The energy barriers for the breakage of the C-C bonds (C(1')-C(2'), C(2')-C(3'), C(3')-C(4'), and C(4')-C(5')) of the ribose group of the radical anion of thymidine are found to be high (ca. 42-57 kcal/mol). The total energies of the C-C bond-broken products are significantly higher than that of the radical anion dT(*-). The decomposition of dT(*-) through the C-C bond rupture is unlikely to take place. The rupture of the C(1')-O(4') bond of dT(*-) needs an activation energy as low as 10.4 kcal/mol. However, the reversed reaction (C(1')-O(4') bond formation) needs the activation energy low as 0.3 kcal/mol. Therefore, the intermediate product LM1(C1')-(O4') is unlikely to be stable and the C(1')-O(4') bond-broken is not favored. The activation energy of the C(4')-O(4') bond rupture process amounts to 20.5 kcal/mol. The total energy of the C(4')-O(4') bond broken product is about 6.5 kcal/mol lower than that of the reactant dT(*-). The subsequent N1-glycosidic bond breaking process is found to have a very low energy barrier. Therefore, the LEE-induced base release through the C(4')-O(4') bond rupture might be a possible pathway.     

Are radical cation states delocalized over GG and GGG hole traps in DNA?

Recently, Kurnikov et al. (J. Phys. Chem. B 2002, 106, 7) have shown that solvation of DNA duplexes destabilizes holes of sizes larger than three base pairs. In this paper, we consider the effects of solvation and internal reorganization on the hole charge distribution in sequences 5'-X-GG-Y-3'. Radical cation states in DNA are found to be localized to a single guanine site independent of the nature of adjacent base pairs.     

Conserved patterns in backbone torsional changes allow for single base flipping from duplex DNA with minimal distortion of the double helix

Base flipping is a structural mechanism common to many DNA processing and repair enzymes. Changes in the local backbone torsions that occur during base flipping and the effect of environment on their behavior are of particular interest in understanding different base flipping mechanisms. In the present study, structures sampled during umbrella sampling molecular dynamics (MD) simulations of base flipping in aqueous and protein-bound environments, carried out with two different MD simulation strategies, are analyzed to find the most significant phosphodiester backbone distortions in the vicinity of the flipping base. Torsional sampling on the 5' side of the flipping base during flipping through the major groove shows similarities to the torsional sampling on the 3' side during flipping through the minor groove and vice versa. In differing environments, this behavior varies only marginally. These compensating torsional changes in the DNA backbone on 5' and 3' sides of the flipping base limit overall distortion of the DNA double helix during single base flipping. Rotameric intermediate states observed during base flipping are identified and postulated to be metastable states implicated in both large-scale structural changes and functional effects of chemical modifications in DNA.     

Sequence dependence in base flipping: experimental and computational studies

Base flipping is the movement of a DNA base from an intrahelical, base-stacked position to an extrahelical, solvent-exposed position. As there are favorable interactions for an intrahelical base, both hydrogen bonding and base stacking, base flipping is expected to be energetically prohibitive for an undamaged DNA duplex. For damaged DNA bases, however, the energetic cost of base flipping may be considerably lower. Using a selective, non-covalent assay for base flipping, the sequence dependence of base flipping in DNA sequences containing an abasic site has been studied. The dissociation constants of the zinc-cyclen complex to small molecules and single strands of DNA as well as the equilibrium constants for base flipping have been determined for these sequences. Molecular dynamics simulations of the zinc-cyclen complex bound to both single- and double-stranded DNA have been performed in an attempt to rationalize the differences in the dissociation constants obtained for the two systems. The results are compared to previous studies of base flipping in DNA containing an abasic site.     

Free energy simulations of uncatalyzed DNA replication fidelity: structure and stability of T.G and dTTP.G terminal DNA mismatches flanked by a single dangling nucleotide

A reference system for DNA replication fidelity was studied by free energy perturbation (FEP) and linear interaction energy (LIE) methods. The studied system included a hydrated duplex DNA with the 5'-CG dangling end of the templating strand, and dCTP4-.Mg2+ or dTTP4-.Mg2+ inserted opposite the dangling G to form a correct (i.e., Watson-Crick) or incorrect (i.e., wobble) base pair, respectively. The average distance between the 3'-terminal oxygen of the primer strand and the alpha-phosphorus of dNTP was found to be 0.2 A shorter for the correct base pair than for the incorrect base pair. Binding of the incorrect dNTP was found to be disfavored by 0.4 kcal/mol relative to the correct dNTP. We estimated that improved binding and more near-attack configurations sampled by the correct base pair should translate in aqueous solution and in the absence of DNA polymerase into a six times faster rate for the incorporation of the correct dNTP into DNA. The accuracy of the calculated binding free energy difference was verified by examining the relative free energy for melting duplex DNA containing GC and GT terminal base pairs flanked by a 5' dangling C. The calculated LIE and FEP free energies of 1.7 and 1.1 kcal/mol, respectively, compared favorably with the experimental estimate of 1.4 kcal/mol obtained using the nearest neighbor parameters. To decompose the calculated free energies into additive electrostatic and van der Waals contributions and to provide a set of rigorous theoretical data for the parametrization of the LIE method, we suggested a variant of the FEP approach, for which we coined a binding-relevant free energy (BRFE) acronym. BRFE approach is characterized by its unique perturbation pathway and by its exclusion of the intramolecular energy of a rigid part of the ligand from the total potential energy.     

6MAP, a fluorescent adenine analogue, is a probe of base flipping by DNA photolyase

Cyclobutylpyrimidine dimers (CPDs) are formed between adjacent pyrimidines in DNA when it absorbs ultraviolet light. CPDs can be directly repaired by DNA photolyase (PL) in the presence of visible light. How PL recognizes and binds its substrate is still not well understood. Fluorescent nucleic acid base analogues are powerful probes of DNA structure. We have used the fluorescent adenine analogue 6MAP, a pteridone, to probe the local double helical structure of the CPD substrate when bound by photolyase. Duplex melting temperatures were obtained by both UV-vis absorption and fluorescence spectroscopies to ascertain the effect of the probe and the CPD on DNA stability. Steady-state fluorescence measurements of 6MAP-containing single-stranded and doubled-stranded oligos with and without protein show that the local region around the CPD is significantly disrupted. 6MAP shows a different quenching pattern compared to 2-aminopurine, another important adenine analogue, although both probes show that the structure of the complementary strand opposing the 5'-side of the CPD lesion is more destacked than that opposing the 3'-side in substrate/protein complexes. We also show that 6MAP/CPD duplexes are substrates for PL. Vertical excitation energies and transition dipole moment directions for 6MAP were calculated using time-dependent density functional theory. Using these results, the F?rster resonance energy transfer efficiency between the individual adenine analogues and the oxidized flavin cofactor was calculated to account for the observed intensity pattern. These calculations suggest that energy transfer is highly efficient for the 6MAP probe and less so for the 2Ap probe. However, no experimental evidence for this process was observed in the steady-state emission spectra.     

Intermolecular Interactions in Photodamaged DNA from Density Functional Theory Symmetry‐Adapted Perturbation Theory

The intermolecular interactions of the photodamaged cyclobutane pyrimidine dimer (CPD) lesion with adjacent nucleobases in the native intrahelical DNA double strand are investigated at the level of density functional theory symmetry‐adapted perturbation theory (DFT‐SAPT) and compared to the original (or repaired) case with pyrimidines (TpT) instead of CPD. The CPD aggregation is on average destabilized by about 6 kcal mol^?1 relative to that involving TpT. The effect of destabilization is asymmetric, that is, it involves a single H‐bonding (Watson–Crick (WC) type) base‐pair interaction.     

Base pairing and replicative processing of the formamidopyrimidine-dG DNA lesion

The 2,6-diamino-4-hydroxy-5-formamidopyrimidine of 2'-deoxyguanosine (FaPydG) is one of the major DNA lesions found after oxidative stress in cells. To clarify the base pairing and coding potential of this major DNA lesion with the aim to estimate its mutagenic effect, we prepared oligonucleotides containing a cyclopentane based analogue of the DNA lesion (cFaPydG). In addition, oligonucleotides containing the cyclopentane analogue of 2'-deoxyguanosine (cdG), and oligonucleotides containing 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG) were synthesized. The thermodynamic stability of duplexes containing these building blocks and all canonical counterbases were determined by concentration dependent melting-point measurements (van't Hoff plots). The data reveal that cFaPydG greatly destabilizes a DNA duplex (DeltaDeltaG degrees (298K) approximately 2-4 kcal mol(-1)). The optimal base pairing partner for the cFaPydG lesion is dC. Investigation of duplexes containing dG and cdG shows that the effect of substituting the deoxyribose by a cyclopentane moiety is marginal. The data also provide strong evidence that the FaPydG lesion is unable to form a base pair with dA. Our computational studies indicate that the syn-conformation required for base pairing with dA is energetically unfavorable. This is in contrast to 8-oxodG for which the syn-conformation represents the energetic minimum. Kinetic primer extension studies using S. cerevisiae Pol eta reveal that cFaPydG is replicated in an error-free fashion. dC is inserted 2-3 orders of magnitude more efficiently than dT or dA, showing that FaPydG is a lesion which retains the coding potential of dG. This is also in contrast to 8-oxodG, for which base pairing with dC and dA was established.     

A theoretical study of the rotational isomers of nitrosomethanol by semiempirical (AM1) and ab initio methods

The conformational potential energy surface as a function of the two internal torsion angles in C-nitrosomethanol has been obtained using the semiempirical AM1 method. Optimized geometries are reported for the local minima on this surface and also for the corresponding points on the HF/6-31G, 6-31G*, and 6-31G** surfaces. All methods predict cis and trans minima which occur in degenerate pairs, each pair being connected by a transition state of C_s symmetry. The AM1 structures are found to compare well with the corresponding ab initio structures. Ab initio HF/6-31G and HF/6-31G* harmonic vibrational frequencies are reported for the cis and trans forms of nitrosomethanol. When scaled appropriately the calculated frequencies are found to compare well with experimental frequencies. The ab initio calculations predict the energy barrier for cis → trans isomerization to be between 5.8 and 6.5 kcal/mol with the trans → cis isomerization barrier lying between 2.3 and 6.5 kcal/mol. The corresponding AM1 energy barriers are around 1 kcal/mol lower in energy. The ab initio calculations predict the barrier to conversion between the two cis rotamers to be very small with the AM1 value being around 1 kcal/mol. Both AM1 and ab initio calculations predict interconversion between trans rotamers to require between 1.2 and 1.4 kcal/mol.     

On the molecular mechanism of drug intercalation into DNA: a simulation study of the intercalation pathway, free energy, and DNA structural changes

Intercalation into DNA (insertion between a pair of base pairs) is a critical step in the function of many anticancer drugs. Despite its importance, a detailed mechanistic understanding of this process at the molecular level is lacking. We have constructed, using extensive atomistic computer simulations and umbrella sampling techniques, a free energy landscape for the intercalation of the anticancer drug daunomycin into a twelve base pair B-DNA. A similar free energy landscape has been constructed for a probable intermediate DNA minor groove-bound state. These allow a molecular level understanding of aspects of the thermodynamics, DNA structural changes, and kinetic pathways of the intercalation process. Key DNA structural changes involve opening the future intercalation site base pairs toward the minor groove (positive roll), followed by an increase in the rise, accompanied by hydrogen bonding changes of the minor groove waters. The calculated intercalation free energy change is -12.3 kcal/mol, in reasonable agreement with the experimental estimate -9.4 kcal/mol. The results point to a mechanism in which the drug first binds to the minor groove and then intercalates into the DNA in an activated process, which is found to be in general agreement with experimental kinetic results.     

Consecutive GC base pairs determine the energy barrier of DNA duplex formation under molecularly crowded conditions

The kinetics of DNA duplex formation was affected by the addition of PEGs with different masses (MW = 200-8000) to an aqueous solution; for each condition, two duplexes (5'-TAGGTTATAA-3'/5'-TTATAACCTA-3' and 5'-CAGGTCACAG-3'/5'-CTGTGACCTG-3') with different stabilities were formed after overcoming the same association activation energy barrier, suggesting that the formation of consecutive GC base pairs in the helices rather than the helix terminus is the initiation nucleus for DNA duplex formation not only in the absence, but also in the presence of PEGs.     

Polarizable model potential function for nucleic acid bases

A polarizable model potential (PMP) function for adenine (A), cytosine (C), guanine (G), thymine (T), and uracil (U) is developed on the basis of ab initio molecular orbital calculations at the MP2/6-31+G* level. The PMP function consists of Coulomb, van der Waals, and polarization terms. The permanent atomic charges of the Coulomb term are determined by using electrostatic potential (ESP) optimization. The multicenter polarizabilities of the polarization term are determined by using polarized one-electron potential (POP) optimization in which the electron density changes induced by a test charge are target. Isotropic and anisotropic polarizabilities are adopted as the multicenter polarizabilities. In the PMP calculations using the optimized parameters, the interaction energies of Watson-Crick type A-T and C-G base pairs were -15.6 and -29.4 kcal/mol, respectively. The interaction energy of Hoogsteen type A-T base pair was -17.8 kcal/mol. These results reproduce well the quantum chemistry calculations at the MP2/6-311++G(3df,2pd) level within the differences of 0.6 kcal/mol. The stacking energies of A-T and C-G were -9.7 and -10.9 kcal/mol. These reproduce well the calculation results at the MP2/6-311++G (2d,2p) level within the differences of 1.3 kcal/mol. The potential energy surfaces of the system in which a sodium ion or a chloride ion is adjacent to the nucleic acid base are calculated. The interaction energies of the PMP function reproduced well the calculation results at the MP2/6-31+G* or MP2/6-311++G(2d,2p) level. The reason why the PMP function reproduces well the high-level quantum mechanical interaction energies is addressed from the viewpoint of each energy terms.     

Neomycin-neomycin dimer: an all-carbohydrate scaffold with high affinity for AT-rich DNA duplexes

A dimeric neomycin-neomycin conjugate 3 with a flexible linker, 2,2'-(ethylenedioxy)bis(ethylamine), has been synthesized and characterized. Dimer 3 can selectively bind to AT-rich DNA duplexes with high affinity. Biophysical studies have been performed between 3 and different nucleic acids with varying base composition and conformation by using ITC (isothermal calorimetry), CD (circular dichroism), FID (fluorescent intercalator displacement), and UV (ultraviolet) thermal denaturation experiments. A few conclusions can be drawn from this study: (1) FID assay with 3 and polynucleotides demonstrates the preference of 3 toward AT-rich sequences over GC-rich sequences. (2) FID assay and UV thermal denaturation experiments show that 3 has a higher affinity for the poly(dA)·poly(dT) DNA duplex than for the poly(dA)·2poly(dT) DNA triplex. Contrary to neomycin, 3 destabilizes poly(dA)·2poly(dT) triplex but stabilizes poly(dA)·poly(dT) duplex, suggesting the major groove as the binding site. (3) UV thermal denaturation studies and ITC experiments show that 3 stabilizes continuous AT-tract DNA better than DNA duplexes with alternating AT bases. (4) CD and FID titration studies show a DNA binding site size of 10-12 base pairs/drug, depending upon the structure/sequence of the duplex for AT-rich DNA duplexes. (5) FID and ITC titration between 3 and an intramolecular DNA duplex [d(5'-A(12)-x-T(12)-3'), x = hexaethylene glycol linker] results in a binding stoichiometry of 1:1 with a binding constant ～10(8) M(-1) at 100 mM KCl. (6) FID assay using 3 and 512 hairpin DNA sequences that vary in their AT base content and placement also show a higher binding selectivity of 3 toward continuous AT-rich than toward DNA duplexes with alternate AT base pairs. (7) Salt-dependent studies indicate the formation of three ion pairs during binding of the DNA duplex d[5'-A(12)-x-T(12)-3'] and 3. (8) ITC-derived binding constants between 3 and DNA duplexes have the following order: AT continuous, d[5'-G(3)A(5)T(5)C(3)-3'] > AT alternate, d[5'-G(3)(AT)(5)C(3)-3'] > GC-rich d[5'-A(3)G(5)C(5)T(3)-3']. (9) 3 binds to the AT-tract-containing DNA duplex (B* DNA, d[5'-G(3)A(5)T(5)C(3)-3']) with 1 order of magnitude higher affinity than to a DNA duplex with alternating AT base pairs (B DNA, d[5'-G(3)(AT)(5)C(3)-3']) and with almost 3 orders of magnitude higher affinity than a GC-rich DNA (A-form, d[5'-A(3)G(5)C(5)T(3)-3']).