We develop an integral equation theory based on the partial wave expansion of the molecular Ornstein-Zernike (OZ) equation. The theory provides a rigorous and transparent framework for multiple site treatments of molecular fluids. We examine free-energy functional and closure expressions with pilot calculations of homonuclear diatomic Lennard-Jones liquids.