A generalized convolution for the Fourier sine transform and the Kontorovich - Lebedev integral transform is introduced. Several properties of this generalized convolution and its application to solve certain systems of integral equations are considered.