Bending of DNA by asymmetric charge neutralization: all-atom energy simulations

DNA dodecamers of the alternating d(CG).d(CG) sequence with six phosphate groups either charge-neutralized or substituted by neutral methylphosphonates across the major or minor groove have been subjected to energy minimization to determine the conformational effect of the asymmetric elimination of phosphate charge. We report bending angles, directions of bending, and detailed structural characteristics such as groove widths and local base-pair parameters. Our principal results are that charge neutralization on one face of the DNA induces significant bending toward the neutralized face, in agreement with theoretical predictions on a simplified model and experimental data on a similar base-pair sequence, and that the DNA conformation averaged over all stereospecific methylphosphonate substitutions is nearly the same as the conformation produced by charge neutralization of the phosphates. Individual isomers, however, cover a wide range of structures, with the magnitude and direction of overall bending sensitive to the precise stereochemical pattern of neutralization. Our simulation does not explicitly contain counterions, and the results therefore suggest that counterions can influence DNA structure by neutralizing the phosphate charge. These data provide new hints into the molecular mechanisms which underlie the deformations of DNA structure induced by the binding of positively charged proteins and other tightly associated cationic species.