G protein-coupled receptors (GPCRs) are important targets in drug discovery because of their roles in physiological and pathological processes. Orphan GPCRs are GPCR proteins for which endogenous ligands have not yet been identified and they present interesting avenues for therapeutic intervention. This study focuses on GPR78, an orphan GPCR that is expressed in the central nervous system and linked to neurological disorders. GPR78 has no reported crystal structure and there is limited research. In this study, we have predicted the three dimensional model of GPR78 and its probable binding pocket. Structure-based virtual screening was carried out using the ChemDiv and Enamine REAL databases, followed by induced-fit docking studies to identify potential lead compounds with favorable interactions. These lead compounds were then embedded into a POPC lipid bilayer for a 200 ns molecular dynamics simulation. Free energy landscapes and MM-PBSA analyses were performed to assess the binding energies and conformational dynamics. The results highlight the dynamic nature of GPR78 in the presence of lead compounds and show favorable binding interactions. This study aims to predict a reliable three dimensional model of GPR78 and identify novel lead compounds through a comprehensive in silico approach. The identification of these potential GPR78 agonists represents a significant step in the development of new therapeutics for neurological disorders, highlighting the therapeutic potential of orphan GPR78 in CNS disorders.