Sulfite oxidase (SOX) is a crucial molybdenum cofactor-containing enzyme in plants that re-oxidizes the sulfite back to sulfate in sulfite assimilation pathway. However, studies of this crucial enzyme are quite limited hence this work was attempted to understand the SOXs in four plant species namely, Arabidopsis thaliana, Solanum lycopersicum, Populus trichocarpa and Brachypodium distachyon. Herein studied SOX enzyme was characterized with both oxidoreductase molybdopterin binding and Moco oxidoreductase dimerization domains. The alignment and motif analyses revealed the highly conserved primary structure of SOXs. The phylogeny constructed with additional species demonstrated a clear divergence of monocots, dicots and lower plants. In addition, to further understand the phylogenetic relationship and make a functional inference, a structure-based phylogeny was constructed using normalized RMSD values in five superposed models from four modelled plant SOXs herein and one previously characterized chicken SOX structure. The plant and animal SOXs showed a clear divergence and also implicated their functional divergences. Based on tree topology, monocot B. distachyon appeared to be diverged from other dicots, pointing out a possible monocot-dicot split. The expression patterns of sulfite scavengers including SOX were differentially modulated under cold, heat, salt and high light stresses. Particularly, they tend to be upregulated under high light and heat while being downregulated under cold and salt stresses. The presence of cisregulatory motifs associated with different stresses in upstream regions of SOX genes was thus justified. The protein-protein interaction network of AtSOX and network enrichment with gene ontology (GO) terms showed that most predicted proteins, including sulfite reductase, ATP sulfurylases and APS reductases were among prime enzymes involved in sulfite pathway. Finally, SOX-sulfite docked structures indicated that arginine residues particularly Arg374 is crucial for SOX-sulfite binding and additional two other residues such as Arg51 and Arg103 may be important for SOX-sulfite bindings in plants.