The circular dichroism (CD) spectrum of a large biological molecule represented as an array of interacting chromophores is investigated. The configuration-averaged Green function formalism is developed to describe the CD and absorption spectra. The perturbation theory expansion is derived for absorption and CD spectra in the case of strong interaction of chromophores with their environment (solvent and/or internal dynamics) compared to their interaction with each other. We apply this formalism to study CD spectra of DNA hairpins.