This article describes paramagnetic relaxation enhancement (PRE) in systems with high electron spin, S, where there is molecular interaction between a paramagnetic ion and a ligand outside of the first coordination sphere. The new feature of our treatment is an improved handling of the electron-spin relaxation, making use of the Redfield theory. Following a common approach, a well-defined second coordination sphere is assumed, and the PRE contribution from these more distant and shorter-lived ligands is treated in a way similar to that used for the first coordination sphere. This model is called "ordered second sphere," OSS. In addition, we develop here a formalism similar to that of Hwang and Freed [J. Chem. Phys. 63, 4017 (1975)], but accounting for the electron-spin relaxation effects. We denote this formalism "diffuse second sphere," DSS. The description of the dynamics of the intermolecular dipole-dipole interaction is based on the Smoluchowski equation, with a potential of mean force related to the radial distribution function. We have used a finite-difference method to calculate numerically a correlation function for translational motion, taking into account the intermolecular forces leading to an arbitrary radial distribution of the ligand protons. The OSS and DSS models, including the Redfield description of the electron-spin relaxation, were used to interpret the PRE in an aqueous solution of a slowly rotating gadolinium (III) complex (S=7/2) bound to a protein.