Effect of loop sequence on unzipping of short DNA hairpins

The dependence of stability on the sequence of a DNA hairpin has been investigated through atomistic simulations. For this, a sequence of 16 bases of a hairpin, which consists of a loop of four bases and a stem of six base pairs, has been considered. We have taken eight different sequences, where the first five base pairs were kept fixed in all sequences, whereas the loop sequence and the identity of the duplex base pair closing the loop have been varied. For these hairpin structures, force-induced melting (unzipping) studies were carried out to investigate the effect of the variables on the stability of hairpin. The temperature at which half of the base pairs are open is termed the melting temperature. We defined the unzipping force F_{h} (half of the base pairs are open) and showed that it may not provide the effect of closing the base pair or loop sequence on the stability of the DNA hairpin. In order to have a better understanding of the stability of a DNA hairpin, the closing base pair or hairpin loop must be open. This requires complete opening of the stem. We defined a force F_{c} at which all base pairs of the stem are open, and we showed that the F_{c} gives better understanding of DNA hairpin stability.