The Color Glass Condensate (CGC) is an effective field theory describing partonic interactions inside the nucleus at high energies (small-x). In this regime, partons with small longitudinal momentum fraction exhibit rapid growth, which is tamed by non-linear QCD effects, leading to a saturated state of the nucleus. In this talk, I will present our progress towards the computation of dijet production in pA scattering at small-x at next-to-leading order (NLO). This computation is essential to unveil the effects of saturation and disentangle them from other competing mechanisms. In particular, I will present results on trijet production in the gluon-initiated channel in pA collisions, and on NLO real corrections for an integrated (anti)quark in the back-to-back limit, where saturation effects are expected to be stronger.