An X- (9.7 GHz and W-band (94 GHz) pulsed electron-nuclear double resonance (ENDOR) study of the flavin cofactor of Escherichia coli DNA photolyase in its neutral radical form is presented. Through proton and deuteron ENDOR measurements at T = 80 K, we detect and characterize the full anisotropy of the hyperfine coupling (hfc) tensor of the proton or deuteron bound to N(5) of the isoalloxazine ring. Scaling of the anisotropic proton hfc components by multiplication with the quotient of the magnetogyric ratio of a deuteron and a proton, chiD/chiH, reveals subtle differences compared to the respective deuteron couplings obtained by 95-GHz deuterium ENDOR spectroscopy on an H-->D buffer-exchanged sample. These differences can be attributed to the different lengths of N(5)-H and N(5)-D bonds arising from the different masses of protons and deuterons. From the R(-3) dependence of the dipolar hyperfine splitting, we estimated that the N(5)-D bond is about 2.5% shorter than the respective N(5)-H bond. That such subtle bond-length differences can be resolved by pulsed ENDOR spectroscopy suggests that this method may be favorably used to probe the geometry of hydrogen bonds between the H(5) of the paramagnetic flavin and the protein backbone. Such information is only obtained with difficulty by other types of spectroscopy.