Protein stability imposes limits on organism complexity and speed of molecular evolution

Distribution of specific prokaryotic immune systems correlates with host optimal growth temperature

Prokaryotes encode an arsenal of highly diverse immune systems to protect themselves against invading nucleic acids such as viruses, plasmids and transposons. This includes invader-interfering systems that neutralize invaders to protect their host, and abortive-infection systems, which trigger dormancy or cell death in their host to offer population-level immunity. Most prokaryotic immune systems are found across different environments and prokaryotic phyla, but their distribution appears biased and the factors that influence their distribution are largely unknown. Here, we compared and combined the prokaryotic immune system identification tools DefenseFinder and PADLOC to obtain an expanded view of the immune system arsenal. Our results show that the number of immune systems encoded is positively correlated with genome size and that the distribution of specific immune systems is linked to phylogeny. Furthermore, we reveal that certain invader-interfering systems are more frequently encoded by hosts with a relatively high optimum growth temperature, while abortive-infection systems are generally more frequently encoded by hosts with a relatively low optimum growth temperature. Combined, our study reveals several factors that correlate with differences in the distribution of prokaryotic immune systems and extends our understanding of how prokaryotes protect themselves from invaders in different environments.

Phylogenetic variations in a novel family of hyperstable apple snail egg proteins: insights into structural stability and functional trends

The relationship between protein stability and functional evolution is little explored in proteins purified from natural sources. Here, we investigated a novel family of egg proteins (Perivitellin-1, PV1) from Pomacea snails. Their remarkable stability and clade-related functions in most derived clades (Canaliculata and Bridgesii) make them excellent candidates for exploring this issue. To that aim, we studied PV1 (PpaPV1) from the most basal lineage, Flagellata. PpaPV1 displays unparalleled structural and kinetic stability, surpassing PV1s from derived clades, ranking among the most hyperstable proteins documented in nature. Its spectral features contribute to a pale egg coloration, exhibiting a milder glycan binding lectin activity with a narrower specificity than PV1s from the closely related Bridgesii clade. These findings provide evidence for substantial structural and functional changes throughout the genus' PV1 evolution. We observed that structural and kinetic stability decreased in a clade-related fashion and was associated with large variations in defensive traits. For instance, pale PpaPV1 lectin turns potent in the Bridgesii clade, adversely affecting gut morphology, while giving rise to brightly colored PV1s providing eggs with a conspicuous, probably warning signal in the Canaliculata clade. This work provides a comprehensive comparative analysis of PV1s from various apple snail species within a phylogenetic framework, offering insights into the interplay among their structural features, stability profiles and functional roles. More broadly, our work provides one of the first examples from natural evolution showing the crucial link among protein structure, stability and evolution of new functions.

Analysis of proteins in the light of mutations

Proteins have evolved through mutations-amino acid substitutions-since life appeared on Earth, some 109 years ago. The study of these phenomena has been of particular significance because of their impact on protein stability, function, and structure. This study offers a new viewpoint on how the most recent findings in these areas can be used to explore the impact of mutations on protein sequence, stability, and evolvability. Preliminary results indicate that: (1) mutations can be viewed as sensitive probes to identify 'typos' in the amino-acid sequence, and also to assess the resistance of naturally occurring proteins to unwanted sequence alterations; (2) the presence of 'typos' in the amino acid sequence, rather than being an evolutionary obstacle, could promote faster evolvability and, in turn, increase the likelihood of higher protein stability; (3) the mutation site is far more important than the substituted amino acid in terms of the marginal stability changes of the protein, and (4) the unpredictability of protein evolution at the molecular level-by mutations-exists even in the absence of epistasis effects. Finally, the Darwinian concept of evolution "descent with modification" and experimental evidence endorse one of the results of this study, which suggests that some regions of any protein sequence are susceptible to mutations while others are not. This work contributes to our general understanding of protein responses to mutations and may spur significant progress in our efforts to develop methods to accurately forecast changes in protein stability, their propensity for metamorphism, and their ability to evolve.

Divergent folding-mediated epistasis among unstable membrane protein variants

Many membrane proteins are prone to misfolding, which compromises their functional expression at the plasma membrane. This is particularly true for the mammalian gonadotropin-releasing hormone receptor GPCRs (GnRHR). We recently demonstrated that evolutionary GnRHR modifications appear to have coincided with adaptive changes in cotranslational folding efficiency. Though protein stability is known to shape evolution, it is unclear how cotranslational folding constraints modulate the synergistic, epistatic interactions between mutations. We therefore compared the pairwise interactions formed by mutations that disrupt the membrane topology (V276T) or tertiary structure (W107A) of GnRHR. Using deep mutational scanning, we evaluated how the plasma membrane expression of these variants is modified by hundreds of secondary mutations. An analysis of 251 mutants in three genetic backgrounds reveals that V276T and W107A form distinct epistatic interactions that depend on both the severity and the mechanism of destabilization. V276T forms predominantly negative epistatic interactions with destabilizing mutations in soluble loops. In contrast, W107A forms positive interactions with mutations in both loops and transmembrane domains that reflect the diminishing impacts of the destabilizing mutations in variants that are already unstable. These findings reveal how epistasis is remodeled by conformational defects in membrane proteins and in unstable proteins more generally.

Enhancing predictions of protein stability changes induced by single mutations using MSA-based Language Models

MOTIVATION

Protein Language Models offer a new perspective for addressing challenges in structural biology, while relying solely on sequence information. Recent studies have investigated their effectiveness in forecasting shifts in thermodynamic stability caused by single amino acid mutations, a task known for its complexity due to the sparse availability of data, constrained by experimental limitations. To tackle this problem, we introduce two key novelties: leveraging a Protein Language Model that incorporates Multiple Sequence Alignments to capture evolutionary information, and using a recently released mega-scale dataset with rigorous data pre-processing to mitigate overfitting.

RESULTS

We ensure comprehensive comparisons by fine-tuning various pre-trained models, taking advantage of analyses such as ablation studies and baselines evaluation. Our methodology introduces a stringent policy to reduce the widespread issue of data leakage, rigorously removing sequences from the training set when they exhibit significant similarity with the test set. The MSA Transformer emerges as the most accurate among the models under investigation, given its capability to leverage co-evolution signals encoded in aligned homologous sequences. Moreover, the optimized MSA Transformer outperforms existing methods and exhibits enhanced generalization power, leading to a notable improvement in predicting changes in protein stability resulting from point mutations.

AVAILABILITY AND IMPLEMENTATION

Code and data at https://github.com/RitAreaSciencePark/PLM4Muts.

SUPPLEMENTARY INFORMATION

Supplementary Information is available at Bioinformatics online.

Concepts and methods for predicting viral evolution

The seasonal human influenza virus undergoes rapid evolution, leading to significant changes in circulating viral strains from year to year. These changes are typically driven by adaptive mutations, particularly in the antigenic epitopes, the regions of the viral surface protein haemagglutinin targeted by human antibodies. Here we describe a consistent set of methods for data-driven predictive analysis of viral evolution. Our pipeline integrates four types of data: (1) sequence data of viral isolates collected on a worldwide scale, (2) epidemiological data on incidences, (3) antigenic characterization of circulating viruses, and (4) intrinsic viral phenotypes. From the combined analysis of these data, we obtain estimates of relative fitness for circulating strains and predictions of clade frequencies for periods of up to one year. Furthermore, we obtain comparative estimates of protection against future viral populations for candidate vaccine strains, providing a basis for pre-emptive vaccine strain selection. Continuously updated predictions obtained from the prediction pipeline for influenza and SARS-CoV-2 are available on the website previr.app .

Concepts and methods for predicting viral evolution

The seasonal human influenza virus undergoes rapid evolution, leading to significant changes in circulating viral strains from year to year. These changes are typically driven by adaptive mutations, particularly in the antigenic epitopes, the regions of the viral surface protein haemagglutinin targeted by human antibodies. Here we describe a consistent set of methods for data-driven predictive analysis of viral evolution. Our pipeline integrates four types of data: (1) sequence data of viral isolates collected on a worldwide scale, (2) epidemiological data on incidences, (3) antigenic characterization of circulating viruses, and (4) intrinsic viral phenotypes. From the combined analysis of these data, we obtain estimates of relative fitness for circulating strains and predictions of clade frequencies for periods of up to one year. Furthermore, we obtain comparative estimates of protection against future viral populations for candidate vaccine strains, providing a basis for pre-emptive vaccine strain selection. Continuously updated predictions obtained from the prediction pipeline for influenza and SARS-CoV-2 are available on the website previr.app.

Decoding the Role of the Global Proteome Dynamics for Cellular Thermal Stability

Molecular mechanisms underlying the thermal response of cells remain elusive. On the basis of the recent result that the short-time diffusive dynamics of the Escherichia coli proteome is an excellent indicator of temperature-dependent bacterial metabolism and death, we used neutron scattering (NS) spectroscopy and molecular dynamics (MD) simulations to investigate the sub-nanosecond proteome mobility in psychro-, meso-, and hyperthermophilic bacteria over a wide temperature range. The magnitude of thermal fluctuations, measured by atomic mean square displacements, is similar among all studied bacteria at their respective thermal cell death. Global roto-translational motions turn out to be the main factor distinguishing the bacterial dynamical properties. We ascribe this behavior to the difference in the average proteome net charge, which becomes less negative for increasing bacterial thermal stability. We propose that the chemical-physical properties of the cytoplasm and the global dynamics of the resulting proteome are fine-tuned by evolution to uphold optimal thermal stability conditions.

Identification of the SIRT1 gene's most harmful non-synonymous SNPs and their effects on functional and structural features-an in silico analysis

Introduction

The sirtuin (Silent mating type information regulation 2 homolog)1(SIRT1) protein plays a vital role in many disorders such as diabetes, cancer, obesity, inflammation, and neurodegenerative and cardiovascular diseases. The objective of this in silico analysis of SIRT1's functional single nucleotide polymorphisms (SNPs) was to gain valuable insight into the harmful effects of non-synonymous SNPs (nsSNPs) on the protein. The objective of the study was to use bioinformatics methods to investigate the genetic variations and modifications that may have an impact on the SIRT1 gene's expression and function.

Methods

nsSNPs of SIRT1 protein were collected from the dbSNP site, from its three (3) different protein accession IDs. These were then fed to various bioinformatic tools such as SIFT, Provean, and I- Mutant to find the most deleterious ones. Functional and structural effects were examined using the HOPE server and I-Tasser. Gene interactions were predicted by STRING software. The SIFT, Provean, and I-Mutant tools detected the most deleterious three nsSNPs (rs769519031, rs778184510, and rs199983221).

Results

Out of 252 nsSNPs, SIFT analysis showed that 94 were deleterious, Provean listed 67 dangerous, and I-Mutant found 58 nsSNPs resulting in lowered stability of proteins. HOPE modelling of rs199983221 and rs769519031 suggested reduced hydrophobicity due to Ile 4Thr and Ile223Ser resulting in decreased hydrophobic interactions. In contrast, on modelling rs778184510, the mutant protein had a higher hydrophobicity than the wild type.

Conclusions

Our study reports that three nsSNPs (D357A, I223S, I4T) are the most damaging mutations of the SIRT1 gene. Mutations may result in altered protein structure and functions. Such altered protein may be the basis for various disorders. Our findings may be a crucial guide in establishing the pathogenesis of various disorders.

Phylogenetic variations in a novel family of hyperstable apple snail egg proteins: insights into structural stability and functional trends

The relationship between protein stability and function evolution has not been explored in proteins from natural sources. Here, we investigate the phylogenetic differences of Perivitellin-1 (PV1) a novel family of hyperstable egg carotenoproteins crucial to the reproductive success of Pomacea snails, as they have evolved clade-specific protective functions. We studied P. patula PV1 (PpaPV1) from Flagellata clade eggs, the most basal of Pomacea and compared it with PV1s orthologs from derived clades. PpaPV1 stands as the most stable, with longer unfolding half-life, resistance to detergent unfolding, and therefore higher kinetic stability than PV1s from derived clades. In fact, PpaPV1 is among the most hyperstable proteins described in nature. In addition, its spectral characteristics providing a pale egg coloration, mild lectin activity and glycan specificity are narrower than derived clades. Our results provide evidence indicating large structural and functional changes along the evolution of the genus. Notably, the lectin binding of PpaPV1 is less pronounced, and its glycan specificity is narrower compared to PV1s in the sister Bridgesii clade. Our findings underscore the phylogenetic disparities in terms of structural and kinetic stability, as well as defensive traits like a potent lectin activity affecting the gut morphology of potential predators within the Bridgesii clade or a conspicuous, likely warning coloration, within the Canaliculata clade. This work provides a comprehensive comparison of the structural attributes, stability profiles, and functional roles of apple snail egg PV1s from multiple species within a phylogenetic context. Furthermore, it proposes an evolutionary hypothesis suggesting a trade-off between structural stability and the functional aspects of apple snail's major egg defense protein.

GENERALIST: A latent space based generative model for protein sequence families

Generative models of protein sequence families are an important tool in the repertoire of protein scientists and engineers alike. However, state-of-the-art generative approaches face inference, accuracy, and overfitting- related obstacles when modeling moderately sized to large proteins and/or protein families with low sequence coverage. Here, we present a simple to learn, tunable, and accurate generative model, GENERALIST: GENERAtive nonLInear tenSor-factorizaTion for protein sequences. GENERALIST accurately captures several high order summary statistics of amino acid covariation. GENERALIST also predicts conservative local optimal sequences which are likely to fold in stable 3D structure. Importantly, unlike current methods, the density of sequences in GENERALIST-modeled sequence ensembles closely resembles the corresponding natural ensembles. Finally, GENERALIST embeds protein sequences in an informative latent space. GENERALIST will be an important tool to study protein sequence variability.

Viral Evolution Shaped by Host Proteostasis Networks

Understanding the factors that shape viral evolution is critical for developing effective antiviral strategies, accurately predicting viral evolution, and preventing pandemics. One fundamental determinant of viral evolution is the interplay between viral protein biophysics and the host machineries that regulate protein folding and quality control. Most adaptive mutations in viruses are biophysically deleterious, resulting in a viral protein product with folding defects. In cells, protein folding is assisted by a dynamic system of chaperones and quality control processes known as the proteostasis network. Host proteostasis networks can determine the fates of viral proteins with biophysical defects, either by assisting with folding or by targeting them for degradation. In this review, we discuss and analyze new discoveries revealing that host proteostasis factors can profoundly shape the sequence space accessible to evolving viral proteins. We also discuss the many opportunities for research progress proffered by the proteostasis perspective on viral evolution and adaptation.

Protein folding rate evolution upon mutations

Despite the spectacular success of cutting-edge protein fold prediction methods, many critical questions remain unanswered, including why proteins can reach their native state in a biologically reasonable time. A satisfactory answer to this simple question could shed light on the slowest folding rate of proteins as well as how mutations-amino-acid substitutions and/or post-translational modifications-might affect it. Preliminary results indicate that (i) Anfinsen's dogma validity ensures that proteins reach their native state on a reasonable timescale regardless of their sequence or length, and (ii) it is feasible to determine the evolution of protein folding rates without accounting for epistasis effects or the mutational trajectories between the starting and target sequences. These results have direct implications for evolutionary biology because they lay the groundwork for a better understanding of why, and to what extent, mutations-a crucial element of evolution and a factor influencing it-affect protein evolvability. Furthermore, they may spur significant progress in our efforts to solve crucial structural biology problems, such as how a sequence encodes its folding.

Diffusive Dynamics of Bacterial Proteome as a Proxy of Cell Death

Temperature variations have a big impact on bacterial metabolism and death, yet an exhaustive molecular picture of these processes is still missing. For instance, whether thermal death is determined by the deterioration of the whole or a specific part of the proteome is hotly debated. Here, by monitoring the proteome dynamics of E. coli, we clearly show that only a minor fraction of the proteome unfolds at the cell death. First, we prove that the dynamical state of the E. coli proteome is an excellent proxy for temperature-dependent bacterial metabolism and death. The proteome diffusive dynamics peaks at about the bacterial optimal growth temperature, then a dramatic dynamical slowdown is observed that starts just below the cell's death temperature. Next, we show that this slowdown is caused by the unfolding of just a small fraction of proteins that establish an entangling interprotein network, dominated by hydrophobic interactions, across the cytoplasm. Finally, the deduced progress of the proteome unfolding and its diffusive dynamics are both key to correctly reproduce the E. coli growth rate.

Design of an artificial phage-display library based on a new scaffold improved for average stability of the randomized proteins

Scaffold-based protein libraries are designed to be both diverse and rich in functional/folded proteins. However, introducing an extended diversity while preserving stability of the initial scaffold remains a challenge. Here we developed an original approach to select the ensemble of folded proteins from an initial library. The thermostable CheY protein from Thermotoga maritima was chosen as scaffold. Four loops of CheY were diversified to create a new binding surface. The subset of the library giving rise to folded proteins was first selected using a natural protein partner of the template scaffold. Then, a gene shuffling approach based on a single restriction enzyme was used to recombine DNA sequences encoding these filtrated variants. Taken together, the filtration strategy and the shuffling of the filtrated sequences were shown to enrich the library in folded and stable sequences while maintaining a large diversity in the final library (Lib-Cheytins 2.1). Binders of the Oplophorus luciferase Kaz domain were then selected by phage display from the final library, showing affinities in the μM range. One of the best variants induced a loss of 92% of luminescent activity, suggesting that this Cheytin preferentially binds to the Kaz active site.

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.

The evolution of the HIV-1 protease folding stability

The evolution of structural proteins is generally constrained by the folding stability. However, little is known about the particular capacity of viral proteins to accommodate mutations that can potentially affect the protein stability and, in general, the evolution of the protein stability over time. As an illustrative model case, here, we investigated the evolution of the stability of the human immunodeficiency virus (HIV-1) protease (PR), which is a common HIV-1 drug target, under diverse evolutionary scenarios that include (1) intra-host virus evolution in a cohort of seventy-five patients sampled over time, (2) intra-host virus evolution sampled before and after specific PR-based treatments, and (3) inter-host evolution considering extant and ancestral (reconstructed) PR sequences from diverse HIV-1 subtypes. We also investigated the specific influence of currently known HIV-1 PR resistance mutations on the PR folding stability. We found that the HIV-1 PR stability fluctuated over time within a constant and wide range in any studied evolutionary scenario, accommodating multiple mutations that partially affected the stability while maintaining activity. We did not identify relationships between change of PR stability and diverse clinical parameters such as viral load, CD4⁺ T-cell counts, and a surrogate of time from infection. Counterintuitively, we predicted that nearly half of the studied HIV-1 PR resistance mutations do not significantly decrease stability, which, together with compensatory mutations, would allow the protein to adapt without requiring dramatic stability changes. We conclude that the HIV-1 PR presents a wide structural plasticity to acquire molecular adaptations without affecting the overall evolution of stability.

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.

Increasing protein stability by inferring substitution effects from high-throughput experiments

We apply a computational model, global multi-mutant analysis (GMMA), to inform on effects of most amino acid substitutions from a randomly mutated gene library. Using a high mutation frequency, the method can determine mutations that increase the stability of even very stable proteins for which conventional selection systems have reached their limit. As a demonstration of this, we screened a mutant library of a highly stable and computationally redesigned model protein using an in vivo genetic sensor for folding and assigned a stability effect to 374 of 912 possible single amino acid substitutions. Combining the top 9 substitutions increased the unfolding energy 47 to 69 kJ/mol in a single engineering step. Crystal structures of stabilized variants showed small perturbations in helices 1 and 2, which rendered them closer in structure to the redesign template. This case study illustrates the capability of the method, which is applicable to any screen for protein function.

Stability and expression of SARS-CoV-2 spike-protein mutations

Protein fold stability likely plays a role in SARS-CoV-2 S-protein evolution, together with ACE2 binding and antibody evasion. While few thermodynamic stability data are available for S-protein mutants, many systematic experimental data exist for their expression. In this paper, we explore whether such expression levels relate to the thermodynamic stability of the mutants. We studied mutation-induced SARS-CoV-2 S-protein fold stability, as computed by three very distinct methods and eight different protein structures to account for method- and structure-dependencies. For all methods and structures used (24 comparisons), computed stability changes correlate significantly (99% confidence level) with experimental yeast expression from the literature, such that higher expression is associated with relatively higher fold stability. Also significant, albeit weaker, correlations were seen between stability and ACE2 binding effects. The effect of thermodynamic fold stability may be direct or a correlate of amino acid or site properties, notably the solvent exposure of the site. Correlation between computed stability and experimental expression and ACE2 binding suggests that functional properties of the SARS-CoV-2 S-protein mutant space are largely determined by a few simple features, due to underlying correlations. Our study lends promise to the development of computational tools that may ideally aid in understanding and predicting SARS-CoV-2 S-protein evolution.

Evolution avoids a pathological stabilizing interaction in the immune protein S100A9

Stability constrains evolution. While much is known about constraints on destabilizing mutations, less is known about the constraints on stabilizing mutations. We recently identified a mutation in the innate immune protein S100A9 that provides insight into such constraints. When introduced into human S100A9, M63F simultaneously increases the stability of the protein and disrupts its natural ability to activate Toll-like receptor 4. Using chemical denaturation, we found that M63F stabilizes a calcium-bound conformation of hS100A9. We then used NMR to solve the structure of the mutant protein, revealing that the mutation distorts the hydrophobic binding surface of hS100A9, explaining its deleterious effect on function. Hydrogen-deuterium exchange (HDX) experiments revealed stabilization of the region around M63F in the structure, notably Phe37. In the structure of the M63F mutant, the Phe37 and Phe63 sidechains are in contact, plausibly forming an edge-face π-stack. Mutating Phe37 to Leu abolished the stabilizing effect of M63F as probed by both chemical denaturation and HDX. It also restored the biological activity of S100A9 disrupted by M63F. These findings reveal that Phe63 creates a molecular staple with Phe37 that stabilizes a nonfunctional conformation of the protein, thus disrupting function. Using a bioinformatic analysis, we found that S100A9 proteins from different organisms rarely have Phe at both positions 37 and 63, suggesting that avoiding a pathological stabilizing interaction indeed constrains S100A9 evolution. This work highlights an important evolutionary constraint on stabilizing mutations, namely, that they must avoid inappropriately stabilizing nonfunctional protein conformations.

In Silico Analysis of PORD Mutations on the 3D Structure of P450 Oxidoreductase

Cytochrome P450 oxidoreductase (POR) is a membrane-bound flavoprotein that helps in transferring electrons from its NADPH domain to all cytochrome P450 (CYP450) enzymes. Mutations in the POR gene could severely affect the metabolism of steroid hormones and the development of skeletal muscles, a condition known as Cytochrome P450 oxidoreductase deficiency (PORD). PORD is associated with clinical presentations of disorders of sex development, Antley and Bixler’s syndrome (ABS), as well as an abnormal steroid hormone profile. We have performed an in silico analysis of POR 3D X-ray protein crystal structure to study the effects of reported mutations on the POR enzyme structure. A total of 32 missense mutations were identified, from 170 PORD patients, and mapped on the 3D crystal structure of the POR enzyme. In addition, five of the missense mutations (R457H, A287P, D210G, Y181D and Y607C) were further selected for an in-depth in silico analysis to correlate the observed changes in POR protein structure with the clinical phenotypes observed in PORD patients. Overall, missense mutations found in the binding sites of POR cofactors could lead to a severe form of PORD, emphasizing the importance of POR cofactor binding domains in transferring electrons to the CYP450 enzyme family.

Proteins' Evolution upon Point Mutations

As the reader must be already aware, state-of-the-art protein folding prediction methods have reached a smashing success in their goal of accurately determining the three-dimensional structures of proteins. Yet, a solution to simple problems such as the effects of protein point mutations on their (i) native conformation; (ii) marginal stability; (iii) ensemble of high-energy nativelike conformations; and (iv) metamorphism propensity and, hence, their evolvability, remains as an unsolved problem. As a plausible solution to the latter, some properties of the amide hydrogen-deuterium exchange, a highly sensitive probe of the structure, stability, and folding of proteins, are assessed from a new perspective. The preliminary results indicate that the protein marginal stability change upon point mutations provides the necessary and sufficient information to estimate, through a Boltzmann factor, the evolution of the amide hydrogen exchange protection factors and, consequently, that of the ensemble of folded conformations coexisting with the native state. This work contributes to our general understanding of the effects of point mutations on proteins and may spur significant progress in our efforts to develop methods to determine the appearance of new folds and functions accurately.

Protein folding stabilities are a major determinant of oxidation rates for buried methionine residues

The oxidation of protein-bound methionines to form methionine sulfoxides has a broad range of biological ramifications, making it important to delineate factors that influence methionine oxidation rates within a given protein. This is especially important for biopharmaceuticals, where oxidation can lead to deactivation and degradation. Previously, neighboring residue effects and solvent accessibility have been shown to impact the susceptibility of methionine residues to oxidation. In this study, we provide proteome-wide evidence that oxidation rates of buried methionine residues are also strongly influenced by the thermodynamic folding stability of proteins. We surveyed the Escherichia coli proteome using several proteomic methodologies and globally measured oxidation rates of methionine residues in the presence and absence of tertiary structure, as well as the folding stabilities of methionine-containing domains. These data indicated that buried methionines have a wide range of protection factors against oxidation that correlate strongly with folding stabilities. Consistent with this, we show that in comparison to E. coli, the proteome of the thermophile Thermus thermophilus is significantly more stable and thus more resistant to methionine oxidation. To demonstrate the utility of this correlation, we used native methionine oxidation rates to survey the folding stabilities of E. coli and T. thermophilus proteomes at various temperatures and propose a model that relates the temperature dependence of the folding stabilities of these two species to their optimal growth temperatures. Overall, these results indicate that oxidation rates of buried methionines from the native state of proteins can be used as a metric of folding stability.

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.

Main protease mutants of SARS-CoV-2 variants remain susceptible to nirmatrelvir

The COVID-19 pandemic continues to be a public health threat. Multiple mutations in the spike protein of emerging variants of SARS-CoV-2 appear to impact on the effectiveness of available vaccines. Specific antiviral agents are keenly anticipated but their efficacy may also be compromised in emerging variants. One of the most attractive coronaviral drug targets is the main protease (M^pro). A promising M^pro inhibitor of clinical relevance is the peptidomimetic nirmatrelvir (PF-07321332). We expressed M^pro of six SARS-CoV-2 lineages (C.37 Lambda, B.1.1.318, B.1.2, B.1.351 Beta, B.1.1.529 Omicron, P.2 Zeta), each of which carries a strongly prevalent missense mutation (G15S, T21I, L89F, K90R, P132H, L205V). Enzyme kinetics reveal that these M^pro variants are catalytically competent to a similar degree as the wildtype. We show that nirmatrelvir has similar potency against the variants as the wildtype. Our in vitro data suggest that the efficacy of the specific M^pro inhibitor nirmatrelvir is not compromised in current COVID-19 variants.

Inferring the distribution of fitness effects in patient-sampled and experimental virus populations: two case studies

We here propose an analysis pipeline for inferring the distribution of fitness effects (DFE) from either patient-sampled or experimentally-evolved viral populations, that explicitly accounts for non-Wright-Fisher and non-equilibrium population dynamics inherent to pathogens. We examine the performance of this approach via extensive power and performance analyses, and highlight two illustrative applications - one from an experimentally-passaged RNA virus, and the other from a clinically-sampled DNA virus. Finally, we discuss how such DFE inference may shed light on major research questions in virus evolution, ranging from a quantification of the population genetic processes governing genome size, to the role of Hill-Robertson interference in dictating adaptive outcomes, to the potential design of novel therapeutic approaches to eradicate within-patient viral populations via induced mutational meltdown.

Data set and fitting dependencies when estimating protein mutant stability: Toward simple, balanced, and interpretable models

Accurate prediction of protein stability changes upon mutation (ΔΔG) is increasingly important to evolution studies, protein engineering, and screening of disease-causing gene variants but is challenged by biases in training data. We investigated 45 linear regression models trained on data sets that account systematically for destabilization bias and mutation-type bias B_M . The models were externally validated on three test data sets probing different pathologies and for internal consistency (symmetry and neutrality). Model structure and performance substantially depended on training data and even fitting method. We developed two final models: SimBa-IB for typical natural mutations and SimBa-SYM for situations where stabilizing and destabilizing mutations occur to a similar extent. SimBa-SYM, despite is simplicity, is essentially non-biased (vs. the S^sym data set) while still performing well for all data sets (R ~ 0.46-0.54, MAE = 1.16-1.24 kcal/mol). The simple models provide advantage in terms of interpretability, use and future improvement, and are freely available on GitHub.

Structure and Mutations of SARS-CoV-2 Spike Protein: A Focused Overview

The spike protein (S-protein) of SARS-CoV-2, the protein that enables the virus to infect human cells, is the basis for many vaccines and a hotspot of concerning virus evolution. Here, we discuss the outstanding progress in structural characterization of the S-protein and how these structures facilitate analysis of virus function and evolution. We emphasize the differences in reported structures and that analysis of structure-function relationships is sensitive to the structure used. We show that the average residue solvent exposure in nearly complete structures is a good descriptor of open vs closed conformation states. Because of structural heterogeneity of functionally important surface-exposed residues, we recommend using averages of a group of high-quality protein structures rather than a single structure before reaching conclusions on specific structure-function relationships. To illustrate these points, we analyze some significant chemical tendencies of prominent S-protein mutations in the context of the available structures. In the discussion of new variants, we emphasize the selectivity of binding to ACE2 vs prominent antibodies rather than simply the antibody escape or ACE2 affinity separately. We note that larger chemical changes, in particular increased electrostatic charge or side-chain volume of exposed surface residues, are recurring in mutations of concern, plausibly related to adaptation to the negative surface potential of human ACE2. We also find indications that the fixated mutations of the S-protein in the main variants are less destabilizing than would be expected on average, possibly pointing toward a selection pressure on the S-protein. The richness of available structures for all of these situations provides an enormously valuable basis for future research into these structure-function relationships.

Is Protein Folding a Thermodynamically Unfavorable, Active, Energy-Dependent Process?

The prevailing current view of protein folding is the thermodynamic hypothesis, under which the native folded conformation of a protein corresponds to the global minimum of Gibbs free energy G. We question this concept and show that the empirical evidence behind the thermodynamic hypothesis of folding is far from strong. Furthermore, physical theory-based approaches to the prediction of protein folds and their folding pathways so far have invariably failed except for some very small proteins, despite decades of intensive theory development and the enormous increase of computer power. The recent spectacular successes in protein structure prediction owe to evolutionary modeling of amino acid sequence substitutions enhanced by deep learning methods, but even these breakthroughs provide no information on the protein folding mechanisms and pathways. We discuss an alternative view of protein folding, under which the native state of most proteins does not occupy the global free energy minimum, but rather, a local minimum on a fluctuating free energy landscape. We further argue that ΔG of folding is likely to be positive for the majority of proteins, which therefore fold into their native conformations only through interactions with the energy-dependent molecular machinery of living cells, in particular, the translation system and chaperones. Accordingly, protein folding should be modeled as it occurs in vivo, that is, as a non-equilibrium, active, energy-dependent process.

Surface residues and non-additive interactions stabilize a consensus homeodomain protein

Despite the widely reported success of consensus design in producing highly stabilized proteins, little is known about the physical mechanisms underlying this stabilization. Here we explore the potential sources of stabilization by performing a systematic analysis of the 29 substitutions that we previously found to collectively stabilize a consensus homeodomain compared to an extant homeodomain. By separately introducing groups of consensus substitutions that alter or preserve charge state, occur at varying degrees of residue burial, and occur at positions of varying degrees of conservation, we determine the extent to which these three features contribute to the consensus stability enhancement. Surprisingly, we find that the largest total contribution to stability comes from consensus substitutions on the protein surface and that the largest per-substitution contributions come from substitutions that maintain charge state. This finding suggests that although consensus proteins are often enriched in charged residues, consensus stabilization does not result primarily from interactions involving charged residues. Although consensus substitutions at strongly conserved positions also contribute disproportionately to stabilization, significant stabilization is also contributed from substitutions at weakly conserved positions. Furthermore, we find that identical consensus substitutions show larger stabilizing effects when introduced into the consensus background than when introduced into an extant homeodomain, indicating that synergistic, stabilizing interactions among the consensus residues contribute to consensus stability enhancement of the homeodomain. By measuring DNA binding affinity for the same set of variants, we find that although consensus design of the homeodomain increases both affinity and folding stability, it does so using a largely non-overlapping set of substitutions.

Quantifying the impact of changes in effective population size and expression level on the rate of coding sequence evolution

Molecular sequences are shaped by selection, where the strength of selection relative to drift is determined by effective population size (N_e). Populations with high N_e are expected to undergo stronger purifying selection, and consequently to show a lower substitution rate for selected mutations relative to the substitution rate for neutral mutations (ω). However, computational models based on biophysics of protein stability have suggested that ω can also be independent of N_e. Together, the response of ω to changes in N_e depends on the specific mapping from sequence to fitness. Importantly, an increase in protein expression level has been found empirically to result in decrease of ω, an observation predicted by theoretical models assuming selection for protein stability. Here, we derive a theoretical approximation for the response of ω to changes in N_e and expression level, under an explicit genotype-phenotype-fitness map. The method is generally valid for additive traits and log-concave fitness functions. We applied these results to protein undergoing selection for their conformational stability and corroborate out findings with simulations under more complex models. We predict a weak response of ω to changes in either N_e or expression level, which are interchangeable. Based on empirical data, we propose that fitness based on the conformational stability may not be a sufficient mechanism to explain the empirically observed variation in ω across species. Other aspects of protein biophysics might be explored, such as protein-protein interactions, which can lead to a stronger response of ω to changes in N_e.

Epistatic interactions promote persistence of NS3-Q80K in HCV infection by compensating for protein folding instability

The Q80K polymorphism in the NS3-4A protease of the hepatitis C virus is associated with treatment failure of direct-acting antiviral agents. This polymorphism is highly prevalent in genotype 1a infections and stably transmitted between hosts. Here, we investigated the underlying molecular mechanisms of evolutionarily conserved coevolving amino acids in NS3-Q80K and revealed potential implications of epistatic interactions in immune escape and variants persistence. Using purified protein, we characterized the impact of epistatic amino acid substitutions on the physicochemical properties and peptide cleavage kinetics of the NS3-Q80K protease. We found that Q80K destabilized the protease protein fold (p < 0.0001). Although NS3-Q80K showed reduced peptide substrate turnover (p < 0.0002), replicative fitness in an H77S.3 cell culture model of infection was not significantly inferior to the WT virus. Epistatic substitutions at residues 91 and 174 in NS3-Q80K stabilized the protein fold (p < 0.0001) and leveraged the WT protease stability. However, changes in protease stability inversely correlated with enzymatic activity. In infectious cell culture, these secondary substitutions were not associated with a gain of replicative fitness in NS3-Q80K variants. Using molecular dynamics, we observed that the total number of residue contacts in NS3-Q80K mutants correlated with protein folding stability. Changes in the number of contacts reflected the compensatory effect on protein folding instability by epistatic substitutions. In summary, epistatic substitutions in NS3-Q80K contribute to viral fitness by mechanisms not directly related to RNA replication. By compensating for protein-folding instability, epistatic interactions likely protect NS3-Q80K variants from immune cell recognition.

Intra-Population Competition during Adaptation to Increased Temperature in an RNA Bacteriophage

Evolution of RNA bacteriophages of the family Leviviridae is governed by the high error rates of their RNA-dependent RNA polymerases. This fact, together with their large population sizes, leads to the generation of highly heterogeneous populations that adapt rapidly to most changes in the environment. Throughout adaptation, the different mutants that make up a viral population compete with each other in a non-trivial process in which their selective values change over time due to the generation of new mutations. In this work we have characterised the intra-population dynamics of a well-studied levivirus, Qβ, when it is propagated at a higher-than-optimal temperature. Our results show that adapting populations experienced rapid changes that involved the ascent of particular genotypes and the loss of some beneficial mutations of early generation. Artificially reconstructed populations, containing a fraction of the diversity present in actual populations, fixed mutations more rapidly, illustrating how population bottlenecks may guide the adaptive pathways. The conclusion is that, when the availability of beneficial mutations under a particular selective condition is elevated, the final outcome of adaptation depends more on the occasional occurrence of population bottlenecks and how mutations combine in genomes than on the selective value of particular mutations.

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.

From genotypes to organisms: State-of-the-art and perspectives of a cornerstone in evolutionary dynamics

Understanding how genotypes map onto phenotypes, fitness, and eventually organisms is arguably the next major missing piece in a fully predictive theory of evolution. We refer to this generally as the problem of the genotype-phenotype map. Though we are still far from achieving a complete picture of these relationships, our current understanding of simpler questions, such as the structure induced in the space of genotypes by sequences mapped to molecular structures, has revealed important facts that deeply affect the dynamical description of evolutionary processes. Empirical evidence supporting the fundamental relevance of features such as phenotypic bias is mounting as well, while the synthesis of conceptual and experimental progress leads to questioning current assumptions on the nature of evolutionary dynamics-cancer progression models or synthetic biology approaches being notable examples. This work delves with a critical and constructive attitude into our current knowledge of how genotypes map onto molecular phenotypes and organismal functions, and discusses theoretical and empirical avenues to broaden and improve this comprehension. As a final goal, this community should aim at deriving an updated picture of evolutionary processes soundly relying on the structural properties of genotype spaces, as revealed by modern techniques of molecular and functional analysis.

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.

A meta-analysis of the activity, stability, and mutational characteristics of temperature-adapted enzymes

Understanding the characteristics that define temperature-adapted enzymes has been a major goal of extremophile enzymology in recent decades. In the present study, we explore these characteristics by comparing psychrophilic, mesophilic, and thermophilic enzymes. Through a meta-analysis of existing data, we show that psychrophilic enzymes exhibit a significantly larger gap (T_g) between their optimum and melting temperatures compared with mesophilic and thermophilic enzymes. These results suggest that T_g may be a useful indicator as to whether an enzyme is psychrophilic or not and that models of psychrophilic enzyme catalysis need to account for this gap. Additionally, by using predictive protein stability software, HoTMuSiC and PoPMuSiC, we show that the deleterious nature of amino acid substitutions to protein stability increases from psychrophiles to thermophiles. How this ultimately affects the mutational tolerance and evolutionary rate of temperature adapted organisms is currently unknown.

Three Simple Properties Explain Protein Stability Change upon Mutation

Accurate prediction of protein stability upon mutation enables rational engineering of new proteins and insights into protein evolution and monogenetic diseases caused by single-point amino acid substitutions. Many tools have been developed to this aim, ranging from energy-based models to machine-learning methods that use large amounts of experimental data. However, as the methods become more complex, the interpretation of the chemistry underlying the protein stability effects becomes obscure. It is thus of interest to identify the simplest prediction model that retains complete amino acid specific interpretation; for a given number of input descriptors, we expect such a model to be almost universal. In this study, we identify such a limiting model, SimBa, a simple multilinear regression model trained on a substitution-type-balanced experimental data set. The model accounts only for the solvent accessibility of the site, volume difference, and polarity difference caused by mutation. Our results show that this very simple and directly applicable model performs comparably to other much more complex, widely used protein stability prediction methods. This suggests that a hard limit of ∼1 kcal/mol numerical accuracy and an R ∼ 0.5 trend accuracy exists and that new features, such as account of unfolded states, water colocalization, and amino acid correlations, are required to improve accuracy to, e.g., 1/2 kcal/mol.

Nothing in Evolution Makes Sense Except in the Light of Biology

A key question in biology is the predictability of the evolutionary process. If we can correctly predict the outcome of evolution, we may be better equipped to anticipate and manage species' adaptation to climate change, habitat loss, invasive species, or emerging infectious diseases, as well as improve our basic understanding of the history of life on Earth. In the present article, we ask the questions when, why, and if the outcome of future evolution is predictable. We first define predictable and then discuss two conflicting views: that evolution is inherently unpredictable and that evolution is predictable given the ability to collect the right data. We identify factors that generate unpredictability, the data that might be required to make predictions at some level of precision or at a specific timescale, and the intellectual and translational value of understanding when prediction is or is not possible.

Thermal stability and structure of glyceraldehyde-3-phosphate dehydrogenase from the coral Acropora millepora

Corals are vulnerable to increasing ocean temperatures. It is known that elevated temperatures lead to the breakdown of an essential mutualistic relationship with photosynthetic algae. The molecular mechanisms of this temperature-dependent loss of symbiosis are less well understood. Here, the thermal stability of a critical metabolic enzyme, glyceraldehyde-3-phosphate dehydrogenase, from the stony coral Acropora millepora was found to increase significantly in the presence of its cofactor NAD+. Determination of the structure of the cofactor-enzyme complex (PDB ID 6PX2) revealed variable NAD+ occupancy across the four monomers of the tetrameric enzyme. The structure of the fully occupied monomers was compared to those with partial cofactor occupancy, identifying regions of difference that may account for the increased thermal stability.

A base measure of precision for protein stability predictors: structural sensitivity

BACKGROUND

Prediction of the change in fold stability (ΔΔG) of a protein upon mutation is of major importance to protein engineering and screening of disease-causing variants. Many prediction methods can use 3D structural information to predict ΔΔG. While the performance of these methods has been extensively studied, a new problem has arisen due to the abundance of crystal structures: How precise are these methods in terms of structure input used, which structure should be used, and how much does it matter? Thus, there is a need to quantify the structural sensitivity of protein stability prediction methods.

RESULTS

We computed the structural sensitivity of six widely-used prediction methods by use of saturated computational mutagenesis on a diverse set of 87 structures of 25 proteins. Our results show that structural sensitivity varies massively and surprisingly falls into two very distinct groups, with methods that take detailed account of the local environment showing a sensitivity of ~ 0.6 to 0.8 kcal/mol, whereas machine-learning methods display much lower sensitivity (~ 0.1 kcal/mol). We also observe that the precision correlates with the accuracy for mutation-type-balanced data sets but not generally reported accuracy of the methods, indicating the importance of mutation-type balance in both contexts.

CONCLUSIONS

The structural sensitivity of stability prediction methods varies greatly and is caused mainly by the models and less by the actual protein structural differences. As a new recommended standard, we therefore suggest that ΔΔG values are evaluated on three protein structures when available and the associated standard deviation reported, to emphasize not just the accuracy but also the precision of the method in a specific study. Our observation that machine-learning methods deemphasize structure may indicate that folded wild-type structures alone, without the folded mutant and unfolded structures, only add modest value for assessing protein stability effects, and that side-chain-sensitive methods overstate the significance of the folded wild-type structure.

Evolution of frustrated and stabilising contacts in reconstructed ancient proteins

Energetic properties of a protein are a major determinant of its evolutionary fitness. Using a reconstruction algorithm, dating the reconstructed proteins and calculating the interaction network between their amino acids through a coevolutionary approach, we studied how the interactions that stabilise 890 proteins, belonging to five families, evolved for billions of years. In particular, we focused our attention on the network of most strongly attractive contacts and on that of poorly optimised, frustrated contacts. Our results support the idea that the cluster of most attractive interactions extends its size along evolutionary time, but from the data, we cannot conclude that protein stability or that the degree of frustration tends always to decrease.

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.

About the Protein Space Vastness

An accurate estimation of the Protein Space size, in light of the factors that govern it, is a long-standing problem and of paramount importance in evolutionary biology, since it determines the nature of protein evolvability. A simple analysis will enable us to, firstly, reduce an unrealistic Protein Space size of ~ 10¹³⁰ sequences, for a 100-residues polypeptide chain, to ~ 10⁹ functional proteins and, secondly, estimate a robust average-mutation rate per amino acid (ξ ~ 1.23) and infer from it, in light of the protein marginal stability, that only a fraction of the sequence will be available at any one time for a functional protein to evolve. Although this result does not solve the Protein Space vastness problem frames it in a more rational one and illustrates the impact of the marginal stability on protein evolvability.

Exploring the sequence fitness landscape of a bridge between protein folds

Most foldable protein sequences adopt only a single native fold. Recent protein design studies have, however, created protein sequences which fold into different structures apon changes of environment, or single point mutation, the best characterized example being the switch between the folds of the GA and GB binding domains of streptococcal protein G. To obtain further insight into the design of sequences which can switch folds, we have used a computational model for the fitness landscape of a single fold, built from the observed sequence variation of protein homologues. We have recently shown that such coevolutionary models can be used to design novel foldable sequences. By appropriately combining two of these models to describe the joint fitness landscape of GA and GB, we are able to describe the propensity of a given sequence for each of the two folds. We have successfully tested the combined model against the known series of designed GA/GB hybrids. Using Monte Carlo simulations on this landscape, we are able to identify pathways of mutations connecting the two folds. In the absence of a requirement for domain stability, the most frequent paths go via sequences in which neither domain is stably folded, reminiscent of the propensity for certain intrinsically disordered proteins to fold into different structures according to context. Even if the folded state is required to be stable, we find that there is nonetheless still a wide range of sequences which are close to the transition region and therefore likely fold switches, consistent with recent estimates that fold switching may be more widespread than had been thought.

Sparse Epistatic Patterns in the Evolution of Terpene Synthases

We explore sequence determinants of enzyme activity and specificity in a major enzyme family of terpene synthases. Most enzymes in this family catalyze reactions that produce cyclic terpenes-complex hydrocarbons widely used by plants and insects in diverse biological processes such as defense, communication, and symbiosis. To analyze the molecular mechanisms of emergence of terpene cyclization, we have carried out in-depth examination of mutational space around (E)-β-farnesene synthase, an Artemisia annua enzyme which catalyzes production of a linear hydrocarbon chain. Each mutant enzyme in our synthetic libraries was characterized biochemically, and the resulting reaction rate data were used as input to the Michaelis-Menten model of enzyme kinetics, in which free energies were represented as sums of one-amino-acid contributions and two-amino-acid couplings. Our model predicts measured reaction rates with high accuracy and yields free energy landscapes characterized by relatively few coupling terms. As a result, the Michaelis-Menten free energy landscapes have simple, interpretable structure and exhibit little epistasis. We have also developed biophysical fitness models based on the assumption that highly fit enzymes have evolved to maximize the output of correct products, such as cyclic products or a specific product of interest, while minimizing the output of byproducts. This approach results in nonlinear fitness landscapes that are considerably more epistatic. Overall, our experimental and computational framework provides focused characterization of evolutionary emergence of novel enzymatic functions in the context of microevolutionary exploration of sequence space around naturally occurring enzymes.

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.

Effects of Single Mutations on Protein Stability Are Gaussian Distributed

The distribution of protein stability effects is known to be well approximated by a Gaussian distribution from previous empirical fits. Starting from first-principles statistical mechanics, we more rigorously motivate this empirical observation by deriving per-residue-position protein stability effects to be Gaussian. Our derivation requires the number of amino acids to be large, which is satisfied by the standard set of 20 amino acids found in nature. No assumption is needed on the number of residues in close proximity in space, in contrast to previous applications of the central limit theorem to protein energetics. We support our derivation results with computational and experimental data on mutant protein stabilities across all types of protein residues.