Assessment of Solvated Interaction Energy Function for Ranking Antibody-Antigen Binding Affinities

Improved Structure-Based Histidine pKa Prediction for pH-Responsive Protein Design

The near neutral pKa of histidine is commonly exploited to engineer pH-sensitive biomolecules. For example, histidine mutations introduced in the complementarity-determining region (CDR) of therapeutic antibodies can enhance selectivity for antigens in the acidic microenvironment of solid tumors or increase dissociation rates in the acidic early endosomes of cells. While solvent-exposed histidines typically have a pKa near 6.5, interacting histidines can experience pKa shifts of up to 4 pH units in either direction, making histidine one of the most variable titratable residues. To assist in selecting potential histidine mutation sites, pKa prediction software should achieve an accuracy significantly better than the current standard of around 1.0 pH unit. However, the limited availability of experimental histidine pKa measurements hinders the use of AI-based methods. This study evaluates histidine pKa predictions using Amber force field electrostatics combined with a continuum solvent model, previously calibrated in the solvated interaction energy (SIE) function for binding affinity predictions. By incorporating limited rotameric sampling, proton optimization, and an empirical correction for buried side-chains, the method achieves a mean unsigned error of 0.4 pH units across a diverse set of 91 histidines from 38 distinct protein structures obtained from the PKAD database. This approach should improve the in-silico design of pH-responsive mutations. The method is implemented in the software program JustHISpKa (https://mm.nrc-cnrc.gc.ca/software/JustHISpKa).

A benchmark for evaluation of structure-based online tools for antibody-antigen binding affinity

The prediction of binding affinity changes caused by missense mutations can elucidate antigen-antibody interactions. A few accessible structure-based online computational tools have been proposed. However, selecting suitable software for particular research is challenging, especially research on the SARS-CoV-2 spike protein with antibodies. Therefore, benchmarking of the mutation-diverse SARS-CoV-2 datasets is critical. Here, we collected the datasets including 1216 variants about the changes in binding affinity of antigens from 22 complexes for SARS-CoV-2 S proteins and 22 monoclonal antibodies as well as applied them to evaluate the performance of seven binding affinity prediction tools. The tested tools' Pearson correlations between predicted and measured changes in binding affinity were between -0.158 and 0.657, while accuracy in classification tasks on predicting increasing or decreasing affinity ranged from 0.444 to 0.834. These tools performed relatively better on predicting single mutations, especially at epitope sites, whereas poor performance on extremely decreasing affinity. The tested tools were relatively insensitive to the experimental techniques used to obtain structures of complexes. In summary, we constructed a list of datasets and evaluated a range of structure-based online prediction tools that will explicate relevant processes of antigen-antibody interactions and enhance the computational design of therapeutic monoclonal antibodies.

Sequence-based engineering of pH-sensitive antibodies for tumor targeting or endosomal recycling applications

The engineering of pH-sensitive therapeutic antibodies, particularly for improving effectiveness and specificity in acidic solid-tumor microenvironments, has recently gained traction. While there is a justified need for pH-dependent immunotherapies, current engineering techniques are tedious and laborious, requiring repeated rounds of experiments under different pH conditions. Inexpensive computational techniques to predict the effectiveness of His pH-switches require antibody-antigen complex structures, but these are lacking in most cases. To circumvent these requirements, we introduce a sequence-based in silico method for predicting His mutations in the variable region of antibodies, which could lead to pH-biased antigen binding. This method, called Sequence-based Identification of pH-sensitive Antibody Binding (SIpHAB), was trained on 3D-structure-based calculations of 3,490 antibody-antigen complexes with solved experimental structures. SIpHAB was parametrized to enhance preferential binding either toward or against the acidic pH, for selective targeting of solid tumors or for antigen release in the endosome, respectively. Applications to nine antibody-antigen systems with previously reported binding preferences at different pHs demonstrated the utility and enrichment capabilities of this high-throughput computational tool. SIpHAB, which only requires knowledge of the antibody primary amino-acid sequence, could enable a more efficient triage of pH-sensitive antibody candidates than could be achieved conventionally. An online webserver for running SipHAB is available freely at https://mm.nrc-cnrc.gc.ca/software/siphab/runner/.

Solvated interaction energy: from small-molecule to antibody drug design

Scoring functions are ubiquitous in structure-based drug design as an aid to predicting binding modes and estimating binding affinities. Ideally, a scoring function should be broadly applicable, obviating the need to recalibrate and refit its parameters for every new target and class of ligands. Traditionally, drugs have been small molecules, but in recent years biologics, particularly antibodies, have become an increasingly important if not dominant class of therapeutics. This makes the goal of having a transferable scoring function, i.e., one that spans the range of small-molecule to protein ligands, even more challenging. One such broadly applicable scoring function is the Solvated Interaction Energy (SIE), which has been developed and applied in our lab for the last 15 years, leading to several important applications. This physics-based method arose from efforts to understand the physics governing binding events, with particular care given to the role played by solvation. SIE has been used by us and many independent labs worldwide for virtual screening and discovery of novel small-molecule binders or optimization of known drugs. Moreover, without any retraining, it is found to be transferrable to predictions of antibody-antigen relative binding affinities and as accurate as functions trained on protein-protein binding affinities. SIE has been incorporated in conjunction with other scoring functions into ADAPT (Assisted Design of Antibody and Protein Therapeutics), our platform for affinity modulation of antibodies. Application of ADAPT resulted in the optimization of several antibodies with 10-to-100-fold improvements in binding affinity. Further applications included broadening the specificity of a single-domain antibody to be cross-reactive with virus variants of both SARS-CoV-1 and SARS-CoV-2, and the design of safer antibodies by engineering of a pH switch to make them more selective towards acidic tumors while sparing normal tissues at physiological pH.

In Silico Analysis of a GH3 β-Glucosidase from Microcystis aeruginosa CACIAM 03

Cyanobacteria are rich sources of secondary metabolites and have the potential to be excellent industrial enzyme producers. β-glucosidases are extensively employed in processing biomass degradation as they mediate the most crucial step of bioconversion of cellobiose (CBI), hence controlling the efficiency and global rate of biomass hydrolysis. However, the production and availability of these enzymes derived from cyanobacteria remains limited. In this study, we evaluated the β-glucosidase from Microcystis aeruginosa CACIAM 03 (MaBgl3) and its potential for bioconversion of cellulosic biomass by analyzing primary/secondary structures, predicting physicochemical properties, homology modeling, molecular docking, and simulations of molecular dynamics (MD). The results showed that MaBgl3 derives from an N-terminal domain folded as a distorted β-barrel, which contains the conserved His-Asp catalytic dyad often found in glycosylases of the GH3 family. The molecular docking results showed relevant interactions with Asp81, Ala271 and Arg444 residues that contribute to the binding process during MD simulation. Moreover, the MD simulation of the MaBgl3 was stable, shown by analyzing the root mean square deviation (RMSD) values and observing favorable binding free energy in both complexes. In addition, experimental data suggest that MaBgl3 could be a potential enzyme for cellobiose-hydrolyzing degradation.

Binding affinity prediction for antibody-protein antigen complexes: A machine learning analysis based on interface and surface areas

Specific antibodies can bind to protein antigens with high affinity and specificity, and this property makes them one of the best protein-based therapeutics. Accurate prediction of antibody‒protein antigen binding affinity is crucial for designing effective antibodies. The current predictive methods for protein‒protein binding affinity usually fail to predict the binding affinity of an antibody‒protein antigen complex with a comparable level of accuracy. Here, new models specific for antibody‒antigen binding affinity prediction are developed according to the different types of interface and surface areas present in antibody‒antigen complex. The contacts-based descriptors are also employed to construct or train different models specific for antibody‒protein antigen binding affinity prediction. The results of this study show that (i) the area-based descriptors are slightly better than the contacts-based descriptors in terms of the predictive power; (ii) the new models specific for antibody‒protein antigen binding affinity prediction are superior to the previously-used general models for predicting the protein‒protein binding affinities; (iii) the performances of the best area-based and contacts-based models developed in this work are better than the performances of a recently-developed graph-based model (i.e., CSM-AB) specific for antibody‒protein antigen binding affinity prediction. The new models developed in this work would not only help understand the mechanisms underlying antibody‒protein antigen interactions, but would also be of some applicable utility in the design and virtual screening of antibody-based therapeutics.

Optimizing Antibody-Antigen Binding Affinities with the ADAPT Platform

The ADAPT (Assisted Design of Antibody and Protein Therapeutics) platform guides the selection of mutants that improve/modulate the affinity of antibodies and other biologics. Predicted affinities are based on a consensus z-score from three scoring functions. Computational predictions are interleaved with experimental validation, significantly enhancing the robustness of the design and selection of mutants. A key step is an initial exhaustive virtual single-mutant scan that identifies hot spots and the mutations predicted to improve affinity. A small number of proposed single mutants are then produced and assayed. Only the validated single mutants (i.e., having improved affinity) are used to design double and higher-order mutants in subsequent rounds of design, avoiding the combinatorial explosion that arises from random mutagenesis. Typically, with a total of about 30-50 designed single, double, and triple mutants, affinity improvements of 10- to 100-fold are obtained.

Structure-based dual affinity optimization of a SARS-CoV-1/2 cross-reactive single-domain antibody

The SARS coronavirus 2 (SARS-CoV-2) spike (S) protein binding to the human ACE2 receptor is the molecular event that initiates viral entry into host cells and leads to infection and virus replication. There is a need for agents blocking viral entry into host cells that are cross-reactive with emerging virus variants. VHH-72 is an anti-SARS-CoV-1 single-domain antibody that also exhibits cross-specificity with SARS-CoV-2 but with decreased binding affinity. Here we applied a structure-based approach to affinity-mature VHH-72 for the SARS-CoV-2 spike protein while retaining the original affinity for SARS-CoV-1. This was achieved by employing the computational platform ADAPT in a constrained dual-affinity optimization mode as a means of broadening specificity. Select mutants designed by ADAPT were formatted as fusions with a human IgG1-Fc fragment. These mutants demonstrated improved binding to the SARS-CoV-2 spike protein due to decreased dissociation rates. Functional testing for virus neutralization revealed improvements relative to the parental VHH72-Fc up to 10-fold using a SARS-CoV-2 pseudotyped lentivirus and 20-fold against the SARS-CoV-2 authentic live virus (Wuhan variant). Binding and neutralization improvements were maintained for some other SARS-CoV-2 variants currently in circulation. These improved VHH-72 mutants are predicted to establish novel interactions with the S antigen. They will be useful, alone or as fusions with other functional modules, in the global quest for treatments of COVID-19 infections.

Progress and challenges for the machine learning-based design of fit-for-purpose monoclonal antibodies

Although the therapeutic efficacy and commercial success of monoclonal antibodies (mAbs) are tremendous, the design and discovery of new candidates remain a time and cost-intensive endeavor. In this regard, progress in the generation of data describing antigen binding and developability, computational methodology, and artificial intelligence may pave the way for a new era of in silico on-demand immunotherapeutics design and discovery. Here, we argue that the main necessary machine learning (ML) components for an in silico mAb sequence generator are: understanding of the rules of mAb-antigen binding, capacity to modularly combine mAb design parameters, and algorithms for unconstrained parameter-driven in silico mAb sequence synthesis. We review the current progress toward the realization of these necessary components and discuss the challenges that must be overcome to allow the on-demand ML-based discovery and design of fit-for-purpose mAb therapeutic candidates.

Redesigning an antibody H3 loop by virtual screening of a small library of human germline-derived sequences

The design of superior biologic therapeutics, including antibodies and engineered proteins, involves optimizing their specific ability to bind to disease-related molecular targets. Previously, we developed and applied the Assisted Design of Antibody and Protein Therapeutics (ADAPT) platform for virtual affinity maturation of antibodies (Vivcharuk et al. in PLoS One 12(7):e0181490, 10.1371/journal.pone.0181490, 2017). However, ADAPT is limited to point mutations of hot-spot residues in existing CDR loops. In this study, we explore the possibility of wholesale replacement of the entire H3 loop with no restriction to maintain the parental loop length. This complements other currently published studies that sample replacements for the CDR loops L1, L2, L3, H1 and H2. Given the immense sequence space theoretically available to H3, we focused on the virtual grafting of over 5000 human germline-derived H3 sequences from the IGMT/LIGM database increasing the diversity of the sequence space when compared to using crystalized H3 loop sequences. H3 loop conformations are generated and scored to identify optimized H3 sequences. Experimental testing of high-ranking H3 sequences grafted into the framework of the bH1 antibody against human VEGF-A led to the discovery of multiple hits, some of which had similar or better affinities relative to the parental antibody. In over 75% of the tested designs, the re-designed H3 loop contributed favorably to overall binding affinity. The hits also demonstrated good developability attributes such as high thermal stability and no aggregation. Crystal structures of select re-designed H3 variants were solved and indicated that although some deviations from predicted structures were seen in the more solvent accessible regions of the H3 loop, they did not significantly affect predicted affinity scores.

Application of an integrated computational antibody engineering platform to design SARS-CoV-2 neutralizers

As the COVID-19 pandemic continues to spread, hundreds of new initiatives including studies on existing medicines are running to fight the disease. To deliver a potentially immediate and lasting treatment to current and emerging SARS-CoV-2 variants, new collaborations and ways of sharing are required to create as many paths forward as possible. Here, we leverage our expertise in computational antibody engineering to rationally design/engineer three previously reported SARS-CoV neutralizing antibodies and share our proposal towards anti-SARS-CoV-2 biologics therapeutics. SARS-CoV neutralizing antibodies, m396, 80R and CR-3022 were chosen as templates due to their diversified epitopes and confirmed neutralization potency against SARS-CoV (but not SARS-CoV-2 except for CR3022). Structures of variable fragment (Fv) in complex with receptor binding domain (RBD) from SARS-CoV or SARS-CoV-2 were subjected to our established in silico antibody engineering platform to improve their binding affinity to SARS-CoV-2 and developability profiles. The selected top mutations were ensembled into a focused library for each antibody for further screening. In addition, we convert the selected binders with different epitopes into the trispecific format, aiming to increase potency and to prevent mutational escape. Lastly, to avoid antibody-induced virus activation or enhancement, we suggest application of NNAS and DQ mutations to the Fc region to eliminate effector functions and extend half-life.

An expanded benchmark for antibody-antigen docking and affinity prediction reveals insights into antibody recognition determinants

Accurate predictive modeling of antibody-antigen complex structures and structure-based antibody design remain major challenges in computational biology, with implications for biotherapeutics, immunity, and vaccines. Through a systematic search for high-resolution structures of antibody-antigen complexes and unbound antibody and antigen structures, in conjunction with identification of experimentally determined binding affinities, we have assembled a non-redundant set of test cases for antibody-antigen docking and affinity prediction. This benchmark more than doubles the number of antibody-antigen complexes and corresponding affinities available in our previous benchmarks, providing an unprecedented view of the determinants of antibody recognition and insights into molecular flexibility. Initial assessments of docking and affinity prediction tools highlight the challenges posed by this diverse set of cases, which includes camelid nanobodies, therapeutic monoclonal antibodies, and broadly neutralizing antibodies targeting viral glycoproteins. This dataset will enable development of advanced predictive modeling and design methods for this therapeutically relevant class of protein-protein interactions.

Structure-based engineering of pH-dependent antibody binding for selective targeting of solid-tumor microenvironment

Recent development of monoclonal antibodies as mainstream anticancer agents demands further optimization of their safety for use in humans. Potent targeting and/or effector activities on normal tissues is an obvious toxicity concern. Optimization of specific tumor targeting could be achieved by taking advantage of the extracellular acidity of solid tumors relative to normal tissues. Here, we applied a structure-based computational approach to engineer anti-human epidermal growth factor receptor 2 (Her2) antibodies with selective binding in the acidic tumor microenvironment. We used an affinity maturation platform in which dual-pH histidine-scanning mutagenesis was implemented for pH selectivity optimization. Testing of a small set of designs for binding to the recombinant Her2 ectodomain led to the identification of antigen-binding fragment (Fab) variants with the desired pH-dependent binding behavior. Binding selectivity toward acidic pH was improved by as much as 25-fold relative to the parental bH1-Fab. In vitro experiments on cells expressing intact Her2 confirmed that designed variants formatted as IgG1/k full-size antibodies have high affinity and inhibit the growth of tumor spheroids at a level comparable to that of the benchmark anti-Her2 antibody trastuzumab (Herceptin®) at acidic pH, whereas these effects were significantly reduced at physiological pH. In contrast, both Herceptin and the parental bH1 antibody exhibited strong cell binding and growth inhibition irrespective of pH. This work demonstrates the feasibility of computational optimization of antibodies for selective targeting of the acidic environment such as that found in many solid tumors.

Predicting antibody affinity changes upon mutations by combining multiple predictors

Antibodies are proteins working in our immune system with high affinity and specificity for target antigens, making them excellent tools for both biotherapeutic and bioengineering applications. The prediction of antibody affinity changes upon mutations (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\Delta \Delta {\mathrm{G}}}}_{\mathrm{binding}}$$\end{document}ΔΔGbinding) is important for antibody engineering. Numerous computational methods have been proposed based on different approaches including molecular mechanics and machine learning. However, the accuracy by each individual predictor is not enough for efficient antibody development. In this study, we develop a new prediction method by combining multiple predictors based on machine learning. Our method was tested on the SiPMAB database, evaluating the Pearson’s correlation coefficient between predicted and experimental \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\Delta \Delta {\mathrm{G}}}}_{\mathrm{binding}}$$\end{document}ΔΔGbinding. Our method achieved higher accuracy (R = 0.69) than previous molecular mechanics or machine-learning based methods (R = 0.59) and the previous method using the average of multiple predictors (R = 0.64). Feature importance analysis indicated that the improved accuracy was obtained by combining predictors with different importance, which have different protocols for calculating energies and for generating mutant and unbound state structures. This study demonstrates that machine learning is a powerful framework for combining different approaches to predict antibody affinity changes.

Binding and functional profiling of antibody mutants guides selection of optimal candidates as antibody drug conjugates

An increasingly appreciated conundrum in the discovery of antibody drug conjugates (ADCs) is that an antibody that was selected primarily for strong binding to its cancer target may not serve as an optimal ADC. In this study, we performed mechanistic cell-based experiments to determine the correlation between antibody affinity, avidity, internalization and ADC efficacy. We used structure-guided design to assemble a panel of antibody mutants with predicted Her2 affinities ranging from higher to lower relative to the parent antibody, Herceptin. These antibodies were ranked for binding via SPR and via flow-cytometry on high-Her2 SKOV3 cells and low-Her2 MCF7 cells, the latter acting as a surrogate for low-Her2 normal cells. A subpanel of variants, representative of different Her2-binding affinities (2 strong, 2 moderate and 3 weak), were further screened via high-content imaging for internalization efficacies in high versus low-Her2 cells. Finally, these antibodies were evaluated in ADC cytotoxicity screening assays (using DM1 and MMAE secondary antibodies) and as antibody-drug conjugates (DM1 and PNU159682). Our results identified specific but weak Her2-binding variants as optimal candidates for developing DM1 and PNU ADCs since they exhibited high potencies (low to sub-nM) in high-Her2 SKOV3 cells and low toxicities in low-Her2 cells. The 2 strong-affinity variants were highly potent in SKOV3 cells but also showed significant toxicities in low-Her2 cells and therefore are predicted to be toxic in normal tissues. Our findings show that pharmacological profiling of an antibody library in multiple binding and functional assays allows for selection of optimal ADCs.

Binding symmetry and surface flexibility mediate antibody self-association

Solution stability is an important factor in the optimization of engineered biotherapeutic candidates such as monoclonal antibodies because of its possible effects on manufacturability, pharmacology, efficacy and safety. A detailed atomic understanding of the mechanisms governing self-association of natively folded protein monomers is required to devise predictive tools to guide screening and re-engineering along the drug development pipeline. We investigated pairs of affinity-matured full-size antibodies and observed drastically different propensities to aggregate from variants differing by a single amino-acid. Biophysical testing showed that antigen-binding fragments (Fabs) from the aggregating antibodies also reversibly associated with equilibrium dissociation constants in the low-micromolar range. Crystal structures (PDB accession codes 6MXR, 6MXS, 6MY4, 6MY5) and bottom-up hydrogen-exchange mass spectrometry revealed that Fab self-association occurs in a symmetric mode that involves the antigen complementarity-determining regions. Subtle local conformational changes incurred upon point mutation of monomeric variants foster formation of complementary polar interactions and hydrophobic contacts to generate a dimeric Fab interface. Testing of popular in silico tools generally indicated low reliabilities for predicting the aggregation propensities observed. A structure-aggregation data set is provided here in order to stimulate further improvements of in silico tools for prediction of native aggregation. Incorporation of intermolecular docking, conformational flexibility, and short-range packing interactions may all be necessary features of the ideal algorithm.

Genomic screening of new putative antiviral lectins from Amazonian cyanobacteria based on a bioinformatics approach

Lectins are proteins of nonimmune origin, which are capable of recognizing and binding to glycoconjugate moieties. Some of them can block the interaction of viral glycoproteins to the host cell receptors acting as antiviral agents. Although cyanobacterial lectins have presented broad biotechnological potential, little research has been directed to Amazonian Cyanobacterial diversity. In order to identify new antiviral lectins, we performed genomic analysis in seven cyanobacterial strains from Coleção Amazônica de Cianobactérias e Microalgas (CACIAM). We found 75 unique CDS presenting one or more lectin domains. Since almost all were annotated as hypothetical proteins, we used homology modeling and molecular dynamics simulations to evaluate the structural and functional properties of three CDS that were more similar to known antiviral lectins. Nostoc sp. CACIAM 19 as well as Tolypothrix sp. CACIAM 22 strains presented cyanovirin‐N homologues whose function was confirmed by binding free energy calculations. Asn, Glu, Thr, Lys, Leu, and Gly, which were described as binding residues for cyanovirin, were also observed on those structures. As for other known cyanovirins, those residues in both our models also made favorable interactions with dimannose. Finally, Alkalinema sp. CACIAM 70d presented one CDS, which was identified as a seven‐bladed beta‐propeller structure with binding sites predicted for sialic acid and N‐acetylglucosamine. Despite its singular structure, our analysis suggested this molecule as a new putative antiviral lectin. Overall, the identification and the characterization of new lectins and their homologues are a promising area in antiviral research, and Amazonian cyanobacteria present biotechnological potential to be explored in this regard.

Application of Assisted Design of Antibody and Protein Therapeutics (ADAPT) improves efficacy of a Clostridium difficile toxin A single-domain antibody

Assisted Design of Antibody and Protein Therapeutics (ADAPT) is an affinity maturation platform interleaving predictions and testing that was previously validated on monoclonal antibodies (mAbs). This study expands the applicability of ADAPT to single-domain antibodies (sdAbs), a promising class of recombinant antibody-based biologics. As a test case, we used the camelid sdAb A26.8, a V_HH that binds Clostridium difficile toxin A (TcdA) relatively weakly but displays a reasonable level of TcdA neutralization. ADAPT-guided A26.8 affinity maturation resulted in an improvement of one order of magnitude by point mutations only, reaching an equilibrium dissociation constant (K_D) of 2 nM, with the best binding mutants having similar or improved stabilities relative to the parent sdAb. This affinity improvement generated a 6-fold enhancement of efficacy at the cellular level; the A26.8 double-mutant T56R,T103R neutralizes TcdA cytotoxicity with an IC₅₀ of 12 nM. The designed mutants with increased affinities are predicted to establish novel electrostatic interactions with the antigen. Almost full additivity of mutation effects is observed, except for positively charged residues introduced at adjacent positions. Furthermore, analysis of false-positive predictions points to general directions for improving the ADAPT platform. ADAPT guided the efficacy enhancement of an anti-toxin sdAb, an alternative therapeutic modality for C. difficile.

Antibody Affinity Maturation by Computational Design

The immune systems protect our bodies from foreign molecules or antigens, where antibodies play important roles. Antibodies evolve over time upon antigen encounter by somatically mutating their genome sequences. The end result is a series of antibodies that display higher affinities and specificities to specific antigens. This process is called affinity maturation. Recent improvements in computer hardware and modeling algorithms now enable the rational design of protein structures and functions, and several works on computer-aided antibody design have been published. In this chapter, we briefly describe computational methods for antibody affinity maturation, focusing on methods for sampling antibody conformations and for scoring designed antibody variants. We also discuss lessons learned from the successful computer-aided design of antibodies.

Assisted Design of Antibody and Protein Therapeutics (ADAPT)

Effective biologic therapeutics require binding affinities that are fine-tuned to their disease-related molecular target. The ADAPT (Assisted Design of Antibody and Protein Therapeutics) platform aids in the selection of mutants that improve/modulate the affinity of antibodies and other biologics. It uses a consensus z-score from three scoring functions and interleaves computational predictions with experimental validation, significantly enhancing the robustness of the design and selection of mutants. The platform was tested on three antibody Fab-antigen systems that spanned a wide range of initial binding affinities: bH1-VEGF-A (44 nM), bH1-HER2 (3.6 nM) and Herceptin-HER2 (0.058 nM). Novel triple mutants were obtained that exhibited 104-, 46- and 32-fold improvements in binding affinity for each system, respectively. Moreover, for all three antibody-antigen systems over 90% of all the intermediate single and double mutants that were designed and tested showed higher affinities than the parent sequence. The contributions of the individual mutants to the change in binding affinity appear to be roughly additive when combined to form double and triple mutants. The new interactions introduced by the affinity-enhancing mutants included long-range electrostatics as well as short-range nonpolar interactions. This diversity in the types of new interactions formed by the mutants was reflected in SPR kinetics that showed that the enhancements in affinities arose from increasing on-rates, decreasing off-rates or a combination of the two effects, depending on the mutation. ADAPT is a very focused search of sequence space and required only 20-30 mutants for each system to be made and tested to achieve the affinity enhancements mentioned above.

Anti-dengue virus activity of scytovirin and evaluation of point mutation effects by molecular dynamics and binding free energy calculations

The absence of a specific treatment against DENV has led to intensive research into developing strategies for curing the infection. One lectin with high antiviral activity is scytovirin, which was isolated from the cyanobacterium Scytonema varium and has proven activity against HIV and Zaire Ebola Virus. To achieve the results presented here, we tested the affinity of full-length scytovirin, SD1 and SD2 separately, and six SD1 mutants for DENV glycoprotein E carbohydrate by Molecular Dynamics (MD) simulations and binding free energy calculations. It was possible to identify the key residues for protein-ligand interaction such as Glu10, Ala11, Pro17, Ans18, Arg30, Thr41, Ser42 and Arg43, which also has importance action against HIV. All binding free energy calculations showed negative values to ΔGbind of protein-DENV carbohydrate complexation. Additionally, these results are similar to the values of scytovirin and HIV gp120 carbohydrate complexation (-32.20 kcal/mol). Furthermore, we found that SD1 individually has more affinity to the carbohydrate and the Asn9, Glu10, Asn18, Arg30 and Arg43 demonstrated an important role in this matter. We also found that mutant G48R has better affinity (-34.10 kcal/mol) for the DENV carbohydrate than the wild type protein (-27.15 kcal/mol).

Determination of equilibrium dissociation constants for recombinant antibodies by high-throughput affinity electrophoresis

High-quality immunoreagents enhance the performance and reproducibility of immunoassays and, in turn, the quality of both biological and clinical measurements. High quality recombinant immunoreagents are generated using antibody-phage display. One metric of antibody quality – the binding affinity – is quantified through the dissociation constant (K_D) of each recombinant antibody and the target antigen. To characterize the K_D of recombinant antibodies and target antigen, we introduce affinity electrophoretic mobility shift assays (EMSAs) in a high-throughput format suitable for small volume samples. A microfluidic card comprised of free-standing polyacrylamide gel (fsPAG) separation lanes supports 384 concurrent EMSAs in 30 s using a single power source. Sample is dispensed onto the microfluidic EMSA card by acoustic droplet ejection (ADE), which reduces EMSA variability compared to sample dispensing using manual or pin tools. The K_D for each of a six-member fragment antigen-binding fragment library is reported using ~25-fold less sample mass and ~5-fold less time than conventional heterogeneous assays. Given the form factor and performance of this micro- and mesofluidic workflow, we have developed a sample-sparing, high-throughput, solution-phase alternative for biomolecular affinity characterization.

A Rational Engineering Strategy for Designing Protein A-Binding Camelid Single-Domain Antibodies

Staphylococcal protein A (SpA) and streptococcal protein G (SpG) affinity chromatography are the gold standards for purifying monoclonal antibodies (mAbs) in therapeutic applications. However, camelid VHH single-domain Abs (sdAbs or VHHs) are not bound by SpG and only sporadically bound by SpA. Currently, VHHs require affinity tag-based purification, which limits their therapeutic potential and adds considerable complexity and cost to their production. Here we describe a simple and rapid mutagenesis-based approach designed to confer SpA binding upon a priori non-SpA-binding VHHs. We show that SpA binding of VHHs is determined primarily by the same set of residues as in human mAbs, albeit with an unexpected degree of tolerance to substitutions at certain core and non-core positions and some limited dependence on at least one residue outside the SpA interface, and that SpA binding could be successfully introduced into five VHHs against three different targets with no adverse effects on expression yield or antigen binding. Next-generation sequencing of llama, alpaca and dromedary VHH repertoires suggested that species differences in SpA binding may result from frequency variation in specific deleterious polymorphisms, especially Ile57. Thus, the SpA binding phenotype of camelid VHHs can be easily modulated to take advantage of tag-less purification techniques, although the frequency with which this is required may depend on the source species.