Investigating the heat capacity of metal-organic frameworks via path integral molecular dynamics simulations