Using liquid integral equation theory, we calculate the pair correlations of particles that interact via a smooth repulsive pair potential in d = 4 spatial dimensions. We discuss the performance of different closures for the Ornstein-Zernike equation, by comparing the results to computer simulation data. Our results are of relevance to understand crystal and glass formation in high-dimensional systems.