Extant fold-switching proteins are widespread

Data-driven probabilistic definition of the low energy conformational states of protein residues

Protein dynamics and related conformational changes are essential for their function but difficult to characterise and interpret. Amino acids in a protein behave according to their local energy landscape, which is determined by their local structural context and environmental conditions. The lowest energy state for a given residue can correspond to sharply defined conformations, e.g. in a stable helix, or can cover a wide range of conformations, e.g. in intrinsically disordered regions. A good definition of such low energy states is therefore important to describe the behaviour of a residue and how it changes with its environment. We propose a data-driven probabilistic definition of six low energy conformational states typically accessible for amino acid residues in proteins. This definition is based on solution NMR information of 1322 proteins through a combined analysis of structure ensembles with interpreted chemical shifts. We further introduce a conformational state variability parameter that captures, based on an ensemble of protein structures from molecular dynamics or other methods, how often a residue moves between these conformational states. The approach enables a different perspective on the local conformational behaviour of proteins that is complementary to their static interpretation from single structure models.

Metamorphic proteins and how to find them

In the last two decades, our existing notion that most foldable proteins have a unique native state has been challenged by the discovery of metamorphic proteins, which reversibly interconvert between multiple, sometimes highly dissimilar, native states. As the number of known metamorphic proteins increases, several computational and experimental strategies have emerged for gaining insights about their refolding processes and identifying unknown metamorphic proteins amongst the known proteome. In this review, we describe the current advances in biophysically and functionally ascertaining the structural interconversions of metamorphic proteins and how coevolution can be harnessed to identify novel metamorphic proteins from sequence information. We also discuss the challenges and ongoing efforts in using artificial intelligence-based protein structure prediction methods to discover metamorphic proteins and predict their corresponding three-dimensional structures.

Heritable epigenetic changes are constrained by the dynamics of regulatory architectures

Interacting molecules create regulatory architectures that can persist despite turnover of molecules. Although epigenetic changes occur within the context of such architectures, there is limited understanding of how they can influence the heritability of changes. Here, I develop criteria for the heritability of regulatory architectures and use quantitative simulations of interacting regulators parsed as entities, their sensors, and the sensed properties to analyze how architectures influence heritable epigenetic changes. Information contained in regulatory architectures grows rapidly with the number of interacting molecules and its transmission requires positive feedback loops. While these architectures can recover after many epigenetic perturbations, some resulting changes can become permanently heritable. Architectures that are otherwise unstable can become heritable through periodic interactions with external regulators, which suggests that mortal somatic lineages with cells that reproducibly interact with the immortal germ lineage could make a wider variety of architectures heritable. Differential inhibition of the positive feedback loops that transmit regulatory architectures across generations can explain the gene-specific differences in heritable RNA silencing observed in the nematode Caenorhabditis elegans. More broadly, these results provide a foundation for analyzing the inheritance of epigenetic changes within the context of the regulatory architectures implemented using diverse molecules in different living systems.

Identifying protein conformational states in the Protein Data Bank: Toward unlocking the potential of integrative dynamics studies

Studying protein dynamics and conformational heterogeneity is crucial for understanding biomolecular systems and treating disease. Despite the deposition of over 215 000 macromolecular structures in the Protein Data Bank and the advent of AI-based structure prediction tools such as AlphaFold2, RoseTTAFold, and ESMFold, static representations are typically produced, which fail to fully capture macromolecular motion. Here, we discuss the importance of integrating experimental structures with computational clustering to explore the conformational landscapes that manifest protein function. We describe the method developed by the Protein Data Bank in Europe - Knowledge Base to identify distinct conformational states, demonstrate the resource's primary use cases, through examples, and discuss the need for further efforts to annotate protein conformations with functional information. Such initiatives will be crucial in unlocking the potential of protein dynamics data, expediting drug discovery research, and deepening our understanding of macromolecular mechanisms.

Can Protein Structure Prediction Methods Capture Alternative Conformations of Membrane Transporters?

Understanding the conformational dynamics of proteins, such as the inward-facing (IF) and outward-facing (OF) transition observed in transporters, is vital for elucidating their functional mechanisms. Despite significant advances in protein structure prediction (PSP) over the past three decades, most efforts have been focused on single-state prediction, leaving multistate or alternative conformation prediction (ACP) relatively unexplored. This discrepancy has led to the development of highly accurate PSP methods such as AlphaFold, yet their capabilities for ACP remain limited. To investigate the performance of current PSP methods in ACP, we curated a data set, named IOMemP, consisting of 32 experimentally determined high-resolution IF and OF structures of 16 membrane proteins with substantial conformational changes. We benchmarked 12 representative PSP methods, along with two recent multistate methods based on AlphaFold, against this data set. Our findings reveal a remarkably consistent preference for specific states across various PSP methods. We elucidated how coevolution information in MSAs influences state preference. Moreover, we showed that AlphaFold, when excluding coevolution information, estimated similar energies between the experimental IF and OF conformations, indicating that the energy model learned by AlphaFold is not biased toward any particular state. Our IOMemP data set and benchmark results are anticipated to advance the development of robust ACP methods.

High-throughput prediction of protein conformational distributions with subsampled AlphaFold2

This paper presents an innovative approach for predicting the relative populations of protein conformations using AlphaFold 2, an AI-powered method that has revolutionized biology by enabling the accurate prediction of protein structures. While AlphaFold 2 has shown exceptional accuracy and speed, it is designed to predict proteins' ground state conformations and is limited in its ability to predict conformational landscapes. Here, we demonstrate how AlphaFold 2 can directly predict the relative populations of different protein conformations by subsampling multiple sequence alignments. We tested our method against nuclear magnetic resonance experiments on two proteins with drastically different amounts of available sequence data, Abl1 kinase and the granulocyte-macrophage colony-stimulating factor, and predicted changes in their relative state populations with more than 80% accuracy. Our subsampling approach worked best when used to qualitatively predict the effects of mutations or evolution on the conformational landscape and well-populated states of proteins. It thus offers a fast and cost-effective way to predict the relative populations of protein conformations at even single-point mutation resolution, making it a useful tool for pharmacology, analysis of experimental results, and predicting evolution.

Proteostatic tuning underpins the evolution of novel multicellular traits

The evolution of multicellularity paved the way for the origin of complex life on Earth, but little is known about the mechanistic basis of early multicellular evolution. Here, we examine the molecular basis of multicellular adaptation in the multicellularity long-term evolution experiment (MuLTEE). We demonstrate that cellular elongation, a key adaptation underpinning increased biophysical toughness and organismal size, is convergently driven by down-regulation of the chaperone Hsp90. Mechanistically, Hsp90-mediated morphogenesis operates by destabilizing the cyclin-dependent kinase Cdc28, resulting in delayed mitosis and prolonged polarized growth. Reinstatement of Hsp90 or Cdc28 expression resulted in shortened cells that formed smaller groups with reduced multicellular fitness. Together, our results show how ancient protein folding systems can be tuned to drive rapid evolution at a new level of biological individuality by revealing novel developmental phenotypes.

The Cryptic Bacterial Microproteome

Microproteins encoded by small open reading frames (smORFs) comprise the "dark matter" of proteomes. Although functional microproteins were identified in diverse organisms from all three domains of life, bacterial smORFs remain poorly characterized. In this comprehensive study of intergenic smORFs (ismORFs, 15-70 codons) in 5,668 bacterial genomes of the family Enterobacteriaceae, we identified 67,297 clusters of ismORFs subject to purifying selection. The ismORFs mainly code for hydrophobic, potentially transmembrane, unstructured, or minimally structured microproteins. Using AlphaFold Multimer, we predicted interactions of some of the predicted microproteins encoded by transcribed ismORFs with proteins encoded by neighboring genes, revealing the potential of microproteins to regulate the activity of various proteins, particularly, under stress. We compiled a catalog of predicted microprotein families with different levels of evidence from synteny analysis, structure prediction, and transcription and translation data. This study offers a resource for investigation of biological functions of microproteins.

Proteostatic tuning underpins the evolution of novel multicellular traits

The evolution of multicellularity paved the way for the origin of complex life on Earth, but little is known about the mechanistic basis of early multicellular evolution. Here, we examine the molecular basis of multicellular adaptation in the Multicellularity Long Term Evolution Experiment (MuLTEE). We demonstrate that cellular elongation, a key adaptation underpinning increased biophysical toughness and organismal size, is convergently driven by downregulation of the chaperone Hsp90. Mechanistically, Hsp90-mediated morphogenesis operates by destabilizing the cyclin-dependent kinase Cdc28, resulting in delayed mitosis and prolonged polarized growth. Reinstatement of Hsp90 or Cdc28 expression resulted in shortened cells that formed smaller groups with reduced multicellular fitness. Together, our results show how ancient protein folding systems can be tuned to drive rapid evolution at a new level of biological individuality by revealing novel developmental phenotypes.

Sequence clustering confounds AlphaFold2

Though typically associated with a single folded state, some globular proteins remodel their secondary and/or tertiary structures in response to cellular stimuli. AlphaFold21 (AF2) readily generates one dominant protein structure for these fold-switching (a.k.a. metamorphic) proteins2, but it often fails to predict their alternative experimentally observed structures3,4. Wayment-Steele, et al. steered AF2 to predict alternative structures of a few metamorphic proteins using a method they call AF-cluster5. However, their Paper lacks some essential controls needed to assess AF-cluster's reliability. We find that these controls show AF-cluster to be a poor predictor of metamorphic proteins. First, closer examination of the Paper's results reveals that random sequence sampling outperforms sequence clustering, challenging the claim that AF-cluster works by "deconvolving conflicting sets of couplings." Further, we observe that AF-cluster mistakes some single-folding KaiB homologs for fold switchers, a critical flaw bound to mislead users. Finally, proper error analysis reveals that AF-cluster predicts many correct structures with low confidence and some experimentally unobserved conformations with confidences similar to experimentally observed ones. For these reasons, we suggest using ColabFold6-based random sequence sampling7-augmented by other predictive approaches-as a more accurate and less computationally intense alternative to AF-cluster.

Predicting multiple conformations via sequence clustering and AlphaFold2

AlphaFold2 (ref. 1) has revolutionized structural biology by accurately predicting single structures of proteins. However, a protein's biological function often depends on multiple conformational substates2, and disease-causing point mutations often cause population changes within these substates3,4. We demonstrate that clustering a multiple-sequence alignment by sequence similarity enables AlphaFold2 to sample alternative states of known metamorphic proteins with high confidence. Using this method, named AF-Cluster, we investigated the evolutionary distribution of predicted structures for the metamorphic protein KaiB5 and found that predictions of both conformations were distributed in clusters across the KaiB family. We used nuclear magnetic resonance spectroscopy to confirm an AF-Cluster prediction: a cyanobacteria KaiB variant is stabilized in the opposite state compared with the more widely studied variant. To test AF-Cluster's sensitivity to point mutations, we designed and experimentally verified a set of three mutations predicted to flip KaiB from Rhodobacter sphaeroides from the ground to the fold-switched state. Finally, screening for alternative states in protein families without known fold switching identified a putative alternative state for the oxidoreductase Mpt53 in Mycobacterium tuberculosis. Further development of such bioinformatic methods in tandem with experiments will probably have a considerable impact on predicting protein energy landscapes, essential for illuminating biological function.

Predicting Relative Populations of Protein Conformations without a Physics Engine Using AlphaFold 2

This paper presents a novel approach for predicting the relative populations of protein conformations using AlphaFold 2, an AI-powered method that has revolutionized biology by enabling the accurate prediction of protein structures. While AlphaFold 2 has shown exceptional accuracy and speed, it is designed to predict proteins' ground state conformations and is limited in its ability to predict conformational landscapes. Here, we demonstrate how AlphaFold 2 can directly predict the relative populations of different protein conformations by subsampling multiple sequence alignments. We tested our method against NMR experiments on two proteins with drastically different amounts of available sequence data, Abl1 kinase and the granulocyte-macrophage colony-stimulating factor, and predicted changes in their relative state populations with more than 80% accuracy. Our subsampling approach worked best when used to qualitatively predict the effects of mutations or evolution on the conformational landscape and well-populated states of proteins. It thus offers a fast and cost-effective way to predict the relative populations of protein conformations at even single-point mutation resolution, making it a useful tool for pharmacology, NMR analysis, and evolution.

AlphaFold2 has more to learn about protein energy landscapes

Recent work suggests that AlphaFold2 (AF2)-a deep learning-based model that can accurately infer protein structure from sequence-may discern important features of folded protein energy landscapes, defined by the diversity and frequency of different conformations in the folded state. Here, we test the limits of its predictive power on fold-switching proteins, which assume two structures with regions of distinct secondary and/or tertiary structure. Using several implementations of AF2, including two published enhanced sampling approaches, we generated >280,000 models of 93 fold-switching proteins whose experimentally determined conformations were likely in AF2's training set. Combining all models, AF2 predicted fold switching with a modest success rate of ~25%, indicating that it does not readily sample both experimentally characterized conformations of most fold switchers. Further, AF2's confidence metrics selected against models consistent with experimentally determined fold-switching conformations in favor of inconsistent models. Accordingly, these confidence metrics-though suggested to evaluate protein energetics reliably-did not discriminate between low and high energy states of fold-switching proteins. We then evaluated AF2's performance on seven fold-switching proteins outside of its training set, generating >159,000 models in total. Fold switching was accurately predicted in one of seven targets with moderate confidence. Further, AF2 demonstrated no ability to predict alternative conformations of two newly discovered targets without homologs in the set of 93 fold switchers. These results indicate that AF2 has more to learn about the underlying energetics of protein ensembles and highlight the need for further developments of methods that readily predict multiple protein conformations.

Heritable epigenetic changes are constrained by the dynamics of regulatory architectures

Interacting molecules create regulatory architectures that can persist despite turnover of molecules. Although epigenetic changes occur within the context of such architectures, there is limited understanding of how they can influence the heritability of changes. Here I develop criteria for the heritability of regulatory architectures and use quantitative simulations of interacting regulators parsed as entities, their sensors and the sensed properties to analyze how architectures influence heritable epigenetic changes. Information contained in regulatory architectures grows rapidly with the number of interacting molecules and its transmission requires positive feedback loops. While these architectures can recover after many epigenetic perturbations, some resulting changes can become permanently heritable. Such stable changes can (1) alter steady-state levels while preserving the architecture, (2) induce different architectures that persist for many generations, or (3) collapse the entire architecture. Architectures that are otherwise unstable can become heritable through periodic interactions with external regulators, which suggests that the evolution of mortal somatic lineages with cells that reproducibly interact with the immortal germ lineage could make a wider variety of regulatory architectures heritable. Differential inhibition of the positive feedback loops that transmit regulatory architectures across generations can explain the gene-specific differences in heritable RNA silencing observed in the nematode C. elegans, which range from permanent silencing to recovery from silencing within a few generations and subsequent resistance to silencing. More broadly, these results provide a foundation for analyzing the inheritance of epigenetic changes within the context of the regulatory architectures implemented using diverse molecules in different living systems.

SSDraw: Software for generating comparative protein secondary structure diagrams

The program SSDraw generates publication-quality protein secondary structure diagrams from three-dimensional protein structures. To depict relationships between secondary structure and other protein features, diagrams can be colored by conservation score, B-factor, or custom scoring. Diagrams of homologous proteins can be registered according to an input multiple sequence alignment. Linear visualization allows the user to stack registered diagrams, facilitating comparison of secondary structure and other properties among homologous proteins. SSDraw can be used to compare secondary structures of homologous proteins with both conserved and divergent folds. It can also generate one secondary structure diagram from an input protein structure of interest. The source code can be downloaded (https://github.com/ncbi/SSDraw) and run locally for rapid structure generation, while a Google Colab notebook allows easy use.

ColabFold predicts alternative protein structures from single sequences, coevolution unnecessary for AF-cluster

Though typically associated with a single folded state, globular proteins are dynamic and often assume alternative or transient structures important for their functions1,2. Wayment-Steele, et al. steered ColabFold3 to predict alternative structures of several proteins using a method they call AF-cluster4. They propose that AF-cluster "enables ColabFold to sample alternate states of known metamorphic proteins with high confidence" by first clustering multiple sequence alignments (MSAs) in a way that "deconvolves" coevolutionary information specific to different conformations and then using these clusters as input for ColabFold. Contrary to this Coevolution Assumption, clustered MSAs are not needed to make these predictions. Rather, these alternative structures can be predicted from single sequences and/or sequence similarity, indicating that coevolutionary information is unnecessary for predictive success and may not be used at all. These results suggest that AF-cluster's predictive scope is likely limited to sequences with distinct-yet-homologous structures within ColabFold's training set.

The Boltzmann distributions of molecular structures predict likely changes through random mutations

New folded molecular structures can only evolve after arising through mutations. This aspect is modeled using genotype-phenotype maps, which connect sequence changes through mutations to changes in molecular structures. Previous work has shown that the likelihood of appearing through mutations can differ by orders of magnitude from structure to structure and that this can affect the outcomes of evolutionary processes. Thus, we focus on the phenotypic mutation probabilities φqp, i.e., the likelihood that a random mutation changes structure p into structure q. For both RNA secondary structures and the HP protein model, we show that a simple biophysical principle can explain and predict how this likelihood depends on the new structure q: φqp is high if sequences that fold into p as the minimum-free-energy structure are likely to have q as an alternative structure with high Boltzmann frequency. This generalizes the existing concept of plastogenetic congruence from individual sequences to the entire neutral spaces of structures. Our result helps us understand why some structural changes are more likely than others, may be useful for estimating these likelihoods via sampling and makes a connection to alternative structures with high Boltzmann frequency, which could be relevant in evolutionary processes.

SSDraw: software for generating comparative protein secondary structure diagrams

The program SSDraw generates publication-quality protein secondary structure diagrams from three-dimensional protein structures. To depict relationships between secondary structure and other protein features, diagrams can be colored by conservation score, B-factor, or custom scoring. Diagrams of homologous proteins can be registered according to an input multiple sequence alignment. Linear visualization allows the user to stack registered diagrams, facilitating comparison of secondary structure and other properties among homologous proteins. SSDraw can be used to compare secondary structures of homologous proteins with both conserved and divergent folds. It can also generate one secondary structure diagram from an input protein structure of interest. The source code can be downloaded (https://github.com/ethanchen1301/SSDraw) and run locally for rapid structure generation, while a Google Colab notebook allows easy use.

Global analysis of aging-related protein structural changes uncovers enzyme-polymerization-based control of longevity

Aging is associated with progressive phenotypic changes. Virtually all cellular phenotypes are produced by proteins, and their structural alterations can lead to age-related diseases. However, we still lack comprehensive knowledge of proteins undergoing structural-functional changes during cellular aging and their contributions to age-related phenotypes. Here, we conducted proteome-wide analysis of early age-related protein structural changes in budding yeast using limited proteolysis-mass spectrometry (LiP-MS). The results, compiled in online ProtAge catalog, unraveled age-related functional changes in regulators of translation, protein folding, and amino acid metabolism. Mechanistically, we found that folded glutamate synthase Glt1 polymerizes into supramolecular self-assemblies during aging, causing breakdown of cellular amino acid homeostasis. Inhibiting Glt1 polymerization by mutating the polymerization interface restored amino acid levels in aged cells, attenuated mitochondrial dysfunction, and led to lifespan extension. Altogether, this comprehensive map of protein structural changes enables identifying mechanisms of age-related phenotypes and offers opportunities for their reversal.

Simulations of a protein fold switch reveal crowding-induced population shifts driven by disordered regions

Macromolecular crowding effects on globular proteins, which usually adopt a single stable fold, have been widely studied. However, little is known about crowding effects on fold-switching proteins, which reversibly switch between distinct folds. Here we study the mutationally driven switch between the folds of GA and GB, the two 56-amino acid binding domains of protein G, using a structure-based dual-basin model. We show that, in the absence of crowders, the fold populations PA and PB can be controlled by the strengths of contacts in the two folds, κA and κB. A population balance, PA ≈ PB, is obtained for κB/κA = 0.92. The resulting model protein is subject to crowding at different packing fractions, ϕc. We find that crowding increases the GB population and reduces the GA population, reaching PB/PA ≈ 4 at ϕc = 0.44. We analyze the ϕc-dependence of the crowding-induced GA-to-GB switch using scaled particle theory, which provides a qualitative, but not quantitative, fit of our data, suggesting effects beyond a spherical description of the folds. We show that the terminal regions of the protein chain, which are intrinsically disordered only in GA, play a dominant role in the response of the fold switch to crowding effects.

Evolutionary selection of proteins with two folds

Although most globular proteins fold into a single stable structure, an increasing number have been shown to remodel their secondary and tertiary structures in response to cellular stimuli. State-of-the-art algorithms predict that these fold-switching proteins adopt only one stable structure, missing their functionally critical alternative folds. Why these algorithms predict a single fold is unclear, but all of them infer protein structure from coevolved amino acid pairs. Here, we hypothesize that coevolutionary signatures are being missed. Suspecting that single-fold variants could be masking these signatures, we developed an approach, called Alternative Contact Enhancement (ACE), to search both highly diverse protein superfamilies-composed of single-fold and fold-switching variants-and protein subfamilies with more fold-switching variants. ACE successfully revealed coevolution of amino acid pairs uniquely corresponding to both conformations of 56/56 fold-switching proteins from distinct families. Then, we used ACE-derived contacts to (1) predict two experimentally consistent conformations of a candidate protein with unsolved structure and (2) develop a blind prediction pipeline for fold-switching proteins. The discovery of widespread dual-fold coevolution indicates that fold-switching sequences have been preserved by natural selection, implying that their functionalities provide evolutionary advantage and paving the way for predictions of diverse protein structures from single sequences.

Fluid protein fold space and its implications

Fold-switching proteins, which remodel their secondary and tertiary structures in response to cellular stimuli, suggest a new view of protein fold space. For decades, experimental evidence has indicated that protein fold space is discrete: dissimilar folds are encoded by dissimilar amino acid sequences. Challenging this assumption, fold-switching proteins interconnect discrete groups of dissimilar protein folds, making protein fold space fluid. Three recent observations support the concept of fluid fold space: (1) some amino acid sequences interconvert between folds with distinct secondary structures, (2) some naturally occurring sequences have switched folds by stepwise mutation, and (3) fold switching is evolutionarily selected and likely confers advantage. These observations indicate that minor amino acid sequence modifications can transform protein structure and function. Consequently, proteomic structural and functional diversity may be expanded by alternative splicing, small nucleotide polymorphisms, post-translational modifications, and modified translation rates.

A geometrical framework for thinking about proteins

We present a model, based on symmetry and geometry, for proteins. Using elementary ideas from mathematics and physics, we derive the geometries of discrete helices and sheets. We postulate a compatible solvent-mediated emergent pairwise attraction that assembles these building blocks, while respecting their individual symmetries. Instead of seeking to mimic the complexity of proteins, we look for a simple abstraction of reality that yet captures the essence of proteins. We employ analytic calculations and detailed Monte Carlo simulations to explore some consequences of our theory. The predictions of our approach are in accord with experimental data. Our framework provides a rationalization for understanding the common characteristics of proteins. Our results show that the free energy landscape of a globular protein is pre-sculpted at the backbone level, sequences and functionalities evolve in the fixed backdrop of the folds determined by geometry and symmetry, and that protein structures are unique in being simultaneously characterized by stability, diversity, and sensitivity.

The non-deterministic genotype-phenotype map of RNA secondary structure

Selection and variation are both key aspects in the evolutionary process. Previous research on the mapping between molecular sequence (genotype) and molecular fold (phenotype) has shown the presence of several structural properties in different biological contexts, implying that these might be universal in evolutionary spaces. The deterministic genotype-phenotype (GP) map that links short RNA sequences to minimum free energy secondary structures has been studied extensively because of its computational tractability and biologically realistic nature. However, this mapping ignores the phenotypic plasticity of RNA. We define a GP map that incorporates non-deterministic (ND) phenotypes, and take RNA as a case study; we use the Boltzmann probability distribution of folded structures and examine the structural properties of ND GP maps for RNA sequences of length 12 and coarse-grained RNA structures of length 30 (RNAshapes30). A framework is presented to study robustness, evolvability and neutral spaces in the ND map. This framework is validated by demonstrating close correspondence between the ND quantities and sample averages of their deterministic counterparts. When using the ND framework we observe the same structural properties as in the deterministic GP map, such as bias, negative correlation between genotypic robustness and evolvability, and positive correlation between phenotypic robustness and evolvability.

Metamorphic protein folding as evolutionary adaptation

Metamorphic proteins switch reversibly between multiple distinct, stable structures, often with different functions. It was previously hypothesized that metamorphic proteins arose as intermediates in the evolution of a new fold - rare and transient exceptions to the 'one sequence, one fold' paradigm. However, as described herein, mounting evidence suggests that metamorphic folding is an adaptive feature, preserved and optimized over evolutionary time as exemplified by the NusG family and the chemokine XCL1. Analysis of extant protein families and resurrected protein ancestors demonstrates that large regions of sequence space are compatible with metamorphic folding. As a category that enhances biological fitness, metamorphic proteins are likely to employ fold switching to perform important biological functions and may be more common than previously thought.

Current successes and remaining challenges in protein function prediction

In recent years, improvements in protein function prediction methods have led to increased success in annotating protein sequences. However, the functions of over 30% of protein-coding genes remain unknown for many sequenced genomes. Protein functions vary widely, from catalyzing chemical reactions to binding DNA or RNA or forming structures in the cell, and some types of functions are challenging to predict due to the physical features associated with those functions. Other complications in understanding protein functions arise due to the fact that many proteins have more than one function or very small differences in sequence or structure that correspond to different functions. We will discuss some of the recent developments in predicting protein functions and some of the remaining challenges.

Intrinsically disordered proteins and conformational noise: The hypothesis a decade later

Phenotypic plasticity is the ability of individual genotypes to produce different phenotypes in response to environmental perturbations. We previously postulated how conformational noise emanating from conformational dynamics of intrinsically disordered proteins (IDPs) which is distinct from transcriptional noise, can contribute to phenotypic switching by rewiring the cellular protein interaction network. Since most transcription factors are IDPs, we posited that conformational noise is an integral component of transcriptional noise implying that IDPs may amplify total noise in the system either stochastically or in response to environmental changes. Here, we review progress in elucidating the details of the hypothesis. We highlight empirical evidence supporting the hypothesis, discuss conceptual advances that underscore its fundamental importance and implications, and identify areas for future investigations.

Crowded environments tune the fold-switching in metamorphic proteins

Metamorphic proteins such as circadian clock protein KaiB and human chemokine XCL1 play vital roles in regulating biological processes, including gene expression, circadian clock and innate immune responses, and perform distinct functions in living cell by switching different structures in response to cellular environment stimuli. However, it is unclear how complex and crowded intracellular environments affect conformational rearrangement of metamorphic proteins. Here, the kinetics and thermodynamics of two well-characterized metamorphic proteins, circadian clock protein KaiB and human chemokine XCL1, were quantified in physiologically relevant environments by using NMR spectroscopy, indicating that crowded agents shift equilibrium towards the inactive form (ground-state KaiB and Ltn10-like state XCL1) without disturbing the corresponding structures, and crowded agents have predominantly impact on the exchange rate of XCL1 that switches folds on timescales of seconds, but have slightly impact on the exchange rate of KaiB that switches folds on timescales of hours. Our data shed light on how metamorphic proteins can respond immediately to the changed crowded intracellular conditions that induced by environmental cues and then execute different functions in living cell, and it also enhances our understanding of how environments enrich the sequence-structure-function paradigm.

Identification of a covert evolutionary pathway between two protein folds

Although homologous protein sequences are expected to adopt similar structures, some amino acid substitutions can interconvert α-helices and β-sheets. Such fold switching may have occurred over evolutionary history, but supporting evidence has been limited by the: (1) abundance and diversity of sequenced genes, (2) quantity of experimentally determined protein structures, and (3) assumptions underlying the statistical methods used to infer homology. Here, we overcome these barriers by applying multiple statistical methods to a family of ~600,000 bacterial response regulator proteins. We find that their homologous DNA-binding subunits assume divergent structures: helix-turn-helix versus α-helix + β-sheet (winged helix). Phylogenetic analyses, ancestral sequence reconstruction, and AlphaFold2 models indicate that amino acid substitutions facilitated a switch from helix-turn-helix into winged helix. This structural transformation likely expanded DNA-binding specificity. Our approach uncovers an evolutionary pathway between two protein folds and provides a methodology to identify secondary structure switching in other protein families.

Systematic conformation-to-phenotype mapping via limited deep sequencing of proteins

Non-native conformations drive protein-misfolding diseases, complicate bioengineering efforts, and fuel molecular evolution. No current experimental technique is well suited for elucidating them and their phenotypic effects. Especially intractable are the transient conformations populated by intrinsically disordered proteins. We describe an approach to systematically discover, stabilize, and purify native and non-native conformations, generated in vitro or in vivo, and directly link conformations to molecular, organismal, or evolutionary phenotypes. This approach involves high-throughput disulfide scanning (HTDS) of the entire protein. To reveal which disulfides trap which chromatographically resolvable conformers, we devised a deep-sequencing method for double-Cys variant libraries of proteins that precisely and simultaneously locates both Cys residues within each polypeptide. HTDS of the abundant E. coli periplasmic chaperone HdeA revealed distinct classes of disordered hydrophobic conformers with variable cytotoxicity depending on where the backbone was cross-linked. HTDS can bridge conformational and phenotypic landscapes for many proteins that function in disulfide-permissive environments.

Moonlighting enzymes: when cellular context defines specificity

It is not often realized that the absolute protein specificity is an exception rather than a rule. Two major kinds of protein multi-specificities are promiscuity and moonlighting. This review discusses the idea of enzyme specificity and then focusses on moonlighting. Some important examples of protein moonlighting, such as crystallins, ceruloplasmin, metallothioniens, macrophage migration inhibitory factor, and enzymes of carbohydrate metabolism are discussed. How protein plasticity and intrinsic disorder enable the removing the distinction between enzymes and other biologically active proteins are outlined. Finally, information on important roles of moonlighting in human diseases is updated.

Distinguishing features of fold-switching proteins

Though many folded proteins assume one stable structure that performs one function, a small-but-increasing number remodel their secondary and tertiary structures and change their functions in response to cellular stimuli. These fold-switching proteins regulate biological processes and are associated with autoimmune dysfunction, severe acute respiratory syndrome coronavirus-2 infection, and more. Despite their biological importance, it is difficult to computationally predict fold switching. With the aim of advancing computational prediction and experimental characterization of fold switchers, this review discusses several features that distinguish fold-switching proteins from their single-fold and intrinsically disordered counterparts. First, the isolated structures of fold switchers are less stable and more heterogeneous than single folders but more stable and less heterogeneous than intrinsically disordered proteins (IDPs). Second, the sequences of single fold, fold switching, and intrinsically disordered proteins can evolve at distinct rates. Third, proteins from these three classes are best predicted using different computational techniques. Finally, late-breaking results suggest that single folders, fold switchers, and IDPs have distinct patterns of residue-residue coevolution. The review closes by discussing high-throughput and medium-throughput experimental approaches that might be used to identify new fold-switching proteins.

Genetically encoded non-canonical amino acids reveal asynchronous dark reversion of chromophore, backbone and side-chains in EL222

Photoreceptors containing the light-oxygen-voltage (LOV) domain elicit biological responses upon excitation of their flavin mononucleotide (FMN) chromophore by blue light. The mechanism and kinetics of dark-state recovery are not well understood. Here we incorporated the non-canonical amino acid p-cyanophenylalanine (CNF) by genetic code expansion technology at forty-five positions of the bacterial transcription factor EL222. Screening of light-induced changes in infrared (IR) absorption frequency, electric field and hydration of the nitrile groups identified residues CNF31 and CNF35 as reporters of monomer/oligomer and caged/decaged equilibria, respectively. Time-resolved multi-probe UV/Visible and IR spectroscopy experiments of the lit-to-dark transition revealed four dynamical events. Predominantly, rearrangements around the A'α helix interface (CNF31 and CNF35) precede FMN-cysteinyl adduct scission, folding of α-helices (amide bands), and relaxation of residue CNF151. This study illustrates the importance of characterizing all parts of a protein and suggests a key role for the N-terminal A'α extension of the LOV domain in controlling EL222 photocycle length. This article is protected by copyright. All rights reserved.

A review of visualisations of protein fold networks and their relationship with sequence and function

Proteins form arguably the most significant link between genotype and phenotype. Understanding the relationship between protein sequence and structure, and applying this knowledge to predict function, is difficult. One way to investigate these relationships is by considering the space of protein folds and how one might move from fold to fold through similarity, or potential evolutionary relationships. The many individual characterisations of fold space presented in the literature can tell us a lot about how well the current Protein Data Bank represents protein fold space, how convergence and divergence may affect protein evolution, how proteins affect the whole of which they are part, and how proteins themselves function. A synthesis of these different approaches and viewpoints seems the most likely way to further our knowledge of protein structure evolution and thus, facilitate improved protein structure design and prediction.

Systematic conformation-to-phenotype mapping via limited deep-sequencing of proteins

Non-native conformations drive protein misfolding diseases, complicate bioengineering efforts, and fuel molecular evolution. No current experimental technique is well-suited for elucidating them and their phenotypic effects. Especially intractable are the transient conformations populated by intrinsically disordered proteins. We describe an approach to systematically discover, stabilize, and purify native and non-native conformations, generated in vitro or in vivo, and directly link conformations to molecular, organismal, or evolutionary phenotypes. This approach involves high-throughput disulfide scanning (HTDS) of the entire protein. To reveal which disulfides trap which chromatographically resolvable conformers, we devised a deep-sequencing method for double-Cys variant libraries of proteins that precisely and simultaneously locates both Cys residues within each polypeptide. HTDS of the abundant E. coli periplasmic chaperone HdeA revealed distinct classes of disordered hydrophobic conformers with variable cytotoxicity depending on where the backbone was cross-linked. HTDS can bridge conformational and phenotypic landscapes for many proteins that function in disulfide-permissive environments.

Design and characterization of a protein fold switching network

To better understand how amino acid sequence encodes protein structure, we engineered mutational pathways that connect three common folds (3α, β-grasp, and α/β-plait). The structures of proteins at high sequence-identity intersections in the pathways (nodes) were determined using NMR spectroscopy and analyzed for stability and function. To generate nodes, the amino acid sequence encoding a smaller fold is embedded in the structure of an ~50% larger fold and a new sequence compatible with two sets of native interactions is designed. This generates protein pairs with a 3α or β-grasp fold in the smaller form but an α/β-plait fold in the larger form. Further, embedding smaller antagonistic folds creates critical states in the larger folds such that single amino acid substitutions can switch both their fold and function. The results help explain the underlying ambiguity in the protein folding code and show that new protein structures can evolve via abrupt fold switching.

Reversible switching between two common protein folds in a designed system using only temperature

Naturally occurring metamorphic proteins have the ability to interconvert from one folded state to another through either a limited set of mutations or by way of a change in the local environment. Here, we show in a designed system that it is possible to switch reversibly between two of the most common monomeric folds employing only temperature changes. We demonstrate that a latent 3α state can be unmasked from an α/β-plait topology with a single V90T amino acid substitution, populating both forms simultaneously. The equilibrium between these two states exhibits temperature dependence, such that the 3α state is predominant (>90%) at 5 °C, while the α/β-plait fold is the major species (>90%) at 30 °C. We describe the structure and dynamics of these topologies, how mutational changes affect the temperature dependence, and the energetics and kinetics of interconversion. Additionally, we demonstrate how ligand-binding function can be tightly regulated by large amplitude changes in protein structure over a relatively narrow temperature range that is relevant to biology. The 3α/αβ switch thus represents a potentially useful approach for designing proteins that alter their fold topologies in response to environmental triggers. It may also serve as a model for computational studies of temperature-dependent protein stability and fold switching.

Evolutionary selection of proteins with two folds

Although most globular proteins fold into a single stable structure 1 , an increasing number have been shown to remodel their secondary and tertiary structures in response to cellular stimuli 2 . State-of-the-art algorithms 3-5 predict that these fold-switching proteins assume only one stable structure 6,7 , missing their functionally critical alternative folds. Why these algorithms predict a single fold is unclear, but all of them infer protein structure from coevolved amino acid pairs. Here, we hypothesize that coevolutionary signatures are being missed. Suspecting that over-represented single-fold sequences may be masking these signatures, we developed an approach to search both highly diverse protein superfamilies-composed of single-fold and fold-switching variants-and protein subfamilies with more fold-switching variants. This approach successfully revealed coevolution of amino acid pairs uniquely corresponding to both conformations of 56/58 fold-switching proteins from distinct families. Then, using a set of coevolved amino acid pairs predicted by our approach, we successfully biased AlphaFold2 5 to predict two experimentally consistent conformations of a candidate protein with unsolved structure. The discovery of widespread dual-fold coevolution indicates that fold-switching sequences have been preserved by natural selection, implying that their functionalities provide evolutionary advantage and paving the way for predictions of diverse protein structures from single sequences.

Structural characterization of SARS-CoV-2 dimeric ORF9b reveals potential fold-switching trigger mechanism

The constant emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants indicates the evolution and adaptation of the virus. Enhanced innate immune evasion through increased expression of viral antagonist proteins, including ORF9b, contributes to the improved transmission of the Alpha variant; hence, more attention should be paid to these viral proteins. ORF9b is an accessory protein that suppresses innate immunity via a monomer conformation by binding to Tom70. Here, we solved the dimeric structure of SARS-CoV-2 ORF9b with a long hydrophobic tunnel containing a lipid molecule that is crucial for the dimeric conformation and determined the specific lipid ligands as monoglycerides by conducting a liquid chromatography with tandem mass spectrometry analysis, suggesting an important role in the viral life cycle. Notably, a long intertwined loop accessible for host factor binding was observed in the structure. Eight phosphorylated residues in ORF9b were identified, and residues S50 and S53 were found to contribute to the stabilization of dimeric ORF9b. Additionally, we proposed a model of multifunctional ORF9b with a distinct conformation, suggesting that ORF9b is a fold-switching protein, while both lipids and phosphorylation contribute to the switching. Specifically, the ORF9b monomer interacts with Tom70 to suppress the innate immune response, whereas the ORF9b dimer binds to the membrane involving mature virion assembly. Our results provide a better understanding of the multiple functions of ORF9b.

Creative destruction: New protein folds from old

Mechanisms of emergence and divergence of protein folds pose central questions in biological sciences. Incremental mutation and stepwise adaptation explain relationships between topologically similar protein folds. However, the universe of folds is diverse and riotous, suggesting more potent and creative forces are at play. Sequence and structure similarity are observed between distinct folds, indicating that proteins with distinct folds may share common ancestry. We found evidence of common ancestry between three distinct β-barrel folds: Scr kinase family homology (SH3), oligonucleotide/oligosaccharide-binding (OB), and cradle loop barrel (CLB). The data suggest a mechanism of fold evolution that interconverts SH3, OB, and CLB. This mechanism, which we call creative destruction, can be generalized to explain many examples of fold evolution including circular permutation. In creative destruction, an open reading frame duplicates or otherwise merges with another to produce a fused polypeptide. A merger forces two ancestral domains into a new sequence and spatial context. The fused polypeptide can explore folding landscapes that are inaccessible to either of the independent ancestral domains. However, the folding landscapes of the fused polypeptide are not fully independent of those of the ancestral domains. Creative destruction is thus partially conservative; a daughter fold inherits some motifs from ancestral folds. After merger and refolding, adaptive processes such as mutation and loss of extraneous segments optimize the new daughter fold. This model has application in disease states characterized by genetic instability. Fused proteins observed in cancer cells are likely to experience remodeled folding landscapes and realize altered folds, conferring new or altered functions.

Understanding the co-evolutionary molecular mechanisms of resistance in the HIV-1 Gag and protease

Due to high human immunodeficiency virus type 1 (HIV-1) subtype C infections coupled with increasing antiretroviral treatment failure, the elucidation of complex drug resistance mutational patterns arising through protein co-evolution is required. Despite the inclusion of potent protease inhibitors Lopinavir (LPV) and Darunavir (DRV) in second- and third-line therapies, many patients still fail treatment due to the accumulation of mutations in protease (PR) and recently, Gag. To understand the co-evolutionary molecular mechanisms of resistance in the HIV-1 PR and Gag, we performed 100 ns molecular dynamic simulations on multidrug resistant PR's when bound to LPV, DRV or a mutated A431V NC|p1 Gag cleavage site (CS). Here we showed that distinct changes in PR's active site, flap and elbow regions due to several PR resistance mutations (L10F, M46I, I54V, L76V, V82A) were found to alter LPV and DRV drug binding. However, binding was significantly exacerbated when the mutant PRs were bound to the NC|p1 Gag CS. Although A431V was shown to coordinate several residues in PR, the L76V PR mutation was found to have a significant role in substrate recognition. Consequently, a greater binding affinity was observed when the mutated substrate was bound to an L76V-inclusive PR mutant (G_bind: -62.46 ± 5.75 kcal/mol) than without (Gbind: -50.34 ± 6.28 kcal/mol). These data showed that the co-selection of resistance mutations in the enzyme and substrate can simultaneously constrict regions in PR's active site whilst flexing the flaps to allow flexible movement of the substrate and multiple, complex mechanisms of resistance to occur. Communicated by Ramaswamy H. Sarma.

Simple mechanisms for the evolution of protein complexity

Proteins are tiny models of biological complexity: specific interactions among their many amino acids cause proteins to fold into elaborate structures, assemble with other proteins into higher-order complexes, and change their functions and structures upon binding other molecules. These complex features are classically thought to evolve via long and gradual trajectories driven by persistent natural selection. But a growing body of evidence from biochemistry, protein engineering, and molecular evolution shows that naturally occurring proteins often exist at or near the genetic edge of multimerization, allostery, and even new folds, so just one or a few mutations can trigger acquisition of these properties. These sudden transitions can occur because many of the physical properties that underlie these features are present in simpler proteins as fortuitous by-products of their architecture. Moreover, complex features of proteins can be encoded by huge arrays of sequences, so they are accessible from many different starting points via many possible paths. Because the bridges to these features are both short and numerous, random chance can join selection as a key factor in explaining the evolution of molecular complexity.

Metamorphic proteins under a computational microscope: Lessons from a fold-switching RfaH protein

Metamorphic proteins constitute unexpected paradigms of the protein folding problem, as their sequences encode two alternative folds, which reversibly interconvert within biologically relevant timescales to trigger different cellular responses. Once considered a rare aberration, metamorphism may be common among proteins that must respond to rapidly changing environments, exemplified by NusG-like proteins, the only transcription factors present in every domain of life. RfaH, a specialized paralog of bacterial NusG, undergoes an all-α to all-β domain switch to activate expression of virulence and conjugation genes in many animal and plant pathogens and is the quintessential example of a metamorphic protein. The dramatic nature of RfaH structural transformation and the richness of its evolutionary history makes for an excellent model for studying how metamorphic proteins switch folds. Here, we summarize the structural and functional evidence that sparked the discovery of RfaH as a metamorphic protein, the experimental and computational approaches that enabled the description of the molecular mechanism and refolding pathways of its structural interconversion, and the ongoing efforts to find signatures and general properties to ultimately describe the protein metamorphome.

Many dissimilar NusG protein domains switch between α-helix and β-sheet folds

Folded proteins are assumed to be built upon fixed scaffolds of secondary structure, α-helices and β-sheets. Experimentally determined structures of >58,000 non-redundant proteins support this assumption, though it has recently been challenged by ~100 fold-switching proteins. Though ostensibly rare, these proteins raise the question of how many uncharacterized proteins have shapeshifting–rather than fixed–secondary structures. Here, we use a comparative sequence-based approach to predict fold switching in the universally conserved NusG transcription factor family, one member of which has a 50-residue regulatory subunit experimentally shown to switch between α-helical and β-sheet folds. Our approach predicts that 24% of sequences in this family undergo similar α-helix ⇌ β-sheet transitions. While these predictions cannot be reproduced by other state-of-the-art computational methods, they are confirmed by circular dichroism and nuclear magnetic resonance spectroscopy for 10 out of 10 sequence-diverse variants. This work suggests that fold switching may be a pervasive mechanism of transcriptional regulation in all kingdoms of life.

Folded proteins are composed of secondary structures, α-helices and β-sheets, that are generally assumed to be stable. Here, the authors combine computational prediction with experimental validation to show that many sequence-diverse NusG protein domains switch completely from α-helix to β-sheet folds.

Coevolution-derived native and non-native contacts determine the emergence of a novel fold in a universally conserved family of transcription factors

The NusG protein family is structurally and functionally conserved in all domains of life. Its members directly bind RNA polymerases and regulate transcription processivity and termination. RfaH, a divergent sub-family in its evolutionary history, is known for displaying distinct features than those in NusG proteins, which allows them to regulate the expression of virulence factors in enterobacteria in a DNA sequence-dependent manner. A striking feature is its structural interconversion between an active fold, which is the canonical NusG three-dimensional structure, and an autoinhibited fold, which is distinctively novel. How this novel fold is encoded within RfaH sequence to encode a metamorphic protein remains elusive. In this work, we used publicly available genomic RfaH protein sequences to construct a complete multiple sequence alignment, which was further augmented with metagenomic sequences and curated by predicting their secondary structure propensities using JPred. Coevolving pairs of residues were calculated from these sequences using plmDCA and GREMLIN, which allowed us to detect the enrichment of key metamorphic contacts after sequence filtering. Finally, we combined our coevolutionary predictions with molecular dynamics to demonstrate that these interactions are sufficient to predict the structures of both native folds, where coevolutionary-derived non-native contacts may play a key role in achieving the compact RfaH novel fold. All in all, emergent coevolutionary signals found within RfaH sequences encode the autoinhibited and active folds of this protein, shedding light on the key interactions responsible for the action of this metamorphic protein.

Design and discovery of metamorphic proteins

Metamorphic proteins are single amino acid sequences that reversibly interconvert between multiple, dramatically different native structures, often with distinct functions. Since the discovery of the first metamorphic proteins in the early 2000s, several additional metamorphic proteins have been identified, and it was suggested that up to 4% of proteins in the PDB may switch folds. Metamorphic proteins have been found to share common features such as marginal thermostability and inconsistencies in predicted secondary structures. Outstanding challenges in the field include the search for more metamorphic proteins and the design of new proteins that switch folds. Identification of novel metamorphic proteins in nature will improve therapeutic targeting of fold-switching proteins involved in human pathology and will enhance the design of protein-based therapies. Designed fold switching proteins have applications as biosensors, molecular switches, molecular machines, and self-assembling systems.

A reappraisal of the form: function problem-theory and phenomenology

This paper is aimed at demonstrating that some geometrical and topological transformations and operations serve not only as promoters of many specific genetic and cellular events in multicellular living organisms, but also as initiators of the organization and regulation of their functions. Thus, changes in the form and structure of macromolecular and cellular systems must be directly associated to their functions. There are specific classes of enzymes that manipulate the geometry and topology of complex DNA-protein structures, and thereby they perform many important cellular processes, including segregation of daughter chromosomes, gene regulation, and DNA repair. We argue that form has an organizing power, hence a causal action, in the sense that it enables to induce functional events during different biological processes, at the supramolecular, cellular, and organismal levels of organization. Clearly, topological forms must be matched with specific kinetic and dynamical parameters to have a functional effectiveness in living systems. This effectiveness is remarkably apparent, to give an example, in the regulation of the genome functions and in cell activity. In more general terms, we try to show that the conformational plasticity of biological systems depends on different kinds of topological manipulations performed by specific families of enzymes. In doing so, they catalyze all those spatial and dynamical changes of biological structures that are suitable for the functions to be acted by the organism.

Biomolecular and Biological Applications of Solid-State NMR with Dynamic Nuclear Polarization Enhancement

Solid-state NMR spectroscopy (ssNMR) with magic-angle spinning (MAS) enables the investigation of biological systems within their native context, such as lipid membranes, viral capsid assemblies, and cells. However, such ambitious investigations often suffer from low sensitivity due to the presence of significant amounts of other molecular species, which reduces the effective concentration of the biomolecule or interaction of interest. Certain investigations requiring the detection of very low concentration species remain unfeasible even with increasing experimental time for signal averaging. By applying dynamic nuclear polarization (DNP) to overcome the sensitivity challenge, the experimental time required can be reduced by orders of magnitude, broadening the feasible scope of applications for biological solid-state NMR. In this review, we outline strategies commonly adopted for biological applications of DNP, indicate ongoing challenges, and present a comprehensive overview of biological investigations where MAS-DNP has led to unique insights.