The binding interactions of the surfactants: anionic sodium dodecyl sulphate (SDS), cationic cetyltrimethylammonium bromide (CTAB), non-ionic octyl glucoside (OG), and zwitterionic 3-[Hexadecyl(dimethyl)ammonio]-1-propanesulfonate (HPS), with bovine serum albumin (BSA) were investigated by computer simulation. The results disclosed that the surfactants bound stably between hydrophobic subdomain IIA and IIIA where tryptophan-213 residue, an important intrinsic fluorophore in BSA is housed. The interactions of the surfactants with the BSA were electrostatic and hydrophobic interactions. The head-groups of SDS, HPS and OG formed hydrogen bonds with the BSA, while that of CTAB was shielded from intermolecular hydrogen-bonding due to intervening methyl groups. Subsequently, molecular dynamics (MD) simulation of the protein-surfactant complexes revealed that hydrogen bonds formed by OG were stronger than those of SDS and HPS. However, the decomposed force-field energies showed that OG had the least interaction energy with the BSA. In addition to MD simulation, it was found by density functional theory (DFT) that the differences in the coulomb interaction energies can be attributed to charge distribution in the surfactants. Overall, free energies calculated by linear interaction energy (LIE) proved that the binding of each surfactant was dominated by differences between van der Waals interactions in bound and free states.