Modeling the coupled return-spread high frequency dynamics of large tick assets