  We present exploratory results from dynamical simulations of QCD in
isolation, as well as QCD coupled to QED, with C* boundary conditions. In
finite volume, the use of C* boundary conditions allows for a gauge invariant
and local formulation of QED without zero modes. In particular we show that the
simulations reproduce known results and that masses of charged mesons can be
extracted in a completely gauge invariant way. For the simulations we use a
modified version of the HiRep code. The primary features of the simulation code
are presented and we discuss some details regarding the implementation of C*
boundary conditions and the simulated lattice action.
