PDE9 inhibitors have been studied to validate their potential to treat diabetes, neurodegenerative disorders, cardiovascular diseases, and erectile dysfunction. In this report, we have selected highly potent previously reported selective PDE9 inhibitors BAY73-6691R, BAY73-6691S, 28r, 28s, 3r, 3s, PF-0447943, PF-4181366, and 4r to elucidate the differences in their interaction patterns in the presence of different metal systems such as Zn/Mg, Mg/Mg, and Zn/Zn. The initial complexes were generated by molecular docking followed by molecular dynamics simulation for 100 ns in triplicate for each system to understand the interactions’ stability. The results were carefully analyzed, focusing on the ligands’ non-bonded interactions with PDE9 in different metal systems.