Pure quantum mechanics can be formulated as a Hamiltonian system in terms of the density matrix. Dissipative effects are modeled via coupling to a macroscopic system, where the coupling operators act via commutators. Following Öttinger (2010) we use the GENERIC framework (General Equations for Non-Equilibrium Reversible Irreversible Coupling) to construct thermodynamically consistent evolution equations as a sum of a Hamiltonian and a gradient-flow contribution, which satisfy a particular non-interaction condition. One of our models couples a quantum system to a finite number of heat baths each of which is described by a time-dependent temperature. The dissipation mechanism is modeled via the canonical correlation operator, which is the inverse of the Kubo-Mori metric for density matrices and which is strongly linked to the von Neumann entropy for quantum systems. Thus, one recovers the dissipative double-bracket operators of the Lindblad equations but encounters a correction term for the consistent coupling to the dissipative dynamics. For the finite-dimensional and isothermal case we provide a general existence result and discuss sufficient conditions that guarantee that all solutions converge to the unique thermal equilibrium state. Finally, we compare of our gradient flow formulation for quantum systems with the Wasserstein gradient flow formulation for the Fokker-Planck equation and the entropy gradient flow formulation for reversible Markov chains.