We investigate the problem of heat conduction across molecular junctions connecting two nanoparticles, both in vacuum and in a liquid environment, using classical molecular dynamics simulations. In vacuum, the well-known result of a length independent conductance is recovered; its precise value, however, is found to depend sensitively on the overlap between the vibrational spectrum of the junction and the density of states of the nanoparticles that act as thermal contacts. In a liquid environment, the conductance is constant up to a crossover length, above which a standard Fourier regime is recovered.