In the cyclic guanosine monophosphate (cGMP) signaling pathway, phosphodiesterase 6 (PDE6) maintains a critical balance of the intracellular concentration of cGMP by catalyzing it to 5' guanosine monophosphate (5'-GMP). To gain insight into the mechanistic impacts of the PDE6 somatic mutations that are implicated in cancer and retinitis pigmentosa, we first defined the structure and organization of the human PDE6 heterodimer using computational comparative modelling. Each subunit of PDE6αβ possesses three domains connected through long α-helices. The heterodimer model indicates that the two chains are likely related by a pseudo two-fold axis. The N-terminal region of each subunit is comprised of two allosteric cGMP-binding domains (Gaf-A & Gaf-B), oriented in the same way and interacting with the catalytic domain present at the C-terminal in a way that would allow the allosteric cGMP-binding domains to influence catalytic activity. Subsequently, we applied an integrated knowledge-driven in silico mutation analysis approach to understand the structural and functional implications of experimentally identified mutations that cause various cancers and retinitis pigmentosa, as well as computational saturation mutagenesis of the dimer interface and cGMP-binding residues of both Gaf-A, and the catalytic domains. We studied the impact of mutations on the stability of PDE6αβ structure, subunit-interfaces and Gaf-cGMP interactions. Further, we discussed the changes in interatomic interactions of mutations that are destabilizing in Gaf-A (R93L, V141 M, F162 L), catalytic domain (D600N, F742 L, F776 L) and at the dimer interface (F426A, F248G, F424 N). This study establishes a possible link of change in PDE6αβ structural stability to the experimentally observed disease phenotypes.