For the ordered and random GO models we found that there is no significant structural variation upon solvation. The oxygen functions (i.e., epoxide and hydroxide groups) form 30-70 % of the possible HBs with the surrounding water molecules, with smaller values for the respective ordered GO models. Based on geometrical criteria, most of GO-water HBs are of moderate strength.