In this paper, we consider the existence of a pullback attractor for the random dynamical system generated by stochastic two-compartment Gray-Scott equation for a multiplicative noise with the homogeneous Neumann boundary condition on a bounded domain of space dimension n ≤ 3. We first show that the stochastic Gray-Scott equation generates a random dynamical system by transforming this stochastic equation into a random one. We also show that the existence of a random attractor for the stochastic equation follows from the conjugation relation between systems. Then, we prove pullback asymptotical compactness of solutions through the uniform estimate on the solutions. Finally, we obtain the existence of a pullback attractor.