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

  1. Julia M Flynn  Is a corresponding author
  2. Neha Samant
  3. Gily Schneider-Nachum
  4. David T Barkan
  5. Nese Kurt Yilmaz
  6. Celia A Schiffer
  7. Stephanie A Moquin
  8. Dustin Dovala
  9. Daniel NA Bolon  Is a corresponding author
  1. Department of Biochemistry and Molecular Biotechnology, University of Massachusetts Chan Medical School, United States
  2. Novartis Institutes for Biomedical Research, United States
7 figures, 1 table and 2 additional files

Figures

Figure 1 with 1 supplement
Experimental strategy to measure the function of all individual mutations of main protease (Mpro).

(A) Fluorescence resonance energy transfer (FRET)-based reporter screen. The Mpro variants were sorted based on their ability to cleave at the Mpro cut-site, separating the YFP-CFP FRET pair. Cells …

Figure 1—figure supplement 1
Main protease (Mpro) expression in cells harboring the LexA-UbMpro plasmid construct.

(A) Yeast cells transformed with a plasmid expressing C145A Ub-Mpro-his6 under the LexA promoter were grown to exponential phase followed by the addition of 2 µM β-estradiol to induce expression for …

Figure 2 with 1 supplement
Main protease (Mpro) functional scores are reproducible, and variants can be clearly distinguished based on function.

(A) Correlation between biological replicates of functional scores of all Mpro variants for each screen. Red line indicates best fit. (B) Distribution of functional scores for all variants (gray), …

Figure 2—source data 1

Sequencing counts and functional scores for each amino acid of main protease (Mpro) in both replicates of all three screens.

For each dataset, the sequencing counts, unnormalized functional scores, and normalized functional scores (normalized to average stop = 0, average wild-type barcode = 1) are reported. For the growth screens, the selection coefficients are also reported. All figures in this paper use the data from replicate 1 of each screen.

https://cdn.elifesciences.org/articles/77433/elife-77433-fig2-data1-v2.xlsx
Figure 2—figure supplement 1
Cumulative frequency distributions for all variants (gray), stops (red), and wild-type (WT) barcodes (blue) for all three screens.
Figure 3 with 2 supplements
Heatmap representation of the main protease (Mpro) functional scores measured in the fluorescence resonance energy transfer (FRET) screen (replicate 1).

Arrows represent positions that form β-sheets, coils represent α-helices, and red triangles indicate the catalytic dyad residues H41 and C145.

Figure 3—figure supplement 1
Heatmap representation of scores from the transcription factor (TF) screen (replicate 1).

Arrows represent positions that form beta sheets, coils represent α-helices, and red triangles indicate the catalytic dyad residues H41 and C145.

Figure 3—figure supplement 2
Heatmap representation of scores from the growth screen (replicate 1).

Arrows represent positions that form beta sheets, coils represent α-helices, and red triangles indicate the catalytic dyad residues H41 and C145.

Functional scores reflect fundamental biophysical constraints of main protease (Mpro).

(A) Heatmap representation of the average functional score at each position (excluding stops) in replicate 1 of each screen (see Figure 4—source data 1). (B) The average functional score at each …

Figure 4—source data 1

Average functional score (excluding stops) at each position of main protease (Mpro) in replicate 1 of each screen.

https://cdn.elifesciences.org/articles/77433/elife-77433-fig4-data1-v2.xlsx
Figure 4—source data 2

Comparison of previously measured relative catalytic rates of individual mutations to functional scores.

https://cdn.elifesciences.org/articles/77433/elife-77433-fig4-data2-v2.docx
Figure 5 with 1 supplement
Functional scores indicate that natural amino acid variants of main protease (Mpro) are generally fit.

(A) Comparison of functional scores in the FRET screen (left panel) and growth screen (right panel) to the number of observations among clinical samples. All missense mutations excluding stops are …

Figure 5—source data 1

Frequency at which the clinical variants of the main protease (Mpro) gene have been observed.

https://cdn.elifesciences.org/articles/77433/elife-77433-fig5-data1-v2.xlsx
Figure 5—figure supplement 1
Functional scores indicate that natural amino acid variants of main protease (Mpro) are generally fit.

(A) Comparison of functional scores in the transcription factor (TF) screen to the number of observations among clinical samples. All missense mutations excluding stops are indicated with black …

Figure 6 with 1 supplement
Structural distribution of main protease (Mpro) positions that are intolerant to mutation.

(A) Mpro positions that are intolerant of mutations with 17 or more substitutions having null-like function are represented by red spheres on chain A (shown in gray) and pink spheres on chain B …

Figure 6—figure supplement 1
Comparison of the average transcription factor (TF) functional score of each position to conservation observed in a broad sample of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) main protease (Mpro) homologs.
Figure 7 with 1 supplement
Substrate and inhibitor binding sites are variably sensitive to mutation.

(A) All main protease (Mpro) positions that contact the Nsp4/5 substrate peptide are represented in spheres and colored by their average fluorescence resonance energy transfer (FRET) functional …

Figure 7—figure supplement 1
Substrate and inhibitor binding sites are variably sensitive to mutation.

(A) All main protease (Mpro) positions that contact the Nsp4/5 substrate peptide are represented in spheres and colored by their average transcription factor (TF) functional score (PDB 7T70). The …

Tables

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Gene
(SARS-CoV-2)
ORF1ab/ nsp5A-BNIH GenBankNC_045512Mpro
Strain, Strain background (Saccharomyces cerevisiae)W303Saccharomyces Genome DatabaseGenBank JRIU00000000
Antibodyanti-his tag HRP-labelled (Mouse monoclonal)R&D systemsCAT#: MAB050HWB (1:4000)
Recombinant DNA reagentBarcoded UbMpro plasmid libraryThis paperp416LexA-
UbMpro(lib)-N18
See Materials and Methods section “Generating mutant libraries”
Recombinant DNA reagentBarcoded WT UbMpro plasmidThis paperp416LexA-
UbMpro(WT)-N18
See Materials and Methods section “Construction of WT Ub-Mpro vector”
Recombinant DNA reagentC145A-Mpro-his6 plasmidThis paperp416LexA-UbMpro
(C145A)-his
See Materials and Methods section “Analysis of Mpro expression”
Recombinant DNA reagentpCyPet-HisAddgene#14,040
Recombinant DNA reagentpYPet-HisAddgene#14,031
Recombinant DNA reagentCyPet-MproCS-YPet fusion geneThis paperSee Materials and Methods section “Generating FRET strain”
Recombinant DNA reagentpDK-ATCPMID:28660202Integrative bidirectional plasmid with TEF and CUP promoters
Recombinant DNA reagentpDK-ATGPMID:28660202Integrative bidirectional plasmid with TEF and GPD promoters
Recombinant DNA reagentDBD-MproCS-AD fusion geneThis paperSee Materials and Methods section “Generating split TF strain”
Commercial assay or kitKAPA SYBR FAST qPCR Master MixKapa BiosystemsKK4600
Commercial assay or kitBCA protein assay kitPierceCAT#
23,225
Chemical compound, drugβ-EstradiolSigma AldrichE2768
Software, algorithmScripts to tabulate variant countsThis paperhttps://github.com/JuliaFlynn/BolonLab,
(copy archived at swh:1:rev:b54d80818c2681fb89533ae330c18a3d39f32ab6)
See Materials and Methods section “Analysis of Illumina sequencing data”
Software, algorithmScripts to associate barcodes with variantsThis paperhttps://github.com/JuliaFlynn/PacBio_barcode_assocation,
(copy archived at swh:1:rev:29eac92475a9ff8e24fb390986c865b504c03f51)
See Materials and Methods section “Barcode Association”
Software, algorithmGraphPad Prism 9Graphpad.comRRID:SCR_008520
Software, algorithmFlowjo v.10.8.0BD BiosciencesRRID:SCR_008520
Software, algorithmPymol v. 2.5.2SchrödingerRRID:SCR_000305
Software, algorithmMatPlotLibhttp://matplotlib.sourceforge.netRRID:SCR_008624
Sequence-based reagentSequencing primersThis paperSee Supplementary file 1
Sequence-based reagentSite-directed mutagenesis primersThis paperSee Supplementary file 1

Additional files

Download links