Yes, of course, cubic symmetry is a hugely constraining condition on the coefficients $A_{nlm}$ - I suppose that you meant that the coefficients depend on $l$ as well.
The cubic symmetry is generated by a few generators such as some rotations by 90 degrees. For example, if the state is symmetric with respect to the 90-degree rotation around the $z$ axis, it means that $\exp(\pi i J_z/2)$ has to keep the state invariant. It implies that $m$ i.e. the eigenvalue of $J_z$ has to be a multiple of four. All coefficients $A_{nlm}$ for which $m$ is not a multiple of four have to vanish.
If the state is also symmetric with respect to the $z\to(-z)$ reflection (parity), it means that $l$ has to be even: all coefficients with $l$ odd vanish.
Finally, if the state is invariant under the rotation by 90 degrees around the $x$ axis, there is a similar condition for the coefficients that is slightly more complicated than that $m$ has to be a multiple of four (the condition will constrain some linear combinations of the coefficients), but it is essentially the same thing. If you impose this extra condition, you will find out that only 1/48 of the coefficients may be nonzero because 48 is the order of the group. (Replace 48 by 24 if you don't include the parity-odd transformations.)