Fold Switching Research Articles

Impact of N-Terminal Domain Conformation and Domain Interactions on RfaH Fold Switching.

RfaH is a two-domain metamorphic protein involved in transcription regulation and translation initiation. To carry out its dual functions, RfaH relies on two coupled structural changes: Domain dissociation and fold switching. In the free state, the C-terminal domain (CTD) of RfaH adopts an all-α fold and is tightly associated with the N-terminal domain (NTD). Upon binding to RNA polymerase (RNAP), the domains dissociate and the CTD transforms into an all-β fold while the NTD remains largely, but not entirely, unchanged. We test the idea that a change in the conformation of an extended β-hairpin (β3-β4) located on the NTD, helps trigger domain dissociation. To this end, we use homology modeling to construct a structure, H1, which is similar to free RfaH but with a remodeled β3-β4 hairpin. We then use an all-atom physics-based model enhanced with a dual basin structure-based potential to simulate domain separation driven by the thermal unfolding of the CTD with NTD in a fixed, folded conformation. We apply our model to both free RfaH and H1. For H1 we find, in line with our hypothesis, that the CTD exhibits lower stability and the domains dissociate at a lower temperature T, as compared to free RfaH. We do not, however, observe complete refolding to the all-β state in these simulations, suggesting that a change in β3-β4 orientation aids in, but is not sufficient for, domain dissociation. In addition, we study the reverse fold switch in which RfaH returns from a domain-open all-β state to its domain-closed all-α state. We observe a T-dependent transition rate; fold switching is slow at low T, where the CTD tends to be kinetically trapped in its all-β state, and at high-T, where the all-α state becomes unstable. Consequently, our simulations suggest an optimal T at which fold switching is most rapid. At this T, the stabilities of both folds are reduced. Overall, our study suggests that both inter-domain interactions and conformational changes within NTD may be important for the proper functioning of RfaH.

AlphaFold predictions of fold-switched conformations are driven by structure memorization

Recent work suggests that AlphaFold (AF)–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. We find that (1) AF is a weak predictor of fold switching and (2) some of its successes result from memorization of training-set structures rather than learned protein energetics. Combining >280,000 models from several implementations of AF2 and AF3, a 35% success rate was achieved for fold switchers likely in AF’s training sets. AF2’s confidence metrics selected against models consistent with experimentally determined fold-switching structures and failed to discriminate between low and high energy conformations. Further, AF captured only one out of seven experimentally confirmed fold switchers outside of its training sets despite extensive sampling of an additional ~280,000 models. Several observations indicate that AF2 has memorized structural information during training, and AF3 misassigns coevolutionary restraints. These limitations constrain the scope of successful predictions, highlighting the need for physically based methods that readily predict multiple protein conformations.

A Transient, Excited Species of the Autoinhibited α-State of the Bacterial Transcription Factor RfaH May Represent an Early Intermediate on the Fold-Switching Pathway.

RfaH is a two-domain transcription factor in which the C-terminal domain switches fold from an α-helical hairpin to a β-roll upon binding the ops-paused RNA polymerase. To ascertain the presence of a sparsely populated excited state that may prime the autoinhibited resting state of RfaH for binding ops-paused RNA polymerase, we carried out a series of NMR-based exchange experiments to probe for conformational exchange on the millisecond time scale. Quantitative analysis of these data reveals exchange between major ground (∼95%) and sparsely populated excited (∼5%) states with an exchange lifetime of ∼3 ms involving residues at the interface between the N-terminal and C-terminal domains formed by the β3/β4 hairpin and helix α3 of the N-terminal domain and helices α4 and α5 of the C-terminal domain. The largest 15N backbone chemical shift differences are associated with the β3/β4 hairpin, leading us to suggest that the excited state may involve a rigid body lateral displacement/rotation away from the C-terminal domain to adopt a position similar to that seen in the active RNA polymerase-bound state. Such a rigid body reorientation would result in a reduction in the interface between the N- and C-terminal domains with the possible introduction of a cavity or cavities. This hypothesis is supported by the observation that the population of the excited species and the exchange rate of interconversion between ground and excited states are reduced at a high (2.5 kbar) pressure. Mechanistic implications for fold switching of the C-terminal domain in the context of RNA polymerase binding are discussed.

A Conserved Ribosomal Protein Has Entirely Dissimilar Structures in Different Organisms.

Ribosomes from different species can markedly differ in their composition by including dozens of ribosomal proteins that are unique to specific lineages but absent in others. However, it remains unknown how ribosomes acquire new proteins throughout evolution. Here, to help answer this question, we describe the evolution of the ribosomal protein msL1/msL2 that was recently found in ribosomes from the parasitic microorganism clade, microsporidia. We show that this protein has a conserved location in the ribosome but entirely dissimilar structures in different organisms: in each of the analyzed species, msL1/msL2 exhibits an altered secondary structure, an inverted orientation of the N-termini and C-termini on the ribosomal binding surface, and a completely transformed 3D fold. We then show that this fold switching is likely caused by changes in the ribosomal msL1/msL2-binding site, specifically, by variations in rRNA. These observations allow us to infer an evolutionary scenario in which a small, positively charged, de novo-born unfolded protein was first captured by rRNA to become part of the ribosome and subsequently underwent complete fold switching to optimize its binding to its evolving ribosomal binding site. Overall, our work provides a striking example of how a protein can switch its fold in the context of a complex biological assembly, while retaining its specificity for its molecular partner. This finding will help us better understand the origin and evolution of new protein components of complex molecular assemblies-thereby enhancing our ability to engineer biological molecules, identify protein homologs, and peer into the history of life on Earth.

In silico protein dynamics in the human cytoplasm: Partial folding, misfolding, fold switching, and non-native interactions.

We examine the influence of cellular interactions in all-atom models of a section of the Homo sapiens cytoplasm on the early folding events of the three-helix bundle protein B (PB). While genetically engineered PB is known to fold in dilute water box simulations in three microseconds, the three initially unfolded PB copies in our two cytoplasm models using a similar force field did not reach the native state during 30-microsecond simulations. We did however capture the formation of all three helices in a compact native-like topology. Folding in vivo is delayed because intramolecular contact formation within PB is in direct competition with intermolecular contacts between PB and surrounding macromolecules. In extreme cases, intermolecular beta-sheets are formed. Interactions with other macromolecules are also observed to promote structure formation, for example when a PB helix in our simulations is shielded from solvent by macromolecular crowding. Sticking and crowding in our models initiate sampling of helix/sheet structural plasticity of PB. Relatedly, in past in vitro experiments, similar GA domains were shown to switch between two different folds. Finally, we also observed that stickiness between PB and the cellular environment can be modulated in our simulations through the reduction in protein hydrophobicity when we reversed PB back to the wild-type sequence. This study demonstrates that even fast-folding proteins can get stuck in non-native states in the cell, making them useful models for protein-chaperone interactions and early stages of aggregate formation relevant to cellular disease.

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.

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.

Reframing prosegment-dependent folding and limits on natural protein folding landscapes from an evolutionary perspective.

In this perspective, we propose that the folding energy landscapes of model proteases including pepsin and alpha-lytic protease (αLP), which lack thermodynamic stability and fold on the order of months to millennia, respectively, should be viewed as not evolved and fundamentally distinct from their extended zymogen forms. These proteases have evolved to fold with prosegment domains and robustly self-assemble as expected. In this manner, general protein folding principles are strengthened. In support of our view, αLP and pepsin exhibit hallmarks of frustration associated with unevolved folding landscapes, such as non-cooperativity, memory effects, and substantial kinetic trapping. The evolutionary implications of this folding strategy are considered in detail. Direct applications of this folding strategy on enzyme design, finding new drug targets, and constructing tunable folding landscapes are also discussed. Together with certain proteases, growing examples of other folding "exceptions"-including protein fold switching, functional misfolding, and prevalent inability to refold-suggests a paradigm shift in which proteins may evolve to exist in a wide range of energy landscapes and structures traditionally thought to be avoided in nature.

Empirical Study of Protein Feature Representation on Deep Belief Networks Trained With Small Data for Secondary Structure Prediction.

Protein secondary structure (SS) prediction is a classic problem of computational biology and is widely used in structural characterization and to infer homology. While most SS predictors have been trained on thousands of sequences, a previous approach had developed a compact model of training proteins that used a C-Alpha, C-Beta Side Chain (CABS)-algorithm derived energy based feature representation. Here, the previous approach is extended to Deep Belief Networks (DBN). Deep learning methods are notorious for requiring large datasets and there is a wide consensus that training deep models from scratch on small datasets, works poorly. By contrast, we demonstrate a simple DBN architecture containing a single hidden layer, trained only on the CB513 dataset. Testing on an independent set of G Switch proteins improved the Q 3 score of the previous compact model by almost 3%. The findings are further confirmed by comparison to several deep learning models which are trained on thousands of proteins. Finally, the DBN performance is also compared with Position Specific Scoring Matrix (PSSM)-profile based feature representation. The importance of (i) structural information in protein feature representation and (ii) complementary small dataset learning approaches for detection of structural fold switching are demonstrated.

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.

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.

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.

AlphaFold2 fails to predict protein fold switching.

AlphaFold2 has revolutionized protein structure prediction by leveraging sequence information to rapidly model protein folds with atomic-level accuracy. Nevertheless, previous work has shown that these predictions tend to be inaccurate for structurally heterogeneous proteins. To systematically assess factors that contribute to this inaccuracy, we tested AlphaFold2's performance on 98-fold-switching proteins, which assume at least two distinct-yet-stable secondary and tertiary structures. Topological similarities were quantified between five predicted and two experimentally determined structures of each fold-switching protein. Overall, 94% of AlphaFold2 predictions captured one experimentally determined conformation but not the other. Despite these biased results, AlphaFold2's estimated confidences were moderate-to-high for 74% of fold-switching residues, a result that contrasts with overall low confidences for intrinsically disordered proteins, which are also structurally heterogeneous. To investigate factors contributing to this disparity, we quantified sequence variation within the multiple sequence alignments used to generate AlphaFold2's predictions of fold-switching and intrinsically disordered proteins. Unlike intrinsically disordered regions, whose sequence alignments show low conservation, fold-switching regions had conservation rates statistically similar to canonical single-fold proteins. Furthermore, intrinsically disordered regions had systematically lower prediction confidences than either fold-switching or single-fold proteins, regardless of sequence conservation. AlphaFold2's high prediction confidences for fold switchers indicate that it uses sophisticated pattern recognition to search for one most probable conformer rather than protein biophysics to model a protein's structural ensemble. Thus, it is not surprising that its predictions often fail for proteins whose properties are not fully apparent from solved protein structures. Our results emphasize the need to look at protein structure as an ensemble and suggest that systematic examination of fold-switching sequences may reveal propensities for multiple stable secondary and tertiary structures.

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.

Cover Image

BiopolymersVolume 112, Issue 10 e23378 COVER IMAGEFree Access Cover Image First published: 25 October 2021 https://doi.org/10.1002/bip.23378AboutPDF ToolsRequest permissionExport citationAdd to favoritesTrack citation ShareShare Give accessShare full text accessShare full-text accessPlease review our Terms and Conditions of Use and check box below to share full-text version of article.I have read and accept the Wiley Online Library Terms and Conditions of UseShareable LinkUse the link below to share a full-text version of this article with your friends and colleagues. Learn more.Copy URL Graphical Abstract This special issue of Biopolymers highlights fold-switching proteins, which remodel their secondary and tertiary structures, and change their functions in response to cellular stimuli. The cover highlights six distinct fold switchers involved in mitosis (Mad2), proteolysis (Selecase), transcription (RfaH), oxidative stress response (CLIC1), circadian rhythm (KaiB) and the human immune system (β2-microglobulin). Topics addressed in the issue include experimental characterization of triggers and conditions that foster fold switching, as well as computational approaches that predict and characterize these shapeshifting proteins. Volume112, Issue10Special Issue: Special Issue: Fold‐Switching ProteinsOctober 2021e23378 RelatedInformation

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.

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.

Evolutionary insights into shape-shifting proteins

Over millions of years a protein that now folds into two stable structures likely favored first one configuration, then the other, before settling on both.

Folding, Fold Switching and Energy Landscape of Regulatory Protein RfaH

Folding, Fold Switching and Energy Landscape of Regulatory Protein RfaH