Structure and dynamics of double helical DNA in torsion angle hyperspace: a molecular mechanics approach

Abstract

Analysis of the conformational space populated by the torsion angles and the correlation between the conformational energy and the sequence of DNA are important for fully understanding DNA structure and function. Presence of seven variable torsion angles about single covalent bonds in DNA main chain puts a big challenge for such analysis. We have carried out restrained energy minimization studies for four representative dinucleosides, namely d(ApA):d(TpT), d(CpG):d(CpG), d(GpC):d(GpC) and d(CpA):d(TpG) to determine the energy hyperspace of DNA in context to the values of the torsion angles and the structural properties of the DNA conformations populating the favorable regions of this energy hyperspace. The torsion angles were manipulated by constraining their values at the reference points and then performing energy minimization. The energy minima obtained on the potential energy contour plots mostly correspond to the conformations populated in crystal structures of DNA. Some novel favorable conformations that are not present in crystal structure data are also found. The plots also suggest few low energy routes for conformational transitions or the associated energy barrier heights. Analyses of base pairing and stacking possibility reveal structural changes accompanying these transitions as well as the flexibility of different base steps towards variations in different torsion angles.