Computer simulations of many-body quantum dynamics of indistinguishable particles is a challenging task for computational physics. In this paper we demonstrate that the method of coupled coherent states (CCS) developed previously for multidimensional quantum dynamics of distinguishable particles can be used to study indistinguishable bosons in the second quantisation formalism. To prove its validity, the technique termed here coupled coherent states for indistinguishable bosons (CCSB) is tested on two model problems. The first is a system-bath problem consisting of a tunnelling mode coupled to a harmonic bath, previously studied by CCS and other methods in distinguishable representation in 20 dimensions. The harmonic bath is comprised of identical oscillators, and may be second quantised for use with CCSB, so that this problem may be thought of as a bosonic bath with an impurity. The cross-correlation function for the dynamics of the system and Fourier transform spectrum compare extremely well with a benchmark calculation, which none of the prior methods of studying the problem achieved. The second model problem involves 100 bosons in a shifted harmonic trap. Breathing oscillations in the 1-body density are calculated and shown to compare favourably to a multiconfigurational time-dependent Hartree for bosons calculation, demonstrating the applicability of the method as a new formally exact way to study the quantum dynamics of Bose-Einstein condensates.