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

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

Contributions of Hyperactive Mutations in Mpro from SARS-CoV-2 to Drug Resistance

The appearance and spread of mutations that cause drug resistance in rapidly evolving diseases, including infections by the SARS-CoV-2 virus, are major concerns for human health. Many drugs target enzymes, and resistance-conferring mutations impact inhibitor binding or enzyme activity. Nirmatrelvir, the most widely used inhibitor currently used to treat SARS-CoV-2 infections, targets the main protease (Mpro) preventing it from processing the viral polyprotein into active subunits. Our previous work systematically analyzed resistance mutations in Mpro that reduce binding to inhibitors; here, we investigate mutations that affect enzyme function. Hyperactive mutations that increase Mpro activity can contribute to drug resistance but have not been thoroughly studied. To explore how hyperactive mutations contribute to resistance, we comprehensively assessed how all possible individual mutations in Mpro affect enzyme function using a mutational scanning approach with a fluorescence resonance energy transfer (FRET)-based yeast readout. We identified hundreds of mutations that significantly increased the Mpro activity. Hyperactive mutations occurred both proximal and distal to the active site, consistent with protein stability and/or dynamics impacting activity. Hyperactive mutations were observed 3 times more than mutations which reduced apparent binding to nirmatrelvir in recent studies of laboratory-grown viruses selected for drug resistance. Hyperactive mutations were also about three times more prevalent than nirmatrelvir binding mutations in sequenced isolates from circulating SARS-CoV-2. Our findings indicate that hyperactive mutations are likely to contribute to the natural evolution of drug resistance in Mpro and provide a comprehensive list for future surveillance efforts.

Systematic profiling of dominant ubiquitin variants reveals key functional nodes contributing to evolutionary selection

Dominant-negative mutations can help to investigate the biological mechanisms and to understand the selective pressures for multifunctional proteins. However, most studies have focused on recessive mutant effects that occur in the absence of a second functional gene copy, which overlooks the fact that most eukaryotic genomes contain more than one copy of many genes. We have identified dominant effects on yeast growth rate among all possible point mutations in ubiquitin expressed alongside a wild-type allele. Our results reveal more than 400 dominant-negative mutations, indicating that dominant-negative effects make a sizable contribution to selection acting on ubiquitin. Cellular and biochemical analyses of individual ubiquitin variants show that dominant-negative effects are explained by varied accumulation of polyubiquitinated cellular proteins and/or defects in conjugation of ubiquitin variants to ubiquitin ligases. Our approach to identify dominant-negative mutations is general and can be applied to other proteins of interest.

Systematic Analyses of the Resistance Potential of Drugs Targeting SARS-CoV-2 Main Protease

Drugs that target the main protease (Mpro) of SARS-CoV-2 are effective therapeutics that have entered clinical use. Wide-scale use of these drugs will apply selection pressure for the evolution of resistance mutations. To understand resistance potential in Mpro, we performed comprehensive surveys of amino acid changes that can cause resistance to nirmatrelvir (Pfizer), and ensitrelvir (Xocova) in a yeast screen. We identified 142 resistance mutations for nirmatrelvir and 177 for ensitrelvir, many of which have not been previously reported. Ninety-nine mutations caused apparent resistance to both inhibitors, suggesting likelihood for the evolution of cross-resistance. The mutation with the strongest drug resistance score against nirmatrelvir in our study (E166V) was the most impactful resistance mutation recently reported in multiple viral passaging studies. Many mutations that exhibited inhibitor-specific resistance were consistent with the distinct interactions of each inhibitor in the substrate binding site. In addition, mutants with strong drug resistance scores tended to have reduced function. Our results indicate that strong pressure from nirmatrelvir or ensitrelvir will select for multiple distinct-resistant lineages that will include both primary resistance mutations that weaken interactions with drug while decreasing enzyme function and compensatory mutations that increase enzyme activity. The comprehensive identification of resistance mutations enables the design of inhibitors with reduced potential of developing resistance and aids in the surveillance of drug resistance in circulating viral populations.

A complete allosteric map of a GTPase switch in its native cellular network

Allosteric regulation is central to protein function in cellular networks. A fundamental open question is whether cellular regulation of allosteric proteins occurs only at a few defined positions or at many sites distributed throughout the structure. Here, we probe the regulation of GTPases-protein switches that control signaling through regulated conformational cycling-at residue-level resolution by deep mutagenesis in the native biological network. For the GTPase Gsp1/Ran, we find that 28% of the 4,315 assayed mutations show pronounced gain-of-function responses. Twenty of the sixty positions enriched for gain-of-function mutations are outside the canonical GTPase active site switch regions. Kinetic analysis shows that these distal sites are allosterically coupled to the active site. We conclude that the GTPase switch mechanism is broadly sensitive to cellular allosteric regulation. Our systematic discovery of new regulatory sites provides a functional map to interrogate and target GTPases controlling many essential biological processes.

Mutational fitness landscape and drug resistance

Robust technology has been developed to systematically quantify fitness landscapes that provide valuable opportunities to improve our understanding of drug resistance and define new avenues to develop drugs with reduced resistance susceptibility. We outline the critical importance of drug resistance studies and the potential for fitness landscape approaches to contribute to this effort. We describe the major technical advancements in mutational scanning, which is the primary approach used to quantify protein fitness landscapes. There are many complex steps to consider in planning and executing mutational scanning projects including developing a selection scheme, generating mutant libraries, tracking the frequency of variants using next-generation sequencing, and processing and interpreting the data. Key experimental parameters impacting each of these steps are discussed to aid in planning fitness landscape studies. There is a strong need for improved understanding of drug resistance, and fitness landscapes provide a promising new approach.

Sequence dependencies and biophysical features both govern cleavage of diverse cut-sites by HIV protease

The infectivity of HIV-1 requires its protease (PR) cleave multiple cut-sites with low sequence similarity. The diversity of cleavage sites has made it challenging to investigate the underlying sequence properties that determine binding and turnover of substrates by PR. We engineered a mutational scanning approach utilizing yeast display, flow cytometry, and deep sequencing to systematically measure the impacts of all individual amino acid changes at 12 positions in three different cut-sites (MA/CA, NC/p1, and p1/p6). The resulting fitness landscapes revealed common physical features that underlie cutting of all three cut-sites at the amino acid positions closest to the scissile bond. In contrast, positions more than two amino acids away from the scissile bond exhibited a strong dependence on the sequence background of the rest of the cut-site. We observed multiple amino acid changes in cut-sites that led to faster cleavage rates, including a preference for negative charge five and six amino acids away from the scissile bond at locations where the surface of protease is positively charged. Analysis of individual cut sites using full-length matrix-capsid proteins indicate that long-distance sequence context can contribute to cutting efficiency such that analyses of peptides or shorter engineered constructs including those in this work should be considered carefully. This work provides a framework for understanding how diverse substrates interact with HIV-1 PR and can be extended to investigate other viral PRs with similar properties.

Defining the substrate envelope of SARS-CoV-2 main protease to predict and avoid drug resistance

Coronaviruses can evolve and spread rapidly to cause severe disease morbidity and mortality, as exemplified by SARS-CoV-2 variants of the COVID-19 pandemic. Although currently available vaccines remain mostly effective against SARS-CoV-2 variants, additional treatment strategies are needed. Inhibitors that target essential viral enzymes, such as proteases and polymerases, represent key classes of antivirals. However, clinical use of antiviral therapies inevitably leads to emergence of drug resistance. In this study we implemented a strategy to pre-emptively address drug resistance to protease inhibitors targeting the main protease (M^pro) of SARS-CoV-2, an essential enzyme that promotes viral maturation. We solved nine high-resolution cocrystal structures of SARS-CoV-2 M^pro bound to substrate peptides and six structures with cleavage products. These structures enabled us to define the substrate envelope of M^pro, map the critical recognition elements, and identify evolutionarily vulnerable sites that may be susceptible to resistance mutations that would compromise binding of the newly developed M^pro inhibitors. Our results suggest strategies for developing robust inhibitors against SARS-CoV-2 that will retain longer-lasting efficacy against this evolving viral pathogen.

Comprehensive fitness landscape of SARS-CoV-2 Mpro reveals insights into viral resistance mechanisms

With the continual evolution of new strains of SARS-CoV-2 that are more virulent, transmissible, and able to evade current vaccines, there is an urgent need for effective anti-viral drugs SARS-CoV-2 main protease (M^pro) is a leading target for drug design due to its conserved and indispensable role in the viral life cycle. Drugs targeting M^pro appear promising but will elicit selection pressure for resistance. To understand resistance potential in M^pro, we performed a comprehensive mutational scan of the protease that analyzed the function of all possible single amino acid changes. We developed three separate high-throughput assays of M^pro function in yeast, based on either the ability of M^pro variants to cleave at a defined cut-site or on the toxicity of their expression to yeast. We used deep sequencing to quantify the functional effects of each variant in each screen. The protein fitness landscapes from all three screens were strongly correlated, indicating that they captured the biophysical properties critical to M^pro function. The fitness landscapes revealed a non-active site location on the surface that is extremely sensitive to mutation making it a favorable location to target with inhibitors. In addition, we found a network of critical amino acids that physically bridge the two active sites of the M^pro dimer. The clinical variants of M^pro were predominantly functional in our screens, indicating that M^pro is under strong selection pressure in the human population. Our results provide predictions of mutations that will be readily accessible to M^pro evolution and that are likely to contribute to drug resistance. This complete mutational guide of M^pro can be used in the design of inhibitors with reduced potential of evolving viral resistance.

Quantitative structural analysis of influenza virus by cryo-electron tomography and convolutional neural networks

Influenza viruses pose severe public health threats globally. Influenza viruses are extensively pleomorphic, in shape, size, and organization of viral proteins. Analysis of influenza morphology and ultrastructure can help elucidate viral structure-function relationships and aid in therapeutics and vaccine development. While cryo-electron tomography (cryoET) can depict the 3D organization of pleomorphic influenza, the low signal-to-noise ratio inherent to cryoET and viral heterogeneity have precluded detailed characterization of influenza viruses. In this report, we leveraged convolutional neural networks and cryoET to characterize the morphological architecture of the A/Puerto Rico/8/34 (H1N1) influenza strain. Our pipeline improved the throughput of cryoET analysis and accurately identified viral components within tomograms. Using this approach, we successfully characterized influenza morphology, glycoprotein density, and conducted subtomogram averaging of influenza glycoproteins. Application of this processing pipeline can aid in the structural characterization of not only influenza viruses, but other pleomorphic viruses and infected cells.

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

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

The Adaptive Potential of the Middle Domain of Yeast Hsp90

The distribution of fitness effects (DFEs) of new mutations across different environments quantifies the potential for adaptation in a given environment and its cost in others. So far, results regarding the cost of adaptation across environments have been mixed, and most studies have sampled random mutations across different genes. Here, we quantify systematically how costs of adaptation vary along a large stretch of protein sequence by studying the distribution of fitness effects of the same ≈2,300 amino-acid changing mutations obtained from deep mutational scanning of 119 amino acids in the middle domain of the heat shock protein Hsp90 in five environments. This region is known to be important for client binding, stabilization of the Hsp90 dimer, stabilization of the N-terminal-Middle and Middle-C-terminal interdomains, and regulation of ATPase–chaperone activity. Interestingly, we find that fitness correlates well across diverse stressful environments, with the exception of one environment, diamide. Consistent with this result, we find little cost of adaptation; on average only one in seven beneficial mutations is deleterious in another environment. We identify a hotspot of beneficial mutations in a region of the protein that is located within an allosteric center. The identified protein regions that are enriched in beneficial, deleterious, and costly mutations coincide with residues that are involved in the stabilization of Hsp90 interdomains and stabilization of client-binding interfaces, or residues that are involved in ATPase–chaperone activity of Hsp90. Thus, our study yields information regarding the role and adaptive potential of a protein sequence that complements and extends known structural information.

Constrained Mutational Sampling of Amino Acids in HIV-1 Protease Evolution

The evolution of HIV-1 protein sequences should be governed by a combination of factors including nucleotide mutational probabilities, the genetic code, and fitness. The impact of these factors on protein sequence evolution is interdependent, making it challenging to infer the individual contribution of each factor from phylogenetic analyses alone. We investigated the protein sequence evolution of HIV-1 by determining an experimental fitness landscape of all individual amino acid changes in protease. We compared our experimental results to the frequency of protease variants in a publicly available data set of 32,163 sequenced isolates from drug-naïve individuals. The most common amino acids in sequenced isolates supported robust experimental fitness, indicating that the experimental fitness landscape captured key features of selection acting on protease during viral infections of hosts. Amino acid changes requiring multiple mutations from the likely ancestor were slightly less likely to support robust experimental fitness than single mutations, consistent with the genetic code favoring chemically conservative amino acid changes. Amino acids that were common in sequenced isolates were predominantly accessible by single mutations from the likely protease ancestor. Multiple mutations commonly observed in isolates were accessible by mutational walks with highly fit single mutation intermediates. Our results indicate that the prevalence of multiple-base mutations in HIV-1 protease is strongly influenced by mutational sampling.

Synonymous Mutations at the Beginning of the Influenza A Virus Hemagglutinin Gene Impact Experimental Fitness

The fitness effects of synonymous mutations can provide insights into biological and evolutionary mechanisms. We analyzed the experimental fitness effects of all single-nucleotide mutations, including synonymous substitutions, at the beginning of the influenza A virus hemagglutinin (HA) gene. Many synonymous substitutions were deleterious both in bulk competition and for individually isolated clones. Investigating protein and RNA levels of a subset of individually expressed HA variants revealed that multiple biochemical properties contribute to the observed experimental fitness effects. Our results indicate that a structural element in the HA segment viral RNA may influence fitness. Examination of naturally evolved sequences in human hosts indicates a preference for the unfolded state of this structural element compared to that found in swine hosts. Our overall results reveal that synonymous mutations may have greater fitness consequences than indicated by simple models of sequence conservation, and we discuss the implications of this finding for commonly used evolutionary tests and analyses.

An experimental evaluation of drug-induced mutational meltdown as an antiviral treatment strategy

The rapid evolution of drug resistance remains a critical public health concern. The treatment of influenza A virus (IAV) has proven particularly challenging, due to the ability of the virus to develop resistance against current antivirals and vaccines. Here, we evaluate a novel antiviral drug therapy, favipiravir, for which the mechanism of action in IAV involves an interaction with the viral RNA-dependent RNA polymerase resulting in an effective increase in the viral mutation rate. We used an experimental evolution framework, combined with novel population genetic method development for inference from time-sampled data, to evaluate the effectiveness of favipiravir against IAV. Evaluating whole genome polymorphism data across 15 time points under multiple drug concentrations and in controls, we present the first evidence for the ability of IAV populations to effectively adapt to low concentrations of favipiravir. In contrast, under high concentrations, we observe population extinction, indicative of mutational meltdown. We discuss the observed dynamics with respect to the evolutionary forces at play and emphasize the utility of evolutionary theory to inform drug development.

Systematic Mutant Analyses Elucidate General and Client-Specific Aspects of Hsp90 Function

To probe the mechanism of the Hsp90 chaperone that is required for the maturation of many signaling proteins in eukaryotes, we analyzed the effects of all individual amino acid changes in the ATPase domain on yeast growth rate. The sensitivity of a position to mutation was strongly influenced by proximity to the phosphates of ATP, indicating that ATPase-driven conformational changes impose stringent physical constraints on Hsp90. To investigate how these constraints may vary for different clients, we performed biochemical analyses on a panel of Hsp90 mutants spanning the full range of observed fitness effects. We observed distinct effects of nine Hsp90 mutations on activation of v-src and glucocorticoid receptor (GR), indicating that different chaperone mechanisms can be utilized for these clients. These results provide a detailed guide for understanding Hsp90 mechanism and highlight the potential for inhibitors of Hsp90 that target a subset of clients.

Quantifying and understanding the fitness effects of protein mutations: Laboratory versus nature

The last decade has seen a growing number of experiments aimed at systematically mapping the effects of mutations in different proteins, and of attempting to correlate their biophysical and biochemical effects with organismal fitness. While insightful, systematic laboratory measurements of fitness effects present challenges and difficulties. Here, we discuss the limitations associated with such measurements, and in particular the challenge of correlating the effects of mutations at the single protein level ("protein fitness") with their effects on organismal fitness. A variety of experimental setups are used, with some measuring the direct effects on protein function and others monitoring the growth rate of a model organism carrying the protein mutants. The manners by which fitness effects are calculated and presented also vary, and the conclusions, including the derived distributions of fitness effects of mutations, vary accordingly. The comparison of the effects of mutations in the laboratory to the natural protein diversity, namely to amino acid changes that have fixed in the course of millions of years of evolution, is also debatable. The results of laboratory experiments may, therefore, be less relevant to understanding long-term inter-species variations yet insightful with regard to short-term polymorphism, for example, in the study of the effects of human SNPs.

Mechanistic Asymmetry in Hsp90 Dimers

Hsp90 is a molecular chaperone that facilitates the maturation of signaling proteins including many kinases and steroid hormone receptors. Through these client proteins, Hsp90 is a key mediator of many physiological processes and has emerged as a promising drug target in cancer. Additionally, Hsp90 can mask or potentiate the impact of mutations in clients with remarkable influence on evolutionary adaptations. The influential roles of Hsp90 in biology and disease have stimulated extensive research into the molecular mechanism of this chaperone. These studies have shown that Hsp90 is a homodimeric protein that requires ATP hydrolysis and a host of accessory proteins termed co-chaperones to facilitate the maturation of clients to their active states. Flexible hinge regions between its three structured domains enable Hsp90 to sample dramatically distinct conformations that are influenced by nucleotide, client, and co-chaperone binding. While it is clear that Hsp90 can exist in symmetrical conformations, recent studies have indicated that this homodimeric chaperone can also assume a variety of asymmetric conformations and complexes that are important for client maturation. The visualization of Hsp90-client complexes at high resolution together with tools to independently manipulate each subunit in the Hsp90 dimer are providing new insights into the asymmetric function of each subunit during client maturation.

A systematic survey of an intragenic epistatic landscape

Mutations are the source of evolutionary variation. The interactions of multiple mutations can have important effects on fitness and evolutionary trajectories. We have recently described the distribution of fitness effects of all single mutations for a nine-amino-acid region of yeast Hsp90 (Hsp82) implicated in substrate binding. Here, we report and discuss the distribution of intragenic epistatic effects within this region in seven Hsp90 point mutant backgrounds of neutral to slightly deleterious effect, resulting in an analysis of more than 1,000 double mutants. We find negative epistasis between substitutions to be common, and positive epistasis to be rare--resulting in a pattern that indicates a drastic change in the distribution of fitness effects one step away from the wild type. This can be well explained by a concave relationship between phenotype and genotype (i.e., a concave shape of the local fitness landscape), suggesting mutational robustness intrinsic to the local sequence space. Structural analyses indicate that, in this region, epistatic effects are most pronounced when a solvent-inaccessible position is involved in the interaction. In contrast, all 18 observations of positive epistasis involved at least one mutation at a solvent-exposed position. By combining the analysis of evolutionary and biophysical properties of an epistatic landscape, these results contribute to a more detailed understanding of the complexity of protein evolution.

Designed Hsp90 heterodimers reveal an asymmetric ATPase-driven mechanism in vivo

Hsp90 is a homodimeric ATPase that is essential in eukaryotes for the maturation of client proteins frequently involved in signal transduction, including many kinases and nuclear steroid hormone receptors. Competitive inhibitors of ATP binding to Hsp90 prevent client maturation and show promise as anticancer agents in clinical trials. However, the role of ATP binding and hydrolysis in each subunit of the Hsp90 dimer has been difficult to investigate because of an inability to assemble and study dimers of defined composition. We used protein engineering to generate functional Hsp90 subunits that preferentially assemble as heterodimers. We analyzed dimers wherein one subunit harbors a disruptive mutation and observed that ATP binding by both subunits is essential for function in yeast, whereas ATP hydrolysis is only required in one subunit. These findings demonstrate important functional contributions from both symmetric and asymmetric Hsp90 dimers and provide valuable reagents for future investigations of Hsp90 mechanism.

Alanine scan of core positions in ubiquitin reveals links between dynamics, stability, and function

Mutations at solvent-inaccessible core positions in proteins can impact function through many biophysical mechanisms including alterations to thermodynamic stability and protein dynamics. As these properties of proteins are difficult to investigate, the impacts of core mutations on protein function are poorly understood for most systems. Here, we determined the effects of alanine mutations at all 15 core positions in ubiquitin on function in yeast. The majority (13 of 15) of alanine substitutions supported yeast growth as the sole ubiquitin. Both the two null mutants (I30A and L43A) were less stable to temperature-induced unfolding in vitro than wild type (WT) but were well folded at physiological temperatures. Heteronuclear NMR studies indicated that the L43A mutation reduces temperature stability while retaining a ground-state structure similar to WT. This structure enables L43A to bind to common ubiquitin receptors in vitro. Many of the core alanine ubiquitin mutants, including one of the null variants (I30A), exhibited an increased accumulation of high-molecular-weight species, suggesting that these mutants caused a defect in the processing of ubiquitin-substrate conjugates. In contrast, L43A exhibited a unique accumulation pattern with reduced levels of high-molecular-weight species and undetectable levels of free ubiquitin. When conjugation to other proteins was blocked, L43A ubiquitin accumulated as free ubiquitin in yeast. Based on these findings, we speculate that ubiquitin's stability to unfolding may be required for efficient recycling during proteasome-mediated substrate degradation.

Shifting fitness landscapes in response to altered environments

The role of adaptation in molecular evolution has been contentious for decades. Here, we shed light on the adaptive potential in Saccharomyces cerevisiae by presenting systematic fitness measurements for all possible point mutations in a region of Hsp90 under four environmental conditions. Under elevated salinity, we observe numerous beneficial mutations with growth advantages up to 7% relative to the wild type. All of these beneficial mutations were observed to be associated with high costs of adaptation. We thus demonstrate that an essential protein can harbor adaptive potential upon an environmental challenge, and report a remarkable fit of the data to a version of Fisher's geometric model that focuses on the fitness trade-offs between mutations in different environments.

Latent effects of Hsp90 mutants revealed at reduced expression levels

In natural systems, selection acts on both protein sequence and expression level, but it is unclear how selection integrates over these two dimensions. We recently developed the EMPIRIC approach to systematically determine the fitness effects of all possible point mutants for important regions of essential genes in yeast. Here, we systematically investigated the fitness effects of point mutations in a putative substrate binding loop of yeast Hsp90 (Hsp82) over a broad range of expression strengths. Negative epistasis between reduced expression strength and amino acid substitutions was common, and the endogenous expression strength frequently obscured mutant defects. By analyzing fitness effects at varied expression strengths, we were able to uncover all mutant effects on function. The majority of mutants caused partial functional defects, consistent with this region of Hsp90 contributing to a mutation sensitive and critical process. These results demonstrate that important functional regions of proteins can tolerate mutational defects without experimentally observable impacts on fitness.

Analyses of the effects of all ubiquitin point mutants on yeast growth rate

The amino acid sequence of a protein governs its function. We used bulk competition and focused deep sequencing to investigate the effects of all ubiquitin point mutants on yeast growth rate. Many aspects of ubiquitin function have been carefully studied, which enabled interpretation of our growth analyses in light of a rich structural, biophysical and biochemical knowledge base. In one highly sensitive cluster on the surface of ubiquitin, almost every amino acid substitution caused growth defects. In contrast, the opposite face tolerated virtually all possible substitutions. Surface locations between these two faces exhibited intermediate mutational tolerance. The sensitive face corresponds to the known interface for many binding partners. Across all surface positions, we observe a strong correlation between burial at structurally characterized interfaces and the number of amino acid substitutions compatible with robust growth. This result indicates that binding is a dominant determinant of ubiquitin function. In the solvent-inaccessible core of ubiquitin, all positions tolerated a limited number of substitutions, with hydrophobic amino acids especially interchangeable. Some mutations null for yeast growth were previously shown to populate folded conformations indicating that, for these mutants, subtle changes to conformation caused functional defects. The most sensitive region to mutation within the core was located near the C-terminus that is a focal binding site for many critical binding partners. These results indicate that core mutations may frequently cause functional defects through subtle disturbances to structure or dynamics.

Hydrophobic core flexibility modulates enzyme activity in HIV-1 protease

Human immunodeficiency virus Type-1 (HIV-1) protease is crucial for viral maturation and infectivity. Studies of protease dynamics suggest that the rearrangement of the hydrophobic core is essential for enzyme activity. Many mutations in the hydrophobic core are also associated with drug resistance and may modulate the core flexibility. To test the role of flexibility in protease activity, pairs of cysteines were introduced at the interfaces of flexible regions remote from the active site. Disulfide bond formation was confirmed by crystal structures and by alkylation of free cysteines and mass spectrometry. Oxidized and reduced crystal structures of these variants show the overall structure of the protease is retained. However, cross-linking the cysteines led to drastic loss in enzyme activity, which was regained upon reducing the disulfide cross-links. Molecular dynamics simulations showed that altered dynamics propagated throughout the enzyme from the engineered disulfide. Thus, altered flexibility within the hydrophobic core can modulate HIV-1 protease activity, supporting the hypothesis that drug resistant mutations distal from the active site can alter the balance between substrate turnover and inhibitor binding by modulating enzyme activity.