Using a new Bayesian Optimization algorithm to guide the design of mechanical metamaterials, we design nonhomogeneous 3D structures possessing the Cauchy symmetry, which dictates the relationship between continuum and atomic deformations. Recent efforts to merge optimization techniques with the design of mechanical metamaterials has resulted in a concentrated effort to tailor their elastic and post elastic properties. Even though these properties of either individual unit cells or homogenized continua can be simulated using multi-physics solvers and well established optimization schemes, they are often computationally expensive and require many design iterations, rendering the validation stage a significant obstacle in the design of new metamaterial designs. This study aims to provide a framework on how to utilize miniscule computational cost to control the elastic properties of metamaterials such that specific symmetries can be accomplished. Using the Cauchy symmetry as a design objective, we engineer structures through the strategic arrangement of 5 different unit cells in a 5×5×5 cubic symmetric microlattice structure. This lattice design, despite constituting a design space with 510 3D lattice configurations, can converge to an effective solution in only 69 function calls as a result of the efficiency of the new Bayesian optimization scheme. To validate the mechanical behavior of the design, the lattice structures were fabricated using multiphoton lithography and mechanically tested, revealing a close correlation between experiments and simulated results in the elastic regime. Ultimately, a similar methodology can be utilized to design metamaterials with other material properties, aspiring to control properties at different length scales, an endeavor that requires inordinate computation cost.