Lineshape simulations are presented for the multiple, overlapping X-band electron paramagnetic resonance (EPR) spectra in two non-heme, high-spin iron proteins: phenylalanine hydroxylase (PAH) and diferric transferrin. The aim of the calculations is to determine the fraction of iron contributing to each of the sites visible by EPR. The simulations are limited to the experimentally accessible transitions occurring at g-values greater than 1.7. In both PAH and transferrin, at least one of the iron sites is characterized by the ratio of zero-field splitting parameters, E/D, near 1/3 and a broad, asymmetric lineshape. A distribution in E/D-values is used in the simulations to account for this breadth and asymmetry. To test the E/D-distribution model, experimental X-band spectra of diferric transferrin at several salt concentrations are fit by simulation. In this test, first the low-field features arising from transitions between the lowest Kramers doublet levels are simulated using E/D-distributions for two sites. Second, parameters that provide a good fit for the lowest doublet transitions are shown also to fit the resonance near an effective g-value of 4.3 from the middle Kramers doublet transition. When applied to spectra of PAH in the resting state, the E/D-distribution approach accounts for the intensity of one of the two major species of iron. The other species is characterized by E/D = 0.032, and the spectrum of this portion of the resting enzyme may be simulated using a frequency-swept Gaussian lineshape. Spectra for the enzyme in an inhibitor-saturated state are also simulated. The simulations are consistent with previous biochemical studies that indicate that only the E/D = 0.032 form of iron participates in catalysis.