Microbial redox reactions of inorganic sulfur compounds, mainly the sulfur anions, are one of the vital reactions responsible for the environmental sulfur balance. These reactions are mediated by phylogenetically diverse prokaryotes, which also take part in the extraction of metal ions from their sulfur containing ores. The sulfur oxidizing gene cluster (sox) of alpha-Proteobacteria comprises of at least 16 genes, forming two transcriptional units, viz., soxSRT and soxVWXYZABCDEFGH. SoxY is known to be a sulfur covalently binding protein, which binds sulfur anions (such as sulfate) to form SoxY-thiocysteine-S-sulfate, the first covalently bound sulfur adduct in the novel global sulfur anion oxidation cycle. SoxZ, a sulfur compound chelating protein, binds to SoxY forming a complex to which SoxB, a sulfate thiol-esterase, binds and ultimately cleaves the sulfur adduct. We employed homology modeling to construct the three-dimensional structures of the SoxY, SoxZ, and SoxB from Paracoccus pantotrophus. With the help of docking and molecular dynamics studies we have identified the residues of SoxY, SoxZ, and SoxB involved in the interaction. The probable mechanisms of the binding of SoxY with Sulfate as well as the removal of sulfate from the SoxYZ complex are also established. Our study provides a rational basis to illustrate the molecular mechanism of the biochemistry of sulfur anion oxidation reactions by these industrially important organisms. (c) 2005 Elsevier Inc. All rights reserved.