Typical MRI simulators are either based on the signal equations or on the solution to the Bloch equations for each voxel, with appropriate image encoding incorporated. In this work, a novel modular approach to simulate MR sequences, termed quantitative voxel-wise simulator (QVS), is proposed. This simulator employs quantitative parametric maps as inputs, models the evolved MR signal as a multi-dimensional filter output and reconstructs the image by using a voxel-by-voxel modulated k-space approach. As a proof-of-concept, 3D MP-RAGE sequence is simulated and the simulated images are compared with in-vivo data. QVS is designed for multi-centre brain MRI harmonization.