Misfolded proteins impose a dosage-dependent fitness cost and trigger a cytosolic unfolded protein response in yeast

Unraveling the mechanism of thermotolerance by Set302 in Cryptococcus neoformans

The underlying mechanism of thermotolerance, which is a key virulence factor essential for pathogenic fungi such as Cryptococcus neoformans, is largely unexplored. In this study, our findings suggest that Set302, a homolog of Set3 and a subunit of histone deacetylase complex Set3C, contributes to thermotolerance in C. neoformans. Specifically, the deletion of the predicted Set3C core subunit, Set302, resulted in further reduction in the growth of C. neoformans at 39°C, and survival of transient incubation at 50°C. Transcriptomics analysis revealed that the expression levels of numerous heat stress-responsive genes altered at both 30°C and 39°C due to the lack of Set302. Notably, at 39°C, the absence of Set302 led to the downregulation of gene expression related to the ubiquitin-proteasome system (UPS). Based on the GFP-α-synuclein overexpression model to characterize misfolded proteins, we observed a pronounced accumulation of misfolded GFP-α-synuclein at 39°C, consequently inhibiting C. neoformans thermotolerance. Furthermore, the loss of Set302 exacerbated the accumulation of misfolded GFP-α-synuclein during heat stress. Interestingly, the set302∆ strain exhibited a similar phenotype under proteasome stress as it did at 39°C. Moreover, the absence of Set302 led to reduced production of capsule and melanin. set302∆ strain also displayed significantly reduced pathogenicity and colonization ability compared to the wild-type strain in the murine infection model. Collectively, our findings suggest that Set302 modulates thermotolerance by affecting the degradation of misfolded proteins and multiple virulence factors to mediate the pathogenicity of C. neoformans.IMPORTANCECryptococcus neoformans is a pathogenic fungus that poses a potential and significant threat to public health. Thermotolerance plays a crucial role in the wide distribution in natural environments and host colonization of this fungus. Herein, Set302, a critical core subunit for the integrity of histone deacetylase complex Set3C and widely distributed in various fungi and mammals, governs thermotolerance and affects survival at extreme temperatures as well as the formation of capsule and melanin in C. neoformans. Additionally, Set302 participates in regulating the expression of multiple genes associated with the ubiquitin-proteasome system (UPS). By eliminating misfolded proteins under heat stress, Set302 significantly contributes to the thermotolerance of C. neoformans. Moreover, Set302 regulates the pathogenicity and colonization ability of C. neoformans in a murine model. Overall, this study provides new insight into the mechanism of thermotolerance in C. neoformans.

Using yeasts for the studies of nonfunctional factors in protein evolution

The evolution of protein sequence is driven not only by factors directly related to protein function and shape but also by nonfunctional factors. Such factors in protein evolution might be categorized as those connected to energetic costs, synthesis efficiency, and avoidance of misfolding and toxicity. A common approach to studying them is correlational analysis contrasting them with some characteristics of the protein, like amino acid composition, but these features are interdependent. To avoid possible bias, empirical studies are needed, and not enough work has been done to date. In this review, we describe the role of nonfunctional factors in protein evolution and present an experimental approach using yeast as a suitable model organism. The focus of the proposed approach is on the potential negative impact on the fitness of mutations that change protein properties not related to function and the frequency of mutations that change these properties. Experimental results of testing the misfolding avoidance hypothesis as an explanation for why highly expressed proteins evolve slowly are inconsistent with correlational research results. Therefore, more efforts should be made to empirically test the effects of nonfunctional factors in protein evolution and to contrast these results with the results of the correlational analysis approach.

Distinguishing mutants that resist drugs via different mechanisms by examining fitness tradeoffs

There is growing interest in designing multidrug therapies that leverage tradeoffs to combat resistance. Tradeoffs are common in evolution and occur when, for example, resistance to one drug results in sensitivity to another. Major questions remain about the extent to which tradeoffs are reliable, specifically, whether the mutants that provide resistance to a given drug all suffer similar tradeoffs. This question is difficult because the drug-resistant mutants observed in the clinic, and even those evolved in controlled laboratory settings, are often biased towards those that provide large fitness benefits. Thus, the mutations (and mechanisms) that provide drug resistance may be more diverse than current data suggests. Here, we perform evolution experiments utilizing lineage-tracking to capture a fuller spectrum of mutations that give yeast cells a fitness advantage in fluconazole, a common antifungal drug. We then quantify fitness tradeoffs for each of 774 evolved mutants across 12 environments, finding these mutants group into 6 classes with characteristically different tradeoffs. Their unique tradeoffs may imply that each group of mutants affects fitness through different underlying mechanisms. Some of the groupings we find are surprising. For example, we find some mutants that resist single drugs do not resist their combination, while others do. And some mutants to the same gene have different tradeoffs than others. These findings, on one hand, demonstrate the difficulty in relying on consistent or intuitive tradeoffs when designing multidrug treatments. On the other hand, by demonstrating that hundreds of adaptive mutations can be reduced to a few groups with characteristic tradeoffs, our findings may yet empower multidrug strategies that leverage tradeoffs to combat resistance. More generally speaking, by grouping mutants that likely affect fitness through similar underlying mechanisms, our work guides efforts to map the phenotypic effects of mutation.

Toxicity of the model protein 3×GFP arises from degradation overload, not from aggregate formation

Although protein aggregation can cause cytotoxicity, such aggregates can also form to mitigate cytotoxicity from misfolded proteins, although the nature of these contrasting aggregates remains unclear. We previously found that overproduction (op) of a three green fluorescent protein-linked protein (3×GFP) induces giant aggregates and is detrimental to growth. Here, we investigated the mechanism of growth inhibition by 3×GFP-op using non-aggregative 3×MOX-op as a control in Saccharomyces cerevisiae. The 3×GFP aggregates were induced by misfolding, and 3×GFP-op had higher cytotoxicity than 3×MOX-op because it perturbed the ubiquitin-proteasome system. Static aggregates formed by 3×GFP-op dynamically trapped Hsp70 family proteins (Ssa1 and Ssa2 in yeast), causing the heat-shock response. Systematic analysis of mutants deficient in the protein quality control suggested that 3×GFP-op did not cause a critical Hsp70 depletion and aggregation functioned in the direction of mitigating toxicity. Artificial trapping of essential cell cycle regulators into 3×GFP aggregates caused abnormalities in the cell cycle. In conclusion, the formation of the giant 3×GFP aggregates itself is not cytotoxic, as it does not entrap and deplete essential proteins. Rather, it is productive, inducing the heat-shock response while preventing an overload to the degradation system.

Dominant negative mutations in yeast Hsp90 reveal triage decision mechanism targeting client proteins for degradation

Most of the fundamental processes of cells are mediated by proteins. However, the biologically-relevant mechanism of most proteins are poorly understood. Dominant negative mutations have provided a valuable tool for investigating protein mechanisms but can be difficult to isolate because of their toxic effects. We used a mutational scanning approach to identify dominant negative mutations in yeast Hsp90. Hsp90 is a chaperone that forms dynamic complexes with many co-chaperones and client proteins. In vitro analyses have elucidated some key biochemical states and structures of Hsp90, co-chaperones, and clients; however, the biological mechanism of Hsp90 remains unclear. For example, high throughput studies have found that many E3 ubiquitin ligases bind to Hsp90, but it is unclear if these are primarily clients or acting to tag other clients for degradation. We introduced a library of all point mutations in the ATPase domain of Hsp90 into yeast and noticed that 176 were more than 10-fold depleted at the earliest point that we could analyze. There were two hot spot regions of the depleted mutations that were located at the hinges of a loop that closes over ATP. We quantified the dominant negative growth effects of mutations in the hinge regions using a library of mutations driven by an inducible promoter. We analyzed individual dominant negative mutations in detail and found that addition of the E33A mutation that prevents ATP hydrolysis by Hsp90 abrogated the dominant negative phenotype. Pull-down experiments did not reveal any stable binding partners, indicating that the dominant effects were mediated by dynamic complexes. DN Hsp90 decreased the expression level of two model Hsp90 clients, glucocorticoid receptor (GR) and v-src kinase. Using MG132, we found that GR was rapidly destabilized in a proteasome-dependent fashion. These findings provide evidence that the binding of E3 ligases to Hsp90 may serve a quality control function fundamental to eukaryotes.

An adaptive biomolecular condensation response is conserved across environmentally divergent species

Cells must sense and respond to sudden maladaptive environmental changes-stresses-to survive and thrive. Across eukaryotes, stresses such as heat shock trigger conserved responses: growth arrest, a specific transcriptional response, and biomolecular condensation of protein and mRNA into structures known as stress granules under severe stress. The composition, formation mechanism, adaptive significance, and even evolutionary conservation of these condensed structures remain enigmatic. Here we provide a remarkable view into stress-triggered condensation, its evolutionary conservation and tuning, and its integration into other well-studied aspects of the stress response. Using three morphologically near-identical budding yeast species adapted to different thermal environments and diverged by up to 100 million years, we show that proteome-scale biomolecular condensation is tuned to species-specific thermal niches, closely tracking corresponding growth and transcriptional responses. In each species, poly(A)-binding protein-a core marker of stress granules-condenses in isolation at species-specific temperatures, with conserved molecular features and conformational changes modulating condensation. From the ecological to the molecular scale, our results reveal previously unappreciated levels of evolutionary selection in the eukaryotic stress response, while establishing a rich, tractable system for further inquiry.

Agglomeration: when folded proteins clump together

Protein self-association is a widespread phenomenon that results in the formation of multimeric protein structures with critical roles in cellular processes. Protein self-association can lead to finite protein complexes or open-ended, and potentially, infinite structures. This review explores the concept of protein agglomeration, a process that results from the infinite self-assembly of folded proteins. We highlight its differences from other better-described processes with similar macroscopic features, such as aggregation and liquid-liquid phase separation. We review the sequence, structural, and biophysical factors influencing protein agglomeration. Lastly, we briefly discuss the implications of agglomeration in evolution, disease, and aging. Overall, this review highlights the need to study protein agglomeration for a better understanding of cellular processes.

Biosynthesis of natural and halogenated plant monoterpene indole alkaloids in yeast

Monoterpenoid indole alkaloids (MIAs) represent a large class of plant natural products with marketed pharmaceutical activities against a wide range of indications, including cancer, malaria and hypertension. Halogenated MIAs have shown improved pharmaceutical properties; however, synthesis of new-to-nature halogenated MIAs remains a challenge. Here we demonstrate a platform for de novo biosynthesis of two MIAs, serpentine and alstonine, in baker's yeast Saccharomyces cerevisiae and deploy it to systematically explore the biocatalytic potential of refactored MIA pathways for the production of halogenated MIAs. From this, we demonstrate conversion of individual haloindole derivatives to a total of 19 different new-to-nature haloserpentine and haloalstonine analogs. Furthermore, by process optimization and heterologous expression of a modified halogenase in the microbial MIA platform, we document de novo halogenation and biosynthesis of chloroalstonine. Together, this study highlights a microbial platform for enzymatic exploration and production of complex natural and new-to-nature MIAs with therapeutic potential.

Facile Method for High-throughput Identification of Stabilizing Mutations

Characterizing the effects of mutations on stability is critical for understanding the function and evolution of proteins and improving their biophysical properties. High throughput folding and abundance assays have been successfully used to characterize missense mutations associated with reduced stability. However, screening for increased thermodynamic stability is more challenging since such mutations are rarer and their impact on assay readout is more subtle. Here, a multiplex assay for high throughput screening of protein folding was developed by combining deep mutational scanning, fluorescence-activated cell sorting, and deep sequencing. By analyzing a library of 2000 variants of Adenylate kinase we demonstrate that the readout of the method correlates with stability and that mutants with up to 13 °C increase in thermal melting temperature could be identified with low false positive rate. The discovery of many stabilizing mutations also enabled the analysis of general substitution patterns associated with increased stability in Adenylate kinase. This high throughput method to identify stabilizing mutations can be combined with functional screens to identify mutations that improve both stability and activity.

A Zpr1 co-chaperone mediates folding of eukaryotic translation elongation factor 1A via a GTPase cycle

General protein folding is mediated by chaperones that utilize ATP hydrolysis to regulate client binding and release. Zinc-finger protein 1 (Zpr1) is an essential ATP-independent chaperone dedicated to the biogenesis of eukaryotic translation elongation factor 1A (eEF1A), a highly abundant GTP-binding protein. How Zpr1-mediated folding is regulated to ensure rapid Zpr1 recycling remains an unanswered question. Here, we use yeast genetics and microscopy analysis, biochemical reconstitution, and structural modeling to reveal that folding of eEF1A by Zpr1 requires GTP hydrolysis. Furthermore, we identify the highly conserved altered inheritance of mitochondria 29 (Aim29) protein as a Zpr1 co-chaperone that recognizes eEF1A in the GTP-bound, pre-hydrolysis conformation. This interaction dampens Zpr1⋅eEF1A GTPase activity and facilitates client exit from the folding cycle. Our work reveals that a bespoke ATP-independent chaperone system has mechanistic similarity to ATPase chaperones but unexpectedly relies on client GTP hydrolysis to regulate the chaperone-client interaction.

Intra-FCY1: a novel system to identify mutations that cause protein misfolding

Protein misfolding is a common intracellular occurrence. Most mutations to coding sequences increase the propensity of the encoded protein to misfold. These misfolded molecules can have devastating effects on cells. Despite the importance of protein misfolding in human disease and protein evolution, there are fundamental questions that remain unanswered, such as, which mutations cause the most misfolding? These questions are difficult to answer partially because we lack high-throughput methods to compare the destabilizing effects of different mutations. Commonly used systems to assess the stability of mutant proteins in vivo often rely upon essential proteins as sensors, but misfolded proteins can disrupt the function of the essential protein enough to kill the cell. This makes it difficult to identify and compare mutations that cause protein misfolding using these systems. Here, we present a novel in vivo system named Intra-FCY1 that we use to identify mutations that cause misfolding of a model protein [yellow fluorescent protein (YFP)] in Saccharomyces cerevisiae. The Intra-FCY1 system utilizes two complementary fragments of the yeast cytosine deaminase Fcy1, a toxic protein, into which YFP is inserted. When YFP folds, the Fcy1 fragments associate together to reconstitute their function, conferring toxicity in media containing 5-fluorocytosine and hindering growth. But mutations that make YFP misfold abrogate Fcy1 toxicity, thus strains possessing misfolded YFP variants rise to high frequency in growth competition experiments. This makes such strains easier to study. The Intra-FCY1 system cancels localization of the protein of interest, thus can be applied to study the relative stability of mutant versions of diverse cellular proteins. Here, we confirm this method can identify novel mutations that cause misfolding, highlighting the potential for Intra-FCY1 to illuminate the relationship between protein sequence and stability.

Ncs2* mediates in vivo virulence of pathogenic yeast through sulphur modification of cytoplasmic transfer RNA

Fungal pathogens threaten ecosystems and human health. Understanding the molecular basis of their virulence is key to develop new treatment strategies. Here, we characterize NCS2*, a point mutation identified in a clinical baker's yeast isolate. Ncs2 is essential for 2-thiolation of tRNA and the NCS2* mutation leads to increased thiolation at body temperature. NCS2* yeast exhibits enhanced fitness when grown at elevated temperatures or when exposed to oxidative stress, inhibition of nutrient signalling, and cell-wall stress. Importantly, Ncs2* alters the interaction and stability of the thiolase complex likely mediated by nucleotide binding. The absence of 2-thiolation abrogates the in vivo virulence of pathogenic baker's yeast in infected mice. Finally, hypomodification triggers changes in colony morphology and hyphae formation in the common commensal pathogen Candida albicans resulting in decreased virulence in a human cell culture model. These findings demonstrate that 2-thiolation of tRNA acts as a key mediator of fungal virulence and reveal new mechanistic insights into the function of the highly conserved tRNA-thiolase complex.

Natural variation in codon bias and mRNA folding strength interact synergistically to modify protein expression in Saccharomyces cerevisiae

Codon bias and mRNA folding strength (mF) are hypothesized molecular mechanisms by which polymorphisms in genes modify protein expression. Natural patterns of codon bias and mF across genes as well as effects of altering codon bias and mF suggest that the influence of these 2 mechanisms may vary depending on the specific location of polymorphisms within a transcript. Despite the central role codon bias and mF may play in natural trait variation within populations, systematic studies of how polymorphic codon bias and mF relate to protein expression variation are lacking. To address this need, we analyzed genomic, transcriptomic, and proteomic data for 22 Saccharomyces cerevisiae isolates, estimated protein accumulation for each allele of 1,620 genes as the log of protein molecules per RNA molecule (logPPR), and built linear mixed-effects models associating allelic variation in codon bias and mF with allelic variation in logPPR. We found that codon bias and mF interact synergistically in a positive association with logPPR, and this interaction explains almost all the effects of codon bias and mF. We examined how the locations of polymorphisms within transcripts influence their effects and found that codon bias primarily acts through polymorphisms in domain-encoding and 3' coding sequences, while mF acts most significantly through coding sequences with weaker effects from untranslated regions. Our results present the most comprehensive characterization to date of how polymorphisms in transcripts influence protein expression.

An adaptive biomolecular condensation response is conserved across environmentally divergent species

Cells must sense and respond to sudden maladaptive environmental changes-stresses-to survive and thrive. Across eukaryotes, stresses such as heat shock trigger conserved responses: growth arrest, a specific transcriptional response, and biomolecular condensation of protein and mRNA into structures known as stress granules under severe stress. The composition, formation mechanism, adaptive significance, and even evolutionary conservation of these condensed structures remain enigmatic. Here we provide an unprecedented view into stress-triggered condensation, its evolutionary conservation and tuning, and its integration into other well-studied aspects of the stress response. Using three morphologically near-identical budding yeast species adapted to different thermal environments and diverged by up to 100 million years, we show that proteome-scale biomolecular condensation is tuned to species-specific thermal niches, closely tracking corresponding growth and transcriptional responses. In each species, poly(A)-binding protein-a core marker of stress granules-condenses in isolation at species-specific temperatures, with conserved molecular features and conformational changes modulating condensation. From the ecological to the molecular scale, our results reveal previously unappreciated levels of evolutionary selection in the eukaryotic stress response, while establishing a rich, tractable system for further inquiry.

Modelling genetic stability in engineered cell populations

Predicting the evolution of engineered cell populations is a highly sought-after goal in biotechnology. While models of evolutionary dynamics are far from new, their application to synthetic systems is scarce where the vast combination of genetic parts and regulatory elements creates a unique challenge. To address this gap, we here-in present a framework that allows one to connect the DNA design of varied genetic devices with mutation spread in a growing cell population. Users can specify the functional parts of their system and the degree of mutation heterogeneity to explore, after which our model generates host-aware transition dynamics between different mutation phenotypes over time. We show how our framework can be used to generate insightful hypotheses across broad applications, from how a device's components can be tweaked to optimise long-term protein yield and genetic shelf life, to generating new design paradigms for gene regulatory networks that improve their functionality.

Extreme Sensitivity of Fitness to Environmental Conditions: Lessons from #1BigBatch

The phrase "survival of the fittest" has become an iconic descriptor of how natural selection works. And yet, precisely measuring fitness, even for single-celled microbial populations growing in controlled laboratory conditions, remains a challenge. While numerous methods exist to perform these measurements, including recently developed methods utilizing DNA barcodes, all methods are limited in their precision to differentiate strains with small fitness differences. In this study, we rule out some major sources of imprecision, but still find that fitness measurements vary substantially from replicate to replicate. Our data suggest that very subtle and difficult to avoid environmental differences between replicates create systematic variation across fitness measurements. We conclude by discussing how fitness measurements should be interpreted given their extreme environment dependence. This work was inspired by the scientific community who followed us and gave us tips as we live tweeted a high-replicate fitness measurement experiment at #1BigBatch.

Systematic conformation-to-phenotype mapping via limited deep sequencing of proteins

Non-native conformations drive protein-misfolding diseases, complicate bioengineering efforts, and fuel molecular evolution. No current experimental technique is well suited for elucidating them and their phenotypic effects. Especially intractable are the transient conformations populated by intrinsically disordered proteins. We describe an approach to systematically discover, stabilize, and purify native and non-native conformations, generated in vitro or in vivo, and directly link conformations to molecular, organismal, or evolutionary phenotypes. This approach involves high-throughput disulfide scanning (HTDS) of the entire protein. To reveal which disulfides trap which chromatographically resolvable conformers, we devised a deep-sequencing method for double-Cys variant libraries of proteins that precisely and simultaneously locates both Cys residues within each polypeptide. HTDS of the abundant E. coli periplasmic chaperone HdeA revealed distinct classes of disordered hydrophobic conformers with variable cytotoxicity depending on where the backbone was cross-linked. HTDS can bridge conformational and phenotypic landscapes for many proteins that function in disulfide-permissive environments.

Genome-wide analysis reveals Hsf1 maintains high transcript abundance of target genes controlled by strong constitutive promoter in Saccharomyces cerevisiae

BACKGROUND

In synthetic biology, the strength of promoter elements is the basis for precise regulation of target gene transcription levels, which in turn increases the yield of the target product. However, the results of many researches proved that excessive transcription levels of target genes actually reduced the yield of the target product. This phenomenon has been found in studies using different microorganisms as chassis cells, thus, it becomes a bottleneck problem to improve the yield of the target product.

RESULTS

In this study, promoters PGK1p and TDH3p with different strengths were used to regulate the transcription level of alcohol acetyl transferase encoding gene ATF1. The results demonstrated that the strong promoter TDH3p decreased the production of ethyl acetate. The results of Real-time PCR proved that the transcription level of ATF1 decreased rapidly under the control of TDH3p, and the unfolded protein reaction was activated, which may be the reason for the abnormal production caused by the strong promoter. RNA-sequencing analysis showed that the overexpression of differential gene HSP30 increased the transcriptional abundance of ATF1 gene and production of ethyl acetate. Interestingly, deletion of the heat shock protein family (e.g., Hsp26, Hsp78, Hsp82) decreased the production of ethyl acetate, suggesting that the Hsp family was also involved in the regulation of ATF1 gene transcription. Furthermore, the results proved that the Hsf1, an upstream transcription factor of Hsps, had a positive effect on alleviating the unfolded protein response and that overexpression of Hsf1 reprogramed the pattern of ATF1 gene transcript levels. The combined overexpression of Hsf1 and Hsps further increased the production of ethyl acetate. In addition, kinase Rim15 may be involved in this regulatory pathway. Finally, the regulation effect of Hsf1 on recombinant strains constructed by other promoters was verified, which confirmed the universality of the strategy.

CONCLUSIONS

Our results elucidated the mechanism by which Rim15-Hsf1-Hsps pathway reconstructed the repression of high transcription level stress and increased the production of target products, thereby providing new insights and application strategies for the construction of recombinant strains in synthetic biology.

Genes Vary Greatly in Their Propensity for Collateral Fitness Effects of Mutations

Mutations can have deleterious fitness effects when they decrease protein specific activity or decrease active protein abundance. Mutations will also be deleterious when they cause misfolding or misinteractions that are toxic to the cell (i.e., independent of whether the mutations affect specific activity and abundance). The extent to which protein evolution is shaped by these and other collateral fitness effects is unclear in part because little is known of their frequency and magnitude. Using deep mutational scanning (DMS), we previously found at least 42% of missense mutations in the TEM-1 β-lactamase antibiotic resistance gene cause deleterious collateral fitness effects. Here, we used DMS to comprehensively determine the collateral fitness effects of missense mutations in three genes encoding the antibiotic resistance proteins New Delhi metallo-β-lactamase (NDM-1), chloramphenicol acetyltransferase I (CAT-I), and 2″-aminoglycoside nucleotidyltransferase (AadB). AadB (20%), CAT-I (0.9%), and NDM-1 (0.2%) were less susceptible to deleterious collateral fitness effects than TEM-1 (42%) indicating that genes have different propensities for these effects. As was observed with TEM-1, all the studied deleterious aadB mutants increased aggregation. However, aggregation did not correlate with collateral fitness effects for many of the deleterious mutants of CAT-I and NDM-1. Select deleterious mutants caused unexpected phenotypes to emerge. The introduction of internal start codons in CAT-1 caused loss of the episome and a mutation in aadB made its cognate antibiotic essential for growth. Our study illustrates how the complexity of the cell provides a rich environment for collateral fitness effects and new phenotypes to emerge.

Systematic conformation-to-phenotype mapping via limited deep-sequencing of proteins

Non-native conformations drive protein misfolding diseases, complicate bioengineering efforts, and fuel molecular evolution. No current experimental technique is well-suited for elucidating them and their phenotypic effects. Especially intractable are the transient conformations populated by intrinsically disordered proteins. We describe an approach to systematically discover, stabilize, and purify native and non-native conformations, generated in vitro or in vivo, and directly link conformations to molecular, organismal, or evolutionary phenotypes. This approach involves high-throughput disulfide scanning (HTDS) of the entire protein. To reveal which disulfides trap which chromatographically resolvable conformers, we devised a deep-sequencing method for double-Cys variant libraries of proteins that precisely and simultaneously locates both Cys residues within each polypeptide. HTDS of the abundant E. coli periplasmic chaperone HdeA revealed distinct classes of disordered hydrophobic conformers with variable cytotoxicity depending on where the backbone was cross-linked. HTDS can bridge conformational and phenotypic landscapes for many proteins that function in disulfide-permissive environments.

Synonymous mutations in the phosphoglycerate kinase 1 gene induce an altered response to protein misfolding in Schizosaccharomyces pombe

Background

Proteostasis refers to the processes that regulate the biogenesis, folding, trafficking, and degradation of proteins. Any alteration in these processes can lead to cell malfunction. Protein synthesis, a key proteostatic process, is highly-regulated at multiple levels to ensure adequate adaptation to environmental and physiological challenges such as different stressors, proteotoxic conditions and aging, among other factors. Because alterations in protein translation can lead to protein misfolding, examining how protein translation is regulated may also help to elucidate in part how proteostasis is controlled. Codon usage bias has been implicated in the fine-tuning of translation rate, as more-frequent codons might be read faster than their less-frequent counterparts. Thus, alterations in codon usage due to synonymous mutations may alter translation kinetics and thereby affect the folding of the nascent polypeptide, without altering its primary structure. To date, it has been difficult to predict the effect of synonymous mutations on protein folding and cellular fitness due to a scarcity of relevant data. Thus, the purpose of this work was to assess the effect of synonymous mutations in discrete regions of the gene that encodes the highly-expressed enzyme 3-phosphoglycerate kinase 1 (pgk1) in the fission yeast Schizosaccharomyces pombe.

Results

By means of systematic replacement of synonymous codons along pgk1, we found slightly-altered protein folding and activity in a region-specific manner. However, alterations in protein aggregation, heat stress as well as changes in proteasome activity occurred independently of the mutated region. Concomitantly, reduced mRNA levels of the chaperones Hsp9 and Hsp16 were observed.

Conclusion

Taken together, these data suggest that codon usage bias of the gene encoding this highly-expressed protein is an important regulator of protein function and proteostasis.

Purifying selection enduringly acts on the sequence evolution of highly expressed proteins in Escherichia coli

The evolutionary speed of a protein sequence is constrained by its expression level, with highly expressed proteins evolving relatively slowly. This negative correlation between expression levels and evolutionary rates (known as the E–R anticorrelation) has already been widely observed in past macroevolution between species from bacteria to animals. However, it remains unclear whether this seemingly general law also governs recent evolution, including past and de novo, within a species. However, the advent of genomic sequencing and high-throughput phenotyping, particularly for bacteria, has revealed fundamental gaps between the 2 evolutionary processes and has provided empirical data opposing the possible underlying mechanisms which are widely believed. These conflicts raise questions about the generalization of the E–R anticorrelation and the relevance of plausible mechanisms. To explore the ubiquitous impact of expression levels on molecular evolution and test the relevance of the possible underlying mechanisms, we analyzed the genome sequences of 99 strains of Escherichia coli for evolution within species in nature. We also analyzed genomic mutations accumulated under laboratory conditions as a model of de novo evolution within species. Here, we show that E–R anticorrelation is significant in both past and de novo evolution within species in E. coli. Our data also confirmed ongoing purifying selection on highly expressed genes. Ongoing selection included codon-level purifying selection, supporting the relevance of the underlying mechanisms. However, the impact of codon-level purifying selection on the constraints in evolution within species might be smaller than previously expected from evolution between species.

Nondiploid cancer cells: Stress, tolerance and therapeutic inspirations

Aberrant ploidy status is a prominent characteristic in malignant neoplasms. Approximately 90% of solid tumors and 75% of haematopoietic malignancies contain aneuploidy cells, and 30%-60% of tumors undergo whole-genome doubling, indicating that nondiploidy might be a prevalent genomic aberration in cancer. Although the role of aneuploid and polyploid cells in cancer remains to be elucidated, recent studies have suggested that nondiploid cells might be a dangerous minority that severely challenges cancer management. Ploidy shifts cause multiple fitness coasts for cancer cells, mainly including genomic, proteotoxic, metabolic and immune stresses. However, nondiploid comprises a well-adopted subpopulation, with many tolerance mechanisms evident in cells along with ploidy shifts. Aneuploid and polyploid cells elegantly maintain an autonomous balance between the stress and tolerance during adaptive evolution in cancer. Breaking the balance might provide some inspiration for ploidy-selective cancer therapy and alleviation of ploidy-related chemoresistance. To understand of the complex role and therapeutic potential of nondiploid cells better, we reviewed the survival stresses and adaptive tolerances within nondiploid cancer cells and summarized therapeutic ploidy-selective alterations for potential use in developing future cancer therapy.

Genetics of adaptation

The rediscovery of Mendel’s work showing that the heredity of phenotypes is controlled by discrete genes was followed by the reconciliation of Mendelian genetics with evolution by natural selection in the middle of the last century with the Modern Synthesis. In the past two decades, dramatic advances in genomic methods have facilitated the identification of the loci, genes, and even individual mutations that underlie phenotypic variants that are the putative targets of natural selection. Moreover, these methods have also changed how we can study adaptation by flipping the problem around, allowing us to first examine what loci show evidence of having been under selection, and then connecting these genetic variants to phenotypic variation. As a result, we now have an expanding list of actual genetic changes that underlie potentially adaptive phenotypic variation. Here, we synthesize how considering the effects of these adaptive loci in the context of cellular environments, genomes, organisms, and populations has provided new insights to the genetic architecture of adaptation.

A limit on the evolutionary rescue of an Antarctic bacterium from rising temperatures

Climate change is gradual, but it can also cause brief extreme heat waves that can exceed the upper thermal limit of any one organism. To study the evolutionary potential of upper thermal tolerance, we evolved the cold-adapted Antarctic bacterium Pseudoalteromonas haloplanktis to survive at 30°C, beyond its ancestral thermal limit. This high-temperature adaptation occurred rapidly and in multiple populations. It involved genomic changes that occurred in a highly parallel fashion and mitigated the effects of protein misfolding. However, it also confronted a physiological limit, because populations failed to grow beyond 30°C. Our experiments aimed to facilitate evolutionary rescue by using a small organism with large populations living at temperatures several degrees below their upper thermal limit. Larger organisms with smaller populations and living at temperatures closer to their upper thermal tolerances are even more likely to go extinct during extreme heat waves.

DnaK response to expression of protein mutants is dependent on translation rate and stability

Chaperones play a central part in the quality control system in cells by clearing misfolded and aggregated proteins. The chaperone DnaK acts as a sensor for molecular stress by recognising short hydrophobic stretches of misfolded proteins. As the level of unfolded protein is a function of protein stability, we hypothesised that the level of DnaK response upon overexpression of recombinant proteins would be correlated to stability. Using a set of mutants of the λ-repressor with varying thermal stabilities and a fluorescent reporter system, the effect of stability on DnaK response and protein abundance was investigated. Our results demonstrate that the initial DnaK response is largely dependent on protein synthesis rate but as the recombinantly expressed protein accumulates and homeostasis is approached the response correlates strongly with stability. Furthermore, we observe a large degree of cell-cell variation in protein abundance and DnaK response in more stable proteins.

Massive expression of cysteine-containing proteins causes abnormal elongation of yeast cells by perturbing the proteasome

The enhanced green fluorescent protein (EGFP) is considered to be a harmless protein because the critical expression level that causes growth defects is higher than that of other proteins. Here, we found that overexpression of EGFP, but not a glycolytic protein Gpm1, triggered the cell elongation phenotype in the budding yeast Saccharomyces cerevisiae. By the morphological analysis of the cell overexpressing fluorescent protein and glycolytic enzyme variants, we revealed that cysteine content was associated with the cell elongation phenotype. The abnormal cell morphology triggered by overexpression of EGFP was also observed in the fission yeast Schizosaccharomyces pombe. Overexpression of cysteine-containing protein was toxic, especially at high-temperature, while the toxicity could be modulated by additional protein characteristics. Investigation of protein aggregate formation, morphological abnormalities in mutants, and transcriptomic changes that occur upon overexpression of EGFP variants suggested that perturbation of the proteasome by the exposed cysteine of the overexpressed protein causes cell elongation. Overexpression of proteins with relatively low folding properties, such as EGFP, was also found to promote the formation of SHOTA (Seventy kDa Heat shock protein-containing, Overexpression-Triggered Aggregates), an intracellular aggregate that incorporates Hsp70/Ssa1, which induces a heat shock response, while it was unrelated to cell elongation. Evolutionary analysis of duplicated genes showed that cysteine toxicity may be an evolutionary bias to exclude cysteine from highly expressed proteins. The overexpression of cysteine-less moxGFP, the least toxic protein revealed in this study, would be a good model system to understand the physiological state of protein burden triggered by ultimate overexpression of harmless proteins.

Quantifying shifts in natural selection on codon usage between protein regions: a population genetics approach

Background

Codon usage bias (CUB), the non-uniform usage of synonymous codons, occurs across all domains of life. Adaptive CUB is hypothesized to result from various selective pressures, including selection for efficient ribosome elongation, accurate translation, mRNA secondary structure, and/or protein folding. Given the critical link between protein folding and protein function, numerous studies have analyzed the relationship between codon usage and protein structure. The results from these studies have often been contradictory, likely reflecting the differing methods used for measuring codon usage and the failure to appropriately control for confounding factors, such as differences in amino acid usage between protein structures and changes in the frequency of different structures with gene expression.

Results

Here we take an explicit population genetics approach to quantify codon-specific shifts in natural selection related to protein structure in S. cerevisiae and E. coli. Unlike other metrics of codon usage, our approach explicitly separates the effects of natural selection, scaled by gene expression, and mutation bias while naturally accounting for a region’s amino acid usage. Bayesian model comparisons suggest selection on codon usage varies only slightly between helix, sheet, and coil secondary structures and, similarly, between structured and intrinsically-disordered regions. Similarly, in contrast to prevous findings, we find selection on codon usage only varies slightly at the termini of helices in E. coli. Using simulated data, we show this previous work indicating “non-optimal” codons are enriched at the beginning of helices in S. cerevisiae was due to failure to control for various confounding factors (e.g. amino acid biases, gene expression, etc.), and rather than selection to modulate cotranslational folding.

Conclusions

Our results reveal a weak relationship between codon usage and protein structure, indicating that differences in selection on codon usage between structures are slight. In addition to the magnitude of differences in selection between protein structures being slight, the observed shifts appear to be idiosyncratic and largely codon-specific rather than systematic reversals in the nature of selection. Overall, our work demonstrates the statistical power and benefits of studying selective shifts on codon usage or other genomic features from an explicitly evolutionary approach. Limitations of this approach and future potential research avenues are discussed.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08635-0).

Chaperones directly and efficiently disperse stress-triggered biomolecular condensates

Stresses such as heat shock trigger the formation of protein aggregates and the induction of a disaggregation system composed of molecular chaperones. Recent work reveals that several cases of apparent heat-induced aggregation, long thought to be the result of toxic misfolding, instead reflect evolved, adaptive biomolecular condensation, with chaperone activity contributing to condensate regulation. Here we show that the yeast disaggregation system directly disperses heat-induced biomolecular condensates of endogenous poly(A)-binding protein (Pab1) orders of magnitude more rapidly than aggregates of the most commonly used misfolded model substrate, firefly luciferase. Beyond its efficiency, heat-induced condensate dispersal differs from heat-induced aggregate dispersal in its molecular requirements and mechanistic behavior. Our work establishes a bona fide endogenous heat-induced substrate for long-studied heat shock proteins, isolates a specific example of chaperone regulation of condensates, and underscores needed expansion of the proteotoxic interpretation of the heat shock response to encompass adaptive, chaperone-mediated regulation.

The Unfolded Protein Response at the Tumor-Immune Interface

The tumor-immune interface has surged to primary relevance in an effort to understand the hurdles facing immune surveillance and cancer immunotherapy. Reports over the past decades have indicated a role for the unfolded protein response (UPR) in modulating not only tumor cell fitness and drug resistance, but also local immunity, with emphasis on the phenotype and altered function of immune cells such as myeloid cells and T cells. Emerging evidence also suggests that aneuploidy correlates with local immune dysregulation. Recently, we reported that the UPR serves as a link between aneuploidy and immune cell dysregulation in a cell nonautonomous way. These new findings add considerable complexity to the organization of the tumor microenvironment (TME) and the origin of its altered function. In this review, we summarize these data and also discuss the role of aneuploidy as a negative regulator of local immunity.

Proteome evaluation of homolog abundance patterns in Arachis hypogaea cv. Tifrunner

BACKGROUND

Cultivated peanut (Arachis hypogaea, AABB genome), an allotetraploid from a cross between A. duranensis (AA genome) and A. ipaensis (BB genome), is an important oil and protein crop with released genome and RNA-seq sequence datasets. These datasets provide the molecular foundation for studying gene expression and evolutionary patterns. However, there are no reports on the proteomic data of A. hypogaea cv. Tifrunner, which limits understanding of its gene function and protein level evolution.

RESULTS

This study sequenced the A. hypogaea cv. Tifrunner leaf and root proteome using the tandem mass tag technology. A total of 4803 abundant proteins were identified. The 364 differentially abundant proteins were estimated by comparing protein abundances between leaf and root proteomes. The differentially abundant proteins enriched the photosystem process. The number of biased abundant homeologs between the two sub-genomes A (87 homeologs in leaf and root) and B (69 and 68 homeologs in leaf and root, respectively) was not significantly different. However, homeologous proteins with biased abundances in different sub-genomes enriched different biological processes. In the leaf, homeologs biased to sub-genome A enriched biosynthetic and metabolic process, while homeologs biased to sub-genome B enriched iron ion homeostasis process. In the root, homeologs with biased abundance in sub-genome A enriched inorganic biosynthesis and metabolism process, while homeologs with biased abundance in sub-genome B enriched organic biosynthesis and metabolism process. Purifying selection mainly acted on paralogs and homeologs. The selective pressure values were negatively correlated with paralogous protein abundance. About 77.42% (24/31) homeologous and 80% (48/60) paralogous protein pairs had asymmetric abundance, and several protein pairs had conserved abundances in the leaf and root tissues.

CONCLUSIONS

This study sequenced the proteome of A. hypogaea cv. Tifrunner using the leaf and root tissues. Differentially abundant proteins were identified, and revealed functions. Paralog abundance divergence and homeolog bias abundance was elucidated. These results indicate that divergent abundance caused retention of homologs in A. hypogaea cv. Tifrunner.

Systematic investigation of the link between enzyme catalysis and cold adaptation

Cold temperature is prevalent across the biosphere and slows the rates of chemical reactions. Increased catalysis has been predicted to be a dominant adaptive trait of enzymes to reduced temperature, and this expectation has informed physical models for enzyme catalysis and influenced bioprospecting strategies. To systematically test rate enhancement as an adaptive trait to cold, we paired kinetic constants of 2223 enzyme reactions with their organism's optimal growth temperature (T_Growth) and analyzed trends of rate constants as a function of T_Growth. These data do not support a general increase in rate enhancement in cold adaptation. In the model enzyme ketosteroid isomerase (KSI), there is prior evidence for temperature adaptation from a change in an active site residue that results in a tradeoff between activity and stability. Nevertheless, we found that little of the rate constant variation for 20 KSI variants was accounted for by T_Growth. In contrast, and consistent with prior expectations, we observed a correlation between stability and T_Growth across 433 proteins. These results suggest that temperature exerts a weaker selection pressure on enzyme rate constants than stability and that evolutionary forces other than temperature are responsible for the majority of enzymatic rate constant variation.

The heat shock response and small molecule regulators

The heat shock response (HSR) is a highly conserved cellular pathway that is responsible for stress relief and the refolding of denatured proteins [1]. When a host cell is exposed to conditions such as heat shock, ischemia, or toxic substances, heat shock factor-1 (HSF-1), a transcription factor, activates the genes that encode for the heat shock proteins (Hsps), which are a family of proteins that work alongside other chaperones to relieve stress and refold proteins that have been denatured (Burdon, 1986) [2]. Along with the refolding of denatured proteins, Hsps facilitate the removal of misfolded proteins by escorting them to degradation pathways, thereby preventing the accumulation of misfolded proteins [3]. Research has indicated that many pathological conditions, such as diabetes, cancer, neuropathy, cardiovascular disease, and aging have a negative impact on HSR function and are commonly associated with misfolded protein aggregation [4,5]. Studies indicate an interplay between mitochondrial homeostasis and HSF-1 levels can impact stress resistance, proteostasis, and malignant cell growth, which further support the role of Hsps in pathological and metabolic functions [6]. On the other hand, Hsp activation by specific small molecules can induce the heat shock response, which can afford neuroprotection and other benefits [7]. This review will focus on the modulation of Hsps and the HSR as therapeutic options to treat these conditions.

Analyses of HIV proteases variants at the threshold of viability reveals relationships between processing efficiency and fitness

Investigating the relationships between protein function and fitness provides keys for understanding biochemical mechanisms that underly evolution. Mutations with partial fitness defects can delineate the threshold of biochemical function required for viability. We utilized a previous deep mutational scan of HIV-1 protease (PR) to identify variants with 15–45 per cent defects in replication and analysed the biochemical function of eight variants (L10M, L10S, V32C, V32I, A71V, A71S, Q92I, Q92N). We purified each variant and assessed the efficiency of peptide cleavage for three cut sites (MA-CA, TF-PR, and PR-RT) as well as gel-based analyses of processing of purified Gag. The cutting activity of at least one site was perturbed relative to WT protease for all variants, consistent with cutting activity being a primary determinant of fitness effects. We examined the correlation of fitness defects with cutting activity of different sites. MA-CA showed the weakest correlation (R² = 0.02) with fitness, suggesting relatively weak coupling with viral replication. In contrast, cutting of the TF-PR site showed the strongest correlation with fitness (R² = 0.53). Cutting at the TF-PR site creates a new PR protein with a free N-terminus that is critical for activity. Our findings indicate that increasing the pool of active PR is rate limiting for viral replication, making this an ideal step to target with inhibitors.

The unfolded protein response links tumor aneuploidy to local immune dysregulation

Aneuploidy is a chromosomal abnormality associated with poor prognosis in many cancer types. Here, we tested the hypothesis that the unfolded protein response (UPR) mechanistically links aneuploidy and local immune dysregulation. Using a single somatic copy number alteration (SCNA) score inclusive of whole‐chromosome, chromosome arm, and focal alterations in a pan‐cancer analysis of 9,375 samples in The Cancer Genome Atlas (TCGA) database, we found an inverse correlation with a cytotoxicity (CYT) score across disease stages. Co‐expression patterns of UPR genes changed substantially between SCNA^low and SCNA^high groups. Pathway activity scores showed increased activity of multiple branches of the UPR in response to aneuploidy. The PERK branch showed the strongest association with a reduced CYT score. The conditioned medium of aneuploid cells transmitted XBP1 splicing and caused IL‐6 and arginase 1 transcription in receiver bone marrow‐derived macrophages and markedly diminished the production of IFN‐γ and granzyme B in activated human T cells. We propose the UPR as a mechanistic link between aneuploidy and immune dysregulation in the tumor microenvironment.

Cell adaptation to aneuploidy by the environmental stress response dampens induction of the cytosolic unfolded-protein response

Aneuploid yeast cells are in a chronic state of proteotoxicity, yet do not constitutively induce the cytosolic unfolded protein response, or heat shock response (HSR) by heat shock factor 1 (Hsf1). Here, we demonstrate that an active environmental stress response (ESR), a hallmark of aneuploidy across different models, suppresses Hsf1 induction in models of single-chromosome gain. Furthermore, engineered activation of the ESR in the absence of stress was sufficient to suppress Hsf1 activation in euploid cells by subsequent heat shock while increasing thermotolerance and blocking formation of heat-induced protein aggregates. Suppression of the ESR in aneuploid cells resulted in longer cell doubling times and decreased viability in the presence of additional proteotoxicity. Last, we show that in euploids, Hsf1 induction by heat shock is curbed by the ESR. Strikingly, we found a similar relationship between the ESR and the HSR using an inducible model of aneuploidy. Our work explains a long-standing paradox in the field and provides new insights into conserved mechanisms of proteostasis with potential relevance to cancers associated with aneuploidy.

An Overexpression Experiment Does Not Support the Hypothesis That Avoidance of Toxicity Determines the Rate of Protein Evolution

The misfolding avoidance hypothesis postulates that sequence mutations render proteins cytotoxic and therefore the higher the gene expression, the stronger the operation of selection against substitutions. This translates into prediction that relative toxicity of extant proteins is higher for those evolving faster. In the present experiment, we selected pairs of yeast genes which were paralogous but evolving at different rates. We expressed them artificially to high levels. We expected that toxicity would be higher for ones bearing more mutations, especially that overcrowding should rather exacerbate than reverse the already existing differences in misfolding rates. We did find that the applied mode of overexpression caused a considerable decrease in fitness and that the decrease was proportional to the amount of excessive protein. However, it was not higher for proteins which are normally expressed at lower levels (and have less conserved sequence). This result was obtained consistently, regardless whether the rate of growth or ability to compete in common cultures was used as a proxy for fitness. In additional experiments, we applied factors that reduce accuracy of translation or enhance structural instability of proteins. It did not change a consistent pattern of independence between the fitness cost caused by overexpression of a protein and the rate of its sequence evolution.

Hsf1 activation by proteotoxic stress requires concurrent protein synthesis

Heat shock factor 1 (Hsf1) activation is responsible for increasing the abundance of protein-folding chaperones and degradation machinery in response to proteotoxic conditions that give rise to misfolded or aggregated proteins. Here we systematically explored the link between concurrent protein synthesis and proteotoxic stress in the budding yeast, Saccharomyces cerevisiae. Consistent with prior work, inhibiting protein synthesis before inducing proteotoxic stress prevents Hsf1 activation, which we demonstrated across a broad array of stresses and validate using orthogonal means of blocking protein synthesis. However, other stress-dependent transcription pathways remained activatable under conditions of translation inhibition. Titrating the protein denaturant ethanol to a higher concentration results in Hsf1 activation in the absence of translation, suggesting extreme protein-folding stress can induce proteotoxicity independent of protein synthesis. Furthermore, we demonstrate this connection under physiological conditions where protein synthesis occurs naturally at reduced rates. We find that disrupting the assembly or subcellular localization of newly synthesized proteins is sufficient to activate Hsf1. Thus, new proteins appear to be especially sensitive to proteotoxic conditions, and we propose that their aggregation may represent the bulk of the signal that activates Hsf1 in the wake of these insults.

Avoidance of protein unfolding constrains protein stability in long-term evolution

Every amino acid residue can influence a protein's overall stability, making stability highly susceptible to change throughout evolution. We consider the distribution of protein stabilities evolutionarily permittable under two previously reported protein fitness functions: flux dynamics and misfolding avoidance. We develop an evolutionary dynamics theory and find that it agrees better with an extensive protein stability data set for dihydrofolate reductase orthologs under the misfolding avoidance fitness function rather than the flux dynamics fitness function. Further investigation with ribonuclease H data demonstrates that not any misfolded state is avoided; rather, it is only the unfolded state. At the end, we discuss how our work pertains to the universal protein abundance-evolutionary rate correlation seen across organisms' proteomes. We derive a closed-form expression relating protein abundance to evolutionary rate that captures Escherichia coli, Saccharomyces cerevisiae, and Homo sapiens experimental trends without fitted parameters.

Human RGR Gene and Associated Features of Age-Related Macular Degeneration in Models of Retina-Choriocapillaris Atrophy

Age-related macular degeneration (AMD) is a progressive eye disease and the most common cause of blindness among the elderly. AMD is characterized by early atrophy of the choriocapillaris and retinal pigment epithelium (RPE). Although AMD is a multifactorial disease with many environmental and genetic risk factors, a hallmark of the disease is the origination of extracellular deposits, or drusen, between the RPE and Bruch membrane. Human retinal G-protein-coupled receptor (RGR) gene generates an exon-skipping splice variant of RGR-opsin (RGR-d; NP_001012740) that is a persistent component of small and large drusen. Herein, the findings show that abnormal RGR proteins, including RGR-d, are pathogenic in an animal retina with degeneration of the choriocapillaris, RPE, and photoreceptors. A frameshift truncating mutation resulted in severe retinal degeneration with a continuous band of basal deposits along the Bruch membrane. RGR-d produced less severe disease with choriocapillaris and RPE atrophy, including focal accumulation of abnormal RGR-d protein at the basal boundary of the RPE. Degeneration of the choriocapillaris was marked by a decrease in endothelial CD31 protein and choriocapillaris breakdown at the ultrastructural level. Fundus lesions with patchy depigmentation were characteristic of old RGR-d mice. RGR-d was mislocalized in cultured cells and caused a strong cell growth defect. These results uphold the notion of a potential hidden link between AMD and a high-frequency RGR allele.

Endoplasmic reticulum stress regulates the intestinal stem cell state through CtBP2

Enforcing differentiation of cancer stem cells is considered as a potential strategy to sensitize colorectal cancer cells to irradiation and chemotherapy. Activation of the unfolded protein response, due to endoplasmic reticulum (ER) stress, causes rapid stem cell differentiation in normal intestinal and colon cancer cells. We previously found that stem cell differentiation was mediated by a Protein kinase R-like ER kinase (PERK) dependent arrest of mRNA translation, resulting in rapid protein depletion of WNT-dependent transcription factor c-MYC. We hypothesize that ER stress dependent stem cell differentiation may rely on the depletion of additional transcriptional regulators with a short protein half-life that are rapidly depleted due to a PERK-dependent translational pause. Using a novel screening method, we identify novel transcription factors that regulate the intestinal stem cell fate upon ER stress. ER stress was induced in LS174T cells with thapsigargin or subtilase cytotoxin (SubAB) and immediate alterations in nuclear transcription factor activity were assessed by the CatTFRE assay in which transcription factors present in nuclear lysate are bound to plasmid DNA, co-extracted and quantified using mass-spectrometry. The role of altered activity of transcription factor CtBP2 was further examined by modification of its expression levels using CAG-rtTA3-CtBP2 overexpression in small intestinal organoids, shCtBP2 knockdown in LS174T cells, and familial adenomatous polyposis patient-derived organoids. CtBP2 overexpression organoids were challenged by ER stress and ionizing irradiation. We identified a unique set of transcription factors with altered activation upon ER stress. Gene ontology analysis showed that transcription factors with diminished binding were involved in cellular differentiation processes. ER stress decreased CtBP2 protein expression in mouse small intestine. ER stress induced loss of CtBP2 expression which was rescued by inhibition of PERK signaling. CtBP2 was overexpressed in mouse and human colorectal adenomas. Inducible CtBP2 overexpression in organoids conferred higher clonogenic potential, resilience to irradiation-induced damage and a partial rescue of ER stress-induced loss of stemness. Using an unbiased proteomics approach, we identified a unique set of transcription factors for which DNA-binding activity is lost directly upon ER stress. We continued investigating the function of co-regulator CtBP2, and show that CtBP2 mediates ER stress-induced loss of stemness which supports the intestinal stem cell state in homeostatic stem cells and colorectal cancer cells.

Resource uptake and the evolution of moderately efficient enzymes

Enzymes speed up reactions that would otherwise be too slow to sustain the metabolism of selfreplicators. Yet, most enzymes seem only moderately efficient, exhibiting kinetic parameters orders of magnitude lower than their expected physically achievable maxima and spanning over surprisingly large ranges of values. Here, we question how these parameters evolve using a mechanistic model where enzyme efficiency is a key component of individual competition for resources. We show that kinetic parameters are under strong directional selection only up to a point, above which enzymes appear to evolve under near-neutrality, thereby confirming the qualitative observation of other modeling approaches. While the existence of a large fitness plateau could potentially explain the extensive variation in enzyme features reported, we show using a population genetics model that such a widespread distribution is an unlikely outcome of evolution on a common landscape, as mutation–selection–drift balance occupy a narrow area even when very moderate biases towards lower efficiency are considered. Instead, differences in the evolutionary context encountered by each enzyme should be involved, such that each evolves on an individual, unique landscape. Our results point to drift and effective population size playing an important role, along with the kinetics of nutrient transporters, the tolerance to high concentrations of intermediate metabolites, and the reversibility of reactions. Enzyme concentration also shapes selection on kinetic parameters, but we show that the joint evolution of concentration and efficiency does not yield extensive variance in evolutionary outcomes when documented costs to protein expression are applied.

Investigating the role of ERAD on antibody processing in glycoengineered Saccharomyces cerevisiae

N-glycosylation plays an important role in the endoplasmic reticulum quality control (ERQC). N-glycan biosynthesis pathways have been engineered in yeasts and fungi to enable the production of therapeutic glycoproteins with human-compatible N-glycosylation, and some glycoengineering approaches alter the synthesis of the lipid-linked oligosaccharide (LLO). Because the effects of LLO engineering on ERQC are currently unknown, we characterized intracellular processing of IgG in glycoengineered Δalg3 Δalg11 Saccharomyces cerevisiae strain and analyzed how altered LLO structures affect endoplasmic reticulum-associated degradation (ERAD). Intracellular IgG light and heavy chain molecules expressed in Δalg3 Δalg11 strain are ERAD substrates and targeted to ERAD independently of Yos9p and Htm1p, whereas in the presence of ALG3 ERAD targeting is dependent on Yos9p but does not require Htm1p. Blocking of ERAD accumulated ER and post-Golgi forms of IgG and increased glycosylation of matα secretion signal but did not improve IgG secretion. Our results show ERAD targeting of a heterologous glycoprotein in yeast, and suggest that proteins in the ER can be targeted to ERAD via other mechanisms than the Htm1p-Yos9p-dependent route when the LLO biosynthesis is altered.

Consequences of mitotic failure - The penalties and the rewards

Eukaryotic cells are usually diploid, meaning they contain two copies of each chromosome. However, aberrant chromosome numbers due to both, chromosome gains and losses, are often observed in nature. They can occur as a planned developmental step, but are more often an uninvited result of mitotic failure. Recent discoveries have improved our understanding of the cellular effects of aneuploidy - uneven chromosome numbers, and polyploidy - multiplication of entire sets of chromosomes - in eukaryotic cells. The results show that mitotic errors lead to rapid and extensive modifications of many cellular processes and affect proliferation, proteome balance, genome stability and more. The findings picture the cellular response to aneuploidy and polyploidy as a complex, tissue and context dependent network of events. Here I review the latest discoveries, with an emphasis on pathological aspects of aneuploidy and polyploidy in human cells.

A Reporter System for Cytosolic Protein Aggregates in Yeast

Protein misfolding and aggregation are linked to neurodegenerative diseases of mammals and suboptimal protein expression within biotechnology. Tools for monitoring protein aggregates are therefore useful for studying disease-related aggregation and for improving soluble protein expression in heterologous hosts for biotechnology purposes. In this work, we developed a promoter-reporter system for aggregated protein on the basis of the yeast native response to misfolded protein. To this end, we first studied the proteome of yeast in response to the expression of folded soluble and aggregation-prone protein baits and identified genes encoding proteins related to protein folding and the response to heat stress as well as the ubiquitin-proteasome system that are over-represented in cells expressing an aggregation-prone protein. From these data, we created and validated promoter-reporter constructs and further engineered the best performing promoters by increasing the copy number of upstream activating sequences and optimization of culture conditions. Our best promoter-reporter has an output dynamic range of approximately 12-fold upon expression of the aggregation-prone protein and responded to increasing levels of aggregated protein. Finally, we demonstrate that the system can discriminate between yeast cells expressing different prion precursor proteins and select the cells expressing folded soluble protein from mixed populations. Our reporter system is thus a simple tool for diagnosing protein aggregates in living cells and should be applicable for the health and biotechnology industries.

Adaptability of the ubiquitin-proteasome system to proteolytic and folding stressors

Aging, disease, and environmental stressors are associated with failures in the ubiquitin-proteasome system (UPS), yet a quantitative understanding of how stressors affect the proteome and how the UPS responds is lacking. Here we assessed UPS performance and adaptability in yeast under stressors using quantitative measurements of misfolded substrate stability and stress-dependent UPS regulation by the transcription factor Rpn4. We found that impairing degradation rates (proteolytic stress) and generating misfolded proteins (folding stress) elicited distinct effects on the proteome and on UPS adaptation. Folding stressors stabilized proteins via aggregation rather than overburdening the proteasome, as occurred under proteolytic stress. Still, the UPS productively adapted to both stressors using separate mechanisms: proteolytic stressors caused Rpn4 stabilization while folding stressors increased RPN4 transcription. In some cases, adaptation completely prevented loss of UPS substrate degradation. Our work reveals the distinct effects of proteotoxic stressors and the versatility of cells in adapting the UPS.

Structure and function of naturally evolved de novo proteins

Comparative evolutionary genomics has revealed that novel protein coding genes can emerge randomly from non-coding DNA. While most of the myriad of transcripts which continuously emerge vanish rapidly, some attain regulatory regions, become translated and survive. More surprisingly, sequence properties of de novo proteins are almost indistinguishable from randomly obtained sequences, yet de novo proteins may gain functions and integrate into eukaryotic cellular networks quite easily. We here discuss current knowledge on de novo proteins, their structures, functions and evolution. Since the existence of de novo proteins seems at odds with decade-long attempts to construct proteins with novel structures and functions from scratch, we suggest that a better understanding of de novo protein evolution may fuel new strategies for protein design.

The Relationship between the Misfolding Avoidance Hypothesis and Protein Evolutionary Rates in the Light of Empirical Evidence

For more than a decade, the misfolding avoidance hypothesis (MAH) and related theories have dominated evolutionary discussions aimed at explaining the variance of the molecular clock across cellular proteins. In this study, we use various experimental data to further investigate the consistency of the MAH predictions with empirical evidence. We also critically discuss experimental results that motivated the MAH development and that are often viewed as evidence of its major contribution to the variability of protein evolutionary rates. We demonstrate, in Escherichia coli and Homo sapiens, the lack of a substantial negative correlation between protein evolutionary rates and Gibbs free energies of unfolding, a direct measure of protein stability. We then analyze multiple new genome-scale data sets characterizing protein aggregation and interaction propensities, the properties that are likely optimized in evolution to alleviate deleterious effects associated with toxic protein misfolding and misinteractions. Our results demonstrate that the propensity of proteins to aggregate, the fraction of charged amino acids, and protein stickiness do correlate with protein abundances. Nevertheless, across multiple organisms and various data sets we do not observe substantial correlations between proteins’ aggregation- and stability-related properties and evolutionary rates. Therefore, diverse empirical data support the conclusion that the MAH and similar hypotheses do not play a major role in mediating a strong negative correlation between protein expression and the molecular clock, and thus in explaining the variability of evolutionary rates across cellular proteins.

Validation of DBFOLD: An efficient algorithm for computing folding pathways of complex proteins

Atomistic simulations can provide valuable, experimentally-verifiable insights into protein folding mechanisms, but existing ab initio simulation methods are restricted to only the smallest proteins due to severe computational speed limits. The folding of larger proteins has been studied using native-centric potential functions, but such models omit the potentially crucial role of non-native interactions. Here, we present an algorithm, entitled DBFOLD, which can predict folding pathways for a wide range of proteins while accounting for the effects of non-native contacts. In addition, DBFOLD can predict the relative rates of different transitions within a protein’s folding pathway. To accomplish this, rather than directly simulating folding, our method combines equilibrium Monte-Carlo simulations, which deploy enhanced sampling, with unfolding simulations at high temperatures. We show that under certain conditions, trajectories from these two types of simulations can be jointly analyzed to compute unknown folding rates from detailed balance. This requires inferring free energies from the equilibrium simulations, and extrapolating transition rates from the unfolding simulations to lower, physiologically-reasonable temperatures at which the native state is marginally stable. As a proof of principle, we show that our method can accurately predict folding pathways and Monte-Carlo rates for the well-characterized Streptococcal protein G. We then show that our method significantly reduces the amount of computation time required to compute the folding pathways of large, misfolding-prone proteins that lie beyond the reach of existing direct simulation. Our algorithm, which is available online, can generate detailed atomistic models of protein folding mechanisms while shedding light on the role of non-native intermediates which may crucially affect organismal fitness and are frequently implicated in disease.

Many proteins must adopt a specific structure in order to function. Computational simulations have been used to shed light on the mechanisms of protein folding, but unfortunately, realistic simulations can typically only be run for small proteins, due to severe limits in computational speed. Here, we present a method to solve this problem, whereby instead of directly simulating folding from an unfolded state, we run simulations that allow for computation of equilibrium folding free energies, alongside high temperature simulations to compute unfolding rates. From these quantities, folding rates can be computed using detailed balance. Importantly, our method can account for the effects of nonnative contacts which transiently form during folding and must be broken prior to adoption of the native state. Such contacts, which are often excluded from simple models of folding, may crucially affect real protein folding pathways and are often observed in folding intermediates implicated in disease.

Yeast Cellular Stress: Impacts on Bioethanol Production

Bioethanol is the largest biotechnology product and the most dominant biofuel globally. Saccharomyces cerevisiae is the most favored microorganism employed for its industrial production. However, obtaining maximum yields from an ethanol fermentation remains a technical challenge, since cellular stresses detrimentally impact on the efficiency of yeast cell growth and metabolism. Ethanol fermentation stresses potentially include osmotic, chaotropic, oxidative, and heat stress, as well as shifts in pH. Well-developed stress responses and tolerance mechanisms make S. cerevisiae industrious, with bioprocessing techniques also being deployed at industrial scale for the optimization of fermentation parameters and the effective management of inhibition issues. Overlap exists between yeast responses to different forms of stress. This review outlines yeast fermentation stresses and known mechanisms conferring stress tolerance, with their further elucidation and improvement possessing the potential to improve fermentation efficiency.