This contribution presents a two-scale computational homogenization framework for the micro-macro simulation of magnetoelectro-mechanically coupled materials. Energetically consistent micro-macro transition conditions will be derived from a generalized form of the classical Hill-Mandel condition. A focus of the work is on the computation of effective magnetoelectric moduli which are derived on the basis of an algorithmically consistent linearization of the macroscopic field equations. The method will be applied to the homogenization of magneto-electric composites which are composed of piezomagnetic and piezoelectric phases. The effective magneto-electric moduli of two-phase composites will be computed.