The dynamics of sheared wet foams is investigated computationally in this study. For this purpose, an established 3D (parallel) implementation of a level-set method for incompressible two-phase flow has been extended to account for the presence of surfactants that are soluble in the liquid, and the associated modified stress conditions at interfaces. In particular, we account for surface rheological behavior (such as surface viscosity) beyond minimal surfactant models, to describe realistic systems. The results of 2D and 3D tests will be demonstrated to compare favourably with the literature. We shall report on the role of surfactants and their properties on T1 events, wherein adjacent bubbles are sheared past each other, with implications for foam instability.