Computational Evaluation of Intermolecular Interactions of a Universal Base 3-Nitropyrrole in Stacked Dimers and DNA Duplexes

Strict Interactions of Fifth Letters, Hydrophobic Unnatural Bases, in XenoAptamers with Target Proteins

XenoAptamers are DNA fragments containing additional letters (unnatural bases, UBs) that bind specifically to their target proteins with high affinities (sub-nanomolar KD values). One of the UBs is the highly hydrophobic 7-(2-thienyl)imidazo[4,5-b]pyridine (Ds), which significantly increases XenoAptamers' affinities to targets. Originally, Ds was developed as a third base pair with a complementary UB, 2-nitro-4-propynylpyrrole (Px), for replication, and thus it can be used for aptamer generation by an evolutional engineering method involving PCR amplification. However, it is unclear whether the Ds base is the best component as the hydrophobic fifth-letter ligand for interactions with target proteins. To optimize the ligand structure of the fifth letter, we prepared 13 Ds variants and examined the affinities of XenoAptamers containing these variants to target proteins. The results obtained using four XenoAptamers prepared by the replacement of Ds bases with variants indicated that subtle changes in the chemical structure of Ds significantly affect the XenoAptamer affinities. Among the variants, placing either 4-(2-thienyl)pyrrolo[2,3-b]pyridine (Ys) or 4-(2-thienyl)benzimidazole (Bs) at specific Ds positions in each original XenoAptamer greatly improved their affinities to targets. The Ys and Bs bases are variants derived by replacing only one nitrogen with a carbon in the Ds base. These results demonstrate the strict intramolecular interactions, which are not simple hydrophobic contacts between UBs and targets, thus providing a method to mature XenoAptamers' affinities to targets.

A computational proposal for the experimentally observed discriminatory behavior of hypoxanthine, a weak universal nucleobase

A computational model composed of six nucleobases was used to investigate why hypoxanthine does not yield duplexes of equal stability when paired opposite each of the natural DNA nucleobases. The magnitudes of all nearest-neighbor interactions in a DNA helix were calculated, including hydrogen-bonding, intra- and interstrand stacking interactions, as well as 1-3 intrastrand stacking interactions. Although the stacking interactions in DNA relevant arrangements are significant and account for at least one third of the total stabilization energy in our nucleobase complexes, the trends in the magnitude of the stacking interactions cannot explain the relative experimental melting temperatures previously reported in the literature. Furthermore, although the total hydrogen-bonding interactions explain why hypoxanthine preferentially pairs with cytosine, the experimental trend for the remaining nucleobases (A, T, G) is not explained. In fact, the calculated pairing preference of hypoxanthine matches that determined experimentally only when the sum of all types of nearest-neighbor interactions is considered. This finding highlights a strong correlation between the relative magnitude of the total nucleobase-nucleobase interactions and measured melting temperatures for DNA strands containing hypoxanthine despite the potential role of other factors (including hydration, temperature, sugar-phosphate backbone). By considering a large range of sequence combinations, we reveal that the binding preference of hypoxanthine is strongly dependent on the nucleobase sequence, which may explain the varied ability of hypoxanthine to universally bind to the natural bases. As a result, we propose that future work should closely examine the interplay between the dominant nucleobase-nucleobase interactions and the overall strand stability to fully understand how sequence context affects the universal binding properties of modified bases and to aid the design of new molecules with ambiguous pairing properties.

Comparison of intrinsic stacking energies of ten unique dinucleotide steps in A-RNA and B-DNA duplexes. Can we determine correct order of stability by quantum-chemical calculations?

High level ab initio methods have been used to study stacking interactions in ten unique base pair steps both in A-RNA and in B-DNA duplexes. The protocol for selection of geometries based on molecular dynamics (MD) simulations is proposed, and its suitability is demonstrated by comparison with stacking in steps at fiber diffraction geometries. It is shown that fiber diffraction geometries are not sufficiently accurate for interaction energy calculations. In addition, the protocol for selection of geometries based on MD simulations allows for the evaluation of the variability of the intrinsic stacking energies along the MD trajectories. The uncertainty in stacking energies (difference between the most and least stable geometry) due to the dynamical nature of systems can be, in some cases, as large as 3.0 kcal x mol(-1), which is almost 50% of the actual sequence dependence of base stacking energies (the energy difference between the most and least stable sequences). Thus, assessing the relative magnitude of the gas phase stacking energy using a single geometry for each sequence is insufficient to obtain an unambiguous order of gas phase stacking energies in canonical double helices. Though the ordering of ten unique dinucleotide steps cannot be definitive, some general conclusions were drawn. The stacking energies of base pair steps in A-RNA are more evenly separated compared to B-DNA, and their ordering is less sensitive to the dynamics of the system compared to be B-DNA. The most stable step both in B-DNA and A-RNA is the GC/GC [corrected] step that is well separated from the second most stable step CG/CG. [corrected] Also the least stable step (the CC/GG step) is well separated from the rest of the structures. The calculations further show that B-DNA stacking is favorable only marginally (on average by 1.14 kcal x mol(-1) per base pair step) over A-RNA stacking, and this difference vanishes after subtracting the stabilizing van der Waals effect of the thymine 5-methyl group that is absent in RNA. Basically, no correlation between the sequence dependence of gas phase stacking energies and the sequence dependence of DeltaG degrees(37) free energies used in nearest-neighbor models was found either for B-DNA or for A-RNA. This reflects the complexity of the balance of forces that are responsible for the sequence dependence of thermodynamics stability of nucleic acids, which masks the effect of the intrinsic interactions between the stacked base pairs.

Nature and magnitude of aromatic stacking of nucleic acid bases

This review summarises recent advances in quantum chemical calculations of base-stacking forces in nucleic acids. We explain in detail the very complex relationship between the gas-phase base-stacking energies, as revealed by quantum chemical (QM) calculations, and the highly variable roles of these interactions in nucleic acids. This issue is rarely discussed in quantum chemical and physical chemistry literature. We further extensively discuss methods that are available for base-stacking studies, complexity of comparison of stacking calculations with gas phase experiments, balance of forces in stacked complexes of nucleic acid bases, and the relation between QM and force field descriptions. We also review all recent calculations on base-stacking systems, including details analysis of the B-DNA stacking. Specific attention is paid to the highest accuracy QM calculations, to the decomposition of the interactions, and development of dispersion-balanced DFT methods. Future prospects of computational studies of base stacking are discussed.

Physical nature of ethidium and proflavine interactions with nucleic acid bases in the intercalation plane

On the basis of the crystallographic structures of three nucleic acid intercalation complexes involving ethidium and proflavine, we have analyzed the interaction energies between intercalator chromophores and their four nearest bases, using a hybrid variation-perturbation method at the second-order Møller-Plesset theory level (MP2) with a 6-31G(d,p) basis set. A total MP2 interaction energy minimum precisely reproduces the crystallographic position of the ethidium chromophore in the intercalation plane between UA/AU bases. The electrostatic component constitutes the same fraction of the total energy for all three studied structures. The multipole electrostatic interaction energy, calculated from cumulative atomic multipole moments (CAMMs), was found to converge only after including components above the fifth order. CAMM interaction surfaces, calculated on grids in the intercalation planes of these structures, reasonably reproduce the alignment of intercalators in crystal structures; they exhibit additional minima in the direction of the DNA grooves, however, which also need to be examined at higher theory levels if no crystallographic data are given.

How Do Size-Expanded DNA Nucleobases Enhance Duplex Stability? Computational Analysis of the Hydrogen-Bonding and Stacking Ability of xDNA Bases

Computational chemistry (B3LYP, MP2) is used to study the properties of size-expanded DNA nucleobases generated by inserting a benzene spacer into the natural nucleobases. Although the addition of the spacer does not significantly affect the hydrogen-bonding properties of natural nucleobases, the orientation of the base about the glycosidic bond necessary for Watson-Crick binding is destabilized, which could have implications for the selectivity of expanded bases, as well as the stability of expanded duplexes. Consideration of the (stacked) binding energies in the preferred relative orientation of natural and expanded nucleobases aligned according to their centers of mass reveals that the stacking within natural dimers can be increased by up to 50% upon expansion of one nucleobase and up to 90% upon expansion of two nucleobases. The implications of these findings to the stability of expanded duplexes were revealed by considering simplified models of natural and mixed duplexes composed of four nucleobases. Although intra- and interstrand interactions within double helices are typically less than those predicted when nucleobases are stacked according to their centers of mass, some nucleobases utilize their full stacking potential within double helices, where both intra- and interstrand interactions can be significant. Most importantly, increasing the size of nucleobases within the duplex significantly increases both intra- and interstrand stacking interactions. Specifically, some interactions are double the magnitude of the corresponding intrastrand interactions in natural helices, and even greater increases in interstrand interactions are sometimes found. Thus, our work suggests that mixed duplexes composed of natural bases hydrogen bound to expanded bases may exploit the increase in the inherent stacking ability of the expanded bases in more than one way and thereby afford duplexes with greater stability than natural DNA.

Environmental Effects on the Enhancement in Natural and Damaged DNA Nucleobase Acidity Because of Discrete Hydrogen-Bonding Interactions

The present study uses density functional theory to carefully consider the effects of the environment on the enhancement in (natural and damaged) DNA nucleobase acidities because of multiple hydrogen-bonding interactions. Although interactions with one small molecule can increase the acidity of the nucleobases by up to 60 kJ mol-1 in the gas phase, the maximum increase in enzymatic-like environments is expected to be approximately 40 kJ mol-1, which reduces to approximately 30 kJ mol-1 in water. Furthermore, the calculated (simultaneous) effects of two, three, or four molecules are increasingly less than the sum of the individual (additive) effects with an increase in the number and acidity of the small molecules bound or the dielectric constant of the solvent. Regardless of these trends, our calculations reveal that additional hydrogen-bonding interactions will have a significant effect on nucleobase acidity in a variety of environments, where the exact magnitude of the effect depends on the properties of the small molecule bound, the nucleobase binding site, and the solvent. The maximum increase in nucleobase acidity because of interactions with up to four small molecules is approximately 80 kJ mol-1 in enzymatic-like environments (or 65 kJ mol-1 in water). These results suggest that hydrogen-bonding interactions likely play an important role in many biological processes by changing the physical and chemical properties of the nucleobases.

Effects of Hydrogen-Bonding and Stacking Interactions with Amino Acids on the Acidity of Uracil

The effects of hydrogen-bonding interactions with amino acids on the (N1) acidity of uracil are evaluated using (B3LYP) density functional theory. Many different binding arrangements of each amino acid to three uracil binding sites are considered. The effects on the uracil acidity are found to significantly depend upon the nature of the amino acid and the binding orientation, but weakly depend on the binding site. Our results reveal that in some instances small models for the amino acids can be used, while for other amino acids larger models are required to properly describe the binding to uracil. The gas-phase acidity of uracil is found to increase by up to approximately 60 kJ mol(-1) due to discrete hydrogen-bonding interactions. Although (MP2) stacking interactions with aromatic amino acids decrease the acidity of uracil, unexpected increases in the acidity are found when any of the aromatic amino acids, or the backbone, hydrogen bond to uracil. Consideration of enzymatic and aqueous environments leads to decreases in the effects of the amino acids on the acidity of uracil. However, we find that the magnitude of the decrease varies with the nature of the molecule bound, as well as the (gas-phase) binding orientations and strengths, and therefore solvation effects should be considered on a case-by-case basis in future work. Nevertheless, the effects of amino acid interactions within enzymatic environments are as much as approximately 35 kJ mol(-1). The present study has general implications for understanding the nature of active site amino acids in enzymes, such as DNA repair enzymes, that catalyze reactions involving anionic nucleobase intermediates.

Synthesis of 2'-O-methyl-RNAs incorporating a 3-deazaguanine, and UV melting and computational studies on its hybridization properties

2′-O-Methyl-RNAs incorporating 3-deazaguanine (c³G) were synthesized by use of N,N-diphenylcarbamoyl and N,N-dimethylaminomethylene as its base protecting groups to suppress sheared-type 5′-GA-3′/5′-GA-3′ tandem mismatched base pairing which requires the N³ atom. These modified RNAs hybridized more weakly with the complementary and single mismatch-containing RNAs than the unmodified RNAs. The T_m experiments were performed to clarify the effects of replacement of the fifth G with c³G on stabilization of 2′-O-methyl-(5′-CGGCGAGGAG-3′)/5′-CUCCGAGCCG-3′ and 2′-O-methyl-(5′-CGGGGACGAG-3′)/5′-CUCGGACCCG-3′duplexes, which form sheared-type and face-to-face type 5′-GA-3′/5′-GA-3′ tandem mismatched base pairs, respectively. Consequently, this replacement led to more pronounced destabilization of the former duplex that needs the N³ atom for the sheared-type base pair than the latter that does not need it for the face-to-face type base pair. A similar tendency was observed for 2′-O-methyl-RNA/DNA duplexes. These results suggest that the N³ atom of G plays an important role in stabilization of the canonical G/C base pair as well as the base discrimination and its loss suppressed formation of the undesired sheared-type mismatched base pair. Computational studies based on ab initio calculations suggest that the weaker hydrogen bonding ability and larger dipole moment of c³G can be the origin of the lower T_m.