Abstract
The ruggedness of a protein’s fitness landscape constrains how it can evolve and how effectively it can be engineered. What determines whether a fitness landscape is rugged, and where in sequence space that ruggedness concentrates, remains poorly understood. Existing measures of ruggedness rely on complete combinatorial datasets and cannot be applied to landscapes in which only a fraction of possible sequences have been characterized. Here, we introduce a graph-based measure of ruggedness based on how quickly fitness differences diffuse across a landscape, which can be applied to sparse datasets. Applying this measure to 64 protein domains and several experimental systems with characterized folding and binding constraints, we show that ruggedness is not uniformly distributed but concentrates at regions where the dominant constraint on fitness changes. These results establish that the composition of individually simple constraints can generate localized ruggedness, providing a mechanistic basis for the complexity of protein fitness landscapes.
Full text
1,870 characters
· extracted from
oa-doi-fallback
· click to expand
Abstract
Epistasis, in which the effect of a mutation depends on the genetic background it is introduced into, drives protein evolution and design by shaping the fitness landscape. Quantifying how constrained by epistasis (and thus rugged) a fitness landscape is when data is combinatorially incomplete, however, remains challenging. Here, we propose a spectral graph theory approach to measure fitness landscape ruggedness as the characteristic time step of a heat diffusion process. Our approach is conceptually analogous to modelling a fitness signal as heat and measuring the ruggedness of the fitness landscape as the length of time that signal has diffused over a network graph for. Through theory and simulation, we show that our diffusion-based approach describes the global ruggedness of a fitness landscape, even when only a fraction of the fitness map is experimentally characterized. Applying this approach to empirical fitness landscapes that encompass deep mutational scans, combinatorially complete datasets, and sparse sampling of homologous sequences, we find large variation in ruggedness both among and within different protein folds and functions. Finally, we use spectral graph theory to demonstrate that ruggedness is intimately related to the fraction of sequence space that maintains sufficient fitness to overcome purifying selection, which we term the solution set of the fitness function. These findings offer a theoretically grounded explanation for why complex biophysical traits often exhibit strong epistasis, while simpler traits, such as fold stability, evolve over smooth fitness landscapes. Our results underscore how spectral methods can help unify and quantify the interplay between fitness constraints and evolutionary outcomes in protein sequence space.
Competing Interest Statement
The authors have declared no competing interest.
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.