The CG-MD simulations of the M model and of the BTA model were carried out using the GROMACS software 61 (versions 2018.6 and 2020.2). All the MD simulations were performed in NVT conditions, using a constant volume for the simulation box of 22.056 × 22.056 × 22.056 nm , with periodic boundary conditions, a constant number of molecules N = 500, and a constant temperature. This corresponds to a monomer density of 0.0466 nm -3 , equivalent to a concentration of ∼ 77 mM. The temperature was set to T = 300 K, in the M model, and at different values (see main text) for the BTA model. The systems were simulated in implicit-solvent via Langevin dynamics, accounting for the friction of the solvent and thermal fluctuations of the system. We used the stochastic dynamics (sd) integrator, setting the inverse of the friction constant to tau-t= 0.1 ps. tau-t also sets the coupling with the random force term, that determines the temperature of the system. The During this second part of the CG-MD we sampled the conformations every 300 ps.