We investigate, via Brownian dynamics simulations, the reaction dynamics of a generic, nonlinear chemical network under spatial confinement and crowding conditions. In detail, the Willamowski-Rossler chemical reaction system has been "extended" and considered as a prototype reaction-diffusion system. Our results are potentially relevant to a number of open problems in biophysics and biochemistry, such as the synthesis of primitive cellular units (protocells) and the definition of their role in the chemical origin of life and the characterization of vesicle-mediated drug delivery processes. More generally, the computational approach presented in this work makes the case for the use of spatial stochastic simulation methods for the study of biochemical networks in vivo where the "well-mixed" approximation is invalid and both thermal and intrinsic fluctuations linked to the possible presence of molecular species in low number copies cannot be averaged out.