In an effort to determine the favorable type and optimum concentration of functional groups in GO structure, we performed molecular dynamics (MD) simulations of a wide range of GO structures, with purely epoxy and hydroxyl functional groups, in contact with S82- polysulfides solvated in a standard electrolyte solvent DME – DOL in 1:1 v/v ratio. Only physical interactions between the polysulfides and the functional groups in GO were considered. Optimized potential for liquid simulations (OPLS) parameters were utilized to model GO and polysulfides along with solvent. Partial charge on carbon atoms in GO with epoxy and hydroxyl groups were evaluated from the optimized structures obtained using Density Functional Theory (DFT). Bond, angle and dihedral parameters along with partial charges on sulfur particles in polysulfide ion S82- were evaluated from the optimized structure using DFT. The calculated density of equilibrated solvent was within 5% of the experimental value. The time averaged distribution of polysulfides, normal to the surface of graphene oxide, was evaluated and showed that polysulfides tend to be closer to GO substrate with hydroxyl functional group as compared to ones with epoxy functional groups. The radial distribution (RDF) plots between H atom of hydroxyl functional group and S atoms of polysulfide ion suggested that a strong Coulombic interaction is responsible for physisorption of polysulfide onto the GO substrate. The diffusion co-efficients of polysulfides were also calculated and were observed to be significantly lower for GO as compared to pure graphene thus indicating a relatively stable deposition. For GO with epoxy functional groups, as the concentration of functional groups was reduced, the physisorption of polysulfides was observed to increase. The polysulfides tended to stabilize at a greater distance from the GO surface with epoxy functional groups. The study also analyzed the role of solvent molecules (especially DOL) in stabilizing polysulfides away from the GO substrate with epoxy functional groups. Overall, the results obtained in the present study show that GO with hydroxyl functional groups have a greater tendency to adsorb polysulfides onto the surface as compared to GO with epoxy groups. These results can potentially pave the way for design of molecularly-tailored cathode supports to mitigate polysulfide shuttle and therefore improve performance of Li-S batteries.