Protein evolution speed depends on its stability and abundance and on chaperone concentrations

Interactions between chaperone and energy storage networks during the evolution of Legionella pneumophila under heat shock

Waterborne transmission of the bacterium Legionella pneumophila has emerged as a major cause of severe nosocomial infections of major public health impact. The major route of transmission involves the uptake of aerosolized bacteria, often from the contaminated hot water systems of large buildings. Public health regulations aimed at controlling the mesophilic pathogen are generally concerned with acute pasteurization and maintaining high temperatures at the heating systems and throughout the plumbing of hot water systems, but L. pneumophila is often able to survive these treatments due to both bacterium-intrinsic and environmental factors. Previous work has established an experimental evolution system to model the observations of increased heat resistance in repeatedly but unsuccessfully pasteurized L. pneumophila populations. Here, we show rapid fixation of novel alleles in lineages selected for resistance to heat shock and shifts in mutational profile related to increases in the temperature of selection. Gene-level and nucleotide-level parallelisms between independently-evolving lineages show the centrality of the DnaJ/DnaK chaperone system in the heat resistance of L. pneumophila. Inference of epistatic interactions through reverse genetics shows an unexpected interaction between DnaJ/DnaK and the polyhydroxybutyrate-accumulation energy storage mechanism used by the species to survive long-term starvation in low-nutrient environments.

High-throughput Selection of Human de novo-emerged sORFs with High Folding Potential

De novo genes emerge from previously noncoding stretches of the genome. Their encoded de novo proteins are generally expected to be similar to random sequences and, accordingly, with no stable tertiary fold and high predicted disorder. However, structural properties of de novo proteins and whether they differ during the stages of emergence and fixation have not been studied in depth and rely heavily on predictions. Here we generated a library of short human putative de novo proteins of varying lengths and ages and sorted the candidates according to their structural compactness and disorder propensity. Using Förster resonance energy transfer combined with Fluorescence-activated cell sorting, we were able to screen the library for most compact protein structures, as well as most elongated and flexible structures. We find that compact de novo proteins are on average slightly shorter and contain lower predicted disorder than less compact ones. The predicted structures for most and least compact de novo proteins correspond to expectations in that they contain more secondary structure content or higher disorder content, respectively. Our experiments indicate that older de novo proteins have higher compactness and structural propensity compared with young ones. We discuss possible evolutionary scenarios and their implications underlying the age-dependencies of compactness and structural content of putative de novo proteins.

High temperature delays and low temperature accelerates evolution of a new protein phenotype

Since the origin of life, temperatures on earth have fluctuated both on short and long time scales. How such changes affect the rate at which Darwinian evolution can bring forth new phenotypes remains unclear. On the one hand, high temperature may accelerate phenotypic evolution because it accelerates most biological processes. On the other hand, it may slow phenotypic evolution, because proteins are usually less stable at high temperatures and therefore less evolvable. Here, to test these hypotheses experimentally, we evolved a green fluorescent protein in E. coli towards the new phenotype of yellow fluorescence at different temperatures. Yellow fluorescence evolved most slowly at high temperature and most rapidly at low temperature, in contradiction to the first hypothesis. Using high-throughput population sequencing, protein engineering, and biochemical assays, we determined that this is due to the protein-destabilizing effect of neofunctionalizing mutations. Destabilization is highly detrimental at high temperature, where neofunctionalizing mutations cannot be tolerated. Their detrimental effects can be mitigated through excess stability at low temperature, leading to accelerated adaptive evolution. By modifying protein folding stability, temperature alters the accessibility of mutational paths towards high-fitness genotypes. Our observations have broad implications for our understanding of how temperature changes affect evolutionary adaptations and innovations.

Engineered odorant receptors illuminate structural principles of odor discrimination

A central challenge in olfaction is understanding how the olfactory system detects and distinguishes odorants with diverse physicochemical properties and molecular configurations. Vertebrate animals perceive odors via G protein-coupled odorant receptors (ORs). In humans, ~400 ORs enable the sense of smell. The OR family is composed of two major classes: Class I ORs are tuned to carboxylic acids while Class II ORs, representing the vast majority of the human repertoire, respond to a wide variety of odorants. How ORs recognize chemically diverse odorants remains poorly understood. A fundamental bottleneck is the inability to visualize odorant binding to ORs. Here, we uncover fundamental molecular properties of odorant-OR interactions by employing engineered ORs crafted using a consensus protein design strategy. Because such consensus ORs (consORs) are derived from the 17 major subfamilies of human ORs, they provide a template for modeling individual native ORs with high sequence and structural homology. The biochemical tractability of consORs enabled four cryoEM structures of distinct consORs with unique ligand recognition properties. The structure of a Class I consOR, consOR51, showed high structural similarity to the native human receptor OR51E2 and yielded a homology model of a related member of the human OR51 family with high predictive power. Structures of three Class II consORs revealed distinct modes of odorant-binding and activation mechanisms between Class I and Class II ORs. Thus, the structures of consORs lay the groundwork for understanding molecular recognition of odorants by the OR superfamily.

Mega-scale experimental analysis of protein folding stability in biology and design

Advances in DNA sequencing and machine learning are providing insights into protein sequences and structures on an enormous scale1. However, the energetics driving folding are invisible in these structures and remain largely unknown2. The hidden thermodynamics of folding can drive disease3,4, shape protein evolution5-7 and guide protein engineering8-10, and new approaches are needed to reveal these thermodynamics for every sequence and structure. Here we present cDNA display proteolysis, a method for measuring thermodynamic folding stability for up to 900,000 protein domains in a one-week experiment. From 1.8 million measurements in total, we curated a set of around 776,000 high-quality folding stabilities covering all single amino acid variants and selected double mutants of 331 natural and 148 de novo designed protein domains 40-72 amino acids in length. Using this extensive dataset, we quantified (1) environmental factors influencing amino acid fitness, (2) thermodynamic couplings (including unexpected interactions) between protein sites, and (3) the global divergence between evolutionary amino acid usage and protein folding stability. We also examined how our approach could identify stability determinants in designed proteins and evaluate design methods. The cDNA display proteolysis method is fast, accurate and uniquely scalable, and promises to reveal the quantitative rules for how amino acid sequences encode folding stability.

Underlying Role of Hydrophobic Environments in Tuning Metal Elements for Efficient Enzyme Catalysis

The catalytic functions of metalloenzymes are often strongly correlated with metal elements in the active sites. However, dioxygen-activating nonheme quercetin dioxygenases (QueD) are found with various first-row transition-metal ions when metal swapping inactivates their innate catalytic activity. To unveil the molecular basis of this seemingly promiscuous yet metal-specific enzyme, we transformed manganese-dependent QueD into a nickel-dependent enzyme by sequence- and structure-based directed evolution. Although the net effect of acquired mutations was primarily to rearrange hydrophobic residues in the active site pocket, biochemical, kinetic, X-ray crystallographic, spectroscopic, and computational studies suggest that these modifications in the secondary coordination spheres can adjust the electronic structure of the enzyme-substrate complex to counteract the effects induced by the metal substitution. These results explicitly demonstrate that such noncovalent interactions encrypt metal specificity in a finely modulated manner, revealing the underestimated chemical power of the hydrophobic sequence network in enzyme catalysis.

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.

Global Analysis of Multi-Mutants to Improve Protein Function

The identification of amino acid substitutions that both enhance the stability and function of a protein is a key challenge in protein engineering. Technological advances have enabled assaying thousands of protein variants in a single high-throughput experiment, and more recent studies use such data in protein engineering. We present a Global Multi-Mutant Analysis (GMMA) that exploits the presence of multiply-substituted variants to identify individual amino acid substitutions that are beneficial for the stability and function across a large library of protein variants. We have applied GMMA to a previously published experiment reporting on >54,000 variants of green fluorescent protein (GFP), each with known fluorescence output, and each carrying 1-15 amino acid substitutions (Sarkisyan et al., 2016). The GMMA method achieves a good fit to this dataset while being analytically transparent. We show experimentally that the six top-ranking substitutions progressively enhance GFP. More broadly, using only a single experiment as input our analysis recovers nearly all the substitutions previously reported to be beneficial for GFP folding and function. In conclusion, we suggest that large libraries of multiply-substituted variants may provide a unique source of information for protein engineering.

Phylogenetic and Structural Analyses of VPS13 Proteins in Archaeplastida Reveal Their Complex Evolutionary History in Viridiplantae

VPS13 is a lipid transfer protein family conserved among Eukaryotes and playing roles in fundamental processes involving vesicular transport and membrane expansion including autophagy and organelle biogenesis. VPS13 folds into a long hydrophobic tunnel, allowing lipid transport, decorated by distinct domains involved in protein localization and regulation. Whereas VPS13 organization and function have been extensively studied in yeast and mammals, information in organisms originating from primary endosymbiosis is scarce. In the higher plant Arabidopsis thaliana, four paralogs, AtVPS13S, X, M1, and M2, were identified, AtVPS13S playing a role in the regulation of root growth, cell patterning, and reproduction. In this work, we performed phylogenetic, as well as domain and structural modeling of VPS13 proteins in Archaeplastida in order to understand their general organization and evolutionary history. We confirmed the presence of human VPS13B orthologues in some phyla and described two new VPS13 families presenting a particular domain arrangement: VPS13R in Rhodophytes and VPS13Y in Chlorophytes and Streptophytes. By focusing on Viridiplantae, we were able to draw the evolutionary history of these proteins made by multiple gene gains and duplications as well as domain rearrangements. We showed that some Chlorophytes have only three (AtVPS13M, S, Y) whereas some Charophytes have up to six VPS13 paralogs (AtVPS13M1, M2, S, Y, X, B). We also highlighted specific structural features of VPS13M and X paralogs. This study reveals the complex evolution of VPS13 family and opens important perspectives for their functional characterization in photosynthetic organisms.

Protein Stability: Enhancement and Measurement

This chapter defines protein stability, emphasizes its importance, and surveys the field of protein stabilization, with summary reference to a selection of 2014-2021 publications. One can enhance stability, particularly by protein engineering strategies but also by chemical modification and by other means. General protocols are set out on how to measure a given protein's (i) kinetic thermal stability and (ii) oxidative stability and (iii) how to undertake chemical modification of a protein in solution.

Bacterial Hsp90 predominantly buffers but does not potentiate the phenotypic effects of deleterious mutations during fluorescent protein evolution

Chaperones facilitate the folding of other ("client") proteins and can thus affect the adaptive evolution of these clients. Specifically, chaperones affect the phenotype of proteins via two opposing mechanisms. On the one hand, they can buffer the effects of mutations in proteins and thus help preserve an ancestral, premutation phenotype. On the other hand, they can potentiate the effects of mutations and thus enhance the phenotypic changes caused by a mutation. We study that how the bacterial Hsp90 chaperone (HtpG) affects the evolution of green fluorescent protein. To this end, we performed directed evolution of green fluorescent protein under low and high cellular concentrations of Hsp90. Specifically, we evolved green fluorescent protein under both stabilizing selection for its ancestral (green) phenotype and directional selection toward a new (cyan) phenotype. While Hsp90 did only affect the rate of adaptive evolution transiently, it did affect the phenotypic effects of mutations that occurred during adaptive evolution. Specifically, Hsp90 allowed strongly deleterious mutations to accumulate in evolving populations by buffering their effects. Our observations show that the role of a chaperone for adaptive evolution depends on the organism and the trait being studied.

Lessons Learned from Two Decades of Modeling the Heat-Shock Response

The Heat Shock Response (HSR) is a highly conserved genetic system charged with protecting the proteome in a wide range of organisms and species. Experiments since the early 1980s have elucidated key elements in these pathways and revealed a canonical mode of regulation, which relies on a titration feedback. This system has been subject to substantial modeling work, addressing questions about resilience, design and control. The compact core regulatory circuit, as well as its apparent conservation, make this system an ideal ‘hydrogen atom’ model for the regulation of stress response. Here we take a broad view of the models of the HSR, focusing on the different questions asked and the approaches taken. After 20 years of modeling work, we ask what lessons had been learned that would have been hard to discover without mathematical models. We find that while existing models lay strong foundations, many important questions that can benefit from quantitative modeling are still awaiting investigation.

The Statistical Trends of Protein Evolution: A Lesson from AlphaFold Database

The recent development of artificial intelligence provides us with new and powerful tools for studying the mysterious relationship between organism evolution and protein evolution. In this work, based on the AlphaFold Protein Structure Database (AlphaFold DB), we perform comparative analyses of the proteins of different organisms. The statistics of AlphaFold-predicted structures show that, for organisms with higher complexity, their constituent proteins will have larger radii of gyration, higher coil fractions, and slower vibrations, statistically. By conducting normal mode analysis and scaling analyses, we demonstrate that higher organismal complexity correlates with lower fractal dimensions in both the structure and dynamics of the constituent proteins, suggesting that higher functional specialization is associated with higher organismal complexity. We also uncover the topology and sequence bases of these correlations. As the organismal complexity increases, the residue contact networks of the constituent proteins will be more assortative, and these proteins will have a higher degree of hydrophilic-hydrophobic segregation in the sequences. Furthermore, by comparing the statistical structural proximity across the proteomes with the phylogenetic tree of homologous proteins, we show that, statistical structural proximity across the proteomes may indirectly reflect the phylogenetic proximity, indicating a statistical trend of protein evolution in parallel with organism evolution. This study provides new insights into how the diversity in the functionality of proteins increases and how the dimensionality of the manifold of protein dynamics reduces during evolution, contributing to the understanding of the origin and evolution of lives.

Pulse labeling reveals the tail end of protein folding by proteome profiling

Accurate and efficient folding of nascent protein sequences into their native states requires support from the protein homeostasis network. Herein we probe which newly translated proteins are thermo-sensitive, making them susceptible to misfolding and aggregation under heat stress using pulse-SILAC mass spectrometry. We find a distinct group of proteins that is highly sensitive to this perturbation when newly synthesized but not once matured. These proteins are abundant and highly structured. Notably, they display a tendency to form β sheet secondary structures, have more complex folding topology, and are enriched for chaperone-binding motifs, suggesting a higher demand for chaperone-assisted folding. These polypeptides are also more often components of stable protein complexes in comparison with other proteins. Combining these findings suggests the existence of a specific subset of proteins in the cell that is particularly vulnerable to misfolding and aggregation following synthesis before reaching the native state.

Mechanisms of protein evolution

How do proteins evolve? How do changes in sequence mediate changes in protein structure, and in turn in function? This question has multiple angles, ranging from biochemistry and biophysics to evolutionary biology. This review provides a brief integrated view of some key mechanistic aspects of protein evolution. First, we explain how protein evolution is primarily driven by randomly acquired genetic mutations and selection for function, and how these mutations can even give rise to completely new folds. Then, we also comment on how phenotypic protein variability, including promiscuity, transcriptional and translational errors, may also accelerate this process, possibly via "plasticity-first" mechanisms. Finally, we highlight open questions in the field of protein evolution, with respect to the emergence of more sophisticated protein systems such as protein complexes, pathways, and the emergence of pre-LUCA enzymes.

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.

Protein lifetimes in aged brains reveal a proteostatic adaptation linking physiological aging to neurodegeneration

Aging is a prominent risk factor for neurodegenerative disorders (NDDs); however, the molecular mechanisms rendering the aged brain particularly susceptible to neurodegeneration remain unclear. Here, we aim to determine the link between physiological aging and NDDs by exploring protein turnover using metabolic labeling and quantitative pulse-SILAC proteomics. By comparing protein lifetimes between physiologically aged and young adult mice, we found that in aged brains protein lifetimes are increased by ~20% and that aging affects distinct pathways linked to NDDs. Specifically, a set of neuroprotective proteins are longer-lived in aged brains, while some mitochondrial proteins linked to neurodegeneration are shorter-lived. Strikingly, we observed a previously unknown alteration in proteostasis that correlates to parsimonious turnover of proteins with high biosynthetic costs, revealing an overall metabolic adaptation that preludes neurodegeneration. Our findings suggest that future therapeutic paradigms, aimed at addressing these metabolic adaptations, might be able to delay NDD onset.

In vitro evolution of uracil glycosylase towards DnaKJ and GroEL binding evolves different misfolded states

Natural evolution is driven by random mutations that improve fitness. In vitro evolution mimics this process, however, on a short time-scale and is driven by the given bait. Here, we used directed in vitro evolution of a random mutant library of Uracil glycosylase (eUNG) displayed on yeast surface to select for binding to chaperones GroEL, DnaK+DnaJ+ATP (DnaKJ) or E.coli cell extract (CE), using binding to the eUNG inhibitor Ugi as probe for native fold. The CE selected population was further divided to Ugi binders (+U) or non-binders (-U). The aim here was to evaluate the sequence space and physical state of the evolved protein binding the different baits. We found that GroEL, DnaKJ and CE-U select and enrich for mutations causing eUNG to misfold, with the three being enriched in mutations in buried and conserved positions, with a tendency to increase positive charge. Still, each selection had its own trajectory, with GroEL and CE-U selecting mutants highly sensitive to protease cleavage while DnaKJ selected partially structured misfolded species with a tendency to refold, making them less sensitive to proteases. More general, our results show that GroEL has a higher tendency to purge promiscuous misfolded protein mutants from the system, while DnaKJ binds misfolding-prone mutant species that are, upon chaperone release, more likely to natively refold. CE-U shares some of the properties of GroEL- and DnaKJ-selected populations, while harboring also unique properties that can be explained by the presence of additional chaperones in CE, such as Trigger factor, HtpG and ClpB.

GroEL/S overexpression helps to purge deleterious mutations and reduce genetic diversity during adaptive protein evolution

Chaperones are proteins that help other proteins fold. They also affect the adaptive evolution of their client proteins by buffering the effect of deleterious mutations and increasing the genetic diversity of evolving proteins. We study how the bacterial chaperone GroE (GroEL + GroES) affects the evolution of green fluorescent protein (GFP). To this end we subjected GFP to multiple rounds of mutation and selection for its color phenotype in four replicate E. coli populations, and studied its evolutionary dynamics through high-throughput sequencing and mutant engineering. We evolved GFP both under stabilizing selection for its ancestral (green) phenotype, and to directional selection for a new (cyan) phenotype. We did so both under low and high expression of the chaperone GroE. In contrast to previous work, we observe that GroE does not just buffer but also helps purge deleterious (fluorescence reducing) mutations from evolving populations. In doing so, GroE helps reduce the genetic diversity of evolving populations. In addition, it causes phenotypic heterogeneity in mutants with the same genotype, helping to enhance their fluorescence in some cells, and reducing it in others. Our observations show that chaperones can affect adaptive evolution in more than one way.

In vivo, in vitro and in silico: an open space for the development of microbe-based applications of synthetic biology

Living systems are studied using three complementary approaches: living cells, cell-free systems and computer-mediated modelling. Progresses in understanding, allowing researchers to create novel chassis and industrial processes rest on a cycle that combines in vivo, in vitro and in silico studies. This design-build-test-learn iteration loop cycle between experiments and analyses combines together physiology, genetics, biochemistry and bioinformatics in a way that keeps going forward. Because computer-aided approaches are not directly constrained by the material nature of the entities of interest, we illustrate here how this virtuous cycle allows researchers to explore chemistry which is foreign to that present in extant life, from whole chassis to novel metabolic cycles. Particular emphasis is placed on the importance of evolution.

Histone H3K27 methylation perturbs transcriptional robustness and underpins dispensability of highly conserved genes in fungi

Epigenetic modifications are key regulators of gene expression and underpin genome integrity. Yet, how epigenetic changes affect the evolution and transcriptional robustness of genes remains largely unknown. Here, we show how the repressive histone mark H3K27me3 underpins the trajectory of highly conserved genes in fungi. We first performed transcriptomic profiling on closely related species of the plant pathogen Fusarium graminearum species complex. We determined transcriptional responsiveness of genes across environmental conditions to determine expression robustness. To infer evolutionary conservation, we used a framework of 23 species across the Fusarium genus including three species covered with histone methylation data. Gene expression variation is negatively correlated with gene conservation confirming that highly conserved genes show higher expression robustness. In contrast, genes marked by H3K27me3 do not show such associations. Furthermore, highly conserved genes marked by H3K27me3 encode smaller proteins, exhibit weaker codon usage bias, higher levels of hydrophobicity, show lower intrinsically disordered regions, and are enriched for functions related to regulation and membrane transport. The evolutionary age of conserved genes with H3K27me3 histone marks falls typically within the origins of the Fusarium genus. We show that highly conserved genes marked by H3K27me3 are more likely to be dispensable for survival during host infection. Lastly, we show that conserved genes exposed to repressive H3K27me3 marks across distantly related Fusarium fungi are associated with transcriptional perturbation at the microevolutionary scale. In conclusion, we show how repressive histone marks are entangled in the evolutionary fate of highly conserved genes across evolutionary timescales.

The interactions of molecular chaperones with client proteins: why are they so weak?

The major classes of molecular chaperones have highly variable sequences, sizes, and shapes, yet they all bind to unfolded proteins, limit their aggregation, and assist in their folding. Despite the central importance of this process to protein homeostasis, it has not been clear exactly how chaperones guide this process or whether the diverse families of chaperones use similar mechanisms. For the first time, recent advances in NMR spectroscopy have enabled detailed studies of how unfolded, "client" proteins interact with both ATP-dependent and ATP-independent classes of chaperones. Here, we review examples from four distinct chaperones, Spy, Trigger Factor, DnaK, and HscA-HscB, highlighting the similarities and differences between their mechanisms. One striking similarity is that the chaperones all bind weakly to their clients, such that the chaperone-client interactions are readily outcompeted by stronger, intra- and intermolecular contacts in the folded state. Thus, the relatively weak affinity of these interactions seems to provide directionality to the folding process. However, there are also key differences, especially in the details of how the chaperones release clients and how ATP cycling impacts that process. For example, Spy releases clients in a largely folded state, while clients seem to be unfolded upon release from Trigger Factor or DnaK. Together, these studies are beginning to uncover the similarities and differences in how chaperones use weak interactions to guide protein folding.

Heterogeneity in conservation of multifunctional partner enzymes with meiotic importance, CDK2 kinase and BRCA1 ubiquitin ligase

The evolution of proteins can be accompanied by changes not only to their amino acid sequences, but also their structural and spatial molecular organization. Comparison of the protein conservation within different taxonomic groups (multifunctional, or highly specific) allows to clarify their specificity and the direction of evolution. Two multifunctional enzymes, cyclin-dependent kinase 2 (CDK2) and BRCA1 ubiquitin ligase, that are partners in some mitotic and meiotic processes were investigated in the present work. Two research methods, bioinformatics and immunocytochemical, were combined to examine the conservation levels of the two enzymes. It has been established that CDK2 is a highly conserved protein in different taxonomic lineages of the eukaryotic tree. Immunocytochemically, a conserved CDK2 pattern was revealed in the meiotic autosomes of five rodent species and partially in domestic turkey and clawed frog. Nevertheless, variable CDK2 distribution was detected at the unsynapsed segments of the rodent X chromosomes. BRCA1 was shown to be highly conserved only within certain mammalian taxa. It was also noted that in those rodent nuclei, where BRCA1 specifically binds to antigens, asynaptic regions of sex chromosomes were positive. BRCA1 staining was not always accompanied by specific binding, and a high nonspecificity in the nucleoplasm was observed. Thus, the studies revealed different conservation of the two enzymes at the level of protein structure as well as at the level of chromosome behavior. This suggests variable rates of evolution due to both size and configuration of the protein molecules and their multifunctionality.

The Protein Folding Problem: The Role of Theory

The protein folding problem was first articulated as question of how order arose from disorder in proteins: How did the various native structures of proteins arise from interatomic driving forces encoded within their amino acid sequences, and how did they fold so fast? These matters have now been largely resolved by theory and statistical mechanics combined with experiments. There are general principles. Chain randomness is overcome by solvation-based codes. And in the needle-in-a-haystack metaphor, native states are found efficiently because protein haystacks (conformational ensembles) are funnel-shaped. Order-disorder theory has now grown to encompass a large swath of protein physical science across biology.

Spatially Enriched Paralog Rearrangements Argue Functionally Diverse Ribosomes Arise during Cold Acclimation in Arabidopsis

Ribosome biogenesis is essential for plants to successfully acclimate to low temperature. Without dedicated steps supervising the 60S large subunits (LSUs) maturation in the cytosol, e.g., Rei-like (REIL) factors, plants fail to accumulate dry weight and fail to grow at suboptimal low temperatures. Around REIL, the final 60S cytosolic maturation steps include proofreading and assembly of functional ribosomal centers such as the polypeptide exit tunnel and the P-Stalk, respectively. In consequence, these ribosomal substructures and their assembly, especially during low temperatures, might be changed and provoke the need for dedicated quality controls. To test this, we blocked ribosome maturation during cold acclimation using two independent reil double mutant genotypes and tested changes in their ribosomal proteomes. Additionally, we normalized our mutant datasets using as a blank the cold responsiveness of a wild-type Arabidopsis genotype. This allowed us to neglect any reil-specific effects that may happen due to the presence or absence of the factor during LSU cytosolic maturation, thus allowing us to test for cold-induced changes that happen in the early nucleolar biogenesis. As a result, we report that cold acclimation triggers a reprogramming in the structural ribosomal proteome. The reprogramming alters the abundance of specific RP families and/or paralogs in non-translational LSU and translational polysome fractions, a phenomenon known as substoichiometry. Next, we tested whether the cold-substoichiometry was spatially confined to specific regions of the complex. In terms of RP proteoforms, we report that remodeling of ribosomes after a cold stimulus is significantly constrained to the polypeptide exit tunnel (PET), i.e., REIL factor binding and functional site. In terms of RP transcripts, cold acclimation induces changes in RP families or paralogs that are significantly constrained to the P-Stalk and the ribosomal head. The three modulated substructures represent possible targets of mechanisms that may constrain translation by controlled ribosome heterogeneity. We propose that non-random ribosome heterogeneity controlled by specialized biogenesis mechanisms may contribute to a preferential or ultimately even rigorous selection of transcripts needed for rapid proteome shifts and successful acclimation.

Consequences of mutation accumulation for growth performance are more likely to be resource-dependent at higher temperatures

Background

Mutation accumulation (MA) has profound ecological and evolutionary consequences. One example is that accumulation of conditionally neutral mutations leads to fitness trade-offs among heterogenous habitats which cause population divergence. Here we suggest that temperature, which controls the rates of all biochemical and biophysical processes, should play a crucial role for determining mutational effects. Particularly, warmer temperatures may mitigate the effects of some, not all, deleterious mutations and cause stronger environmental dependence in MA effects.

Results

We experimentally tested the above hypothesis by measuring the growth performance of ten Escherichia coli genotypes on six carbon resources across ten temperatures, where the ten genotypes were derived from a single ancestral strain and accumulated spontaneous mutations. We analyzed resource dependence of MA consequences for growth yields. The MA genotypes typically showed reduced growth yields relative to the ancestral type; and the magnitude of reduction was smaller at intermediate temperatures. Stronger resource dependence in MA consequences for growth performance was observed at higher temperatures. Specifically, the MA genotypes were more likely to show impaired growth performance on all the six carbon resources when grown at lower temperatures; but suffered growth performance loss only on some, not all the six, carbon substrates at higher temperatures.

Conclusions

Higher temperatures increase the chance that MA causes conditionally neutral fitness effects while MA is more likely to cause fitness loss regardless of available resources at lower temperatures. This finding has implications for understanding how geographic patterns in population divergence may emerge, and how conservation practices, particularly protection of diverse microhabitats, may mitigate the impacts of global warming.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01846-1.

Self-assembled peptide and protein nanostructures for anti-cancer therapy: Targeted delivery, stimuli-responsive devices and immunotherapy

Self-assembled peptides and proteins possess tremendous potential as targeted drug delivery systems and key applications of these well-defined nanostructures reside in anti-cancer therapy. Peptides and proteins can self-assemble into nanostructures of diverse sizes and shapes in response to changing environmental conditions such as pH, temperature, ionic strength, as well as host and guest molecular interactions; their countless benefits include good biocompatibility and high loading capacity for hydrophobic and hydrophilic drugs. These self-assembled nanomaterials can be adorned with functional moieties to specifically target tumor cells. Stimuli-responsive features can also be incorporated with respect to the tumor microenvironment. This review sheds light on the growing interest in self-assembled peptides and proteins and their burgeoning applications in cancer treatment and immunotherapy.

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.

Engineering promiscuity of chloramphenicol acetyltransferase for microbial designer ester biosynthesis

Robust and efficient enzymes are essential modules for metabolic engineering and synthetic biology strategies across biological systems to engineer whole-cell biocatalysts. By condensing an acyl-CoA and an alcohol, alcohol acyltransferases (AATs) can serve as interchangeable metabolic modules for microbial biosynthesis of a diverse class of ester molecules with broad applications as flavors, fragrances, solvents, and drop-in biofuels. However, the current lack of robust and efficient AATs significantly limits their compatibility with heterologous precursor pathways and microbial hosts. Through bioprospecting and rational protein engineering, we identified and engineered promiscuity of chloramphenicol acetyltransferases (CATs) from mesophilic prokaryotes to function as robust and efficient AATs compatible with at least 21 alcohol and 8 acyl-CoA substrates for microbial biosynthesis of linear, branched, saturated, unsaturated and/or aromatic esters. By plugging the best engineered CAT (CATec3 Y20F) into the gram-negative mesophilic bacterium Escherichia coli, we demonstrated that the recombinant strain could effectively convert various alcohols into desirable esters, for instance, achieving a titer of 13.9 g/L isoamyl acetate with 95% conversion by fed-batch fermentation. The recombinant E. coli was also capable of simulating the ester profile of roses with high conversion (>97%) and titer (>1 g/L) from fermentable sugars at 37 °C. Likewise, a recombinant gram-positive, cellulolytic, thermophilic bacterium Clostridium thermocellum harboring CATec3 Y20F could produce many of these esters from recalcitrant cellulosic biomass at elevated temperatures (>50 °C) due to the engineered enzyme's remarkable thermostability. Overall, the engineered CATs can serve as a robust and efficient platform for designer ester biosynthesis from renewable and sustainable feedstocks.

Fitness effects of spontaneous mutations in a warming world

Spontaneous mutations fuel evolutionary processes and differ in consequence, but the consequences depend on the environment. Biophysical considerations of protein thermostability predict that warm temperatures may systematically increase the deleteriousness of mutation. We sought to test whether mutation reduced fitness more when measured in an environment that reflected climate change projections for temperature. We investigated the effects of spontaneous mutations on life history, size, and fitness in 21 mutation accumulation lines and 12 control lines of Daphnia pulex at standard and elevated (+4℃) temperatures. Warmer temperature accelerated life history and reduced body length and clutch sizes. Mutation led to reduced mean clutch sizes and fitness estimates at both temperatures. We found no evidence of a systematic temperature–mutation interaction on trait means, although some lines showed evidence of beneficial mutation at one temperature and deleterious mutation at the other. However, trait variances are also influenced by mutation, and we observed increased variances due to mutation for most traits. For variance of the intrinsic rate of increase and some reproductive traits, we found significant temperature–mutation interactions, with a larger increase due to mutation in the warmer environment. This suggests that selection on new mutations will be more efficient at elevated temperatures.

Dynamical spectroscopy and microscopy of proteins in cells

With a strong understanding of how proteins fold in hand, it is now possible to ask how in-cell environments modulate their folding, binding and function. Studies accessing fast (ns to s) in-cell dynamics have accelerated over the past few years through a combination of in-cell NMR spectroscopy and time-resolved fluorescence microscopies. Here, we discuss this recent work and the emerging picture of protein surfaces as not just hydrophilic coats interfacing the solvent to the protein's core and functional regions, but as critical components in cells controlling protein mobility, function and communication with post-translational modifications.

Structural-functional diversity of malaria parasite's PfHSP70-1 and PfHSP40 chaperone pair gives an edge over human orthologs in chaperone-assisted protein folding

Plasmodium falciparum, the human malaria parasite harbors a metastable proteome which is vulnerable to proteotoxic stress conditions encountered during its lifecycle. How parasite's chaperone machinery is able to maintain its aggregation-prone proteome in functional state, is poorly understood. As HSP70-40 system forms the central hub in cellular proteostasis, we investigated the protein folding capacity of PfHSP70-1 and PfHSP40 chaperone pair and compared it with human orthologs (HSPA1A and DNAJA1). Despite the structural similarity, we observed that parasite chaperones and their human orthologs exhibit striking differences in conformational dynamics. Comprehensive biochemical investigations revealed that PfHSP70-1 and PfHSP40 chaperone pair has better protein folding, aggregation inhibition, oligomer remodeling and disaggregase activities than their human orthologs. Chaperone-swapping experiments suggest that PfHSP40 can also efficiently cooperate with human HSP70 to facilitate the folding of client-substrate. SPR-derived kinetic parameters reveal that PfHSP40 has higher binding affinity towards unfolded substrate than DNAJA1. Interestingly, the observed slow dissociation rate of PfHSP40-substrate interaction allows PfHSP40 to maintain the substrate in folding-competent state to minimize its misfolding. Structural investigation through small angle x-ray scattering gave insights into the conformational architecture of PfHSP70-1 (monomer), PfHSP40 (dimer) and their complex. Overall, our data suggest that the parasite has evolved functionally diverged and efficient chaperone machinery which allows the human malaria parasite to survive in hostile conditions. The distinct allosteric landscapes and interaction kinetics of plasmodial chaperones open avenues for the exploration of small-molecule based antimalarial interventions.

Elevated temperature increases genome-wide selection on de novo mutations

Adaptation in new environments depends on the amount of genetic variation available for evolution, and the efficacy by which natural selection discriminates among this variation. However, whether some ecological factors reveal more genetic variation, or impose stronger selection pressures than others, is typically not known. Here, we apply the enzyme kinetic theory to show that rising global temperatures are predicted to intensify natural selection throughout the genome by increasing the effects of DNA sequence variation on protein stability. We test this prediction by (i) estimating temperature-dependent fitness effects of induced mutations in seed beetles adapted to ancestral or elevated temperature, and (ii) calculate 100 paired selection estimates on mutations in benign versus stressful environments from unicellular and multicellular organisms. Environmental stress per se did not increase mean selection on de novo mutation, suggesting that the cost of adaptation does not generally increase in new ecological settings to which the organism is maladapted. However, elevated temperature increased the mean strength of selection on genome-wide polymorphism, signified by increases in both mutation load and mutational variance in fitness. These results have important implications for genetic diversity gradients and the rate and repeatability of evolution under climate change.

Incorporating antagonistic pleiotropy into models for molecular replicators

In modern cells, chromosomal genes composed of DNA encode multi-subunit protein/RNA complexes that catalyze the replication of the chromosome and cell. One prevailing theory for the origin of life posits an early stage involving self-replicating macromolecules called replicators, which can be considered genes capable of self-replication. One prevailing theory for the genetics of aging in humans and other organisms is antagonistic pleiotropy, which posits that a gene can be beneficial in one context, and detrimental in another context. We previously reported that the conceptual simplicity of molecular replicators facilitates the generation of two simple models involving antagonistic pleiotropy. Here a third model is proposed, and each of the three models is presented with improved definition of the time variable. Computer simulations were used to calculate the proliferation of a hypothetical two-subunit replicator (AB), when one of the two subunits (B) exhibits antagonistic pleiotropy, leading to an advantage for B to be unstable. In model 1, instability of B yields free A subunits, which in turn stimulate the activity of other AB replicators. In model 2, B is lost and sometimes replaced by a more active mutant form, B'. In model 3, B becomes damaged and loses activity, and its instability allows it to be replaced by a new B. For each model, conditions were identified where instability of B was detrimental, and where instability of B was beneficial. The results are consistent with the hypothesis that antagonistic pleiotropy can promote molecular instability and system complexity, and provide further support for a model linking aging and evolution.

An experimental test of temperature-dependent selection on mitochondrial haplotypes in Callosobruchus maculatus seed beetles

Mitochondrial DNA (mtDNA) consists of few but vital maternally inherited genes that interact closely with nuclear genes to produce cellular energy. How important mtDNA polymorphism is for adaptation is still unclear. The assumption in population genetic studies is often that segregating mtDNA variation is selectively neutral. This contrasts with empirical observations of mtDNA haplotypes affecting fitness-related traits and thermal sensitivity, and latitudinal clines in mtDNA haplotype frequencies. Here, we experimentally test whether ambient temperature affects selection on mtDNA variation, and whether such thermal effects are influenced by intergenomic epistasis due to interactions between mitochondrial and nuclear genes, using replicated experimental evolution in Callosobruchus maculatus seed beetle populations seeded with a mixture of different mtDNA haplotypes. We also test for sex-specific consequences of mtDNA evolution on reproductive success, given that mtDNA mutations can have sexually antagonistic fitness effects. Our results demonstrate natural selection on mtDNA haplotypes, with some support for thermal environment influencing mtDNA evolution through mitonuclear epistasis. The changes in male and female reproductive fitness were both aligned with changes in mtDNA haplotype frequencies, suggesting that natural selection on mtDNA is sexually concordant in stressful thermal environments. We discuss the implications of our findings for the evolution of mtDNA.

How Do Cells Adapt? Stories Told in Landscapes

Cells adapt to changing environments. Perturb a cell and it returns to a point of homeostasis. Perturb a population and it evolves toward a fitness peak. We review quantitative models of the forces of adaptation and their visualizations on landscapes. While some adaptations result from single mutations or few-gene effects, others are more cooperative, more delocalized in the genome, and more universal and physical. For example, homeostasis and evolution depend on protein folding and aggregation, energy and protein production, protein diffusion, molecular motor speeds and efficiencies, and protein expression levels. Models provide a way to learn about the fitness of cells and cell populations by making and testing hypotheses.

The Role of Conformational Dynamics and Allostery in Modulating Protein Evolution

Advances in sequencing techniques and statistical methods have made it possible not only to predict sequences of ancestral proteins but also to identify thousands of mutations in the human exome, some of which are disease associated. These developments have motivated numerous theories and raised many questions regarding the fundamental principles behind protein evolution, which have been traditionally investigated horizontally using the tip of the phylogenetic tree through comparative studies of extant proteins within a family. In this article, we review a vertical comparison of the modern and resurrected ancestral proteins. We focus mainly on the dynamical properties responsible for a protein's ability to adapt new functions in response to environmental changes. Using the Dynamic Flexibility Index and the Dynamic Coupling Index to quantify the relative flexibility and dynamic coupling at a site-specific, single-amino-acid level, we provide evidence that the migration of hinges, which are often functionally critical rigid sites, is a mechanism through which proteins can rapidly evolve. Additionally, we show that disease-associated mutations in proteins often result in flexibility changes even at positions distal from mutational sites, particularly in the modulation of active site dynamics.

Opposing Pressures of Speed and Efficiency Guide the Evolution of Molecular Machines

Many biomolecular machines need to be both fast and efficient. How has evolution optimized these machines along the tradeoff between speed and efficiency? We explore this question using optimizable dynamical models along coordinates that are plausible evolutionary degrees of freedom. Data on 11 motors and ion pumps are consistent with the hypothesis that evolution seeks an optimal balance of speed and efficiency, where any further small increase in one of these quantities would come at great expense to the other. For F_oF₁-ATPases in different species, we also find apparent optimization of the number of subunits in the c-ring, which determines the number of protons pumped per ATP synthesized. Interestingly, these ATPases appear to more optimized for efficiency than for speed, which can be rationalized through their key role as energy transducers in biology. The present modeling shows how the dynamical performance properties of biomolecular motors and pumps may have evolved to suit their corresponding biological actions.

Evolution of interface binding strengths in simplified model of protein quaternary structure

The self-assembly of proteins into protein quaternary structures is of fundamental importance to many biological processes, and protein misassembly is responsible for a wide range of proteopathic diseases. In recent years, abstract lattice models of protein self-assembly have been used to simulate the evolution and assembly of protein quaternary structure, and to provide a tractable way to study the genotype-phenotype map of such systems. Here we generalize these models by representing the interfaces as mutable binary strings. This simple change enables us to model the evolution of interface strengths, interface symmetry, and deterministic assembly pathways. Using the generalized model we are able to reproduce two important results established for real protein complexes: The first is that protein assembly pathways are under evolutionary selection to minimize misassembly. The second is that the assembly pathway of a complex mirrors its evolutionary history, and that both can be derived from the relative strengths of interfaces. These results demonstrate that the generalized lattice model offers a powerful new idealized framework to facilitate the study of protein self-assembly processes and their evolution.