A new approach to carry out molecular dynamics simulations of chemical reactions in solution using combined density functional theory/molecular mechanics potentials is presented. We focus our attention on the analysis of reactive trajectories, dynamic solvent effects and transmission coefficient rather than on the evaluation of free energy which is another important topic that will be examined elsewhere. In a previous paper we have described the generalities of this hybrid molecular dynamics method and it has been employed to investigate low energy barrier proton transfer process in water. The study of processes with activation energies larger than a few kT requires the use of specific techniques adapted to "rare events" simulations. We describe here a method that consists in the simulation of shea trajectories starting from an equilibrated transition state in solution, the structure of which has been approximately established. This calculation is particularly efficient when carried out with parallel computers since the study of a reactive process is decomposed in a set of short time trajectories that are completely independent. The procedure is close to that used by other authors in the context of classical molecular dynamics but present the advantage of describing the chemical system with rigorous quantum mechanical calculations. It is illustrated through the study of the first reaction step ir; electrophilic bromination of ethylene in water. This elementary process is representative of many charge separation reactions for which static and dynamic solvent effects play a fundamental role.