Single nucleotide mutation rates have critical implications for human evolution and genetic diseases. Accurate modeling of these mutation rates has long remained an open problem since the rates vary substantially across the human genome. A recent model, however, explained much of the variation by considering higher order nucleotide interactions in the local (7-mer) sequence context around mutated nucleotides. Despite this model’s predictive value, we still lack a clear understanding of the biophysical mechanisms underlying the variations in genome-wide mutation rates. DNA shape features are geometric measurements of DNA structural properties, such as helical twist and tilt, and are known to capture information on interactions between neighboring nucleotides within a local context. Motivated by this characteristic of DNA shape features, we used them to model mutation rates in the human genome. These DNA shape feature based models improved both the accuracy (up to 14%) and the interpretability over the current nucleotide sequence-based models. The models also discovered the specific shape features that capture the most variability in mutation rates, and distinguished between the most and the least mutated sequence contexts, thus characterizing mutation promoting properties of the genomic DNA. To our knowledge, this is the first attempt that demonstrates the structural underpinnings of nucleotide mutations in the human genome and lays the groundwork for future studies to incorporate DNA shape information in modeling genetic variations.