We note here that to model a disaccharide using our procedure a bead that includes the glycosidic bond is combined into one of the hexose units. This is required to represent the glycosidic link between the monosaccharides and leads to a slightly different mapping compared to the original individual monosaccharides reported by Grunewald et al . Consequently, the bonded terms are also slightly different. We have compared the bonded terms from our trehalose model to that reported in the work of Grunewald et al , and achieve good agreement (more details are in the result section, Supplementary Table and Supplementary Figure ). The next part of the procedure involves performing a set of standard analyses to assess the new parameter set including evaluation of distribution of the bond lengths, angles and dihedrals and solvent accessible surface area (much of this is automated). Given space constraints, the data for 3 out of the 28 disaccharides that were parameterized using this fully automated procedure are shown in Figure , all of the other data is available in the SI.