We report on the results obtained in the study of the parton saturation effects, taken into account through the multiscattering Glauber-Mueller approach, applied to the Drell-Yan (DY) process described in the color dipole picture. As a main result, one shows that those effects play an important role in the estimates of the DY differential cross section at RHIC energies.The Drell-Yan (DY) process, i.e. the production of massive lepton pairs in hadronic collisions, in conjunction with deep inelastic scattering (DIS), has been one of the important processes probing strong interaction physics. Recently, in connection with the availability of high energy accelerators, a great attention has been focused on the small-x region of QCD. There, the parton densities become high and the limits of the perturbative methods are tested. Such a region presents the onset of the saturation phenomenon (at scale Q 2 s ), i.e. the taming of the parton (mostly gluon) distribution due to nonlinear dynamics associated with unitarity effects 1 . In this contribution, we study the high energy DY cross section in the target rest frame. In this case, the relevant degrees of freedom are the projectile wavefunction and the dipole-proton effective cross section. The underlying process is the scattering of a parton from the projectile structure function off the target color field. This parton radiates (bremsstrahlung) a massive photon, which subsequently decays into a lepton pair. The interaction with the target can occur before or after the photon emission. A remarkable feature emerging is that the γ * q-proton (or hadron) interaction can be described by the same qq (color dipole)-proton cross section as in DIS. Although diagramatically no dipole is present, the interference among graphs results in a product of two quark amplitudes in the DY cross section testing the external gluonic field at two different transverse positions (impact parameters), in a similar way to DIS 2 . In such a representation, the photoabsorbtion cross section on deep inelastic is described by the convolution of the wavefunctions, Ψ γ * , from the virtual photon and the interaction dipole cross section, σ dip . The wavefunctions are considered taking into account the simplest photon Fock state configuration,