The behaviors of annular sector plates made of functionally graded graphene origami-enabled auxetic metamaterials (FG-GOEAMs) under the action of buckling and bending loads are investigated herein using a numerical approach. Multiple GOEAM layers with graphene origami (GOri) content are considered for the plate which changes in layer-wise patterns. Moreover, the material properties of the plate are calculated using novel genetic programming-assisted micromechanical models. In order to derive the governing equations, the Mindlin plate theory in conjunction with Hamilton’s principle is utilized within the framework of the variational differential quadrature (VDQ) method. The vector-matrix representation of the presented formulation can be beneficial from the viewpoint of implementing numerical approaches. The governing equations are then discretized and solved based on VDQ matrix differential and integral operators so as to obtain the critical buckling load and maximum lateral deflection of plates with different edge conditions. The effects of GOri content, folding degree and distribution pattern on the results are studied. It is revealed that increasing the GOri content leads to the negative Poisson’s ratio (NPR) effect which respectively reduces and increases the lateral deflection and the critical buckling load.