Software for the analysis and visualization of deep mutational scanning data

popDMS infers mutation effects from deep mutational scanning data

SUMMARY

Deep mutational scanning (DMS) experiments provide a powerful method to measure the functional effects of genetic mutations at massive scales. However, the data generated from these experiments can be difficult to analyze, with significant variation between experimental replicates. To overcome this challenge, we developed popDMS, a computational method based on population genetics theory, to infer the functional effects of mutations from DMS data. Through extensive tests, we found that the functional effects of single mutations and epistasis inferred by popDMS are highly consistent across replicates, comparing favorably with existing methods. Our approach is flexible and can be widely applied to DMS data that includes multiple time points, multiple replicates, and different experimental conditions.

AVAILABILITY AND IMPLEMENTATION

popDMS is implemented in Python and Julia, and is freely available on GitHub at https://github.com/bartonlab/popDMS.

Mapping mutational fitness effects across the coxsackievirus B3 proteome reveals distinct profiles of mutation tolerability

RNA viruses have notoriously high mutation rates due to error-prone replication by their RNA polymerase. However, natural selection concentrates variability in a few key viral proteins. To test whether this stems from different mutation tolerance profiles among viral proteins, we measured the effect of >40,000 non-synonymous mutations across the full proteome of coxsackievirus B3 as well as >97% of all possible codon deletions in the nonstructural proteins. We find significant variation in mutational tolerance within and between individual viral proteins, which correlated with both general and protein-specific structural and functional attributes. Furthermore, mutational fitness effects remained stable across cell lines, suggesting selection pressures are mostly conserved across environments. In addition to providing a rich dataset for understanding virus biology and evolution, our results illustrate that incorporation of mutational tolerance data into druggable pocket discovery can aid in selecting targets with high barriers to drug resistance.

Natural variation in neuraminidase activity influences the evolutionary potential of the seasonal H1N1 lineage hemagglutinin

The antigenic evolution of the influenza A virus hemagglutinin (HA) gene poses a major challenge for the development of vaccines capable of eliciting long-term protection. Prior efforts to understand the mechanisms that govern viral antigenic evolution mainly focus on HA in isolation, ignoring the fact that HA must act in concert with the viral neuraminidase (NA) during replication and spread. Numerous studies have demonstrated that the degree to which the receptor-binding avidity of HA and receptor-cleaving activity of NA are balanced with each other influences overall viral fitness. We recently showed that changes in NA activity can significantly alter the mutational fitness landscape of HA in the context of a lab-adapted virus strain. Here, we test whether natural variation in relative NA activity can influence the evolutionary potential of HA in the context of the seasonal H1N1 lineage (pdmH1N1) that has circulated in humans since the 2009 pandemic. We observed substantial variation in the relative activities of NA proteins encoded by a panel of H1N1 vaccine strains isolated between 2009 and 2019. We comprehensively assessed the effect of NA background on the HA mutational fitness landscape in the circulating pdmH1N1 lineage using deep mutational scanning and observed pronounced epistasis between NA and residues in or near the receptor-binding site of HA. To determine whether NA variation could influence the antigenic evolution of HA, we performed neutralizing antibody selection experiments using a panel of monoclonal antibodies targeting different HA epitopes. We found that the specific antibody escape profiles of HA were highly contingent upon NA background. Overall, our results indicate that natural variation in NA activity plays a significant role in governing the evolutionary potential of HA in the currently circulating pdmH1N1 lineage.

Deep mutational scanning of hepatitis B virus reveals a mechanism for cis-preferential reverse transcription

Hepatitis B virus (HBV) is a small double-stranded DNA virus that chronically infects 296 million people. Over half of its compact genome encodes proteins in two overlapping reading frames, and during evolution, multiple selective pressures can act on shared nucleotides. This study combines an RNA-based HBV cell culture system with deep mutational scanning (DMS) to uncouple cis- and trans-acting sequence requirements in the HBV genome. The results support a leaky ribosome scanning model for polymerase translation, provide a fitness map of the HBV polymerase at single-nucleotide resolution, and identify conserved prolines adjacent to the HBV polymerase termination codon that stall ribosomes. Further experiments indicated that stalled ribosomes tether the nascent polymerase to its template RNA, ensuring cis-preferential RNA packaging and reverse transcription of the HBV genome.

Minimum information and guidelines for reporting a multiplexed assay of variant effect

Multiplexed assays of variant effect (MAVEs) have emerged as a powerful approach for interrogating thousands of genetic variants in a single experiment. The flexibility and widespread adoption of these techniques across diverse disciplines have led to a heterogeneous mix of data formats and descriptions, which complicates the downstream use of the resulting datasets. To address these issues and promote reproducibility and reuse of MAVE data, we define a set of minimum information standards for MAVE data and metadata and outline a controlled vocabulary aligned with established biomedical ontologies for describing these experimental designs.

Natural variation in neuraminidase activity influences the evolutionary potential of the seasonal H1N1 lineage hemagglutinin

The antigenic evolution of the influenza A virus hemagglutinin (HA) gene poses a major challenge for the development of vaccines capable of eliciting long-term protection. Prior efforts to understand the mechanisms that govern viral antigenic evolution mainly focus on HA in isolation, ignoring the fact that HA must act in concert with the viral neuraminidase (NA) during replication and spread. Numerous studies have demonstrated that the degree to which the receptor binding avidity of HA and receptor cleaving activity of NA are balanced with each other influences overall viral fitness. We recently showed that changes in NA activity can significantly alter the mutational fitness landscape of HA in the context of a lab-adapted virus strain. Here, we test whether natural variation in relative NA activity can influence the evolutionary potential of HA in the context of the seasonal H1N1 lineage (pdmH1N1) that has circulated in humans since the 2009 pandemic. We observed substantial variation in the relative activities of NA proteins encoded by a panel of H1N1 vaccine strains isolated between 2009 and 2019. We comprehensively assessed the effect of NA background on the HA mutational fitness landscape in the circulating pdmH1N1 lineage using deep mutational scanning and observed pronounced epistasis between NA and residues in or near the receptor binding site of HA. To determine whether NA variation could influence the antigenic evolution of HA, we performed neutralizing antibody selection experiments using a panel of monoclonal antibodies targeting different HA epitopes. We found that the specific antibody escape profiles of HA were highly contingent upon NA background. Overall, our results indicate that natural variation in NA activity plays a significant role in governing the evolutionary potential of HA in the currently circulating pdmH1N1 lineage.

popDMS infers mutation effects from deep mutational scanning data

Deep mutational scanning (DMS) experiments provide a powerful method to measure the functional effects of genetic mutations at massive scales. However, the data generated from these experiments can be difficult to analyze, with significant variation between experimental replicates. To overcome this challenge, we developed popDMS, a computational method based on population genetics theory, to infer the functional effects of mutations from DMS data. Through extensive tests, we found that the functional effects of single mutations and epistasis inferred by popDMS are highly consistent across replicates, comparing favorably with existing methods. Our approach is flexible and can be widely applied to DMS data that includes multiple time points, multiple replicates, and different experimental conditions.

The effect of single mutations in Zika virus envelope on escape from broadly neutralizing antibodies

IMPORTANCE

The wide endemic range of mosquito-vectored flaviviruses-such as Zika virus and dengue virus serotypes 1-4-places hundreds of millions of people at risk of infection every year. Despite this, there are no widely available vaccines, and treatment of severe cases is limited to supportive care. An avenue toward development of more widely applicable vaccines and targeted therapies is the characterization of monoclonal antibodies that broadly neutralize all these viruses. Here, we measure how single amino acid mutations in viral envelope protein affect neutralizing antibodies with both broad and narrow specificities. We find that broadly neutralizing antibodies with potential as vaccine prototypes or biological therapeutics are quantifiably more difficult to escape than narrow, virus-specific neutralizing antibodies.

Deep mutational scanning reveals the functional constraints and evolutionary potential of the influenza A virus PB1 protein

IMPORTANCE

The influenza virus polymerase is important for adaptation to new hosts and, as a determinant of mutation rate, for the process of adaptation itself. We performed a deep mutational scan of the polymerase basic 1 (PB1) protein to gain insights into the structural and functional constraints on the influenza RNA-dependent RNA polymerase. We find that PB1 is highly constrained at specific sites that are only moderately predicted by the global structure or larger domain. We identified a number of beneficial mutations, many of which have been shown to be functionally important or observed in influenza virus' natural evolution. Overall, our atlas of PB1 mutations and their fitness impacts serves as an important resource for future studies of influenza replication and evolution.

Deep mutational scanning of proteins in mammalian cells

Protein mutagenesis is essential for unveiling the molecular mechanisms underlying protein function in health, disease, and evolution. In the past decade, deep mutational scanning methods have evolved to support the functional analysis of nearly all possible single-amino acid changes in a protein of interest. While historically these methods were developed in lower organisms such as E. coli and yeast, recent technological advancements have resulted in the increased use of mammalian cells, particularly for studying proteins involved in human disease. These advancements will aid significantly in the classification and interpretation of variants of unknown significance, which are being discovered at large scale due to the current surge in the use of whole-genome sequencing in clinical contexts. Here, we explore the experimental aspects of deep mutational scanning studies in mammalian cells and report the different methods used in each step of the workflow, ultimately providing a useful guide toward the design of such studies.

Key Amino Acid Residues That Determine the Antigenic Properties of Highly Pathogenic H5 Influenza Viruses Bearing the Clade 2.3.4.4 Hemagglutinin Gene

The H5 subtype highly pathogenic avian influenza viruses bearing the clade 2.3.4.4 HA gene have been pervasive among domestic poultry and wild birds worldwide since 2014, presenting substantial risks to human and animal health. Continued circulation of clade 2.3.4.4 viruses has resulted in the emergence of eight subclades (2.3.4.4a-h) and multiple distinct antigenic groups. However, the key antigenic substitutions responsible for the antigenic change of these viruses remain unknown. In this study, we analyzed the HA gene sequences of 5713 clade 2.3.4.4 viruses obtained from a public database and found that 23 amino acid residues were highly variable among these strains. We then generated a series of single-amino-acid mutants based on the H5-Re8 (a vaccine seed virus) background and tested their reactivity with a panel of eight monoclonal antibodies (mAbs). Six mutants bearing amino acid substitutions at positions 120, 126, 141, 156, 185, or 189 (H5 numbering) led to reduced or lost reactivity to these mAbs. Further antigenic cartography analysis revealed that the amino acid residues at positions 126, 156, and 189 acted as immunodominant epitopes of H5 viruses. Collectively, our findings offer valuable guidance for the surveillance and early detection of emerging antigenic variants.

Dynamics and durability of HIV-1 neutralization are determined by viral replication

Human immunodeficiency virus type 1 (HIV-1)-neutralizing antibodies (nAbs) that prevent infection are the main goal of HIV vaccine discovery. But as no nAb-eliciting vaccines are yet available, only data from HIV-1 neutralizers-persons with HIV-1 who naturally develop broad and potent nAbs-can inform about the dynamics and durability of nAb responses in humans, knowledge which is crucial for the design of future HIV-1 vaccine regimens. To address this, we assessed HIV-1-neutralizing immunoglobulin G (IgG) from 2,354 persons with HIV-1 on or off antiretroviral therapy (ART). Infection with non-clade B viruses, CD4+ T cell counts <200 µl-1, being off ART and a longer time off ART were independent predictors of a more potent and broad neutralization. In longitudinal analyses, we found nAb half-lives of 9.3 and 16.9 years in individuals with no- or low-level viremia, respectively, and 4.0 years in persons who newly initiated ART. Finally, in a potent HIV-1 neutralizer, we identified lower fractions of serum nAbs and of nAb-encoding memory B cells after ART initiation, suggesting that a decreasing neutralizing serum activity after antigen withdrawal is due to lower levels of nAbs. These results collectively show that HIV-1-neutralizing responses can persist for several years, even at low antigen levels, suggesting that an HIV-1 vaccine may elicit a durable nAb response.

Comprehensive profiling of neutralizing polyclonal sera targeting coxsackievirus B3

Despite their fundamental role in resolving viral infections, our understanding of how polyclonal neutralizing antibody responses target non-enveloped viruses remains limited. To define these responses, we obtained the full antigenic profile of multiple human and mouse polyclonal sera targeting the capsid of a prototypical picornavirus, coxsackievirus B3. Our results uncover significant variation in the breadth and strength of neutralization sites targeted by individual human polyclonal responses, which contrasted with homogenous responses observed in experimentally infected mice. We further use these comprehensive antigenic profiles to define key structural and evolutionary parameters that are predictive of escape, assess epitope dominance at the population level, and reveal a need for at least two mutations to achieve significant escape from multiple sera. Overall, our data provide a comprehensive analysis of how polyclonal sera target a non-enveloped viral capsid and help define both immune dominance and escape at the population level.

MBE: model-based enrichment estimation and prediction for differential sequencing data

Characterizing differences in sequences between two conditions, such as with and without drug exposure, using high-throughput sequencing data is a prevalent problem involving quantifying changes in sequence abundances, and predicting such differences for unobserved sequences. A key shortcoming of current approaches is their extremely limited ability to share information across related but non-identical reads. Consequently, they cannot use sequencing data effectively, nor be directly applied in many settings of interest. We introduce model-based enrichment (MBE) to overcome this shortcoming. We evaluate MBE using both simulated and real data. Overall, MBE improves accuracy compared to current differential analysis methods.

The effect of single mutations in Zika virus envelope on escape from broadly neutralizing antibodies

Zika virus and dengue virus are co-circulating flaviviruses with a widespread endemic range. Eliciting broad and potent neutralizing antibodies is an attractive goal for developing a vaccine to simultaneously protect against these viruses. However, the capacity of viral mutations to confer escape from broadly neutralizing antibodies remains undescribed, due in part to limited throughput and scope of traditional approaches. Here, we use deep mutational scanning to map how all possible single amino acid mutations in Zika virus envelope protein affect neutralization by antibodies of varying breadth and potency. While all antibodies selected viral escape mutations, the mutations selected by broadly neutralizing antibodies conferred less escape relative to those selected by narrow, virus-specific antibodies. Surprisingly, even for broadly neutralizing antibodies with similar binding footprints, different single mutations led to escape, indicating distinct functional requirements for neutralization not captured by existing structures. Additionally, the antigenic effects of mutations selected by broadly neutralizing antibodies were conserved across divergent, albeit related, flaviviruses. Our approach identifies residues critical for antibody neutralization, thus comprehensively defining the as-yet-unknown functional epitopes of antibodies with clinical potential.

Facile Method for High-throughput Identification of Stabilizing Mutations

Characterizing the effects of mutations on stability is critical for understanding the function and evolution of proteins and improving their biophysical properties. High throughput folding and abundance assays have been successfully used to characterize missense mutations associated with reduced stability. However, screening for increased thermodynamic stability is more challenging since such mutations are rarer and their impact on assay readout is more subtle. Here, a multiplex assay for high throughput screening of protein folding was developed by combining deep mutational scanning, fluorescence-activated cell sorting, and deep sequencing. By analyzing a library of 2000 variants of Adenylate kinase we demonstrate that the readout of the method correlates with stability and that mutants with up to 13 °C increase in thermal melting temperature could be identified with low false positive rate. The discovery of many stabilizing mutations also enabled the analysis of general substitution patterns associated with increased stability in Adenylate kinase. This high throughput method to identify stabilizing mutations can be combined with functional screens to identify mutations that improve both stability and activity.

Position-Specific Enrichment Ratio Matrix scores predict antibody variant properties from deep sequencing data

MOTIVATION

Deep sequencing of antibody and related protein libraries after phage or yeast-surface display sorting is widely used to identify variants with increased affinity, specificity, and/or improvements in key biophysical properties. Conventional approaches for identifying optimal variants typically use the frequencies of observation in enriched libraries or the corresponding enrichment ratios. However, these approaches disregard the vast majority of deep sequencing data and often fail to identify the best variants in the libraries.

RESULTS

Here, we present a method, Position-Specific Enrichment Ratio Matrix (PSERM) scoring, that uses entire deep sequencing datasets from pre- and post-selections to score each observed protein variant. The PSERM scores are the sum of the site-specific enrichment ratios observed at each mutated position. We find that PSERM scores are much more reproducible and correlate more strongly with experimentally measured properties than frequencies or enrichment ratios, including for multiple antibody properties (affinity and non-specific binding) for a clinical-stage antibody (emibetuzumab). We expect that this method will be broadly applicable to diverse protein engineering campaigns.

AVAILABILITY AND IMPLEMENTATION

All deep sequencing datasets and code to perform the analyses presented within are available via https://github.com/Tessier-Lab-UMich/PSERM_paper.

Position-Specific Enrichment Ratio Matrix scores predict antibody variant properties from deep sequencing data

Motivation

Deep sequencing of antibody and related protein libraries after phage or yeast-surface display sorting is widely used to identify variants with increased affinity, specificity and/or improvements in key biophysical properties. Conventional approaches for identifying optimal variants typically use the frequencies of observation in enriched libraries or the corresponding enrichment ratios. However, these approaches disregard the vast majority of deep sequencing data and often fail to identify the best variants in the libraries.

Results

Here, we present a method, Position-Specific Enrichment Ratio Matrix (PSERM) scoring, that uses entire deep sequencing datasets from pre- and post-selections to score each observed protein variant. The PSERM scores are the sum of the site-specific enrichment ratios observed at each mutated position. We find that PSERM scores are much more reproducible and correlate more strongly with experimentally measured properties than frequencies or enrichment ratios, including for multiple antibody properties (affinity and non-specific binding) for a clinical-stage antibody (emibetuzumab). We expect that this method will be broadly applicable to diverse protein engineering campaigns.

Availability

All deep sequencing datasets and code to do the analyses presented within are available via GitHub.

Contact

Peter Tessier, ptessier@umich.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

satmut_utils: a simulation and variant calling package for multiplexed assays of variant effect

The impact of millions of individual genetic variants on molecular phenotypes in coding sequences remains unknown. Multiplexed assays of variant effect (MAVEs) are scalable methods to annotate relevant variants, but existing software lacks standardization, requires cumbersome configuration, and does not scale to large targets. We present satmut_utils as a flexible solution for simulation and variant quantification. We then benchmark MAVE software using simulated and real MAVE data. We finally determine mRNA abundance for thousands of cystathionine beta-synthase variants using two experimental methods. The satmut_utils package enables high-performance analysis of MAVEs and reveals the capability of variants to alter mRNA abundance.

Identification of ADAR1 p150 and p110 Associated Edit Sites

Adenosine deaminase acting on RNA 1 (ADAR1) catalyzes adenosine-to-inosine editing on double-stranded RNA molecules and is involved in regulating cellular responses to endogenous and exogenous RNA. ADAR1 is the primary A-to-I editor of RNA in humans, and the majority of edit sites are found in a class of short interspersed nuclear elements called Alu elements, many of which are located in introns and 3' untranslated regions. Two ADAR1 protein isoforms, p110 (110 kDa) and p150 (150 kDa), are known to be coupled in expression, and decoupling the expression of these isoforms has revealed that the p150 isoform edits a broader range of targets compared to p110. Numerous methods for identification of ADAR1-associated edits have been developed, and we present here a specific method for identification of edit sites associated with individual ADAR1 isoforms.

Scalable Functional Assays for the Interpretation of Human Genetic Variation

Scalable sequence-function studies have enabled the systematic analysis and cataloging of hundreds of thousands of coding and noncoding genetic variants in the human genome. This has improved clinical variant interpretation and provided insights into the molecular, biophysical, and cellular effects of genetic variants at an astonishing scale and resolution across the spectrum of allele frequencies. In this review, we explore current applications and prospects for the field and outline the principles underlying scalable functional assay design, with a focus on the study of single-nucleotide coding and noncoding variants.

Intragenic compensation through the lens of deep mutational scanning

A significant fraction of mutations in proteins are deleterious and result in adverse consequences for protein function, stability, or interaction with other molecules. Intragenic compensation is a specific case of positive epistasis when a neutral missense mutation cancels effect of a deleterious mutation in the same protein. Permissive compensatory mutations facilitate protein evolution, since without them all sequences would be extremely conserved. Understanding compensatory mechanisms is an important scientific challenge at the intersection of protein biophysics and evolution. In human genetics, intragenic compensatory interactions are important since they may result in variable penetrance of pathogenic mutations or fixation of pathogenic human alleles in orthologous proteins from related species. The latter phenomenon complicates computational and clinical inference of an allele's pathogenicity. Deep mutational scanning is a relatively new technique that enables experimental studies of functional effects of thousands of mutations in proteins. We review the important aspects of the field and discuss existing limitations of current datasets. We reviewed ten published DMS datasets with quantified functional effects of single and double mutations and described rates and patterns of intragenic compensation in eight of them.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-022-01005-w.

Co-evolution of interacting proteins through non-contacting and non-specific mutations

Proteins often accumulate neutral mutations that do not affect current functions but can profoundly influence future mutational possibilities and functions. Understanding such hidden potential has major implications for protein design and evolutionary forecasting but has been limited by a lack of systematic efforts to identify potentiating mutations. Here, through the comprehensive analysis of a bacterial toxin-antitoxin system, we identified all possible single substitutions in the toxin that enable it to tolerate otherwise interface-disrupting mutations in its antitoxin. Strikingly, the majority of enabling mutations in the toxin do not contact and promote tolerance non-specifically to many different antitoxin mutations, despite covariation in homologues occurring primarily between specific pairs of contacting residues across the interface. In addition, the enabling mutations we identified expand future mutational paths that both maintain old toxin-antitoxin interactions and form new ones. These non-specific mutations are missed by widely used covariation and machine learning methods. Identifying such enabling mutations will be critical for ensuring continued binding of therapeutically relevant proteins, such as antibodies, aimed at evolving targets.

MAVE-NN: learning genotype-phenotype maps from multiplex assays of variant effect

Multiplex assays of variant effect (MAVEs) are a family of methods that includes deep mutational scanning experiments on proteins and massively parallel reporter assays on gene regulatory sequences. Despite their increasing popularity, a general strategy for inferring quantitative models of genotype-phenotype maps from MAVE data is lacking. Here we introduce MAVE-NN, a neural-network-based Python package that implements a broadly applicable information-theoretic framework for learning genotype-phenotype maps-including biophysically interpretable models-from MAVE datasets. We demonstrate MAVE-NN in multiple biological contexts, and highlight the ability of our approach to deconvolve mutational effects from otherwise confounding experimental nonlinearities and noise.

The endoplasmic reticulum proteostasis network profoundly shapes the protein sequence space accessible to HIV envelope

The sequence space accessible to evolving proteins can be enhanced by cellular chaperones that assist biophysically defective clients in navigating complex folding landscapes. It is also possible, at least in theory, for proteostasis mechanisms that promote strict quality control to greatly constrain accessible protein sequence space. Unfortunately, most efforts to understand how proteostasis mechanisms influence evolution rely on artificial inhibition or genetic knockdown of specific chaperones. The few experiments that perturb quality control pathways also generally modulate the levels of only individual quality control factors. Here, we use chemical genetic strategies to tune proteostasis networks via natural stress response pathways that regulate the levels of entire suites of chaperones and quality control mechanisms. Specifically, we upregulate the unfolded protein response (UPR) to test the hypothesis that the host endoplasmic reticulum (ER) proteostasis network shapes the sequence space accessible to human immunodeficiency virus-1 (HIV-1) envelope (Env) protein. Elucidating factors that enhance or constrain Env sequence space is critical because Env evolves extremely rapidly, yielding HIV strains with antibody- and drug-escape mutations. We find that UPR-mediated upregulation of ER proteostasis factors, particularly those controlled by the IRE1-XBP1s UPR arm, globally reduces Env mutational tolerance. Conserved, functionally important Env regions exhibit the largest decreases in mutational tolerance upon XBP1s induction. Our data indicate that this phenomenon likely reflects strict quality control endowed by XBP1s-mediated remodeling of the ER proteostasis environment. Intriguingly, and in contrast, specific regions of Env, including regions targeted by broadly neutralizing antibodies, display enhanced mutational tolerance when XBP1s is induced, hinting at a role for host proteostasis network hijacking in potentiating antibody escape. These observations reveal a key function for proteostasis networks in decreasing instead of expanding the sequence space accessible to client proteins, while also demonstrating that the host ER proteostasis network profoundly shapes the mutational tolerance of Env in ways that could have important consequences for HIV adaptation.

The host cell’s endoplasmic reticulum proteostasis network has a profound, constraining impact on the protein sequence space accessible to HIV’s envelope protein, which is a major target of the host’s adaptive immune system; in particular, upregulation of stringent quality control pathways appears to restrict the viability of destabilizing envelope variants.

Broadly neutralizing antibodies target a haemagglutinin anchor epitope

Broadly neutralizing antibodies that target epitopes of haemagglutinin on the influenza virus have the potential to provide near universal protection against influenza virus infection1. However, viral mutants that escape broadly neutralizing antibodies have been reported2,3. The identification of broadly neutralizing antibody classes that can neutralize viral escape mutants is critical for universal influenza virus vaccine design. Here we report a distinct class of broadly neutralizing antibodies that target a discrete membrane-proximal anchor epitope of the haemagglutinin stalk domain. Anchor epitope-targeting antibodies are broadly neutralizing across H1 viruses and can cross-react with H2 and H5 viruses that are a pandemic threat. Antibodies that target this anchor epitope utilize a highly restricted repertoire, which encodes two public binding motifs that make extensive contacts with conserved residues in the fusion peptide. Moreover, anchor epitope-targeting B cells are common in the human memory B cell repertoire and were recalled in humans by an oil-in-water adjuvanted chimeric haemagglutinin vaccine4,5, which is a potential universal influenza virus vaccine. To maximize protection against seasonal and pandemic influenza viruses, vaccines should aim to boost this previously untapped source of broadly neutralizing antibodies that are widespread in the human memory B cell pool.

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.

Deep mutational scanning for therapeutic antibody engineering

The biophysical and functional properties of monoclonal antibody (mAb) drug candidates are often improved by protein engineering methods to increase the probability of clinical efficacy. One emerging method is deep mutational scanning (DMS) which combines the power of exhaustive protein mutagenesis and functional screening with deep sequencing and bioinformatics. The application of DMS has yielded significant improvements to the affinity, specificity, and stability of several preclinical antibodies alongside novel applications such as introducing multi-specific binding properties. DMS has also been applied directly on target antigens to precisely map antibody-binding epitopes and notably to profile the mutational escape potential of viral targets (e.g., SARS-CoV-2 variants). Finally, DMS combined with machine learning is enabling advances in the computational screening and engineering of therapeutic antibodies.

Intelligent host engineering for metabolic flux optimisation in biotechnology

Optimising the function of a protein of length N amino acids by directed evolution involves navigating a 'search space' of possible sequences of some 20N. Optimising the expression levels of P proteins that materially affect host performance, each of which might also take 20 (logarithmically spaced) values, implies a similar search space of 20P. In this combinatorial sense, then, the problems of directed protein evolution and of host engineering are broadly equivalent. In practice, however, they have different means for avoiding the inevitable difficulties of implementation. The spare capacity exhibited in metabolic networks implies that host engineering may admit substantial increases in flux to targets of interest. Thus, we rehearse the relevant issues for those wishing to understand and exploit those modern genome-wide host engineering tools and thinking that have been designed and developed to optimise fluxes towards desirable products in biotechnological processes, with a focus on microbial systems. The aim throughput is 'making such biology predictable'. Strategies have been aimed at both transcription and translation, especially for regulatory processes that can affect multiple targets. However, because there is a limit on how much protein a cell can produce, increasing kcat in selected targets may be a better strategy than increasing protein expression levels for optimal host engineering.

Shifts in amino acid preferences as proteins evolve: A synthesis of experimental and theoretical work

Amino acid preferences vary across sites and time. While variation across sites is widely accepted, the extent and frequency of temporal shifts are contentious. Our understanding of the drivers of amino acid preference change is incomplete: To what extent are temporal shifts driven by adaptive versus nonadaptive evolutionary processes? We review phenomena that cause preferences to vary (e.g., evolutionary Stokes shift, contingency, and entrenchment) and clarify how they differ. To determine the extent and prevalence of shifted preferences, we review experimental and theoretical studies. Analyses of natural sequence alignments often detect decreases in homoplasy (convergence and reversions) rates, and variation in replacement rates with time-signals that are consistent with temporally changing preferences. While approaches inferring shifts in preferences from patterns in natural alignments are valuable, they are indirect since multiple mechanisms (both adaptive and nonadaptive) could lead to the observed signal. Alternatively, site-directed mutagenesis experiments allow for a more direct assessment of shifted preferences. They corroborate evidence from multiple sequence alignments, revealing that the preference for an amino acid at a site varies depending on the background sequence. However, shifts in preferences are usually minor in magnitude and sites with significantly shifted preferences are low in frequency. The small yet consistent perturbations in preferences could, nevertheless, jeopardize the accuracy of inference procedures, which assume constant preferences. We conclude by discussing if and how such shifts in preferences might influence widely used time-homogenous inference procedures and potential ways to mitigate such effects.

Linking genome variants to disease: scalable approaches to test the functional impact of human mutations

The application of genomics to medicine has accelerated the discovery of mutations underlying disease and has enhanced our knowledge of the molecular underpinnings of diverse pathologies. As the amount of human genetic material queried via sequencing has grown exponentially in recent years, so too has the number of rare variants observed. Despite progress, our ability to distinguish which rare variants have clinical significance remains limited. Over the last decade, however, powerful experimental approaches have emerged to characterize variant effects orders of magnitude faster than before. Fueled by improved DNA synthesis and sequencing and, more recently, by CRISPR/Cas9 genome editing, multiplex functional assays provide a means of generating variant effect data in wide-ranging experimental systems. Here, I review recent applications of multiplex assays that link human variants to disease phenotypes and I describe emerging strategies that will enhance their clinical utility in coming years.

Highly-potent, synthetic APOBEC3s restrict HIV-1 through deamination-independent mechanisms

The APOBEC3 (A3) genes encode cytidine deaminase proteins with potent antiviral and anti-retroelement activity. This locus is characterized by duplication, recombination, and deletion events that gave rise to the seven A3s found in primates. These include three single deaminase domain A3s (A3A, A3C, and A3H) and four double deaminase domain A3s (A3B, A3D, A3F, and A3G). The most potent of the A3 proteins against HIV-1 is A3G. However, it is not clear if double deaminase domain A3s have a generalized functional advantage to restrict HIV-1. In order to test whether superior restriction factors could be created by genetically linking single A3 domains into synthetic double domains, we linked A3C and A3H single domains in novel combinations. We found that A3C/A3H double domains acquired enhanced antiviral activity that is at least as potent, if not better than, A3G. Although these synthetic double domain A3s package into budding virions more efficiently than their respective single domains, this does not fully explain their gain of antiviral potency. The antiviral activity is conferred both by cytidine-deaminase dependent and independent mechanisms, with the latter correlating to an increase in RNA binding affinity. T cell lines expressing this A3C-A3H super restriction factor are able to control replicating HIV-1ΔVif infection to similar levels as A3G. Together, these data show that novel combinations of A3 domains are capable of gaining potent antiviral activity to levels similar to the most potent genome-encoded A3s, via a primarily non-catalytic mechanism.

Applications of Deep Mutational Scanning in Virology

Several recently developed high-throughput techniques have changed the field of molecular virology. For example, proteomics studies reveal complete interactomes of a viral protein, genome-wide CRISPR knockout and activation screens probe the importance of every single human gene in aiding or fighting a virus, and ChIP-seq experiments reveal genome-wide epigenetic changes in response to infection. Deep mutational scanning is a relatively novel form of protein science which allows the in-depth functional analysis of every nucleotide within a viral gene or genome, revealing regions of importance, flexibility, and mutational potential. In this review, we discuss the application of this technique to RNA viruses including members of the Flaviviridae family, Influenza A Virus and Severe Acute Respiratory Syndrome Coronavirus 2. We also briefly discuss the reverse genetics systems which allow for analysis of viral replication cycles, next-generation sequencing technologies and the bioinformatics tools that facilitate this research.

High-resolution mapping of the neutralizing and binding specificities of polyclonal sera post-HIV Env trimer vaccination

Mapping polyclonal serum responses is critical to rational vaccine design. However, most high-resolution mapping approaches involve isolating and characterizing individual antibodies, which incompletely defines the polyclonal response. Here we use two complementary approaches to directly map the specificities of the neutralizing and binding antibodies of polyclonal anti-HIV-1 sera from rabbits immunized with BG505 Env SOSIP trimers. We used mutational antigenic profiling to determine how all mutations in Env affected viral neutralization and electron microscopy polyclonal epitope mapping (EMPEM) to directly visualize serum Fabs bound to Env trimers. The dominant neutralizing specificities were generally only a subset of the more diverse binding specificities. Additional differences between binding and neutralization reflected antigenicity differences between virus and soluble Env trimer. Furthermore, we refined residue-level epitope specificity directly from sera, revealing subtle differences across sera. Together, mutational antigenic profiling and EMPEM yield a holistic view of the binding and neutralizing specificity of polyclonal sera.

Globally defining the effects of mutations in a picornavirus capsid

The capsids of non-enveloped viruses are highly multimeric and multifunctional protein assemblies that play key roles in viral biology and pathogenesis. Despite their importance, a comprehensive understanding of how mutations affect viral fitness across different structural and functional attributes of the capsid is lacking. To address this limitation, we globally define the effects of mutations across the capsid of a human picornavirus. Using this resource, we identify structural and sequence determinants that accurately predict mutational fitness effects, refine evolutionary analyses, and define the sequence specificity of key capsid-encoded motifs. Furthermore, capitalizing on the derived sequence requirements for capsid-encoded protease cleavage sites, we implement a bioinformatic approach for identifying novel host proteins targeted by viral proteases. Our findings represent the most comprehensive investigation of mutational fitness effects in a picornavirus capsid to date and illuminate important aspects of viral biology, evolution, and host interactions.

A virus is made up of genetic material that is encased with a protective protein coat called the capsid. The capsid also helps the virus to infect host cells by binding to the host receptor proteins and releasing its genetic material. Inside the cell, the virus hitchhikes the infected cell’s machinery to grow or replicate its own genetic material.

Viral capsids are the main target of the host’s defence system, and therefore, continuously change in an attempt to escape the immune system by introducing alterations (known as mutations) into the genes encoding viral capsid proteins. Mutations occur randomly, and so while some changes to the viral capsid might confer an advantage, others may have no effect at all, or even weaken the virus.

To better understand the effect of capsid mutations on the virus’ ability to infect host cells, Mattenberger et al. studied the Coxsackievirus B3, which is linked to heart problems and acute heart failure in humans. The researchers analysed around 90% of possible amino acid mutations (over 14,800 mutations) and correlated each mutation to how it influenced the virus’ ability to replicate in human cells grown in the laboratory.

Based on these results, Mattenberger et al. developed a computer model to predict how a particular mutation might affect the virus. The analysis also identified specific amino acid sequences of capsid proteins that are essential for certain tasks, such as building the capsid. It also included an analysis of sequences in the capsid that allow it to be recognized by another viral protein, which cuts the capsid proteins into the right size from a larger precursor. By looking for similar sequences in human genes, the researchers identified several ones that the virus may attack and inactivate to support its own replication.

These findings may help identify potential drug targets to develop new antiviral therapies. For example, proteins of the capsid that are less likely to mutate will provide a better target as they lower the possibility of the virus to become resistant to the treatment. They also highlight new proteins in human cells that could potentially block the virus in cells.

Inferring Protein Sequence-Function Relationships with Large-Scale Positive-Unlabeled Learning

Machine learning can infer how protein sequence maps to function without requiring a detailed understanding of the underlying physical or biological mechanisms. It is challenging to apply existing supervised learning frameworks to large-scale experimental data generated by deep mutational scanning (DMS) and related methods. DMS data often contain high-dimensional and correlated sequence variables, experimental sampling error and bias, and the presence of missing data. Notably, most DMS data do not contain examples of negative sequences, making it challenging to directly estimate how sequence affects function. Here, we develop a positive-unlabeled (PU) learning framework to infer sequence-function relationships from large-scale DMS data. Our PU learning method displays excellent predictive performance across ten large-scale sequence-function datasets, representing proteins of different folds, functions, and library types. The estimated parameters pinpoint key residues that dictate protein structure and function. Finally, we apply our statistical sequence-function model to design highly stabilized enzymes.

Structural and functional characterization of G protein-coupled receptors with deep mutational scanning

The >800 human G protein–coupled receptors (GPCRs) are responsible for transducing diverse chemical stimuli to alter cell state- and are the largest class of drug targets. Their myriad structural conformations and various modes of signaling make it challenging to understand their structure and function. Here, we developed a platform to characterize large libraries of GPCR variants in human cell lines with a barcoded transcriptional reporter of G protein signal transduction. We tested 7800 of 7828 possible single amino acid substitutions to the beta-2 adrenergic receptor (β₂AR) at four concentrations of the agonist isoproterenol. We identified residues specifically important for β₂AR signaling, mutations in the human population that are potentially loss of function, and residues that modulate basal activity. Using unsupervised learning, we identify residues critical for signaling, including all major structural motifs and molecular interfaces. We also find a previously uncharacterized structural latch spanning the first two extracellular loops that is highly conserved across Class A GPCRs and is conformationally rigid in both the inactive and active states of the receptor. More broadly, by linking deep mutational scanning with engineered transcriptional reporters, we establish a generalizable method for exploring pharmacogenomics, structure and function across broad classes of drug receptors.

Identification of Antibodies Targeting the H3N2 Hemagglutinin Receptor Binding Site following Vaccination of Humans

Antibodies targeting the receptor binding site (RBS) of the influenza virus hemagglutinin (HA) protein are usually not broadly reactive because their footprints are typically large and extend to nearby variable HA residues. Here, we identify several human H3N2 HA RBS-targeting monoclonal antibodies (mAbs) that are sensitive to substitutions in conventional antigenic sites and are therefore not broadly reactive. However, we also identify an H3N2 HA RBS-targeting mAb that is exceptionally broadly reactive despite being sensitive to substitutions in residues outside of the RBS. We show that similar antibodies are present at measurable levels in the sera of some individuals but that they are inefficiently elicited by conventional vaccines. Our data indicate that HA RBS-targeting antibodies can be effective against variable viral strains even when they are somewhat sensitive to substitutions in HA residues adjacent to the RBS.

Zost et al. show that most antibodies targeting the RBS of the H3N2 HAs are not broadly reactive. They identify one broadly reactive H3 HA RBS antibody that is tolerant of substitutions in adjacent antigenic sites but show that these types of antibodies are not efficiently elicited by vaccination.

DiMSum: an error model and pipeline for analyzing deep mutational scanning data and diagnosing common experimental pathologies

Deep mutational scanning (DMS) enables multiplexed measurement of the effects of thousands of variants of proteins, RNAs, and regulatory elements. Here, we present a customizable pipeline, DiMSum, that represents an end-to-end solution for obtaining variant fitness and error estimates from raw sequencing data. A key innovation of DiMSum is the use of an interpretable error model that captures the main sources of variability arising in DMS workflows, outperforming previous methods. DiMSum is available as an R/Bioconda package and provides summary reports to help researchers diagnose common DMS pathologies and take remedial steps in their analyses.

Implications of localized charge for human influenza A H1N1 hemagglutinin evolution: Insights from deep mutational scans

Seasonal influenza A viruses of humans evolve rapidly due to strong selection pressures from host immune responses, principally on the hemagglutinin (HA) viral surface protein. Based on mouse transmission experiments, a proposed mechanism for immune evasion consists of increased avidity to host cellular receptors, mediated by electrostatic charge interactions with negatively charged cell surfaces. In support of this, the HA charge of the globally circulating H3N2 has increased over time since its pandemic. However, the same trend was not seen in H1N1 HA sequences. This is counter-intuitive, since immune escape due to increased avidity (due itself to an increase in charge) was determined experimentally. Here, we explore whether patterns of local charge of H1N1 HA can explain this discrepancy and thus further associate electrostatic charge with immune escape and viral evolutionary dynamics. Measures of site-wise functional selection and expected charge computed from deep mutational scan data on an early H1N1 HA yield a striking division of residues into three groups, separated by charge. We then explored evolutionary dynamics of these groups from 1918 to 2008. In particular, one group increases in net charge over time and consists of sites that are evolving the fastest, that are closest to the receptor binding site (RBS), and that are exposed to solvent (i.e., on the surface). By contrast, another group decreases in net charge and consists of sites that are further away from the RBS and evolving slower, but also exposed to solvent. The last group consists of those sites in the HA core, with no change in net charge and that evolve very slowly. Thus, there is a group of residues that follows the same trend as seen for the entire H3N2 HA. It is possible that the H1N1 HA is under other biophysical constraints that result in compensatory decreases in charge elsewhere on the protein. Our results implicate localized charge in HA interactions with host cells, and highlight how deep mutational scan data can inform evolutionary hypotheses.

The hemagglutinin (HA) surface protein of influenza A viruses evolves rapidly to evade host immunity, leading to sizable yearly epidemics in human populations. Previous transmission experiments with H1N1 in mice have tied immune escape to an increase in HA avidity for cellular receptors, mediated by electrostatic charge. Furthermore, retrospective sequence analyses from a previous study confirmed that the HA of circulating global H3N2 has increased in net charge, yet surprisingly, that of H1N1 has not varied significantly. How is a stable net charge related to local patterns of H1N1 HA charge in response to selection? To elucidate the role of local electrostatic charge in host-virus interactions, we investigate characteristics of local charge on the H1N1 HA using functional data from deep mutational scan experiments. Combining measures of functional selection and expected charge at each site on DMS data from a 1933 H1N1 HA yields a striking visual pattern that identifies three groups of sites that have different biophysical properties and that prove to have distinct evolutionary patterns in natural human sequences. Essentially, we find evidence for an increase in charge near the receptor binding site and for compensatory changes elsewhere on the protein. Thus, our findings may reconcile disparate results from transmission experiments and natural sequence analyses, and highlight the importance of local properties of the HA protein. Overall, our findings further support the hypothesis of immune escape due to increased HA avidity to cellular receptors mediated by electrostatic charge. More generally, our work illustrates a novel method of leveraging deep mutational scans in conjunction with natural sequences to shed light on virus-host interactions.

Different genetic barriers for resistance to HA stem antibodies in influenza H3 and H1 viruses

The discovery and characterization of broadly neutralizing human antibodies (bnAbs) to the highly conserved stem region of influenza hemagglutinin (HA) have contributed to considerations of a universal influenza vaccine. However, the potential for resistance to stem bnAbs also needs to be more thoroughly evaluated. Using deep mutational scanning, with a focus on epitope residues, we found that the genetic barrier to resistance to stem bnAbs is low for the H3 subtype but substantially higher for the H1 subtype owing to structural differences in the HA stem. Several strong resistance mutations in H3 can be observed in naturally circulating strains and do not reduce in vitro viral fitness and in vivo pathogenicity. This study highlights a potential challenge for development of a truly universal influenza vaccine.

Improved mutant function prediction via PACT: Protein Analysis and Classifier Toolkit

MOTIVATION

Deep mutational scanning experiments have enabled the measurement of the sequence-function relationship for thousands of mutations in a single experiment. The Protein Analysis and Classifier Toolkit (PACT) is a Python software package that marries the fitness metric of a given mutation within these experiments to sequence and structural features enabling downstream analyses. PACT enables the easy development of user sharable protocols for custom deep mutational scanning experiments as all code is modular and reusable between protocols. Protocols for mutational libraries with single or multiple mutations are included. To exemplify its utility, PACT assessed two deep mutational scanning datasets that measured the tradeoff of enzyme activity and enzyme stability.

RESULTS

PACT efficiently evaluated classifiers that predict protein mutant function tested on deep mutational scanning screens. We found that the classifiers with the lowest false positive and highest true positive rate assesses sequence homology, contact number and if mutation involves proline.

AVAILABILITY AND IMPLEMENTATION

PACT and the processed datasets are distributed freely under the terms of the GPL-3 license. The source code is available at GitHub (https://github.com/JKlesmith/PACT).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

APOBEC3C Tandem Domain Proteins Create Super Restriction Factors against HIV-1

As a part of the innate immune system, humans encode proteins that inhibit viruses such as HIV-1. These broadly acting antiviral proteins do not protect humans from viral infections because viruses encode proteins that antagonize the host antiviral proteins to evade the innate immune system. One such example of a host antiviral protein is APOBEC3C (A3C), which weakly inhibits HIV-1. Here, we show that we can improve the antiviral activity of A3C by duplicating the DNA sequence to create a synthetic tandem domain and, furthermore, that the proteins thus generated are relatively resistant to the viral antagonist Vif. Together, these data give insights about how nature has evolved a defense against viral pathogens such as HIV.

Humans encode proteins, called restriction factors, that inhibit replication of viruses such as HIV-1. The members of one family of antiviral proteins, apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3; shortened here to A3), act by deaminating cytidines to uridines during the reverse transcription reaction of HIV-1. The A3 locus encodes seven genes, named A3A to A3H. These genes have either one or two cytidine deaminase domains, and several of these A3s potently restrict HIV-1. A3C, which has only a single cytidine deaminase domain, however, inhibits HIV-1 only very weakly. We tested novel double domain protein combinations by genetically linking two A3C genes to make a synthetic tandem domain protein. This protein created a “super restriction factor” that had more potent antiviral activity than the native A3C protein, which correlated with increased packaging into virions. Furthermore, disabling one of the active sites of the synthetic tandem domain protein resulted in an even greater increase in the antiviral activity—recapitulating a similar evolution seen in A3F and A3G (double domain A3s that use only a single catalytically active deaminase domain). These A3C tandem domain proteins do not have an increase in mutational activity but instead inhibit formation of reverse transcription products, which correlates with their ability to form large higher-order complexes in cells. Finally, the A3C-A3C super restriction factor largely escaped antagonism by the HIV-1 viral protein Vif.

Antibody Lineages with Vaccine-Induced Antigen-Binding Hotspots Develop Broad HIV Neutralization

The vaccine-mediated elicitation of antibodies (Abs) capable of neutralizing diverse HIV-1 strains has been a long-standing goal. To understand how broadly neutralizing antibodies (bNAbs) can be elicited, we identified, characterized, and tracked five neutralizing Ab lineages targeting the HIV-1-fusion peptide (FP) in vaccinated macaques over time. Genetic and structural analyses revealed two of these lineages to belong to a reproducible class capable of neutralizing up to 59% of 208 diverse viral strains. B cell analysis indicated each of the five lineages to have been initiated and expanded by FP-carrier priming, with envelope (Env)-trimer boosts inducing cross-reactive neutralization. These Abs had binding-energy hotspots focused on FP, whereas several FP-directed Abs induced by immunization with Env trimer-only were less FP-focused and less broadly neutralizing. Priming with a conserved subregion, such as FP, can thus induce Abs with binding-energy hotspots coincident with the target subregion and capable of broad neutralization.

Identification of HIV-1 Envelope Mutations that Enhance Entry Using Macaque CD4 and CCR5

Although Rhesus macaques are an important animal model for HIV-1 vaccine development research, most transmitted HIV-1 strains replicate poorly in macaque cells. A major genetic determinant of this species-specific restriction is a non-synonymous mutation in macaque CD4 that results in reduced HIV-1 Envelope (Env)-mediated viral entry compared to human CD4. Recent research efforts employing either laboratory evolution or structure-guided design strategies have uncovered several mutations in Env’s gp120 subunit that enhance binding of macaque CD4 by transmitted/founder HIV-1 viruses. In order to identify additional Env mutations that promote infection of macaque cells, we utilized deep mutational scanning to screen thousands of Env point mutants for those that enhance HIV-1 entry via macaque receptors. We identified many uncharacterized amino acid mutations in the N-terminal heptad repeat (NHR) and C-terminal heptad repeat (CHR) regions of gp41 that increased entry into cells bearing macaque receptors up to 9-fold. Many of these mutations also modestly increased infection of cells bearing human CD4 and CCR5 (up to 1.5-fold). NHR/CHR mutations identified by deep mutational scanning that enhanced entry also increased sensitivity to neutralizing antibodies targeting the MPER epitope, and to inactivation by cold-incubation, suggesting that they promote sampling of an intermediate trimer conformation between closed and receptor bound states. Identification of this set of mutations can inform future macaque model studies, and also further our understanding of the relationship between Env structure and function.

Restriction of HIV-1 Escape by a Highly Broad and Potent Neutralizing Antibody

Broadly neutralizing antibodies (bNAbs) represent a promising approach to prevent and treat HIV-1 infection. However, viral escape through mutation of the HIV-1 envelope glycoprotein (Env) limits clinical applications. Here we describe 1-18, a new V_H1-46-encoded CD4 binding site (CD4bs) bNAb with outstanding breadth (97%) and potency (GeoMean IC₅₀ = 0.048 μg/mL). Notably, 1-18 is not susceptible to typical CD4bs escape mutations and effectively overcomes HIV-1 resistance to other CD4bs bNAbs. Moreover, mutational antigenic profiling uncovered restricted pathways of HIV-1 escape. Of most promise for therapeutic use, even 1-18 alone fully suppressed viremia in HIV-1-infected humanized mice without selecting for resistant viral variants. A 2.5-Å cryo-EM structure of a 1-18-BG505_SOSIP.664 Env complex revealed that these characteristics are likely facilitated by a heavy-chain insertion and increased inter-protomer contacts. The ability of 1-18 to effectively restrict HIV-1 escape pathways provides a new option to successfully prevent and treat HIV-1 infection.

Protein tolerance to random circular permutation correlates with thermostability and local energetics of residue-residue contacts

Adenylate kinase (AK) orthologs with a range of thermostabilities were subjected to random circular permutation, and deep mutational scanning was used to evaluate where new protein termini were nondisruptive to activity. The fraction of circularly permuted variants that retained function in each library correlated with AK thermostability. In addition, analysis of the positional tolerance to new termini, which increase local conformational flexibility, showed that bonds were either functionally sensitive to cleavage across all homologs, differentially sensitive, or uniformly tolerant. The mobile AMP-binding domain, which displays the highest calculated contact energies, presented the greatest tolerance to new termini across all AKs. In contrast, retention of function in the lid and core domains was more dependent upon AK melting temperature. These results show that family permutation profiling identifies primary structure that has been selected by evolution for dynamics that are critical to activity within an enzyme family. These findings also illustrate how deep mutational scanning can be applied to protein homologs in parallel to differentiate how topology, stability, and local energetics govern mutational tolerance.

Deep Mutational Scanning Comprehensively Maps How Zika Envelope Protein Mutations Affect Viral Growth and Antibody Escape

Functional constraints on viral proteins are often assessed by examining sequence conservation among natural strains, but this approach is relatively ineffective for Zika virus because all known sequences are highly similar. Here, we take an alternative approach to map functional constraints on Zika virus's envelope (E) protein by using deep mutational scanning to measure how all amino acid mutations to the E protein affect viral growth in cell culture. The resulting sequence-function map is consistent with existing knowledge about E protein structure and function but also provides insight into mutation-level constraints in many regions of the protein that have not been well characterized in prior functional work. In addition, we extend our approach to completely map how mutations affect viral neutralization by two monoclonal antibodies, thereby precisely defining their functional epitopes. Overall, our study provides a valuable resource for understanding the effects of mutations to this important viral protein and also offers a roadmap for future work to map functional and antigenic selection to Zika virus at high resolution.IMPORTANCE Zika virus has recently been shown to be associated with severe birth defects. The virus's E protein mediates its ability to infect cells and is also the primary target of the antibodies that are elicited by natural infection and vaccines that are being developed against the virus. Therefore, determining the effects of mutations to this protein is important for understanding its function, its susceptibility to vaccine-mediated immunity, and its potential for future evolution. We completely mapped how amino acid mutations to the E protein affected the virus's ability to grow in cells in the laboratory and escape from several antibodies. The resulting maps relate changes in the E protein's sequence to changes in viral function and therefore provide a valuable complement to existing maps of the physical structure of the protein.