Thermal state of transiently accreting neutron stars