'Heavy' (isotopically-labeled) enzyme isotope effects offer a direct experimental probe of the role of protein vibrations on enzyme-catalyzed reactions. Here we have developed a strategy to generate isotopologues of the flavoenzyme pentaerythritol tetranitrate reductase (PETNR) where the protein and/or intrinsic flavin mononucleotide (FMN) cofactor are isotopically labeled with 2 H, 15 N and 13 C. Both the protein and cofactor contribute to the enzyme isotope effect on the reductive hydride transfer reaction, but their contributions are not additive and may partially cancel each other out. However, the isotope effect specifically arising from the FMN suggests that vibrations local to the active site play a role in the hydride transfer chemistry, while the protein-only 'heavy enzyme' effect demonstrates that protein vibrations contribute to catalysis in PETNR. In all cases, enthalpy-entropy compensation plays a major role in minimizing the magnitude of 'heavy enzyme' isotope effects. Fluorescence lifetime measurements of the intrinsic flavin mononucleotide show marked differences between 'light' and 'heavy' enzymes on the ns-ps timescale, suggesting relevant timescale(s) for those vibrations implicated in the 'heavy enzyme' isotope effect on the PETNR reaction.