Rattlesnake venom can differ in composition and in metalloproteinase-associated activities. The molecular basis for this intra-species variation in Crotalus scutulatus scutulatus (Mojave rattlesnake) remains an enigma. To understand the molecular basis for intra-species variation of metalloproteinase-associated activities, we modeled the threedimensional structures of four metalloproteinases based on the amino acid sequence of four variations of the proteinase domain of the C. s. scutulatus metalloproteinase gene (GP1, GP2, GP3, and GP4). For comparative purposes, we modeled the atrolysin metalloproteinases of C. atrox as well. All molecular models shared the same topology. While the atrolysin metalloproteinase molecular models contained highly conserved substrate binding sites, the Mojave rattlesnake metalloproteinases showed higher structural divergence when superimposed onto each other. The highest structural divergence among the four C. s. scutulatus molecular models was located at the northern cleft wall and the S' 1 -pocket of the substrate binding site, molecular regions that modulate substrate selectivity. Molecular dynamics and field potential maps for each C. s. scutulatus metalloproteinase model demonstrated that the non-hemorrhagic metalloproteinases (GP2 and GP3) contain highly basic molecular and field potential surfaces while the hemorrhagic metalloproteinases GP1 and atrolysin C showed extensive acidic field potential maps and shallow but less dynamic active site pockets. Hence, differences in the spatial arrangement of the northern cleft wall, the S' 1 -pocket, and the physico-chemical environment surrounding the catalytic site contribute to differences in metalloproteinase activities in the Mojave rattlesnake. Our results provide a structural basis for variation of metalloproteinase-associated activities in the rattlesnake venom of the Mojave rattlesnake.