The effect of non-denaturing concentrations of three different organic solvents, formamide, acetone and isopropanol, on the structure of haloalkane dehalogenases DhaA, LinB, and DbjA at the protein-solvent interface was studied using molecular dynamics simulations. Analysis of B-factors revealed that the presence of a given organic solvent mainly affects the dynamical behavior of the specificity-determining cap domain, with the exception of DbjA in acetone. Orientation of organic solvent molecules on the protein surface during the simulations was clearly dependent on their interaction with hydrophobic or hydrophilic surface patches, and the simulations suggest that the behavior of studied organic solvents in the vicinity of hyrophobic patches on the surface is similar to the air/water interface. DbjA was the only dimeric enzyme among studied haloalkane dehalogenases and provided an opportunity to explore effects of organic solvents on the quaternary structure. Penetration and trapping of organic solvents in the network of interactions between both monomers depends on the physico-chemical properties of the organic solvents. Consequently, both monomers of this enzyme oscillate differently in different organic solvents. With the exception of LinB in acetone, the structures of studied enzymes were stabilized in water-miscible organic solvents.