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.

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 (2) = 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 (2) = 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.

Interactions among mutations within a protein have the potential to make molecular evolution contingent and irreversible, but the extent to which epistasis actually shaped historical evolutionary trajectories is unclear. We addressed this question by identifying all amino acid substitutions that occurred during the billion-year evolutionary history of the heat shock protein 90 (Hsp90) ATPase domain beginning from a deep eukaryotic ancestor to modern Saccharomyces cerevisiae and then precisely measuring their fitness effects when introduced into both extant and reconstructed ancestral Hsp90 proteins. We find a pervasive influence of epistasis: of 98 derived states that evolved during history, most were deleterious at times before they happened, and the vast majority also became subsequently entrenched, with the ancestral state becoming deleterious after its substitution. This epistasis was primarily caused by specific interactions among sites rather than a general permissive or restrictive effect on the protein's tolerance to mutation. Our results show that epistasis continually opens and closes windows of mutational opportunity over evolutionary timescales, producing histories and biological states that reflect the transient internal constraints imposed by a protein's fleeting sequence states.

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.

High-throughput sequencing has enabled many powerful approaches in biological research. Here, we review sequencing approaches to measure frequency changes within engineered mutational libraries subject to selection. These analyses can provide direct estimates of biochemical and fitness effects for all individual mutations across entire genes (and likely compact genomes in the near future) in genetically tractable systems such as microbes, viruses, and mammalian cells. The effects of mutations on experimental fitness can be assessed using sequencing to monitor time-dependent changes in mutant frequency during bulk competitions. The impact of mutations on biochemical functions can be determined using reporters or other means of separating variants based on individual activities (e.g., binding affinity for a partner molecule can be interrogated using surface display of libraries of mutant proteins and isolation of bound and unbound populations). The comprehensive investigation of mutant effects on both biochemical function and experimental fitness provide promising new avenues to investigate the connections between biochemistry, cell physiology, and evolution. We summarize recent findings from systematic mutational analyses; describe how they relate to a field rich in both theory and experimentation; and highlight how they may contribute to ongoing and future research into protein structure-function relationships, systems-level descriptions of cell physiology, and population-genetic inferences on the relative contributions of selection and drift.