Abstract
Proteins can be studied through their sequence statistics or structural properties. These represent complementary views that are useful but lack a quantitative framework to tell, family by family, which is most informative and how to combine them. We introduce a framework that builds both views in parallel: amino acid (AA) alignments are translated into parallel alignments over a 3D interaction (3Di) structure-informed alphabet. Variational autoencoders compress each into a two-dimensional map, and direct coupling analysis places a shared coevolutionary energy on both maps, turning them into latent generative landscapes. On these landscapes, we define information-theoretic distance metrics that quantify how sequence changes drive structural and functional variation in protein families. We demonstrate the framework on five families: in malate dehydrogenases, the 3Di landscape identifies the structurally conserved scaffold that this family uses to encode thermal adaptation via sequence variability revealed in the AA landscape. In globins and transient receptor potential melastatin (TRPM), the 3Di landscape recovers known functional subfamilies. In the Flaviviridae E1 and E2 glycoproteins, structure reveals evolutionary relationships invisible at the sequence level. Because many sequences encode the same fold, our framework lets us disentangle family-sequence variability from structural and functional variation. These generative landscapes allow sampling near functional regions, and we show they can help us gain mechanistic insight into the evolutionary forces shaping sequence–structure–function variation and guide the design of new proteins.