Predicting mutation frequencies in stem-loop structures of derepressed genes: implications for evolution.

This work provides evidence that, during transcription, the mutability (propensity to mutate) of a base in a DNA secondary structure depends both on the stability of the structure and on the extent to which the base is unpaired. Zuker's DNA folding computer program reveals the most probable stem-loop structures (SLSs) and negative energies of folding (-DeltaG) for any given nucleotide sequence. We developed an interfacing program that calculates (i) the percentage of folds in which each base is unpaired during transcription; and (ii) the mutability index (MI) for each base, expressed as an absolute value and defined as -follows: MI = (% total folds in which the base is unpaired) x (highest -DeltaG of all folds in which it is unpaired). Thus, MIs predict the relative mutation or reversion frequencies of unpaired bases in SLSs. MIs for 16 mutable bases in auxotrophs, selected during starvation in derepressed genes, are compared with 70 background mutations in lacI and ebgR that were not derepressed during mutant selection. All the results are consistent with the location of known mutable bases in SLSs. Specific conclusions are: (i) Of 16 mutable bases in transcribing genes, 87% have higher MIs than the average base of the sequence analysed, compared with 50% for the 70 background mutations. (ii) In 15 of the mutable bases of transcribing genes, the correlation between MIs and relative mutation frequencies determined experimentally is good. There is no correlation for 35 mutable bases in the lacI gene. (iii) In derepressed auxotrophs, 100% of the codons containing the mutable bases are within one codon's length of a stem, compared with 53% for the background mutable bases in lacI. (iv) The data suggest that environmental stressors may cause as well as select mutations in derepressed genes. The implications of these results for evolution are discussed.