Schindler disease is a lysosomal storage disorder caused due to deficiency or defective activity of alpha-N-acetylgalactosaminidase (α-NAGA). Mutations in gene encoding α-NAGA cause wide range of diseases, characterized with mild to severe clinical features. Molecular effects of these mutations are yet to be explored in detail. Therefore, this study was focused on four missense mutations of α-NAGA namely, S160C, E325K, R329Q and R329W. Native and mutant structures of α-NAGA were analysed to determine geometrical deviations such as the contours of root mean square deviation, root mean square fluctuation, percentage of residues in allowed regions of Ramachandran plot and solvent accessible surface area, using conformational sampling technique. Additionally, global energy-minimized structures of native and mutants were further analysed to compute their intra-molecular interactions, hydrogen bond dilution and distribution of secondary structure. In addition, docking studies were also performed to determine variations in binding energies between native and mutants. The deleterious effects of mutants were evident due to variations in their active site residues pertaining to spatial conformation and flexibility, comparatively. Hence, variations exhibited by mutants, namely S160C, E325K, R329Q and R329W to that of native, consequently, lead to the detrimental effects causing Schindler disease. This study computationally explains the underlying reasons for the pathogenesis of the disease, thereby aiding future researchers in drug development and disease management. © 2016, International Association of Scientists in the Interdisciplinary Areas and Springer-Verlag Berlin Heidelberg.