Exploring the sequence fitness landscape of a bridge between protein folds

Most foldable protein sequences adopt only a single native fold. Recent protein design studies have, however, created protein sequences which fold into different structures apon changes of environment, or single point mutation, the best characterized example being the switch between the folds of the GA and GB binding domains of streptococcal protein G. To obtain further insight into the design of sequences which can switch folds, we have used a computational model for the fitness landscape of a single fold, built from the observed sequence variation of protein homologues. We have recently shown that such coevolutionary models can be used to design novel foldable sequences. By appropriately combining two of these models to describe the joint fitness landscape of GA and GB, we are able to describe the propensity of a given sequence for each of the two folds. We have successfully tested the combined model against the known series of designed GA/GB hybrids. Using Monte Carlo simulations on this landscape, we are able to identify pathways of mutations connecting the two folds. In the absence of a requirement for domain stability, the most frequent paths go via sequences in which neither domain is stably folded, reminiscent of the propensity for certain intrinsically disordered proteins to fold into different structures according to context. Even if the folded state is required to be stable, we find that there is nonetheless still a wide range of sequences which are close to the transition region and therefore likely fold switches, consistent with recent estimates that fold switching may be more widespread than had been thought.Author summary: While most proteins self-assemble (or “fold”) to a unique three-dimensional structure, a few have been identified that can fold into two distinct structures. These so-called “metamorphic” proteins that can switch folds have attracted a lot of recent interest, and it has been suggested that they may be much more widespread than currently appreciated. We have developed a computational model that captures the propensity of a given protein sequence to fold into either one of two specific structures (GA and GB), in order to investigate which sequences are able to fold to both GA and GB (“switch sequences”), versus just one of them. Our model predicts that there is a large number of switch sequences that could fold into both structures, but also that the most likely such sequences are those for which the folded structures have low stability, in agreement with available experimental data. This also suggests that intrinsically disordered proteins which can fold into different structures on binding may provide an evolutionary path in sequence space between protein folds.