It is known that binding (aka association) constants are in fact unitless, as has been discussed here already. However, I'm not a chemist and am confused about when one should or should not use units when working with association constants. One source says:
$K_\text{eq}$ for a reaction with unequal numbers of reactants and products is always given with units, even in published papers.
But why is that? Why not always use unitless values? Is there something inherently wrong with never using units for binding constants?
Consider this example. Let $\ce{A}$ bind to $\ce{X}$ as an n-mer (e.g. $\ce{A}$ can be a transcriptional activator binding to gene promoter). This results in active state, $\ce{X_{A}}$ (e.g. a state leading to gene transcription):
$$\ce{X + nA ->[k_{\text{on}}][k_{\text{off}}] X_{A}}$$
Association constant: $K_\text{A} = \frac{k_\text{on}}{k_\text{off}} $
Assuming equilibrium and law of mass action:
$$K_\text{A} \cdot \text{X} \cdot \mathrm{A^n} = \mathrm{X_A}$$
Now, the fractional occupancy (active states to all states ratio) is:
$$y=\mathrm{\frac{X_A}{X+X_A}} = \frac{K_\text{A}\cdot \mathrm{A^n}}{1 + K_\text{A}\cdot \mathrm{A^n}}$$
As $y$ must be unitless (i.e. it can be interpreted as a probability), for this particular equation, $K_\text{A}$ must have units of $\mathrm{M^{-n}}$ (in general, $\text{concentration}^{-n}$) for the equation to work out, correct (assuming concentration of $\ce{A}$ has units $\mathrm{M}$)? So in this particular case, using unitless $K_\text{A}$ seems wrong, but is it really, or is there something I'm missing?
Now let's extend this equation so that it includes a half-saturation constant $h$, i.e. concentration of $\ce{A}$ required for $y=0.5$ (50% activation). If I'm doing this correctly, we get:
$$y=\frac{K_\text{A}\cdot \mathrm{A^n}}{K_\text{A}\cdot h^n + K_\text{A}\cdot \mathrm{A^n}} = \frac{\mathrm{A^n}}{h^n + \mathrm{A^n}}$$
Note that this is equivalent to the Hill equation for an activator. This generalized equation, unlike the previous one, works just fine regardless of whether or not $K_\text{A}$ is unitless.
Is my understanding of this correct, and does the choice of unitless vs non-unitless binding constant indeed depend on the formulation of a specific equation?