The malleability of intrinsically disordered proteins (IDPs) has generated great interest in understanding how their conformations respond to crowded cellular environments. Experiments can report gross properties such as fluorescence resonance energy transfer (FRET) efficiency but cannot resolve the conformational ensembles of IDPs and their interactions with macromolecular crowders. Computation can in principle provide the latter information but in practice has been hampered by the enormous expense for realistic modeling of IDPs and crowders and for sufficient conformational sampling. Here, taking advantage of a powerful method called FMAP (fast Fourier transform-based modeling of atomistic protein-crowder interactions), we computed how the conformational ensembles of three IDPs are modified in concentrated polyethylene glycol (PEG) 6000 solutions. We represented the IDPs at the all-atom level and the PEG molecules at a coarse-grained level and calculated the experimental observable, i.e., FRET efficiency. Whereas accounting for only steric repulsion of PEG led to overestimation of crowding effects, quantitative agreement with experimental data was obtained upon including mild IDP-PEG attraction. The present work demonstrates that realistic modeling of IDPs under crowded conditions for direct comparison with experiments is now achievable.