We present an improved simulation methodology to describe non-photochemical hole-burned (NPHB) spectra. The model, which includes both frequency-dependent excitation energy transfer (EET) rate distributions and burning following EET, provides reasonable fits of various optical spectra including resonant and nonresonant holes in the case of FMO complex. A qualitative description of the NPHB process in light of a very complex protein energy landscape is briefly discussed. As an example, we show that both resonant and nonresonant HB spectra obtained for the 825 nm band of the trimeric FMO of C. tepidum are consistent with the presence of a relatively slow EET between the lowest energy states of the monomers of the trimer (mostly localized on BChl a 3), with a weak (similar to 1 cm(-1)) coupling between these states revealed via calculated emission spectra. We argue that the nature of the so-called 825 nm absorption band of the FMO trimer, contrary to the presently accepted Consensus, cannot be explained by a single transition.