  The weak decays $K^\pm\to\pi^\pm a$ offer a powerful probe of axion-like
particles (ALPs). In this work, we provide a comprehensive analysis of these
processes within chiral perturbation theory, extending existing calculations by
including complete next-to-leading order (NLO) contributions and
isospin-breaking corrections at first order in $(m_d-m_u)$. We show that the
consistent incorporation of ALPs in the QCD and weak chiral Lagrangians
requires a non-trivial extension of the corresponding operator bases, which we
describe in detail. Furthermore, we show that in the presence of an ALP the
so-called ``weak mass term'', which is unobservable in the Standard Model, is
non-redundant already at leading order. We find that NLO corrections associated
with flavor-violating ALP couplings modify the leading-order result by a few
percent, with negligible uncertainties. NLO corrections proportional to
flavor-conserving ALP couplings lead to potentially larger corrections, which,
however, are accompanied by sizable uncertainties mainly due to the currently
limited knowledge of various low-energy constants. We study how these
corrections impact bounds on the ALP couplings, first model independently, and
then specializing to the case of an ALP with flavor-universal couplings in the
UV. Our findings confirm that the decays $K^\pm\to\pi^\pm a$ provide the
strongest particle-physics constraints for $m_a\lesssim 300$\,MeV. In addition,
we point out that these bounds have interesting implications for the ALP
couplings to nucleons, which were so far only constrained by astrophysical
measurements and non-accelerator experiments.
