A sequential nonlinear multiscale method for the simulation of elastic metamaterials subject to large deformations and instabilities is proposed. For the finite strain homogenization of cubic beam lattice unit cells, a stochastic perturbation approach is applied to induce buckling. Then, three variants of anisotropic effective constitutive models built upon artificial neural networks are trained on the homogenization data and investigated: one is hyperelastic and fulfills the material symmetry conditions by construction, while the other two are hyperelastic and elastic, respectively, and approximate the material symmetry through data augmentation based on strain energy densities and stresses. Finally, macroscopic nonlinear finite element simulations are conducted and compared to fully resolved simulations of a lattice structure. The good agreement between both approaches in tension and compression scenarios shows that the sequential multiscale approach based on anisotropic constitutive models can accurately reproduce the highly nonlinear behavior of buckling-driven 3D metamaterials at lesser computational effort.