For all CG simulations, the Martini 3 force field was used. C8S3 monomers and aggregates were solvated in Martini 3 water beads (W, representing 4 H 2 O molecules), and Na + (TQ5 bead with +1 charge) were used to neutralise the charge of each system. In all systems, the Berendsen barostat 39 was used to maintain the pressure constant at 1 bar with a time constant of 1 ps and compressibility of 3x10 -4 bar -1 . The temperature was kept constant at 300 K by using the v-rescale algorithm with a time constant of 0.1 ps . Electrostatic and van der Waals interactions were calculated using the Verlet list scheme with the cut-off set at 1.1 nm and a buffer tolerance of 0.01 kJ mol -1 ps -1 . Long range interactions were treated using the Reaction Field method and the bond lengths were constrained by using the LINCS method .