Theoretical study of the geometrical arrangement of GT and GA wobble pairs in two short duplexes, Proton NMR chemical shifts and interaction energy calculations

NMR shielding constants are calculated for the base protons of duplexes formed by the dodecamer d(CGTGAATTCGCG) and the decamer d(CCAAGATTGG). A good agreement with experimental data is obtained for B-DNA helices in which the wobble GT and GA pairs are in the plane of the corresponding GC pairs of the parent duplexes formed by d(CGCGAATTCGCG) and d(CCAAGCTTGG), if the glycosyl bonds of T and G or A and G are symmetrical with respect to the dyad axis of the Watson-Crick GC pair. Interaction energy calculations show that this type of geometrical arrangement, which implies a distortion of the ribonphosphate backbone of both strands of the duplexes are more stable than those in which only one strand has its conformation modified by the presence of the wobble pair. For the duplex containing the GA pair, NMR chemical shifts as well as interaction energy computations favour the Watson-Crick hydrogen bonding scheme. The variation of the different contributions (intrastrand, interstrand, pair-pair) to the interaction energy between the bases of the duplexes, with the geometrical arrangement of the wobble pairs, is reported.