1 H and 31 P-NMR spectra were acquired using a Bruker Ascend Avance Neo operating at 500 MHz and 202.4 MHz, respectively. Samples consisting of H 2 O:D 2 O (9:1, v/v) mixtures were analysed using HOD suppression to collect 1 H-NMR data. Chemical shifts (δ) are given in ppm, the notations s, d, t, q, and quint represent the multiplicities of singlet, doublet, triplet, quartet, and quintet, respectively, and coupling constants (J) are given in Hertz (Hz). Yields of conversion were determined by relative integration of the signals in the 1 H or 31 P-NMR spectra. Data analysis was performed using MestReNova (version 14.3.0) and GraphPad Prism (version 9.5.1).