The chemical feedback between building blocks in templated polymerization of diblock copolymers and their consecutive micellization was studied for the first time by means of coarse-grained molecular dynamics simulations. Using a stochastic polymerization model, we were able to reproduce the experimental findings on the effect of chemical feedback on the polymerization rates at low and high solution concentrations. The size and shape of micelles were computed using a newly developed software in Python conjugated with graph theory. In full agreement with the experiments, our simulations revealed that micelles formed by the templated micellization are more spherical and have a lower radius of gyration than those formed by the traditional two-step micellization method. The advantage of molecular simulation over the traditional kinetic models is that with the simulation, one studies in detail the heterogeneous polymerization in the presence of the oppositely charged template while also accounting for the incompatibility between reacted species, which significantly influences the reaction process.