We apply the integral equation theory based on the distributed partial wave expansion to molecular liquids. Improvements in the present study consist of: (1) an orientational-dependent closure relation using the two-dimensional Lebedev quadrature, (2) small- k expansions of distributed partial functions and their projections, and (3) a renormalized formula for charged molecular systems. We examine the method in the calculations of tetrahedral model and non-polar/polar diatomic fluids.