In this paper we report on the design and analyis of a multilevel method for the Ornstein-Zernike equations and related systems of integro-algebraic equations from computational chemistry. Our approach is based on an extension of the Atkinson-Brakhage method, with Newton-GMRES used as the coarse mesh solver. We report on several numerical experiments to illustrate the effectiveness of the method.