Understanding how species adapt to different temperatures is crucial to predict their response to global warming, and thermal performance curves (TPCs) have been employed recurrently to study this topic. Nevertheless, fundamental questions regarding how thermodynamic constraints and evolution interact to shape TPCs in lineages inhabiting different environments remain unanswered. Here, we study
Drosophila simulans
along a latitudinal gradient spanning 3000 km to test opposing hypotheses based on thermodynamic constrains (
hotter-is-better
) versus biochemical adaptation (
jack-of-all-temperatures
) as primary determinants of TPCs variation across populations. We compare thermal responses in metabolic rate and the egg-to-adult survival as descriptors of organismal performance and fitness, respectively, and show that different descriptors of TPCs vary in tandem with mean environmental temperatures, providing strong support to
hotter-is-better
. Thermodynamic constraints also resulted in a strong negative association between maximum performance and thermal breadth. Lastly, we show that descriptors of TPCs for metabolism and egg-to-adult survival are highly correlated, providing evidence of co-adaptation, and that curves for egg-to-adult survival are systematically narrower and displaced toward lower temperatures. Taken together, our results support the pervasive role of thermodynamics constraining thermal responses in
Drosophila
populations along a latitudinal gradient, that are only partly compensated by evolutionary adaptation.