To account for the relative loop motion, we also defined a hydrogen bond parameter between the phosphate and amide between the bases G (residue) and A (residue) at 5 and 8 positions in the loop structure (Fig. ). The switching function used here ensures that the calculated hydrogen bond values are continuous derivatives.