A contact-based analysis of local energetic frustration dynamics identifies key residues enabling RfaH fold-switch.

Fold-switching enables metamorphic proteins to reversibly interconvert between two highly dissimilar native states to regulate their protein functions. While about 100 proteins have been identified to undergo fold-switching, unveiling the key residues behind this mechanism for each protein remains challenging. Reasoning that fold-switching in proteins is driven by dynamic changes in local energetic frustration, we combined fold-switching simulations generated using simplified structure-based models with frustration analysis to identify key residues involved in this process based on the change in the density of minimally frustrated contacts during refolding. Using this approach to analyze the fold-switch of the bacterial transcription factor RfaH, we identified 20 residues that significantly change their frustration during its fold-switch, some of which have been experimentally and computationally reported in previous works. Our approach, which we developed as an additional module for the FrustratometeR package, highlights the role of local frustration dynamics in protein fold-switching and offers a robust tool to enhance our understanding of other proteins with significant conformational shifts.

Dissecting the structural and functional consequences of the evolutionary proline-glycine deletion in the wing 1 region of the forkhead domain of human FoxP1.

The human FoxP transcription factors dimerize via three-dimensional domain swapping, a unique feature among the human Fox family, as result of evolutionary sequence adaptations in the forkhead domain. This is the case for the conserved glycine and proline residues in the wing 1 region, which are absent in FoxP proteins but present in most of the Fox family. In this work, we engineered both glycine (G) and proline-glycine (PG) insertion mutants to evaluate the deletion events in FoxP proteins in their dimerization, stability, flexibility, and DNA-binding ability. We show that the PG insertion only increases protein stability, whereas the single glycine insertion decreases the association rate and protein stability and promotes affinity to the DNA ligand.

The Light Chain Allosterically Enhances the Protease Activity of Murine Urokinase-Type Plasminogen Activator.

The active form of the murine urokinase-type plasminogen activator (muPA) is formed by a 27-residue disordered light chain connecting the amino-terminal fragment (ATF) with the serine protease domain. The two chains are tethered by a disulfide bond between C1CT in the disordered light chain and C122CT in the protease domain. Previous work showed that the presence of the disordered light chain affected the inhibition of the protease domain by antibodies. Here we show that the disordered light chain induced a 3.7-fold increase in kcat of the protease domain of muPA. In addition, hydrogen-deuterium exchange mass spectrometry (HDX-MS) and accelerated molecular dynamics (AMD) were performed to identify the interactions between the disordered light chain and the protease domain. HDX-MS revealed that the light chain is contacting the 110s, the turn between the ß10- and ß11-strand, and the ß7-strand. A reduction in deuterium uptake was also observed in the activation loop, the 140s loop and the 220s loop, which forms the S1-specificty pocket where the substrate binds. These loops are further away from where the light chain seems to be interacting with the protease domain. Our results suggest that the light chain most likely increases the activity of muPA by allosterically favoring conformations in which the specificity pocket is formed. We propose a model by which the allostery would be transmitted through the ß-strands of the ß-barrels to the loops on the other side of the protease domain.

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.

DNA controls the dimerization of the human FoxP1 forkhead domain.

Transcription factors (TFs) regulate gene expression by binding to specific DNA sequences and gating access to genes. Even when the binding of TFs and their cofactors to DNA is reversible, indicating a reversible control of gene expression, there is little knowledge about the molecular effect DNA has on TFs. Using single-molecule multiparameter fluorescence spectroscopy, molecular dynamics simulations, and biochemical assays, we find that the monomeric form of the forkhead (FKH) domain of the human FoxP1 behaves as a disordered protein and increases its folded population when it dimerizes. Notably, DNA binding promotes a disordered FKH dimer bound to DNA, negatively controlling the stability of the dimeric FoxP1:DNA complex. The DNA-mediated reversible regulation on FKH dimers suggests that FoxP1-dependent gene suppression is unstable, and it must require the presence of other dimerization domains or cofactors to revert the negative impact exerted by the DNA.

Metamorphic proteins and how to find them.

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

An Open One-Step RT-qPCR for SARS-CoV-2 detection.

The COVID-19 pandemic has resulted in millions of deaths globally, and while several diagnostic systems were proposed, real-time reverse transcription polymerase chain reaction (RT-PCR) remains the gold standard. However, diagnostic reagents, including enzymes used in RT-PCR, are subject to centralized production models and intellectual property restrictions, which present a challenge for less developed countries. With the aim of generating a standardized One-Step open RT-qPCR protocol to detect SARS-CoV-2 RNA in clinical samples, we purified and tested recombinant enzymes and a non-proprietary buffer. The protocol utilized M-MLV RT and Taq DNA pol enzymes to perform a Taqman probe-based assay. Synthetic RNA samples were used to validate the One-Step RT-qPCR components, demonstrating sensitivity comparable to a commercial kit routinely employed in clinical settings for patient diagnosis. Further evaluation on 40 clinical samples (20 positive and 20 negative) confirmed its comparable diagnostic accuracy. This study represents a proof of concept for an open approach to developing diagnostic kits for viral infections and diseases, which could provide a cost-effective and accessible solution for less developed countries.

Local energetic frustration conservation in protein families and superfamilies.

Energetic local frustration offers a biophysical perspective to interpret the effects of sequence variability on protein families. Here we present a methodology to analyze local frustration patterns within protein families and superfamilies that allows us to uncover constraints related to stability and function, and identify differential frustration patterns in families with a common ancestry. We analyze these signals in very well studied protein families such as PDZ, SH3, ɑ and ß globins and RAS families. Recent advances in protein structure prediction make it possible to analyze a vast majority of the protein space. An automatic and unsupervised proteome-wide analysis on the SARS-CoV-2 virus demonstrates the potential of our approach to enhance our understanding of the natural phenotypic diversity of protein families beyond single protein instances. We apply our method to modify biophysical properties of natural proteins based on their family properties, as well as perform unsupervised analysis of large datasets to shed light on the physicochemical signatures of poorly characterized proteins such as the ones belonging to emergent pathogens.

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.

Engineering the catalytic activity of an Antarctic PET-degrading enzyme by loop exchange.

Several hydrolases have been described to degrade polyethylene terephthalate (PET) at moderate temperatures ranging from 25°C to 40°C. These mesophilic PET hydrolases (PETases) are less efficient in degrading this plastic polymer than their thermophilic homologs and have, therefore, been the subject of many protein engineering campaigns. However, enhancing their enzymatic activity through rational design or directed evolution poses a formidable challenge due to the need for exploring a large number of mutations. Additionally, evaluating the improvements in both activity and stability requires screening numerous variants, either individually or using high-throughput screening methods. Here, we utilize instead the design of chimeras as a protein engineering strategy to increase the activity and stability of Mors1, an Antarctic PETase active at 25°C. First, we obtained the crystal structure of Mors1 at 1.6 Å resolution, which we used as a scaffold for structure- and sequence-based chimeric design. Then, we designed a Mors1 chimera via loop exchange of a highly divergent active site loop from the thermophilic leaf-branch compost cutinase (LCC) into the equivalent region in Mors1. After restitution of an active site disulfide bond into this chimera, the enzyme exhibited a shift in optimal temperature for activity to 45°C and an increase in fivefold in PET hydrolysis when compared with wild-type Mors1 at 25°C. Our results serve as a proof of concept of the utility of chimeric design to further improve the activity and stability of PETases active at moderate temperatures.

DNA facilitates heterodimerization between human transcription factors FoxP1 and FoxP2 by increasing their conformational flexibility.

Transcription factors regulate gene expression by binding to DNA. They have disordered regions and specific DNA-binding domains. Binding to DNA causes structural changes, including folding and interactions with other molecules. The FoxP subfamily of transcription factors in humans is unique because they can form heterotypic interactions without DNA. However, it is unclear how they form heterodimers and how DNA binding affects their function. We used computational and experimental methods to study the structural changes in FoxP1's DNA-binding domain when it forms a heterodimer with FoxP2. We found that FoxP1 has complex and diverse conformational dynamics, transitioning between compact and extended states. Surprisingly, DNA binding increases the flexibility of FoxP1, contrary to the typical folding-upon-binding mechanism. In addition, we observed a 3-fold increase in the rate of heterodimerization after FoxP1 binds to DNA. These findings emphasize the importance of structural flexibility in promoting heterodimerization to form transcriptional complexes.

Effect of temperature and nucleotide on the binding of BiP chaperone to a protein substrate.

BiP (immunoglobulin heavy-chain binding protein) is a Hsp70 monomeric ATPase motor that plays broad and crucial roles in maintaining proteostasis inside the cell. Structurally, BiP is formed by two domains, a nucleotide-binding domain (NBD) with ATPase activity connected by a flexible hydrophobic linker to the substrate-binding domain. While the ATPase and substrate binding activities of BiP are allosterically coupled, the latter is also dependent on nucleotide binding. Recent structural studies have provided new insights into BiP's allostery; however, the influence of temperature on the coupling between substrate and nucleotide binding to BiP remains unexplored. Here, we study BiP's binding to its substrate at the single molecule level using thermo-regulated optical tweezers which allows us to mechanically unfold the client protein and explore the effect of temperature and different nucleotides on BiP binding. Our results confirm that the affinity of BiP for its protein substrate relies on nucleotide binding, by mainly regulating the binding kinetics between BiP and its substrate. Interestingly, our findings also showed that the apparent affinity of BiP for its protein substrate in the presence of nucleotides remains invariable over a wide range of temperatures, suggesting that BiP may interact with its client proteins with similar affinities even when the temperature is not optimal. Thus, BiP could play a role as a "thermal buffer" in proteostasis.

Unraveling the folding and dimerization properties of the human FoxP subfamily of transcription factors.

Human FoxP proteins share a highly conserved DNA-binding domain that dimerizes via three-dimensional domain swapping, although showing varying oligomerization propensities among its members. Here, we present an experimental and computational characterization of all human FoxP proteins to unravel how their amino acid substitutions impact their folding and dimerization mechanism. We solved the crystal structure of the forkhead domain of FoxP4 to then perform a comparison across all members, finding that their sequence changes impact not only the structural heterogeneity of their forkhead domains but also the protein-protein association energy barrier. Lastly, we demonstrate that the accumulation of a monomeric intermediate is an oligomerization-dependent feature rather than a common aspect of monomers and dimers in this protein subfamily.

Domain tethering impacts dimerization and DNA-mediated allostery in the human transcription factor FoxP1.

Transcription factors are multidomain proteins with specific DNA binding and regulatory domains. In the human FoxP subfamily (FoxP1, FoxP2, FoxP3, and FoxP4) of transcription factors, a 90 residue-long disordered region links a Leucine Zipper (ZIP)-known to form coiled-coil dimers-and a Forkhead (FKH) domain-known to form domain swapping dimers. We used replica exchange discrete molecular dynamics simulations, single-molecule fluorescence experiments, and other biophysical tools to understand how domain tethering in FoxP1 impacts dimerization at ZIP and FKH domains and how DNA binding allosterically regulates their dimerization. We found that domain tethering promotes FoxP1 dimerization but inhibits a FKH domain-swapped structure. Furthermore, our findings indicate that the linker mediates the mutual organization and dynamics of ZIP and FKH domains, forming closed and open states with and without interdomain contacts, thus highlighting the role of the linkers in multidomain proteins. Finally, we found that DNA allosterically promotes structural changes that decrease the dimerization propensity of FoxP1. We postulate that, upon DNA binding, the interdomain linker plays a crucial role in the gene regulatory function of FoxP1.

Conformational Selection of a Tryptophan Side Chain Drives the Generalized Increase in Activity of PET Hydrolases through a Ser/Ile Double Mutation.

Poly(ethylene terephthalate) (PET) is the most common polyester plastic in the packaging industry and a major source of environmental pollution due to its single use. Several enzymes, termed PET hydrolases, have been found to hydrolyze this polymer at different temperatures, with the enzyme from Ideonella sakaiensis (IsPETase) having optimal catalytic activity at 30-35 °C. Crystal structures of IsPETase have revealed that the side chain of a conserved tryptophan residue within an active site loop (W185) shifts between three conformations to enable substrate binding and product release. This is facilitated by two residues unique to IsPETase, S214 and I218. When these residues are inserted into other PET hydrolases in place of the otherwise strictly conserved histidine and phenylalanine residues found at their respective positions, they enhance activity and decrease T opt. Herein, we combine molecular dynamics and well-tempered metadynamics simulations to investigate dynamic changes of the S214/I218 and H214/F218 variants of IsPETase, as well as three other mesophilic and thermophilic PET hydrolases, at their respective temperature and pH optima. Our simulations show that the S214/I218 insertion both increases the flexibility of active site loop regions harboring key catalytic residues and the conserved tryptophan and expands the conformational plasticity of this tryptophan side chain, enabling the conformational transitions that allow for substrate binding and product release in IsPETase. The observed catalytic enhancement caused by this substitution in other PET hydrolases appears to be due to conformational selection, by capturing the conformational ensemble observed in IsPETase.

An Open One-Step RT-qPCR for SARS-CoV-2 detection.

The COVID-19 pandemic has resulted in millions of deaths globally, and while several diagnostic systems were proposed, real-time reverse transcription polymerase chain reaction (RT-PCR) remains the gold standard. However, diagnostic reagents, including enzymes used in RT-PCR, are subject to centralized production models and intellectual property restrictions, which present a challenge for less developed countries. With the aim of generating a standardized One-Step open RT-qPCR protocol to detect SARS-CoV-2 RNA in clinical samples, we purified and tested recombinant enzymes and a non-proprietary buffer. The protocol utilized M-MLV RT and Taq DNA pol enzymes to perform a Taqman probe-based assay. Synthetic RNA samples were used to validate the One-Step RT-qPCR components, and the kit showed comparable sensitivity to approved commercial kits. The One-Step RT-qPCR was then tested on clinical samples and demonstrated similar performance to commercial kits in terms of positive and negative calls. This study represents a proof of concept for an open approach to developing diagnostic kits for viral infections and diseases, which could provide a cost-effective and accessible solution for less developed countries.

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.

Discovery of New Phenylacetone Monooxygenase Variants for the Development of Substituted Indigoids through Biocatalysis.

Indigoids are natural pigments obtained from plants by ancient cultures. Romans used them mainly as dyes, whereas Asian cultures applied these compounds as treatment agents for several diseases. In the modern era, the chemical industry has made it possible to identify and develop synthetic routes to obtain them from petroleum derivatives. However, these processes require high temperatures and pressures and large amounts of solvents, acids, and alkali agents. Thus, enzyme engineering and the development of bacteria as whole-cell biocatalysts emerges as a promising green alternative to avoid the use of these hazardous materials and consequently prevent toxic waste generation. In this research, we obtained two novel variants of phenylacetone monooxygenase (PAMO) by iterative saturation mutagenesis. Heterologous expression of these two enzymes, called PAMOHPCD and PAMOHPED, in E. coli was serendipitously found to produce indigoids. These interesting results encourage us to characterize the thermal stability and enzyme kinetics of these new variants and to evaluate indigo and indirubin production in a whole-cell system by HPLC. The highest yields were obtained with PAMOHPCD supplemented with L-tryptophan, producing ~3000 mg/L indigo and ~130.0 mg/L indirubin. Additionally, both enzymes could oxidize and produce several indigo derivatives from substituted indoles, with PAMOHPCD being able to produce the well-known Tyrian purple. Our results indicate that the PAMO variants described herein have potential application in the textile, pharmaceutics, and semiconductors industries, prompting the use of environmentally friendly strategies to obtain a diverse variety of indigoids.

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.

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.