To quantify the solvent exposure of these sites, we used the double cubic lattice method (DCLM) to calculate the SASA of the CONH atoms, which were the atoms sharing the same plane with the scissile bond, between PheB25 and TyrB26, and between TyrB26 and ThrB27. In DCLM, the accessible surface is defined by tracing the center of the probe sphere (with the radius as the van der Waals radius of water) as it rolls along the van der Waals surface of the solute. The obtained SASA time series was truncated to discard the non-equilibrium region and was then decorrelated . The mean and the standard deviation of the time series were calculated for each glycoform. Finally, the reported value for each glycoform is the value averaged across the five different wild-type bases, with the uncertainty being the standard deviation of the 5 runs. We then plotted the 𝛼-chymotrypsin half-life measured in the experimental work against these averaged SASA values in a correlation plot, with the Kendall's tau correlation coefficient (𝜏) and its uncertainty annotated. The uncertainty of the correlation coefficient was estimated from bootstrapping, as detailed in section 2.3. Note that Kendall's tau was adopted in favor of Pearson correlation coefficient because it was not clear if the relationship between the two variables should be linear, so we assume only a monotonic relationship.