One problem with your formula that $T$ factorizes as the product of a unitary operator $U$ and complex conjugation $K$ is that 'complex conjugation' is meaningless in a Hilbert space a priori.
Let me be more formal about this; consider a vector $\psi \in H$, with $H$ a $n$-dimensional Hilbert space, that is, isomorphic to $\mathbb{C}^n$. How is an (abstract) vector space isomorphic to the canonical Hilbert space $\mathbb{C}^n$? In the following way: pick a basis $\{e_i\}_{i=1,...,n}$ of $H$. Then decompose $\psi$ in this basis:
$$ \psi = \sum_i \psi_i e_i \ .$$
Finally map $\psi$ on the $n$-tupel $(\psi_1,...,\psi_n)$.
As you can see, this isomorphism, lets denote it by $E$, very strongly depends on the chosen basis.
Let's ignore this anyway. On $\mathbb{C}^n$ we can define complex conjugation $K$, it is simply the map
$$ K(\psi_1,...,\psi_n) = (\overline{\psi_1},...,\overline{\psi_n}) \ . $$
Hence we can define a complex conjugation $K_E$ on $H$ simply by
$$ K_E = E^{-1} \circ K \circ E \ .$$
Now let's see what happens if we change a basis; consider another basis $E' = \{e'_i\}_{i=1,...,n}$ with $e_i = \sum_j M_{ij} e'_j$. Then considering $M_{ij}$ as the matrix of an operation on $\mathbb{C}^n$:
$$ K_{E} = E^{-1} \circ K \circ E = E'^{-1} \circ M^{-1} \circ K \circ M \circ E' = E'^{-1} \circ \overline{M}M^{-1} K \circ E' \ ,$$
i.e. unless $\overline{M} = M$, $K_E \neq K_{E'}$. Hence there is no invariant notion of complex conjugation in a complex Hilbert space.
The statement is true that if we fix a basis $E$, we may write complex conjugation w.r.t. any other basis as
$$ K_{E'} = U_{E,E'} K_E \ ,$$
where $U_{E,E'}$ is a unitary. Note that complex conjugation depends on $E$ only through the equivalence class $[E]$ of basis connected to $E$ through real transformations.
In this language one could say that a time-reversal operation is a choice of preferred equivalence class of basis. That is, at least, if $T^2 = 1$. If $T^2 = -1$, we should map $H$ to a $\mathbb{H}^{n/2}$, where $\mathbb{H}$ are the quaternions.
This leaves you with a recipe for computing the time-reversed operator, for $T^2=1$: simply represent $O$ in a real basis (a basis invariant under $T$), then take the complex conjugate of this matrix.
Notice that neither of these manipulations depend on $T$ being a symmetry.