Direct imaging of exoplanets greatly enhances our ability to spectroscopically detect molecules in their atmospheres, but requires both excellent angular and spectral resolution. A highly promising approach uses multiplexed transmission gratings to achieve compressed sensing, but the published theoretical calculations are not realistic enough to determine whether the approach is feasible. We aim to determine the performance of recovering exoplanet signals of an instrument with multiplexed volume Bragg gratings in a full, detailed numerical simulation, specifically examining any effects caused by multiplexing several gratings. Our end-to-end simulation includes realistic stellar and planetary spectra, a closed-loop adaptive optics simulation, and rigorous coupled wave analysis to model the multiplexed gratings. We chose 20 passbands around expected methane absorption features that were optically combined using 20 gratings. We find that exoplanet signals can be recovered down to contrasts of $10^ $, without the addition of a coronagraph. Multiplexing a larger number of gratings improves the deepest recovered contrast, if photon noise is the dominant noise source. When residuals from our simple post-processing approach are the largest noise source, a larger number of multiplexed gratings decreases the performance, due to the stacking of the PSFs at different wavelengths. This is an artifact of our data reduction approach. We conclude that multiplexed Bragg gratings are a viable method to look for exoplanets with a compressed sensing approach and additional gains may be made in post-processing.