Our Galaxy and the nearby Andromeda Galaxy (M31) form a bound system, even though the relative velocity vector of M31 is currently not well constrained. Their orbital motion is highly dependent on the initial conditions, but all the reliable scenarios imply a first close approach in the next 3–5 Gyrs. In our study, we simulate this interaction via direct N-body integration, using the HiGPUs code. Our aim is to investigate the dependence of the time of the merger on the physical and dynamical properties of the system. Finally, we study the dynamical evolution of the two Supermassive Black Holes placed in the two galactic centers, with the future aim to achieve a proper resolution to follow their motion until they form a tight binary system.