Whilst the initial disaccharide positions were equivalent between each pair of replicate CG and AA simulations, the distribution in disaccharide aggregate sizes varied considerably between the two resolutions. There was a strong preference towards the formation of dimers across all CG disaccharide models, followed by trimers and tetramers. The largest aggregate was observed in one replica of the Glcβ(1→3)Glcβ model, comprising 9 individual disaccharides. This unimodal distribution was consistent between disaccharides that were parameterized using either the fully automated or semi-automated protocols, with the second largest aggregates observed in the Manα(1→2)Manα and Glcβ(1→6)Glcβ models, comprising 8 individual disaccharides. In contrast to this, the AA GLYCAM-06j simulations each showed a bimodal distribution in aggregate sizes, with the primary (largest) peak corresponding to aggregate sizes between 18 -20, and the secondary peak corresponding to aggregate sizes of 2. These data highlight that while disaccharide dimer formation was still present within the AA GLYCAM-06j simulations, each system exhibited a preference towards the formation of larger aggregates comprising most, if not all, disaccharide moieties within the system, contrary to what was expected according to the experimentally derived solubility of disaccharides at this concentration . For this reason, we decided to perform the same all-atom simulations using CHARMM36. Interestingly, for all of the disaccharides studied in the aggregation assays, Martini 3 showed excellent agreement with the CHARMM36 results, for both, average number of aggregates per frame and aggregate dynamics (Figure , Supplementary Figures ). This is in line with the aggregation behaviour of trehalose observed in experiments (i.e. soluble at the same concentration).