A pseudovirus system enables deep mutational scanning of the full SARS-CoV-2 spike.

A major challenge in understanding SARS-CoV-2 evolution is interpreting the antigenic and functional effects of emerging mutations in the viral spike protein. Here, we describe a deep mutational scanning platform based on non-replicative pseudotyped lentiviruses that directly quantifies how large numbers of spike mutations impact antibody neutralization and pseudovirus infection. We apply this platform to produce libraries of the Omicron BA.1 and Delta spikes. These libraries each contain ∼7,000 distinct amino acid mutations in the context of up to ∼135,000 unique mutation combinations. We use these libraries to map escape mutations from neutralizing antibodies targeting the receptor-binding domain, N-terminal domain, and S2 subunit of spike. Overall, this work establishes a high-throughput and safe approach to measure how ∼105 combinations of mutations affect antibody neutralization and spike-mediated infection. Notably, the platform described here can be extended to the entry proteins of many other viruses.

Functional and antigenic landscape of the Nipah virus receptor binding protein.

Nipah virus recurrently spills over to humans, causing fatal infections. The viral receptor-binding protein (RBP or G) attaches to host receptors and is a major target of neutralizing antibodies. Here we use deep mutational scanning to measure how all amino-acid mutations to the RBP affect cell entry, receptor binding, and escape from neutralizing antibodies. We identify functionally constrained regions of the RBP, including sites involved in oligomerization, along with mutations that differentially modulate RBP binding to its two ephrin receptors. We map escape mutations for six anti-RBP antibodies, and find that few antigenic mutations are present in natural Nipah strains. Our findings offer insights into the potential for functional and antigenic evolution of the RBP that can inform the development of antibody therapies and vaccines.

Mapping the neutralizing specificity of human anti-HIV serum by deep mutational scanning.

Understanding the specificities of human serum antibodies that broadly neutralize HIV can inform prevention and treatment strategies. Here we describe a deep mutational scanning system that can measure the effects of combinations of mutations to HIV envelope (Env) on neutralization by antibodies and polyclonal serum. We first show that this system can accurately map how all functionally tolerated mutations to Env affect neutralization by monoclonal antibodies. We then comprehensively map Env mutations that affect neutralization by a set of human polyclonal sera known to target the CD4-binding site that neutralize diverse strains of HIV. The neutralizing activities of these sera target different epitopes, with most sera having specificities reminiscent of individual characterized monoclonal antibodies, but one sera targeting two epitopes within the CD4 binding site. Mapping the specificity of the neutralizing activity in polyclonal human serum will aid in assessing anti-HIV immune responses to inform prevention strategies.

Mapping the neutralizing specificity of human anti-HIV serum by deep mutational scanning.

Understanding the specificities of human serum antibodies that broadly neutralize HIV can inform prevention and treatment strategies. Here, we describe a deep mutational scanning system that can measure the effects of combinations of mutations to HIV envelope (Env) on neutralization by antibodies and polyclonal serum. We first show that this system can accurately map how all functionally tolerated mutations to Env affect neutralization by monoclonal antibodies. We then comprehensively map Env mutations that affect neutralization by a set of human polyclonal sera that neutralize diverse strains of HIV and target the site engaging the host receptor CD4. The neutralizing activities of these sera target different epitopes, with most sera having specificities reminiscent of individual characterized monoclonal antibodies, but one serum targeting two epitopes within the CD4-binding site. Mapping the specificity of the neutralizing activity in polyclonal human serum will aid in assessing anti-HIV immune responses to inform prevention strategies.

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.

A biophysical model of viral escape from polyclonal antibodies.

A challenge in studying viral immune escape is determining how mutations combine to escape polyclonal antibodies, which can potentially target multiple distinct viral epitopes. Here we introduce a biophysical model of this process that partitions the total polyclonal antibody activity by epitope and then quantifies how each viral mutation affects the antibody activity against each epitope. We develop software that can use deep mutational scanning data to infer these properties for polyclonal antibody mixtures. We validate this software using a computationally simulated deep mutational scanning experiment and demonstrate that it enables the prediction of escape by arbitrary combinations of mutations. The software described in this paper is available at https://jbloomlab.github.io/polyclonal.

A pseudovirus system enables deep mutational scanning of the full SARS-CoV-2 spike.

A major challenge in understanding SARS-CoV-2 evolution is interpreting the antigenic and functional effects of emerging mutations in the viral spike protein. Here we describe a new deep mutational scanning platform based on non-replicative pseudotyped lentiviruses that directly quantifies how large numbers of spike mutations impact antibody neutralization and pseudovirus infection. We demonstrate this new platform by making libraries of the Omicron BA.1 and Delta spikes. These libraries each contain ~7000 distinct amino-acid mutations in the context of up to ~135,000 unique mutation combinations. We use these libraries to map escape mutations from neutralizing antibodies targeting the receptor binding domain, N-terminal domain, and S2 subunit of spike. Overall, this work establishes a high-throughput and safe approach to measure how ~10 5 combinations of mutations affect antibody neutralization and spike-mediated infection. Notably, the platform described here can be extended to the entry proteins of many other viruses.

Identification of coagulation factor IX variants with enhanced activity through ancestral sequence reconstruction.

Orthologous proteins contain sequence disparity guided by natural selection. In certain cases, species-specific protein functionality predicts pharmacological enhancement, such as greater specific activity or stability. However, immunological barriers generally preclude use of nonhuman proteins as therapeutics, and difficulty exists in the identification of individual sequence determinants among the overall sequence disparity. Ancestral sequence reconstruction (ASR) represents a platform for the prediction and resurrection of ancient gene and protein sequences. Recently, we demonstrated that ASR can be used as a platform to facilitate the identification of therapeutic protein variants with enhanced properties. Specifically, we identified coagulation factor VIII (FVIII) variants with improved specific activity, biosynthesis, stability, and resistance to anti-human FVIII antibody-based inhibition. In the current study, we resurrected a panel of ancient mammalian coagulation factor IX (FIX) variants with the goal of identifying improved pharmaceutical candidates. One variant (An96) demonstrated 12-fold greater FIX activity production than human FIX. Addition of the R338L Padua substitution further increased An96 activity, suggesting independent but additive mechanisms. after adeno-associated virus 2 (AAV2)/8-FIX gene therapy, 10-fold greater plasma FIX activity was observed in hemophilia B mice administered AAV2/8-An96-Padua as compared with AAV2/8-human FIX-Padua. Furthermore, phenotypic correction conferred by the ancestral variant was confirmed using a saphenous vein bleeding challenge and thromboelastography. Collectively, these findings validate the ASR drug discovery platform as well as identify an ancient FIX candidate for pharmaceutical development.

Molecular coevolution of coagulation factor VIII and von Willebrand factor.

Ancestral sequence reconstruction provides a unique platform for investigating the molecular evolution of single gene products and recently has shown success in engineering advanced biological therapeutics. To date, the coevolution of proteins within complexes and protein-protein interactions is mostly investigated in silico via proteomics and/or within single-celled systems. Herein, ancestral sequence reconstruction is used to investigate the molecular evolution of 2 proteins linked not only by stabilizing association in circulation but also by their independent roles within the primary and secondary hemostatic systems of mammals. Using sequence analysis and biochemical characterization of recombinant ancestral von Willebrand factor (VWF) and coagulation factor VIII (FVIII), we investigated the evolution of the essential macromolecular FVIII/VWF complex. Our data support the hypothesis that these coagulation proteins coevolved throughout mammalian diversification, maintaining strong binding affinities while modulating independent and distinct hemostatic activities in diverse lineages.

An experimental phylogeny to benchmark ancestral sequence reconstruction.

Ancestral sequence reconstruction (ASR) is a still-burgeoning method that has revealed many key mechanisms of molecular evolution. One criticism of the approach is an inability to validate its algorithms within a biological context as opposed to a computer simulation. Here we build an experimental phylogeny using the gene of a single red fluorescent protein to address this criticism. The evolved phylogeny consists of 19 operational taxonomic units (leaves) and 17 ancestral bifurcations (nodes) that display a wide variety of fluorescent phenotypes. The 19 leaves then serve as 'modern' sequences that we subject to ASR analyses using various algorithms and to benchmark against the known ancestral genotypes and ancestral phenotypes. We confirm computer simulations that show all algorithms infer ancient sequences with high accuracy, yet we also reveal wide variation in the phenotypes encoded by incorrectly inferred sequences. Specifically, Bayesian methods incorporating rate variation significantly outperform the maximum parsimony criterion in phenotypic accuracy. Subsampling of extant sequences had minor effect on the inference of ancestral sequences.