In vitro facsimiles of biomolecular condensates are formed by different types of intrinsically disordered proteins including prion-like low complexity domains (PLCDs). PLCD condensates are viscoelastic materials defined by time-dependent, sequence-specific complex shear moduli. Here, we show that viscoelastic moduli can be computed directly using a generalization of the Rouse model and information regarding intra- and inter-chain contacts that is extracted from equilibrium configurations of lattice-based Metropolis Monte Carlo (MMC) simulations. The key ingredient of the generalized Rouse model is the Zimm matrix that we compute from equilibrium MMC simulations. We compute two flavors of Zimm matrices, one referred to as the single-chain model that accounts only for intra-chain contacts, and the other referred to as a collective model, that accounts for inter-chain interactions. The single-chain model systematically overestimates the storage and loss moduli, whereas the collective model reproduces the measured moduli with greater fidelity. However, in the long time, low-frequency domain, a mixture of the two models proves to be most accurate. In line with the theory of Rouse, we find that a continuous distribution of relaxation times exists in condensates. The single crossover frequency between dominantly elastic versus dominantly viscous behaviors is influenced by the totality of the relaxation modes. Hence, our analysis suggests that viscoelastic fluid-like condensates are best described as generalized Maxwell fluids. Finally, we show that the complex shear moduli can be used to solve an inverse problem to obtain distributions of relaxation times that underlie the dynamics within condensates.