All the analysis were carried out by means of python scripts employing the MD-Analysis package , with the exception of the analysis reported in Figures , S8, S9, S10 and the analysis reported in Figure , which were performed using the clustsize and rdf GROMACS tools, respectively. The snapshots in Figures -b and S11 were rendered with VMD. For the analyses of monomer aggregation, we proceeded as follows: two M monomers were considered in contact (i.e., they belonged to the same assembly) if their core beads lied within a distance of r cut = 0.6 nm; we chose this value for r cut as this radius includes the first peak of the radial distribution function (g(r)) (Figure ). For the BTA model, the contact radius r cut has been fixed at 0.6 nm, obtained following the same criterion. We also tested how a different choice of the cutoff can impact the analysis (vide infra).