Protein evolution has produced enzymes that maintain stability and function across various thermal environments. While sequence variation, structural dynamics, and intermolecular interactions are known to influence an enzyme's thermal adaptation, how these factors collectively govern stability and function across diverse temperatures remains unresolved. Cytosolic malate dehydrogenase (cMDH), a citric acid cycle enzyme, is an ideal model for studying these mechanisms due to its temperature-sensitive flexibility and broad presence in species from diverse thermal environments. In this study, we employ techniques inspired by deep learning and statistical mechanics to uncover how sequence variation and structural dynamics shape patterns of cMDH's thermal adaptation. By integrating coevolutionary models with variational autoencoders (VAE), we generate a latent generative landscape (LGL) of cMDH sequence space, enabling us to explore evolutionary pathways and predict fitness using direct coupling analysis (DCA). Structural predictions via AlphaFold and molecular dynamics simulations further illuminate how variations in hydrophobic interactions and conformational flexibility contribute to the thermal stability of warm- and cold-adapted cMDH orthologs. The integrative computational framework employed in this study provides powerful insights into protein adaptation at both sequence and structural levels, offering new perspectives on the evolution of thermal stability and creating avenues for the rational design of proteins with optimized thermal properties for biotechnological applications.