We present an efficient implementation of analytical non-adiabatic derivative coupling elements for the coupled cluster singles and doubles model. The derivative coupling elements are evaluated in a biorthonormal formulation in which the nuclear derivative acts on the right electronic state, where this state is biorthonormal with respect to the set of left states. This stands in contrast to earlier implementations based on normalized states and a gradient formula for the derivative coupling. As an illustration of the implementation, we determine a minimum energy conical intersection between the nπ* and ππ* states in the nucleobase thymine.