Advances in deep learning and AlphaFold2 have enabled the large-scale prediction of protein structures across species, opening avenues for studying protein function and evolution. Here we analyse 11,269 predicted and experimentally determined enzyme structures that catalyse 361 metabolic reactions across 225 pathways to investigate metabolic evolution over 400 million years in the Saccharomycotina subphylum. By linking sequence divergence in structurally conserved regions to a variety of metabolic properties of the enzymes, we reveal that metabolism shapes structural evolution across multiple scales, from species-wide metabolic specialization to network organization and the molecular properties of the enzymes. Although positively selected residues are distributed across various structural elements, enzyme evolution is constrained by reaction mechanisms, interactions with metal ions and inhibitors, metabolic flux variability and biosynthetic cost. Our findings uncover hierarchical patterns of structural evolution, in which structural context dictates amino acid substitution rates, with surface residues evolving most rapidly and small-molecule-binding sites evolving under selective constraints without cost optimization. By integrating structural biology with evolutionary genomics, we establish a model in which enzyme evolution is intrinsically governed by catalytic function and shaped by metabolic niche, network architecture, cost and molecular interactions.