A protocol for explicitly constructing the exact time-evolution operators generated by $2 times 2$ time-dependent $PT$-symmetry Hamiltonians is reported. Its mathematical applicability is illustrated with the help of appropriate examples. The physical relevance of the proposed approach within gain-loss system scenarios, like two-coupled wave-guides, is discussed in detail.