RADACK, a stochastic simulation of hydroxyl radical attack to DNA.

RADACK was conceived to simulate the radiation-induced attack to different DNA forms and complexes. It allows to separately calculate the probability of attack to each reactive atom of the sugar and of the base and takes into account the sequence-dependent structure of DNA as known from crystallographic or NMR studies or resulting from molecular modelling. The calculations are aimed to assess sequence-, structure- and ligand-dependent modulation of damages of sugar and bases, leading to single strand breaks (frank strand breaks, FSB) and alkali-labile base modifications (alkali-revealed breaks, ARB), respectively. The modelling procedure and the results of simulations for some representative structures (B, Z and quadruplex forms) are here described and discussed. The calculated relative probabilities of OH* radical attack to all reaction sites are compared to experimental FSB and ARB values. By a fitting procedure, the relative efficiencies of conversion of the C4' and C5'-centred radicals into FSB, epsilon (C4'): epsilon (C5'), and the relative efficiencies of base radicals- to- ARB conversion, epsilon(T) : epsilon(A) : epsilon(C) : epsilon(G), are then deduced for each DNA form. The ability of the model to account for the distribution of damages in DNA-ligand complexes is proven by its successful application to two DNA-protein systems : the lac repressor-lac operator complex and the nuclcosome core.