Deep mutational scanning of Pneumocystis jirovecii dihydrofolate reductase reveals allosteric mechanism of resistance to an antifolate

Pneumocystis jirovecii is a fungal pathogen that causes pneumocystis pneumonia, a disease that mainly affects immunocompromised individuals. This fungus has historically been hard to study because of our inability to grow it in vitro. One of the main drug targets in P. jirovecii is its dihydrofolate reductase (PjDHFR). Here, by using functional complementation of the baker's yeast ortholog, we show that PjDHFR can be inhibited by the antifolate methotrexate in a dose-dependent manner. Using deep mutational scanning of PjDHFR, we identify mutations conferring resistance to methotrexate. Thirty-one sites spanning the protein have at least one mutation that leads to resistance, for a total of 355 high-confidence resistance mutations. Most resistance-inducing mutations are found inside the active site, and many are structurally equivalent to mutations known to lead to resistance to different antifolates in other organisms. Some sites show specific resistance mutations, where only a single substitution confers resistance, whereas others are more permissive, as several substitutions at these sites confer resistance. Surprisingly, one of the permissive sites (F199) is without direct contact to either ligand or cofactor, suggesting that it acts through an allosteric mechanism. Modeling changes in binding energy between F199 mutants and drug shows that most mutations destabilize interactions between the protein and the drug. This evidence points towards a more important role of this position in resistance than previously estimated and highlights potential unknown allosteric mechanisms of resistance to antifolate in DHFRs. Our results offer unprecedented resources for the interpretation of mutation effects in the main drug target of an uncultivable fungal pathogen.

Compensatory mutations potentiate constructive neutral evolution by gene duplication

Protein functions generally depend on their assembly into complexes. During evolution, some complexes have transitioned from homomers encoded by a single gene to heteromers encoded by duplicate genes. This transition could occur without adaptive evolution through intermolecular compensatory mutations. Here, we experimentally duplicate and evolve an homodimeric enzyme to examine if and how this could happen. We identify hundreds of deleterious mutations that inactivate individual homodimers but produce functional enzymes when co-expressed as duplicated proteins that heterodimerize. The structure of one such heteromer reveals how both losses of function are buffered through the introduction of asymmetry in the complex that allows them to subfunctionalize. Constructive neutral evolution can thus occur by gene duplication followed by only one deleterious mutation per duplicate.

Epistasis between promoter activity and coding mutations shapes gene evolvability

The evolution of protein-coding genes proceeds as mutations act on two main dimensions: regulation of transcription level and the coding sequence. The extent and impact of the connection between these two dimensions are largely unknown because they have generally been studied independently. By measuring the fitness effects of all possible mutations on a protein complex at various levels of promoter activity, we show that promoter activity at the optimal level for the wild-type protein masks the effects of both deleterious and beneficial coding mutations. Mutations that are deleterious at low activity but masked at optimal activity are slightly destabilizing for individual subunits and binding interfaces. Coding mutations that increase protein abundance are beneficial at low expression but could potentially incur a cost at high promoter activity. We thereby demonstrate that promoter activity in interaction with protein properties can dictate which coding mutations are beneficial, neutral, or deleterious.

Asymmetrical dose responses shape the evolutionary trade-off between antifungal resistance and nutrient use

Antimicrobial resistance is an emerging threat for public health. The success of resistance mutations depends on the trade-off between the benefits and costs they incur. This trade-off is largely unknown and uncharacterized for antifungals. Here, we systematically measure the effect of all amino acid substitutions in the yeast cytosine deaminase Fcy1, the target of the antifungal 5-fluorocytosine (5-FC, flucytosine). We identify over 900 missense mutations granting resistance to 5-FC, a large fraction of which appear to act through destabilization of the protein. The relationship between 5-FC resistance and growth sustained by cytosine deamination is characterized by a sharp trade-off, such that small gains in resistance universally lead to large losses in canonical enzyme function. We show that this steep relationship can be explained by differences in the dose-response functions of 5-FC and cytosine. Finally, we observe the same trade-off shape for the orthologue of FCY1 in Cryptoccocus neoformans, a human pathogen. Our results provide a powerful resource and platform for interpreting drug target variants in fungal pathogens as well as unprecedented insights into resistance-function trade-offs.

Deep Mutational Scanning of Protein-Protein Interactions Between Partners Expressed from Their Endogenous Loci In Vivo

Deep mutational scanning (DMS) generates mutants of a protein of interest in a comprehensive manner. CRISPR-Cas9 technology enables large-scale genome editing with high efficiency. Using both DMS and CRISPR-Cas9 therefore allows us to investigate the effects of thousands of mutations inserted directly in the genome. Combined with protein-fragment complementation assay (PCA), which enables the quantitative measurement of protein-protein interactions (PPIs) in vivo, these methods allow for the systematic assessment of the effects of mutations on PPIs in living cells. Here, we describe a method leveraging DMS, CRISPR-Cas9, and PCA to study the effect of point mutations on PPIs mediated by protein domains in yeast.

High-Throughput Gene Mutagenesis Screening Using Base Editing

Base editing is a CRISPR-Cas9 genome engineering tool that allows programmable mutagenesis without the creation of double-stranded breaks. Here, we describe the design and execution of large-scale base editing screens using the Target-AID base editor in yeast. Using this approach, thousands of sites can be mutated simultaneously. The effects of these mutations on fitness can be measured using a pooled growth competition assay followed by DNA sequencing of gRNAs as barcodes.

Perturbing proteomes at single residue resolution using base editing

Base editors derived from CRISPR-Cas9 systems and DNA editing enzymes offer an unprecedented opportunity for the precise modification of genes, but have yet to be used at a genome-scale throughput. Here, we test the ability of the Target-AID base editor to systematically modify genes genome-wide by targeting yeast essential genes. We mutate around 17,000 individual sites in parallel across more than 1500 genes. We identify over 700 sites at which mutations have a significant impact on fitness. Using previously determined and preferred Target-AID mutational outcomes, we find that gRNAs with significant effects on fitness are enriched in variants predicted to be deleterious based on residue conservation and predicted protein destabilization. We identify key features influencing effective gRNAs in the context of base editing. Our results show that base editing is a powerful tool to identify key amino acid residues at the scale of proteomes.

Extended Linkers Improve the Detection of Protein-protein Interactions (PPIs) by Dihydrofolate Reductase Protein-fragment Complementation Assay (DHFR PCA) in Living Cells

Understanding the function of cellular systems requires describing how proteins assemble with each other into transient and stable complexes and to determine their spatial relationships. Among the tools available to perform these analyses on a large scale is Protein-fragment Complementation Assay based on the dihydrofolate reductase (DHFR PCA). Here we test how longer linkers between the fusion proteins and the reporter fragments affect the performance of this assay. We investigate the architecture of the RNA polymerases, the proteasome and the conserved oligomeric Golgi (COG) complexes in living cells and performed large-scale screens with these extended linkers. We show that longer linkers significantly improve the detection of protein-protein interactions and allow to measure interactions further in space than the standard ones. We identify new interactions, for instance between the retromer complex and proteins related to autophagy and endocytosis. Longer linkers thus contribute an enhanced additional tool to the existing toolsets for the detection and measurements of protein-protein interactions and protein proximity in living cells.