As I understand it, your question which was marked as a duplicate not only was not a duplicate of any previous question, but also seems to contain erroneous and/or misleading answers.
Let's start with the velocity distribution. Wikipedia says the following: "The Maxwell–Boltzmann distribution applies to the classical ideal gas, which is an idealization of real gases. In real gases, there are various effects (e.g., van der Waals interactions, relativistic speed limits, and quantum exchange interactions) that make their speed distribution sometimes very different from the Maxwell–Boltzmann form."
Contrast this with (my previous comment about a derivation and) the passage found in Landau Lifshitz Statistical Physics (pp. 79-80): "Let us consider the probability distribution for the momenta, and once again emphasize the very important fact that in classical statistics this distribution does not depend on the nature of the interaction of particles within the system or on the nature of the external field, so can be expressed in a form applicable to all bodies. (In quantum statistics this statement is not quite true in general.)" After a derivation of the distribution, the book then goes on explaining that it works also in molecular systems and in systems with Brownian motion. They also show the difference between quantum statistics and the classical description of a harmonic oscillator and how one gets different velocity distributions of the particle (in the latter case accordingly recovers the Maxwell-Boltzmann distribution).
Supposing that we can characterize everything by classical mechanics (a rather fair assumption given the application, certainly one made routinely in molecular dynamics simulations of biological matter), the velocity distributions ought to follow Maxwell-Boltzmann. Now is this to say that Gillespie's formulation is without any fault? No. You are entirely correct in being suspicious of the assumption of no velocity correlation between collisions. However, due to the abundant aqueous solution (and relatively dilute reactants), the non-reactive collisions between the chemical compound and water quickly thermalize and destroy any previous velocity correlations between the two reagants (you might want to see if hydrodynamical interactions exist; I doubt that they are of any relevance in 3D, but for diffusion on top of the cell membrane they might play a minor role. An in-depth analysis would probably even make a decent publication.). The spatial correlation, however, probably takes a longer time to vanish, and it is here that Gillespie's formulation is at its weakest (the Gillespie method can be supplanted by a true continuum correction, i.e. Green's Function Reaction Dynamics, or by dividing space into small boxes and having chemical constants for diffusing from one box to another). For a recent application where spatial heterogeneity is important see for example PNAS 110, 5927 (2013) and references contained therein.