The purpose of this work is to evaluate the homology modeling, in silico prediction, and characterisation of Cytochrome c oxidase from Cyprinus carpio and Tubifex tubifex, as well as molecular docking experiments between the modelled protein and three frequently used surfactants. Using the template crystal structure of bovine heart Cytochrome c oxidase, homology modeling of Cytochrome c oxidase (Subunit 2) of Cyprinus carpio (Accession # P24985) and Cytochrome c oxidase (Subunit 1) of Tubifex tubifex (Accession # Q7YAA6) was conducted. The model structure was improved further with 3Drefine, and the final 3D structure was verified with PROCHEK and ERRATA. The physiochemical, as well as the stereochemical parameters of the modelled protein, were evaluated using various tools like ExPASy’s ProtParam, Hydropathy Analysis and EMBOSS pepwheel. The projected model was then docked with toxic ligands, Sodium dodecyl sulfate (SDS), Cetylpyridinium chloride (CPC), and Sodium laureth sulfate (SLES), whose 3D structures were obtained from the Uniprot database. CPC interacted best with Cytochrome c oxidase subunit 2 of Cyprinus carpio and Cytochrome c oxidase subunit 1 of Tubifex tubifex, according to our findings. Furthermore, in the case of all surfactants, hydrophobic interactions with the active site amino acid residues of the modelled protein were observed to be more common than hydrogen bonds and salt bridges. Molecular simulation studies exhibited that the surfactants alter the structural flexibility of the predicted proteins. Hence it may be inferred that the surfactants might alter the structure and dynamics of Cytochrome c oxidase of both worm and fish.