We motivate the need to include collisional dissipation in gyrokinetic turbulence simulations, and constract criteria for a physically vaUd model of such dissipation. A new analytically manageable operator satisfying those criteria is presented and transformed into gyrokinetic variables. The form of conservation laws for collision operators in gyrokinetic variables is explained. The numerical implementation of our new operator in the code GS 2 is outlined and successful tests from this code presented.