I don't know a reference, but, as you've stated it, this is a trivial result:
If $p$ is a point on a smooth surface $S\subset\mathbb{R}^3$ at which the Gauss curvature is negative, then $p$ is non-umbilic, so, in an open $p$-neighborhood in $S$, both the principal curves and the asymptotic curves define transverse foliations of $S$. In particular, there exist $p$-centered local coordinates $(u,v):V\to\mathbb{R}^2$ such that the $u$-level curves and $v$-level curves are the principal curves and there exist $p$-centered local coordinates $(x,y):V\to\mathbb{R}^2$ such that the $x$-level curves and the $y$-level curves are the asymptotic curves. Hence there is an open $p$-neighborhood $U\subset V$ on which there exists a (unique) diffeomorphism $T:V\to U$ fixing $p$ such that $x = u\circ T$ and $y = v\circ T$. Since $T$ takes each $x$-level curve to a $u$-level curve and each $y$-level curve to a $v$-level curve, it carries every asymptotic curve in $V$ to a principal curve.
I suspect that you may be remembering a paper in which some specific method is given for constructing such a $T$ with more particular properties than just taking asymptotic curves to principal curves, but I confess that I don't remember seeing a discussion about such a particular construction in any of the classical literature that I have read.