The vibrational spectral density is an important physical parameter needed to describe both linear and non-linear spectra of multi-chromophore systems such as photosynthetic complexes. Low-temperature techniques such as hole burning (HB) and fluorescence line narrowing are commonly used to extract the spectral density for a given electronic transition from experimental data. We report here that the lineshape function formula reported by Hayes et al. [J. Phys. Chem. 98, 7337 (1994)] in the mean-phonon approximation and frequently applied to analyzing HB data contains inconsistencies in notation, leading to essentially incorrect expressions in cases of moderate and strong electron-phonon (el-ph) coupling strengths. A corrected lineshape function L(ω) is given that retains the computational and intuitive advantages of the expression of Hayes et al. [J. Phys. Chem. 98, 7337 (1994)]. Although the corrected lineshape function could be used in modeling studies of various optical spectra, we suggest that it is better to calculate the lineshape function numerically, without introducing the mean-phonon approximation. New theoretical fits of the P870 and P960 absorption bands and frequency-dependent resonant HB spectra of Rb. sphaeroides and Rps. viridis reaction centers are provided as examples to demonstrate the importance of correct lineshape expressions. Comparison with the previously determined el-ph coupling parameters [Johnson et al., J. Phys. Chem. 94, 5849 (1990); Lyle et al., ibid. 97, 6924 (1993); Reddy et al., ibid. 97, 6934 (1993)] is also provided. The new fits lead to modified el-ph coupling strengths and different frequencies of the special pair marker mode, ωsp, for Rb. sphaeroides that could be used in the future for more advanced calculations of absorption and HB spectra obtained for various bacterial reaction centers.