In this paper, a numerical method is presented for simulating the 3D interfacial flows with insoluble surfactant. The numerical scheme consists of a 3D immersed interface method (IIM) for solving Stokes equations with jumps across the interface and a 3D level-set method for solving the surfactant convection-diffusion equation along a moving and deforming interface. The 3D IIM Poisson solver modifies the one in the literature by assuming that the jump conditions of the solution and the flux are implicitly given at the grid points in a small neighborhood of the interface. This assumption is convenient in conjunction with the level-set techniques. It allows standard Lagrangian interpolation for quantities at the projection points on the interface. The interface jump relations are re-derived accordingly. A novel rotational procedure is given to generate smooth local coordinate systems and make effective interpolation. Numerical examples demonstrate that the IIM Poisson solver and the Stokes solver achieve second-order accuracy. A 3D drop with insoluble surfactant under shear flow is investigated numerically by studying the influences of different physical parameters on the drop deformation.