I think you’re in the same case as this post, where you want to do a convolution of binomials. This results in a Poisson-Binomial distribution.
Unlike the linked post, you have small n, so you should be able to implement the distribution by implementing the logp on the Wikipedia page and then passing it to a pm.CustomDist