Background: Heat-related mortality is an increasingly important public health burden that is expected to worsen with climate change. In addition to long-term trends, there are also interannual variations in heat-related mortality that are of interest for efficient planning of health services. Large-scale climate patterns have an important influence on summer weather and therefore constitute important tools to understand and predict the variations in heat-related mortality. Methods: In this article, we propose to model summer heat-related mortality using seven climate indices through a two-stage analysis using data covering the period 1981–2018 in two metropolitan areas of the province of Québec (Canada): Montréal and Québec. In the first stage, heat attributable fractions are estimated through a time series regression design and distributed lag nonlinear specification. We consider different definitions of heat. In the second stage, estimated attributable fractions are predicted using climate index curves through a functional linear regression model. Results: Results indicate that the Atlantic Multidecadal Oscillation is the best predictor of heat-related mortality in both Montréal and Québec and that it can predict up to 20% of the interannual variability. Conclusion: We found evidence that one climate index is predictive of summer heat-related mortality. More research is needed with longer time series and in different spatial contexts. The proposed analysis and the results may nonetheless help public health authorities plan for future mortality related to summer heat.