Abstract : Recent tensor network techniques for simulating system-environment wavefunctions have provided profound insights into non-Markovian dissipation and decoherence in open quantum systems. Here, we propose a dynamically adaptive one-site Time-Dependent-Variational-Principle (A1TDVP) method for matrix product states in which local bond dimensions grow to capture developing systembath entanglement. This avoids the need for multiple convergence runs w.r.t. bond dimensions and the unfavourable local Hilbert space scaling of two-site methods. A1TDVP is thus ideally suited for open quantum dynamics in finite-temperature bosonic environments, as the initial states typically have low bond dimension but require very large local physical dimensions. We demonstrate this with simulations of non-equilibrium heat flows through a qubit spin, finding a 30x and 10x speed-up over 2TDVP and 1TDVP, respectively.