We model the cooling of hybrid neutron stars combining a microscopic nuclear equation of state in the Brueckner-Hartree-Fock approach with different quark models. We then analyze the neutron star cooling curves predicted by the different models and single out the preferred ones. We find that the possibility of neutron p-wave pairing can be excluded in our scenario.