Time-dependent density-matrix functional theory for biexcitonic phenomena