Evolution of fold switching in a metamorphic protein

Protein Fold Usages in Ribosomes: Another Glance to the Past

The analysis of protein fold usage, similar to codon usage, offers profound insights into the evolution of biological systems and the origins of modern proteomes. While previous studies have examined fold distribution in modern genomes, our study focuses on the comparative distribution and usage of protein folds in ribosomes across bacteria, archaea, and eukaryotes. We identify the prevalence of certain 'super-ribosome folds,' such as the OB fold in bacteria and the SH3 domain in archaea and eukaryotes. The observed protein fold distribution in the ribosomes announces the future power-law distribution where only a few folds are highly prevalent, and most are rare. Additionally, we highlight the presence of three copies of proto-Rossmann folds in ribosomes across all kingdoms, showing its ancient and fundamental role in ribosomal structure and function. Our study also explores early mechanisms of molecular convergence, where different protein folds bind equivalent ribosomal RNA structures in ribosomes across different kingdoms. This comparative analysis enhances our understanding of ribosomal evolution, particularly the distinct evolutionary paths of the large and small subunits, and underscores the complex interplay between RNA and protein components in the transition from the RNA world to modern cellular life. Transcending the concept of folds also makes it possible to group a large number of ribosomal proteins into five categories of urfolds or metafolds, which could attest to their ancestral character and common origins. This work also demonstrates that the gradual acquisition of extensions by simple but ordered folds constitutes an inexorable evolutionary mechanism. This observation supports the idea that simple but structured ribosomal proteins preceded the development of their disordered extensions.

Accurate prediction of protein function using statistics-informed graph networks

Understanding protein function is pivotal in comprehending the intricate mechanisms that underlie many crucial biological activities, with far-reaching implications in the fields of medicine, biotechnology, and drug development. However, more than 200 million proteins remain uncharacterized, and computational efforts heavily rely on protein structural information to predict annotations of varying quality. Here, we present a method that utilizes statistics-informed graph networks to predict protein functions solely from its sequence. Our method inherently characterizes evolutionary signatures, allowing for a quantitative assessment of the significance of residues that carry out specific functions. PhiGnet not only demonstrates superior performance compared to alternative approaches but also narrows the sequence-function gap, even in the absence of structural information. Our findings indicate that applying deep learning to evolutionary data can highlight functional sites at the residue level, providing valuable support for interpreting both existing properties and new functionalities of proteins in research and biomedicine.

Ancestral Reconstruction and the Evolution of Protein Energy Landscapes

A protein's sequence determines its conformational energy landscape. This, in turn, determines the protein's function. Understanding the evolution of new protein functions therefore requires understanding how mutations alter the protein energy landscape. Ancestral sequence reconstruction (ASR) has proven a valuable tool for tackling this problem. In ASR, one phylogenetically infers the sequences of ancient proteins, allowing characterization of their properties. When coupled to biophysical, biochemical, and functional characterization, ASR can reveal how historical mutations altered the energy landscape of ancient proteins, allowing the evolution of enzyme activity, altered conformations, binding specificity, oligomerization, and many other protein features. In this article, we review how ASR studies have been used to dissect the evolution of energy landscapes. We also discuss ASR studies that reveal how energy landscapes have shaped protein evolution. Finally, we propose that thinking about evolution from the perspective of an energy landscape can improve how we approach and interpret ASR studies.

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.

Back in time to the Gly-rich prototype of the phosphate binding elementary function

Binding of nucleotides and their derivatives is one of the most ancient elementary functions dating back to the Origin of Life. We review here the works considering one of the key elements in binding of (di)nucleotide-containing ligands - phosphate binding. We start from a brief discussion of major participants, conditions, and events in prebiotic evolution that resulted in the Origin of Life. Tracing back to the basic functions, including metal and phosphate binding, and, potentially, formation of primitive protein-protein interactions, we focus here on the phosphate binding. Critically assessing works on the structural, functional, and evolutionary aspects of phosphate binding, we perform a simple computational experiment reconstructing its most ancient and generic sequence prototype. The profiles of the phosphate binding signatures have been derived in form of position-specific scoring matrices (PSSMs), their peculiarities depending on the type of the ligands have been analyzed, and evolutionary connections between them have been delineated. Then, the apparent prototype that gave rise to all relevant phosphate-binding signatures had also been reconstructed. We show that two major signatures of the phosphate binding that discriminate between the binding of dinucleotide- and nucleotide-containing ligands are GxGxxG and GxxGxG, respectively. It appears that the signature archetypal for dinucleotide-containing ligands is more generic, and it can frequently bind phosphate groups in nucleotide-containing ligands as well. The reconstructed prototype's key signature GxGGxG underlies the role of glycine residues in providing flexibility and interactions necessary for binding the phosphate groups. The prototype also contains other ancient amino acids, valine, and alanine, showing versatility towards evolutionary design and functional diversification.

Evolution of Phosphorylase Activity in an Ancestral Glycosyltransferase

The reconstruction of ancestral sequences can offer a glimpse into the fascinating process of molecular evolution by exposing the adaptive pathways that shape the proteins found in nature today. Here, we track the evolution of the carbohydrate-active enzymes responsible for the synthesis and turnover of mannogen, a critical carbohydrate reserve in Leishmania parasites. Biochemical characterization of resurrected enzymes demonstrated that mannoside phosphorylase activity emerged in an ancestral bacterial mannosyltransferase, and later disappeared in the process of horizontal gene transfer and gene duplication in Leishmania. By shuffling through plausible historical sequence space in an ancestral mannosyltransferase, we found that mannoside phosphorylase activity could be toggled on through various combinations of mutations at positions outside of the active site. Molecular dynamics simulations showed that such mutations can affect loop rigidity and shield the active site from water molecules that disrupt key interactions, allowing α-mannose 1-phosphate to adopt a catalytically productive conformation. These findings highlight the importance of subtle distal mutations in protein evolution and suggest that the vast collection of natural glycosyltransferases may be a promising source of engineering templates for the design of tailored phosphorylases.

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.

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.

The First Nucleic Acid Strands May Have Grown on Peptides via Primeval Reverse Translation

The central dogma of molecular biology dictates that, with only a few exceptions, information proceeds from DNA to protein through an RNA intermediate. Examining the enigmatic steps from prebiotic to biological chemistry, we take another road suggesting that primordial peptides acted as template for the self-assembly of the first nucleic acids polymers. Arguing in favour of a sort of archaic "reverse translation" from proteins to RNA, our basic premise is a Hadean Earth where key biomolecules such as amino acids, polypeptides, purines, pyrimidines, nucleosides and nucleotides were available under different prebiotically plausible conditions, including meteorites delivery, shallow ponds and hydrothermal vents scenarios. Supporting a protein-first scenario alternative to the RNA world hypothesis, we propose the primeval occurrence of short two-dimensional peptides termed "selective amino acid- and nucleotide-matching oligopeptides" (henceforward SANMAOs) that noncovalently bind at the same time the polymerized amino acids and the single nucleotides dispersed in the prebiotic milieu. In this theoretical paper, we describe the chemical features of this hypothetical oligopeptide, its biological plausibility and its virtues from an evolutionary perspective. We provide a theoretical example of SANMAO's selective pairing between amino acids and nucleosides, simulating a poly-Glycine peptide that acts as a template to build a purinic chain corresponding to the glycine's extant triplet codon GGG. Further, we discuss how SANMAO might have endorsed the formation of low-fidelity RNA's polymerized strains, well before the appearance of the accurate genetic material's transmission ensured by the current translation apparatus.

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.

Impact of ancestral sequence reconstruction on mechanistic and structural enzymology

Ancestral sequence reconstruction (ASR) provides insight into the changes within a protein sequence across evolution. More specifically, it can illustrate how specific amino acid changes give rise to different phenotypes within a protein family. Over the last few decades it has established itself as a powerful technique for revealing molecular common denominators that govern enzyme function. Here, we describe the strength of ASR in unveiling catalytic mechanisms and emerging phenotypes for a range of different proteins, also highlighting biotechnological applications the methodology can provide.

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.

Molecular determinants of protein evolvability

The plethora of biological functions that sustain life is rooted in the remarkable evolvability of proteins. An emerging view highlights the importance of a protein's initial state in dictating evolutionary success. A deeper comprehension of the mechanisms that govern the evolvability of these initial states can provide invaluable insights into protein evolution. In this review, we describe several molecular determinants of protein evolvability, unveiled by experimental evolution and ancestral sequence reconstruction studies. We further discuss how genetic variation and epistasis can promote or constrain functional innovation and suggest putative underlying mechanisms. By establishing a clear framework for these determinants, we provide potential indicators enabling the forecast of suitable evolutionary starting points and delineate molecular mechanisms in need of deeper exploration.

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.

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.

Thermodynamic Evolution of a Metamorphic Protein: A Theoretical-Computational Study of Human Lymphotactin

Metamorphic, or fold-switching, proteins feature different folds that are physiologically relevant. The human chemokine XCL1 (or Lymphotactin) is a metamorphic protein that features two native states, an [Formula: see text] and an all[Formula: see text] fold, which have similar stability at physiological condition. Here, extended molecular dynamics (MD) simulations, principal component analysis of atomic fluctuations and thermodynamic modeling based on both the configurational volume and free energy landscape, are used to obtain a detailed characterization of the conformational thermodynamics of human Lymphotactin and of one of its ancestors (as was previously obtained by genetic reconstruction). Comparison of our computational results with the available experimental data show that the MD-based thermodynamics can explain the experimentally observed variation of the conformational equilibrium between the two proteins. In particular, our computational data provide an interpretation of the thermodynamic evolution in this protein, revealing the relevance of the configurational entropy and of the shape of the free energy landscape within the essential space (i.e., the space defined by the generalized internal coordinates providing the largest, typically non-Gaussian, structural fluctuations).

What defines a chemokine? - The curious case of CXCL17

Chemotactic cytokines (chemokines) are a group of around 40 small proteins which share a similar protein fold and are well known for their ability to direct the migration of leukocytes to a variety of tissue locations. CXCL17 was the last member of the chemokine family to be assigned and was admitted to the family based on theoretical modelling of the CXCL17 structure and chemotactic activity for monocytes and dendritic cells. Of Interest, CXCL17 expression appears to be restricted to mucosal tissues such as the tongue, stomach and lung, suggestive of specific roles at these locations. A putative CXCL17 receptor, GPR35 was reportedly identified and mice deficient in CXCL17 were generated and characterised. More recently, however, some apparent contradictions regarding aspects of CXCL17 biology have been raised by ourselves and others. Notably, GPR35 appears to be a receptor for the serotonin metabolite 5-hydroxyindoleacetic acid rather than for CXCL17 and modelling of CXCL17 using a variety of platforms fails to identify a chemokine-like fold. In this article, we summarize the discovery of CXCL17 and discuss key papers describing the subsequent characterisation of this protein. Ultimately, we pose the question, 'What defines a chemokine?' (185 words).

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.

Genetic Perturbation Alters Functional Substates in Alkaline Phosphatase

Enzymes inherently exhibit molecule-to-molecule heterogeneity in their conformational and functional states, which is considered to be a key to the evolution of new functions. Single-molecule enzyme assays enable us to directly observe such multiple functional states or functional substates. Here, we quantitatively analyzed functional substates in the wild-type and 69 single-point mutants of Escherichia coli alkaline phosphatase by employing a high-throughput single-molecule assay with a femtoliter reactor array device. Interestingly, many mutant enzymes exhibited significantly heterogeneous functional substates with various types, while the wild-type enzyme showed a highly homogeneous substate. We identified a correlation between the degree of functional substates and the level of improvement in promiscuous activities. Our work provides much comprehensive evidence that the functional substates can be easily altered by mutations, and the evolution toward a new catalytic activity may involve the modulation of the functional substates.

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.

Protein dynamics: The future is bright and complicated!

Biological life depends on motion, and this manifests itself in proteins that display motion over a formidable range of time scales spanning from femtoseconds vibrations of atoms at enzymatic transition states, all the way to slow domain motions occurring on micro to milliseconds. An outstanding challenge in contemporary biophysics and structural biology is a quantitative understanding of the linkages among protein structure, dynamics, and function. These linkages are becoming increasingly explorable due to conceptual and methodological advances. In this Perspective article, we will point toward future directions of the field of protein dynamics with an emphasis on enzymes. Research questions in the field are becoming increasingly complex such as the mechanistic understanding of high-order interaction networks in allosteric signal propagation through a protein matrix, or the connection between local and collective motions. In analogy to the solution to the "protein folding problem," we argue that the way forward to understanding these and other important questions lies in the successful integration of experiment and computation, while utilizing the present rapid expansion of sequence and structure space. Looking forward, the future is bright, and we are in a period where we are on the doorstep to, at least in part, comprehend the importance of dynamics for biological function.

Physiology of chemokines in the cancer microenvironment

Chemokines are chemotactic cytokines whose canonical functions govern movement of receptor-expressing cells along chemical gradients. Chemokines are a physiological system that is finely tuned by ligand and receptor expression, ligand or receptor oligomerization, redundancy, expression of atypical receptors, and non-GPCR binding partners that cumulatively influence discrete pharmacological signaling responses and cellular functions. In cancer, chemokines play paradoxical roles in both the directed emigration of metastatic, receptor-expressing cancer cells out of the tumor as well as immigration of tumor-infiltrating immune cells that culminate in a tumor-unique immune microenvironment. In the age of precision oncology, strategies to effectively harness the power of immunotherapy requires consideration of chemokine gradients within the unique spatial topography and temporal influences with heterogeneous tumors. In this article, we review current literature on the diversity of chemokine ligands and their cellular receptors that detect and process chemotactic gradients and illustrate how differences between ligand recognition and receptor activation influence the signaling machinery that drives cellular movement into and out of the tumor microenvironment. Facets of chemokine physiology across discrete cancer immune phenotypes are contrasted to existing chemokine-centered therapies in cancer.

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.

Intrinsically fluorescent polyureas toward conformation-assisted metamorphosis, discoloration and intracellular drug delivery

Peptidomimetic polymers have attracted increasing interest because of the advantages of facile synthesis, high molecular tunability, resistance to degradation, and low immunogenicity. However, the presence of non-native linkages compromises their ability to form higher ordered structures and protein-inspired functions. Here we report a class of amino acid-constructed polyureas with molecular weight- and solvent-dependent helical and sheet-like conformations as well as green fluorescent protein-mimic autofluorescence with aggregation-induced emission characteristics. The copolymers self-assemble into vesicles and nanotubes and exhibit H-bonding-mediated metamorphosis and discoloration behaviors. We show that these polymeric vehicles with ultrahigh stability, superfast responsivity and conformation-assisted cell internalization efficiency could act as an “on-off” switchable nanocarrier for specific intracellular drug delivery and effective cancer theranosis in vitro and in vivo. This work provides insights into the folding and hierarchical assembly of biomacromolecules, and a new generation of bioresponsive polymers and nonconventional luminescent aliphatic materials for diverse applications.

Biomimetic materials are of interest but can often suffer from limitations caused by the non-native linkages used. Here, the authors report on the creation of amino acid constructed polyureas which can self-assemble into vesicles and nanotubes with aggregation induced fluorescence and the potential for drug delivery applications.

Mechanisms of protein evolution

How do proteins evolve? How do changes in sequence mediate changes in protein structure, and in turn in function? This question has multiple angles, ranging from biochemistry and biophysics to evolutionary biology. This review provides a brief integrated view of some key mechanistic aspects of protein evolution. First, we explain how protein evolution is primarily driven by randomly acquired genetic mutations and selection for function, and how these mutations can even give rise to completely new folds. Then, we also comment on how phenotypic protein variability, including promiscuity, transcriptional and translational errors, may also accelerate this process, possibly via "plasticity-first" mechanisms. Finally, we highlight open questions in the field of protein evolution, with respect to the emergence of more sophisticated protein systems such as protein complexes, pathways, and the emergence of pre-LUCA enzymes.

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.

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.

Proteins' Evolution upon Point Mutations

As the reader must be already aware, state-of-the-art protein folding prediction methods have reached a smashing success in their goal of accurately determining the three-dimensional structures of proteins. Yet, a solution to simple problems such as the effects of protein point mutations on their (i) native conformation; (ii) marginal stability; (iii) ensemble of high-energy nativelike conformations; and (iv) metamorphism propensity and, hence, their evolvability, remains as an unsolved problem. As a plausible solution to the latter, some properties of the amide hydrogen-deuterium exchange, a highly sensitive probe of the structure, stability, and folding of proteins, are assessed from a new perspective. The preliminary results indicate that the protein marginal stability change upon point mutations provides the necessary and sufficient information to estimate, through a Boltzmann factor, the evolution of the amide hydrogen exchange protection factors and, consequently, that of the ensemble of folded conformations coexisting with the native state. This work contributes to our general understanding of the effects of point mutations on proteins and may spur significant progress in our efforts to develop methods to determine the appearance of new folds and functions accurately.

Global Fold Switching of the RafH Protein: Diverse Structures with a Conserved Pathway

It is generally believed that a protein's sequence uniquely determines its structure, the basis for a protein to perform biological functions. However, as a representative metamorphic protein, RfaH can be encoded by a single amino acid sequence into two distinct native state structures. Its C-terminal domain (CTD) either takes an all-α-helical configuration to pack tightly with its N-terminal domain (NTD), or the CTD disassociates from the NTD, transforms into an all-β-barrel fold, and further attaches to the ribosome, leaving the NTD exposed to bind RNA polymerases. Therefore, the RfaH protein couples transcription and translation processes. Although previous studies have provided a preliminary understanding of its function, the full course of the conformational change of RfaH-CTD at the atomic level is elusive. We used teDA2, a feature space-based enhanced sampling protocol, to explore the transformation of RfaH-CTD. We found that it undergoes a large-scale structural rearrangement, with characteristic spectra as the fingerprint, and a global unfolding transition with a tighter and energetically moderate molten globule-like nucleus formed in between. The formation of this nucleus limits the possible intermediate conformations, facilitates the formation of secondary and tertiary structures, and thus ensures the efficiency of transformation. The key features along the transition path disclosed from this work are likely associated with the evolution of RfaH, such that encoding a single sequence into multiple folds with distinct biological functions is energetically unhindered.

Curious Case of MAD2 Protein: Diverse Folding Intermediates Leading to Alternate Native States

Anfinsen's dogma postulates that for one sequence there will be only one unique structure that is necessary for the functioning of the protein. However, over the years there have been a number of departures from this postulate. As far as function is considered, there are growing examples of proteins that "moonlight", perform multiple unrelated functions. With the discovery of intrinsically disordered proteins, morpheeins, chameleonic sequences, and metamorphic proteins that can switch folds, we have acquired a more nuanced understanding of protein folding and dynamics. Appearing to apparently contradict the classical folding paradigm, metamorphic proteins are considered exotic species. In this work, we have explored the free energy landscape and folding pathways of the metamorphic protein MAD2 which is an important component of the spindle checkpoint. It coexists in two alternate states: the inactive open state and the active closed state. Using a dual-basin structure-based model approach we have shown that a variety of intermediates and multiple pathways are available to MAD2 to fold into its alternate forms. This approach involves performing molecular dynamics simulations of coarse-grained models of MAD2 where the structural information regarding both of its native conformations is explicitly included in terms of their native contacts in the force field used. Detailed analyses have indicated that some of the contacts within the protein play a key role in determining which folding pathway will be selected and point to a probable long-range communication between the N and the C termini of the protein that seems to control its folding. Finally, our work also provides a rationale for the experimentally observed preference of the ΔC₁₀ variant of MAD2 to exist in the open state.

Design principles of protein switches

Protein switches perform essential roles in many biological processes and are exciting targets for de novo protein design, which aims to produce proteins of arbitrary shape and functionality. However, the biophysical requirements for switch function - multiple conformational states, fine-tuned energetics, and stimuli-responsiveness - pose a formidable challenge for design by computation (or intuition). A variety of methods have been developed toward tackling this challenge, usually taking inspiration from the wealth of sequence and structural information available for naturally occurring protein switches. More recently, modular switches have been designed computationally, and new methods have emerged for sampling unexplored structure space, providing promising new avenues toward the generation of purpose-built switches and de novo signaling systems for cellular engineering.

A sequence-based method for predicting extant fold switchers that undergo α-helix ↔ β-strand transitions

Extant fold-switching proteins remodel their secondary structures and change their functions in response to cellular stimuli, regulating biological processes and affecting human health. Despite their biological importance, these proteins remain understudied. Predictive methods are needed to expedite the process of discovering and characterizing more of these shapeshifting proteins. Most previous approaches require a solved structure or all-atom simulations, greatly constraining their use. Here, we propose a high-throughput sequence-based method for predicting extant fold switchers that transition from α-helix in one conformation to β-strand in the other. This method leverages two previous observations: (a) α-helix ↔ β-strand prediction discrepancies from JPred4 are a robust predictor of fold switching, and (b) the fold-switching regions (FSRs) of some extant fold switchers have different secondary structure propensities when expressed by themselves (isolated FSRs) than when expressed within the context of their parent protein (contextualized FSRs). Combining these two observations, we ran JPred4 on 99-fold-switching proteins and found strong correspondence between predicted and experimentally observed α-helix ↔ β-strand discrepancies. To test the overall robustness of this finding, we randomly selected regions of proteins not expected to switch folds (single-fold proteins) and found significantly fewer predicted α-helix ↔ β-strand discrepancies. Combining these discrepancies with the overall percentage of predicted secondary structure, we developed a classifier to identify extant fold switchers (Matthews correlation coefficient of .71). Although this classifier had a high false-negative rate (7/17), its false-positive rate was very low (2/136), suggesting that it can be used to predict a subset of extant fold switchers from a multitude of available genomic sequences.

APRIL limits atherosclerosis by binding to heparan sulfate proteoglycans

Atherosclerotic cardiovascular disease causes heart attacks and strokes, which are the leading causes of mortality worldwide¹. The formation of atherosclerotic plaques is initiated when low-density lipoproteins bind to heparan-sulfate proteoglycans (HSPGs)² and become trapped in the subendothelial space of large and medium size arteries, which leads to chronic inflammation and remodelling of the artery wall². A proliferation-inducing ligand (APRIL) is a cytokine that binds to HSPGs³, but the physiology of this interaction is largely unknown. Here we show that genetic ablation or antibody-mediated depletion of APRIL aggravates atherosclerosis in mice. Mechanistically, we demonstrate that APRIL confers atheroprotection by binding to heparan sulfate chains of heparan-sulfate proteoglycan 2 (HSPG2), which limits the retention of low-density lipoproteins, accumulation of macrophages and formation of necrotic cores. Indeed, antibody-mediated depletion of APRIL in mice expressing heparan sulfate-deficient HSPG2 had no effect on the development of atherosclerosis. Treatment with a specific anti-APRIL antibody that promotes the binding of APRIL to HSPGs reduced experimental atherosclerosis. Furthermore, the serum levels of a form of human APRIL protein that binds to HSPGs, which we termed non-canonical APRIL (nc-APRIL), are associated independently of traditional risk factors with long-term cardiovascular mortality in patients with atherosclerosis. Our data reveal properties of APRIL that have broad pathophysiological implications for vascular homeostasis.

Predictable fold switching by the SARS-CoV-2 protein ORF9b

Extant fold-switching proteins remodel their secondary structures and change their functions in response to environmental stimuli. These shapeshifting proteins regulate biological processes and are associated with a number of diseases, including tuberculosis, cancer, Alzheimer's, and autoimmune disorders. Thus, predictive methods are needed to identify more fold-switching proteins, especially since all naturally occurring instances have been discovered by chance. In response to this need, two high-throughput predictive methods have recently been developed. Here we test them on ORF9b, a newly discovered fold switcher and potential therapeutic target from the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Promisingly, both methods correctly indicate that ORF9b switches folds. We then tested the same two methods on ORF9b1, the ORF9b homolog from SARS-CoV-1. Again, both methods predict that ORF9b1 switches folds, a finding consistent with experimental binding studies. Together, these results (a) demonstrate that protein fold switching can be predicted using high-throughput computational approaches and (b) suggest that fold switching might be a general characteristic of ORF9b homologs.

Thoughts on the Protein's Native State

The presence of metamorphism in the protein's native state is not yet fully understood. To shed light on this issue, we present an assessment, in terms of the amide hydrogen exchange protection factor, that aims to determine the possible existence of structural fluctuations in the native-state consistent with both the upper-bound marginal stability of proteins and the presence of metamorphism. The preliminary results enable us to conclude that the native-state metamorphism is, indeed, more probable than previously thought.

Metamorphic Protein Folding Encodes Multiple Anti-Candida Mechanisms in XCL1

Candida species cause serious infections requiring prolonged and sometimes toxic therapy. Antimicrobial proteins, such as chemokines, hold great interest as potential additions to the small number of available antifungal drugs. Metamorphic proteins reversibly switch between multiple different folded structures. XCL1 is a metamorphic, antimicrobial chemokine that interconverts between the conserved chemokine fold (an α–β monomer) and an alternate fold (an all-β dimer). Previous work has shown that human XCL1 kills C. albicans but has not assessed whether one or both XCL1 folds perform this activity. Here, we use structurally locked engineered XCL1 variants and Candida killing assays, adenylate kinase release assays, and propidium iodide uptake assays to demonstrate that both XCL1 folds kill Candida, but they do so via different mechanisms. Our results suggest that the alternate fold kills via membrane disruption, consistent with previous work, and the chemokine fold does not. XCL1 fold-switching thus provides a mechanism to regulate the XCL1 mode of antifungal killing, which could protect surrounding tissue from damage associated with fungal membrane disruption and could allow XCL1 to overcome candidal resistance by switching folds. This work provides inspiration for the future design of switchable, multifunctional antifungal therapeutics.

The IDIP framework for assessing protein function and its application to the prion protein

The quest to determine the function of a protein can represent a profound challenge. Although this task is the mandate of countless research groups, a general framework for how it can be approached is conspicuously lacking. Moreover, even expectations for when the function of a protein can be considered to be 'known' are not well defined. In this review, we begin by introducing concepts pertinent to the challenge of protein function assignments. We then propose a framework for inferring a protein's function from four data categories: 'inheritance', 'distribution', 'interactions' and 'phenotypes' (IDIP). We document that the functions of proteins emerge at the intersection of inferences drawn from these data categories and emphasise the benefit of considering them in an evolutionary context. We then apply this approach to the cellular prion protein (PrP^C ), well known for its central role in prion diseases, whose function continues to be considered elusive by many investigators. We document that available data converge on the conclusion that the function of the prion protein is to control a critical post-translational modification of the neural cell adhesion molecule in the context of epithelial-to-mesenchymal transition and related plasticity programmes. Finally, we argue that this proposed function of PrP^C has already passed the test of time and is concordant with the IDIP framework in a way that other functions considered for this protein fail to achieve. We anticipate that the IDIP framework and the concepts analysed herein will aid the investigation of other proteins whose primary functional assignments have thus far been intractable.

Metamorphic proteins: the Janus proteins of structural biology

The structural paradigm that the sequence of a protein encodes for a unique three-dimensional native fold does not acknowledge the intrinsic plasticity encapsulated in conformational free energy landscapes. Metamorphic proteins are a recently discovered class of biomolecules that illustrate this plasticity by folding into at least two distinct native state structures of comparable stability in the absence of ligands or cofactors to facilitate fold-switching. The expanding list of metamorphic proteins clearly shows that these proteins are not mere aberrations in protein evolution, but may have actually been a consequence of distinctive patterns in selection pressure such as those found in virus–host co-evolution. In this review, we describe the structure–function relationships observed in well-studied metamorphic protein systems, with specific focus on how functional residues are sequestered or exposed in the two folds of the protein. We also discuss the implications of metamorphosis for protein evolution and the efforts that are underway to predict metamorphic systems from sequence properties alone.

Resurrected Ancestors Reveal Origins of Metamorphism in XCL1

In a recent study, Dishman et al. resurrected ancestors of the metamorphic chemokine, XCL1, inferred through phylogenetics, and found that metamorphism arose in the XCL1 lineage ~150 million years ago. A zigzagging evolutionary path suggests that the metamorphic properties are adaptive and reveals three design principles that could be used for technological applications.

Fragile protein folds: Sequence and environmental factors affecting the equilibrium of two interconverting, stably folded protein conformations

Recent research on fold-switching metamorphic proteins has revealed some notable exceptions to Anfinsen's hypothesis of protein folding. We have previously described how a single point mutation can enable a well-folded protein domain, one of the two PAS (Per-ARNT-Sim) domains of the human ARNT (aryl hydrocarbon receptor nuclear translocator) protein, to interconvert between two conformers related by a slip of an internal β-strand. Using this protein as a test case, we advance the concept of a "fragile fold," a protein fold that can reversibly rearrange into another fold that differs by a substantial number of hydrogen bonds, entailing reorganization of single secondary structure elements to more drastic changes seen in metamorphic proteins. Here we use a battery of biophysical tests to examine several factors affecting the equilibrium between the two conformations of the switching ARNT PAS-B Y456T protein. Of note, we find that factors which impact the HI loop preceding the shifted Iβ-strand affect both the equilibrium levels of the two conformers and the denatured state which links them in the interconversion process. Finally, we describe small molecules that selectively bind to and stabilize the wildtype conformation of ARNT PAS-B. These studies form a toolkit for studying fragile protein folds and could enable ways to modulate the biological functions of such fragile folds, both in natural and engineered proteins.