The initial frames for the umbrella sampling (US) simulations were extracted from constantvelocity pulling simulations. In these simulations, the systems were pulled from reaction chain coordinante [133, 134] 𝜉ch = 0.25 to 𝜉ch = 1 over 1000 ns using a force constant of 5000 kJ mol -1 , for five replicates. Nineteen umbrella windows were employed, with reference positions ranging from 𝜉ch = 0.1 to 1 in increments of 0.05, where 1 is the fully formed stalk. The force constant was set to 5000 kJ mol -1 . Each window was simulated for 1 𝜇s, with the initial 500 ns excluded for equilibration. An integration time step of 20 fs was utilized. The potential of mean force (PMF) was computed using the weighted histogram analysis method (WHAM) [135].