We investigate a dissipatively driven XYZ spin-1/2 chain in the Zeno limit of strong dissipation, described by Lindblad master equation. The nonequilibrium steady state is expressed in terms of a matrix product ansatz using novel site-dependent Lax operators. The components of Lax operators satisfy a simple set of linear recurrence equations that generalize the defining algebraic relations of the quantum group Uq(sl2). We reveal connection between the nonequilibrium steady state of the nonunitary dynamics and the respective integrable model with edge magnetic fields, described by coherent unitary dynamics.Introduction.-One of the main current efforts of the condensed matter physics society is to understand quantum states of matter far from equilibrium [1,2]. Understandably, simple models with tractable but nontrivial exact solutions are of key importance in this game. The realm of driven dissipative quantum many-body systems [3] provides nice and rich examples of such models, capable of displaying genuinely out-of-equilibrium phenomena, for example, novel types of non-equilibrium phase transitions [4][5][6][7]. While exact treatment of the aforementioned class of models is essentially limited to quasi-free situations, it is remarkable that some exact solutions have been found even in the case of strong interactions, in particular, in quantum integrable spin chains with dissipation and incoherent driving localized at the chain's boundaries [8]. Despite the fact, that the exact matrix product form of these solutions has been found only for a very specific choice of the boundary jump operators [9, 10], this provided a fresh perspective on the effect of local and quasilocal conservation laws on quantum transport and relaxation [11]. It has, however, remained an open question of how these exact steady state solutions fit into the general framework of integrability. For example, except in the special case of dissipatively driven noninteracting models [12], the solvable dissipatively driven boundaries cannot be generated using the solutions of the ubiquitous reflection equations [13], which constitute the standard framework for generating integrable boundaries in the coherent (nondissipative, Hamiltonian) setting.In the present letter we make a step forward in the understanding of integrability of open XYZ spin-1/2 chain, dissipatively driven at the boundaries. We present an exact solution for the nonequilibrium steady state of the model with arbitrary boundary processes, as long as the edge spins are described by pure states in the large dissipation limit. The inhomogeneous matrix product ansatz for the steady state differs drastically from the previously known solutions that describe steady states of the spin chains driven along a particular axis at both edges.Firstly, we introduce the model and describe how the