yDNA versus xDNA pyrimidine nucleobases: computational evidence for dependence of duplex stability on spacer location

The structural and binding properties of the natural and x- and y-pyrimidines were compared using computational methods. Our calculations show that although the x-pyrimidines favor different orientations about the glycosidic bond compared to the natural pyrimidines, which could have implications for the formation and resulting stability of xDNA duplexes and jeopardize the selectivity of expanded nucleobases, y-pyrimidines have rotational profiles more similar to the natural bases. Increasing the pyrimidine size using a benzene spacer leads to relatively minor changes in the hydrogen-bond strength of isolated Watson-Crick base pairs. However, differences in the anomeric carbon distances in pairs composed of x- or y-pyrimidines suggest yDNA may yield a more optimal expanded structure. By stacking two monomers via their centers of mass, we find that the expanded nucleobases stack much stronger than the natural bases. Additionally, although replacing xT by yT changes the stacking energy by less than 5 kJ mol (-1), replacing xC by yC significantly strengthens complexes with the natural nucleobases (by up to 30%). Calculations on larger duplex models composed of four nucleobases reveal that x- and y-pyrimidines can increase duplex stability of natural helices by strengthening both the intra and interstrand stacking interactions. Furthermore, when the total stability (sum of all hydrogen-bonding and (intrastrand and interstrand) stacking interactions) of the larger models is considered, y-pyrimidines lead to more stable complexes than x-pyrimidines for all but three duplex sequences. Thus, through analysis of a variety of properties, our calculations suggest that the location of the benzene spacer affects the properties of expanded nucleobases and the stability of expanded duplexes, and therefore should be carefully considered when designing future expanded analogues.