We present a new unintegrated (transverse momentum dependent, TMD) gluon density in nuclei valid at small $x$ and incorporating the shadowing effect. We start from recently presented TMD gluon distribution in proton (LLM-2024), which shows saturation behavior, and extend it to nuclei using the geometric scaling. The saturation scale is rescaled in a way to describe available data on nuclear structure functions at low $x$ and $Q^2$ using color dipole approach. This method allows to take into account the shadowing effect at low $x$ in the TMD gluon distribution.