Structural Basis for Transcript Elongation Control by NusG Family Universal Regulators

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.

The yin and yang of the universal transcription factor NusG

RNA polymerase (RNAP), the central enzyme of transcription, intermittently pauses during the elongation stage of RNA synthesis. Pausing provides an opportunity for regulatory events such as nascent RNA folding or the recruitment of transregulators. NusG (Spt5 in eukaryotes and archaea) regulates RNAP pausing and is the only transcription factor conserved across all cellular life. NusG is a multifunctional protein: its N-terminal domain (NGN) binds to RNAP, and its C-terminal KOW domain in bacteria interacts with transcription regulators such as ribosomes and termination factors. In Escherichia coli, NusG acts as an antipausing factor. However, recent studies have revealed that NusG has distinct transcriptional regulatory roles specific to bacterial clades with clinical implications. Here, we focus on NusG's dual roles in the regulation of pausing.

Discovery of Inhibitors Targeting Protein-Protein Interaction between Bacterial RNA Polymerase and NusG as Novel Antimicrobials

Bacterial RNA polymerase (RNAP), the core enzyme responsible for bacterial transcription, requires the NusG factor for efficient transcription elongation and termination. As the primary binding site for NusG, the RNAP clamp-helix (CH) domain represents a potential protein-protein interaction (PPI) target for novel antimicrobial agent design and discovery. In this study, we designed a pharmacophore model based on the essential amino acids of the CH for binding to NusG, such as R270, R278, and R281 (Escherichia coli numbering), and identified a hit compound with mild antimicrobial activity. Subsequent rational design and synthesis of this hit compound led to improved antimicrobial activity against Streptococcus pneumoniae, with the minimum inhibitory concentration (MIC) reduced from 128 to 1 μg/mL. Additional characterization of the antimicrobial activity, inhibitory activity against RNAP-NusG interaction, and cell-based transcription and fluorescent assays of the optimized compounds demonstrated their potential for further lead optimization.

Structural basis of archaeal FttA-dependent transcription termination

The ribonuclease FttA (also known as aCPSF and aCPSF1) mediates factor-dependent transcription termination in archaea1-3. Here we report the structure of a Thermococcus kodakarensis transcription pre-termination complex comprising FttA, Spt4, Spt5 and a transcription elongation complex (TEC). The structure shows that FttA interacts with the TEC in a manner that enables RNA to proceed directly from the TEC RNA-exit channel to the FttA catalytic centre and that enables endonucleolytic cleavage of RNA by FttA, followed by 5'→3' exonucleolytic cleavage of RNA by FttA and concomitant 5'→3' translocation of FttA on RNA, to apply mechanical force to the TEC and trigger termination. The structure further reveals that Spt5 bridges FttA and the TEC, explaining how Spt5 stimulates FttA-dependent termination. The results reveal functional analogy between bacterial and archaeal factor-dependent termination, functional homology between archaeal and eukaryotic factor-dependent termination, and fundamental mechanistic similarities in factor-dependent termination in bacteria, archaea, and eukaryotes.

Antiterminator LoaP loads onto RNA to chase a runaway RNA polymerase

In this issue of Structure, Elghondakly et al.1 present the crystal structure of Thermoanaerobacter pseudethanolicus antiterminator LoaP, a member of a ubiquitous family of NusG transcription factors, bound to its target, a dfn RNA hairpin. LoaP uses RNA as a recognition determinant, which is unique among NusG paralogs and makes unusual contacts in the major groove of the RNA.

Major-groove sequence-specific RNA recognition by LoaP, a paralog of transcription elongation factor NusG

LoaP is a member of the universal NusG protein family. Previously, we reported that unlike other characterized homologs, LoaP binds RNA sequence-specifically, recognizing a stem-loop in the 5'-untranslated region of operons it regulates. To elucidate how this NusG homolog acquired this ability, we now determined the co-crystal structure of Thermoanaerobacter pseudethanolicus LoaP bound to its cognate 26-nucleotide dfn RNA element. Our structure reveals that the LoaP C-terminal KOW domain recognizes the helical portion of the RNA by docking into a broadened major groove, while a protruding β-hairpin of the N-terminal NusG-like domain binds the UNCG tetraloop capping the stem-loop. Major-groove RNA recognition is unusual and is made possible by conserved features of the dfn hairpin. Superposition with structures of other NusG proteins implies that LoaP can bind concurrently to the dfn RNA and the transcription elongation complex, suggesting a new level of co-transcriptional regulation by proteins of this conserved family.

Pause Patrol: Negative Elongation Factor's Role in Promoter-Proximal Pausing and Beyond

RNA polymerase (Pol) II is highly regulated to ensure appropriate gene expression. Early transcription elongation is associated with transient pausing of RNA Pol II in the promoter-proximal region. In multicellular organisms, this pausing is stabilized by the association of transcription elongation factors DRB-sensitivity inducing factor (DSIF) and Negative Elongation Factor (NELF). DSIF is a broadly conserved transcription elongation factor whereas NELF is mostly restricted to the metazoan lineage. Mounting evidence suggests that NELF association with RNA Pol II serves as checkpoint for either release into rapid and productive transcription elongation or premature termination at promoter-proximal pause sites. Here we summarize NELF's roles in promoter-proximal pausing, transcription termination, DNA repair, and signaling based on decades of cell biological, biochemical, and structural work and describe areas for future research.

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.

Structural basis of RfaH-mediated transcription-translation coupling

The NusG paralog RfaH mediates bacterial transcription-translation coupling in genes that contain a DNA sequence element, termed an ops site, required for pausing RNA polymerase (RNAP) and for loading RfaH onto the paused RNAP. Here, we report cryo-electron microscopy structures of transcription-translation complexes (TTCs) containing Escherichia coli RfaH. The results show that RfaH bridges RNAP and the ribosome, with the RfaH N-terminal domain interacting with RNAP and the RfaH C-terminal domain interacting with the ribosome. The results show that the distribution of translational and orientational positions of RNAP relative to the ribosome in RfaH-coupled TTCs is more restricted than in NusG-coupled TTCs because of the more restricted flexibility of the RfaH interdomain linker. The results further suggest that the structural organization of RfaH-coupled TTCs in the 'loading state', in which RNAP and RfaH are located at the ops site during formation of the TTC, is the same as the structural organization of RfaH-coupled TTCs in the 'loaded state', in which RNAP and RfaH are located at positions downstream of the ops site during function of the TTC. The results define the structural organization of RfaH-containing TTCs and set the stage for analysis of functions of RfaH during translation initiation and transcription-translation coupling.

Bacteroides expand the functional versatility of a universal transcription factor and transcribed DNA to program capsule diversity

Human gut Bacteroides species encode numerous (eight or more) tightly regulated capsular polysaccharides (CPS). Specialized paralogs of the universal transcription elongation factor NusG, called UpxY (Y), and an anti-Y UpxZ (Z) are encoded by the first two genes of each CPS operon. The Y-Z regulators combine with promoter inversions to limit CPS transcription to a single operon in most cells. Y enhances transcript elongation whereas Z inhibits noncognate Ys. How Y distinguishes among cognate CPS operons and how Z inhibits only noncognate Ys are unknown. Using in-vivo nascent-RNA sequencing and promoter-less in vitro transcription (PIVoT), we establish that Y recognizes a paused RNA polymerase via sequences in both the exposed non-template DNA and the upstream duplex DNA. Y association is aided by novel 'pause-then-escape' nascent RNA hairpins. Z binds non-cognate Ys to directly inhibit Y association. This Y-Z hierarchical regulatory program allows Bacteroides to create CPS subpopulations for optimal fitness.

Structural basis of archaeal RNA polymerase transcription elongation and Spt4/5 recruitment

Archaeal transcription is carried out by a multi-subunit RNA polymerase (RNAP) that is highly homologous in structure and function to eukaryotic RNAP II. Among the set of basal transcription factors, only Spt5 is found in all domains of life, but Spt5 has been shaped during evolution, which is also reflected in the heterodimerization of Spt5 with Spt4 in Archaea and Eukaryotes. To unravel the mechanistic basis of Spt4/5 function in Archaea, we performed structure-function analyses using the archaeal transcriptional machinery of Pyrococcus furiosus (Pfu). We report single-particle cryo-electron microscopy reconstructions of apo RNAP and the archaeal elongation complex (EC) in the absence and presence of Spt4/5. Surprisingly, Pfu Spt4/5 also binds the RNAP in the absence of nucleic acids in a distinct super-contracted conformation. We show that the RNAP clamp/stalk module exhibits conformational flexibility in the apo state of RNAP and that the enzyme contracts upon EC formation or Spt4/5 engagement. We furthermore identified a contact of the Spt5-NGN domain with the DNA duplex that stabilizes the upstream boundary of the transcription bubble and impacts Spt4/5 activity in vitro. This study, therefore, provides the structural basis for Spt4/5 function in archaeal transcription and reveals a potential role beyond the well-described support of elongation.

A nascent riboswitch helix orchestrates robust transcriptional regulation through signal integration

Widespread manganese-sensing transcriptional riboswitches effect the dependable gene regulation needed for bacterial manganese homeostasis in changing environments. Riboswitches - like most structured RNAs - are believed to fold co-transcriptionally, subject to both ligand binding and transcription events; yet how these processes are orchestrated for robust regulation is poorly understood. Through a combination of single-molecule and bulk approaches, we discover how a single Mn2+ ion and the transcribing RNA polymerase (RNAP), paused immediately downstream by a DNA template sequence, are coordinated by the bridging switch helix P1.1 in the representative Lactococcus lactis riboswitch. This coordination achieves a heretofore-overlooked semi-docked global conformation of the nascent RNA, P1.1 base pair stabilization, transcription factor NusA ejection, and RNAP pause extension, thereby enforcing transcription readthrough. Our work demonstrates how a central, adaptable RNA helix functions analogous to a molecular fulcrum of a first-class lever system to integrate disparate signals for finely balanced gene expression control.

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

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

Concerted transformation of a hyper-paused transcription complex and its reinforcing protein

RfaH, a paralog of the universally conserved NusG, binds to RNA polymerases (RNAP) and ribosomes to activate expression of virulence genes. In free, autoinhibited RfaH, an α-helical KOW domain sequesters the RNAP-binding site. Upon recruitment to RNAP paused at an ops site, KOW is released and refolds into a β-barrel, which binds the ribosome. Here, we report structures of ops-paused transcription elongation complexes alone and bound to the autoinhibited and activated RfaH, which reveal swiveled, pre-translocated pause states stabilized by an ops hairpin in the non-template DNA. Autoinhibited RfaH binds and twists the ops hairpin, expanding the RNA:DNA hybrid to 11 base pairs and triggering the KOW release. Once activated, RfaH hyper-stabilizes the pause, which thus requires anti-backtracking factors for escape. Our results suggest that the entire RfaH cycle is solely determined by the ops and RfaH sequences and provide insights into mechanisms of recruitment and metamorphosis of NusG homologs across all life.

Distinct negative elongation factor conformations regulate RNA polymerase II promoter-proximal pausing

Metazoan gene expression regulation involves pausing of RNA polymerase (Pol II) in the promoter-proximal region of genes and is stabilized by DSIF and NELF. Upon depletion of elongation factors, NELF appears to accompany elongating Pol II past pause sites; however, prior work indicates that NELF prevents Pol II elongation. Here, we report cryoelectron microscopy structures of Pol II-DSIF-NELF complexes with NELF in two distinct conformations corresponding to paused and poised states. The paused NELF state supports Pol II stalling, whereas the poised NELF state enables transcription elongation as it does not support a tilted RNA-DNA hybrid. Further, the poised NELF state can accommodate TFIIS binding to Pol II, allowing for Pol II reactivation at paused or backtracking sites. Finally, we observe that the NELF-A tentacle interacts with the RPB2 protrusion and is necessary for pausing. Our results define how NELF can support pausing, reactivation, and elongation by Pol II.

Molecular basis of mRNA delivery to the bacterial ribosome

Protein synthesis begins with the formation of a ribosome-mRNA complex. In bacteria, the 30S ribosomal subunit is recruited to many mRNAs through base pairing with the Shine Dalgarno (SD) sequence and RNA binding by ribosomal protein bS1. Translation can initiate on nascent mRNAs and RNA polymerase (RNAP) can promote recruitment of the pioneering 30S subunit. Here we examined ribosome recruitment to nascent mRNAs using cryo-EM, single-molecule fluorescence co-localization, and in-cell crosslinking mass spectrometry. We show that bS1 delivers the mRNA to the ribosome for SD duplex formation and 30S subunit activation. Additionally, bS1 mediates the stimulation of translation initiation by RNAP. Together, our work provides a mechanistic framework for how the SD duplex, ribosomal proteins and RNAP cooperate in 30S recruitment to mRNAs and establish transcription-translation coupling.

Single-molecule tracking reveals the functional allocation, in vivo interactions, and spatial organization of universal transcription factor NusG

During transcription elongation, NusG aids RNA polymerase by inhibiting pausing, promoting anti-termination on rRNA operons, coupling transcription with translation on mRNA genes, and facilitating Rho-dependent termination. Despite extensive work, the in vivo functional allocation and spatial distribution of NusG remain unknown. Using single-molecule tracking and super-resolution imaging in live E. coli cells, we found NusG predominantly in a chromosome-associated population (binding to RNA polymerase in elongation complexes) and a slowly diffusing population complexed with the 30S ribosomal subunit; the latter provides a "30S-guided" path for NusG into transcription elongation. Only ∼10% of NusG is fast diffusing, with its mobility suggesting non-specific interactions with DNA for >50% of the time. Antibiotic treatments and deletion mutants revealed that chromosome-associated NusG participates mainly in rrn anti-termination within phase-separated transcriptional condensates and in transcription-translation coupling. This study illuminates the multiple roles of NusG and offers a guide on dissecting multi-functional machines via in vivo imaging.

A nascent riboswitch helix orchestrates robust transcriptional regulation through signal integration

Widespread manganese-sensing transcriptional riboswitches effect the dependable gene regulation needed for bacterial manganese homeostasis in changing environments. Riboswitches - like most structured RNAs - are believed to fold co-transcriptionally, subject to both ligand binding and transcription events; yet how these processes are orchestrated for robust regulation is poorly understood. Through a combination of single molecule and bulk approaches, we discovered how a single Mn 2+ ion and the transcribing RNA polymerase (RNAP), paused immediately downstream by a DNA template sequence, are coordinated by the bridging switch helix P1.1 in the paradigmatic Lactococcus lactis riboswitch. This coordination achieves a heretofore-overlooked semi-docked global conformation of the nascent RNA, P1.1 base pair stabilization, transcription factor NusA ejection, and RNAP pause extension, thereby enforcing transcription readthrough. Our work demonstrates how a central, adaptable RNA helix functions analogous to a molecular fulcrum of a first-class lever system to integrate disparate signals for finely balanced gene expression control.

Line Shape Analysis of 19F NMR-Monitored Chemical Denaturation of a Fold-Switching Protein RfaH Reveals Its Slow Folding Dynamics

The recent discovery of metamorphic proteins, which can switch between multiple conformations under native conditions, has challenged the well-established one sequence-one structure paradigm of protein folding. This is exemplified in the C-terminal domain of the multidomain transcription factor RfaH, which converts from an α-helical coiled-coil conformation in its autoinhibited state to a β-barrel conformation upon activation. Here, we use multisite line shape analysis of 19F NMR-monitored equilibrium chemical denaturation measurements of two 19F-labeled variants of full-length RfaH, to show that it folds/unfolds slowly on the NMR time scale, in an apparent all-or-none fashion at physiological pH and room temperature in the 3.3-4.8 M urea concentration range. The significant thermodynamic stability and slow unfolding rate (kinetic stability) are likely responsible for maintaining the closed autoinhibited state of RfaH, preventing spurious interactions with RNA polymerase (RNAP) when not functional. Our results provide a quantitative understanding of the folding-function relationship in the model fold-switching protein RfaH.

AlphaFold2 has more to learn about protein energy landscapes

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

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.

Structural basis of RfaH-mediated transcription-translation coupling

The NusG paralog RfaH mediates bacterial transcription-translation coupling on genes that contain a DNA sequence element, termed an ops site, required for pausing RNA polymerase (RNAP) and for loading RfaH onto the paused RNAP. Here we report cryo-EM structures of transcription-translation complexes (TTCs) containing RfaH. The results show that RfaH bridges RNAP and the ribosome, with the RfaH N-terminal domain interacting with RNAP, and with the RfaH C-terminal domain interacting with the ribosome. The results show that the distribution of translational and orientational positions of RNAP relative to the ribosome in RfaH-coupled TTCs is more restricted than in NusG-coupled TTCs, due to the more restricted flexibility of the RfaH interdomain linker. The results further show that the structural organization of RfaH-coupled TTCs in the "loading state," in which RNAP and RfaH are located at the ops site during formation of the TTC, is the same as the structural organization of RfaH-coupled TTCs in the "loaded state," in which RNAP and RfaH are located at positions downstream of the ops site during function of the TTC. The results define the structural organization of RfaH-containing TTCs and set the stage for analysis of functions of RfaH during translation initiation and transcription-translation coupling.

One sentence summary

Cryo-EM reveals the structural basis of transcription-translation coupling by RfaH.

Computing Free Energies of Fold-Switching Proteins Using MELD x MD

Some proteins are conformational switches, able to transition between relatively different conformations. To understand what drives them requires computing the free-energy difference ΔGAB between their stable states, A and B. Molecular dynamics (MD) simulations alone are often slow because they require a reaction coordinate and must sample many transitions in between. Here, we show that modeling employing limited data (MELD) x MD on known endstates A and B is accurate and efficient because it does not require passing over barriers or knowing reaction coordinates. We validate this method on two problems: (1) it gives correct relative populations of α and β conformers for small designed chameleon sequences of protein G; and (2) it correctly predicts the conformations of the C-terminal domain (CTD) of RfaH. Free-energy methods like MELD x MD can often resolve structures that confuse machine-learning (ML) methods.

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.

Exploring the structural acrobatics of fold-switching proteins using simplified structure-based models

Metamorphic proteins are a paradigm of the protein folding process, by encoding two or more native states, highly dissimilar in terms of their secondary, tertiary, and even quaternary structure, on a single amino acid sequence. Moreover, these proteins structurally interconvert between these native states in a reversible manner at biologically relevant timescales as a result of different environmental cues. The large-scale rearrangements experienced by these proteins, and their sometimes high mass interacting partners that trigger their metamorphosis, makes the computational and experimental study of their structural interconversion challenging. Here, we present our efforts in studying the refolding landscapes of two quintessential metamorphic proteins, RfaH and KaiB, using simplified dual-basin structure-based models (SBMs), rigorously footed on the energy landscape theory of protein folding and the principle of minimal frustration. By using coarse-grained models in which the native contacts and bonded interactions extracted from the available experimental structures of the two native states of RfaH and KaiB are merged into a single Hamiltonian, dual-basin SBM models can be generated and savvily calibrated to explore their fold-switch in a reversible manner in molecular dynamics simulations. We also describe how some of the insights offered by these simulations have driven the design of experiments and the validation of the conformational ensembles and refolding routes observed using this simple and computationally efficient models.

Structural basis for control of bacterial RNA polymerase pausing by a riboswitch and its ligand

Folding of nascent transcripts can be modulated by the RNA polymerase (RNAP) that carries out their transcription, and vice versa. A pause of RNAP during transcription of a preQ1 riboswitch (termed que-PEC) is stabilized by a previously characterized template consensus sequence and the ligand-free conformation of the nascent RNA. Ligand binding to the riboswitch induces RNAP pause release and downstream transcription termination; however, the mechanism by which riboswitch folding modulates pausing is unclear. Here, we report single-particle cryo-electron microscopy reconstructions of que-PEC in ligand-free and ligand-bound states. In the absence of preQ1, the RNA transcript is in an unexpected hyper-translocated state, preventing downstream nucleotide incorporation. Strikingly, on ligand binding, the riboswitch rotates around its helical axis, expanding the surrounding RNAP exit channel and repositioning the transcript for elongation. Our study reveals the tight coupling by which nascent RNA structures and their ligands can functionally regulate the macromolecular transcription machinery.

Robust regulation of transcription pausing in Escherichia coli by the ubiquitous elongation factor NusG

Transcription elongation by multi-subunit RNA polymerases (RNAPs) is regulated by auxiliary factors in all organisms. NusG/Spt5 is the only universally conserved transcription elongation factor shared by all domains of life. NusG is a component of antitermination complexes controlling ribosomal RNA operons, an essential antipausing factor, and a transcription-translation coupling factor in Escherichia coli. We employed RNET-seq for genome-wide mapping of RNAP pause sites in wild-type and NusG-depleted cells. We demonstrate that NusG is a major antipausing factor that suppresses thousands of backtracked and nonbacktracked pauses across the E. coli genome. The NusG-suppressed pauses were enriched immediately downstream from the translation start codon but were also abundant elsewhere in open reading frames, small RNA genes, and antisense transcription units. This finding revealed a strong similarity of NusG to Spt5, which stimulates the elongation rate of many eukaryotic genes. We propose a model in which promoting forward translocation and/or stabilization of RNAP in the posttranslocation register by NusG results in suppression of pausing in E. coli.

When push comes to shove - RNA polymerase and DNA-bound protein roadblocks

In recent years, transcriptional roadblocking has emerged as a crucial regulatory mechanism in gene expression, whereby other DNA-bound obstacles can block the progression of transcribing RNA polymerase (RNAP), leading to RNAP pausing and ultimately dissociation from the DNA template. In this review, we discuss the mechanisms by which transcriptional roadblocks can impede RNAP progression, as well as how RNAP can overcome these obstacles to continue transcription. We examine different DNA-binding proteins involved in transcriptional roadblocking and their biophysical properties that determine their effectiveness in blocking RNAP progression. The catalytically dead CRISPR-Cas (dCas) protein is used as an example of an engineered programmable roadblock, and the current literature in understanding the polarity of dCas roadblocking is also discussed. Finally, we delve into a stochastic model of transcriptional roadblocking and highlight the importance of transcription factor binding kinetics and its resistance to dislodgement by an elongating RNAP in determining the strength of a roadblock.

Structural and functional basis of the universal transcription factor NusG pro-pausing activity in Mycobacterium tuberculosis

Transcriptional pauses mediate regulation of RNA biogenesis. DNA-encoded pause signals trigger pausing by stabilizing RNA polymerase (RNAP) swiveling and inhibiting DNA translocation. The N-terminal domain (NGN) of the only universal transcription factor, NusG/Spt5, modulates pausing through contacts to RNAP and DNA. Pro-pausing NusGs enhance pauses, whereas anti-pausing NusGs suppress pauses. Little is known about pausing and NusG in the human pathogen Mycobacterium tuberculosis (Mtb). We report that MtbNusG is pro-pausing. MtbNusG captures paused, swiveled RNAP by contacts to the RNAP protrusion and nontemplate-DNA wedged between the NGN and RNAP gate loop. In contrast, anti-pausing Escherichia coli (Eco) NGN contacts the MtbRNAP gate loop, inhibiting swiveling and pausing. Using CRISPR-mediated genetics, we show that pro-pausing NGN is required for mycobacterial fitness. Our results define an essential function of mycobacterial NusG and the structural basis of pro- versus anti-pausing NusG activity, with broad implications for the function of all NusG orthologs.

Correlative Escherichia coli Transcription Rate and Bubble Conformation Remodeled by NusA and NusG

Transcription is highly regulated by a variety of transcription factors, among which NusA and NusG act contradictorily in Escherichia coli (E. coli) that NusA stabilizes a paused RNA polymerase (RNAP) and NusG suppresses it. The mechanism of the NusA and NusG regulations on RNAP transcription has been addressed, but their effect on the conformational changes of the transcription bubble correlated with transcription kinetics remains elusive. By using single-molecule magnetic trap, we identify a reduction in the transcription rate of ∼40% events by NusA. Although the rest ∼60% of transcription events exhibit unaffected transcription rates, a NusA-enhanced standard deviation of the transcription rate is observed. NusA remodeling also increases the extent of DNA unwinding in the transcription bubble by 1-2 base pairs, which can be reduced by NusG. The NusG remodeling is more significant on the RNAP molecules with reduced transcription rates rather than those without. Our results provide a quantitative view on the mechanisms of transcriptional regulation by NusA and NusG factors.

A trailing ribosome speeds up RNA polymerase at the expense of transcript fidelity via force and allostery

In prokaryotes, translation can occur on mRNA that is being transcribed in a process called coupling. How the ribosome affects the RNA polymerase (RNAP) during coupling is not well understood. Here, we reconstituted the E. coli coupling system and demonstrated that the ribosome can prevent pausing and termination of RNAP and double the overall transcription rate at the expense of fidelity. Moreover, we monitored single RNAPs coupled to ribosomes and show that coupling increases the pause-free velocity of the polymerase and that a mechanical assisting force is sufficient to explain the majority of the effects of coupling. Also, by cryo-EM, we observed that RNAPs with a terminal mismatch adopt a backtracked conformation, while a coupled ribosome allosterically induces these polymerases toward a catalytically active anti-swiveled state. Finally, we demonstrate that prolonged RNAP pausing is detrimental to cell viability, which could be prevented by polymerase reactivation through a coupled ribosome.

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.

An ensemble of interconverting conformations of the elemental paused transcription complex creates regulatory options

Transcriptional pausing underpins the regulation of cellular RNA synthesis, but its mechanism remains incompletely understood. Sequence-specific interactions of DNA and RNA with the dynamic, multidomain RNA polymerase (RNAP) trigger reversible conformational changes at pause sites that temporarily interrupt the nucleotide addition cycle. These interactions initially rearrange the elongation complex (EC) into an elemental paused EC (ePEC). ePECs can form longer-lived PECs by further rearrangements or interactions of diffusible regulators. For both bacterial and mammalian RNAPs, a half-translocated state in which the next DNA template base fails to load into the active site appears central to the ePEC. Some RNAPs also swivel interconnected modules that may stabilize the ePEC. However, it is unclear whether swiveling and half-translocation are requisite features of a single ePEC state or if multiple ePEC states exist. Here, we use cryo-electron microscopy (cryo-EM) analysis of ePECs with different RNA-DNA sequences combined with biochemical probes of ePEC structure to define an interconverting ensemble of ePEC states. ePECs occupy either pre- or half-translocated states but do not always swivel, indicating that difficulty in forming the posttranslocated state at certain RNA-DNA sequences may be the essence of the ePEC. The existence of multiple ePEC conformations has broad implications for transcriptional regulation.

Allosteric mechanism of transcription inhibition by NusG-dependent pausing of RNA polymerase

NusG is a transcription elongation factor that stimulates transcription pausing in Gram+ bacteria including B. subtilis by sequence-specific interaction with a conserved pause-inducing _-11TTNTTT_-6 motif found in the non-template DNA (ntDNA) strand within the transcription bubble. To reveal the structural basis of NusG-dependent pausing, we determined a cryo-EM structure of a paused transcription complex (PTC) containing RNA polymerase (RNAP), NusG, and the TTNTTT motif in the ntDNA strand. The interaction of NusG with the ntDNA strand rearranges the transcription bubble by positioning three consecutive T residues in a cleft between NusG and the β-lobe domain of RNAP. We revealed that the RNAP swivel module rotation (swiveling), which widens (swiveled state) and narrows (non-swiveled state) a cleft between NusG and the β-lobe, is an intrinsic motion of RNAP and is directly linked to trigger loop (TL) folding, an essential conformational change of all cellular RNAPs for the RNA synthesis reaction. We also determined cryo-EM structures of RNAP escaping from the paused transcription state. These structures revealed the NusG-dependent pausing mechanism by which NusG-ntDNA interaction inhibits the transition from swiveled to non-swiveled states, thereby preventing TL folding and RNA synthesis allosterically. This motion is also reduced by the formation of an RNA hairpin within the RNA exit channel. Thus, the pause half-life can be modulated by the strength of the NusG-ntDNA interaction and/or the stability of the RNA hairpin. NusG residues that interact with the TTNTTT motif are widely conserved in bacteria, suggesting that NusG-dependent pausing is widespread.

Structural basis of the transcription termination factor Rho engagement with transcribing RNA polymerase from Thermus thermophilus

Transcription termination is an essential step in transcription by RNA polymerase (RNAP) and crucial for gene regulation. For many bacterial genes, transcription termination is mediated by the adenosine triphosphate-dependent RNA translocase/helicase Rho, which causes RNA/DNA dissociation from the RNAP elongation complex (EC). However, the structural basis of the interplay between Rho and RNAP remains obscure. Here, we report the cryo-electron microscopy structure of the Thermus thermophilus RNAP EC engaged with Rho. The Rho hexamer binds RNAP through the carboxyl-terminal domains, which surround the RNA exit site of RNAP, directing the nascent RNA seamlessly from the RNA exit to its central channel. The β-flap tip at the RNA exit is critical for the Rho-dependent RNA release, and its deletion causes an alternative Rho-RNAP binding mode, which is irrelevant to termination. The Rho binding site overlaps with the binding sites of other macromolecules, such as ribosomes, providing a general basis of gene regulation.

Structural basis of Rho-dependent transcription termination

Rho is a ring-shaped hexameric ATP-dependent molecular motor. Together with the transcription elongation factor NusG, Rho mediates factor-dependent transcription termination and transcription-translation-coupling quality control in Escherichia coli1-4. Here we report the preparation of complexes that are functional in factor-dependent transcription termination from Rho, NusG, RNA polymerase (RNAP), and synthetic nucleic acid scaffolds, and we report cryogenic electron microscopy structures of the complexes. The structures show that functional factor-dependent pre-termination complexes contain a closed-ring Rho hexamer; have RNA threaded through the central channel of Rho; have 60 nucleotides of RNA interacting sequence-specifically with the exterior of Rho and 6 nucleotides of RNA interacting sequence-specifically with the central channel of Rho; have Rho oriented relative to RNAP such that ATP-dependent translocation by Rho exerts mechanical force on RNAP; and have NusG bridging Rho and RNAP. The results explain five decades of research on Rho and provide a foundation for understanding Rho's function.

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.

Transcription-Translation Coupling in Bacteria

In bacteria, transcription and translation take place in the same cellular compartment. Therefore, a messenger RNA can be translated as it is being transcribed, a process known as transcription-translation coupling. This process was already recognized at the dawn of molecular biology, yet the interplay between the two key players, the RNA polymerase and ribosome, remains elusive. Genetic data indicate that an RNA sequence can be translated shortly after it has been transcribed. The closer both processes are in time, the less accessible the RNA sequence is between the RNA polymerase and ribosome. This temporal coupling has important consequences for gene regulation. Biochemical and structural studies have detailed several complexes between the RNA polymerase and ribosome. The in vivo relevance of this physical coupling has not been formally demonstrated. We discuss how both temporal and physical coupling may mesh to produce the phenomenon we know as transcription-translation coupling.

Structural basis of transcriptional regulation by a nascent RNA element, HK022 putRNA

Transcription, in which RNA polymerases (RNAPs) produce RNA from DNA, is the first step of gene expression. As such, it is highly regulated either by trans-elements like protein factors and/or by cis-elements like specific sequences on the DNA. Lambdoid phage HK022 contains a cis-element, put, which suppresses pausing and termination during transcription of the early phage genes. The putRNA transcript solely performs the anti-pausing/termination activities by interacting directly with the E.coli RNAP elongation complex (EC) by an unknown structural mechanism. In this study, we reconstituted putRNA-associated ECs and determined the structures using cryo-electron microscopy. The determined structures of putRNA-associated EC, putRNA-absent EC, and σ⁷⁰-bound EC suggest that the putRNA interaction with the EC counteracts swiveling, a conformational change previously identified to promote pausing and σ⁷⁰ might modulate putRNA folding via σ⁷⁰-dependent pausing during elongation.

The structure and activities of the archaeal transcription termination factor Eta detail vulnerabilities of the transcription elongation complex

Transcription must be properly regulated to ensure dynamic gene expression underlying growth, development, and response to environmental cues. Regulation is imposed throughout the transcription cycle, and while many efforts have detailed the regulation of transcription initiation and early elongation, the termination phase of transcription also plays critical roles in regulating gene expression. Transcription termination can be driven by only a few proteins in each domain of life. Detailing the mechanism(s) employed provides insight into the vulnerabilities of transcription elongation complexes (TECs) that permit regulated termination to control expression of many genes and operons. Here, we describe the biochemical activities and crystal structure of the superfamily 2 helicase Eta, one of two known factors capable of disrupting archaeal transcription elongation complexes. Eta retains a twin-translocase core domain common to all superfamily 2 helicases and a well-conserved C terminus wherein individual amino acid substitutions can critically abrogate termination activities. Eta variants that perturb ATPase, helicase, single-stranded DNA and double-stranded DNA translocase and termination activities identify key regions of the C terminus of Eta that, when combined with modeling Eta-TEC interactions, provide a structural model of Eta-mediated termination guided in part by structures of Mfd and the bacterial TEC. The susceptibility of TECs to disruption by termination factors that target the upstream surface of RNA polymerase and potentially drive termination through forward translocation and allosteric mechanisms that favor opening of the clamp to release the encapsulated nucleic acids emerges as a common feature of transcription termination mechanisms.

Structural advances in transcription elongation

Transcription is the first step of gene expression and involves RNA polymerases. After transcription initiation, RNA polymerase enters elongation followed by transcription termination at the end of the gene. Only recently, structures of transcription elongation complexes bound to key transcription elongation factors have been determined in bacterial and eukaryotic systems. These structures have revealed numerous insights including the basis for transcriptional pausing, RNA polymerase interaction with large complexes such as the ribosome and the spliceosome, and the transition into productive elongation. Here, we review these structures and describe areas for future research.

The pleiotropic roles of SPT5 in transcription

Initially discovered by genetic screens in budding yeast, SPT5 and its partner SPT4 form a stable complex known as DSIF in metazoa, which plays pleiotropic roles in multiple steps of transcription. SPT5 is the most conserved transcription elongation factor, being found in all three domains of life; however, its structure has evolved to include new domains and associated posttranslational modifications. These gained features have expanded transcriptional functions of SPT5, likely to meet the demand for increasingly complex regulation of transcription in higher organisms. This review discusses the pleiotropic roles of SPT5 in transcription, including RNA polymerase II (Pol II) stabilization, enhancer activation, Pol II pausing and its release, elongation, and termination, with a focus on the most recent progress of SPT5 functions in regulating metazoan transcription.

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.

AtNusG, a chloroplast nucleoid protein of bacterial origin linking chloroplast transcriptional and translational machineries, is required for proper chloroplast gene expression in Arabidopsis thaliana

In Escherichia coli, transcription-translation coupling is mediated by NusG. Although chloroplasts are descendants of endosymbiotic prokaryotes, the mechanism underlying this coupling in chloroplasts remains unclear. Here, we report transcription-translation coupling through AtNusG in chloroplasts. AtNusG is localized in chloroplast nucleoids and is closely associated with the chloroplast PEP complex by interacting with its essential component PAP9. It also comigrates with chloroplast ribosomes and interacts with their two components PRPS5 (uS5c) and PRPS10 (uS10c). These data suggest that the transcription and translation machineries are coupled in chloroplasts. In the atnusg mutant, the accumulation of chloroplast-encoded photosynthetic gene transcripts, such as psbA, psbB, psbC and psbD, was not obviously changed, but that of their proteins was clearly decreased. Chloroplast polysomic analysis indicated that the decrease in these proteins was due to the reduced efficiency of their translation in this mutant, leading to reduced photosynthetic efficiency and enhanced sensitivity to cold stress. These data indicate that AtNusG-mediated coupling between transcription and translation in chloroplasts ensures the rapid establishment of photosynthetic capacity for plant growth and the response to environmental changes. Therefore, our study reveals a conserved mechanism of transcription-translation coupling between chloroplasts and E. coli, which perhaps represents a regulatory mechanism of chloroplast gene expression. This study provides insights into the underlying mechanisms of chloroplast gene expression in higher plants.

Allosteric couplings upon binding of RfaH to transcription elongation complexes

In every domain of life, NusG-like proteins bind to the elongating RNA polymerase (RNAP) to support processive RNA synthesis and to couple transcription to ongoing cellular processes. Structures of factor-bound transcription elongation complexes (TECs) reveal similar contacts to RNAP, consistent with a shared mechanism of action. However, NusG homologs differ in their regulatory roles, modes of recruitment, and effects on RNA synthesis. Some of these differences could be due to conformational changes in RNAP and NusG-like proteins, which cannot be captured in static structures. Here, we employed hydrogen-deuterium exchange mass spectrometry to investigate changes in local and non-local structural dynamics of Escherichia coli NusG and its paralog RfaH, which have opposite effects on expression of xenogenes, upon binding to TEC. We found that NusG and RfaH regions that bind RNAP became solvent-protected in factor-bound TECs, whereas RNAP regions that interact with both factors showed opposite deuterium uptake changes when bound to NusG or RfaH. Additional changes far from the factor-binding site were observed only with RfaH. Our results provide insights into differences in structural dynamics exerted by NusG and RfaH during binding to TEC, which may explain their different functional outcomes and allosteric regulation of transcriptional pausing by RfaH.

Structural and mechanistic basis of σ-dependent transcriptional pausing

In σ-dependent transcriptional pausing, the transcription initiation factor σ, translocating with RNA polymerase (RNAP), makes sequence-specific protein–DNA interactions with a promoter-like sequence element in the transcribed region, inducing pausing. It has been proposed that, in σ-dependent pausing, the RNAP active center can access off-pathway “backtracked” states that are substrates for the transcript-cleavage factors of the Gre family and on-pathway “scrunched” states that mediate pause escape. Here, using site-specific protein–DNA photocrosslinking to define positions of the RNAP trailing and leading edges and of σ relative to DNA at the λPR′ promoter, we show directly that σ-dependent pausing in the absence of GreB in vitro predominantly involves a state backtracked by 2–4 bp, and σ-dependent pausing in the presence of GreB in vitro and in vivo predominantly involves a state scrunched by 2–3 bp. Analogous experiments with a library of 47 (∼16,000) transcribed-region sequences show that the state scrunched by 2–3 bp—and only that state—is associated with the consensus sequence, T−3N−2Y−1G+1, (where −1 corresponds to the position of the RNA 3′ end), which is identical to the consensus for pausing in initial transcription and which is related to the consensus for pausing in transcription elongation. Experiments with heteroduplex templates show that sequence information at position T−3 resides in the DNA nontemplate strand. A cryoelectron microscopy structure of a complex engaged in σ-dependent pausing reveals positions of DNA scrunching on the DNA nontemplate and template strands and suggests that position T−3 of the consensus sequence exerts its effects by facilitating scrunching.

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

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