W = N*1e-21*120/Radian | — | W = work; N = number of pN nm rad^{−1} of torque; 120 = degrees of the catalytic step; |

| | Radian = the number of degrees in a radian. |

DeltaMudiss = −DeltaMuATP − W | — | Δ*μ*_{diss}=−Δ*μ*_{ATP}−*W*=65.2e-21 pN nm − Work. |

rateRatio = exp(DeltaMudiss/(kB*T)) | — | Rate ratio = exp(Δ*μ*_{diss}/(k_{b}T); |

| | k_{b}= Boltzmann constant = 1.3806488e-23 J K^{−1}. |

rb = Vmax/(1 + rateRatio) | — | rb=Vmax/(1+exp(Δμdiss/(kBT)); Vmax=35 s−1. |

rf = rateRatio*rb | — | rf=exp(Δμdiss/(kBT))rb |

obsb = rb + contamination | — | Observed backward stepping rate = true catalytic backward stepping rate + contamination; |

| | contamination is identical in both directions and preset at levels ranging from zero to 17.5 s^{−1}. |

obsf = rf + contamination | — | Observed forward stepping rate = true catalytic forward stepping rate + contamination. |

pobsf = obsf/(obsf + obsb) | — | Probability of observing steps in the forward direction as a fraction of the total forward and backward stepping rates observed. This calculated probability is the quantity plotted in the curves sloping *up* from left to right (synthetic direction) in figure 3. |

pobsb = obsb/(obsf + obsb) | — | Probability of observing steps in the backward direction as a fraction of the total forward and backward stepping rates observed. This calculated probability is the quantity plotted in the curves sloping *down* from left to right (hydrolytic direction) in figure 3. |

Two further lines of code were used to calculate the curves plotted in figure 4, thus: |

DeltaMudiss = kB*T*ln(rateRatio) | — | Δ*μ*_{diss}=k_{b}T ln(r_{f}/r_{b}) This calculation is actually redundant, given the second line of code shown above. |

| | It shows the true Probability Isotherm for the uncontaminated catalytic steps. |

DeltaMuobs = kB*T*ln(obsf/obsb) | — | Δ*μ*_{obs}=k_{b}T ln(observed forward stepping rate/observed backward stepping rate). This shows the apparent Probability Isotherm as contaminated with catalytically fruitless steps |