For the comprehensive simulation of ion trajectories including reactive collisions at elevated pressure conditions, a chemical reaction simulation (RS) extension to the popular SIMION software package was developed, which is based on the Monte Carlo statistical approach. The RS extension is of particular interest to SIMION users who wish to simulate ion trajectories in collision dominated environments such as atmospheric pressure ion sources, ion guides (e.g., funnels, transfer multi poles), chemical reaction chambers (e.g., proton transfer tubes), and/or ion mobility analyzers. It is well known that ion molecule reaction rate constants frequently reach or exceed the collision limit obtained from kinetic gas theory. Thus with a typical dwell time of ions within the above mentioned devices in the ms range, chemical transformation reactions are likely to occur. In other words, individual ions change critical parameters such as mass, mobility, and chemical reactivity en passage to the analyzer, which naturally strongly affects their trajectories. The RS method simulates elementary reaction events of individual ions reflecting the behavior of a large ensemble by a representative set of simulated reacting particles. The simulation of the proton bound water cluster reactant ion peak (RIP) in ion mobility spectrometry (IMS) was chosen as a benchmark problem. For this purpose, the RIP was experimentally determined as a function of the background water concentration present in the IMS drift tube. It is shown that simulation and experimental data are in very good agreement, demonstrating the validity of the method.