where r is the distance between bonded CG beads and σ is the characteristic size of a Lennard-Jones bead. The FENE potential parameters 𝐾 𝑟 and 𝑅 0 are specified to 30.0ε and 1.5 σ , respectively, where ε is the minimal potential energy value or depth of the potential well. A flexible backbone and sidechains are used without any constraint on the bending rigidity. Interactions between non-bonded CG beads are described using the LJ potential with a cutoff distance of 2.5σ. Unless otherwise noted, all MD simulations were carried out with the Largescale Atomic/Molecular Massively Parallel Simulator (LAMMPS). The mechanical properties were evaluated with non-equilibrium molecular dynamics (NEMD) simulations. The FENE potential inadequately accounted for chain scission during the application of external strain in the NEMD simulations; therefore, a quartic potential was employed to describe the interactions among bonded CG beads in the straining MD simulations, as represented by Equation 2.