We propose two different ways to simulate Boolean circuits in the framework of ciliate computations, based on the intramolecular and on the intermolecular gene assembly models with contextual operations applied in a maximal parallel manner. In both cases, the running time and the size of resources of the simulating devices are proportional to the corresponding depth and size of the simulated Boolean circuits.