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.

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.

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.

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.

The C-terminal domain of transcription factor RfaH: Folding, fold switching and energy landscape

We simulate the folding and fold switching of the C-terminal domain (CTD) of the transcription factor RfaH using an all-atom physics-based model augmented with a dual-basin structure-based potential energy term. We show that this hybrid model captures the essential thermodynamic behavior of this metamorphic domain, that is, a change in the global free energy minimum from an α-helical hairpin to a 5-stranded β-barrel upon the dissociation of the CTD from the rest of the protein. Using Monte Carlo sampling techniques, we then analyze the energy landscape of the CTD in terms of progress variables for folding toward the two folds. We find that, below the folding transition, the energy landscape is characterized by a single, dominant funnel to the native β-barrel structure. The absence of a deep funnel to the α-helical hairpin state reflects a negligible population of this fold for the isolated CTD. We observe, however, a higher α-helix structure content in the unfolded state compared to results from a similar but fold switch-incompetent version of our model. Moreover, in folding simulations started from an extended chain conformation we find transiently formed α-helical structure, occurring early in the process and disappearing as the chain progresses toward the thermally stable β-barrel state.

NusG, an Ancient Yet Rapidly Evolving Transcription Factor

Timely and accurate RNA synthesis depends on accessory proteins that instruct RNA polymerase (RNAP) where and when to start and stop transcription. Among thousands of transcription factors, NusG/Spt5 stand out as the only universally conserved family of regulators. These proteins interact with RNAP to promote uninterrupted RNA synthesis and with diverse cellular partners to couple transcription to RNA processing, modification or translation, or to trigger premature termination of aberrant transcription. NusG homologs are present in all cells that utilize bacterial-type RNAP, from endosymbionts to plants, underscoring their ancient and essential function. Yet, in stark contrast to other core RNAP components, NusG family is actively evolving: horizontal gene transfer and sub-functionalization drive emergence of NusG paralogs, such as bacterial LoaP, RfaH, and UpxY. These specialized regulators activate a few (or just one) operons required for expression of antibiotics, capsules, secretion systems, toxins, and other niche-specific macromolecules. Despite their common origin and binding site on the RNAP, NusG homologs differ in their target selection, interacting partners and effects on RNA synthesis. Even among housekeeping NusGs from diverse bacteria, some factors promote pause-free transcription while others slow the RNAP down. Here, we discuss structure, function, and evolution of NusG proteins, focusing on unique mechanisms that determine their effects on gene expression and enable bacterial adaptation to diverse ecological niches.

Structural fluctuations and mechanical stabilities of the metamorphic protein RfaH

RfaH is a compact two-domain bacterial transcription factor that functions both as a regulator of transcription and an enhancer of translation. Underpinning the dual functional roles of RfaH is a partial but dramatic fold switch, which completely transforms the ~50-amino acid C-terminal domain (CTD) from an all-α state to an all-β state. The fold switch of the CTD occurs when RfaH binds to RNA polymerase (RNAP), however, the details of how this structural transformation is triggered is not well understood. Here we use all-atom Monte Carlo simulations to characterize structural fluctuations and mechanical stability properties of the full-length RfaH and the CTD as an isolated fragment. In agreement with experiments, we find that interdomain contacts are crucial for maintaining a stable, all-α CTD in free RfaH. To probe mechanical properties, we use pulling simulations to measure the work required to inflict local deformations at different positions along the chain. The resulting mechanical stability profile reveals that free RfaH can be divided into a "rigid" part and a "soft" part, with a boundary that nearly coincides with the boundary between the two domains. We discuss the potential role of this feature for how fold switching may be triggered by interaction with RNAP.

Differential Local Stability Governs the Metamorphic Fold Switch of Bacterial Virulence Factor RfaH

RfaH, a two-domain protein from a universally conserved NusG/Spt5 family of regulators, is required for the transcription and translation of long virulence and conjugation operons in many Gram-negative bacterial pathogens. Escherichia coli RfaH action is controlled by a unique large-scale structural rearrangement triggered by recruitment to transcription elongation complexes through a specific DNA element. Upon recruitment, the C-terminal domain of RfaH refolds from an α-hairpin, which is bound to RNA polymerase binding site within the N-terminal domain, into an unbound β-barrel that interacts with the ribosome. Although structures of the autoinhibited (α-hairpin) and active (β-barrel) states and plausible refolding pathways have been reported, how this reversible switch is encoded within RfaH sequence and structure is poorly understood. Here, we combined hydrogen-deuterium exchange measurements by mass spectrometry and nuclear magnetic resonance with molecular dynamics to evaluate the differential local stability between both RfaH folds. Deuteron incorporation reveals that the tip of the C-terminal hairpin (residues 125-145) is stably folded in the autoinhibited state (∼20% deuteron incorporation), whereas the rest of this domain is highly flexible (>40% deuteron incorporation), and its flexibility only decreases in the β-folded state. Computationally predicted ΔG agree with these results by displaying similar anisotropic stability within the tip of the α-hairpin and on neighboring N-terminal domain residues. Remarkably, the β-folded state shows comparable structural flexibility than nonmetamorphic homologs. Our findings provide information critical for understanding the metamorphic behavior of RfaH and other chameleon proteins and for devising targeted strategies to combat bacterial infections.

Ancient Transcription Factors in the News

In every cell from bacteria to mammals, NusG-like proteins bind transcribing RNA polymerase to modulate the rate of nascent RNA synthesis and to coordinate it with numerous cotranscriptional processes that ultimately determine the transcript fate. Housekeeping NusG factors regulate expression of the bulk of the genome, whereas their highly specialized paralogs control just a few targets.

In every cell from bacteria to mammals, NusG-like proteins bind transcribing RNA polymerase to modulate the rate of nascent RNA synthesis and to coordinate it with numerous cotranscriptional processes that ultimately determine the transcript fate. Housekeeping NusG factors regulate expression of the bulk of the genome, whereas their highly specialized paralogs control just a few targets. In Escherichia coli, NusG stimulates silencing of horizontally acquired genes, while its paralog RfaH counters NusG action by activating a subset of these genes. Acting alone or as part of regulatory complexes, NusG factors can promote uninterrupted RNA synthesis, bring about transcription pausing or premature termination, modulate RNA processing, and facilitate translation. Recent structural and mechanistic studies of NusG homologs from all domains of life reveal molecular details of multifaceted interactions that underpin their unexpectedly diverse regulatory roles. NusG proteins share conserved binding sites on RNA polymerase and many effects on the transcription elongation complex but differ in their mechanisms of recruitment, interactions with nucleic acids and secondary partners, and regulatory outcomes. Strikingly, some can alternate between autoinhibited and activated states that possess dramatically different secondary structures to achieve exquisite target specificity.

Flipping states: a few key residues decide the winning conformation of the only universally conserved transcription factor

Transcription factors from the NusG family bind to the elongating RNA polymerase to enable synthesis of long RNAs in all domains of life. In bacteria, NusG frequently co-exists with specialized paralogs that regulate expression of a small set of targets, many of which encode virulence factors. Escherichia coli RfaH is the exemplar of this regulatory mechanism. In contrast to NusG, which freely binds to RNA polymerase, RfaH exists in a structurally distinct autoinhibitory state in which the RNA polymerase-binding site is buried at the interface between two RfaH domains. Binding to an ops DNA sequence triggers structural transformation wherein the domains dissociate and RfaH refolds into a NusG-like structure. Formation of the autoinhibitory state, and thus sequence-specific recruitment, represents the decisive step in the evolutionary history of the RfaH subfamily. We used computational and experimental approaches to identify the residues that confer the unique regulatory properties of RfaH. Our analysis highlighted highly conserved Ile and Phe residues at the RfaH interdomain interface. Replacement of these residues with equally conserved Glu and Val counterpart residues in NusG destabilized interactions between the RfaH domains and allowed sequence-independent recruitment to RNA polymerase, suggesting a plausible pathway for diversification of NusG paralogs.