Motivated by significant discrepancies between experimental observations of electron-doped cuprates and numerical results of the Hubbard and models, we investigate the role of intersite interactions by studying the model on square lattices. Based on large-scale density matrix renormalization group simulations, we identify the ground-state phase diagram across varying intersite interactions and doping concentration . We find that the phase diagram with finite intersite interactions offers a more accurate description of electron-doped cuprates than the conventional Hubbard and models. Moreover, we reveal the role of intersite interactions at varying doping levels: at light doping, intersite interactions favor Néel antiferromagnetic order, and suppress both superconductivity and charge density wave; around optimal doping, these interactions support a pseudogap-like phase while suppressing superconductivity, and we further perform the slave-boson mean-field analysis to understand the numerical results microscopically; at higher doping, the effects of intersite interactions become insignificant, with our numerical predictions suggesting the emergence of an incommensurate spin density wave phase. Our specific focus around optimal doping with various intersite interactions identifies successive phases including phase separation, uniform -wave SC, and a pseudogap-like phase, and reveals a relative insensitivity of the charge density wave to superconductivity. Our study suggests the model as the minimal model to capture the essential physics of electron-doped cuprates.