Mutational fitness landscapes reveal genetic and structural improvement pathways for a vaccine-elicited HIV-1 broadly neutralizing antibody

Optimized single-cell gates for yeast display screening

Yeast display is a widely used technology in antibody discovery and protein engineering. The cell size of yeast enables fluorescence-activated cell sorting (FACS) to precisely screen gene libraries, including for multi-parameter selection of protein phenotypes. However, yeast cells show a broader size distribution than mammalian cells that complicates single-cell gate determination for FACS. In this report, we analyze several yeast display gating options in detail and present an optimized strategy to select single yeast cells via flow cytometry. These data reveal optimized single-cell gating strategies to support robust and high-efficiency yeast display studies.

Clonal Variant Analysis of Antibody Engineering Libraries

In vitro antibody evolution is a powerful technique for improving monoclonal antibodies. This can be achieved by generating artificial diversity on an antibody template, which can be done using different in vitro diversification techniques. The resulting libraries consist of single- or multimutant variants of a defined antibody template that are screened for improved function using antibody display. Here, we describe a bioinformatic protocol for tracking synthetic antibody variants using high-throughput sequencing across screening rounds, enabling efficient high-throughput interpretation of the function of individual mutations in sorted antibody display libraries. The protocol enables a user to achieve precision analysis and interpretation of clonal antibody variant data for discovery purposes, especially for high-throughput antibody engineering or optimization against target antigens.

Antibody Data Analysis from Diverse Immune Libraries

Antibody functional screening studies and next-generation sequencing require careful processing and interpretation of sequence data for optimal results. Here, we provide a detailed protocol for the functional analysis of antibody gene data, including antibody repertoire quantification and functional mapping of high-throughput screening data based on enrichment ratio values, which are a simple way to determine the enrichment of each sequenced antibody after sorting a display library against desired antigens. This protocol enables a user to apply a set of simple yet powerful bioinformatic tools for high-throughput analysis and interpretation of antibody data that is especially well suited for display library screening and for antibody discovery applications.

Beyond Single Clones: High-Throughput Sequencing in Antibody Discovery

Antibody repertoire sequencing and display library screening are powerful approaches for antibody discovery and engineering that can connect DNA sequence with antibody function. Antibody display and screening studies have made a tremendous impact on immunology and biotechnology over the last decade, accelerated by technological advances in high-throughput DNA sequencing techniques. Indeed, bioinformatic analysis of antibody DNA library data has now taken a central role in modern antibody drug discovery, and is also critical for many ongoing studies of human immune development. Here, we describe current trends in antibody DNA library screening and analysis, and introduce a selection of protocols describing fundamental bioinformatic techniques to enable scientists to efficiently study antibody DNA libraries.

A library-on-library screen reveals the breadth expansion landscape of a broadly neutralizing betacoronavirus antibody

Broadly neutralizing antibodies (bnAbs) typically evolve cross-reactivity breadth through acquiring somatic hypermutations. While evolution of breadth requires improvement of binding to multiple antigenic variants, most experimental evolution platforms select against only one antigenic variant at a time. In this study, a yeast display library-on-library approach was applied to delineate the affinity maturation of a betacoronavirus bnAb, S2P6, against 27 spike stem helix peptides in a single experiment. Our results revealed that the binding affinity landscape of S2P6 varies among different stem helix peptides. However, somatic hypermutations that confer general improvement in binding affinity across different stem helix peptides could also be identified. We further showed that a key somatic hypermutation for breadth expansion involves long-range interaction. Overall, our work not only provides a proof-of-concept for using a library-on-library approach to analyze the evolution of antibody breadth, but also has important implications for the development of broadly protective vaccines.

Antibody-directed evolution reveals a mechanism for enhanced neutralization at the HIV-1 fusion peptide site

The HIV-1 fusion peptide (FP) represents a promising vaccine target, but global FP sequence diversity among circulating strains has limited anti-FP antibodies to ~60% neutralization breadth. Here we evolve the FP-targeting antibody VRC34.01 in vitro to enhance FP-neutralization using site saturation mutagenesis and yeast display. Successive rounds of directed evolution by iterative selection of antibodies for binding to resistant HIV-1 strains establish a variant, VRC34.01_mm28, as a best-in-class antibody with 10-fold enhanced potency compared to the template antibody and ~80% breadth on a cross-clade 208-strain neutralization panel. Structural analyses demonstrate that the improved paratope expands the FP binding groove to accommodate diverse FP sequences of different lengths while also recognizing the HIV-1 Env backbone. These data reveal critical antibody features for enhanced neutralization breadth and potency against the FP site of vulnerability and accelerate clinical development of broad HIV-1 FP-targeting vaccines and therapeutics.

Insights into next generation sequencing guided antibody selection strategies

Therapeutic antibody discovery often relies on in-vitro display methods to identify lead candidates. Assessing selected output diversity traditionally involves random colony picking and Sanger sequencing, which has limitations. Next-generation sequencing (NGS) offers a cost-effective solution with increased read depth, allowing a comprehensive understanding of diversity. Our study establishes NGS guidelines for antibody drug discovery, demonstrating its advantages in expanding the number of unique HCDR3 clusters, broadening the number of high affinity antibodies, expanding the total number of antibodies recognizing different epitopes, and improving lead prioritization. Surprisingly, our investigation into the correlation between NGS-derived frequencies of CDRs and affinity revealed a lack of association, although this limitation could be moderately mitigated by leveraging NGS clustering, enrichment and/or relative abundance across different regions to enhance lead prioritization. This study highlights NGS benefits, offering insights, recommendations, and the most effective approach to leverage NGS in therapeutic antibody discovery.

Large-scale antibody immune response mapping of splenic B cells and bone marrow plasma cells in a transgenic mouse model

Molecular characterization of antibody immunity and human antibody discovery is mainly carried out using peripheral memory B cells, and occasionally plasmablasts, that express B cell receptors (BCRs) on their cell surface. Despite the importance of plasma cells (PCs) as the dominant source of circulating antibodies in serum, PCs are rarely utilized because they do not express surface BCRs and cannot be analyzed using antigen-based fluorescence-activated cell sorting. Here, we studied the antibodies encoded by the entire mature B cell populations, including PCs, and compared the antibody repertoires of bone marrow and spleen compartments elicited by immunization in a human immunoglobulin transgenic mouse strain. To circumvent prior technical limitations for analysis of plasma cells, we applied single-cell antibody heavy and light chain gene capture from the entire mature B cell repertoires followed by yeast display functional analysis using a cytokine as a model immunogen. We performed affinity-based sorting of antibody yeast display libraries and large-scale next-generation sequencing analyses to follow antibody lineage performance, with experimental validation of 76 monoclonal antibodies against the cytokine antigen that identified three antibodies with exquisite double-digit picomolar binding affinity. We observed that spleen B cell populations generated higher affinity antibodies compared to bone marrow PCs and that antigen-specific splenic B cells had higher average levels of somatic hypermutation. A degree of clonal overlap was also observed between bone marrow and spleen antibody repertoires, indicating common origins of certain clones across lymphoid compartments. These data demonstrate a new capacity to functionally analyze antigen-specific B cell populations of different lymphoid organs, including PCs, for high-affinity antibody discovery and detailed fundamental studies of antibody immunity.

An optimized thermodynamics integration protocol for identifying beneficial mutations in antibody design

Accurate identification of beneficial mutations is central to antibody design. Many knowledge-based (KB) computational approaches have been developed to predict beneficial mutations, but their accuracy leaves room for improvement. Thermodynamic integration (TI) is an alchemical free energy algorithm that offers an alternative technique for identifying beneficial mutations, but its performance has not been evaluated. In this study, we developed an efficient TI protocol with high accuracy for predicting binding free energy changes of antibody mutations. The improved TI method outperforms KB methods at identifying both beneficial and deleterious mutations. We observed that KB methods have higher accuracies in predicting deleterious mutations than beneficial mutations. A pipeline using KB methods to efficiently exclude deleterious mutations and TI to accurately identify beneficial mutations was developed for high-throughput mutation scanning. The pipeline was applied to optimize the binding affinity of a broadly sarbecovirus neutralizing antibody 10-40 against the circulating severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) omicron variant. Three identified beneficial mutations show strong synergy and improve both binding affinity and neutralization potency of antibody 10-40. Molecular dynamics simulation revealed that the three mutations improve the binding affinity of antibody 10-40 through the stabilization of an altered binding mode with increased polar and hydrophobic interactions. Above all, this study presents an accurate and efficient TI-based approach for optimizing antibodies and other biomolecules.

The evolutionary and functional significance of germline immunoglobulin gene variation

The recombination between immunoglobulin (IG) gene segments determines an individual's naïve antibody repertoire and, consequently, (auto)antigen recognition. Emerging evidence suggests that mammalian IG germline variation impacts humoral immune responses associated with vaccination, infection, and autoimmunity - from the molecular level of epitope specificity, up to profound changes in the architecture of antibody repertoires. These links between IG germline variants and immunophenotype raise the question on the evolutionary causes and consequences of diversity within IG loci. We discuss why the extreme diversity in IG loci remains a mystery, why resolving this is important for the design of more effective vaccines and therapeutics, and how recent evidence from multiple lines of inquiry may help us do so.

Mapping monoclonal anti-SARS-CoV-2 antibody repertoires against diverse coronavirus antigens

Variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have emerged continuously, challenging the effectiveness of vaccines, diagnostics, and treatments. Moreover, the possibility of the appearance of a new betacoronavirus with high transmissibility and high fatality is reason for concern. In this study, we used a natively paired yeast display technology, combined with next-generation sequencing (NGS) and massive bioinformatic analysis to perform a comprehensive study of subdomain specificity of natural human antibodies from two convalescent donors. Using this screening technology, we mapped the cross-reactive responses of antibodies generated by the two donors against SARS-CoV-2 variants and other betacoronaviruses. We tested the neutralization potency of a set of the cross-reactive antibodies generated in this study and observed that most of the antibodies produced by these patients were non-neutralizing. We performed a comparison of the specific and non-specific antibodies by somatic hypermutation in a repertoire-scale for the two individuals and observed that the degree of somatic hypermutation was unique for each patient. The data from this study provide functional insights into cross-reactive antibodies that can assist in the development of strategies against emerging SARS-CoV-2 variants and divergent betacoronaviruses.

Highly protective antimalarial antibodies via precision library generation and yeast display screening

The monoclonal antibody CIS43 targets the Plasmodium falciparum circumsporozoite protein (PfCSP) and prevents malaria infection in humans for up to 9 mo following a single intravenous administration. To enhance the potency and clinical utility of CIS43, we used iterative site-saturation mutagenesis and DNA shuffling to screen precise gene-variant yeast display libraries for improved PfCSP antigen recognition. We identified several mutations that improved recognition, predominately in framework regions, and combined these to produce a panel of antibody variants. The most improved antibody, CIS43_Var10, had three mutations and showed approximately sixfold enhanced protective potency in vivo compared to CIS43. Co-crystal and cryo-electron microscopy structures of CIS43_Var10 with the peptide epitope or with PfCSP, respectively, revealed functional roles for each of these mutations. The unbiased site-directed mutagenesis and screening pipeline described here represent a powerful approach to enhance protective potency and to enable broader clinical use of antimalarial antibodies.

Biophysical Fitness Landscape of the SARS-CoV-2 Delta Variant Receptor Binding Domain

Among the five known SARS-CoV-2 variants of concern, Delta is the most virulent leading to severe symptoms and increased mortality among infected people. Our study seeks to examine how the biophysical parameters of the Delta variant correlate to the clinical observations. Receptor binding domain (RBD) is the first point of contact with the human host cells and is the immunodominant form of the spike protein. Delta variant RBD contains two novel mutations L452R and T478K. We examined the effect of single as well as the double mutations on RBD expression in human Expi293 cells, RBD stability using urea and thermal denaturation, and RBD binding to angiotensin converting enzyme 2 (ACE2) receptor and to neutralizing antibodies using isothermal titration calorimetry. Delta variant RBD showed significantly higher expression compared to the wild-type RBD, and the increased expression is due to L452R mutation. Despite their non-conservative nature, none of the mutations significantly affected RBD structure and stability. All mutants showed similar binding affinity to ACE2 and to Class 1 antibodies (CC12.1 and LY-CoV016) as that of the wild-type. Delta double mutant L452R/T478K showed no binding to Class 2 antibodies (P2B-2F6 and LY-CoV555) and a hundred-fold weaker binding to a Class 3 antibody (REGN10987), and the decreased antibody binding is determined by the L452R mutation. These results indicate that the immune escape from neutralizing antibodies, rather than increased receptor binding, is the main biophysical parameter that determined the fitness landscape of the Delta variant RBD.

Molecular probes of spike ectodomain and its subdomains for SARS-CoV-2 variants, Alpha through Omicron

Since the outbreak of the COVID-19 pandemic, widespread infections have allowed SARS-CoV-2 to evolve in human, leading to the emergence of multiple circulating variants. Some of these variants show increased resistance to vaccine-elicited immunity, convalescent plasma, or monoclonal antibodies. In particular, mutations in the SARS-CoV-2 spike have drawn attention. To facilitate the isolation of neutralizing antibodies and the monitoring of vaccine effectiveness against these variants, we designed and produced biotin-labeled molecular probes of variant SARS-CoV-2 spikes and their subdomains, using a structure-based construct design that incorporated an N-terminal purification tag, a specific amino acid sequence for protease cleavage, the variant spike-based region of interest, and a C-terminal sequence targeted by biotin ligase. These probes could be produced by a single step using in-process biotinylation and purification. We characterized the physical properties and antigenicity of these probes, comprising the N-terminal domain (NTD), the receptor-binding domain (RBD), the RBD and subdomain 1 (RBD-SD1), and the prefusion-stabilized spike ectodomain (S2P) with sequences from SARS-CoV-2 variants of concern or of interest, including variants Alpha, Beta, Gamma, Epsilon, Iota, Kappa, Delta, Lambda, Mu, and Omicron. We functionally validated probes by using yeast expressing a panel of nine SARS-CoV-2 spike-binding antibodies and confirmed sorting capabilities of variant probes using yeast displaying libraries of plasma antibodies from COVID-19 convalescent donors. We deposited these constructs to Addgene to enable their dissemination. Overall, this study describes a matrix of SARS-CoV-2 variant molecular probes that allow for assessment of immune responses, identification of serum antibody specificity, and isolation and characterization of neutralizing antibodies.

Idiosyncratic epistasis leads to global fitness-correlated trends

Epistasis can markedly affect evolutionary trajectories. In recent decades, protein-level fitness landscapes have revealed extensive idiosyncratic epistasis among specific mutations. By contrast, other work has found ubiquitous and apparently nonspecific patterns of global diminishing-returns and increasing-costs epistasis among mutations across the genome. Here, we used a hierarchical CRISPR gene drive system to construct all combinations of 10 missense mutations from across the genome in budding yeast and measured their fitness in six environments. We show that the resulting fitness landscapes exhibit global fitness-correlated trends but that these trends emerge from specific idiosyncratic interactions. We thus provide experimental validation of recent theoretical work arguing that fitness-correlated trends can emerge as the generic consequence of idiosyncratic epistasis.

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.

Immortalization and functional screening of natively paired human T cell receptor repertoires

Functional analyses of the T cell receptor (TCR) landscape can reveal critical information about protection from disease and molecular responses to vaccines. However, it has proven difficult to combine advanced next-generation sequencing technologies with methods to decode the peptide-major histocompatibility complex (pMHC) specificity of individual TCRs. We developed a new high-throughput approach to enable repertoire-scale functional evaluations of natively paired TCRs. In particular, we leveraged the immortalized nature of physically linked TCRα:β amplicon libraries to analyze binding against multiple recombinant pMHCs on a repertoire scale, and to exemplify the utility of this approach, we also performed affinity-based functional mapping in conjunction with quantitative next-generation sequencing to track antigen-specific TCRs. These data successfully validated a new immortalization and screening platform to facilitate detailed molecular analyses of disease-relevant antigen interactions with human TCRs.

Molecular probes of spike ectodomain and its subdomains for SARS-CoV-2 variants, Alpha through Omicron

Since the outbreak of the COVID-19 pandemic, widespread infections have allowed SARS-CoV-2 to evolve in human, leading to the emergence of multiple circulating variants. Some of these variants show increased resistance to vaccines, convalescent plasma, or monoclonal antibodies. In particular, mutations in the SARS-CoV-2 spike have drawn attention. To facilitate the isolation of neutralizing antibodies and the monitoring the vaccine effectiveness against these variants, we designed and produced biotin-labeled molecular probes of variant SARS-CoV-2 spikes and their subdomains, using a structure-based construct design that incorporated an N-terminal purification tag, a specific amino acid sequence for protease cleavage, the variant spike-based region of interest, and a C-terminal sequence targeted by biotin ligase. These probes could be produced by a single step using in-process biotinylation and purification. We characterized the physical properties and antigenicity of these probes, comprising the N-terminal domain (NTD), the receptor-binding domain (RBD), the RBD and subdomain 1 (RBD-SD1), and the prefusion-stabilized spike ectodomain (S2P) with sequences from SARS-CoV-2 variants of concern or of interest, including variants Alpha, Beta, Gamma, Epsilon, Iota, Kappa, Delta, Lambda, Mu, and Omicron. We functionally validated probes by using yeast expressing a panel of nine SARS-CoV-2 spike-binding antibodies and confirmed sorting capabilities of variant probes using yeast displaying libraries of plasma antibodies from COVID-19 convalescent donors. We deposited these constructs to Addgene to enable their dissemination. Overall, this study describes a matrix of SARS-CoV-2 variant molecular probes that allow for assessment of immune responses, identification of serum antibody specificity, and isolation and characterization of neutralizing antibodies.

Molecular probes of spike ectodomain and its subdomains for SARS-CoV-2 variants, Alpha through Omicron

Since the outbreak of the COVID-19 pandemic, widespread infections have allowed SARS-CoV-2 to evolve in human, leading to the emergence of multiple circulating variants. Some of these variants show increased resistance to vaccines, convalescent plasma, or monoclonal antibodies. In particular, mutations in the SARS-CoV-2 spike have drawn attention. To facilitate the isolation of neutralizing antibodies and the monitoring the vaccine effectiveness against these variants, we designed and produced biotin-labeled molecular probes of variant SARS-CoV-2 spikes and their subdomains, using a structure-based construct design that incorporated an N-terminal purification tag, a specific amino acid sequence for protease cleavage, the variant spike-based region of interest, and a C-terminal sequence targeted by biotin ligase. These probes could be produced by a single step using in-process biotinylation and purification. We characterized the physical properties and antigenicity of these probes, comprising the N-terminal domain (NTD), the receptor-binding domain (RBD), the RBD and subdomain 1 (RBD-SD1), and the prefusion-stabilized spike ectodomain (S2P) with sequences from SARS-CoV-2 variants of concern or of interest, including variants Alpha, Beta, Gamma, Epsilon, Iota, Kappa, Delta, Lambda, Mu, and Omicron. We functionally validated probes by using yeast expressing a panel of nine SARS-CoV-2 spike-binding antibodies and confirmed sorting capabilities of variant probes using yeast displaying libraries of plasma antibodies from COVID-19 convalescent donors. We deposited these constructs to Addgene to enable their dissemination. Overall, this study describes a matrix of SARS-CoV-2 variant molecular probes that allow for assessment of immune responses, identification of serum antibody specificity, and isolation and characterization of neutralizing antibodies.

Antibody screening at reduced pH enables preferential selection of potently neutralizing antibodies targeting SARS-CoV-2

Antiviral monoclonal antibody (mAb) discovery enables the development of antibody-based antiviral therapeutics. Traditional antiviral mAb discovery relies on affinity between antibody and a viral antigen to discover potent neutralizing antibodies, but these approaches are inefficient because many high affinity mAbs have no neutralizing activity. We sought to determine whether screening for anti-SARS-CoV-2 mAbs at reduced pH could provide more efficient neutralizing antibody discovery. We mined the antibody response of a convalescent COVID-19 patient at both physiological pH (7.4) and reduced pH (4.5), revealing that SARS-CoV-2 neutralizing antibodies were preferentially enriched in pH 4.5 yeast display sorts. Structural analysis revealed that a potent new antibody called LP5 targets the SARS-CoV-2 N-terminal domain supersite via a unique binding recognition mode. Our data combine with evidence from prior studies to support antibody screening at pH 4.5 to accelerate antiviral neutralizing antibody discovery.

Binding affinity landscapes constrain the evolution of broadly neutralizing anti-influenza antibodies

Over the past two decades, several broadly neutralizing antibodies (bnAbs) that confer protection against diverse influenza strains have been isolated. Structural and biochemical characterization of these bnAbs has provided molecular insight into how they bind distinct antigens. However, our understanding of the evolutionary pathways leading to bnAbs, and thus how best to elicit them, remains limited. Here, we measure equilibrium dissociation constants of combinatorially complete mutational libraries for two naturally isolated influenza bnAbs (CR9114, 16 heavy-chain mutations; CR6261, 11 heavy-chain mutations), reconstructing all possible evolutionary intermediates back to the unmutated germline sequences. We find that these two libraries exhibit strikingly different patterns of breadth: while many variants of CR6261 display moderate affinity to diverse antigens, those of CR9114 display appreciable affinity only in specific, nested combinations. By examining the extensive pairwise and higher order epistasis between mutations, we find key sites with strong synergistic interactions that are highly similar across antigens for CR6261 and different for CR9114. Together, these features of the binding affinity landscapes strongly favor sequential acquisition of affinity to diverse antigens for CR9114, while the acquisition of breadth to more similar antigens for CR6261 is less constrained. These results, if generalizable to other bnAbs, may explain the molecular basis for the widespread observation that sequential exposure favors greater breadth, and such mechanistic insight will be essential for predicting and eliciting broadly protective immune responses.