The essential biological function of phosphodiesterase (PDE) type enzymes is to regulate the cytoplasmic levels of intracellular second messengers, 3',5'-cyclic guanosine monophosphate (cGMP) and/or 3',5'-cyclic adenosine monophosphate (cAMP). PDE targets have 11 isoenzymes. Of these enzymes, PDE5 has attracted a special attention over the years after its recognition as being the target enzyme in treating erectile dysfunction. Due to the amino acid sequence and the secondary structural similarity of PDE6 and PDE11 with the catalytic domain of PDE5, first-generation PDE5 inhibitors (i. e. sildenafil and vardenafil) are also competitive inhibitors of PDE6 and PDE11. Since the major challenge of designing novel PDE5 inhibitors is to decrease their cross-reactivity with PDE6 and PDE11, in this study, we attempt to identify potent tadalafillike PDE5 inhibitors that have PDE5/PDE6 and PDE5/PDE11 selectivity. For this aim, the similarity-based virtual screening protocol is applied for the "clean drug-like subset of ZINC database" that contains more than 20 million small compounds. Moreover, molecular dynamics (MD) simulations of selected hits complexed with PDE5 and off-targets were performed in order to get insights for structural and dynamical behaviors of the selected molecules as selective PDE5 inhibitors. Since tadalafil blocks hERG1 K channels in concentration dependent manner, the cardiotoxicity prediction of the hit molecules was also tested. Results of this study can be useful for designing of novel, safe and selective PDE5 inhibitors.