Using large-scale quantum Monte Carlo simulations, we exactly solve a model of Fermions hopping on the honeycomb lattice with cluster charge interactions, which has been proposed as an effective model with possible application to twisted bilayer graphene near half-filling. We find an interaction driven semimetal to insulator transition to an insulating phase consisting of a valence bond solid with Kekule pattern. Finite size scaling reveals that the phase transition of the semimetal to Kekule valence bond solid phase is continuous and belongs to chiral XY universality class.