In short, the optimized geometries were clustered as described in Wicker and Cooper 49 using an RMSD threshold of 2 Å. Due to the high sensitivity of Boltzmann distribution exponential errors in relative energies, we opted to weight all conformers (within 10 kcal mol from the global minimum) equally and calculate the probability for each cluster as :